对《固体火箭发动机喷管烧穿故障的试验分析》质疑

１引言吕振中等人署名的“固体火箭发动机喷管烧穿故障的试验分析”一文（刊在《推进技术》１９９７年第１８卷第６期上）,从内容分析,所述对象为我部某大型固体火箭发动机地面试车中的喷管烧穿故障。该文（以下简称吕文或吕课题）认为喷管的内型面半径对轴线的导函数在...     

尾部喷流对飞行器阻力影响的数值模拟分析

底阻在弹类飞行器阻力中占比较大,准确预示底阻对于弹类飞行器飞行性能评估至关重要,而发动机尾部喷流对底阻影响明显。采用基于雷诺平均Navier-Stokes方程的流场仿真方法研究了飞行器底部发动机喷流和外流干扰流场特性,主要分析了喷流对飞行器阻力的影响。飞行器安装了两台推力可调的液体火箭发动机,发动机在不同飞行工况下采用不同推力工作。分别研究了亚音速、跨音速和超音速典型飞行工况下弹体底部无喷流状态、单喷管喷流状态和双喷管喷流状态时,飞行器阻力变化情况。结果表明:不同马赫数下,发动机喷流对底部阻力影响情况基本一致,与无喷流情况相比,当发动机工作时,无论单喷管喷流还是双喷管喷流状态,底部发动机喷流引射效应明显,弹体阻力系数明显增加。     

运载火箭发动机喷管变形机构型综合分析研究

为有效缩短箭体长度,增大喷管扩张比,提高发动机比冲,对运载火箭发动机进行了喷管变形机构型综合分析。首先分析了国内外喷管变形机构发展历程,在此基础上系统总结了喷管变形机构的典型应用,并对典型喷管变形机构进行了方案分析,研究了其组成原理和功能设计特点,最后基于机构拓扑图分析方法给出机构构型方案及设计方法。介绍了国内外相关喷管变形机构各构型方案运动原理及设计特点,并综合分析评价了各方案,提出了喷管变形机构技术发展建议,可为运载火箭发动机变形机构设计提供支撑。     

超声速边界层转捩实验可控扰动的引入方法

通过数值模拟的方法,研究了在二维平板超声速边界层中,由壁面处安装的激振薄膜引入的扰动的演化过程.主要得到了如下结论:由激振薄膜引入的扰动能很快演化成标准的T-S波,但扰动的各物理量的演化为T-S波的速度并不相同,即各物理量并不同时形成T-S波,速度最先形成,而温度最晚.在相同的条件下,由不同激振形式或不同宽度的薄膜引入的扰动的演化速度基本相同.     

长尾喷管两相流流场计算

为了给长尾喷管热防护结构的烧蚀和温度场计算提供准确的热边界条件，从二维轴对称N－S方程出发，采用颗粒轨道模型和有限体积的Jameson格式计算了长层喷管纯气相流场和两相流流场，分析和比较表明：在长尾喷管喉部附近以及扩散段的相同轴向位置，两相流情况下轴线上的燃气马赫数小于纯气相情况下的燃气温度；在进行长尾喷管热防护结构设计时应以两相流流场作为热边界条件。     

基于Mohring声类比发动机喷流噪声数值计算

为了研究非均匀流场中航空发动机喷流噪声特性,满足未来的航空器噪声适航要求,采用计算流体力学（CFD）和计算航空声学（CAA）相结合的混合数值算法,对大涵道比发动机喷管的简化缩尺模型进行气动噪声仿真计算。采用大涡模拟（LES）计算喷管的瞬态喷流流场;在流场计算的基础上使用Mohring声类比进行声源提取,将时均流场插值到声学网格作为背景流,结合有限元和无限元方法对喷流噪声近场以及远场的辐射特性进行数值计算及分析,并通过单通道锯齿形喷管试验验证数值计算方法的可行性。数值结果表明：发动机喷流噪声主要是由内外涵剪切层内的涡环破碎产生的大尺度涡而形成的,噪声辐射峰值主要集中在低频范围内,随着频率升高,各方向角的声压级都在降低,在1000～2500 Hz,从125 dB快速降低到105 dB,之后衰减速度变缓,到100 dB趋于稳定。数值计算方法精确度高,最大计算误差为1.97%。为发动机噪声适航提供了一种噪声预测方法。     

惯性颗粒和壁湍流的相互作用:湍流调制及颗粒分布

直接数值模拟壁面湍流与惯性点粒子方法被用来研究粒子对湍流的调制机理以及粒子在湍流场中的分布.由于在壁面湍流中,内层(约为y+<100,y+为黏性尺度下的壁面垂直方向坐标)与外层(约为y+>100)中存在着两种特征结构,即内层中的大尺度结构(large-scale structures,LSMs)与外层中的超大尺度结构(very-large-scale structures,VLSMs),这两种特征结构控制着壁面湍流的动力学过程.因此研究方向主要是探究粒子与LSMs以及VLSMs之间的相互作用:即通过研究粒子惯性效应对LSMs主导的湍流内层自维持过程来解释不同惯性的粒子导致湍流增强或减弱的原因;通过研究粒子惯性效应对VLSMs主导的湍流外层的调制,提出粒子的惯性效应导致VLSMs增强的直接与间接机理.同时,通过人为控制LSMs和VLSMs与粒子分别耦合来研究这两种特征结构的存在导致的粒子聚集与群聚行为.最后,推导了惯性粒子沉积速度的控制方程,用于分析控制粒子沉积速度的关键因素.     

超音速分离线喷管内流场数值模拟

采用数值方法求解超音速分离线(SSSL)喷管内流场,研究了不同摆角对喷管流场分布的影响,对比分析超音速分离线与亚音速分离线喷管的轴向推力、径向推力及偏转放大因子随喷管摆角的变化规律,为超音速分离线喷管的设计研究提供理论参考。计算结果表明,摆动对超音速分离线喷管内流场影响显著,随着摆角的增大,内流场的非对称性和激波强度均增加;在相同摆管的轴向力分力略有减小,而径向分力则呈现增大的趋势;超音速分离线喷管与亚音速分离线喷管的径向分力比值,即偏转放大因子则随喷管摆角呈先增大、后减小的变化规律,本算例中的最佳放大因子1.36,对应的喷管摆角为2.5°;另外,随着摆角增大,超音速分离线喷管内流场Al2O3粒子分布的非对称特性也逐渐加强,活动体小端局部范围粒子浓度显著增大。     

关于混合气体输运系数计算方案的影响研究

为了研究混合气体分子输运系数的不同近似计算方案对边界层及射流剪切层流场计算的影响,以含双方程k-ωSST湍流模式的质量平均Navier-Stokes方程为控制方程,分别采用两种计算方案:公式计算混合气体导热系数、扩散系数的直接计算方案和由普朗特数和施密特数等关联混合气体导热系数、扩散系数的关联计算方案,对Burrow和Kurkov实验流场进行了数值模拟。结果表明:与采用直接计算方案相比,采用关联计算方案,流场数值模拟在获得相当的计算结果的同时,大大降低了计算量;抛开与控制方程的相容性问题,混合气体扩散系数各直接计算方案对流场参数计算影响不明显;在湍流计算时,边界层/剪切层区域的混合气体湍流输运系数大于分子输运系数,由此弱化了混合气体分子输运系数不同的近似计算方案对流场数值模拟的影响;相对输运系数近似计算方案,湍流可压缩修正模型对组分剖面的影响更为明显。计算表明,在边界层近壁流场、剪切层流场以及层流流场的精细化数值模拟中,混合气体分子输运系数的计算方案有重要的影响。     

简化三方程转捩模型在高超声速流动中的应用

在雷诺平均方法的基础上,通过耦合k_(-ω)SST湍流模型和间歇因子转捩模型,引入湍流模型和转捩模型的可压缩修正方法,对高超声速平板、双楔、尖锥三类模型边界层转捩流动开展了数值模拟研究。与实验结果的对比分析表明基于压力梯度表征参数T_w=R_(TΩ/ω)的简化三方程转捩模型,能够准确捕捉高超声速平板边界层流动的转捩起始位置、转捩区域长度以及湍流区壁面热流。而对于双楔、尖锥模型,改进前的简化三方程转捩模型由于受到流动可压缩效应的影响,边界层转捩后湍流区的壁面热流模拟预测结果明显高于实验值。在添加模型可压缩修正方法后,转捩区域长度和湍流区壁面热流模拟结果得到有效改善,与实验值吻合较好。可见,简化三方程转捩模型在添加可压缩修正方法后具备准确模拟预测高超声速边界层流动转捩的潜力。