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

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

低密度小喷管流场的DSMC直接模拟

用ＤＳＭＣ方法数值模拟了低密度轴对小喷管的流场流动，获得了相应的结果，并将所得结果与用差分求解Ｎ－Ｓ方程方法获得的结果进行了对比，二者吻合较好。     

超燃冲压发动机燃烧室三维化学反应流场的数值模拟

采用ＬＵ隐式方法强耦合求解有限体积法离散的完全Ｎ－Ｓ方程及组分方程，数值模拟超燃冲压发动机机燃烧室机的化学反应流流场。计算中和了ＡＵＳＭ通量分裂格式，两方程的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代化学反应流场。     

双燃式冲压发动机中超燃燃烧室冷态流场数值模拟

采用Ｍａｃ Ｃｏｒｍ ａｃｋ 矢通分裂格式, 求解雷诺平均的Ｎ－Ｓ方程组, 采用Ｂａｌｄｗ ｉｎ－Ｌｏ－ｍ ａｘ 代数湍流模型和混合长度模型, 模拟超燃燃烧室内两股平行超声速来流的冷态流场, 取得了与试验相一致的结果, 计算表明回流区和射流流动参数对超燃燃烧室流场内波系的变化有着较大的影响。     

一种模拟喷管二维流场的新方法

采用ＳＩＭＰＬＥ方法对固体火箭发动机喷管进行了纯气相流场的数值模拟。推导了建立在同位网格之上利用协变物理速度分量为计算变量的离散方程，采用协变物理速度分量推导压力修正方程，对密度采用一阶迎风格式插值，将ＳＩＭＰＬＥ方法扩大到计算可压流动。计算结果表明：本方法计算精度高，收敛速度快，程序编写简单。     

复杂几何形状喷管内外三维流场的数值模拟

研究了球面收敛矩形扩张的二元喷管（ＳＣＦＮ）流场的数值模拟,为了有效地捕获激波与滑流面及较好地模拟分离,采用了Ｒｏｅ格式、ＭＵＳＣＬ方法及Ｖｅｎｄａｔａｋｒｉｓｈｎａｎ限制器、Ｓｐａｌａｒｔ－Ａｌｌｍａｒａｓ方程紊流模型。研究中给出了两个喷管模型四个流动状态的数值模拟结果。结果表明,除了有严重分离的流动情况以外,该算法能得到相当满意的计算结果。     

粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

喷管化学反应流的对角化点隐式法求解

采用弱耦合对角化点隐式方法的ＭａｃＣｏｒｍａｃｋ格式,求解喷管的粘性化学反应流动,解决了控制方程组的刚性问题,又大大减少了常规点隐式方法中求逆矩阵的大量工作。对液氢液氧火箭发动机,采用６种组分、８个反应有限速率的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流场参数在喷管中的分布。计算结果表明,采用对角化点隐式方法求解化学反应喷管粘性流动能提高计算效率。     

一维摩擦喷管流动

本文阐述了一维摩擦喷管流动控制方程的数值求解方法,分析研究了摩擦对喷管流场及发动机性能的影响。确定边界条件时考虑了燃烧室工作过程,并给出了计算例题。     

次流喷射控制推力矢量喷管的流场数值模拟

为了探讨次流喷射控制二维推力矢量喷管的性能,用时间推进求解Ｎ－Ｓ方程的方法数值模拟了喷管的流场,数值格式用中心有限体积格式,借助于当地时间步长和隐式残差光顺技术加速收敛,得到了在不同二次喷流压力下推力矢量喷管内流场的变化特征。计算结果与实验结果比较有满意的一致性。研究表明,应用次流喷射控制主流流动可以实现较大的推力矢量转折,但是,二次喷流必须具有足够的压力值。     

连续旋转爆轰发动机参数特性的三维数值模拟

采用一步化学反应模型,基于任意坐标系下Euler控制方程,对连续旋转爆轰发动机(CRDE)进行了三维数值模拟研究, 详细分析了来流总压对CRDE参数特性的影响.研究发现:随着来流总压的增大,爆轰波峰值随之增大,但是燃烧室头部Laval型喷注段的壅塞比基本不变.不同来流总压下,可燃气体在燃烧室头部均以亚声速入射.随着来流总压的增大,燃烧室内流场的平均压强增大,但是其平均轴向流速基本不变.CRDE的流量、推力和比冲均随着来流总压的增加而变大;但是,不同来流总压下,CRDE的流量与Laval型喷注段最大流量之比不变,并且此比值约等Laval型喷注段的壅塞比.Laval型喷注段最大流量与来流总压成正比,因此以上分析从理论上进一步解释了CRDE的流量与来流总压成正比的原因.      

喷水射流预冷对发动机进气温度及流场变化影响的数值研究

为减小整个预旋系统的流动损失,首先对带预旋集气腔进气孔、预旋集气腔、预旋喷嘴的冷气预旋流路进行了分析,发现进气孔和集气腔会导致预旋喷嘴进口流场不均匀,相较于进口均匀条件,预旋喷嘴总压损失系数增大0.026。在此基础上提出了一种将预旋集气腔进气孔、预旋集气腔和预旋喷嘴融合设计的低损失融合式预旋喷嘴设计方案,分析表明：融合式预旋喷嘴能有效减小冷气在预旋系统内的流动损失,在设计工况总压损失系数减小0.032,并使冷气在预旋系统内流动更加均匀,提升了预旋系统的整体性能。     

低损失融合式预旋喷嘴设计与研究

为减小整个预旋系统的流动损失,首先对带预旋集气腔进气孔、预旋集气腔、预旋喷嘴的冷气预旋流路进行了分析,发现进气孔和集气腔会导致预旋喷嘴进口流场不均匀,相较于进口均匀条件,预旋喷嘴总压损失系数增大0.026。在此基础上提出了一种将预旋集气腔进气孔、预旋集气腔和预旋喷嘴融合设计的低损失融合式预旋喷嘴设计方案,分析表明:融合式预旋喷嘴能有效减小冷气在预旋系统内的流动损失,在设计工况总压损失系数减小0.032,并使冷气在预旋系统内流动更加均匀,提升了预旋系统的整体性能。     

旋流对轴对称喷管性能影响的数值研究

为了研究旋流进气对喷管性能参数的影响,使用数值模拟的方法对轴对称喷管直流进气与旋流进气两种工作状态进行了对比研究,分析了不同空气流量、不同旋流强度下的喷管流量系数、真空推力系数以及比冲的变化规律.研究发现,使用喷管入口轴向质量平均总压计算所得喷管性能参数不随旋流强度的变化而变化,且与直流状态近似相等;旋流的切向动压随旋流数的增大而增大,但对喷管推力没有贡献,应将其归入总压损失.      

波瓣喷管引射-混合器的数值研究与验证

对基于Navier-Stokes方程理论预测波瓣喷管引射一混合器引射流量比的计算方法进行了探讨，在计算过程中．主流进口采用速度边界条件，二次流进口采用总压压力边界条件，混合流出口采用静压压力边界条件，两者均设置为环境大气压力，与相关实验数据的对比验证表明计算结果与实验结果仅相差10％左右；同时通过改变混合管结构参数．得到了该参数对引射一混合特性的影响规律，进一步揭示了波瓣喷管有利于强化引射一混合的内在机理，计算结果符合物理过程本质。研究表明，本文的计算方法可以有效地预测引射流量比和揭示混合流场特性。     

抛物线喷管型面参数对流动分离影响的数值模拟

为了研究抛物线喷管型面的3个参数——初始膨胀圆弧半径、入口角和出口角的变化对喷管流动分离过程的影响规律，模拟了不同参数组合下的抛物线喷管流场参数，获得了喷管流动位置、分离模态随集气室压强增加的变化过程.结果表明:这3个参数能够影响发动机起动过程中喷管流场的发展过程.在一定范围内分离模态转换喷管压比（NPR）与初始膨胀圆弧半径呈正相关；通过合理设定入口角和出口角，可以推迟和缩短受限激波分离过程，为解决彭冠侧向载荷问题提供了可能的方向.      

波瓣喷管引射-混合器的数值研究与验证

对基于CFD数值预测波瓣喷管引射 混合器引射流量比的计算方法进行了探讨,在计算过程中,主流进口采用速度边界条件,二次流进口采用总压压力边界条件,混合流出口采用静压压力边界条件,两者均设置为环境大气压力,与相关实验数据的对比验证表明计算结果与实验结果仅相差10%左右;同时通过改变混合管结构参数,得到了混合管结构参数对引射 混合特性的影响规律,进一步揭示了波瓣喷管有利于强化引射 混合的内在机理。在波瓣喷管出口对应于波谷区域存在一个相对低的静压区,对于引射 混合器系统存在一个最优的性能设计点,计算结果符合物理过程本质。     

气动塞式喷管底部二次流特性的数值模拟

发展了一套可以对底部有二次流的气动塞式喷管的性能进行预示的数值方法，采用二阶精度的NND格式求解二维层流N-S方程，对底部二次流的入口边界条件进行了特殊处理。以一种实验用的直排式气动塞式喷管为对象，针对两种不同的背压条件，对二次流流量分别为主流的0%，1%，3%和5%的几种工况研究了塞式喷管底部二次流流量变化对喷管性能影响。计算结果表明，底部二次流的加入使得气动塞式喷管性能有比较明显的提高。数值方法可以用于气动塞式喷管的设计和性能预报。     

进口为均匀剪切流的二元亚音扩压器

本文在主流为均匀流的二元亚音扩压器流动计算方法的基础上,发展了一种进口为均匀剪切流的二元亚音扩压器流动特性计算方法。为验证本方法,计算了三种实例,并和试验测定结果进行了比较,二者符合情况良好。     

斜切反喷管性能计算

本文提出了斜切反喷管流场计算方法,内部点用MacCormack二步显格式数值解控制方程,入口边界和固壁边界的参数用物理边界条件和可应用的沿双特征线的相容性方程计算.固壁边界计算中采用了激波装配技术和开关格式,有效地捕获了流场内的激波系列,计算中得到的激波情况和物理分析一致.本文在流场计算的基础上,进行了反喷管的性能计算.本文提出的方法可应用于反喷管的性能预估.