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

Unexpected results from interFoam

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

LinkBack Thread Tools Display Modes
Old   February 23, 2010, 17:39
Default Unexpected results from interFoam
New Member
Join Date: Jan 2010
Posts: 4
Rep Power: 9
mconroy is on a distinguished road
I am running the interFoam solver to do what I hoped would be a simple mulitphase simulation: filling a one-dimensional column with a mixture of water and air (without gravity).

When I run the simulation with the input liquid volume fraction equal to one, I get results that I expect (see filling tube plot).

However, when I changed the volume fraction to 0.75, I get weird oscillations in the liquid volume fraction as a function of height (see weird plot). I have tried changing the cell size and time steps, but I still had similar problems.

Boundary conditions:

Top of column (RHS in plots):
pressure = 1 atm
velocity -> zeroGradient
liquid volume fraction -> zeroGradient

Bottom of column (LHS in plots):
pressure -> zeroGradient
velocity = 10 cm/s
liquid volume fraction = {1 seems to work fine, <1 does not}

Initial conditions:
pressure = 1 atm
velocity = zero
liquid volume fraction = 0

I would like to understand what I am doing wrong and to hopefully get some help on how to get it right. Please let me know if any more information is needed.

Attached Files
File Type: pdf filling_tube.pdf (2.1 KB, 39 views)
File Type: pdf weird.pdf (3.7 KB, 32 views)
mconroy is offline   Reply With Quote

Old   March 30, 2011, 22:31
New Member
Dave West
Join Date: Mar 2011
Location: US
Posts: 12
Rep Power: 8
Terp is on a distinguished road
Mike -
I'm no expert on VOF methods, but my understanding is that it is looking for the interface between 2 fluids that don't mix. If you put in a fluid that is a uniform 75% phase 1 (and 25% phase 2), they are mixed to the sub-grid scale, and the code can't find a surface to separate them on. I wouldn't expect a good answer. If you bring in some of each phase, at 100% and 0%, separated to the scale of the mesh, I expect you would see better results.

Terp is offline   Reply With Quote

Old   March 31, 2011, 11:05
New Member
Join Date: Jan 2010
Posts: 4
Rep Power: 9
mconroy is on a distinguished road

Thanks for the reply. I need to understand the tools before I know what I can do with them

mconroy 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
Different Results from Fluent 5.5 and Fluent 6.0 Rajeev Kumar Singh FLUENT 6 December 19, 2010 12:33
Are the VOF results of interFoam grid independent asaha OpenFOAM Running, Solving & CFD 1 June 26, 2009 04:27
Is that a circle Results from interFoam sega OpenFOAM Paraview & paraFoam 2 May 8, 2008 02:49
Strange results from interFoam solution converges but sum of all forces not equal to zero nicasch OpenFOAM Running, Solving & CFD 0 April 15, 2008 02:01
unexpected results Kevin Rogers Main CFD Forum 7 October 6, 2005 11:34

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