Thanks! Now I managed to get the expected result. It took me some time to realize that since ρ is time-dependent I had to use hasOld=True and updateOld(). I wonder if it's possible for fipy to detect that the TransientTerm coefficient has changed and issue an informative warning.
Thanks again for fipy, Noam On 10 July 2013 16:32, Daniel Wheeler <[email protected]> wrote: > > > > On Wed, Jul 10, 2013 at 9:02 AM, Noam Yorav-Raphael <[email protected]>wrote: > >> I'm sorry, it turns out I was wrong about the equation. The correct >> equation is ∂(e/ρ)/∂t = 1/3 ∂²e/∂m², which is solvable by using >> TransientTerm(coeff=1/rho). >> Thanks! >> > > I hope that things work out for you and thanks for the kind words > regarding fipy. > > -- > Daniel Wheeler > > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > >
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
