CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

interFoam - mass conservation

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

Like Tree3Likes
  • 3 Post By carsten

LinkBack Thread Tools Display Modes
Old   August 5, 2009, 16:00
Exclamation interFoam - mass conservation
New Member
Wolfgang Betz
Join Date: Mar 2009
Posts: 6
Rep Power: 8
wbetz is on a distinguished road

Why isn't the boundary condidition 'zeroGradient' for pressure on walls not working anymore in OpenFOAM 1.6 (interFoam)? (if you do, there is a problem with mass conservation ...)

type zeroGradient;

... instead you have to use: (why ???)

type buoyantPressure;
value uniform 0;

wbetz is offline   Reply With Quote

Old   August 6, 2009, 01:40
Senior Member
sega's Avatar
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
I suppose this is due to the replacement of pd with static pressure p.
Have a look into the Release Notes for version 1.6 (Section 5.1 Solver restructuring):

pd replaced by static pressure p. All solvers in which buoyancy affects might be strong have been converted from using pd to p with improved numerics to give equally good accuracy and stability. This change is prompted by the need to remove the confusion surrounding the meaning and purpose of pd.
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   August 6, 2009, 03:51
New Member
Wolfgang Betz
Join Date: Mar 2009
Posts: 6
Rep Power: 8
wbetz is on a distinguished road
I agree with that - however, does somebody has more information on that behaviour ('zeroGradient' for pressure on walls ist not working anymore)?

And are there any other BCs you have to be aware of (as to a similar behaviour) in OF-version 1.6?
wbetz is offline   Reply With Quote

Old   September 21, 2009, 10:40
Carsten Thorenz
Join Date: Mar 2009
Location: Germany
Posts: 34
Rep Power: 8
carsten is on a distinguished road
Dear Wolfgang,

due to the nature of the PISO-algorithm used in interFoam, it is necessary to have a very "good" pressure field, as the pressure field is required in the velocity correction step of the scheme in order to reach mass conservation (i.e.div(U)=0).

If you look at the PISO-algorithm you will see that this means if the flux is zero (as it is orthogonal to the wall) it is necessary that the second derivative of pressure (orthogonal to the wall) must be zero, too (in areas where the density is constant).

buoyantPressure fulfills this requirement in areas of constant VOF-function / density. I'm not sure about mass conservation if sloped walls touch areas with VOF-values 0<VOF<1.

Bye and good luck,

Carsten Thorenz
maddalena, mgg and Pirlu like this.
carsten is offline   Reply With Quote

Old   December 7, 2010, 16:22
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 9
idrama is on a distinguished road
Has anybody ever done a mass conservation analysis in a pipe flow within VOF approach?
idrama is offline   Reply With Quote


buoyantpressure, interfoam, walls, zerogradient

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
InterFoam mass conservation question msrinath80 OpenFOAM Running, Solving & CFD 5 December 9, 2013 02:19
Mass Conservation in LES Sas Main CFD Forum 3 July 6, 2007 10:37
Mass Conservation in LES Jaswant Main CFD Forum 8 July 4, 2005 22:32
Mass conservation Diego Nogueira Main CFD Forum 1 July 30, 2004 15:50
Two-phase mass conservation ? Michael R. Rasmussen CFX 3 September 1, 2000 08:36

All times are GMT -4. The time now is 21:32.