Thanks! Huge help! So, does phi.getCellVolumeAverage() * l act just like integration under FiPy then? If so, for 2-D systems, would I simply do phi.getCellVolumeAverage() * l * l ?
On Wed, Jun 20, 2012 at 9:09 AM, Jonathan Guyer <[email protected]> wrote: > > On Jun 20, 2012, at 7:43 AM, Yun Tao wrote: > > > However, though the initial distribution integrates to 1 over the finite > domain, the area under the curve gradually decays as the function converges > on the origin (see attachment). This can be witnessed from the numerical > outputs starting at step 413 printed on the terminal. *The integration > function used is numpy.trapz.* Given the enclosure by the zero-flux > boundaries, I can't fathom how this can happen, whether it's a resolvable > issues or is conservation not guaranteed under certain conditions. > > x runs from the center of the leftmost cell to the center of the rightmost > cell, so you are not integrating the material between those cell centers > and the boundaries of the domain. If you print phi.getCellVolumeAverage() * > l, you will see that it doesn't leak. > > > _______________________________________________ > 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
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
