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?
On Wed, Jun 20, 2012 at 1:36 PM, Jonathan Guyer <[email protected]> wrote: > > On Jun 20, 2012, at 3:25 PM, Yun Tao wrote: > > > Thanks! Huge help! So, does phi.getCellVolumeAverage() * l act just like > integration under FiPy then? > > Yes. > > > > If so, for 2-D systems, would I simply do phi.getCellVolumeAverage() * l > * l ? > > phi.getCellVolumeAverage() * mesh.getCellVolumes().sum() is what you > should use in general. > > > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > -- Graduate Group of Ecology Doctoral Candidate Department of Environmental Science and Policy Center for Population Biology University of California, Davis
s-cd1dx.py
Description: Binary data
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
