改进的二维三阶半离散中心迎风格式

对二维三阶半离散中心迎风格式中的权函数给出了简化改进. 在保持格式精度的基础上, 改进后的权函数在二维情况下具有更加简单直接的结构而且严格非负. 该改进方法得到的格式仍然具有半离散中心迎风格式的优点, 同时保持了重构函数的非振性. 时间离散采用保持强稳定性的三阶Runge-Kutta方法, 并利用四阶Lax-Wendroff (L-W) 格式计算磁流体算例中的磁场散度. 用该修正格式计算了二维磁流体数值算例, 得到高精度无振荡的结果, 验证了此方法的有效性.      

超声速进气道可压及不可压流动数值模拟

采用Chio-Merkle预处理矩阵对可压NS(Navier-Stokes)方程组进行时间预处理,分析了预处理方法的物理和数学背景.用有限体积方法,结合LU-SGS(Lower-Upper-Symmetric-Gauss-Seidel)隐式时间积分和AUSM+(P)(Precondition Advection Upstream Splitting Method)格式、中心差分2种空间离散格式,求解该预处理NS方程组.通过圆弧凸包流动和二维方腔顶盖驱动流动的数值试验表明,该方法克服了传统时间迭代方法模拟低速流动时的刚性问题,加速了收敛过程,可以同时有效地模拟可压及不可压流动.用此方法与块结构化网格技术相结合,进行了在不同飞行马赫数、攻角和出口反压条件下,混压式轴对称超声速进气道三维流场数值模拟.计算结果表明:该方法能准确地捕获复杂激波系,清楚地揭示进气道三维流动现象;进气道性能随工作条件的变化规律,与理论分析一致.      

广义拉格朗日乘子MHD方程组的两类逆风分裂格式

开发高性能的磁流体力学数值模拟方法是提高空间天气数值预报研究的一个重要方面.有限体积法的逆风分裂格式具有良好的间断捕获能力，Steger-Warming和AUSM（Advection Upstream Splitting Method）是逆风分裂格式FVS（Flux Vector Splitting）方法中具有代表性的两种格式.采用这两种格式求解具有伽利略不变性的扩展型广义拉格朗日乘子磁流体力学（EGLM-MHD）方程组，对Orszag-Tang涡流问题和三维爆炸波问题进行数值模拟，结果表明两种格式均能得到稳定精确的数值结果.与Steger-Warming格式相比，AUSM格式产生的磁场散度误差更小，计算速度更快.      

交替LU分裂算法及其在CFD中的应用

在CFD(Computational Fluid Dynamics)时间相关算法中,为了保证计算的稳定性,时间步长的取值通常会很小,这将导致计算过程收敛缓慢.针对这一问题,提出了一种新的迭代算法—交替LU分裂(ALUS,Alternating Lower-Upper Splitting)算法,可以有效加速收敛,提高计算效率.ALUS算法将系数矩阵分裂成上、下三角矩阵,因此仅需要利用追赶法求解两个三角矩阵,计算量较小,容易实现.给出了ALUS算法收敛的定理,并且通过线性问题以及CFD圆柱绕流的数值模拟对ALUS算法进行了检验.理论分析和数值实验的结果均表明:ALUS算法计算量小,大大节省了计算时间,而且该算法是鲁棒的.因此ALUS算法是高效的、稳定的算法,适用于CFD数值模拟.     

数值模拟高超音速化学非平衡流的隐式格式

针对高超音速化学非平衡流提出了两种隐式格式.一种是近似因子分解算法对角化形式在隐式TVD(全变差递减)格式上的推广,另一种可视为隐式TVD格式的LU分解形式.数值试验表明,本文给出的算法保持了TVD格式的高分辨率特性,同时计算量大大减少,是高超音速化学非平衡流数值模拟的有效方法.     

高超音速化学非平衡流的隐式格式迭代解法

针对高超音速化学非平衡流计算中的隐式TVD格式给出了一种迭代解法.该方法是在原有的Harten-Yee格式基础上,在每一时间步对无粘部分进行迭代运算,目的是为了加速收敛.借助对模型方程的收敛性分析,证明该迭代算法是收敛的.二维数值实验表明该方法对于高马赫数时的化学非平衡流能明显提高隐式格式的计算效率,其计算结果是满意的.     

高阶熵条件格式下Euler方程与 N-S方程的混合算法

针对在工程问题中使用N-S方程计算量大的问题,提出只在粘性作用显著的局部区域求解N-S方程,在流场的其余部分采用Euler方程计算的办法来模拟具有复杂外形的流动问题,并提出了一种基于高阶熵条件格式的算子分裂算法.应用所构造的算法对绕迫击炮弹的亚、跨超声速流动进行了模拟计算.计算结果与试验数据的对比表明,所提出的方法是切实可行的.      

多单元直排塞式喷管发动机性能

为了了解优化设计的塞式喷管的性能及燃气流动中热力学参数变化对性能的影响,对比钟型喷管与塞式喷管的高度特性,从曲线坐标下的三维平均雷诺N-S方程和Euler方程出发,采用LU时间隐式格式、MUSCL空间离散方法,发展了模拟塞式喷管三维流场的数值程序.计算了从喉部圆转方内喷管的性能,比较了冻结和变化热力学参数对塞式喷管性能的影响及塞式喷管与相同面积比钟型喷管的高度特性曲线.计算结果可为塞式喷管的设计研制提供参考.      

三维抛物化N-S方程全隐格式矢通量分裂求解

采用全隐格式的矢通量分裂算法,在通用曲线坐标系内求解三维无因次抛物化强守恒型N-S(Navier-Stokes)方程.所建立的计算程序,可计算较为复杂边界的三维流场,如求解亚、跨、超音速共存的流场,分区模拟内外流共存流动及其混合剪切尾流.给出了一轴对称亚跨超音速喷管内流场的计算结果,包括计算网格,流矢分布,压力、温度和马赫数等值线分布.也给出了一喷管内外流及其混合尾流的计算结果,包括计算网格,流矢分布,无因次压力、密度和温度等值线分布.     

绕AOTV高超音速流场的数值计算

本文采用实质上为二阶精度的隐式无波动、无自由参数耗散(NND)差分格式结合Steger-Warming矢通量分裂方法求解N-S方程,对AOTV外型的高超音速绕流进行了模拟。在计算中,对脱体激波采用了捕获法。作为验证,首先计算了二维圆柱绕流,并同文献[1]的数据进行了比较,取得了满意的结果,算例的计算结果还表明,采用本文的方法对激波的模拟具有较高的分辨能力。     

层流平板摩擦阻力的数值计算

为提高层流摩擦阻力数值计算的精度,采用数值求解N-S方程的方法,对层流平板摩擦阻力进行了计算.计算中采用4种先进上风格式,同时配合以14种不同的壁面法向网格分布和5种常用限制器,并将计算结果与布拉修斯解进行对比.结果表明,van Leer的FVS(Flux Vector Splitting)格式数值耗散大,对网格分布和限制器敏感,不适合用于摩擦阻力计算.限制器和计算格式的粘性分辨率越高,摩擦阻力的计算精度越高.在计算量允许的情况下,可以减小壁面法向第1层网格高度,以提高计算精度,但应同时兼顾网格长宽比.对于摩擦阻力计算,边界层内应至少布置4个网格点.      

快速大时间步长熵条件格式的分辨率研究

提出了高分辨率快速大时间步长熵条件格式的构造方法.用激波管问题对一族熵条件格式进行研究.在精度、步长、限制器方面进行了详细的数值实验,研究了同样计算量下各种格式的表现品质.从理论上保证了大时间步长格式的无振荡性质,从具体的数值实验分析中确定了大时间步长格式的分辨率问题.      

非定常N－S方程的线化及其差分求解

在微幅简谐振荡的假设下，对计算平面内非定常Ｎ－Ｓ方程进行了线比，并在ＭａｃＣｏｒｍａｒｋ于１９８０年提出的半隐格式的基础上，发展了一种适用于求解计算平面内线化Ｎ－Ｓ方程的差分格式．用该差分格式求解计算平面内的线化Ｎ－Ｓ方程时，可以避免求解五对角方程；且时间步长可以取得较大，计算效率高．     

基于RoeM格式思想的激波稳定格式构造

Roe格式有着较高的黏性分辨率以及间断分辨率,但它在实际流动模拟时很容易出现激波不稳定现象.此外,它在无黏定常计算时无法保证总焓守恒.为了在保留Roe格式分辨率高的优点基础上改善其上述缺陷,将RoeM格式与高阶激波判别方法相结合,提出了RoeMW1,RoeMW2格式.计算结果表明,在流动较为简单的情况下提出的RoeMW1, RoeMW2 格式可以有效避免激波不稳定现象的出现且保持总焓守恒.而在一些较为复杂的流动结构计算中,虽然比之于RoeM格式均有所改善,但在RoeMW1基础上改进的RoeMW2格式在间断模拟时鲁棒性更强.