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