Re: Solving convection only problem

2016-10-20 Thread Zhekai Deng
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

RE: Question on time-marching scheme in Fipy

2016-10-20 Thread Campbell, Ian
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,

Re: Solving convection only problem

2016-10-20 Thread Daniel Wheeler
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

sweep with more than one variable

2016-10-20 Thread Francisco Vega Reyes
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).