CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Why results differs between upwind and vanLeerV (http://www.cfd-online.com/Forums/openfoam-solving/66741-why-results-differs-between-upwind-vanleerv.html)

kjetil July 23, 2009 04:58

Why results differs between upwind and vanLeerV
 
2 Attachment(s)
Dear community,
I have been experimenting with different schemes for a pipe case, and have been running upwind and vanLeerV for the initial phase for two almost identical cases. A two phase flow w/ alpha = 0.5 with different velocities is used.

In setup 1 I used this setup for div schemes:
divSchemes
{
default none;
div(phia,Ua) Gauss upwind;
div(phib,Ub) Gauss upwind;
div(phib,k) Gauss upwind;
div(phib,epsilon) Gauss upwind;
div(phi,alpha) Gauss limitedLinear01 1;
div((-nuEffa*grad(Ua).T())) Gauss linear;
div((-nuEffb*grad(Ub).T())) Gauss linear;
}
... and for setup 2:
divSchemes
{
default none;
div(phia,Ua) Gauss vanLeerV 1;
div(phib,Ub) Gauss vanLeerV 1;
div(phib,k) Gauss upwind;
div(phib,epsilon) Gauss upwind;
div(phi,alpha) Gauss limitedLinear01 1;
div((-nuEffa*grad(Ua).T())) Gauss linear;
div((-nuEffb*grad(Ub).T())) Gauss linear;
}
But the results from the two setups are quite different, and I believe the only difference in my setup is the different schemes. After 10 seconds, the pipe exit looks very different in the two setups - see attached photos. And it is the vanLeer-case which is "delayed". Does anybody have an explanation to this? I suppose the two schemes should yield the same results?

idrama November 10, 2012 05:55

Hallo kjetil,

the difference between upwind and vanLeer is, that vanLeer is a scheme composed of upwind and central differencing. Of course, you can have different results, because upwind is alway first order accurate whereby vanLeer can have second order accuracy - when no step gradients occur. For a better understanding, I recommend to read Versteeg (Chapter about upwind, QUCIK etc.). He does a very good analysis about the different schemes.

Cheers


All times are GMT -4. The time now is 01:18.