Difference between revisions of "Hysing benchmark by Lionel Gamet"
Jozsef Nagy (talk | contribs) |
Jozsef Nagy (talk | contribs) |
||
Line 13: | Line 13: | ||
[[File:gamet3D.png|400px|right|isosurface alpha=0.5 colored by U t=3.5s]] | [[File:gamet3D.png|400px|right|isosurface alpha=0.5 colored by U t=3.5s]] | ||
− | + | ==Introduction== | |
− | + | This case is a reference test case for VoF simulations. It is the case number 26 of the 3D quantitative benchmark configuration by [1]. The case is also available as an example from the Basilisk website [http://basilisk.fr] and is detailed in the article of Cano-Lozano [1]. It consists in a single rising bubble in a large tank. In the physical conditions of this case, the rising bubble undergoes a spiralling path. | |
− | |||
− | + | ==Setting up the test case in OpenFOAM== | |
− | The | + | The bubble rises along +z direction and is initialized as a sphere at z_0/D_0 = 3:5, where D_0 is the bubble initial diameter. The fluid domain is of size 32x32x128 D_0. |
− | The case | + | The density ratio ρ_1/rho;_2 between the fluids is 1000 and the dynamic viscosity ratio μ_1/mu;_2 is 100. Index 1 refers to the continuous liquid phase while index 2 refers to the gas phase. The chosen Bond/Eotvos number Bo = ρ_1 g D_0^2= 10 and Galilei number Ga = (ρ_1 g^{1/2} D_0^{3/2})/μ_1 = 100.25 classify the current bubble in the oscillatory dynamics regime, with dominant inertial forces [2]. In the simulations, the gravity g and �first phase density ρ_1 are taken as unity, which gives a surface tension σ = 0.1Nm^{-1} and a rise velocity of the order of unity. |
+ | |||
+ | |||
+ | A video of the case can be downloaded [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/Bubble26_CanoLozano_video_isoAdvectorPlicRDF5_level4.avi here]. | ||
+ | |||
+ | The starting case case can be downloaded [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/Bubble26_CanoLozano.tgz here]. | ||
+ | |||
+ | The results can be found [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/Bubble26_CanoLozano_Results.tgz here]. | ||
− | |||
− | |||
==References== | ==References== | ||
− | [1] | + | [1] J. Cano-Lozano, C. Mart��nez-Baz�an, J. Magnaudet, and J. Tchoufag, Paths and wakes of deformable nearly spheroidal rising bubbles close to the transition to path instability," Physical Review Fluids, vol. 1, no. 5, 2016. |
− | + | [2] M. K. Tripathi, K. C. Sahu, and R. Govindarajan, Dynamics of an initially spherical bubble rising in quiescent liquid," Nature Communications, vol. 6, 2015. | |
− | [ | + | [3] H. Scheuer and J. Roenby, Accurate and effcient surface reconstruction from volume fraction data on general meshes," J. Comp. Phys., vol. 383, pp.1-23, 2019. |
+ | [4] L. Gamet, M. Scala, J. Roenby, H. Scheuer, and J.-L. Pierson, Validation of volume-of-fluid OpenFOAM isoAdvector solvers using single bubble benchmarks, Submitted to Computers and Fluids, 2020. | ||
[[File:gamet2D.png|540px|center|liquid phase fraction t=3s]] | [[File:gamet2D.png|540px|center|liquid phase fraction t=3s]] |
Revision as of 06:36, 3 August 2020
- contributor: Lionel Gamet
- affiliation: IFP Energies nouvelles, France
- contact: click here for email address
- OpenFOAM version: v1812
- Published under: CC BY-NC-SA license (creative commons licenses)
Go back to Multiphase modeling.
Contents
Hysing benchmark
Introduction
This case is a reference test case for VoF simulations. It is the case number 26 of the 3D quantitative benchmark configuration by [1]. The case is also available as an example from the Basilisk website [1] and is detailed in the article of Cano-Lozano [1]. It consists in a single rising bubble in a large tank. In the physical conditions of this case, the rising bubble undergoes a spiralling path.
Setting up the test case in OpenFOAM
The bubble rises along +z direction and is initialized as a sphere at z_0/D_0 = 3:5, where D_0 is the bubble initial diameter. The fluid domain is of size 32x32x128 D_0.
The density ratio ρ_1/rho;_2 between the fluids is 1000 and the dynamic viscosity ratio μ_1/mu;_2 is 100. Index 1 refers to the continuous liquid phase while index 2 refers to the gas phase. The chosen Bond/Eotvos number Bo = ρ_1 g D_0^2= 10 and Galilei number Ga = (ρ_1 g^{1/2} D_0^{3/2})/μ_1 = 100.25 classify the current bubble in the oscillatory dynamics regime, with dominant inertial forces [2]. In the simulations, the gravity g and �first phase density ρ_1 are taken as unity, which gives a surface tension σ = 0.1Nm^{-1} and a rise velocity of the order of unity.
A video of the case can be downloaded here.
The starting case case can be downloaded here.
The results can be found here.
References
[1] J. Cano-Lozano, C. Mart��nez-Baz�an, J. Magnaudet, and J. Tchoufag, Paths and wakes of deformable nearly spheroidal rising bubbles close to the transition to path instability," Physical Review Fluids, vol. 1, no. 5, 2016. [2] M. K. Tripathi, K. C. Sahu, and R. Govindarajan, Dynamics of an initially spherical bubble rising in quiescent liquid," Nature Communications, vol. 6, 2015. [3] H. Scheuer and J. Roenby, Accurate and effcient surface reconstruction from volume fraction data on general meshes," J. Comp. Phys., vol. 383, pp.1-23, 2019. [4] L. Gamet, M. Scala, J. Roenby, H. Scheuer, and J.-L. Pierson, Validation of volume-of-fluid OpenFOAM isoAdvector solvers using single bubble benchmarks, Submitted to Computers and Fluids, 2020.