火箭推进

为研究驻波压力场速度波腹位置射流的雾化特性,基于八叉树结构形式的网格自适应方法与多尺度仿真算法构建了射流雾化的数值模拟方案,通过在计算域边界施加扰动构建了一阶横向驻波压力场。在此基础上实现了振荡压力场与雾化的多物理场耦合计算,研究了驻波压力场速度波腹位置射流雾化的响应特性,基于气体动力学理论阐述了速度波腹位置射流雾化的响应机理。结果表明:基于自适应网格与多尺度仿真算法建立的雾化数值模拟方案可以实现雾化过程较为准确的求解,是研究雾化过程的有力工具; 反压振荡与雾化的耦合求解算法可以研究振荡压力场下的雾化特性,对于进一步认识非定常雾化特性并揭示热声耦合机理起到重要的作用; 射流处于速度波腹位置时,圆柱射流发生变形变成扁平液膜,并伴随气流的运动发生周期性摆动,破碎长度减小,破碎程度加剧。

To investigate the characteristics of liquid jet atomization in the standing wave pressure field, a numerical scheme of liquid jet atomization was established based on octree adaptive mesh refinement and multiscale algorithm. A first-order transverse standing wave pressure field was established by exerting perturbations at the boundary of the computational zone. On this basis of the above methods, the multi-physical coupling simulation of pressure oscillation and atomization was carried out. The response characteristics of liquid jet atomization on velocity antinode position at the standing wave pressure field was studied. The response mechanism of liquid jet atomization was interpreted based on gas dynamics theory. The results reveal that the atomization process can be accurately solved by the established numerical scheme based on adaptive mesh refinement and multiscale algorithm, and it is a powerful tool to investigate the atomization process. The atomization characteristics coupled with pressure oscillations can be studied based on the coupled simulation algorithm between pressure oscillations and atomization. It will play an important role to further investigate the unsteady atomization process and the mechanism of thermo-acoustic positive feedback. When the jet is in the position of the velocity antinode, the liquid column jet is flattened to liquid sheet, acompanied by the periotic swing of liquid jet with the gas flow. The breakup length is shortened, and the breakup degree is enhanced.

引言
1 雾化过程的多尺度仿真算法
2 驻波压力场的激励方法
3 驻波压力场速度波腹位置射流雾化的响应分析
4 结论

表1 数值模拟参数设置<br/>Tab.1 Setup of simulation parameters

表1 数值模拟参数设置
Tab.1 Setup of simulation parameters

图1 计算域及边界条件的设定<br/>Fig.1 Computational zone and setup of boundary conditions

图1 计算域及边界条件的设定
Fig.1 Computational zone and setup of boundary conditions

图2 计算域的初始网格<br/>Fig.2 Initial meshes of computational zone

图2 计算域的初始网格
Fig.2 Initial meshes of computational zone

图3 PLIC-VOF方法捕捉到的球形液体<br/>Fig.3 Spherical droplet captured by PLIC-VOF method

图3 PLIC-VOF方法捕捉到的球形液体
Fig.3 Spherical droplet captured by PLIC-VOF method

图4 不同网格尺寸下的射流形态<br/>Fig.4 Jet morphology of different mesh dimensions

图4 不同网格尺寸下的射流形态
Fig.4 Jet morphology of different mesh dimensions

图5 宏观射流形态<br/>Fig.5 Macroscopic jet morphology

图5 宏观射流形态
Fig.5 Macroscopic jet morphology

图6 本文计算结果与文献[18]中仿真结果的对比<br/>Fig.6 Comparison of simulation results between this paper andRef.[18]

图6 本文计算结果与文献[18]中仿真结果的对比
Fig.6 Comparison of simulation results between this paper andRef.[18]

图7 发动机截面的一阶切向振型<br/>Fig.7 The first order tangential mode of engine cross-section

图7 发动机截面的一阶切向振型
Fig.7 The first order tangential mode of engine cross-section

图8 xy截面的压力分布云图<br/>Fig.8 Pressure distribution nephogram of xy-section

图8 xy截面的压力分布云图
Fig.8 Pressure distribution nephogram of xy-section

图9 圆柱形燃烧室中的一阶切向振型<br/>Fig.9 The first order tangential mode in cylindrical combustion chamber

图9 圆柱形燃烧室中的一阶切向振型
Fig.9 The first order tangential mode in cylindrical combustion chamber

图10 各监测平面平均压力随时间的变化<br/>Fig.10 Variation of mean pressure for each monitoring plane

图10 各监测平面平均压力随时间的变化
Fig.10 Variation of mean pressure for each monitoring plane

图11 0.12 ms时刻VOF捕捉到的射流形态 <br/>Fig.11 Jet morphology captured by VOF method at 0.12 ms

图11 0.12 ms时刻VOF捕捉到的射流形态
Fig.11 Jet morphology captured by VOF method at 0.12 ms

$图12 0.12 ms时刻不同截面的体积分数分布 <br/>Fig.12 Distribution of volume fraction on different sections at 0.12 ms$

图12 0.12 ms时刻不同截面的体积分数分布
Fig.12 Distribution of volume fraction on different sections at 0.12 ms

图13 不同时刻射流的动态响应过程<br/>Fig.13 Dynamic response process of the jet at different times

图13 不同时刻射流的动态响应过程
Fig.13 Dynamic response process of the jet at different times

$图14 不同时刻射流xy截面体积分数分布云图<br/>Fig.14 Volume fraction distribution of xy-section at different times$

图14 不同时刻射流xy截面体积分数分布云图
Fig.14 Volume fraction distribution of xy-section at different times

图15 0.095 ms时刻xy截面的流线图<br/>Fig.15 Streamline diagram of xy-section at 0.095 ms

图15 0.095 ms时刻xy截面的流线图
Fig.15 Streamline diagram of xy-section at 0.095 ms

图16 0.095 ms时刻射流轴线不同位置的横截面<br/>Fig.16 Cross sections of the jet for different positions at 0.095 ms

图16 0.095 ms时刻射流轴线不同位置的横截面
Fig.16 Cross sections of the jet for different positions at 0.095 ms

图17 0.12 ms时的射流雾场 <br/>Fig.17 Jet atomization field at 0.12 ms

图17 0.12 ms时的射流雾场
Fig.17 Jet atomization field at 0.12 ms

图18 射流在横向气流中变形破碎示意图<br/>Fig.18 Schematic diagram of jet deformation and fragmentation in transverse airflow

图18 射流在横向气流中变形破碎示意图
Fig.18 Schematic diagram of jet deformation and fragmentation in transverse airflow

[1] VOROB'EV A P. Effect of acoustic oscillations on the stability of a plane jet[J]. Fluid Dynamics, 1991, 26(4): 521-526.
[2]CARPENTIER J B, BAILLOT F, BLAISOT J B, et al. Behavior of cylindrical liquid jets evolving in a transverse acoustic field[J]. Physics of Fluids, 2009, 21(2): 23601.
[3]JU D H, SUN X C, JIA X X, et al. Experimental investigation of the atomization behavior of ethanol and kerosene in acoustic fields[J]. Fuel, 2017, 202: 613-619.
[4]JIA X X, HUANG Z, JU D H, et al. Effect of high frequency acoustic field on atomization behavior of ethanol and kerosene[Z]. 2017.
[5]贾晓旭. 高频声振荡对航空煤油和乙醇喷雾宏观特性的影响[D]. 上海: 上海交通大学, 2018.
[6]SRINIVASAN V, SALAZAR A J, SAITO K. Modeling the disintegration of modulated liquid jets using volume-of-fluid(VOF)methodology[J]. Applied Mathematical Modelling, 2011, 35(8): 3710-3730.
[7]SRINIVASAN V, SALAZAR A J, SAITO K. Numerical simulation of the disintegration of forced liquid jets using volume-of-fluid method[J]. International Journal of Computational Fluid Dynamics, 2010, 24(8): 317-333.
[8]YANG X C, TURAN A L. Simulation of liquid jet atomization coupled with forced perturbation[J]. Physics of Fluids, 2017, 29(2): 022103.
[9]POPINET S. Gerris: a tree-based adaptive solver for the incompressible Euler equations in complex geometries[J]. Journal of Computational Physics, 2003, 190(2): 572-600.
[10]POPINET S. An accurate adaptive solver for surface-tension-driven interfacial flows[J]. Journal of Computational Physics, 2009, 228(16): 5838-5866.
[11]谢远, 聂万胜, 高玉超, 等. 流量脉动对针栓式喷嘴雾化特性影响的仿真分析[J]. 火箭推进, 2023, 49(1): 54-64.
[12]CHEN X D, MA D J, YANG V, et al. High-fidelity simulations of impinging jet atomization[J]. Atomization and Sprays, 2013, 23(12): 1079-1101.
[13]WEI Q, LIANG G Z. Coupled Lagrangian impingement spray model for doublet impinging injectors under liquid rocket engine operating conditions[J]. Chinese Journal of Aeronautics, 2017, 30(4): 1391-1406.
[14]TOMAR G, FUSTER D, ZALESKI S, et al. Multiscale simulations of primary atomization[J]. Computers & Fluids, 2010, 39(10): 1864-1874.
[15]刘昌波, 雷凡培, 周立新. 两股湍流射流撞击雾化过程的数值研究[J]. 推进技术, 2014, 35(12): 1669-1678.
[16]丁群. 基于OpenFOAM的VOF-DDM模型对雾化现象的研究[D]. 哈尔滨: 哈尔滨工程大学, 2018.
[17]HIRT C W, NICHOLS B D. Volume of fluid(VOF)method for the dynamics of free boundaries[J]. Journal of Computational Physics, 1981, 39(1): 201-225.
[18]SHINJO J, UMEMURA A. Simulation of liquid jet primary breakup: dynamics of ligament and droplet formation[J]. International Journal of Multiphase Flow, 2010, 36(7): 513-532.
[19]QUINLAN J M, ZINN B. Simulating full-scale transverse combustion instabilities in a lab-scale facility[C]//49th AIAA/ASME/SAE/ASEE Joint Propulsion Conference. Reston, Virginia: AIAA, 2013.
[20]DUCRUIX S, REY C, CANDEL S. A method for the transverse modulation of reactive flows with application to combustion instability[J]. Combustion Theory and Modelling, 2005, 9(1): 5-22.
[21]DUCRUIX S, CANDEL S. External modulation of confined flows in computation fluid dynamics[J]. AIAA Journal, 2004, 42(8):1550-1558.
[22]HAKIM L, SCHMITT T, DUCRUIX S, et al. Dynamics of a transcritical coaxial flame under a high-frequency transverse acoustic forcing: influence of the modulation frequency on the flame response[J]. Combustion and Flame, 2015, 162(10): 3482-3502.
[23]GONZALEZ-FLESCA M, SCHMITT T, DUCRUIX S, et al. Large eddy simulations of a transcritical round jet submitted to transverse acoustic modulation[J]. Physics of Fluids, 2016, 28(5): 55-106.
[24]SHARIFI V, KEMPF A M, BECK C. Large-eddy simulation of acoustic flame response to high-frequency transverse excitations[J]. AIAA Journal, 2019, 57(1): 327-340.
[25]SCHMITT T, RODRIGUEZ J, LEYVA I A, et al. Experiments and numerical simulation of mixing under supercritical conditions[J]. Physics of Fluids, 2012, 24: 55-104.
[26]HARDI J S. Experimental investigation of high frequency combustion instability in cryogenic oxygen-hydrogen rocket engines[D]. Adelaide: The University of Adelaide, 2012.
[27] CLARK B J. Breakup of a liquid jet in a transverse flow of gas[R]. NASA TN D-2424, 1964.

备注

引言

1 雾化过程的多尺度仿真算法

2 驻波压力场的激励方法

3 驻波压力场速度波腹位置射流雾化的响应分析

4 结论

期刊信息