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

Slug two-phase flow with InterFoam

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

LinkBack Thread Tools Display Modes
Old   October 23, 2015, 04:16
Default Slug two-phase flow with InterFoam
New Member
Francesco Tocci
Join Date: Sep 2015
Location: Delft, The Netherlands
Posts: 6
Rep Power: 2
Frenk_T is on a distinguished road
Hello Everybody,

I'm working with the solver interFoam (OpenFOAM-2.4.0), trying to reproduce the results shown in "Comparison of experimental and Computational Fluid Dynamics (CFD) studies of slug flow in a vertical riser" by M. Abdulkadir, V. Hernandez-Perez, S. Lo, I.S. Lowndes, B.J. Azzopardi.

I have test several parameters (mesh size, time step, maxCo, cAlpha, fvSchemes) but so far I'm not able to reproduce the results in terms of the rise velocity of Taylor bubble.

Maintaining a sharp and accurate interface is far from trivial. As you can see in the attached figures it seems that I can not avoid the effects of numerical smearing of the interface. The solution depends also on the initial condition (internal field) of k and epsilon (Fig. "instability" is obtained changing only the internal field for k,epsilon ---> Rayleigh-Taylor instability?).

Is there any problem with my settings?

Attached are an image of the mesh, snapshots and a tar.gz file with the 0, System folders and checkMesh output.

Any help is appreciated.


Attached Images
File Type: jpg Mesh.jpg (146.2 KB, 23 views)
File Type: jpg t_0_44s.jpg (112.9 KB, 21 views)
File Type: jpg t_1s.jpg (99.1 KB, 18 views)
File Type: jpg instability.jpg (32.9 KB, 18 views)
Attached Files
File Type: gz Riser.tar.gz (3.7 KB, 2 views)

Last edited by Frenk_T; October 23, 2015 at 04:20. Reason: Missing attachments
Frenk_T is online now   Reply With Quote

Old   October 23, 2015, 05:15
Senior Member
akidess's Avatar
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 983
Rep Power: 17
akidess will become famous soon enough
I think your mesh is far too coarse. In interFoam (VOF) you should strive to resolve every bubble, otherwise you need to use an Euler-Euler approach. Also, limit your timestep to very small Courant number (see Using standard k-epsilon for a multiphase simulation might not be ideal either, perhaps you can use LES.
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
*Help define the OpenFOAM stackexchange Q&A site:
akidess is offline   Reply With Quote

Old   October 23, 2015, 05:31
New Member
Francesco Tocci
Join Date: Sep 2015
Location: Delft, The Netherlands
Posts: 6
Rep Power: 2
Frenk_T is on a distinguished road
Dear Anton, thank you for your reply.
You are right about the mesh: however in the article that I attached, a similar mesh is employed with good results (they use Star-CCM+). So I was trying to replicate exactly the same case with OpenFoam. For the same reason I am using standard k-epsilon. Anyway I will try your advices.
Frenk_T is online now   Reply With Quote


interfoam bubble, multiphase, slug flow

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
Discrete Phase Model, outlet mass flow rate does not fit edu_aero FLUENT 28 September 18, 2015 06:53
Channel flow using InterFOAM DanM OpenFOAM Running, Solving & CFD 37 June 13, 2015 14:50
gas-liquid 2 phase flow in microchannels using VoF LyngHoo FLUENT 1 November 13, 2012 20:01
Open Channel Flow using InterFoam type solver sxhdhi OpenFOAM Running, Solving & CFD 3 May 5, 2009 21:58
Multiphase flow problem icedou FLUENT 6 July 10, 2005 02:52

All times are GMT -4. The time now is 07:45.