煤油液滴直径对两相旋转爆轰发动机流场的影响

杨帆; 姜春雪; 王宇辉; 李世全; 王健平; 张国庆

doi:10.11883/bzycj-2022-0068

煤油液滴直径对两相旋转爆轰发动机流场的影响

doi: 10.11883/bzycj-2022-0068

1.
北京化工大学机电工程学院，北京 100029
2.
北京大学工学院，北京 100871
3.
北京理工大学宇航学院，北京 100081

基金项目: 国家自然科学基金（52076003）；中央高校基本科研业务费专项资金（buctrc201913）

详细信息

作者简介:
杨　帆（1996－　），男，硕士研究生，yfhbxt@163.com

通讯作者:
王宇辉（1986－　），男，副教授，aowuki@163.com

中图分类号: O382; V231
计量
- 文章访问数: 661
- HTML全文浏览量: 213
- PDF下载量: 78
- 被引次数: 4
出版历程
- 收稿日期: 2022-02-25
- 修回日期: 2022-10-12
- 网络出版日期: 2022-10-13
- 刊出日期: 2023-02-05

Influence of kerosene droplet diameters on the flow field of a two-phase rotating detonation engine

1.
College of Mechanical and Electrical Engineering, Beijing University of Chemical Technology, Beijing 100029, China
2.
College of Engineering, Peking University, Beijing 100871, China
3.
School of Aerospace Engineering, Beijing Institute of Technology, Beijing 100081, China

摘要

摘要: 为探究煤油液滴不同初始直径对气液两相旋转爆轰发动机流场的影响，假设初始注入的煤油液滴具有均匀直径，考虑雾化破碎、蒸发等过程，建立了非定常两相爆轰的Eulerian-Lagrangian模型，进行了液态煤油/高温空气爆轰的非预混二维数值模拟。结果表明：在初始液滴直径为1~70 μm的工况范围，燃烧室内均形成了单个稳定传播的旋转爆轰波；全局当量比为1时，爆轰波前的空气区域大于液滴煤油的蒸气区域，导致波前燃料空气混合不均匀，波前均存在富油区和贫油区，两相速度差导致分离出的空气形成低温条带；当煤油液滴的初始直径较小时，波前的反应物混合过程主要受蒸发的影响，爆轰波可稳定传播；当直径减小至1 μm时，煤油液滴在入口处即蒸发，旋转爆轰波表现为气相传播的特性，爆轰波结构平整；当煤油液滴的初始直径较大时，波前的反应物混合过程主要受液滴破碎的影响；对于相同的燃料质量流量，在不同初始煤油液滴直径工况下，煤油液滴最大的停留时间均占爆轰波传播时间尺度的80%以上；爆轰波前燃料预蒸发为气相的占比越高，爆轰波的传播速度越高；初始液滴直径为10~70 μm的工况范围内，爆轰波的速度随初始直径的增大先升高后降低。
- 旋转爆轰发动机 /
- 气液两相流 /
- 燃烧 /
- 液滴直径
Abstract: To investigate the influence of the initial droplet diameter on the flow field of gas-liquid two-phase rotating detonation engine, an Eulerian-Lagrangian model of unsteady two-phase detonation is established based on the assumption of an initially uniform droplet diameter and considering atomization and evaporation processes. Non-premixed two-dimensional numerical simulations of detonation for liquid kerosene and high temperature air mixture are conducted. The results show that a single stable rotating detonation wave is formed in the initial droplet diameter range of 1–70 μm. For the global equivalent ratio of 1, the air area before the detonation wave front is larger than the vapor area of kerosene droplets, resulting in inhomogeneous mixing before the wave front. Both oil-rich and oil-poor areas form before the wave front. Due to the speed difference between two phases of the gas and droplets, the air is separated to form a low-temperature strip. When the initial diameter of kerosene droplets is small, the mixing process of reactants is mainly affected by evaporation and the detonation wave propagates stably. When the initial droplet diameter is reduced to 1 μm, evaporation occurs at the entrance, and the rotating detonation flow field shows the characteristics of gas phase propagation, and the structure of the detonation wave is smooth. When the initial diameter of kerosene droplets is relatively large, the mixing process of reactants before the wave front is mainly affected by droplet break-up. For the same fuel mass flow rate with different initial droplet diameters, the maximum residence time of kerosene droplets accounts for more than 80% of the detonation wave propagation time and the detonation velocity increases with the increased ratio of gaseous part of the fuel. The velocity of the detonation wave increases first and then decreases with the increased initial droplet diameter in the range of 10–70 μm.
- rotating detonation engine /
- gas-liquid two-phase flow /
- combustion /
- droplet diameters

HTML全文

爆轰发动机作为传统爆燃推进系统的替代方案，在理论上可以实现发动机整体性能的提高^[1-3]。旋转爆轰发动机（rotating detonation engine, RDE）通常使用预爆轰管切向点火，产生一个或多个旋转爆轰波（rotating detonation wave, RDW），爆轰波在接近燃烧室入口处周向传播消耗燃料，产生连续推力^[4-8]。对于气相旋转爆轰发动机，已有学者开展了针对其旋转爆轰波的模态跃迁^[9-10]、推力性能测量^[11-13]、光学测量^[14-16]及热测量^[17-18]等的实验研究；或通过数值模拟方法研究了几何尺寸^[19]，几何结构（如环形^[20-22]、空心^[23]、圆盘形^[24]），燃料和氧化剂喷注方式^[25-26]，不同排气喷管^[27]等对RDE流场的影响。而实际上，发动机多使用液态燃料，相比于气相爆轰，气液两相爆轰需要考虑燃料的雾化、液滴的破碎、燃料部分蒸发及燃料与氧化剂不完全混合等过程。液滴燃料的大小、液滴在时间和空间上的分布等参数对气液两相爆轰过程起到了主要作用^[28]。

液态煤油作为航空动力的常用燃料，在室温条件下很难与空气发生爆轰。由于燃料液滴必须在形成开始放热所需的自由基之前蒸发，因此部分能量用于蒸发液滴，导致用于驱动爆轰反应的热量较少^[29]。实验中通常采用掺氢^[30]、富氧^[31]、煤油部分汽化^[32]等方式来实现旋转爆轰波的传播。由于仅依赖实验难以深入探究RDE内部流场结构，而对于气液两相旋转爆轰等复杂过程的数值模拟，通常也需要对实验条件进行一定的简化处理。已有研究人员在数值模拟中考虑了液滴的蒸发、雾化过程，对初始液滴采用均匀直径假设^[33]和球状假设^[34]。煤油本身作为由链烷烃、芳香烃、环烷烃组成的混合物，数值模拟中通常将其视为单一组分平均分子式或代替燃油组分，然后再给定相应的反应机理。

以往的文献中，涉及RDE流场中液滴分布的分析较少，气液两相RDE的流场研究有待进一步开展。此外，大多数描述液滴破碎的模型理论都是基于无限范围的均匀气流中的单个液滴，研究表明，喷雾液滴的破碎时间可能是孤立液滴的55%～90%，这取决于液滴间距、直径和流量^[35]。泰勒类比破碎（Taylor analogue breakup, TAB）模型由内燃机雾化燃烧过程发展而来，除了应用于低速流动和雾化燃烧^[36-38]外，在涉及激波液滴相互作用的超声速流动中也有广泛应用^[39-40]。由于气液两相介质的爆轰效果取决于液滴的破碎程度和液滴的尺寸分布，而液滴破碎的程度决定了冲击波后化学能的释放速度^[41]。因此，液滴的大小必然对气液两相介质的爆轰参数产生很大的影响。

本文中，以液态煤油C₁₂H₂₃为燃料，高温空气为氧化剂，考虑液滴的雾化破碎、蒸发过程，忽略复杂的燃料裂解过程，采用单步总包反应机理研究旋转爆轰流场，采用成熟的泰勒类比破碎（TAB）模型模拟液滴的雾化破碎过程，以期揭示气液两相RDE的基本流场结构，分析液滴运动，并比较不同初始液滴直径对RDE流场的影响。

1. 研究方案

1.1 数学模型

旋转爆轰流场区域存在气相和液相两相流动。Fluent软件中模拟多相流的模型可以分为欧拉-欧拉多相流模型和欧拉-拉格朗日多相流模型，前者的连续相和离散相均采用欧拉法进行求解，后者的离散相则采用拉格朗日法计算。其中，欧拉多相流模型包括VOF（volume of fluid）模型、混合（mixture）模型和欧拉（Euler）模型。对于离散相体积率超过10%的液滴流动，可选用欧拉多相流模型中的混合模型或者欧拉模型；对于体积分数小于10%的液滴流动，可选用离散相模型（discrete phase model, DPM）。

1.1.1 气相控制方程

气相满足连续方程、动量方程、能量守恒方程和组分方程^[42]：

$\frac{{\partial \rho }}{{\partial t}} + \nabla \cdot \left( {\rho {\boldsymbol{U}}} \right) = \varGamma$

(1)

$\frac{{\partial (\rho u)}}{{\partial t}} + \nabla \cdot \left( {\rho u{\boldsymbol{U}}} \right) = \nabla \cdot \left( {\mu \nabla u} \right) - \frac{{\partial p}}{{\partial x}} + {S _u}$

(2)

$\frac{{\partial \left( {\rho v} \right)}}{{\partial t}} + \nabla \cdot \left( {\rho v{\boldsymbol{U}}} \right) = \nabla \cdot \left( {\mu \nabla v} \right) - \frac{{\partial p}}{{\partial y}} + {S _v}$

(3)

$\frac{{\partial \left( {\rho e} \right)}}{{\partial t}} + \nabla \cdot \left( {\rho e{\boldsymbol{U}}} \right) + p\nabla \cdot {\boldsymbol{U}} = \nabla \cdot \left( {\lambda \nabla T} \right) + \varPhi + {S _{{\rm{in}}}}$

(4)

$\varPhi = 2\mu \left[ {{{\left( {\frac{{\partial u}}{{\partial x}}} \right)}^2} + {{\left( {\frac{{\partial v}}{{\partial y}}} \right)}^2}+ \frac{1}{2}{{\left( {\frac{{\partial u}}{{\partial y}}+ \frac{{\partial v}}{{\partial x}}} \right)}^2}} \right] - \frac{2}{3}\mu {\left( {\nabla {\boldsymbol{U}}} \right)^2}$

(5)

$\frac{{\partial \left( {\rho {\omega_i}} \right)}}{{\partial t}} + \nabla \cdot \left( {\rho {\omega_i}{\boldsymbol{U}}} \right) = \nabla \cdot \left( {{\varGamma _i}\nabla {\omega_i}} \right) + {R_i}$

(6)

式中：ρ为流体密度，U为流体的速度矢量，u和v为速度矢量U的2个分量，Γ为质量源项，μ为动力黏度，p为压力，S_u和S_v为动量方程的广义源项，e为比内能，λ为流体的导热系数，Φ为黏性耗散函数，S_in为内热源，ω_i为组分i的质量分数，Γ_i为组分i的质量扩散系数，R_i为单位时间单位体积因化学反应和液滴蒸发产生的质量。

1.1.2 离散相控制方程

在拉格朗日坐标系下，由牛顿运动定律得出液滴运动方程^[43]：

$\frac{{{\rm{d}}{{\boldsymbol{u}}_{\rm{p}}}}}{{{\rm{d}}t}} = {F_{\rm{D}}}\left( {{\boldsymbol{U}} - {{\boldsymbol{u}}_{\rm{p}}}} \right) + \frac{{{\boldsymbol{g}}\left( {{\rho _{\rm{p}}} - \rho } \right)}}{{{\rho _{\rm{p}}}}} + {\boldsymbol{F}}$

(7)

${F_{\rm{D}}} = \frac{{18\mu }}{{{\rho _{\rm{p}}}{d_{\rm{p}}^2}}}\frac{{{c_{\rm{d}}}Re}}{{24}}$

(8)

$Re = \frac{{\rho {d_{\rm{p}}}}}{\mu }\left| {{{\boldsymbol{u}}_{\rm{p}}} - {\boldsymbol{U}}} \right|$

(9)

式中：u_p为液滴的速度，ρ_p为液滴密度，d_p为液滴直径，Re为相对雷诺数，c_d为阻力系数，g为重力加速度，F为液滴其他作用力。忽略重力。当液滴周围流体的密度远小于液滴密度时，可不考虑虚拟质量力。本文中，液态煤油密度为780 kg/m³，空气密度为1.225 kg/m³，空气密度远小于液态煤油密度，因此不考虑虚拟质量力。

通常，在计算液滴运动阻力时采用球形假设，液滴阻力系数（ $c_{\rm{d,sphere}}$ ）的计算公式^[43]可以表示为：

${c_{{\text{d,sphere}}}} = \left\{ {\begin{array}{*{20}{l}} 0.424 & \;\; \;\; Re {\text{＞}} 1\;000 \\ \dfrac{{24}}{{Re}}\left( {1 + \dfrac{1}{6}R{e^{{2}/{3}}}} \right) & \;\;\;\; Re {\text{≤}} 1\;000 \end{array}} \right.$

(10)

考虑到RDE流动中球形液滴在高速气流作用力下将会产生变形，液滴形状的变化决定了需要动态计算阻力系数，表示为^[43]：

${c_{\rm{d}}} = {c_{{\rm{d,sphere}}}}\left( {1 + 2.632\,y} \right)$

(11)

式中：y为液滴变形因子。

基于液滴震荡、扭曲与弹簧质量系统的比拟，求解控制液滴震荡变形的模型方程组^[44]：

$m\frac{{{{\rm{d}}^2}x}}{{{\rm{d}}{t^2}}} = {F_{\rm{g}}} - kx - d\frac{{{\rm{d}}x}}{{{\rm{d}}t}}$

(12)

$\frac {F_{\rm{g}}}{m} = {c_{\rm{F}}}\frac{\rho_{\rm{g}} }{{{\rho _{\rm{l}}}r}}{{{v_{\rm{r}}^2}}}$

(13)

$\frac{k}{m} = {c_{\rm{k}}}\frac{\sigma }{{{\rho _{\rm{l}}}{r^3}}}$

(14)

$\frac{d}{m} = {c_{\rm{d}}}\frac{{{\mu _{\rm{l}}}}}{{{\rho _{\rm{l}}}{r^2}}}$

(15)

$y = \frac{x}{{{c_{\rm{b}}}r}}$

(16)

式中：m为液滴的质量，x为液滴平衡位置径向变化量，F_g为气动力，ρ_g和ρ_l为气相和液相的密度，d为阻尼系数，k为弹性系数，v_r为液滴与气体的相对速度大小，r为液滴半径，σ和μ_l分别为液滴表面张力系数和黏性系数^[42]，无量纲常数c_F=1/3，c_k=8，c_d=5，c_b=1/2。

通过求解模型方程组，可以确定液滴的震荡变形，液滴震荡增长到临界值时，母液滴将会分裂出许多更小的子液滴。联合式(12)～(16)，可以得出：

$\frac{{{{\rm{d}}^2}y}}{{{\rm{d}}{t^2}}} = \frac{{{c_{\rm{F}}}}}{{{c_{\rm{b}}}}}\frac{\rho_{\rm{g}} }{{{\rho _{\rm{l}}}}}\frac{{{u^2}}}{{{r^2}}} - \frac{{{c_{\rm{k}}}\sigma }}{{{\rho _{\rm{l}}}{r^3}}}y - \frac{{{c_{\rm{d}}}{\mu _{\rm{l}}}}}{{{\rho _{\rm{l}}}{r^2}}}\frac{{{\rm{d}}y}}{{{\rm{d}}t}}$

(17)

当y=0时，将按照球形公式计算液滴的阻力；y=1时，对应液滴临界变形量；y>1时，液滴破碎，式(17)由TAB模型得出。

液滴蒸发过程中，质量守恒方程^[44]为：

$\frac{{{\rm{d}}m}}{{{\rm{d}}t}} = - 4{\text{π}} {r^2}H$

(18)

式中：H为液滴单位面积的蒸发质量速率。

1.2 化学反应模型

化学反应模型为一步反应的层流有限速率模型：

$\rm {C}_{12}{H}_{23}+17.75\left({O}_{2}+3.76{N}_{2}\right)\longrightarrow 11.5{H}_{2}O+12C{O}_{2}+66.74{N}_{2}$

(19)

该模型忽略湍流脉动对化学反应的影响。采用Arrhenius公式计算反应速率^[45]：

$\frac{{{\rm{d}}{c_{{{\rm{C}}_{12}}{{\rm{H}}_{23}}}}}}{{{\rm{d}}t}} = - A{T^b}\exp\left[ { - E/(RT)} \right]{c_{{{\rm{C}}_{12}}{{\rm{H}}_{23}}}^\alpha} {c_{{{\rm{O}}_2}}^\beta}$

(20)

式中：A为指前因子，E为活化能，R为通用气体常数，T为温度， ${c_{{{\rm{C}}_{12}}{{\rm{H}}_{23}}}}$ 为C₁₂H₂₃的物质的量浓度， ${c_{{{\rm{O}}_2}}}$ 为O₂的物质的量浓度，α、β为相应的级数，计算参数如表1所示。

表 1 反应速率计算参数^[46]

Table 1. Parameters used to calculate the reaction rate^[46]

A	b	α	E/(kJ·mol⁻¹)	β	R/(J·mol⁻¹·K⁻¹)
2.587×10⁹	0	0.25	125.6	0.15	8.314

下载: 导出CSV

| 显示表格

1.3 物理模型建立及边界条件

将环形RDE结构展开为二维以降低计算成本，工作原理如图1所示。液态燃料爆轰燃烧反应流场涉及到燃料破碎、剥离、蒸发、扩散、化学反应、相变、湍流两相流、激波等复杂物理化学过程。求解过程参数多、计算量大、收敛难度大，网格质量的好坏直接影响计算收敛。本文中，采用均匀四边形网格来保证网格质量。模型X方向长度为0.11 m，Y方向长度为0.10 m。虚线框设置为3 000 K、2 MPa高温高压点火区，点火区初始X方向速度为2 000 m/s。左、右边界在爆轰波传播第一周时设置为壁面，抑制反向爆轰波的产生。当爆轰波传播至邻近右边界时，将左右边界修改为周期性边界。压力出口边界的背压为0.1 MPa，总温为300 K。空气质量流入口的总温为1 000 K。空气入口采用离散相表面注入方式注入煤油液滴。不同工况的煤油参数如表2所示。

图 1 旋转爆轰发动机工作原理

Figure 1. Operating principle of a RDE

下载: 全尺寸图片幻灯片

表 2 煤油液滴注入参数

Table 2. Injection parameters of kerosene droplets

工况	全局当量比	煤油质量流量/（kg·s⁻¹）	初始液滴直径/μm	液滴注入速度/（m·s⁻¹）	液滴温度/K
1	1	2.056 2	1	50	300
2	1	2.056 2	10	50	300
3	1	2.056 2	20	50	300
4	1	2.056 2	30	50	300
5	1	2.056 2	40	50	300
6	1	2.056 2	50	50	300
7	1	2.056 2	70	50	300

下载: 导出CSV

| 显示表格

1.4 模型与网格有效性验证

根据离散相与空气的质量流量和密度估算出离散相的体积分数为0.011%，远小于10%，因此采用DPM模型。初始液滴粒径为30 μm时RDE流场的离散相体积分数等值线分布如图2所示。图中离散相的体积分数最大为0.05%，验证了模型的合理性。

图 2 初始直径为30 μm的DPM体积分数等值线分布

Figure 2. Contours of DPM volume fraction at an initial diameter of 30 μm

下载: 全尺寸图片幻灯片

真实的喷雾计算中可能有百万液滴，众多液滴轨迹计算成本高昂。两相RDE模拟中引入了粒子包裹的概念，粒子包裹由单个液滴参数表征，用以减小计算量。而网格尺寸影响单次注入的包裹数目，单次注入包裹的数量等于注入边界的网格单元数量。单次注入的每个包裹所包含的液滴数量为：

$N = {\dot m_{\rm{s}}}\frac{{\Delta t}}{{nm}}$

(21)

式中：N为单个包裹的液滴数量， ${\dot m_{\rm{s}}}$ 为液滴流的质量流量，n为注入边界网格单元数量，Δt为液滴时间步。为了保证一定的计算精度，注入的每个包裹所包含的液滴数量不能太多，Δt设置为1 µs。

选取合适的网格尺寸可以合理地控制流场中的包裹数目。对初始直径为30 μm的工况，分别在均匀网格尺寸0.20、0.25、0.40和0.50 mm下进行计算。图3给出了不同网格尺寸下RDE流场稳定后的温度等值线。从图3可以看出，爆轰波、斜激波、缓燃接触面均可捕捉，但在0.20和0.25 mm网格尺寸下，可以更清晰地捕捉到接触间断附近的低温条带。

图 3 不同网格单元尺寸的温度等值线分布

Figure 3. Contours of temperature for different cell sizes

下载: 全尺寸图片幻灯片

在X=0.02 m，Y=0.01 m处设置压力监测点，爆轰波传播稳定后，求解监测点处的爆轰波速度：

$\overline{v} = \frac{{al}}{{\displaystyle\sum\limits_{i = 1}^a {{T_i}} }}$

(22)

式中： $\overline{v}$ 为爆轰波平均速度，l为几何模型的X向长度，T_i为爆轰波压力曲线第i个周期，a为所取的周期个数。不同网格尺寸计算所得的爆轰波平均速度、温度和反应区宽度如表3所示，其中反应区宽度指化学反应速率不为零的区域的宽度。

表 3 不同网格尺寸计算所得的爆轰波平均速度、温度和反应区宽度

Table 3. Average velocity, temperature and reaction zone of detonation waves calculated for different cell sizes

网格单元尺寸/mm	爆轰波平均速度/(m∙s⁻¹)	温度/K	反应区宽度/mm
0.20	1 170	2 509	0.70
0.25	1 160	2 500	0.75
0.40	1 200	2 543	0.85
0.50	1 240	2 482	1.00

下载: 导出CSV

| 显示表格

由表3可以看出，网格尺寸为0.20和0.25 mm时对应的爆轰波参数较接近，结合图3的流场分布，选取0.25 mm网格尺寸进行计算。

2. 流场分析

2.1 气液两相RDE流场结构

以初始液滴直径为50 μm的工况为例，揭示了气液两相旋转爆轰基本流场，如图4所示。图4(a)显示了典型的旋转爆轰波、斜激波与缓燃接触面三波交汇结构，爆轰波运动方向自左向右，红色区域为煤油液滴分布区域，煤油液滴从爆轰波后注入，在爆轰波单个周期内，波前液滴受高温空气影响蒸发时间最长，因此波前存在无液滴区域（图4(a)中白色框区域）。图4(b)中箭头显示了横波的位置与运动方向。此外，温度场中存在低温条带和爆轰波前液滴与低温反应物三角形区域不重合现象。

图 4 初始液滴直径50 μm的温度和压力等值线分布

Figure 4. Contours of temperature and pressure at an initial droplet diameter of 50 μm

下载: 全尺寸图片幻灯片

图5展示了空气Y方向与X方向的速度分布。由图5可知，RDE燃烧室X向速度的最大值集中分布在爆轰波处。轴向Y速度最大值在燃烧室出口处可达千米每秒。由于斜激波的运动导致出口流速不均匀。图5(a)给出了局部放大的爆轰波处Y轴方向速度的等值线分布，入口爆轰波后存在一段空气阻塞区，空气Y轴方向的速度达百米每秒，而液滴的初始注入速度为50 m/s，如图6所示。因此旋转爆轰波后，入射液滴被空气加速，在爆轰波单个旋转周期内，空气的注入高度大于液滴注入高度，导致波前的液滴与空气三角形区域不重合。

图 5 Y方向速度与X方向的速度等值线分布

Figure 5. Contours of Y-velocity and X-velocity

下载: 全尺寸图片幻灯片

图 6 液滴速度等值线分布

Figure 6. Contours of droplet velocity

下载: 全尺寸图片幻灯片

图7为煤油蒸气质量分数与液滴分布。由图7可知，爆轰波前煤油蒸气与液滴高度几乎一致。这是由于爆轰波单个旋转周期较短，而煤油液滴蒸发为气相后未来得及扩散（仅箭头所示，微量煤油蒸气向空气层扩散），而爆轰波的传播需要化学反应提供能量，因此爆轰波强度高的范围也主要集中于波前蒸气区域，约20 mm。

图 7 煤油蒸气质量分数和液滴分布

Figure 7. Mass fraction of kerosene vapor and distribution of droplets

下载: 全尺寸图片幻灯片

图8给出了初始液滴直径为50 µm的O₂、N₂和CO₂的质量分数等值线分布。由图8中O₂、N₂的分布可知，空气三角形最大注入高度约30 mm，且由CO₂分布可知波前反应物在高温空气预蒸发过程中未发生提前燃烧。本文中所有工况的燃料三角形区域均无燃料提前反应现象。由图7推测的爆轰波强度范围低于空气层最大注入高度。因此，爆轰波相对于空气向前运动时，远离燃烧室入口处的空气会分离出去形成低温间断条带，如图4(a)中箭头所示。

图 8 初始液滴直径为50 μm的O₂、N₂和CO₂质量分数等值线分布

Figure 8. Contours of mass fractions of O₂, N₂ and CO₂ at an initial droplet diameter of 50 μm

下载: 全尺寸图片幻灯片

波前的燃料由蒸气和液滴组成，为分析燃料的混合程度，定义当量比^[45]：

$\varphi = \frac{{{{(m_{\rm{air}}/m_{\rm{fule}})}_{{\rm{stoic}}}}}}{{m_{\rm{air}}/m_{\rm{fule}}}} = \frac{{14.59(c {V_{{\rm{cell}}}} + \rho_{\rm{g}} {V_{{\rm{cell}}}} {\omega_{{{\rm{C}}_{{{12}}}}{{\rm{H}}_{{{23}}}}({\rm{g}})}})}}{{\rho_{\rm{g}} {V_{{\rm{cell}}}} ({\omega_{{{\rm{O}}_{{2}}}}} + {\omega_{{{\rm{N}}_2}}})}}$

(23)

式中： $({ m_{\rm{air}}/m_{\rm{fule}}})_{\rm{stoic}}$ 为化学当量的空燃比， $m_{\rm{air}}$ 和 $m_{\rm{fuel}}$ 分别为空气和燃料的质量；c为离散相质量浓度；V_cell为网格单元体积； ${\omega}_{{\mathrm{C}}_{12}{\mathrm{H}}_{23}\left(\mathrm{g}\right)}$ 、 ${\omega}_{{\mathrm{O}}_{2}}$ 、 ${\omega}_{{\mathrm{N}}_{2}}$ 分别为煤油蒸气、O₂和N₂的质量分数。爆轰波后存在未燃烧完全的煤油蒸气。由图8中的O₂和N₂分布可知，氧化剂几乎被爆轰波消耗，式(23)仅分析波前的混合程度。初始液滴直径为50 μm的爆轰波前当量比等值线分布如图9所示。本文中，各工况的注入条件保持全局当量比为1，空气区域大于液滴煤油蒸气区域，导致波前燃料空气混合不均匀，波前均存在富油区及贫油区。

图 9 初始液滴直径50 μm的爆轰波前当量比等值线分布

Figure 9. Contours of equivalence ratios before the detonation wave at an initial droplet diameter of 50 μm

下载: 全尺寸图片幻灯片

2.2 液滴运动分析

2.1节仅展示了液滴的空间分布。为了更清晰地说明气液两相RDE燃烧室内液滴的运动规律，绘制了煤油液滴分布示意图，如图10所示。从图10可以看出，爆轰波后由高压产生了一段空气阻塞区域，由于离散相采用表面注入方式，爆轰波后液滴仍然可以注入，然后与气相相互作用，红色球表示初始注入的、粒径较大的液滴。爆轰波后区域受膨胀波的影响，压力降低，空气开始注入，且液滴被加速，蓝色箭头所指的注入区域为相对低速空气注入区域。爆轰波阻塞区后部分液滴并未进入产物区域，且与气相相对速度较低，此时液滴不易破碎，缓慢蒸发向下游运动，最终形成较大液滴条带，如图10中绿色球形液滴所示。红色注入区域为相对高速空气注入区域，此时液滴在与空气较大的剪切作用下发生剧烈破碎，产生较小粒径的紫色球形液滴。流场粒径大小分布如图11所示。

图 10 液滴分布示意图

Figure 10. Schematic diagram of droplets distribution

下载: 全尺寸图片幻灯片

图 11 初始液滴直径50 μm液滴直径等值线分布

Figure 11. Contours of the droplet diameter at an initial droplet diameter of 50 μm

下载: 全尺寸图片幻灯片

在爆轰波传播的单个循环过程中，液滴由爆轰波后注入，爆轰波前远离燃烧室入口的液滴经历的蒸发时间最长，图12显示了初始液滴直径50 μm工况下流场中煤油液滴的停留时间等值线分布，最大液滴停留时间为86.9 µs。

图 12 初始液滴直径50 μm的液滴停留时间等值线分布

Figure 12. Contours of the droplet residence time at an initial droplet diameter of 50 μm

下载: 全尺寸图片幻灯片

2.3 不同初始液滴直径对流场的影响

2.3.1 初始直径1 μm工况

当初始直径为1 μm时，液滴在入口附近即被高温空气完全蒸发。根据监测点（X=0.02 m，Y=0.01 m）压力曲线计算爆轰波平均速度达1 330 m/s。气液两相旋转爆轰表现出气相爆轰的特性，爆轰波极为平整。图13为初始液滴直径为1 μm的温度和压力等值线分布。

图 13 初始液滴直径1 μm的温度和压力等值线分布

Figure 13. Contours of temperature and pressure at an initial droplet diameter of 1 μm

下载: 全尺寸图片幻灯片

2.3.2 初始直径10～70 μm工况

图14为初始直径为10～70 μm范围工况下旋转爆轰波自持稳定传播后液滴直径等值线分布。

图 14 不同初始液滴直径的液滴直径等值线分布

Figure 14. Droplet diameter distribution for different initial droplet diameters

下载: 全尺寸图片幻灯片

2.2节中分析了初始直径为50 μm工况下旋转爆轰流场的液滴直径分布规律，解释了爆轰波前远离入口处较大液滴条带产生的原因。由图14(a)可知，初始液滴直径为10 μm时，流场中存在大量未破碎的液滴。初始液滴直径增大至20 μm时，流场中仅混杂少量未破碎的液滴，如图14(b)中的箭头所示。初始液滴直径为30 μm时，流场中几乎不存在未破碎的液滴。在爆轰波的一个旋转周期内，爆轰波后初始液滴注入，此时主要经历蒸发混合等过程，受空气剪切破碎的影响较小，蒸发速度较慢。

从较小初始液滴直径范围的工况可知，初始液滴直径越小，液滴越具有保持稳定的趋势，越不容易破碎。相应地，对于较大初始液滴直径，如图14(e)和(f)，初始液滴直径越大，液滴越具有不稳定的趋势，在空气作用下越容易发生破碎。这导致随着初始直径的增大，下游的大液滴条带不明显。

定义流场中平均液滴直径^[47]：

$\overline {D} = \frac{{\displaystyle\int_0^\infty {Df(D){\rm{d}}D} }}{{\displaystyle\int_0^\infty {f(D){\rm{d}}D} }}$

(24)

式中： $f(D)$ 为流场中液滴直径的概率密度分布函数，等式右边分母值为1。

图15为液滴平均直径随初始液滴直径变化的函数分布。可以看出，除10 μm工况外，随着初始直径的增大，流场的液滴平均直径先减小后缓慢增大。40～70 μm的初始液滴直径工况范围内，流场平均液滴直径比10～30 μm工况更小且差别不大。这再次验证了较大初始直径的液滴更易破碎。可以推断出对于较大的初始液滴直径，不同的初始液滴直径受相同的高速空气流剪切破碎后，总是会破碎成差别不大的液滴直径。

图 15 液滴平均直径随初始液滴直径变化的函数分布

Figure 15. Droplet mean diameter as a function of the initial droplet diameter

下载: 全尺寸图片幻灯片

图16展示了在X=0.02 m、Y=0.01 m监测点处计算所得的爆轰波平均速度、压力和温度。可以看出，随着初始直径增大，爆轰波的压力波动变化，而爆轰波的温度变化趋势则与爆轰波速度变化趋势基本一致。

图 16 爆轰波速度、温度和压力随初始液滴直径变化的分布函数

Figure 16. Detonation wave velocity, temperature and pressure as functions of the initial droplet diameter

下载: 全尺寸图片幻灯片

表4给出了不同液滴直径工况流场中液滴最大停留时间和爆轰波周期。在不同工况下，液滴最大停留时间均占爆轰波传播时间尺度的80%以上，保证了爆轰波的自持传播。

表 4 不同初始直径的液滴最大停留时间和爆轰波周期

Table 4. Maximum residence time of droplets and detonation cycle time for different droplet diameters

初始液滴直径/μm	最大停留时间/μs	爆轰波周期/μs	最大停留时间与爆轰波周期的比值/%
10	80.0	96.5	82.9
20	79.0	91.7	86.2
30	85.0	94.8	89.7
40	86.0	94.8	90.7
50	86.9	96.5	90.1
70	84.5	96.5	87.6

下载: 导出CSV

| 显示表格

为了分析爆轰波速度发生变化的原因，定义爆轰波波前不同初始直径工况下液滴的蒸发效率：

$\eta = \frac{{{\rho _{\text{g}}}{V_{{\text{cell}}}}{f_{{{\text{C}}_{12}}{{\text{H}}_{23}}\left( {\text{g}} \right)}}}}{{{\rho _{\text{g}}}{V_{{\text{cell}}}}{f_{{{\text{C}}_{12}}{{\text{H}}_{23}}\left( g \right)}} + c{V_{{\text{Cell}}}}}}$

(25)

式(25)仅计算爆轰波前的蒸发混合过程，爆轰波后无液滴，η恒为1。

初始直径30 μm工况的蒸发效率等值线分布如图17所示。初始直径工况下爆轰波前蒸发效率最高。取工况波前区域求解蒸发效率平均值。工况爆轰波前白色矩形区域划分标准为：Y方向自燃烧室入口沿Y正向延伸，高度为爆轰波高度；X方向宽度均为0.01 m，X坐标自爆轰波的最高点处沿X正向延伸。不同初始直径工况下平均蒸发效率与爆轰波速度的关系如图18所示。

图 17 初始液滴直径为30 μm工况的蒸发效率分布

Figure 17. Distribution of evaporation efficiency for an initial droplet diameter of 30 μm

下载: 全尺寸图片幻灯片

图 18 平均蒸发效率、爆轰波速度与初始液滴直径的关系

Figure 18. Average evaporation efficiency and detonation wave velocity as functions of the initial droplet diameter

下载: 全尺寸图片幻灯片

从图18可以看出，爆轰波波前的蒸发效率与爆轰波速度随初始直径的变化趋势一致。燃料质量流量不变的情况下，改变初始直径，波前蒸发效率较高时，爆轰波前气相燃料与所有燃料的质量比值较大，爆轰波速度较高。

3. 结　论

采用单步总包反应机理，考虑液滴的雾化破碎、碰撞黏合过程，采用初始均匀直径假设，在直径1～70 μm工况范围内计算了二维RDE的气液两相流流场，得出以下主要结论。

(1) 不同初始液滴直径工况下，旋转爆轰燃烧室内均形成了稳定传播的爆轰波，由于波前液滴与煤油蒸气分布不均匀，爆轰波呈现曲面分布。仅当初始直径减小至1 μm时，爆轰波才变得更平整，呈现出气相爆轰特性。

(2) 初始液滴越小，越具有保持稳定的趋势，在爆轰波周期内受高温空气作用，此时主要经历蒸发过程。而初始液滴越大，越易破碎成小液滴，此时受空气破碎影响较大。

(3) 对于相同的燃料质量流量，在不同初始液滴直径工况下，爆轰波前燃料的蒸发效率越高，爆轰波速度越大。初始直径为10～70 μm工况下，随着初始液滴直径的增大，爆轰波速度先升高后降低，爆轰波速度极大值出现在初始液滴直径为20 μm的工况。

图 1 旋转爆轰发动机工作原理

Figure 1. Operating principle of a RDE

下载: 全尺寸图片幻灯片

图 2 初始直径为30 μm的DPM体积分数等值线分布

Figure 2. Contours of DPM volume fraction at an initial diameter of 30 μm

下载: 全尺寸图片幻灯片

图 3 不同网格单元尺寸的温度等值线分布

Figure 3. Contours of temperature for different cell sizes

下载: 全尺寸图片幻灯片

图 4 初始液滴直径50 μm的温度和压力等值线分布

Figure 4. Contours of temperature and pressure at an initial droplet diameter of 50 μm

下载: 全尺寸图片幻灯片

图 5 Y方向速度与X方向的速度等值线分布

Figure 5. Contours of Y-velocity and X-velocity

下载: 全尺寸图片幻灯片

图 6 液滴速度等值线分布

Figure 6. Contours of droplet velocity

下载: 全尺寸图片幻灯片

图 7 煤油蒸气质量分数和液滴分布

Figure 7. Mass fraction of kerosene vapor and distribution of droplets

下载: 全尺寸图片幻灯片

图 8 初始液滴直径为50 μm的O₂、N₂和CO₂质量分数等值线分布

Figure 8. Contours of mass fractions of O₂, N₂ and CO₂ at an initial droplet diameter of 50 μm

下载: 全尺寸图片幻灯片

图 9 初始液滴直径50 μm的爆轰波前当量比等值线分布

Figure 9. Contours of equivalence ratios before the detonation wave at an initial droplet diameter of 50 μm

下载: 全尺寸图片幻灯片

图 10 液滴分布示意图

Figure 10. Schematic diagram of droplets distribution

下载: 全尺寸图片幻灯片

图 11 初始液滴直径50 μm液滴直径等值线分布

Figure 11. Contours of the droplet diameter at an initial droplet diameter of 50 μm

下载: 全尺寸图片幻灯片

图 12 初始液滴直径50 μm的液滴停留时间等值线分布

Figure 12. Contours of the droplet residence time at an initial droplet diameter of 50 μm

下载: 全尺寸图片幻灯片

图 13 初始液滴直径1 μm的温度和压力等值线分布

Figure 13. Contours of temperature and pressure at an initial droplet diameter of 1 μm

下载: 全尺寸图片幻灯片

图 14 不同初始液滴直径的液滴直径等值线分布

Figure 14. Droplet diameter distribution for different initial droplet diameters

下载: 全尺寸图片幻灯片

图 15 液滴平均直径随初始液滴直径变化的函数分布

Figure 15. Droplet mean diameter as a function of the initial droplet diameter

下载: 全尺寸图片幻灯片

图 16 爆轰波速度、温度和压力随初始液滴直径变化的分布函数

Figure 16. Detonation wave velocity, temperature and pressure as functions of the initial droplet diameter

下载: 全尺寸图片幻灯片

图 17 初始液滴直径为30 μm工况的蒸发效率分布

Figure 17. Distribution of evaporation efficiency for an initial droplet diameter of 30 μm

下载: 全尺寸图片幻灯片

图 18 平均蒸发效率、爆轰波速度与初始液滴直径的关系

Figure 18. Average evaporation efficiency and detonation wave velocity as functions of the initial droplet diameter

下载: 全尺寸图片幻灯片

表 1 反应速率计算参数^[46]

Table 1. Parameters used to calculate the reaction rate^[46]

A	b	α	E/(kJ·mol⁻¹)	β	R/(J·mol⁻¹·K⁻¹)
2.587×10⁹	0	0.25	125.6	0.15	8.314

下载: 导出CSV

表 2 煤油液滴注入参数

Table 2. Injection parameters of kerosene droplets

工况	全局当量比	煤油质量流量/（kg·s⁻¹）	初始液滴直径/μm	液滴注入速度/（m·s⁻¹）	液滴温度/K
1	1	2.056 2	1	50	300
2	1	2.056 2	10	50	300
3	1	2.056 2	20	50	300
4	1	2.056 2	30	50	300
5	1	2.056 2	40	50	300
6	1	2.056 2	50	50	300
7	1	2.056 2	70	50	300

下载: 导出CSV

表 3 不同网格尺寸计算所得的爆轰波平均速度、温度和反应区宽度

Table 3. Average velocity, temperature and reaction zone of detonation waves calculated for different cell sizes

网格单元尺寸/mm	爆轰波平均速度/(m∙s⁻¹)	温度/K	反应区宽度/mm
0.20	1 170	2 509	0.70
0.25	1 160	2 500	0.75
0.40	1 200	2 543	0.85
0.50	1 240	2 482	1.00

下载: 导出CSV

表 4 不同初始直径的液滴最大停留时间和爆轰波周期

Table 4. Maximum residence time of droplets and detonation cycle time for different droplet diameters

初始液滴直径/μm	最大停留时间/μs	爆轰波周期/μs	最大停留时间与爆轰波周期的比值/%
10	80.0	96.5	82.9
20	79.0	91.7	86.2
30	85.0	94.8	89.7
40	86.0	94.8	90.7
50	86.9	96.5	90.1
70	84.5	96.5	87.6

下载: 导出CSV

参考文献(47)

[1]	TEASLEY T, WILLIAMS B, LARKEY A, et al. A review towards the design optimization of high-performance additively manufactured rotating detonation rocket engine injectors [C] // AIAA Propulsion and Energy Forum. Virtual Event: AIAA, 2021. DOI: 10.2514/6.2021-3655.
[2]	王顺利, 吴云, 金迪, 等. 不同当量比下喷管对旋转爆震特性的影响研究 [J]. 爆炸与冲击, 2020, 40(10): 102102. DOI: 10.11883/bzycj-2019-0481. WANG S L, WU Y, JIN D, et al. Effects of nozzles on performance of rotating detonation at different equivalence ratios [J]. Explosion and Shock Waves, 2020, 40(10): 102102. DOI: 10.11883/bzycj-2019-0481.
[3]	RAJ P, MEADOWS J. A method for modeling dynamic injector response in a rotating detonation combustor (RDC) [C] // AIAA Propulsion and Energy Forum. Virtual Event: AIAA, 2021. DOI: 10.2514/6.2021-3647.
[4]	葛高杨, 马元, 侯世卓, 等. 当量比对汽油燃料两相旋转爆轰发动机工作特性影响实验研究 [J]. 爆炸与冲击, 2021, 41(11): 112102. DOI: 10.11883/bzycj-2020-0465. GE G Y, MA Y, HOU S Z, et al. Experimental study on the effect of equivalent ratio on working characteristics of gasoline fuel two-phase rotating detonation engine [J]. Explosion and Shock Waves, 2021, 41(11): 112102. DOI: 10.11883/bzycj-2020-0465.
[5]	李宝星, 许桂阳, 翁春生, 等. 燃烧室宽度对液态燃料旋转爆轰发动机影响实验研究 [J]. 推进技术, 2021, 42(2): 372–381. DOI: 10.13675/j.cnki.tjjs.190472. LI B X, XU G Y, WENG C S, et al. Experimental investigation for effects of chamber width on rotating detonation engine with liquid fuel [J]. Journal of Propulsion Technology, 2021, 42(2): 372–381. DOI: 10.13675/j.cnki.tjjs.190472.
[6]	PAXSON D E, FOTIA M L, HOKE J, et al. Comparison of numerically simulated and experimentally measured performance of a rotating detonation engine [C] // 53rd AIAA Aerospace Sciences Meeting. Kissmmee, Florida, USA: AIAA, 2015. DOI: 10.2514/6.2015-1101.
[7]	ZHOU R, WU D, WANG J P. Progress of continuously rotating detonation engines [J]. Chinese Journal of Aeronautics, 2016, 29(1): 15–29. DOI: 10.1016/j.cja.2015.12.006.
[8]	STECHMANN D, LIM D, HEISTER S D. Survey of rotating detonation wave combustor technology and potential rocket vehicle applications [C] // 50th AIAA/ASME/SAE/ASEE Joint Propulsion Conference. Cleveland, OH, USA: AIAA, 2014. DOI: 10.2514/6.2014-3902.
[9]	WANG Y H, QIAO W Y, LE J L. Combustion characteristics in rotating detonation engines [J]. International Journal of Aerospace Engineering, 2021: 8839967. DOI: 10.1155/2021/8839967.
[10]	SUCHOCKI J A, YU S J, HOKE J L, et al. Rotating detonation engine operation [C] // 50th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition. Nashville, Tennessee, USA: AIAA, 2012. DOI: 10.2514/6.2012-119.
[11]	FOTIA M L, KAEMMING T A, HOKE J, et al. Study of the experimental performance of a rotating detonation engine with nozzled exhaust flow [C] // 53rd AIAA Aerospace Sciences Meeting. Kissimmee, Florida, USA: AIAA, 2015. DOI: 10.2514/6.2015-0631.
[12]	ISHIHARA K, KATO Y, MATSUOKA K, et al. Performance evaluation of a rotating detonation engine with conical-shape tail [C] // 53rd AIAA Aerospace Sciences Meeting. Kissimmee, Florida, USA: AIAA, 2015. DOI: 10.2514/6.2015-0630.
[13]	FOTIA M, SCHAUER F, HOKE J. Experimental study of performance scaling in rotating detonation engines operated on hydrogen and gaseous hydrocarbon fuel [C] // 20th AIAA International Space Planes and Hypersonic Systems and Technologies Conference. Glasgow, Scotland: AIAA, 2015. DOI: 10.2514/6.2015-3626.
[14]	WANG Y H, LE J L. Rotating detonation engines with two fuel orifice schemes [J]. Acta Astronautica, 2019, 161: 262–275. DOI: 10.1016/j.actaastro.2019.05.035.
[15]	NAKAGAMI S, MATSUOKA K, KASAHARA J, et al. Visualization of rotating detonation waves in a plane combustor with a cylindrical wall injector [C] // 53rd AIAA Aerospace Sciences Meeting. Kissimmee, Florida, USA: AIAA, 2015. DOI: 10.2514/6.2015-0878.
[16]	RANKIN B A, RICHARDSON D R, CASWELL A W, et al. Imaging of OH* chemiluminescence in an optically accessible nonpremixed rotating detonation engine [C] // 53rd AIAA Aerospace Sciences Meeting. Kissimmee, Florida, USA: AIAA, 2015. DOI: 10.2514/6.2015-1604.
[17]	MEYER S J, POLANKA M D, SCHAUER F R, et al. Parameter impact on heat flux in a rotating detonation engine [C] // 53rd AIAA Aerospace Sciences Meeting. Kissimmee, Florida, USA: AIAA, 2018. DOI: 10.2514/6.2018-0400.
[18]	THEUERKAUF S W, SCHAUER F R, ANTHONY R, et al. Experimental characterization of high-frequency heat flux in a rotating detonation engine [C] // 53rd AIAA Aerospace Sciences Meeting. Kissimmee, Florida, USA: AIAA, 2015. DOI: 10.2514/6.2015-1603.
[19]	李宝星, 许桂阳, 舒慧明, 等. 燃烧室轴向和周向长度对气液两相旋转爆轰特性的影响 [J]. 航空动力学报, 2020, 35(8): 1601–1611. DOI: 10.13224/j.cnki.jasp.2020.08.005. LI B X, XU G Y, SHU H M, et al. Influence of axial and circumferential lengths of combustion chamber on gas-liquid two-phase rotating detonation characteristics [J]. Journal of Aerospace Power, 2020, 35(8): 1601–1611. DOI: 10.13224/j.cnki.jasp.2020.08.005.
[20]	徐灿, 邓利, 马虎, 等. 环缝宽度对旋转爆震发动机工作特性的影响 [J]. 爆炸与冲击, 2019, 39(3): 032102. DOI: 10.11883/bzycj-2017-0248. XU C, DENG L, MA H, et al. Annular gaps width effecting on performance of rotating detonation engine [J]. Explosion and Shock Waves, 2019, 39(3): 032102. DOI: 10.11883/bzycj-2017-0248.
[21]	ZHANG L F, ZHANG S J, WANG J P. A direct simulation of continuous detonation engine with the Navier-Stokes equations [C] //21st AIAA International Space Planes and Hypersonics Technologies Conference. Xiamen, Fujian, China: AIAA, 2017. DOI: 10.2514/6.2017-2323.
[22]	SHAO Y T, LIU M, WANG J P. Continuous detonation engine and effects of different types of nozzle on its propulsion performance [J]. Chinese Journal Aeronautics, 2010, 23(6): 647–652. DOI: 10.1016/S1000-9361(09)60266-1.
[23]	SUN J, ZHOU J, LIU S J, et al. Numerical investigation of a non-premixed hollow rotating detonation engine [J]. International Journal of Hydrogen Energy, 2019, 44(31): 17084–17094. DOI: 10.1016/j.ijhydene.2019.04.168.
[24]	HAYASHI K A, SHIMOMURA K, TSUBOI N, et al. 3D numerical study on flow field in disc-RDE [C] // AIAA Propulsion and Energy Forum. Virtual Event: AIAA, 2021. DOI: 10.25Virtual Event14/6.2021-3665.
[25]	徐雪阳, 卓长飞, 武晓松, 等. 非预混喷注对旋转爆震发动机影响的数值研究 [J]. 航空学报, 2016, 37(4): 1184–1195. DOI: 10.7527/s1000-6893.2015.0195. XU X Y, ZHUO C F, WU X S, et al. Numerical simulation of injection schemes with separate supply of fuel and oxidizer effects on rotating detonation engine [J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(4): 1184–1195. DOI: 10.7527/s1000-6893.2015.0195.
[26]	周蕊, 李晓鹏. 连续旋转爆轰发动机冷流场的混合特性研究 [J]. 航空学报, 2016, 37(12): 3668–3674. DOI: 10.7527/s1000-6893.2016.0108. ZHOU R, LI X P. Numerical investigation of mixing characteristic of cold continuously rotating detonation engine [J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(12): 3668–3674. DOI: 10.7527/s1000-6893.2016.0108.
[27]	SARACOGLU, BAYINDIR H, et al. Unsteady performance of rotating detonation engines with different exhaust nozzles [J]. Journal of Propulsion and Power, 2017, 33(1): 1–10. DOI: 10.2514/1.B36164.
[28]	KAILASANATH K. Liquid-fueled detonations in tubes [J]. Journal of Propulsion and Power, 2006, 22(6): 1261–1268. DOI: 10.2514/1.19624.
[29]	SCHWER D A. Multi-dimensional simulations of liquid-fueled JP-10/oxygen detonations [C] // AIAA Propulsion and Energy Forum. Indianapolis, USA: AIAA, 2019. DOI: 10.2514/6.2019-4042.
[30]	BYKOVSKII F A, ZHDAN S A, VEDERNIKOV E F. Continuous spin detonation of a heterogeneous kerosene-air mixture with addition of hydrogen [J]. Combustion, Explosion, and Shock Waves, 2016, 52(3): 371–373. DOI: 10.1134/S0010508216030187.
[31]	ZHENG Q, MENG H L, WENG C S, et al. Experimental research on the instability propagation characteristics of liquid kerosene rotating detonation wave [J]. Defence Technology, 2020, 16(6): 1106–1115. DOI: 10.1016/j.dt.2020.06.028.
[32]	LI J M, CHANG P H, LI L, et al. Investigation of injection strategy for liquid-fuel rotating detonation engine [C] // AIAA Aerospace Sciences Meeting. Kissimmee, Florida, USA: AIAA, 2018. DOI: 10.2514/6.2018-0403.
[33]	MENG Q Y, ZHAO M J, ZHENG H T, et al. Eulerian-Lagrangian modelling of rotating detonative combustion in partially pre-vaporized n-heptane sprays with hydrogen addition [J]. Fuel, 2021, 290: 119808. DOI: 10.1016/j.fuel.2020.119808.
[34]	徐高, 翁春生, 康楠, 等. 考虑燃料雾化的气液两相连续旋转爆轰数值模拟 [J]. 推进技术, 2022, 43(1): 200249. DOI: 10.13675/j.cnki.tjjs.200249. XU G, WENG C S, KANG N, et al. Numerical simulation of gas-liquid two-phase continuous rotating detonation considering fuel atomization [J]. Journal of Propulsion Technology, 2022, 43(1): 200249. DOI: 10.13675/j.cnki.tjjs.200249.
[35]	FOX G E, DABORA E K. Breakup of liquid drops due to convective flow in shocked sprays [J]. Symposium (International) on Combustion, 1973, 14(1): 1365–1373. DOI: 10.1016/S0082-0784(73)80122-5.
[36]	KHOSLA S, CROCKER D S. CFD modeling of the atomization of plain liquid jets in cross flow for gas turbine applications [C] // Proceedings of ASME Turbo Expo 2004: Power for Land, Sea, and Air. Vienna, Austria: American Society of Mechanical Engineers, 2004. DOI: 10.1115/GT2004-54269.
[37]	KHARE P, WANG S, YANG V. Modeling of finite-size droplets and particles in multiphase flows [J]. Chinese Journal of Aeronautics, 2015, 28(4): 974–982. DOI: 10.1016/j.cja.2015.05.004.
[38]	BROUMAND M, FAROKHI M, BIROUK M. Penetration height of a circular liquid jet in a subsonic gaseous crossflow: an eulerian-lagrangian approach [C] // 54th AIAA Aerospace Sciences Meeting. San Diego, California, USA: AIAA, 2016. DOI: 10.2514/6.2016-1591.
[39]	BHANDARKAR A, MANNA P, CHAKRABORTY D. Assessment of droplet breakup models in high-speed cross-flow [J]. Atomization and Sprays, 2017, 27(1): 61–79. DOI: 10.1615/atomizspr.2016015409.
[40]	HU R, LI Q, LI C, LI C. Effects of an accompanied gas jet on transverse liquid injection in a supersonic crossflow [J]. Acta Astronautica, 2019, 159: 440–451. DOI: 10.1016/j.actaastro.2019.01.040.
[41]	BURCAT A, EIDELMAN A. Evolution of a detonation wave in a cloud of fuel droplets: part Ⅱ. influence of fuel droplets [J]. AIAA Journal, 1980, 18(10): 80–4099. DOI: 10.2514/3.7717.
[42]	吴望一. 流体力学 [M]. 北京: 北京大学出版社, 2001.
[43]	LIU A B, MATHER D, REITZ R D. Modeling the effects of drop drag and breakup on fuel sprays: 930072 [R]. SAE Transactions, Detroit, USA, 1993.
[44]	O’ROURKE P J, AMSDEN A A. The TAB method for numerical calculation of spray droplet breakup [R]. Los Alamos, NM, USA: Los Alamos National Laboratory, 1987.
[45]	STEPHEN R T, DANIEL C H. An introduction to combustion: concepts and applications [M]. New York, USA: McGraw Hill, 2000.
[46]	GUO K K, NIE W S, CAI H H. Numerical simulation of damping capacity between injector-formed baffle and normal blade baffle in a kero/lox liquid rocket engine [C] // Proceedings of the Advances in Materials, Machinery, Electical Engineering (AMMEE 2017). Prais, France: Atlantis Press, 2017: 576–585. DOI: 10.2991/ammee-17.2017.110.
[47]	KUO K K Y. Principles of combustion [M]. New York, USA: John Wiley and Sons, 1986.

施引文献

期刊类型引用(3)

1.	徐鸿飞，王放，武郁文，翁春生. 基于小尺寸模型分析氢氧旋转爆轰波传播的不稳定性机制. 爆炸与冲击. 2025(01): 30-42 . 本站查看
2.	杨志龙，田文斌，张珍，王志英，王一伟. 单孔排气气泡特征识别与参数提取方法. 化工进展. 2024(02): 808-817 . 百度学术
3.	刘秋月，王放，翁春生，赵庆军. 预混煤油/空气两相旋转爆轰传播特性数值研究. 推进技术. 2024(09): 118-128 . 百度学术

其他类型引用(1)

资源附件(0)

访问统计

图(18) / 表(4)

计量

文章访问数: 661
HTML全文浏览量: 213
PDF下载量: 78
被引次数: 4

1. 研究方案
1.1 数学模型
1.2 化学反应模型
1.3 物理模型建立及边界条件
1.4 模型与网格有效性验证
2. 流场分析
2.1 气液两相RDE流场结构
2.2 液滴运动分析
2.3 不同初始液滴直径对流场的影响
3. 结　论

煤油液滴直径对两相旋转爆轰发动机流场的影响

doi: 10.11883/bzycj-2022-0068

作者简介: 杨 帆（1996－ ），男，硕士研究生，yfhbxt@163.com

通讯作者: 王宇辉（1986－ ），男，副教授，aowuki@163.com

计量

出版历程

Influence of kerosene droplet diameters on the flow field of a two-phase rotating detonation engine

1. 研究方案

1.1 数学模型

1.1.1 气相控制方程

1.1.2 离散相控制方程

1.2 化学反应模型

1.3 物理模型建立及边界条件

1.4 模型与网格有效性验证

2. 流场分析

2.1 气液两相RDE流场结构

2.2 液滴运动分析

2.3 不同初始液滴直径对流场的影响

2.3.1 初始直径1 μm工况

2.3.2 初始直径10～70 μm工况

3. 结 论

期刊类型引用(3)

其他类型引用(1)

计量

出版历程

目录

1. 研究方案

1.1 数学模型

1.2 化学反应模型

1.3 物理模型建立及边界条件

1.4 模型与网格有效性验证

2. 流场分析

2.1 气液两相RDE流场结构

2.2 液滴运动分析

2.3 不同初始液滴直径对流场的影响

3. 结 论

作者简介:
杨　帆（1996－　），男，硕士研究生，yfhbxt@163.com

通讯作者:
王宇辉（1986－　），男，副教授，aowuki@163.com

3. 结　论

3. 结　论