In this study, a solver coupling the level-set and volume-of-fluid methods is developed for interface computation of incompressible two-phase flows on the OpenFOAM platform. It corrects the surface tension, density, and viscosity of the fluids with the level-set function synchronously. To examine the effect of density and viscosity corrections, an intermediate solver that only corrects the surface tension force is also developed. Two test cases, that is, a bubble rising freely in stagnant water and a drop impacting on a liquid surface, are performed to validate and compare the capabilities of the new solvers. It is found that the coupling solver which corrects all the fluids’ properties mentioned above has the best prediction accuracy, especially for Rayleigh jet issues. With the application of this solver, the mechanism for the bubble entrainment in drop impact on a liquid pool is explained.
WuZ.Drop “impact” on an airfoil surface. Adv Colloid Interf Sci2018; 256: 23–47.
2.
VargasM. Droplet deformation prediction with the droplet deformation and breakup model (DDB).In: 4th AIAA Atmospheric and Space Environments Conference2012, New Orleans, Louisiana, 25–28 June 2012.
3.
WuZLvBCaoY. Modeling of drop deformation in proximity to an airfoil. Proc Inst Mech Eng C J Mech Eng Sci2018; 233: 1655–1659.
4.
WuZLvBCaoY. Improved Taylor analogy breakup and Clark models for droplet deformation prediction. Proc Inst Mech Eng G J Aerospace Eng2019; 233: 767–775.
5.
LopesP. Free-surface flow interface and air-entrainment modelling using OpenFOAM. PhD Thesis. Coimbra, Portugal: University of Coimbra, 2013.
6.
UbbinkO. Numerical prediction of two fluid systems with sharp interfaces. PhD Thesis. UK: Imperial College of Science, 1997.
7.
DalyBJ. A technique for including surface tension effects in hydrodynamic calculations. J Comput Phys1969; 4: 97–117.
8.
NicholsBDHirtCW. Calculating three-dimensional free surface flows in the vicinity of submerged and exposed structures. J Comput Phys1973; 12: 234–246.
9.
OsherSSethianJA. Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations. J Comput Phys1988; 79: 12–49.
10.
HarlowFHWelchJE. Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface. Phys Fluids1965; 8: 2182–2189.
11.
HirtCWNicholsBD. Volume of fluid (VOF) method for the dynamics of free boundaries. J Comput Phys1981; 39: 201–225.
12.
PopinetS. Numerical models of surface tension. Annu Rev Fluid Mech2018; 50: 49–75.
13.
PeskinCS. Flow patterns around heart valves: a numerical method. J Comput Phys1972; 10: 252–271.
14.
TryggvasonGBunnerBEsmaeeliA, et al.A front-tracking method for the computations of multiphase flow. J Comput Phys2001; 169: 708–759.
15.
PopinetSZaleskiS. A front‐tracking algorithm for accurate representation of surface tension. Int J Numer Methods Fluids1999; 30: 775–793.
16.
AndersonDMMcFaddenGBWheelerAA. Diffuse-interface methods in fluid mechanics. Annu Rev Fluid Mech1998; 30: 139–165.
17.
SussmanMSmerekaPOsherS. A level set approach for computing solutions to incompressible two-phase flow. J Comput Phys1994; 114: 146–159.
18.
BourliouxA. A coupled level-set volume-of-fluid algorithm for tracking material interfaces. In: Proceedings of the 6th International Symposium on Computational Fluid Dynamics, Lake Tahoe, CA, 4‐8, September 1995.
19.
AlbadawiADonoghueDBRobinsonAJ, et al.Influence of surface tension implementation in volume of fluid and coupled volume of fluid with level set methods for bubble growth and detachment. Int J Multiph Flow2013; 53: 11–28.
20.
HaghshenasMWilsonJAKumarR. Algebraic coupled level set-volume of fluid method for surface tension dominant two-phase flows. Int J Multiph Flow2017; 90: 13–28.
21.
LiuFXuYLiY. A coupled level-set and volume-of-fluid method for simulating axi-symmetric incompressible two-phase flows. Appl Math Comput2017; 293: 112–130.
22.
FerrariAMagniniMThomeJR. A flexible coupled level set and volume of fluid (flexCLV) method to simulate microscale two-phase flow in non-uniform and unstructured meshes. Int J Multiph Flow2017; 91: 276–295.
23.
NguyenTDHinoT. An improvement of interface computation of incompressible two-phase flows based on coupling volume of fluid with level-set methods. Int J Comput Fluid Dyn2020; 34: 75–89.
24.
WuZCaoYYaoY. Anatomy of air entrapment in drop impact on a solid surface. Int J Multiph Flow2021; 142: 103724.
GreenshieldsC. OpenFOAM user guide version 4.0. London: OpenFOAM Foundation Ltd, 2016.
27.
WellerHG. A new approach to VOF-based interface capturing methods for incompressible and compressible flow. Technical Report TR/HGW/04, 4. (2008, accessed 19 June 2012)
28.
BerberovićEvan HinsbergNPJakirlićS, et al.Drop impact onto a liquid layer of finite thickness: dynamics of the cavity evolution. Phys Rev E2009; 79: 036306.
29.
HemidaH. OpenFOAM tutorial: free surface tutorial using interFoam and rasInterFoam. Gothenburg: Chalmers University of Technology, 2008.
30.
SunDLTaoWQ. A coupled volume-of-fluid and level set (VOSET) method for computing incompressible two-phase flows. Int J Heat Mass Transf2010; 53: 645–655.
31.
HysingSTurekSKuzminD, et al.Quantitative benchmark computations of two-dimensional bubble dynamics. Int J Numer Methods Fluids2009; 60: 1259–1288.
32.
MortonDRudmanMJong-LengL. An investigation of the flow regimes resulting from splashing drops. Phys Fluids2000; 12: 747–763.
33.
RayBBiswasGSharmaA. Regimes during liquid drop impact on a liquid pool. J Fluid Mech2015; 768: 492–523.
34.
BhagaDWeberME. Bubbles in viscous liquids: shapes, wakes and velocities. J Fluid Mech1981; 105: 61–85.
35.
OguzHNProsperettiA. Bubble entrainment by the impact of drops on liquid surfaces. J Fluid Mech1990; 219: 143–179.
