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!
On 9 July 2013 17:27, Noam Yorav-Raphael <noamr...@gmail.com> wrote: > Hello, > > I'm trying to simulate radiation diffusion in an expanding matter. Since > the matter is expanding and the energy is bound to matter, it makes sense > to use mass coordinates, which are lagrangian. > > The equation turns out to be: ∂e/∂t = ρ/3 ∂²e/∂m² > where e is the energy density per volume, and ρ is mass per volume. The > total energy is ∫e/ρ dm, which is preserved. > > Is it possible to simulate this using fipy? > > By the way, I'm extremely impressed with fipy - I did a project using fipy > and a friend of mine used matlab. fipy was significantly faster than > matlab, produced better results, produced much more legible visualization > and with cleaner code - well done! > > Thank you very much, > Noam >
_______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]