Hysing benchmark by Lionel Gamet - Revision history

Jozsef Nagy at 07:36, 26 January 2024

2024-01-26T07:36:40Z

Jozsef Nagy at 12:06, 30 August 2022

2022-08-30T12:06:03Z

Jozsef Nagy at 06:14, 3 December 2021

2021-12-03T06:14:00Z

Jozsef Nagy at 06:51, 17 September 2020

2020-09-17T06:51:56Z

Jozsef Nagy at 10:19, 24 August 2020

2020-08-24T10:19:00Z

Jozsef Nagy at 10:18, 24 August 2020

2020-08-24T10:18:50Z

Jozsef Nagy: /* Hysing benchmark */

2020-08-24T10:18:32Z

Hysing benchmark

Jozsef Nagy at 06:57, 3 August 2020

2020-08-03T06:57:14Z

Jozsef Nagy at 06:55, 3 August 2020

2020-08-03T06:55:17Z

Jozsef Nagy at 06:53, 3 August 2020

2020-08-03T06:53:34Z

@@ Line 4: / Line 4: @@
 * '''affiliation''': IFP Energies nouvelles, France
 * '''contact''': <mail address='lionel.gamet@ifpen.fr' description='author'>click here for email address</mail>
-* '''OpenFOAM version''': v2206
+* '''OpenFOAM version''': v2312
 * '''Published under''': CC BY-NC-SA license ([https://creativecommons.org/licenses creative commons licenses])

@@ Line 4: / Line 4: @@
 * '''affiliation''': IFP Energies nouvelles, France
 * '''contact''': <mail address='lionel.gamet@ifpen.fr' description='author'>click here for email address</mail>
-* '''OpenFOAM version''': v2106
+* '''OpenFOAM version''': v2206
 * '''Published under''': CC BY-NC-SA license ([https://creativecommons.org/licenses creative commons licenses])

@@ Line 4: / Line 4: @@
 * '''affiliation''': IFP Energies nouvelles, France
 * '''contact''': <mail address='lionel.gamet@ifpen.fr' description='author'>click here for email address</mail>
-* '''OpenFOAM version''': v2006
+* '''OpenFOAM version''': v2106
 * '''Published under''': CC BY-NC-SA license ([https://creativecommons.org/licenses creative commons licenses])

← Older revision		Revision as of 06:51, 17 September 2020
Line 31:		Line 31:
	* For [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/HysingBenchmark_2D_uns.zip 2D]: prismatic grids using either Pointwise or Gmsh		* For [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/HysingBenchmark_2D_uns.zip 2D]: prismatic grids using either Pointwise or Gmsh
	* For [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/HysingBenchmark_3D_uns.zip 3D]: tetrahedral grids using either Pointwise or Gmsh. There are also tetMesh and pMesh options for users willing to play with these kind of grids.		* For [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/HysingBenchmark_3D_uns.zip 3D]: tetrahedral grids using either Pointwise or Gmsh. There are also tetMesh and pMesh options for users willing to play with these kind of grids.
		+
		+	Reference OpenFOAM results of the Hysing benchmark can be found [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/HysingBenchmark_Results.tar.gz here].

@@ Line 4: / Line 4: @@
 * '''affiliation''': IFP Energies nouvelles, France
 * '''contact''': <mail address='lionel.gamet@ifpen.fr' description='author'>click here for email address</mail>
-* '''OpenFOAM version''': 2006
+* '''OpenFOAM version''': v2006
 * '''Published under''': CC BY-NC-SA license ([https://creativecommons.org/licenses creative commons licenses])

@@ Line 27: / Line 27: @@
 * [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==

@@ Line 11: / Line 11: @@
 =[https://github.com/jnmlujnmlu/OpenFOAMTeaching/tree/master/LionelGamet Hysing benchmark]=
-==Introduction==
+[[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 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.
+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.
-A video of the case can be downloaded [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/Bubble26_CanoLozano_video_isoAdvectorPlicRDF5_level4.avi here].
+A video of the 2D case can be downloaded [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/alpha.phase1_2D320x640x001.avi here].
-The starting case case can be downloaded [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/Bubble26_CanoLozano.tgz here].
+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 results can be found [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/LionelGamet/Bubble26_CanoLozano_Results.tgz 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 &rho;_1/rho;_2 between the fluids is 1000 and the dynamic viscosity ratio &mu;_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 = &rho;_1 g D_0^2= 10 and Galilei number Ga = (&rho;_1 g^{1/2} D_0^{3/2})/&mu;_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 &rho;_1 are taken as unity, which gives a surface tension &sigma; = 0.1Nm^{-1} and a rise velocity of the order of unity.
+* [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==
-[1] J. Cano-Lozano, C. Martínez-Bazán, 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.
+[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
-[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.
+[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/
 [[File:gamet2D.png|540px|center|liquid phase fraction t=3s]]

@@ Line 10: / Line 10: @@
 =[https://github.com/jnmlujnmlu/OpenFOAMTeaching/tree/master/LionelGamet Hysing benchmark]=
 ==Introduction==
@@ Line 56: / Line 54: @@
 ==Running the case==
 The Allrun_sensitivity script is at the highest level, above Allrun. It is used to run a grid sensitivity for all solvers, interFoam, interIsoFoam, and also the enhanced version of isoAdvector by H. Scheuer [3], now partially integrated in the v2006 with the plicRDF reconstruction scheme. Allrun_sensitivity should be run first to create running directories for the different VoF solvers. It will generate 6 running directories (3 solvers x 2 levels of grid refinement, as discussed in the former paragraph). The grids are the same for all solvers. Users might thus want to only generate the grid for one of the solvers (e.g. interFoam), and then create links in the others.

@@ Line 57: / Line 57: @@
 ==Running the case==
-The Allrun_sensitivity script is at the highest level, above Allrun. It is used to run a grid sensitivity for all solvers, interFoam, interIsoFoam, and also the enhanced version of isoAdvector by H. Scheuer [3], now partially integrated in the v2006 with the plicRDF reconstruction scheme. Allrun_sensitivity should be run �rst to create running directories for the di�erent VoF solvers. It will generate 6 running directories (3 solvers � 2 levels of grid re�nement, as discussed in the former paragraph). The grids are the same for all solvers. Users might thus want to only generate the grid for one of the solvers (e.g. interFoam), and then create links in the others.
+The Allrun_sensitivity script is at the highest level, above Allrun. It is used to run a grid sensitivity for all solvers, interFoam, interIsoFoam, and also the enhanced version of isoAdvector by H. Scheuer [3], now partially integrated in the v2006 with the plicRDF reconstruction scheme. Allrun_sensitivity should be run first to create running directories for the different VoF solvers. It will generate 6 running directories (3 solvers x 2 levels of grid refinement, as discussed in the former paragraph). The grids are the same for all solvers. Users might thus want to only generate the grid for one of the solvers (e.g. interFoam), and then create links in the others.
-The Allrun.pre script should be run first to create the mesh. The background grid is first constructed by running blockMesh, then snappyHexMesh is used to re�ne near the domain centerline. renumberMesh is used to reduce the matrix bandwith. It is a good practice to always run renumberMesh. Then the setAlphaField utility is run to initialize the bubble as a sphere.
+The Allrun.pre script should be run first to create the mesh. The background grid is first constructed by running blockMesh, then snappyHexMesh is used to refine near the domain centerline. renumberMesh is used to reduce the matrix bandwith. It is a good practice to always run renumberMesh. Then the setAlphaField utility is run to initialize the bubble as a sphere.
-The Allrun script should then be run to call the VoF solver. These computations are very long and might bene�t from AMR in future releases of the test cases. Results are provided for reference in the Wiki. Post-processing quantities
+The Allrun script should then be run to call the VoF solver. These computations are very long and might benefit from AMR in future releases of the test cases. Results are provided for reference in the Wiki. Post-processing quantities
-of interest (bubble volume, area, centroid, velocity and sphericity) are extracted from the solver log �le through a grep command.
+of interest (bubble volume, area, centroid, velocity and sphericity) are extracted from the solver log file through a grep command.
-The bubble trajectory is a 3D spiral, drawing almost a perfect circle in the x-y plane for the �ne grid case, as can be seen from �gure 1 taken from the article in press by Gamet et al. [4], with the standard isoAdvector solver
+The bubble trajectory is a 3D spiral, drawing almost a perfect circle in the x-y plane for the fine grid case, as can be seen from figure 1 taken from the article in press by Gamet et al. [4], with the standard isoAdvector solver
-interIsoFoam. A sample wake visualization, using the �2 vortex criterion, is shown on �gure 2 for the level 4 coarse grid.
+interIsoFoam. A sample wake visualization, using the &labda;2 vortex criterion, is shown on figure 2 for the level 4 coarse grid.
 ==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.
+[1] J. Cano-Lozano, C. Martínez-Bazán, 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]]