Difference between revisions of "Hysing benchmark by Lionel Gamet"

From OpenFOAM Wiki
Jump to navigation Jump to search
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]]
  
This case is a reference test case for VoF simulations. It is an elementary quantitative benchmark configuration that was originally established in 2D by Hysing et al. [1]. Adelsberger et al. [2] have published a 3D equivalent of the same benchmark. The case consists in a single rising bubble in a quiescent liquid. The test case number 2 as described by Hysing et al. [1] is modeled here.
+
==Introduction==
  
A video of the 2D case can be downloaded [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/alpha.phase1_2D320x640x001.avi here].
+
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.  
  
A video of the 3D case can be downloaded [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/alpha.phase1_3D320x640x320_colored_by_Umag.avi here].
 
  
The pdf for the 2D case can be found [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/HysingBenchmark_2D.pdf here].
+
==Setting up the test case in OpenFOAM==
  
The pdf for the 3D case can be found [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/HysingBenchmark_3D.pdf here].
+
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 files can be downloaded here:
+
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].
  
* [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/HysingBenchmark_2D_wikitut.zip 2D case]
 
* [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/HysingBenchmark_3D_wikitut.zip 3D case]
 
  
 
==References==
 
==References==
[1] S. Hysing, S. Turek, D. Kuzmin, N. Parolini, E. Burman, S. Ganesan, and L. Tobiska: "Quantitative benchmark computations of two-dimensional bubble dynamics," I.J.N.M.F., vol. 60, no. 11, pp. 1259{1288, 2009. [Online]. Available: http://www.featow.de/en/benchmarks/cfdbenchmarking/bubble.html
+
[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.
[2] J. Adelsberger, P. Esser, M. Griebel, S. Groß, M. Klitz, and A. Ruttgers: "3D incompressible two-phase flow benchmark computations for rising droplets," 2014, proceedings of the 11th World Congress on Computational Mechanics (WCCM XI), Barcelona, Spain, also available as INS Preprint No. 1401 and as IGPM Preprint No. 393. [Online]. Available: http://wissrech.ins.uni-bonn.de/research/projects/risingbubblebenchmark/
+
[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

Go back to Multiphase modeling.

Hysing benchmark

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 [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.

liquid phase fraction t=3s