Hi all, I am still bothering you with my simple diffusion problem between an outer circle (Boundary condition value=1) and an inner polygon (Boundary flux=0)
I am getting negative values on the concentration field !! a simple diffusion process ( eq = TransientTerm() == DiffusionTerm(coeff=D)) should not allow that , no ? even with a fixed flux boundary ?
