非对称体超音速无粘绕流流场数值计算

 本文给出了非对称再入飞行器的超音速三维无粘流场的数值计算方法。它归结为求解三个自变量—阶定常拟线性双曲型偏微分方程组的初边值问题。采用MacCormack二步显式差分格式。 为了说明计算方法的实际能力,对弯头钝锥和带控制翼的再入体作了计算。其精确度是令人满意的。     

斜切反喷管性能计算

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

超声速翼—身组合体Euler方程数值解

本文采用MacCormack二步显式差分格式,用空向推进法给出了超声速翼-身组合体Euler方程的数值解。应用薄翼假设生成简易的多区网格。采用Kentzer方案处理边界。计算中通过对横流速度型的修正改善了背风面压力分布特性,解决了推进中断问题,同时模拟出背风涡。     

绕弯曲头锥超音速流动差分法数值计算

本文给出了绕弯曲头锥理想气体超音速流动差分数值计算方法。这是一个三维绕流问题,由于弯曲头锥的外形特点,头锥部分的数值计算采用了计算轴移动法。数值计算采用MacCormack二阶精度二步格式,物面上用一组由特征相容条件导得的方程。本方法已经用BCY算法语言编制成计算程序,在109丙电子计算机上进行运算。本文最后给出计算实例。     

高超声速无粘流场的推进计算

本文将文[1]中用于时间相关法计算的NND格式推广到定常超声速流动的空间推进计算,采用二步的预测、校正方法保证了推进方向的二阶精度,可以证明,这种二阶精度的NND格式具有TVD性质,是MacCormack二步显式格式的推广。本文首先将格式应用于二维平板上斜激波反射流场的推进计算,以检验格式捕捉激波的能力,同时研究了不同的通量分裂方法对格式捕捉激波能力的影响,得到了相当满意的结果。在此基础上,计算了航天飞机简化外形的身部超声速流场,给出了M_∞=10,α=0°,和M_∞=5,α=5°两种状态的部分结果,计算结果清楚地描绘了由于气流在机翼附近受到强烈压缩而产生的内嵌激波与外激波相交的复杂流场结构,与文[7]相比,流场结构更为清晰。     

介绍一种适合于计算跨音速流动的混合差分格式

本文介绍MacCormack格式和迎风格式的一种混合差分格式。它适用于跨音速流计算。本文用它计算了绕圆弧翼型的跨音速流场,并详述了处理方法和结果。     

一种新的用于MHD方程逆风格式的Jacobian矩阵分裂方法

针对理想MHD方程,提出了一种新的基于MacCormack算法的雅可比矩阵分裂方法,克服了原有方法稳定性差的问题,并成功地应用于理想MHD方程的求解.控制方程在非结构混合网格上进行空间离散,其中对流项采用本文发展的逆风向量分裂格式,并引入了双曲型磁场散度清除技术,时间推进为显式5步龙格-库塔方法.对MHD激波管流动和带均匀磁场干扰的二维高超声速钝头体绕流流场进行了数值模拟,得到了与参考文献相吻合的数值结果,表明本文发展的数值分裂方法可以有效地捕捉MHD流场的流动特征,并且具有比MacCormack方法更高的稳定性和计算精度.     

关于几个二阶显式差分格式的评述

本文着重讨论了MacCormack格式计算头部混合流场时出现的一些问题,当用不同的边界点算法与之配合时对计算稳定性的影响以及格式的稳定性在计算过程中的表现形式,比较了MacCormack格式与其他二阶格式的计算结果,在相同的边界算法的情况下各种格式的效能。本文也提出了一种新的物面边界的算法,其中用了特征概念,但又避免了传统特征线法的迭代过程。     

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

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

复杂前机身—进气道组合体的流场计算

本文提出了一种能够对复杂前机身-进气道组合体进行流场计算与气动分析的方法。通过求解Euler方程来模拟前机身-进气道组合体在亚、跨、超声速,不同攻角和发动机流量条件下的各类流场。     

软式加油过程中软管鞭打现象数值模拟研究

对空中加油过程中软管锥套系统产生的鞭打现象进行了数值模拟.流场计算中的空间离散采用Osher格式,紊流模型采用S-A一方程模型;运用“刚杆-球铰”模型离散方法对软管锥套系统建模,建立了软管多体系统模型并推导其运动控制方程;采用一类基于相邻单元搜索算法,通过搜寻软管节点所在加油机尾流场网格的宿主单元,建立软管微段气动力与流场的联系.运用方法对软管锥套系统的动态特性进行了仿真,结果与国内外飞行试验结果相吻合.     

基于MBD-CFD的软管-锥套空中加油仿真框架

随着无人加油等新兴应用场景的出现，软管-锥套式空中加油技术近年来重新得到了航空工业界的重点关注。使用计算机数值模拟手段对空中加油系统进行流体-动力学仿真，已成为空中加油系统设计、研发和飞行测试评估中的必备手段。一方面，基于工程简化模型的软管-锥套仿真方法的数值模拟精度难以满足实际需求，另一方面，完全基于非定常计算流体力学（CFD）的耦合仿真方法计算量较大，在有限计算资源条件下难以满足快速设计迭代与多参数优化的计算效率需求。通过将多个求解层次的流固耦合分析工具有机结合，建立了基于多体动力学（MBD）和CFD的加油机软管-锥套装置的耦合仿真框架和方法。对典型加油机-软管-锥套组合体的全尺寸构型和缩比构型进行了数值模拟，通过多个算例验证了本方法的可信性。通过将非定常CFD方法与基于气动建模的尾迹流场动态插值方法结合，可将后者的计算精度提高至与非定常CFD方法相当，同时将计算量减小2个数量级以上，从而实现了计算精度和效率的兼顾。     

超音速钝体无粘绕流的计算

研究用三维特征线相容性关系的迎风差分格式,数值模拟超音速和高超音速钝体无粘绕流流场。用钝体反方法计算头激波及其波后的流场,用轴对称流特征线法确定初值面上的流动参数。算例的计算结果与实验数据比较吻合     

MacCormack格式在计算非定常激波时的特点

前言 MacCormack格式(以下简称为MC格式)近年来在国内外用得较多。在非线性气体力学问题的计算中,它与著名的L-W格式或L-W二步格式相比有许多不同。其主要差别是MC格式的数值解中激波后的过头(Overshoot)较小。对于多维问题,MC格式对时间变量分步(Time Splitting)是十分有效的。 为了更好地了解MC格式,本文讨论MC格式的一些性质,并计算了两个典型问题。     

入射激波与飞行器头部激波干扰流场的计算方法

本文论述了入射激波与飞行器头部激波相交后的二维或轴对称干扰流场的计算方法。它是文[3]建立的具有最小耗散系数、保持解的单调性、一阶与二阶自动调节的混合格式的推广。通过对绕楔、球头和锥三种干扰流场的计算表明,本文的方法是可靠的,具有较高的精度,可用于多重间断问题的流场计算。     

声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析。研究了Steger-Warming通矢量分裂、VanLeer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

高超声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析，研究了Steger-Warming通矢量分裂、Van Leer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

超声速平面钝体流动的快速解法

本文介绍了一种平面流动的快速欧拉方程解法。该方法将原参数非定常欧拉方程组重新组合成以广义称曼变量表示的欧拉方程组，再用二点二步迎风格式离散解。针对钝体流动，本文先建立了动网格下的方程，构造了动网格下的算法。提出了一种简单的激波处理方法。计算结果表明，该方法速度快，稳定性好，对初场不敏感。     

固体火箭发动机一维非定常两相喷管流场计算

本文采用MacCormack预校二步差分格式对一维非定常两相流动喷管进行了数值求解,得到了固体火箭发动机压力建立过程中气相和凝相颗粒参数的变化规律.     

钝体头部“逾越”型转捩的直接数值模拟

通过直接数值模拟,研究了钝体头部边界层“逾越”型转捩的问题。在计算上,利用高阶紧致差分格式和高分辨率的无反射边界条件,求解全N-S方程,并将局部有限幅值扰动直接引入到钝体头部的曲面边界层中。摈弃了目前国际上通常采用的周期边界条件,以及以平板代替曲面的做法。结果表明:流场转捩的发生与一种类似“发卡”状的拉伸结构出现有关,而国外的计算结果中并未发现这一结构。此外,流场的失稳过程和物理图像也存在较大差异。这说明国外目前的做法和结论是值得商榷的。同时,模拟的结果也说明,这种方法对研究大扰动转捩来说是一条行之有效的途径。