一种同位网格上的迎风有限元方法

通过Ｇｒｅｅｎ定理将对流项变量从微分算子中分离出来,从插值函数入手引入迎风格式,是对强对流问题有限元计算中对流项变量的一种新的处理方法。按这种方法采用局部斜迎风格式及速－压同位网格公式,构成了一种对高Ｒｅｙｎｏｌｄｓ数流体流动数值模拟比较有效的有限元方法。     

带座舱飞船高超声速再入非平衡流场的数值研究

本文针对带座舱飞船高超声速再入大气层过程中存在的严重气动加热现象 ,利用简单隐式TVD差分格式和激波捕捉法 ,数值求解三维化学非平衡Navier Stokes方程 ,其中化学模型是 7组分 1 5个化学反应的空气化学模型 ,对带座舱飞船再入高度为 70km和速度为6 7km/s的化学非平衡流场进行了数值模拟 ,给出了带座舱飞船再入迎角为 0°,1 0°和 2 0°情况下的壁面热流、表面压力和电子数密度等参数分布     

双曲型守恒律方程的一个大时间步长二阶TVD差分格式

本文根据Harten,A.的大时间步长差分分格式构造思想,为双曲型守恒律方程弱解计算构造了一个2K 3点大时间步长二阶显式差分格式——LTS-LF格式,得到了其在CFL限制K下为总变差不增差分格式(TVD格式)。文章按照Roe的方法推广格式到方程组情形,并就Burger’s方程和Euler方程组黎曼问题进行数值试验,结果令人满意。     

一种改进的MUSCL格式

给出了一种改进的二阶精度的ＭＵＳＣＬ格式，并用来求解Ｅｕｌｅｒ方程，改进后的ＭＵＳＣＬ格式在线化插值步中，用Ｖａｎ Ｌｅｅｒ的插值斜率对特片变量进行插值，并对与熵波、剪切波对应的特片变量使用了小耗散的插值斜率，在迎风步采用Ｒｏｅ的近似黎曼解求得迎风通量。比较一维和二维的计算结果均可以看出，改进后的ＭＵＳＣＬ格式在消除数值振荡的同时，明显提高了流场间断的分辨率。     

AUSM+方法在黏性内流模拟中的应用

简要阐述了三维雷诺平均Navier-Stokes(RANS)方程求解程序中使用的对流迎风分裂格式改进形式(AUSM+)以及Spalart-Allmaras(S-A)一方程湍流模型,然后对Sajben收缩扩散通道内跨声速流动和高速透平静叶排内部流动进行了数值模拟,并将数值结果与实验结果进行了比较,结果表明自主开发的三维R...     

径向基函数在三维Euler方程数值计算中的应用

基于三维Euler方程,对非结构四面体网格给出了一种基于紧支径向基函数重构的ENO型有限体积格式,方法的主要思想是先对每一个四面体单元构造插值径向基函数,而在计算交界面的流通量采用高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法。最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力。     

燃气自由射流的正格式数值模拟

采用二阶正格式方法对复杂燃气自由射流进行了数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对不同欠膨胀压力比下的燃气射流进行了数值计算。计算结果表明,该方法能够较好捕捉到波胞、滑移线、射流边界以及三波点等复杂射流流场波系结构,与实验照片反映的流动特征及高精度、高分辨率的三阶间断有限元方法的计算结果吻合较好,马赫盘位置与理论预估值和实验测量结果误差较小。证明该方法对激波具有较强的捕捉能力,为此类流动的数值研究提供了一种新手段。     

可压缩混合层中的涡结构和激波

采用高精度高分辨率格式进行直接数值模拟,既得到了对流马赫数为1.2的混合层中线性发展之后的λ涡、发卡涡和涡的破碎,还得到了破碎后发展的多模态主导下的分层涡结构,以及对应动量厚度饱和下的涡结构。其流动结构的发展显示了混合层在自由空间充分发展能力的有限性。此外捕捉和显示了混合层中自然发展起来的外缘激波、对旋激波和单涡激波,激波主要是由于涡结构的高速旋转造成的相对超声速、低压区和逆压梯度区而产生的。激波的存在抑制了混合层从主流吸取质量和动量,也抑制了混合层内部流质混合的均匀性。     

非稳定导热问题有限元算法的守恒条件及一种守恒格式

本文对非稳定导热问题的有限元算法的能量守恒条件进行了讨论。认为有限元之解有时出现反常、振荡的原因是对问题时空域的离散化未能使相应格式满足能量守恒原理。依据文中的分析,提出了一种以单元边长做为集中热容阵权因子的守恒格式。验算表明,其稳定、可靠、精度较高。     

基于Van Leer+AUSM混合格式超声速流场的并行数值算法研究

基于CLUSTER集群并行体系结构与消息传递库MPI并行环境,开展了超声速流场的并行数值算 法研究。计算网格采用多块对接粘性结构网格,控制方程为积分形式的Navier\|Stoke s(N-S)方程,空间离散采用基于Van Leer+AUSM混合格式的有限体积法,粘性项采用中心 格式离散并利用格林定理计算粘性通量中的导数项,时间推进采用五步R-K法,应用一方程S -A模型计算湍流的粘性。并行计算采用 分区并行的思想,设计了满足负载平衡要求的区域分解和通讯方法,编程上采用单控制流多 数据流(SPMD)模型以及Master／Slave并行模式。最后,通过对某弹丸和某制导火箭弹超声 速流动的数值模拟, 验证了并行算法的准确性、高效性和处理复杂几何外形的能力。