Hi Everyone, I have a diffusive field c2 defined by the equation :
c2eq= TransientTerm() == DiffusionTerm(coeff=D) the boundary conditions are defined by a source term (value c20 )in a circular region : and a null value on the outside, X, Y = mesh.getFaceCenters() XX = FaceVariable(mesh=mesh, value=X) YY = FaceVariable(mesh=mesh, value=Y) sourcec2 = ((XX-CX)**2+(YY-CY)**2 <= Rstart**2 ) c2.constrain(c20, where=sourcec2) BCsforc2 = FixedValue(faces=mesh.getExteriorFaces(), value=0) c2eq.solve(var=c2, dt=dt,boundaryConditions=BCsforc2,solver = LinearLUSolver(tolerance = 1.e-15, steps = 2000)) as the equation is just diffusion, I cannot get a value higher than c20! but when I plot c2, I find a value of about 3*c20 is there a problem with constrain ? thanks so much for your help, Julien
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy]
