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

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

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

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

涡喷发动机可靠性增长试验方案的设计

根据对产品可靠性的最低可接收值与幂律模型的可行性增长管理公式、可行性增长管理策略，推导了涡喷发动机可靠性增长试验所需的台数公式，参照某些涡喷发动机的试验数据及上述公式，给出了某B－2型涡喷发动机的可靠性增长试验方案及其理想增长曲线。     

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

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

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

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

混合气体流动的一种数值模拟方法

本文建立一个数值模拟完全气体混合流动的理论模型.该模型首先应用混合气体的Euler方程和每种气体组分的质量分数方程来控制流动.为了消除混合网格内气体组分界面附近出现的非物理振荡,我们假定混合气体的每种组分达到了动力学平衡状态然而尚未达到热力学平衡状态.这种思想导致需要另外给定每种气体组分的总能量方程.为使用高分辨格式来求解这组双曲型偏微分方程并且简化对所需要的Jacobi矩阵的推导,混合气体的压力方程也被耦合起来.Godunov型的波传播方法被采用来离散求解所获得的控制方程.从典型算例结果来看,一维问题的数值解与精确解一致,二维问题的数值解与理论分析一致.这说明本文的理论模型是合理的.     

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

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

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

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

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

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

超声速平头圆柱绕流DSMC/EPSM混合算法研究

本文将EPSM算法与DSMC方法结合，构造了可模拟含近连续流区及过渡流区的DSMC／EPSM混合算法。运用混合算法模拟马赫数等于5．37时超声速平头圆柱绕流，并与DSMC结果进行比较，证明了DSMC／EPSM混合算法的有效性，并对两种算法的计算效率进行了比较。