CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Convective Schemes how are they implemented

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   February 28, 2008, 11:44
Default Dear OF users and developers,
New Member
Join Date: Mar 2009
Posts: 23
Rep Power: 10
nicasch is on a distinguished road
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:

Φf =ΦC +(xf – xC) *ψ(rf) * (dΦ/ dx)f
Φf =ΦC +(xf – xC) / (xD– xC) *ψ(rf) * (ΦD-ΦC)

f ,the face
C,central node
D downwind node
ψ(rf) the flux limiter

rf is based on the following formula
rf=(ΦC-ΦU)/( ΦD-ΦC) = 2*grad(ΦC)*(xD-xC) /( ΦD-Φ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.
nicasch is offline   Reply With Quote

Old   November 4, 2010, 09:47
Senior Member
santiagomarquezd's Avatar
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 438
Rep Power: 17
santiagomarquezd will become famous soon enough
Nicasch, have you finally found an explanation for interfaceCompression scheme? I'm struggling with the same question.

Research Scientist
Research Center for Computational Methods (CIMEC) - CONICET/UNL
Tel: 54-342-4511594 Int. 7032
Colectora Ruta Nac. 168 / Paraje El Pozo
(3000) Santa Fe - Argentina.
santiagomarquezd is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
GammaContactAngle not implemented coastal593 OpenFOAM Running, Solving & CFD 5 November 9, 2010 06:21
Convective schemes for vector fields edoardo OpenFOAM Running, Solving & CFD 6 February 16, 2009 12:05
Help for fourth order accurate convective schemes Z.C.Wang Main CFD Forum 0 January 15, 2009 07:53
Where is gSum implemented cosimobianchini OpenFOAM Running, Solving & CFD 0 September 19, 2007 12:35
Approximation Schemes for convective term Michail CFD-Wiki 10 July 31, 2006 14:33

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