超声速前机身及翼—身组合体的全速势方程数值计算

通过数值求解全速热方程，计算了超声速来流的前机身及翼－身组合体。当流场存在亚声速区时，在此区域内采用中心差分格式，迭代求解，并引入多重网格技术，加快收敛；当流场中某些一区域沿某一方向是超声速时，在此区域内采用沿该方向的推进解法。计算结果表明，本文的方法可靠，结果准确，可以向工程应用方面推广。     

超声速和高超声速横向喷流数值模拟

为了解横向喷流的干扰影响和喷口附近的流场结构,采用分块对接网格和“0”型网格技术,数值求解N—S方程来模拟超声速和高超声速流场中横向喷流的干扰流场。对两种尖拱弹身外形的超声速和高超声速喷流干扰流场进行了数值计算,计算结果与风洞实验数据吻合一致。在此基础上,开展了某型导弹多喷构型干扰流场的数值模拟,得出一些横向喷流数值模拟的结论。     

给定下游边界的超声速流场逆向求解方法

为了实现超声速流场出口参数的可控设计,借鉴特征线方法在古尔赛特与柯西问题两类边值条件中的应用,提出了一种根据下游边界条件求解其依赖域的逆柯西问题求解方法,并在此基础上发展了一种均匀来流条件下的超声速流场逆向求解方法。对轴对称内收缩流场、二维平面流场、轴对称外锥流场三个典型算例校核表明,逆向求解方法得到的计算结果与正向特征线方法一致,由此验证了该逆向求解方法的可行性,并有效拓展了超声速流场的设计思路。     

基于各向异性非结构网格的超声速流动自适应计算

基于各向异性非结构网格,实现了自适应求解技术,将其应用于超声速楔形体绕流及超声速横向喷流问题的计算研究,对各向异性网格自适应计算过程的收敛性和求解精度进行了对比分析,展示了各向异性网格自适应算法在降低问题求解规模。提高各向异性流场分辨率等方面的优势及存在的问题。超声速楔形体无粘绕流场的自适应计算,清晰地捕捉到了激波附近的流场信息,激波前后的马赫数、压力、密度。温度等物理量与理论分析解吻合很好。超声速横向喷流流场存在激波、分离涡、边界层等流场结构的相互干扰。计算研究表明,单纯基于Hes-sian度量张量的各向异性网格生成及自适应算法不能有效模拟边界层内的流动情况,是将来需要进一步开展的研究挑战。     

尖拱弹身横向喷流数值模拟

本文通过数值求解层流Navier—Stakes方程来模拟超声速流场中横向超声速喷流的干扰流场，采用分块对接网格和“O”型网格技术，精确模拟喷口截面形状并生成高质量的贴体计算网格，对尖拱弹身的横向超声速喷流的干扰流场进行了数值模拟，研究了喷口附近的流场结构、涡系结构、波系结构及流动分离等流场特性，得到喷流流场的干扰效应，计算结果与实验值吻合较好。     

强激光与高超声速球锥流场干扰数值模拟研究

本文采用有限差分技术，数值求解带能量源项与空气11组分化学动力学模型的轴对称粘性流体力学方程组，采用类氢原子理论模型计算空气对激光的吸收系数，建立了模拟强激光与超声速钝体流场干扰的计算软件，得到了高能浙江与超声速球锥流场干扰数值模拟的结果。     

外压式超声速进气道流场的数值模拟

采用CFD方法对某型外压式超声速进气道流场进行了二维数值计算,叙述了网格生成、边界条件确定方法,控制方程求解收敛技巧以及计算结果的分析.数值计算结果与理论相符较好.     

超声速气流中液滴破碎的大涡模拟

为了研究超声速气流中的液滴破碎过程,采用基于CLSVOF界面追踪方法的大涡模拟方法。该方法中用可压流和不可压流求解器分别计算超声速气流流场和液体流场,并且在分别计算气相或液相时,另一相作为该相计算时的边界条件参与计算,从而实现可压流求解器和不可压求解的融合。通过去散度化液体速度外推解决界面处大液气密度比所带来的动量误差大的问题。利用该方法对来流Ma=1.358的超声速气流中的水滴破碎进行了数值模拟。模拟结果展示了超声速气流中液滴破碎的几何结构变化细节和基本破碎形态,且与实验吻合较好。数值模拟结果揭示了Rayleigh-Taylor不稳定和气动剪切力对液滴破碎的作用过程。破碎过程中,液滴在垂直来流方向上的最大无量纲宽度(Dmax/D0)约为4.5,无量纲破碎时间(tbU∞/D0)约为29.85。     

超声速/高超声速风洞试验段结构形式对流场性能的影响研究

通过求解轴对称 N-S 方程,对Φ1 m 高超声速风洞马赫数3和6状态下的流场进行了模拟,计算结果与试验数据基本一致,验证了所用数值方法的可信性。在此基础上,对比研究了马赫数3和6状态下采用闭口等直圆截面和开口自由射流两种试验段结构形式的超声速/高超声速风洞在起动条件下的稳态流场性能。结果表明：采用闭口等直圆截面试验段和开口自由射流试验段的流场均匀区内速度场性能指标均满足相关标准要求;马赫数3喷管采用闭口试验段时,沿风洞轴向-300mm-900mm 截面范围内的流场均匀区直径均保持在Φ882mm 以上,均匀区面积较开口试验段增加了约31.57%;马赫数6喷管采用闭口试验段时,均匀区面积比开口试验段仅增加了约8.24%,流场品质略为提高。超声速条件下,闭口试验段的流场均匀区增加明显;但在高超声速条件下,闭口试验段的流场均匀区增加比较有限。     

超声速化学反应流动的LU-SGS伪时间迭代空间推进求解

介绍了求解超声速多组分有限速率化学反应流动的伪时间迭代(lower-upper symmetric-Gauss-Seidel,LU-SGS)方法.空间推进求解多组分抛物化(parabolized Navier-Stokes,PNS)方程时,在一个推进面上采用修正的LU-SGS方法迭代至收敛,把得到的结果作为初值赋给下一个推进面.沿推进面依次迭代直至求解完整个流场.采用伪时间迭代LU-SGS方法求解化学反应PNS方程,计算结果的准确性和时间迭代求解完全Navier-Stokes(N-S)方程相当,求解效率提高一个数量级.      

风洞分流锥及孔板整流的数值模拟研究

采用数值模拟技术研究了2m×2m超声速风洞设计引导试验大开角扩散段配置一个中心倒锥和两层球冠状孔板的内流场,中心倒锥和两层球状孔板的不同组合共有五套。孔板有几百至上千个开孔,是模拟的难点,发展了一种孔板流动CFD边界条件模型。另一个难点在于风洞管道内流入口和出口边界条件的准确处理。计算表明,分流锥可将流动有效抑制在大开角洞壁附近,防止洞壁附近的扩张分离,但在分流锥底部将产生大尺度的分离涡环结构,经过孔板的整流后,总压有较大损失,但流线趋于平直,可达到预期效果。     

亚跨超音速喷管内流场的整体矢通量分裂算法

从N－S方程出发，采用先进的矢通量分裂算法，对亚跨超音速喷管内流场进行了整体模拟计算。文中对几种不同计算格式的计算效率进行了比较，给出了超音速占优混合流动的直接分割求解方法。计算采用代数湍流模型，跨音速段的计算结果同其它实验和计算结果进行了比较。计算中发现对出口边界条件进行亚音速、超音速和有回流分别处理很有必要；边界层内的M数沿壁面法向变化非常快；不同壁温条件和不同的流动Re数均会改变边界层的厚度。     

尾喷管内外超声速流场数值模拟

采用有限元法数值模拟带二次流可调收－扩喷管内外跨声速和超声速流场，在计算中采用四边形等参元确定位移函数，利用加权余量法中的Galerkin法建立有限元方程，数值研究三种型式喷管和四种飞行工况下尾喷管内外流场。计算结果与实验数据符合较好，说明可用此数值模拟技术帮助开展实验研究，进行尾喷管优化设计。     

向量化多重网格的Mac Cormack显格式解轴对称底部干扰流场

本文将多重网格法与MacCormack显格式的算法相结合,计算了轴对称底部流场中斜向喷流与超声速外流粘性干扰流场,提高了计算效率。流场的分析表明分辨率较好,计算合理。在YH-1上将程序向量化,大大节省了计算时间。     

三维任意机身的跨音速全位势流计算

 <正> 准确、有效地估计旋成体或任意截面机身在有迎角、有侧滑情形下的跨音速绕流特性,对飞机机身、导弹等外形设计有重要意义。已趋成熟的全位势流计算方法在翼型、单独机翼及翼身组合体绕流中得到广泛应用,并产生了诸如Holst的TWINGn和TAIR,Jameson的FLO系列等实用程序;但对旋成体有迎角绕流或单独机身的全位势绕流计算在文献中还很少见到。现采用守恒型全位势方程、代数法生成的贴体坐标网格和Holst的AF2高效差分方法对旋成体或单独机身有迎角、有侧滑时的跨音速绕流进行研究;并给出分析机身跨音速特性的计算程序是很有必要的。     

机动飞行器中等攻角高超声速无粘绕流数值模拟

本文采用ＮＮＤ差分格式，利用推进－迭代方法，数值模拟了机动飞行器的中等攻角高超声速无粘绕流流场，发展了一种适用于处理中等攻下背风面流场的技术，大量计算表明，结果准确、可靠。气动力力系数和壁面压力分布与实验数据一致，流场波系结构模拟正确。     

用四阶Runge-Kuta法和TVD格式求解两维跨超音流场

利用坐标转换给出两维任意曲线坐标系下雅可比矩阵特征值、特征向量及Harten型TVD格式的详细表达式, 并结合四阶Runge-Kuta法发展了一种新的、有效的求解超音、跨音流场的算法。数值计算表明, 本文的方法可有效地求解跨音流场。      

直接的矩阵分裂计算跨音速喷管流场

给出了分离系数矩阵求解欧拉方程组时系数矩阵的直接分离公式，减少和简化了应用该方法时的矩阵运算，参照Ｂｅａｍ－Ｗａｒｍｉｎｇ显式格式构造与ＳＣＭ方法一致的内点差分格式，并分析了它的数值特性，对有斜激波的跨音速内流场进行了数值仿真。计算中采用变系数的当地时间步长，加快了流场达到的收敛速度，计算无需引入任何人工参数，计算值与实验结果吻合。     

燃气源开槽导管三维流场的数值模拟

为了解某火箭弹控制系统燃气源开槽导管内的气流流动特性,通过ＳＩＭＰＬＥ方法对其复杂流场进行了数值模拟。采用隐式时间分裂格式、幂函数方案及交错网格修正压力等方法,获得了燃烧室后部和导管内气流的速度、压力、密度和温度的三维分布,计算结果与实验现象的规律一致     

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

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