塞式喷管流场计算

采用NND差分格式，对塞式喷管轴对称流场进行了数值模拟。通过流场分析，可以对塞式喷管火箭发动机有清楚的认识。计算表明，NND格式有很好的分辨率；无粘和有粘的情况流场结构和性能相差悬珠；塞式喷管的自适应能力很好。     

固体火箭发动机内轴对称两相流场的一体化数值模拟

将固体火箭发动机燃烧室和喷管作为一个整体，应用时间相关法进行了轴对称两相流场的一体化数值模拟。对气相采用ＭａｃＣｏｒｍａｃｋ显式差分格式，对颗粒相则采用特征线法进行计算。计算结果表明该方法可以对带潜入式喷管的发动机内复杂的两相流场进行有效的数值模拟。     

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

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

跨音速喷管化学反应流数值模拟

用时间相关的半隐格式有限差分数值方法求解了化学非平衡反应跨音速喷管流场，在喷管收敛段，流动接近化学平衡状态，控制方程的刚性问题严重，数值积分困难。通过对时间差分项隐式离散、对空间差分项显式离散，流场边界采用参考平面上的特征线法计算等，成功地解决了由于化学反应有限速率带来的数值解不稳定问题。该格式简单、需要计算机存贮空间少。本文完成了一维和轴对称非平衡化学反应喷管流动计算，并与化学平衡流和冻结流的计算结果做了比较。     

串装双燃速药柱发动机的内流场计算

对一维的加质燃烧室流场、轴对称的二维喷管内流场和高燃速药柱燃烧完成后的燃烧室二维内流场以及三维药柱几何做了一体化计算。应用轴对称的N－S方程、TVD有限体积法的显式MacCormack格式，对燃烧室后段和长尾喷管跨声速流场进行求解，得到了双燃速内弹道性能和轴对称二维流场中燃气参数的分布。结果具有较高的精度，并为进一步的燃烧室和喷管热结构分析提供数值依据。     

固液混合火箭发动机喷管流动计算

利用二维轴对称 N-S方程和组分方程对选用液氧 /端羟基聚丁二烯推进剂的固液混合火箭发动机喷管流动进行了计算。计算采用 LU时间隐式格式、MUSCL空间离散和 Van Leer矢通量分裂方法 ,8组分 1 0反应的化学反应模型 ,对化学源相进行了点隐式处理。喷管入口条件通过燃速计算和热力计算得到 ,计算了多个氧化剂流率和装药长度情况下的喷管流场 ,分析了真空比冲和推力与氧化剂流率的关系 ,为固液混合火箭发动机喷管设计提供了依据      

高精度ENO格式在射流数值模拟中的应用

本文分析了ENO格式的特点，并应用于全NS方程的迁移项和压力顶，模拟了轴对称和三维射流。首先对轴对称自由射流为算例比较了高阶TVD格式和ENO格式，证实了ENO格式捕获复杂激波波系的能力；其次计算了四喷管湍射流干扰流场，给出了相应的计算结果；最后计算了四股射流流入圆管内的流场，给出了物理量沿壁面和截面的分布。     

轴对称矢量喷管内流特性的不同湍流模型计算

利用二阶迎风格式，在采用三种不同湍流模型分别对轴对称矢量喷管三维内流场进行对比计算后，S-A-方程模型成为文中数值模拟中最有效的湍流模型，并将此模型计算的结果与模型试验结果进行了分析比较。研究结果表明，在矢量状态下，计算结果和试验结果相比，气动矢量角的相对误差不大于6％，流量系数的相对误差不大于l％，推力系数的相对误差不大于l％，所开发的计算方法可用于轴对称矢量喷管的工程设计。     

多重网格加速的LUSGS算法用于喷管流场数值模拟

采用有限体积法,二阶精度插值的Roe格式对三维流场N-S方程进行空间离散;采用LUSGS隐式时间推进法和多重网格技术加快计算收敛速度;采用独特的显式算法克服了两方程湍流模型刚性对流场计算收敛速度的影响,并给出了该方法的理论依据。通过轴对称喷管内流场、二元喷管内外流场的数值模拟验证了采用以上方法在计算效率上的提高。最后将RNG模型计算结果与SA模型计算结果及实验数据做对比分析。     

定常轴对称超音速喷管和超音速射流的计算

本文叙述轴对称超音速喷管和超音速射流的数值计算。内点用MacCormack有限差分格式,分步加粘性;边界用从特征关系列出的二阶精度格式;喷管口用Prandtl-Meyer流公式;射流边界用一个迭代过程确定。对有激波喷管和无激波喷管做了大量的数值实验,并作了比较和讨论,给出了三个喷管的计算结果。     

在烧蚀条件下固体火箭发动机复合喷管的温度预示

本文探讨了由石墨和高硅氧-酚醛组成的固体火箭发动机复合喷管,在考虑热解炭化和烧蚀引起喷管内型面变化条件下的瞬态轴对称二维温度场.借用当量热容的概念和坐标变换的方法,应用单元体能量平衡原理,在内点用显式格式,边界结点用隐式格式进行数值求解,这样的计算比较简捷,所得温度场与实验结果作了比较,说明本文的预示方法是可以符合工程精度要求的.     

轴对称燃气射流的数值计算

本文应用VonMises变换对轴对称射流的各控制方程和相应的初始条件及边界条件进行了变换, 选用普朗特混合长度湍流模型, 用隐式有限差分格式进行数值求解, 成功地计算了某型飞机在巡航状态下燃气由轴对称收缩喷管排入亚音速空气流中的速度场、温度场、浓度场、密度场以及燃气各组分的成分分布。计算结果可以用于预测飞机尾喷管高温燃气射流的红外辐射特性。      

带推力矢量的涡扇发动机实时数学模型研究

利用流场计算结果建立了轴对称矢量喷管的实是数学模型。将发动机部件热力参数间的关系用显式的解析式关系表示，从而去掉部件计算机的迭代过程，以此方法建立了涡扇发动机实时数学模型。在以上两个模型的基础上建立了带轴对称矢量喷管的涡扇发动机数学模型。利用此模型研究了矢量喷管对涡扇发动机工作参数的影响。结果表明，该模型可完成带推力矢量的涡扇发动机静态及动态计算，并可用于推力矢量控制研究。矢量喷管的偏转对涡扇发动机工作具有一定影响。     

轴对称短化喷管型面设计及流动分析

对最短长度喷管（MLN）设计方法进行了改进和优化，在喷管喉部扩张段实现了光滑过渡，采用特征线方法设计轴对称短化喷管型面，采用高阶，高分辨率WENO格式对设计的喷管进行流场数值模拟，结果表明：短长度喷管和喉部圆弧过渡的短化喷管出口流场均匀，能够达到设计要求，经过圆弧过渡的喷管型面出口流场品质更好。     

二元塞式喷管红外特征及壁面降温的红外抑制效果计算

用数值计算的方法研究了涡扇发动机二元塞式喷管在3~5μm波段的红外辐射特性,并与轴对称收扩喷管进行了对比.结果表明:无冷却措施时,相比于轴对称收扩喷管,二元塞式喷管在大部分探测角度上红外辐射增强,只在小部分的探测角度范围内具有红外抑制作用,在窄边探测面内探测角为90°方向上,其红外辐射强度降低了54%.若塞锥后部壁面温度降低300K,二元塞式喷管在大部分探测角度范围内的红外辐射强度明显低于轴对称喷管,在正尾向上的红外辐射强度降低了57%;与轴对称喷管中心锥表面降低相同温度相比,二元塞式喷管能取得更高的红外抑制效果,并且所付出的冷却代价较小.      

塞式喷管流场数值模拟

采用ＮＮＤ２Ｍ差分格式,从考虑к－ε湍流模型的理想气体三维薄层近似Ｎ－Ｓ方程出发,对塞式喷管内外流场进行了数值模拟。计算结果与国外文献十分一致。通过流场分析,对塞式喷管流动机理有了清楚的认识。计算表明ＮＮＤ２Ｍ差分格式及相应计算软件适合于模拟复杂喷流流场,可用于塞式喷管发动机的研究。     

轴对称塞式喷管电磁散射特性数值模拟

为研究塞锥存在对喷管雷达隐身特性的影响,以涡扇发动机轴对称喷管为基础,设计了加有不同锥度塞锥的塞式喷管;运用自主开发的基于物理光学迭代(IPO)和等效边缘电磁流(EEC)方法的程序对各型塞式喷管的雷达散射特性进行了数值计算并与原轴对称喷管特性进行对比分析。结果表明:塞式喷管能够有效降低喷管雷达散射截面(RCS);但0°仰俯角附近会有局部升高,在水平极化和垂直极化方式下RCS值分别比原喷管最大增大8.14%,11.77%;存在最优锥度,在水平和垂直极化方式下都能够最大程度减小喷管总场RCS均值,使得其比原喷管分别减缩13.2%,15.2%。     

轴对称矢量喷管三维传热计算研究

采用封闭腔理论算法，利用矢量喷管的三维内流场数值计算结果，对轴对称矢量喷管的壁温分布进行了数值计算。计算结果表明，在非矢量状态，喷管同一轴向位置处的周向壁温基本相同，壁温分布可按一维处理。在矢量状态，喷管同一轴向位置处的周向壁温相对差值达15％，壁温分布呈现复杂的形式，必须采用三维计算方法计算壁温分布；喷管扩张段壁温明显高于非矢量状态，所开发的计算方法和程序既可用于轴对称矢量喷管的工程设计，也可为试验验证提供理论依据。     

跨声速条件下轴对称收扩喷管内外流场的数值研究

采用有限体积法,二阶精度插值的Roe格式对三维流场N-S方程进行空间离散;采用LUSGS隐式时间推进法和多重网格技术加快计算收敛速度;采用不同湍流模型、壁面函数对轴对称喷管内外流场进行三维数值模拟.计算结果与实验数据对比表明,采用涡粘修正的RNG模型,内流场使用增强壁面函数,外流场使用标准壁面函数可以得到准确的计算结果.研究表明:在亚声速条件下,喷管外壁面压力分布随落压比的变化不大;超声速条件下,喷管外壁面激波随着落压比的增大前移.      

带有分离的二维可压粘性喷管流动计算

采用二维、时变、可压纳维-斯托克斯(N-S)方程计算喷管流动。流动可为无粘、层流或湍流。湍流模型可采用代数混合长模型,一方程和二方程模型。喷管形状为二维(平面或轴对称)。计算网格内点采用MacCormack显式格式。采用局部细网格,快速解法以及时间间隔光滑方法,以加速高雷诺数下收敛。边界点采用参考平面特征线格式计算,粘性项作为源项处理。给出了初步算例结果,与其它计算和实验数据比较,符合良好。