I think FiPy is really badly suited for coupled hyperbolic equations. It is essential to have a "Riemann-type solver" for this type of problem to correctly get the cross-coupling in the waves across the cell faces. I would really recommend using CLAWPACK for this. FiPy is much better suited for convection-diffusion problems and not so good for pure hyperbolic equations. It would be great to have that functionality in FiPy, but unfortunately it isn't available.
On Fri, Nov 28, 2014 at 2:58 PM, Dan Brunner <[email protected]> wrote: > > Thanks for you help so far. I will be more explicit in what I am doing. The > equations are: > > > With n, v, Te, and Ti the variables, x the spatial dimension, t the temporal > dimension, and m a constant. The S-terms are source terms, which will get > more complicated as I progress (including a thermal conductivity with a > T-dependent coefficient). To solve, for now ignoring the terms with m in the > last equation for simplicity, I have tried recasting them in the form: > > I roughly follow the process laid out in the coupled diffusion equations > example: -- Daniel Wheeler _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
