强迫耗散非线性发展方程与完全平方守恒格式

从描述大气和海洋运动的强迫耗散非线性发展方程出发，对强迫耗散非线性大气和海洋方程组显式差分格式的计算稳定性进行了分析，构造了一类强迫耗散非线性发展方程的显式准安全平方守恒差分格式。理论分析和数值试验证明，这类显式准完全平方守恒差分格式是计算稳定的，值得推广应用。     

超声速开式空腔流动的数值模拟

研究了不同长深比、不同雷诺数条件下,超声速开式空腔流动特性的变化。算法采用空间上5阶的WENO格式,时间上具有TVD性质的3阶Runge-Kutta方法,湍流模型选用B－L模型。数值模拟获得了不同的时刻的注场图像,压力波、涡以及空腔结构之间的相互耦合作用,使得空腔内存在典型的周期性压力脉动,但不同长深比空腔的流动性存在较大差别。     

二维非结构网格的一个TVD型有限体积方法

考虑双曲型守恒律方程,对二维非结构三角形网格给出一种 TVD型有限体积方法,主要思想是在一阶单调格式的基础上,在每一个单元上对变量作单调线性重构函数,时间离散采用二阶 Runge-Kutta方法。通过计算分析了该方法的精度,对平面激波反射和空穴流动的计算结果表明该方法是成功的。     

三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

NND格式在多维理想磁流体方程组中的应用

采用修正的四步Runge-Kutta方法求解三维一般曲线坐标系下的理想磁流体方程组,为克服数值振荡,加特征型NND格式进行后处理.特征型NND格式推广到求解三维磁流体(MHD)问题需要知道雅可比通量的左右特征矩阵,在具体计算时需要克服矩阵的奇性.本文用三维程序采用推广的特征NND格式计算了一维MHD激波管和二维(MHD)喷管流动,计算结果表明,特征NND格式保持了TVD格式高精度的优点,又具有计算简单的特点,在包括强弱间断等复杂波系的定常和非定常MHD流场数值模拟中是成功的.     

三维复杂外形横截面拓扑结构随迎角变化研究

通过数值求解抛物化的椭圆型方程生成复杂外形三维网格;采用交替方向隐式分解的隐式NND格式求解全N S方程模拟"类升力体"外形在高超声速下的大迎角流动;采用张涵信发展的流动拓扑结构理论分析给出了"类升力体"外形垂直于体轴的横截面流线随迎角变化的拓扑结构;此外,通过计算发现:大于20°迎角后,在部分横截面背风对称线上出现结构不稳定的鞍点相连现象。     

分区对接网格算法的应用研究

本文应用分区对接网格技术，由于采用了对接面上网格点完全对应的方法，能自动满足通量守恒和计算的稳定性。采用有限体积算法下的原始变量型NND格式，求解NS方程，数值模拟了复杂外形的高超声速层流流场，计算结果与实验测量符合较好。     

一种通量加权型紧致格式

将通量加权的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的加权型差分格式.利用该格式与二阶TVD格式分别计算了一维方波、组合波问题和一些黎曼问题,如Sod问题和Shu问题,以及一维定常激波问题.计算结果的比较表明加权格式无论在捕捉各种间断,还是在分辨各种复杂波系上,都具有较大的优势,并与精确解非常吻合.      

高超声速热化学非平衡钝体绕流数值模拟

对高超声速热化学非平衡钝体绕流进行研究,发展了高温热化学非平衡条件下的高超声速气动热数值模拟方法。采用二阶精度迎风TVD格式的N-S方程有限体积法多块结构网格求解器,考虑了空气5组分及化学非平衡效应的影响,对组分生成源项及振动能源项采用了点隐式的处理,提高了计算效率。对球锥模型和RAM-CII飞行实验等算例进行了数值气动热模拟,分析了不同非平衡模型在热化学非平衡条件下流场的影响。计算结果表明建立的数值模拟方法具有较高的精度。     

可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展.但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题.界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因.为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案...