CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Convective Schemes how are they implemented (http://www.cfd-online.com/Forums/openfoam-solving/59089-convective-schemes-how-they-implemented.html)

 nicasch February 28, 2008 11:44

Dear OF users and developers,

Dear OF users and developers,

I would be very grateful if anyone could please provide some insight into the following divSchemes:

- Gauss limitedLinearV
- Gauss interfaceCompression.

I suppose that they do some kind of blending based on the blending factor (limiter) calculated according to the scheme used. Is this blending between linear UD and CD or something else?

I found a post containing the following:

&Phi;f =&Phi;C +(xf &ndash; xC) *&psi;(rf) * (d&Phi;/ dx)f
Or
&Phi;f =&Phi;C +(xf &ndash; xC) / (xD&ndash; xC) *&psi;(rf) * (&Phi;D-&Phi;C)

f ,the face
C,central node
D downwind node
&psi;(rf) the flux limiter

rf is based on the following formula
rf=(&Phi;C-&Phi;U)/( &Phi;D-&Phi;C) = 2*grad(&Phi;C)*(xD-xC) /( &Phi;D-&Phi;C) -1

I get bit confused while trying to read the code, is there any written text about these two schemes, or can someone please explain how they work?

In limitedLinearV, I suppose that it is a kind of limited CD, but how is the scheme limited? Where is the limiter for this scheme used? Is the limiter function phi(r) used in calculation of the face value using adjacent cell-center values (like normal CD), or is it a blending between linear upwind and CD?

For the interfaceCompression scheme I suppose that it is some blending based on the blending factor calculated according to the either quartic or quadratic formula. Where is the limiter for this scheme used?

Many thanks in advance, best regards.

 santiagomarquezd November 4, 2010 09:47

Nicasch, have you finally found an explanation for interfaceCompression scheme? I'm struggling with the same question.

Regards

 All times are GMT -4. The time now is 09:30.