Thanks for the reply. I searched more and just realized this question has
been answered before. I then switched to VanLeerConvectionTerm, It turns
out this method is able to keep the shock pretty well which has been
mentioned in the previous posts. I might try CLAWPACK in the future, but I
think
Hi Daniel,
With regards to your answer:
"It's fully implicit unless the user sets it up in a different way"
Just to confirm what I understand from this and from your 2009 paper, alpha in
FiPy is 1.0 by default, and so the fully implicit time-stepping scheme that's
used is backward-Euler,
Hi Zhekai,
There is generally a lot of numerical diffusion when solving
convection problems with first order schemes and even some numerical
diffusion when using higher order schemes. There are many different
schemes and a mass of literature on how to preserve square waves,
shocks and hyperbolic
Hello,
Is it possible to express in a simple way the sweep method for more than
one variable?
More concretely, in a similar problem to the Stokes Cavity example I got
an equation with DiffusionTerm both for the X velocity and the Y
velocity (in the same equation, due to tensorial viscosity).