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 fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]