CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

Numerical schemes for free surface flows (VOF)

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

Like Tree5Likes
  • 1 Post By botp
  • 3 Post By santiagomarquezd
  • 1 Post By idrama

LinkBack Thread Tools Search this Thread Display Modes
Old   March 11, 2011, 08:08
Lightbulb Numerical schemes for free surface flows (VOF)
New Member
Bo Terp Paulsen
Join Date: Oct 2010
Posts: 13
Rep Power: 14
botp is on a distinguished road
Dear all,

I am currently working with free surface flows using the interFoam solver.

Unfortunately I have some serious problems with violent instabilities in the air phase and wondered however this is related to the numerical schemes? An Example of the instabilities are shown in this figure. In the figure the velocity field is shown with a contour of the VOF-scalar (gamma=0.5) superimposed. As seen from the figure, vortices are created in the air and the free surface surface is highly distorted. The distortion of the free surface is non-physical and ruin the computations.

So, if anyone has experience with the numerical schemes in relation to the interFoam solver and/or the VOF scheme in general, I will appreciate some advice!

Currently I am using the settings from the damBreak tutorial.

Kind regards,
Bo Terp
ascotdreamer likes this.

Last edited by botp; March 11, 2011 at 09:57. Reason: Additional information added
botp is offline   Reply With Quote

Old   March 11, 2011, 11:28
Senior Member
santiagomarquezd's Avatar
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 452
Rep Power: 23
santiagomarquezd will become famous soon enough
Hi Bo Terp, I had same problems solving an sloshing case. I've been solving other cases too and my suggestion would be run with this parameters:

-upwind difference scheme for all div terms.
-5 PISO correctors
-nAlphaCorr no more than 1
-nAlphaSubCycles, no more than 1
-Courant number below or equal to 0.2 (use adjustTimeStep to do that)

once you have a very diffusive but bounded solution both in alpha and U, then try increasing cAlpha and changing the div schemes to which are used in interFoam tutorials.

Let me know of your findings, I'm learning too.

shipman, randolph and vivek05 like this.
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

Old   March 11, 2011, 15:27
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 17
idrama is on a distinguished road
Properly you have an unstructured mesh which contains non-orthogonal cells. Can you post the current one? I don't know the exact boundary conditions but I would suggest you set something like

nCorrectors 2;
nOrthoCorrectors 4;

div(rho*phi,U) Gauss limitedLinearV 1;
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss interfaceCompression 1;

Alternatively, I can recommend
div(phi,alpha) Gauss limitedLinear 1;
div(phi,alpha) Gauss linearUpwind cellMDLimited Gauss linear 1;

Setting cAlpha 0 means no interface compression, take cAlpha 1. Hasn't that work than increase to 2,3,4 etc.

tikulju likes this.

Last edited by idrama; March 12, 2011 at 02:17.
idrama is offline   Reply With Quote


numerical scheme, openfoam 1.5-dev, vof model

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
free surface of VOF and melting model? wanghong FLUENT 3 March 13, 2006 09:57
Standard for checking and testing numerical schemes? X. Ye Main CFD Forum 7 August 31, 1999 17:05
Free Surface Flows / Ocean Waves / Far-field boundary conditions Sundar Prasad Main CFD Forum 5 May 4, 1999 10:25
free surface flows of submerged objects yogesh amle Main CFD Forum 1 March 16, 1999 02:05
Modeling Free Surface Flows Elliot Schwartz Main CFD Forum 5 August 25, 1998 21:03

All times are GMT -4. The time now is 15:10.