鱼雷燃烧室内衬压力平衡孔数值模拟研究

陈 肯，冯要飞，宗 潇，伊进宝

(中国船舶集团有限公司 第七〇五研究所，陕西 西安 710077)

0 引 言

鱼雷燃烧室是鱼雷热动力装置中最重要的部件之一[1]，HAP 三组元推进剂在燃烧室内燃烧，产生高温高压的燃气作为工质推动发动机做功[2]。为保证燃烧室在工作过程中稳定可靠，燃烧室壳体需进行合理设计[3]。现有结构下的燃烧室，在工作后其内壳体上有明显的烧蚀环槽。

经分析是因公差燃烧室内壳体与内衬在接触位置处存在缝隙，燃气经缝隙泄漏造成冲刷现象，形成烧蚀环槽。为解决这一问题，结构改进的一种方案为在内衬上增开一对压力平衡孔，以解决壳体的烧蚀现象。而压力平衡孔的位置及大小将直接影响结构改进方案的效果。

热流固耦合数值仿真的方法广泛应用于鱼雷、航空等各领域的燃烧室的仿真计算工作当中[4-14]。本文针对具有不同压力平衡孔位置及大小的鱼雷燃烧室结构为研究对象，利用Fluent 软件开展燃烧室高工况工作条件下的热流固耦合数值仿真计算，分析燃烧室内衬压力平衡孔的位置及大小对燃气域流场及燃烧室内壳体应力场的影响规律。

1 燃烧室工作过程中物化过程的数学模型

鱼雷燃烧室工作时，涉及到HAP 三组元推进剂破碎雾化、液滴颗粒及内流场气相的两相流动、液滴颗粒的蒸发燃烧、燃烧室壳体及冷却水与燃气之间的耦合换热这4 个过程。

1.1 液滴破碎控制方程

液体射流在空气中做相对运动，空气动力促使射流破碎；液体表面张力使射流不破碎。在2 个力的作用下液膜破碎成小液滴，如图1 所示。

图1 液滴雾化破碎机理Fig. 1 Droplet atomization and breakage mechanism

由射流扰动机理，对三组元推进剂的破碎状态定义破碎系数R，其数学计算式为：

式中： λ为射流扰动波长；d0为射流起始直径； σ为液体表面张力；d∗为破碎后的液滴的标称直径。

取R-R分布为液滴雾化后的直径分布依据，按中间质量分布准则作为直径的分布准则，可知液滴直径d的分布表达式为：

式中，i为液滴直径分布指数。

1.2 液滴运动及气相两相流动控制方程

破碎后的液滴，在运动过程受到燃气粘性的摩擦阻力，无体积力。离散相燃料液滴的运动控制方程为：

式中：u为流体相速度；up为 颗粒速度； ρ为流体密度； ρp为颗粒密度；d为颗粒直径；Re为颗粒雷诺数；CD为阻力系数。

对于气相部分，连续相流体运动控制方程为：

1）连续性方程

式中：Sm为质量源相，数值上等于注入的燃料在单位时间蒸发的质量流量。

2）动量方程

广义上连续相流体运动的控制方程为[9]：

式中：P为静压； µ 为气相动力粘性；F为外部体积力；SM为动量源相；∇ 为拉普拉斯算子。

对于燃烧室气相，因重力因素忽略，无外部体积力，故气相简化后的动量方程为：

3）湍流模型

燃烧室燃烧过程中燃气处于高湍流状态，选用k-ε湍流模型，其控制方程为：

式中： µ为粘性系数； µt为湍流粘性系数；Gk为速度梯度湍动能；Gb为浮力湍动能，对不可压流体该值为0；YM为湍流脉 动湍动能；C为修正系数常量； σk，σε为湍流普朗特数。

1.3 液滴蒸发燃烧控制方程

1）液滴蒸发控制方程

液滴未达到蒸发温度的温度控制方程[11]为：

达到蒸发温度后的液滴质量控制方程[11]为：

式中：mp，cp，Tp， εp，Ap分别为颗粒的质量、比热、温度、发射率和表面积；h为对流传热系数；T∞为气相温度； σ为波尔兹曼常数； θR为 辐射温度；hfg为潜热。

2）液滴燃烧控制方程

HAP 三组元推进剂完全燃烧的化学反应方程式为：

反应释放的能量由反应物与生成物的总焓差决定。

Eddy-dissipation 涡耗散模型的反应速率Ri,k控制方程为[11]：

式中：mp为所有产物组分和某一特定反应物的质量比重；Mw,j为 物质的分子 量；分别为反应物和生成物的化学当量系数；A和B为经验常数，分别取4 和0.5。

1.4 耦合换热过程能量控制方程

流体微元体微分形式的能量方程为：

2 数值计算几何模型及网格划分

2.1 几何模型预处理及域提取

对燃烧室的几何结构进行简化，以降低网格划分难度、降低网格数量、提高网格质量、帮助模型更快更好的收敛。

利用提取功能处理简化完成的燃烧室固体域，得到研究对象的全体计算域，如图2 所示。

图2 计算域提取Fig. 2 Picking up the compute area

保温层域的缝隙厚度由燃烧室内壳体与内衬的最大加工公差决定，缝隙厚度为0.166 mm，上下两处缝隙位置的几何结构局部放大图如图3 所示。

图3 缝隙处结构Fig. 3 Structure in the gap

2.2 计算域网格划分

采用混合网格的划分方式完成数值计算模型的网格划分，其中使用ICEM 软件的block 功能进行结构化网格的划分，使用ICEM 软件中的八叉树算法进行非结构化网格的划分。在Fluent 中使用make polyhedral 功能将四面体网格转为多面体网格，最终划分的网格数量在千万级左右，网格划分结果如图4 所示。

图4 网格划分结果Fig. 4 The meshing result

2.3 网格无关性验证

验证网格尺寸不会对数值计算结果造成显著影响，开展网格无关性验证。进行网格无关性验证的网格划分方案如表1 所示。

本文针对燃烧室燃气域的流场进行研究，故选取燃烧室燃气域出口段的温度变化情况作为衡量网格无关性的指标，网格无关性验证的结果如图5 所示。

结合表1 的网格划分方案及图5 的网格计算结果，可知网格划分方案2 为最优的网格划分方式，可作为之后数值计算中的网格划分方案。

图5 网格无关性验证结果Fig. 5 Result about grid independence verification

表1 网格无关性验证网格划分方案Tab. 1 Mesh parameters about grid independence verification

3 边界条件及压力平衡孔几何参数设置

3.1 边界条件设置

针对燃烧室热流固耦合数值仿真计算，以燃烧室在高工况下的燃料、冷却水、内压等工况参数作为数值仿真计算的边界参数。具体的边界条件参数设置如表2 所示。

表2 数值计算边界条件设置Tab. 2 The boundary condition setting parameters

3.2 压力平衡孔几何参数

设置9 组不同压力平衡孔位置及尺寸的数值仿真计算case，用于比较压力平衡孔位置及大小对燃烧室燃气域的流场及壳体的应力场的影响，具体结构参数设置如表3 所示。

表3 压力平衡孔几何参数设置Tab. 3 The pressure balance hole geometry parameters setting

4 数值计算结果及分析

4.1 压力平衡孔位置对燃烧室流场及壳体应力场的影响

1）上方缝隙处

不同压力平衡孔位置下，燃气域上方缝隙处的压力、速度、热通量变化图如图6~图8 所示。

由图6 可知，压力平衡孔与内衬上表面距离越小，缝隙两端的压差越小。经分析原因在于燃气腔内的燃气越接近内衬上表面，燃气的压力值与内衬上表面处的压力值越接近，开通压力平衡孔后压力联通，使得缝隙两端的压差减小。

图6 不同压力平衡孔位置燃气域上方缝隙压力图Fig. 6 The pressure diagram in the gas upper gap under the different pressure balance hole location

由图7 可知，压力平衡孔位置的改变不对缝隙处燃气的速度变化规律造成影响，但对速度的大小会产生显著影响，且压力平衡孔与内衬上表面越近，缝隙内的燃气流速越低，这是由于缝隙两端的压差减小导致的。平衡孔位置的改变对缝隙内燃气流速的最大影响值约为3 m/s。

图7 不同压力平衡孔位置燃气域上方缝隙速度图Fig. 7 The velocity diagram in the gas upper gap under the different pressure balance hole location

缝隙内燃气流速的降低，会导致燃气对流换热量的减小，所以压力平衡孔与内衬上表面的距离越小，缝隙处的热通量大小越小，图8 可以印证这一分析。

图8 不同压力平衡孔位置燃气域上方缝隙热通量图Fig. 8 The heat flux diagram in the gas upper gap under the different pressure balance hole location

2）下方缝隙处

不同压力平衡孔位置下，燃气域下方缝隙处的压力、速度、热通量变化图如图9~图11 所示。

图9 不同压力平衡孔位置燃气域下方缝隙压力图Fig. 9 The pressure diagram in the gas bottom gap under the different pressure balance hole location

图11 不同压力平衡孔位置燃气域下方缝隙热通量图Fig. 11 The heat flux diagram in the gas bottom gap under the different pressure balance hole location

可知，改变压力平衡孔的位置对下方缝隙处的燃气的压强、流速、缝隙处的热通量等参数均不产生较大影响。分析认为原因是燃气在外腔内流动，当流至下方缝隙的上方时燃气的流动状态已经趋于稳定，故压力平衡孔对下方缝隙处的燃气流动已不会造成显著影响。

3）壳体应力场

图12 为不同压力平衡孔位置燃烧室壳体的应力场云图。可知，随着压力平衡孔与内衬上表面的距离的减小，壳体缝隙处的最大应力值减小。结合图6~图8，可知原因在于缝隙两端压差降低，降低了燃气流速从而减少了缝隙处的热通量。降低了壳体在缝隙处的温度，从而减小了热应力的大小。

图12 不同压力平衡孔位置壳体应力场图Fig. 12 The stress diagram in the inside shell under the different pressure balance hole location

4.2 压力平衡孔大小对燃烧室流场及壳体应力场的影响

1）上方缝隙处

不同压力平衡孔直径下，燃气域上方缝隙处的压力、速度、热通量变化图如图13~图15 所示。

图10 不同压力平衡孔位置燃气域下方缝隙速度图Fig. 10 The velocity diagram in the gas bottom gap under the different pressure balance hole location

由图13 可知，随着压力平衡孔的直径增大，缝隙两端的压差先变小后增大。分析认为原因在于随着孔径的增大，燃气域外腔内的燃气压力会越来越接近内腔，使得缝隙两端压差减小，但孔径过大会使燃气在孔的末端发展出湍流，在湍流中心处形成低压区，反而使得缝隙两端压差增大。

图13 不同压力平衡孔直径燃气域上方缝隙压力图Fig. 13 The pressure diagram in the gas upper gap under the different pressure balance hole diameter

由图14 可知，压力平衡孔直径的改变不对缝隙处燃气的速度变化规律造成影响，但对速度的大小会产生显著影响。随着压力平衡孔直径的增大，缝隙内的燃气流速先降低后增高，这是由于缝隙两端的压差先变小后变大所导致的。且在孔径大小为4 mm 时，缝隙内的燃气流速最低为8.5 m/s。随着孔径的改变，缝隙内燃气流速的最大差值约为2 m/s。

图14 不同压力平衡孔直径燃气域上方缝隙处速度图Fig. 14 The velocity diagram in the gas upper gap under the different pressure balance hole diameter

由图15 可知，上方缝隙处燃气的对流换热量与燃气的流速具有相同的变化规律，即先变小后变大。

图15 不同压力平衡孔直径燃气域上方缝隙热通量图Fig. 15 The heat flux diagram in the gas upper gap under the different pressure balance hole diameter

2）下方缝隙处

不同压力平衡孔直径下，燃气域下方缝隙处的压力、速度、热通量变化如图16~图18 所示。

图16 不同压力平衡孔直径燃气域下方缝隙两端压力图Fig. 16 The pressure diagram in the gas bottom gap under the different pressure balance hole diameter

图18 不同压力平衡孔直径燃气域下方缝隙热通量图Fig. 18 The heat flux diagram in the gas bottom gap under the different pressure balance hole diameter

可以看出，改变压力平衡孔的直径对燃烧域下方缝隙处的缝隙两端压强、缝隙内的燃气流速、缝隙处的热通量等参数均不产生较大影响，经分析认为原因在于燃气在外腔内流动，当流至下方缝隙上方时燃气的流动状态已经趋于稳定，故压力平衡孔对下方缝隙处的燃气流动不会造成影响。

3）壳体应力场

图19 为不同压力平衡孔直径燃烧室壳体的应力场云图。可知，随着压力平衡孔直径的增大，壳体缝隙处的最大应力值先减小后增大，且当孔径为4 mm时，壳体所受的最大应力值最小。结合图13~图15，可知原因在于缝隙两端压差降低，降低了燃气流速从而减少了缝隙处的热通量。降低了壳体在缝隙处的温度，从而减小了热应力的大小。

图19 不同压力平衡孔直径壳体应力场图Fig. 19 The stress diagram in the inside shell under the different pressure balance hole diameter

图17 不同压力平衡孔直径燃气域下方缝隙处速度图Fig. 17 The velocity diagram in the gas bottom gap under the different pressure balance hole diameter

5 结 语

本文针对压力平衡孔的位置及大小对鱼雷燃烧室的燃气域流场及壳体域应力场的影响规律进行数值计算，经分析可以得出如下结论：

1）对于燃气域上方缝隙处，压力平衡孔的位置距离内衬上表面越近，燃气域缝隙两端的压差越小、燃气的流速越低、缝隙处的热通量越小；对于下方缝隙处，压力平衡孔的位置变化对燃气域无明显影响；对于壳体的最大应力值，压力平衡孔距离内衬上表面越近，最大应力值越小。

2）对于燃气域上方缝隙处，随着压力平衡孔的直径变大，燃气域缝隙两端的压差先变小后变大、燃气的流速先降低后增高、缝隙处的热通量先减小后增大；对于下方缝隙处，压力平衡孔的直径变化对燃气域无明显影响；对于壳体的最大应力值，随着压力平衡孔直径的增大，最大应力值先减小后增大，且当直径为4 mm 时最大应力值最小。