喷管流场和羽流场的N－S数值模拟

利用时间推进的方法，通过求解可压缩的、雷诺平均的Ｎ－Ｓ方程，并利用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数湍流模型和Ｐｒａｎｄｔｌ混合长度模型，数值模拟了飞机和火箭发动机喷管的内流和羽流，获得了喷管内流场和羽流场的参数分布，计算结果与试验结果及国外研究成果一致。网格划分采用了自适应网格技术。     

喷管流场和羽流场的N－S数值模拟

利用时间推进的方法，通过求解可压缩的、雷诺平均的Ｎ－Ｓ方程，并利用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数滴流模型和Ｐｒａｎｄｔｌ混合长度模型，数值模拟了飞机和火箭发动机喷管的内流和羽流，获得了喷管内流场和羽流场的参数分布，计算结果与试验结果及国外研究成果一致。网格划分采和了自适应网格技术。     

固体火箭发动机内流场分区耦合数值计算

采用分区耦合方法计算固体火箭轴对称燃烧室与喷管流场。对于低速的燃烧室流场，选用不可压流的Ｎ－Ｓ方程描述并用ＳＩＭＰＬＥＣ方法数值求解；对于高Ｒｅ数的喷管流场，则采用Ｅｕ－ｌｅｒ方程描述并用ＳＣＭ方法求解。计算时用燃烧室出流为喷管流场提供入口参数，同时用喷管流场压强分布反馈影响燃烧室流动状况。对耦合边界条件处理方法进行了探讨。对典型的侧壁加质燃烧室与喷管流场进行了计算，计算结果揭示了单独喷管流场计算难以反映的喷管收敛段近壁区的低速区域，与已有的燃烧室流场实验结果一致并反映了燃烧室与喷管流场之间的联系，较好地模拟了流动中的物理现象。     

C-H-O-N体系液体火箭发动机喷管流场数值模拟

采用弱耦合点隐式方法的MacCormack格式对C－H－O－N体系的液体火箭发动机喷管粘性流场进行了数值模拟，采用了考虑化学反应和不考虑化学反应的两种模型，其中化学动力学模型为C－H－O－N体系的12组分、14反应有限速率的化学反应模型。数值模拟得到了流场参数在喷管中的分布，结果表明，数值模拟结果与理论分析一致，为发动机的设计及试验提供了参考。     

二元矢量喷管流场的N－S方程计算

对二元矢量喷管在不同工况下的流场进行了数值模拟。求解的控制方程为三维非定常Ｎ－Ｓ方程。空间方向用有限体积离散，时间方向用Ｒｕｎｇｅ－Ｋｕｔｔａ法推进。粘性项的计算上采用了一种有限体积框架下新的离散化方法。湍流模型采用了广泛使用的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模型。计算表明，本文方法能准确地计算二元矢量喷管、非矢量喷管在设计与大设计落压比下的流场。     

真空羽流场的 N-S 和 DSMC 耦合数值模拟

根据调姿发动机真空羽流场的特点,采用差分求解Ｎ－Ｓ方程数值模拟喷管内流场,用ＤＳＭＣ方法数值模拟喷管外羽流场,得到了符合流动规律的数值模拟结果。为姿控发动机喷管设计提供了参考,并为真空羽流污染研究打下了良好的基础。     

喷管分离流动的数值模拟

针对高空喷管地面试验时产生流动分离的现象,用Ｂｅａｍ －Ｗａｒｍ ｉｎｇ 近似因式分解法求解薄层Ｎ－Ｓ方程,对分离流场进行数值模拟。预示了分离点的位置以及喷管壁面压强分布等,可以为高空发动机地面试车提供参考, 并为进一步研究分离的影响因素打下基础     

斜切反喷管粘性流场数值模

为了在级间分离期间提供反推力，许多固体火箭发动机前端都装有一组斜切反喷管，由于反喷管气动型面具有尖点，并在超音速区有台阶，喷管内存在一系列激波，并伴有流动分离现象。从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，利用时间相关法，采用ＭａｃＣｏｒｍａｃｋ两步显格式，结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，数值模拟了斜切反喷管流场。计算得到的壁面压强分布与风洞吹风实验测得的压强分市相当一致。     

斜切反喷管粘性流场数值模拟

为了在级间分离期间提供反推力，许多固体火箭发动机前端都装有一组斜切反管，由于反管气动型面具有尖点，并在超音速区有台阶，喷管内存在一系列激波，并伴有流动分离现象。从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，利用时间相关法，采用ＭａｃＣｏｒｍａｃｋ两步显格式，结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，数值模拟了斜切反喷管流场。     

真空小喷管羽流场的Monte Carlo直接模拟

利用直接模拟的Ｍｏｎｔｅ Ｃａｒｌｏ方法数值模拟了真空喷管羽流场，获得了符合规律的数值模拟结果。由于将随机取样频率法用于轴对称问题中，解决了计算机内存和速度上的困难，同时也证明了Ｒ－Ｓ－Ｆ方法的可靠性。     

喷管中泡状气液两相流影响因素分析

为了获得边界条件对缩放喷管中泡状流的影响,采用包含气泡动力学的连续均质混合流模型,基于求解Rayleigh-Plesset方程的变步长龙格-库塔方法,提出一种稳态泡状流求解的迭代修正计算方法,实现了喷管中流场计算。获得了气泡半径、混合密度在喷管中变化情况。将不同边界条件下的计算结果进行对比分析,得到了气泡平衡半径、喷管入口含气率及入口速度对流动的影响。     

SIMPLE方法的敛散性

 通过对轴对称突扩直管和突扩后缩环管流场的计算考察了SMPL正方法的敛散性。在源项处理、方程的非线性问题和解法等方面进行了分析。在此基础上研究了迭代初场的构成和斜壁网格的处理。提出了一种可行的算法。     

火箭发动机塞式喷管流场的数值研究

用时间相关法和张涵信的无振荡、无自由参数ＮＮＤ显式格式,对火箭发动机塞式喷管二维轴对称定常流场进行了数值模拟。计算得到了有粘和无粘中两种情况下,喷管流场的马赫数压强、温度、流线以及速度矢量分布图,计算网络采用代数方法生成,并通过求解Ｌａｐａｃｅ方程对网格进行修正。计算结果表明：ＮＮＤ格式用于计算塞式喷管流场具有较高的效率和良好的分辨率;无粘和有粘情况得到的流场结构有较大的差别;基础喷管的出口马赫数     

爆震管喷管流动数值研究

为了获得高性能的爆震管,应用时空守恒元和解元(CE/SE)法对带不同形式喷管的爆震管进行二维轴对称数值模拟,分析复杂的爆震波流场,研究喷管对爆震管性能的影响。结果表明,爆震管的推进性能的提高受喷管的长度、夹角和形状的影响,其中长喷管性能优于短喷管,收敛扩张喷管性能优于收敛喷管,钟形喷管性能优于扩张喷管。     

喷管气动参数对推力矢量影响的数值模拟

以二元收扩喷管为对象,开展了基于二次流喷射的流体推力矢量技术研究。基于CFD技术,分析了激波矢量控制技术实现推力矢量的机理。重点分析了喷管落压比NPR,二次流总压比SPR,自由流马赫数Ma∞对流体推力矢量性能的影响。数值模拟结果表明:推力矢量的大小与斜激波的位置、角度以及流动分离区的大小有关。在所计算的参数范围内,推力矢量随喷管主次流的变化规律是:推力矢量角随NPR的增大逐渐减小;随SPR的增加,推力矢量角单调增加;在大落压比时,自由流马赫数对推力矢量的影响是有限的,而在低落压比时,自由流马赫数增加,推力矢量角减小。     

喷管一维流场化学非平衡通用计算软件研制

研制的火箭发动机一维喷管非平衡流通用计算软件，建立在一个统一的化学反应库之上。该通用计算软件可对推进剂元素为碳、氢、氧、氮、氟、氯的火箭发动机喷管进行一维化学非平衡流数值模拟。通过计算推进剂元素为氢、氧和碳、氢、氧、氮的喷管一维化学非平衡流，得到了喷管流场的速度、温度、密度、压力及各组分的质量百分比等参数，为工程设计提供了依据。     

离心式喷嘴内流场特性的数值模拟

基于单流体模型及气液两相混合分数的概念，建立了离心式喷嘴内流场的数学模型，并进行了数值计算，计算结果给出了离心式喷嘴内部液膜与气涡共存的流场结构及气液交界面的几何形状，计算得到的壁面压力分布、出口处液膜内轴向速度分布以及雾角大小都与实验结果吻合的很好，与常用的跟踪气液界面的ALE方法相比，本文建立的数值模拟方法具有简单，计算量小和易于实现的特点。     

多喷管射流气动声学特性的数值研究

采用大涡模拟结合动力Smagorinsky亚格子模型对四喷管超声速高温射流进行了三维非稳态数值模拟。基于非稳态计算的结果,利用FW-H面积分方法对四喷管射流远场的声学特性进行计算,并将结果与单喷管得到的结果进行对比。结果表明,多喷管射流形成了复杂的流场,其产生的气动声场并不仅仅是单喷管简单的叠加,由于四股射流之间的相互影响使得声场也产生了新的特点,不但表现在总声压级数值上,还对改变了射流噪声辐射的指向性,最大声压级的位置由原来的50°左右变化到30°。     

单元塞式喷管结构对性能影响的数值研究

从N－S方程出发，采用二阶精度的NND格式对塞式喷管的流场进行了数值模拟，重点研究了单元喷管的倾角的塞锥的截断长度在偏离设计高度时对性能的影响。结果表明，不管是在低于还是高于设计高度，随着倾角的增加，内喷管对推力的贡献在减小，塞锥和底部对推力的贡献却增大，而侧喷管的最佳倾角却基本相同，都在全长型塞式喷管设计角度附近。截断后的塞式喷管虽然减小了重量，但不管是在低空还是高空，性能均低于全长型塞式喷管，因此增加底部的推力是提高性能的一个重要途径。