[carlocal]

Dear all,

I am opening this thread here since I did not find any solutions elsewhere. For teaching purposes I am trying to implement the QUICK discretization scheme in the convection-diffusion equation - which is then solved iteratively - in MATLAB. However, I have huge stability problems that I cannot identify. I am leaving the entire code here in the hope that someone can take a look at it. Thanks!

[FMDenaro]

Originally Posted by carlocal

Dear all,

I am opening this thread here since I did not find any solutions elsewhere. For teaching purposes I am trying to implement the QUICK discretization scheme in the convection-diffusion equation - which is then solved iteratively - in MATLAB. However, I have huge stability problems that I cannot identify. I am leaving the entire code here in the hope that someone can take a look at it. Thanks!

Before I try to have a look to the code, I need more details.


1) is that a steady 1D problem with the known analytical solution?
2) did you try at low Peclet number?
3) what about your reference textbook for the scheme, the original paper of Leonard?