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

Checkerboarding with interFoam

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

LinkBack Thread Tools Display Modes
Old   July 14, 2015, 08:59
Default Checkerboarding with interFoam
Senior Member
Dr. Fabian Schlegel
Join Date: Apr 2009
Location: Dresden, Germany
Posts: 222
Rep Power: 11
fs82 is on a distinguished road
Dear Foamers,

I have a very annoying problem concerning boundary conditions in multiphase flows. I like to simulate a single twodimensional bubble with interFoam. Boundary conditions are zero velocity condition at top and buttom and slip condition at the left and right wall. Alpha has zeroGradient everywhere and p_rgh fixedFluxPressure.

To generate an appropriate start solution I turned of gravity and gave the simulation some time to generate a steady solution, assuming that the bubble does not move due to missing gravity. However, after a couple of seconds my simulation crashes. The reason seems to be some oscillations (checkerboarding) visible in the horizontal velocity U at the left or right wall (velocity_slip.png). The pressure field looks strange as well (not symmetric, pressue_slip.png). I do not have a clue what is the reason for that. If I change slip to fixedValue everything is fine (velocity_fixedValue.png and pressure_fixedValue.png). I attached a testcase, if somebody wants to try it. I appreciate every help.

Best regards,
Attached Images
File Type: jpg pressure_fixedValue.jpg (27.0 KB, 22 views)
File Type: jpg pressure_slip.jpg (30.3 KB, 23 views)
File Type: jpg velocity_fixedValue.jpg (29.2 KB, 24 views)
File Type: jpg velocity_slip.jpg (33.2 KB, 22 views)
Attached Files
File Type: zip (6.6 KB, 1 views)
fs82 is offline   Reply With Quote

Old   July 15, 2015, 07:19
Senior Member
Dr. Fabian Schlegel
Join Date: Apr 2009
Location: Dresden, Germany
Posts: 222
Rep Power: 11
fs82 is on a distinguished road
Heureka, the solution is written in Ferziger, Computational Methods for Fluid Dynamics, 2008 on page 294 (section 8.8) and is related to the rhie-chow correction to avoid the checkerboarding on staggered grids. There is a minimum criteria for the Courant number of Co>0.01 and I missed it.
fs82 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
interFoam in parallel gooya_kabir OpenFOAM Running, Solving & CFD 0 December 9, 2013 06:09
Problem of InterFoam with LES SpalartAllmarasIDDES keepfit OpenFOAM 3 August 29, 2013 11:21
InterFoam stops after deltaT goes to 1e14 francesco_b OpenFOAM Running, Solving & CFD 8 July 31, 2013 02:29
interFoam vs. simpleFoam channel flow comparison DanM OpenFOAM Running, Solving & CFD 11 January 5, 2013 07:21
Open Channel Flow using InterFoam type solver sxhdhi OpenFOAM Running, Solving & CFD 3 May 5, 2009 21:58

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