弹库导弹注水保护技术关键参数计算方法

李士军，杨 巍，潘树国，汪长波

(中国船舶集团有限公司第七一三研究所，河南郑州450015)

0 前 言

舰船弹库是舰船上储存弹药的舱室，不同弹库储存的弹药类型不同，舰载导弹平时储存在弹库内。舰载导弹自身带有固体火箭发动机，内部携带有大量的固体推进剂。导弹中的固体推进剂处于一种相对稳定的状态，当受到机械强冲击、爆炸冲击波、高温、静电等袭扰时，会发生意外点火，火箭发动机会在短时间内向外喷射高温、高压、高速燃气射流；燃气射流温度超过1000℃，速度超过2000m/s，喷射出燃起流从十几千克到几十千克不等。舰船弹库是一个密闭空间，导弹在短时间内意外点火后，造成周围弹药处于高温高压下，极易发生殉爆，同时会使得弹库内部空间的压力急剧升高，造成弹库发生爆炸。李士军等[1–2]先后研究了导弹意外点火时弹库泄压排气过程、喷洒抑制剂的影响。国外舰船弹库中一般均配置注水系统，例如“鞑靼人”舰空导弹库、M K 41导弹库、MK 48导弹库等。注水系统是通过注水喷头向意外点火的导弹发动机燃气射流中注入高压水，以便降低火箭发动机燃气射流温度，避免导弹发动机在弹库中意外点火时的高温燃气对周围弹药、舰船设备和船体造成损伤，如图1所示。我国国军标GJB 4 000-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）液滴中的温度相同，不存在温度梯度，处于饱和状态即Tsat=393K；

6）气体和水的热力学特性与温度无关；

7）水滴的拖拽力和传热相关，通过公式表示；

8）控制体的轴向长度取值为发动机尾喷管直径的2倍，即n=L/dj=2，确保控制体取在发动机射流中的核心区，发动机射流核心区的长度采用经验公式表示，式中Mje为发动机尾喷管气流的Ma数；

9）液滴为刚性球体，不会受气流作用发生变形；

10）忽略辐射热传导作用和由于水的冲击碰撞、破碎消耗的能量。

2 控制方程

根据图2，建立控制界面1和界面2之间的连续方程、动量方程、能量方程和理想气体状态方程：

式中：m˙W为水沿轴向的质量流率；Fd为液滴拖拽力；up为 液滴速度；η 为蒸发的水质量份数；cpl为液体的比定压热容；m˙W为进入燃气流中的水总质量流率；m˙Wt为水喷嘴喷出的水总质量流率，变量下标1和2代表控制体两侧侧面。

图2 注水控制体示意图Fig.2 Diagram of control body

进入燃气流中的水质量流率和喷水总的质量流率之间的为m˙W=m˙WiN，N为喷头数量。

3 控制体两侧的参数比

将质量守恒方程式（1）代入动量守恒方程式（2），通过变换，可以得到控制体两侧射流速度之比：

式（5）中右侧括号中的第二项可以写作：

式中，CD为拖拽系数，拖拽力Fd表达式为：

将式（5）代入式（4）可以改写为：

式（3）中右边括号内第三项可以改写为雷诺数形式表示：

将式（6）和式（8）代入式（3）中可得：

由能量守恒式（3）和质量守恒式（1），可得到控制体进出口界面上的温度比值：

根据气体等熵方程，可得：

式中，Mj1为 超声速气体的当地马赫数R为气体常数。

将式（6）、式（8）、式（11）和式（12）代入式（10），可得：

4 方程中的参数计算

式（5）中的拖拽系数与拖拽力成正比，直接影响拖拽力的大小。CD是 与Rep数有关的函数，研究表明[8]，当Re数很小(R ep≤0.1)时，粘性力占主导地位，惯性力可以忽略不计，此时的流动为蠕动流，称为Stokes流动，拖拽系数当Rep数增大时，惯性力变得越来越重要而不能被忽略，根据不同的Rep数的使用范围，许多学者提出了相关的经验公式，本项目采用White[7]提出的采用Spalding拖拽系数修正的公式为水的汽化潜热，Tsat为水的饱和温度。

气体传给液滴群的热量为：

qp=η˙mw hfg

注水后液体的汽化部分的吸热量 ，应和式（15）相等，并考虑雷诺数，可得：

式（16）说明在控制体内，水的汽化蒸发比例与发动机喷出的气体质量流率和注水的质量流率无关。

5 控制体内的不变量

6 算例分析

式中： µ0=1.711×10−5N/m2，C=122K，经计算可得µj1=5.444×10−5N/m2。

水的密度为 ρp=1 000 kg/|m3，注水压力在0.6MPa，喷头连接管路的水流速度为5m/s，喷头通径为35 mm，单个喷头的出水量为4.3 kg/s，假设安装2个喷头，喷入水的速度为3 1.6 m/s。 水的汽化潜热hfg=2.256 6×106J/kg ，定压比热容为Cpl=4 220 J/kg，标准大气压下水的饱和温度为Tsat=373K，水的热传导系 数 为k=0.58 W/(m·K)， µpl=1.01×10−3Pa·s，水的温度为Tp=293 K（20℃）。

η=73.08%

在蒸发量 的条件下，喷入水的质量流率和燃气流质量流率比值对燃气温度的影响如图3所示。

图3 注水量对温度的影响曲线Fig.3 Influence curve of water injection on temperature

7 结语

本文采用参数比拟法，充分考虑了水射流形成的液滴及液滴汽化影响，提出注水量与燃气流场温度的计算方法，并通过算例进行了计算分析。该方法可为注水系统的设计提供理论支持和计算方法。