Slug two-phase flow with InterFoam

Frenk_T · October 23, 2015, 04:16

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.

Thanks,

Francesco

akidess · October 23, 2015, 05:15

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 http://dx.doi.org/10.1088/1749-4699/5/1/014016). Using standard k-epsilon for a multiphase simulation might not be ideal either, perhaps you can use LES.

Frenk_T · October 23, 2015, 05:31

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 · January 21, 2016, 13:03

Just to give an update:
I changed solver and with multiphaseEulerFoam I have interesting results. I attach a screenshot; I can see the formation of a Taylor bubble even with a really coarse mesh. The aspect ratio plays an important role (the two meshes have the same number of cells in the z direction).

October 23, 2015, 05:15		#2
akidess Senior Member Anton Kidess Join Date: May 2009 Location: Germany Posts: 1,377 Rep Power: 29	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 http://dx.doi.org/10.1088/1749-4699/5/1/014016). 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.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Channel flow using InterFOAM	DanM	OpenFOAM Running, Solving & CFD	49	July 31, 2020 11:43
Discrete Phase Model, outlet mass flow rate does not fit	edu_aero	FLUENT	29	February 3, 2020 08:38
gas-liquid 2 phase flow in microchannels using VoF	LyngHoo	FLUENT	1	November 13, 2012 19: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

October 23, 2015, 05:31		#3
Frenk_T New Member Join Date: Sep 2015 Location: Germany Posts: 11 Rep Power: 10	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.