On Tue, Jul 3, 2012 at 1:34 AM, Yun Tao <[email protected]> wrote: > Hm. Reducing the mesh range (l=50. -> l=30.) does seem to introduce > leaking even when the initial cell volume = 1 (timestep 485). Why would the > domain matter when the BCs are zero-flux?
The cell volumes averages do change, but very slowly. This could be due round off in the solver. The np.trapz seems to be very different from the cell average. Is that calculation taking into account that the values are cell centered so the end points may not be weighted correctly? -- Daniel Wheeler
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
