A VOF-LPT Coupling Method for Simulating Jet Atomization

doi:10.19596/j.cnki.1001-246x.8460

Abstract

Abstract:

We evaluate several VOF(Volume of Fluid) methods for interfacial reconstruction. A VOF-LPT(Lagrangian Particle Tracking) coupling method is applied in OpenFOAM, an open source CFD package, to achieve both simulation accuracy and computational efficiency. This method simulates accurately from continuous liquid to discrete small droplets for jet atomization. It provides a feasible and efficient simulation scheme for simulation of jet atomization process on large scales.

Key words: jet atomization, volume of fluid, multiphase flow, Lagrangian particle tracking

Zi-feng LI, Ning YANG. A VOF-LPT Coupling Method for Simulating Jet Atomization[J]. Chinese Journal of Computational Physics, 2022, 39(4): 440-452.

Figures/Tables 16

Fig.1 3D coupling algorithm (phase field threshold αt = 0.1, computing time t = ti) (a) initial phase field; (b) labeling liquid grid; (c) droplet identification and numbering

Fig.2 Algorithm logic diagram of VOF-LPT

Fig.3 Initial condition of dambreak simulation

Table 1 CPU time (in unit: s) of interface reconstruction methods (with 6 cores) at different mesh resolution

	粗网格(9070 Cells)	中等网格(36288 Cells)	细网格(145152 Cells)
MULES	20.2	483.2	5 032
Isoadvector	68.3	397.7	5 373
PLIC-RD	43.6	400.2	4 202

Fig.4 Dambreak simulation with different interface reconstruction method (a) medium grid; (b) fine grid

Fig.5 Phase field conservation error parameters in different interface capture methods

Fig.6 Mesh and boundary conditions in jet atomization simulation

Table 2 Computation time for 0.03 s physical time in jet atomization simulation (with 10 cores)

界面重构方法	计算时间/s
MULES	52 667
Isoadvector	39 326
PLIC-RDF	48 539

Fig.7 Jet atomization simulation with different interfacial reconstruction methods (a) phase field ɑ; (b) contours of ɑ = 0.5 colored by velocity Ux

Fig.8 Discrete droplets satisfying transformation threshold (a) VOF results before transformation; (b) droplet distribution after transformation

Fig.9 (a) Number density distribution in (Ux, Uz) plane; (b) Ux velocity distribution of discrete droplets

Fig.10 Initial conditions for the simulation of impact jet flow (D = 635 m, Uj=18.5 m ·s-1, 2α = 60°, We = 2 987, Re = 11 724)

Fig.11 VOF-LPT simulation of impact jet flow (a) pure VOF of Chen et al; (b) VOF-LPT in this work; (c) experiment of Ryan et al

Fig.12 Liquid structure before or after the identification and capture of discrete droplets (a) VOF; (b) VOF-LPT

Fig.13 (a) Droplet probability density distributions of discrete droplet velocity components; (b) Joint density distributions of pairwise combination of velocities of discrete droplets

Fig.14 Scale-dependent coupling of VOF, VOF-LPT and LPT simulation for jet atomization

References 32

1	DIWAN S S , RAVICHANDRAN S , GOVINDARAJAN R , et al. Understanding transmission dynamics of COVID-19-type infections by direct numerical simulations of cough/sneeze flows[J]. Transactions of the Indian National Academy of Engineering, 2020, 5, 255- 261. DOI
2	楼佳明. 气相法聚乙烯超冷凝工艺的开发[D]. 杭州: 浙江大学, 2012.
3	CHEN L , GHONIEM A F . Development of a three-dimensional computational slag flow model for coal combustion and gasification[J]. Fuel, 2013, 113, 357- 366. DOI
4	HOCHHAUS G , MöLLMANN H , DERENDORF H , et al. Pharmacokinetic/pharmacodynamic aspects of aerosol therapy using glucocorticoids as a model[J]. The Journal of Clinical Pharmacology, 1997, 37 (10): 881- 892. DOI
5	GHASEMI A, FUKUDA K, BALACHANDAR R, et al. Numerical investigation of spray characteristics of diesel alternative fuels[R]. SAE Technical Paper, 2012.
6	CAO Qiwei , XIAO Delong , YANG Xianjun , et al. Numerical simulation on nonlinear evolution of magneto-Rayleigh-Taylor instability[J]. Chinese Journal of Computational Physics, 2021, 38 (1): 5- 15.
7	TRYGGVASON G , BUNNER B , ESMAEELI A , et al. A front-tracking method for the computations of multiphase flow[J]. Journal of Computational Physics, 2001, 169 (2): 708- 759. DOI
8	DESJARDINS O , MOUREAU V , PITSCH H . An accurate conservative level set/ghost fluid method for simulating turbulent atomization[J]. Journal of Computational Physics, 2008, 227 (18): 8395- 8416. DOI
9	MÁRQUEZ D S. An extended mixture model for the simultaneous treatment of short and long scale interfaces[D]. Santa Fe: National University of the Littoral, 2013.
10	WAN Decheng . Numerical simulations of the collisions between two solitary waves and the transformation of a solitary wave over backward step by VOF method[J]. Chinese Journal of Computational Physics, 1998, 15 (6): 654- 666.
11	ZHAO Xizeng , YE Zhouteng . Numerical study of droplet splashing in a CIP-based model[J]. Chinese Journal of Computational Physics, 2016, 33 (1): 39- 48.
12	BRACKBILL J U , KOTHE D B , ZEMACH C . A continuum method for modeling surface tension[J]. Journal of Computational Physics, 1992, 100 (2): 335- 354. DOI
13	CHEN X , MA D , YANG V , et al. High-fidelity simulations of impinging jet atomization[J]. Atomization and Sprays, 2013, 23 (12): 1079- 1101. DOI
14	SPITZENBERGER A , NEUMANN S , HEINRICH M , et al. Particle detection in VOF-simulations with OpenFOAM[J]. SoftwareX, 2020, 11, 100382. DOI
15	OpenFOAM: The open source CFD toolbox[EB/OL]. http://www.OpenFOAM.com.
16	VALLIER A. Simulations of cavitation-from the large vapour structures to the small bubble dynamics[D]. Lund: Lund University, 2013.
17	HERRMANN M . A parallel Eulerian interface tracking/Lagrangian point particle multi-scale coupling procedure[J]. Journal of Computational Physics, 2010, 229 (3): 745- 759. DOI
18	HEINRICH M , SCHWARZE R . 3D-coupling of Volume-of-Fluid and Lagrangian particle tracking for spray atomization simulation in OpenFOAM[J]. SoftwareX, 2020, 11, 100483. DOI
19	DESMONS F , COQUERELLE M . A generalized high-order momentum preserving (HOMP)method in the one-fluid model for incompressible two phase flows with high density ratio[J]. Journal of Computational Physics, 2021, 437, 110322. DOI
20	ZHU J . On the higher-order bounded discretization schemes for finite volume computations of incompressible flows[J]. Computer Methods in Applied Mechanics and Engineering, 1992, 98 (3): 345- 360. DOI
21	BRACKBILL J U , KOTHE D B , ZEMACH C . A continuum method for modeling surface tension[J]. Journal of Computational Physics, 1992, 100 (2): 335- 354. DOI
22	WELLER H G , TABOR G , JASAK H , et al. A tensorial approach to computational continuum mechanics using object-oriented techniques[J]. Computers in Physics, 1998, 12 (6): 620- 631. DOI
23	WELLER H G, DERIVATION M. Solution of the conditionally averaged two-phase flow equations[C]. OpenCFD, Ltd, 2005.
24	LOPES P. Free-surface flow interface and air-entrainment modelling using OpenFOAM[D]. Coimbra: University of Coimbra, 2013.
25	ROENBY J , BREDMOSE H , JASAK H . A computational method for sharp interface advection[J]. Royal Society Open Science, 2016, 3 (11): 160405. DOI
26	OLSSON E . A description of isoadvector-a numerical method for improved surface sharpness in two-phase flows[J]. Proceedings of CFD with OpenSource Software, 2017,
27	SCHEUFLER H , ROENBY J . Accurate and efficient surface reconstruction from volume fraction data on general meshes[J]. Journal of Computational Physics, 2019, 383, 1- 23. DOI
28	GAMET L , SCALA M , ROENBY J , et al. Validation of volume-of-fluid OpenFOAM^® isoAdvector solvers using single bubble benchmarks[J]. Computers & Fluids, 2020, 213, 104722.
29	SCHILLER L , NAUMANN Z A . A drag coefficient correlation[J]. Zeitschrift des Vereins Deutscher Ingenieure, 1935, 77, 318- 320.
30	CHOUDHARY A N, THAKUR R. Evaluation of connected component labeling algorithms on shared and distributed memory multiprocessors[C]//Department of Electrical and Computer Engineering 121 Link Hall, Syracuse University Syracuse, 1992: 362-365.
31	ZUZIO D , ESTIVALÈZES J L , DIPIERRO B . An improved multiscale Eulerian-Lagrangian method for simulation of atomization process[J]. Computers & Fluids, 2018, 176, 285- 301.
32	RYAN H M , ANDERSON W E , PAL S , et al. Atomization characteristics of impinging liquid jets[J]. Journal of Propulsion and Power, 1995, 11 (1): 135- 145.

[1]	Xiangxiang WEI, Qihong FENG, Xianmin ZHANG, Yingsong HUANG, Lijie LIU. Distribution of Remaining Oil in Water Flooding at Pore Scale: Volume of Fluid Method [J]. Chinese Journal of Computational Physics, 2021, 38(5): 573-584.
[2]	NIU Xiao, NI Guoxi, MA Wenhua. Application of Adaptive Multi-resolution Method in Numerical Simulation of Reactive Multiphase Flows [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2020, 37(6): 639-652.
[3]	FENG Qihong, ZHAO Yunchang, WANG Sen, ZHANG Yigen, SUN Yeheng, SHI Shubin. Pore-scale Oil-Water Two-phase Flow Simulation Based on Phase Field Method [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2020, 37(4): 439-447.
[4]	LOU Qin, ZANG Chenqiang, WANG Haoyuan, LI Ling. Interfacial Dynamics of Immiscible Gas-Liquid Two-phase Flow for CO₂ in Microchannel: Lattice Boltzmann Method [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2019, 36(2): 153-164.
[5]	XIE Chiyu, ZHANG Jianying, WANG Moran. Lattice Boltzmann Modeling of Non-Newtonian Multiphase Fluid Displacement [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2016, 33(2): 147-155.
[6]	HONG Ning, LI Tiyun, LI Qiuxiang, SHI Baochang. Lattice Boltzmann Simulation of Multiphase Flows in a Square Cavity [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2013, 30(4): 520-526.
[7]	LIANG Xianhong. Second Order Volume of Fluid Interface Reconstruction Method in Three Dimensions [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2013, 30(3): 353-360.
[8]	MING Pingjian, DUAN Wenyang. Numerical Simulation of a Sway Tank with Unstructured Grid High Order Scheme [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2011, 28(4): 507-514.
[9]	ZHOU Lujun, XUAN Yimin, LI Qiang. Multiscale Simulation of Nanofluid Multiphase Flows [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2009, 26(6): 849-856.
[10]	XU Kan, LIU Minghou, LIU Dong, CHEN Yiliang. Simulation of Multiphase Flows in Microscale Hemodialysis [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2009, 26(1): 49-56.
[11]	WANG Ruyun, CHEN Pingping, BAN Changying. Keeping Volume Fraction of Fluid in Reconstructing Moving-interfaces of VOF on Rectangular Meshes [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2008, 25(4): 431-436.
[12]	SONG Yu, WANG Fei, HAO Pengfei, HE Feng. Numerical Simulation of Micro Flows Considering Wall Adhesion [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2008, 25(1): 75-82.