过载对固体发动机凝相粒子运动的影响分析

2020-11-10王立武赵朝坤田维平1郭运强

空气动力学学报 2020年5期

王立武, 赵朝坤, 田维平1,, 郭运强

(1. 西北工业大学燃烧、热结构与内流场重点实验室, 西安 710072; 2. 中国航天科技集团有限公司第四研究院第四十一研究所, 西安 710025; 3. 中国航天科技集团有限公司第四研究院, 西安 710025)

0 引言

随着高机动性能战术导弹技术的发展，对高横向过载下固体发动机的工作稳定性和安全性提出了越来越高的要求[1]。过载会引起燃烧室流场中的凝相粒子局部聚集，对发动机内流动和装药燃烧产生很大的影响，并恶化了发动机绝热层的工作环境[2-5]。对于含高铝粉推进剂的发动机来说，高过载引起燃烧室和喷管内部高速粒子流的局部聚集和冲刷，严重时会导致内绝热防护失效[6-8]。

国外20世纪60年代已经发现过载对火箭飞行性能和安全有重要影响，Corwe[9]和King[10]提出了燃烧物理模型和动力学模型。Greatrix[11]等建立了燃烧室和喷管内稳态和非稳态流动方程的数值模型、侵蚀燃烧模型，研究了横向过载对推进剂燃烧过程的影响。Sabnis[12]采用欧拉-拉格朗日方法模拟了过载条件下粒子轨迹在发动机内部的运动轨迹，认为150 μm的粒子受过载的影响强烈，而50 μm的粒子受气流的影响更为强烈。

国内关于过载条件下固体发动机的研究主要集中在燃烧室绝热层烧蚀方面，并取得了一定研究成果[13-15]，在过载条件对喷管烧蚀的影响方面研究则较少,主要为喷管内流场仿真研究[5,16-19]。何国强等[13]提出了固体发动机过载条件下内流场两相流计算方法，针对内孔装药结构燃烧室进行了纵、横向过载条件下两相流动数值模拟，结果表明，在壁面绝热层上存在粒子冲刷密度最大点，且该点随横向过载的增加而远离发动机后封头。李越森等[14]通过试验验证了横向加速度严重影响局部绝热层烧蚀的事实，认为该现象是由于横向加速度作用下燃气中Al2O3液态粒子偏离发动机中心线,沿离心力方向大量沉积所致。许团委等[15]针对战术发动机的飞行条件，开展了3种典型中、小飞行过载下的流场计算，详细分析了燃烧室内凝相颗粒冲刷参数分布，并初步分析了颗粒聚集状态与绝热层烧蚀之间的关系。姜贵庆等[16]计算了中型喷管中的粒子轨迹，认为入口粒子冲击角对粒子轨迹有明显影响。淡林鹏等[17]利用颗粒轨道模型计算得出了长尾喷管中粒子运动的某些规律。刘洋等[18]对某长尾管发动机内的三维两相流动进行了数值模拟，认为横向过载对颗粒运动轨迹和聚集浓度影响显著，纵向过载影响小。张志峰等[19]应用随机轨道模型对长尾喷管发动机长尾段内流场进行了数值模拟，研究3种不同颗粒直径分布下颗粒的运动轨迹和颗粒在壁面聚集浓度的区别。

固体发动机工作过程中，过载条件对燃烧室绝热层烧蚀和喷管收敛段烧蚀均会产生重要影响，但目前暂未见到综合考虑过载对燃烧室和喷管影响的相关报道。为此，本文采用CFD方法[20]，针对某端面燃烧的固体发动机，开展了过载条件下固体发动机内流场仿真研究，首次同时考虑了过载条件对燃烧室和喷管粒子浓度的综合作用，分析了过载条件对凝相粒子运动规律及壁面粒子浓度分布的影响，并阐明了过载条件引起壁面粒子浓度增大的原因。

1 理论模型和边界条件

1.1 理论模型

与常规固体发动机两相流内流场仿真相比，本文发动机内流场仿真过程中引入了与何国强[9]、许团委[16]等相同的过载条件。计算过程中，把气相作为连续相，采用可压缩黏性平均N-S方程的数值解和湍流模型来求解发动机流场；颗粒相作为离散相进行处理，采用颗粒轨道模型；使用迎风格式离散对流项，中心差分格式离散黏性项，一阶后差格式离散时间项。内流场仿真过程中，作出如下假设：1) 气相满足理想气体状态方程，比热比为常数；2) 不考虑气相化学反应；3) 不考虑颗粒相的燃烧、蒸发、破碎过程；4) 不考虑两相间能量的交换；5) 不考虑气相及颗粒相的湍流效应。

图1 网格划分模型Fig.1 Grid of calculation model

图2 两种网格下压强曲线计算结果Fig.2 Calculation results for two kinds of grid

1.2 边界条件

本文计算中涉三类边界条件：出口边界条件，固体壁面边界条件和加质壁面边界条件。出口边界条件为环境压强和温度；固体壁面采用黏附条件，壁面绝热；在加质入口，设定粒子的质量流率。针对加质壁面边界条件，考虑到本模型具体工况，使用文献[21]的结果，采用颗粒直径的Rosin-Rammler分布，分布指数取3.5，粒度分布范围为10～120 μm，粒径峰值为70 μm。此外，根据仿真需求，给定相应的过载条件。

为了分析过载条件对发动机凝相粒子运动的影响规律，针对几种不同组合过载条件进行了内流场仿真计算，详见表1。其中，横向过载表示为xy平面朝y方向的过载，轴向过载为沿x轴方向的过载。

表1 发动机内流场计算条件Table 1 Conditions of SRM field simulation

2 计算结果及分析

2.1 过载对凝相粒子运动的影响

图3～图5给出了算例1#、5#、9#粒子迹线图，图中颜色表示粒子在发动机中的停留时间。对于有过载和无过载情况，粒子撞击到壁面后若未被壁面吸附，将发生反弹，继续随气流流出喷管；随着过载条件的加入，会引起粒子在发动机内部停留时间的增加，算例5#和算例9#由无过载的0.0481 s分别增加到0.17 s和0.0743 s。

图3 无过载粒子迹线图(单位：s)Fig.3 Particle residence time without overload (unit:s)

图4 横向过载40g时的粒子迹线图(单位：s)Fig.4 Particle residence time under 40g lateral overload (unit:s)

图5 横向过载40g、轴向过载40g时的粒子迹线图(单位：s)Fig.5 Particle residence time under 40g lateral and 40g axial overload (unit:s)

分析认为，无过载时粒子与壁面的撞击主要发生在喷管收敛段，纯横向过载作用下粒子与燃烧室壁面和喷管收敛段壁面均发生撞击，粒子与燃烧室壁面撞击后反弹继续随气流运动，粒子与燃烧室壁面撞击改变了粒子的运动轨迹，撞击和反弹过程中粒子的速度较小，增加了粒子在低速度状态下的运动行程，最终导致粒子在发动机内部停留时间的增加；在横向过载不变的情况下，加入轴向过载后，再次改变了粒子的运动轨迹，使得粒子在燃烧室壁面的撞击区域朝着喷管方向移动，增大了粒子与燃烧室壁面撞击和反弹过程的速度，减小了粒子停留时间。因此，针对文中的计算条件，无过载时粒子停留时间最短，单纯横向过载时粒子停留时间最长，在横向过载不变时引入轴向过载后粒子停留时间在两者之间。

2.2 横向过载对壁面粒子分布的影响

图6给出了不同横向过载条件下壁面粒子浓度分布，表2 给出了横向过载方向壁面最大粒子浓度对比情况，其中，Pc表示粒子浓度，x表示发动机轴向位置，y表示发动机径向位置，粒子浓度是通过Workbench自带的功能提取。从图中可以看出，在纯横向过载作用下，喷管扩张段壁面粒子浓度为零，表明粒子通过喷管喉部后不会撞击到扩张段壁面，而是随中心区主流燃气一起流出喷管，即在喉部后形成文献[14]描述的粒子自由区。

从图6和表2可以看出，随着横向过载的增大，过载方向燃烧室筒段和喷管收敛段局部区域最大粒子浓度总体上呈增大的趋势，且沿着发动机轴线方向燃烧室筒段和喷管收敛段粒子聚集区域逐渐扩大。此外，随着横向过载的增大，燃烧室筒段壁面和喷管收敛段壁面存在粒子浓度局部增大的区域，呈现多峰分布现象，峰数与横向过载大小有关，10g～30g以下呈双峰分布，40g时呈三峰分布。

图6 不同横向过载下壁面粒子浓度分布计算结果Fig.6 Calculation results of wall particle concentration distribution under different lateral overload

表2 横向过载下壁面最大粒子浓度计算结果Table 2 Results of maximum particle concentration under lateral overload

分析认为，在横向过载的作用下，粒子朝着过载方向发生偏转，不利于粒子的随流性，增大了粒子与燃烧室壁面和喷管收敛段壁面局部区域的撞击频率，从而引起粒子在燃烧室和喷管收敛段局部区域聚集，导致该区域粒子浓度的增大，且局部区域粒子浓度随着横向过载的增加而逐渐增大。

2.3 组合过载对壁面粒子分布的影响

图7给出了不同横向过载和轴向过载条件下壁面粒子浓度分布，其中a表示为轴向过载，l表示为横向过载。表3给出了不同过载条件横向过载方向壁面最大粒子浓度对比情况。在横向过载和轴向过载共同作用下，喷管扩张段壁面粒子浓度依然为零，表明粒子通过喷管喉部后仍旧不会撞击到扩装段壁面，同样随中心区主流燃气一起流出喷管。

图7 不同横向和轴向组合过载下的壁面粒子浓度分布计算结果Fig.7 Calculation results of wall particle concentration distributionunder different lateral and axial combined overload

从表3可以看出，保持横向过载40g不变，轴向过载引起燃烧室最大粒子浓度向喷管方向移动，轴向过载20g和30g时燃烧室壁面最大粒子浓度更大，轴向过载40g时最大粒子浓度开始下降，故存在与最大粒子浓度相对应的轴向过载临界点；由图7中可以看出，随着轴向过载的增大，沿着发动机轴线方向燃烧室筒段粒子聚集区域先缩小后扩大，在轴向20g时粒子聚集区域最小。随着轴向过载的增大，喷管收敛段最大粒子浓度呈下降趋势，且沿着发动机轴线方向喷管收敛段粒子聚集区域逐渐缩小。此外，在横向过载和轴向过载共同作用下，燃烧室筒段壁面和喷管收敛段壁面粒子浓度均存在粒子浓度局部增大的区域，依然呈多峰分布现象，以三峰分布为主。

表3 组合过载下壁面最大粒子浓度计算结果Table 3 Results of maximum particle concentration under different overload

分析认为，与单纯横向过载相比，轴向过载的加入再次改变了粒子的运动方向。在燃烧室区域，粒子运动速度较低，轴向过载将增大粒子的运动速度，导致粒子向综合过载方向运动并与燃烧室壁面发生撞击，造成局部粒子浓度的增大，最大粒子浓度与轴向过载直接相关；在喷管收敛段区域，粒子运动速度较高，轴向过载加强了粒子的随流性，导致喷管收敛段壁面最大粒子浓度的下降。

研究表明，发动机绝热层烧蚀受粒子冲刷影响较为严重，喷管烧蚀受热化学烧蚀和机械剥蚀的双重作用，而粒子冲刷是机械剥蚀的主因，发动机在过载条件下工作更是如此，故粒子聚集将加剧局部位置的绝热层烧蚀。因此，在发动机设计过程中，针对燃烧室和喷管粒子浓度聚集部位，应该根据不同的过载大小、过载时间适应性地加厚绝热层，从而有效提高发动机的工作可靠性。