舰船弹库是舰船上储存弹药的舱室,不同弹库储存的弹药类型不同,舰载导弹平时储存在弹库内。舰载导弹自身带有固体火箭发动机,内部携带有大量的固体推进剂。导弹中的固体推进剂处于一种相对稳定的状态,当受到机械强冲击、爆炸冲击波、高温、静电等袭扰时,会发生意外点火,火箭发动机会在短时间内向外喷射高温、高压、高速燃气射流;燃气射流温度超过1000 ℃,速度超过2000 m/s,喷射出燃起流从十几千克到几十千克不等。舰船弹库是一个密闭空间,导弹在短时间内意外点火后,造成周围弹药处于高温高压下,极易发生殉爆,同时会使得弹库内部空间的压力急剧升高,造成弹库发生爆炸。李士军等[1-2]先后研究了导弹意外点火时弹库泄压排气过程、喷洒抑制剂的影响。国外舰船弹库中一般均配置注水系统,例如“鞑靼人”舰空导弹库、MK41导弹库、MK48导弹库等。注水系统是通过注水喷头向意外点火的导弹发动机燃气射流中注入高压水,以便降低火箭发动机燃气射流温度,避免导弹发动机在弹库中意外点火时的高温燃气对周围弹药、舰船设备和船体造成损伤,如图1所示。我国国军标GJB 4000-2000第7组武器发射装置和保障系统721.5.3.4条弹库安全性中规定:发射装置的弹库应根据需要选用喷淋降温、强排风、注水、灭火措施,但是没有给出设计方法和设计指标。
![]() |
图 1 弹库注水示意图 Fig. 1 Diagram of water injection in ammunition depot |
国内外研究向火箭发动机注水的主要目的是为了降低火箭发动机噪声,张彬乾等[3]通过大量试验研究了亚、超音速射流噪声场的分布规律及影响因素,分析了射流噪声源产生的机理,为后面研究射流噪声抑制技术打下了基础;徐悦等[4]分析了火箭燃气射流噪声的气动特性,利用喷水试验探讨了喷水对火箭发动机点火时燃气的喷水降噪效果以及影响喷水降噪的因素;马宏伟等[5]对单兵火箭燃气射流噪声进行了实验研究,获得了燃气射流噪声近场的分布规律,同时应用小波变换方法处理燃气射流噪声数据。张磊等[6]对单兵火箭燃气射流噪声抑制进行了试验研究。Thomas[7]通过试验验证了喷水能有效地抑制燃气射流噪声,并讨论了影响喷水降噪的关键因素;Kandula[8-9]对注水抑制噪声进行了试验研究。
1 注水问题假设条件根据具体问题进行物理假设:
1)发动机尾喷管出口处的压力为常量;
2)射流和射流混合物气体符合理想气体定量,采用理想气体状态方程,气体常数近似一致,且液滴均匀分布在控制体内的气体混合物中;
3)液滴在控制体内的分布均匀;
4)射流中的液滴直径相同;
5)液滴中的温度相同,不存在温度梯度,处于饱和状态即
6)气体和水的热力学特性与温度无关;
7)水滴的拖拽力和传热相关,通过公式表示;
8)控制体的轴向长度取值为发动机尾喷管直径的2倍,即
9)液滴为刚性球体,不会受气流作用发生变形;
10)忽略辐射热传导作用和由于水的冲击碰撞、破碎消耗的能量。
2 控制方程根据图2,建立控制界面1和界面2之间的连续方程、动量方程、能量方程和理想气体状态方程:
${\rho _{j1}}{u_{j1}}{A_{j1}} = {\rho _{j2}}{u_{j2}}{A_{j2}} - \eta {\dot m_W}\text{,}$ | (1) |
${\rho _{j2}}u_{j2}^2{A_{j2}} - {\rho _{j1}}u_{j1}^2{A_{j1}} = - {F_d} + \eta {\dot m_W}{u_p}\text{,}$ | (2) |
$\begin{split}&{\dot m_{j2}}\left( {{c_{pj}}{T_{j2}}+\frac{{u_{j2}^2}}{2}} \right) =\\ &{\dot m_{j1}}\left( {{c_{pj}}{T_{j1}}+\frac{{u_{j1}^2}}{2}} \right) - {F_d}{u_p}+\eta {\dot m_W}{T_p}{c_{pl}}\text{,}\end{split}$ | (3) |
${\rho _{j1}}{R_{j1}}{T_{j1}} = {\rho _{j2}}{R_{j2}}{T_{j2}}\text{。}$ | (4) |
式中:
![]() |
图 2 注水控制体示意图 Fig. 2 Diagram of control body |
进入燃气流中的水质量流率和喷水总的质量流率之间的为
将质量守恒方程式(1)代入动量守恒方程式(2),通过变换,可以得到控制体两侧射流速度之比:
$\frac{{{u_{j2}}}}{{{u_{j1}}}} = \frac{1}{{1 + \eta \dfrac{{{{\dot m}_W}}}{{{{\dot m}_{j1}}}}}} \cdot \left( {1 - \frac{{{F_d}}}{{{\rho _{j1}}u_{j1}^2{A_{j1}}}} + \frac{{\eta {{\dot m}_W}{u_p}}}{{{\rho _{j1}}u_{j1}^2{A_{j1}}}}} \right)\text{,}$ | (5) |
式(5)中右侧括号中的第二项可以写作:
$\varphi = \frac{{{F_d}}}{{{\rho _{j1}}u_{j1}^2{A_{j1}}}} = \frac{1}{2}{N_p}{C_D}\frac{{{A_P}}}{{{A_{j1}}}}{\left( {1 - \frac{{{u_p}}}{{{u_{j1}}}}} \right)^2}\text{,}$ | (6) |
式中,
${F_d} = \frac{1}{2}{\rho _{j1}}{\left( {{u_{j1}} - {u_p}} \right)^2}{C_D}{A_{PT}}\text{,}$ | (7) |
${N_p} = \frac{3}{2}\left( {\frac{{\dot m}}{{{{\dot m}_{j1}}}}} \right) \cdot \left( {\frac{{{\rho _{j1}}}}{{{\rho _p}}}} \right) \cdot {\left( {\frac{{{{\rm d}_{j1}}}}{{{{\rm d}_p}}}} \right)^3} \cdot n\left( {1 - \eta } \right)\text{,}$ | (8) |
将式(5)代入式(4)可以改写为:
$\varphi = \psi \cdot \left( {\frac{{{{\dot m}_W}}}{{{{\dot m}_{j1}}}}} \right)\left( {1 - \eta } \right)\text{,}$ | (9) |
式中,
式(3)中右边括号内第三项可以改写为雷诺数形式表示:
$\frac{{\eta {{\dot m}_W}{u_p}}}{{{\rho _{j1}}u_{j1}^2{A_{j1}}}} = \eta \left( {\frac{{{{\dot m}_W}}}{{{{\dot m}_{j1}}}}} \right)\left[ {1 - \frac{{{{\operatorname{Re} }_p}}}{{{{\operatorname{Re} }_{j1}}}} \cdot \frac{{{{\rm d}_{j1}}}}{{{{\rm d}_p}}}} \right]\text{,}$ | (10) |
将式(6)和式(8)代入式(3)中可得:
$\begin{split}\frac{{{u_{j2}}}}{{{u_{j1}}}} = &\frac{1}{{1 + \eta \dfrac{{{{\dot m}_W}}}{{{{\dot m}_{j1}}}}}}\left[ {1 - \psi \cdot \left( {\frac{{{{\dot m}_W}}}{{{{\dot m}_{j1}}}}} \right)\left( {1 - \eta } \right)} \right.+ \\ &\left. { \eta \left( {\frac{{{{\dot m}_W}}}{{{{\dot m}_{j1}}}}} \right)\left( {1 - \frac{{{{\operatorname{Re} }_p}}}{{{{\operatorname{Re} }_{j1}}}} \cdot \frac{{{{\rm d}_{j1}}}}{{{{\rm d}_p}}}} \right)} \right]\text{,}\end{split}$ | (11) |
由能量守恒式(3)和质量守恒式(1),可得到控制体进出口界面上的温度比值:
$\begin{split}\frac{{{T_{j2}}}}{{{T_{j1}}}} =& \frac{1}{{\left( {1 + \eta \dfrac{{{{\dot m}_W}}}{{{{\dot m}_{j1}}}}} \right)}}\left[ {1 + \frac{{u_{j1}^2}}{{2{c_{pj}}{T_{j1}}}} - \frac{{{F_d}{u_p}}}{{{{\dot m}_{j1}}{c_{pj}}{T_{j1}}}}} \right. +\\ &\left. { \frac{{\eta {{\dot m}_W}{T_p}{c_{pl}}}}{{{{\dot m}_{j1}}{c_{pj}}{T_{j1}}}}} \right] - \frac{{u_{j2}^2}}{{2{c_{pj}}{T_{j1}}}}\text{,}\end{split}$ | (12) |
根据气体等熵方程,可得:
$\frac{{u_{j1}^2}}{{2{c_{pj}}{T_{j1}}}} = \frac{{\gamma - 1}}{2}M_{j1}^2\text{,}$ | (13) |
$\frac{{u_{j2}^2}}{{2{c_{pj}}{T_{j1}}}} = {\left( {\frac{{{u_{j2}}}}{{{u_{j1}}}}} \right)^2}\left( {\frac{{\gamma - 1}}{2}} \right)M_{j1}^2\text{,}$ | (14) |
式中,
将式(6)、式(8)、式(11)和式(12)代入式(10),可得:
$\begin{split}\frac{{{T_{j2}}}}{{{T_{j1}}}} =& \frac{1}{{\left( {1 + \eta \dfrac{{{{\dot m}_W}}}{{{{\dot m}_{j1}}}}} \right)}}\left[ {1 + \frac{{\gamma - 1}}{2}M_{j1}^2 - \psi \cdot \left( {\frac{{{{\dot m}_W}}}{{{{\dot m}_{j1}}}}} \right)\left( {1-\eta } \right)\frac{{{u_{j1}}{u_p}}}{{{c_{pj}}{T_{j1}}}}} \right. +\\ &\left. { \frac{{\eta {{\dot m}_W}{T_p}{c_{pl}}}}{{{{\dot m}_{j1}}{c_{pj}}{T_{j1}}}}} \right] - {\left( {\frac{{{u_{j2}}}}{{{u_{j1}}}}} \right)^2}\left( {\frac{{\gamma - 1}}{2}} \right)M_{j1}^2\text{,}\\[-18pt]\end{split}$ | (15) |
由理想气体状态方程的假设,可得
$\frac{{M{a_{j2}}}}{{M{a_{j1}}}} = \frac{{\dfrac{{{u_{j2}}}}{{{c_{j2}}}}}}{{\dfrac{{{u_{j1}}}}{{{c_{j1}}}}}} = \frac{{{u_{j2}}}}{{{u_{j1}}}} \cdot \sqrt {\frac{{\gamma R{T_{j1}}}}{{\gamma R{T_{j2}}}}} = \frac{{{u_{j2}}}}{{{u_{j1}}}} \cdot \sqrt {\frac{{{T_{j1}}}}{{{T_{j2}}}}} \text{。}$ | (16) |
式(5)中的拖拽系数与拖拽力成正比,直接影响拖拽力的大小。
气体传给液滴群的热量为:
${q_g} = \text{π} {N_p}{N_u}k{d_p}\left( {{T_{je}} - {T_p}} \right)\text{,}$ | (17) |
式中,
注水后液体的汽化部分的吸热量
$\frac{1}{\eta } = 1 + \frac{{Pr }}{{{ N}u}} \cdot \frac{{{{{Re} }_{j1}}}}{{{{\left( {\dfrac{{{d_{j1}}}}{{{d_p}}}} \right)}^2}}} \cdot \frac{1}{{\dfrac{3}{2} \cdot \left( {\dfrac{{{\rho _{j1}}}}{{{\rho _p}}}} \right)n}}\left[ {\frac{{{h_{fg}}}}{{{c_{pj}}\left( {{T_{je}} - {T_{sat}}} \right)}}} \right]\text{,}$ | (18) |
式(16)说明在控制体内,水的汽化蒸发比例与发动机喷出的气体质量流率和注水的质量流率无关。
5 控制体内的不变量根据雷诺数、式(8)和式(16),令
$\frac{d_{j 1}}{d_{p}}=\sqrt{\frac{1}{C_{1}} \cdot \frac{\operatorname{Re}_{j 1}}{\rho_{j 1} / \rho_{p}}}\text{,}$ | (19) |
$\operatorname{Re}_{p}=C_{2} \frac{\operatorname{Re}_{j 1}}{{\rm d}_{j 1} / {\rm d}_{p}}\text{。}$ | (20) |
以某导弹发动机参数为例进行计算分析。发动机燃气射流总温
${\mu _{j1}} = {\mu _0}{\left( {\dfrac{{{T_{j1}}}}{{273}}} \right)^{1.5}}\dfrac{{273 + C}}{{{T_{j1}} + C}}\text{。}$ |
式中:
水的密度为
液滴直径在
本算例中的
一般液滴的直径在
$\begin{split}\frac{{{T_{j2}}}}{{{T_{j1}}}} = &\frac{1}{{\left( {1 + \eta \xi } \right)}}{\left[ {2.156 - 2.46 \cdot \xi \left( {1 - \eta } \right) + 0.453\eta \xi } \right]^2} -\\ & 1.156{\left( {\frac{1}{{1 + \eta \xi }}} \right)^2} \cdot \left[ {1 - 1.13\xi \left( {1 - \eta } \right) + 0.95\eta \xi } \right]\text{。}\end{split}$ | (21) |
在蒸发量
![]() |
图 3 注水量对温度的影响曲线 Fig. 3 Influence curve of water injection on temperature |
本文采用参数比拟法,充分考虑了水射流形成的液滴及液滴汽化影响,提出注水量与燃气流场温度的计算方法,并通过算例进行了计算分析。该方法可为注水系统的设计提供理论支持和计算方法。
[1] |
李士军等. 舰载弹药舱泄压排气理论研究[J]. 舰船科学技术, 2015, 37(4): 227-231. DOI:10.3404/j.issn.1672-7649.2015.04.051 |
[2] |
张起, 李士军. 舰载导弹意外点火时舱室内喷洒抑制的影响研究[J]. 舰船科学技术, 2015, 38(8): 194-197. |
[3] |
张彬乾, 张正科. 射流噪声场特性的试验研究[J]. 西北工业大学学报, 1993, 11(4). |
[4] |
徐悦, 周旭, 张志成, 等. 火箭发动机燃气射流喷水降噪研究[J]. 航空动力学报, 2010, 25(4). |
[5] |
刘礼军. 火箭导弹发射噪声实验与抑制技术研究[D]. 南京: 南京理工大学, 2004.
|
[6] |
张磊, 阮文俊, 王浩. 单兵火箭燃气射流噪声抑制的实验研究[J]. 火力与指挥控制, 2015, 40(7). DOI:10.3969/j.issn.1002-0640.2015.07.045 |
[7] |
THOMAS D. Reductions in multi-component jet noise by water injection[R]. AIAA, 2004: 2976.
|
[8] |
KANDULA M, MICHAEL J. Effective jet properties for the estimation of turbulent mixing noise reduction by water injection[R]. AIAA, 2007: 3654.
|
[9] |
KANDULA M. On the scaling laws for jet noise in subsonicand supersonic flow[R]. AIAA, 2003: 3288.
|