环形火焰引发爆震的数值研究

采用改进的适于多组分反应流的波传播算法,对直管道中甲烷．空气预混气的爆震引发过程进行数值研究,考虑了两种不同构型的初始环形火焰对引发爆震的影响。对这两种构型下爆震引发过程中的压力场、温度场、反应放热率以及组分浓度进行了计算,由此探讨了“热点”形成机制及其对爆震引发的作用。结果表明,火焰发展产生的燃烧压缩波在管壁或管轴心发生汇聚、碰撞所形成的复杂波系与火焰相互作用,导致火焰前锋与压缩波之间的局部区域反应剧烈,反应放热率明显增加,从而促使热点形成。形成的热点迅速放大成为过驱爆震波,并逐渐衰减为稳定爆震波。稳定爆震波阵面参数与CJ理论结果进行了对比,两者符合得较好。     

某型导弹尾喷流形状的数值模拟

从三维薄层近似N—S方程出发，采用EN0差分格式，对某型导弹层流场进行了数值模拟。计算模拟了飞行马赫数Mα0．8—1．6，海拔高度H＝2km—15km，燃烧室总压为20MPa和40MPa条件下导弹的层流场。计算结果表明与飞行速度相比，海拔高度和燃烧室总压对层流宽度有显著影响。通过对流场结构的分析与测量，建立了层流宽度与上述参数的拟合公式，为工程估算这一类型导弹层流影响边界提供参考。     

气动谐振管加热振荡流动过程数值研究

为了了解气动谐振管加热现象的振荡流动过程，利用显式TVD－MacCormark格式数值求解二维轴对称雷诺平均N－S方程，获得了一个谐振循环中完整的喷嘴－谐振管系统的振荡流动湍流流场，结果表明，欠膨胀声速冲击喷流固有的流场不均匀性和不稳定性是谐振管加热现象产生的内在原因，激波耗散是谐振管主要的加热机理。     

稳态等离子体推力器磁场设计与数值分析

首先分析了稳态等离子体推力器性能对磁场的要求，据此介绍了磁路系统的工程设计原则和方法；然后，采用ANSYS大型有限元分析软件对一个具体的稳态等离子体推力器在额定工况下的磁场进行了计算，获得了满意的结果，通过对结果的分析获得了对有关现象的直观深入的认识，为所计算具体推力器的改进设计提供了线索；同时，也证实了以ANSYS为软件平台的稳态等离子体推力器磁路系统计算机数值仿真辅助设计的可行性和有效性。     

椭圆截面喷管排气辐射热流数值分析

文中对直升机红外抑制器椭圆截面喷管排气的辐射热流进行数值分析,并预估排气温度分布、速度和燃气成分分布对辐射热流的影响。为了计算辐射热流和考虑尾焰对周围环境的传热,采用形式简单并可用通用方程计算程序计算的六热流辐射模型。由于排气流场内温度和密度变化较大,故紊流模型采用考虑密度变化的修正双方程k-ε模型。本文中所用的比热C_(pm)是随气流温度和成分而变化,温度用迭代法求得。喷管出口截面为椭圆形,排气为三元自由射流,控制方程为三元抛物形偏微分方程,可采用前进积分求解,这样既可满足计算精度要求,又可节省机时。并根据三元喷流周向边界条件周期性变化的特点,采用CTDMA解法。计算结果与试验结果基本一致,说明本计算方法是可行的。     

外形设计对微喷管性能的影响

利用直接模拟蒙特卡罗法对钟形、喇叭形和锥形等三种不同形状微喷管进行管内气体流动模拟,并对不同温度及入口压力下各形状喷管的性能变化进行模拟分析,模拟中采用与实际分子碰撞过程更为接近的内能松弛模型及CLL物面反射模型。研究结果表明钟形微喷管因扩张段回流区的存在,其推力性能明显低于喇叭形和锥形微喷管,锥形微喷管的推力性能同喇叭形微喷管较为相近。三种形状的微喷管推力值都随着入口温度的升高而降低,而微喷管的推力效率随着入口温度的升高而上升,入口压力的变化对微喷管推力有较大影响,但对微喷管的推力效率影响并不显著。     

聚焦脉冲激光大气传输热晕效应的数值分析

利用激光大气传输四维程序对序列脉冲激光大气传输热晕效应进行了数值计算和分析。对于均匀大气路径上聚焦传输序列脉冲激光进行数值仿真计算得到,风场渡越时间内的脉冲数Np在2~6之间有利于传输;然后对实际大气模式下序列脉冲激光进行了数值计算,得到的平均脉冲数-NP可以在一定条件下对序列脉冲激光进行定量表征的结论。     

强制偏流喷管的数值仿真

采用薄层Ｎ－Ｓ方程成功地求解了强制流喷管内的三维复杂流场，对通量项和压力项用Ｖａｎ Ｌｅｅｒ发展的矢通量分裂法进行了离散，对离散方程采用线超松弛方法迭代求解，在亚音速与超音速气流共同存在的出口边界分情况进行了处理。计算结果验证了强制偏喷管具有自动补偿外界压力变化的能力，与实验结果基本吻合。     

H型叶栅通道流计算网格的改进

提出一种壁面正交、尾迹区网格节点分布逐渐均匀化叶栅通道网格生成方法。为了研究壁面处网格正交性对数值计算结果的影响，推导出壁面附近流动参数一阶导数离散误差与网格线交角的数学关系。采用N-S方程分别对平板紊流附面层流和叶栅通道流进行计算分析。结果表明：壁面处网格正交性对附面层内流动计算预测精度影响较大，而对势流区流动影响很小；尾迹区切向网格节点逐渐均匀化分布能更准确计算尾迹区流动。     

分区和有限容积法计算长尾喷管中的湍流流动

为了给长尾喷管热防护结构设计提供热边界条件，从轴对称二维雷诺平均N－S方程出发，采用分区算法和有限容积法计算了长尾喷管中的湍流流动，分析了流场的结构特性，提出了长尾喷管中燃气温度，压强和马赫数分布的某些规律以及燃气温度沿壁面的分布，计算和分析结果提出了长尾喷管热防护结构设计应注意的事项。