高超声速三维化学非平衡流动的数值模拟

高温非平衡化学反应对高超声速再入飞行器绕流流场的气动热特性有较大影响.通过采用有限体积NND(Non-oscillatory, containing No free parameters and Dissipative scheme)格式求解三维N-S方程,分别计算了小攻角条件下球锥的完全气体粘性流场和化学反应非平衡气体流场.比较并分析了化学非平衡作用对三维高超声速钝头体绕流流场的影响.所得结果与国外文献结果符合的较好,证明了计算方法与程序的有效性和实用性.      

基于简化Boltzmann方程的超声速流稀薄效应分析

从近连续流到近自由分子流,计算分析超声速流场的稀薄效应。模型Boltzmann方程先后经简化分布函数和离散速度坐标变换后,采用一个隐式通量修正二阶迎风TVD格式差分求解。壁面取漫反射气体分子模型。在Knudsen数为0.001,0.01,0.1,1.0,10等5种情形下,数值模拟二维圆柱氩气体绕流,观察到了不同的波系及尾迹结构。总阻力系数计算值同实验数据基本吻合,碰撞项Shakov动力学模型反映的流场细节略好于修正BGK模型。稀薄效应趋于弱化激波等超声速流动结构,物体影响域增大。     

双时间步长加权ENO-强紧致高分辨率格式及在叶轮机械非定常流动中的应用

在非结构网格下,针对非定常三维N-S方程组发展了一种双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间上将r=3的加权ENO格式与强紧致格式相结合去处理N-S方程中的对流项以及离散方程的右端项,并用四阶精度的紧致格式去计算N-S方程中的粘性项。典型的3个算例从不同侧面对格式进行了考核。计算表明:该算法具有高效率与高分辨率的特征,所得的计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下模拟非定常流动的物理过程。      

不同半径圆柱诱导CH4/空气预混燃烧的计算研究

本文研究不同半径圆柱诱导CH4/空气预混燃烧流场。采用保自由流5阶WENO格式求解贴体坐标变换后的多组分Euler方程,用基元反应模型描述CH4/空气燃烧。不同于标准WENO格式通量构造方法,该WENO格式数值通量由方程的解进行WENO插值得到,在曲线坐标系下具有保自由流性质。首先给出了保自由流WENO格式精度和保自由流的数值验证,然后计算圆柱诱导CH4/空气预混气燃烧流场,并考察不同半径圆柱的影响,给出燃烧流场压力与温度等值线和云图、压力和温度沿过驻点线分布。结果表明:在考核计算结果网格无关性基础上,该WENO格式可准确地捕捉激波和火焰阵面形状、激波和火焰面驻点距离,得到的计算结果和文献结果相符。增大圆柱半径,激波和火焰面被推向来流方向,激波和火焰面之间距离也减小。和TVD格式相比,5阶WENO格式采用四分之一的网格数可得到近似相同的计算结果。     

高阶精度方法下的湍流生成项对跨声速流动数值模拟的影响研究

利用五阶空间离散精度的WCNS格式和多块结构网格技术,通过求解雷诺平均NS方程,开展了SST两方程模型不同湍流生成项组合方式对跨声速流动数值模拟影响的计算分析。研究的主要目的是为高阶精度格式在复杂外形上的工程应用提供技术支撑。计算模型采用了RAE2822超临界翼型和DLR-F6翼身组合体构型。研究内容主要包括不同湍流生成项对残差收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布以及模型整体气动力特性的影响。不同湍流生成项组合方式的流场计算结果还与风洞试验数据进行了对比。研究结果表明:对于小迎角不存在明显分离的跨声速流动,不同湍流生成项对流场的高精度计算结果的影响很小,可以不用考虑。     

三阶HWCNS的构造及其在高超声速流动中的应用

对网格质量要求高、计算稳定性差和计算效率低是制约高阶精度格式应用于高超声速复杂流动模拟的重要因素。针对这些问题,发展了三阶精度的混合节点半节点加权紧致非线性格式(HWCNS3),改进其光滑测试因子和非线性权得到了HWCNS3-OP,并给出了它们的频谱特性。利用Lax和Osher-Shu算例测试了格式对间断和高频波的捕捉能力;通过钝锥和航天飞机的高超声速绕流算例,考察了HWCNS3-OP在真实流动模拟中热流和气动力的预测精度及其计算效率。研究结果表明:HWCNS3-OP具有较高的分辨率和良好的间断捕捉能力,高频波捕捉能力相对HWCNS3提高了约3倍,相对守恒律的单调迎风中心格式(MUSCL)提高了约4倍;HWCNS3-OP计算稳定性较好,计算效率相对五阶HWCNS提高了2~3倍,HWCNS3-OP是一种较适合高超声速复杂流动模拟的高阶精度格式。     

消除粘性项高阶离散数值振荡的半结点-结点交错方法

通常的粘性项处理方法是连续两次采用结点型中心差分格式求得.但是,通过两次采用高阶结点型中心格式求得的粘性项易在流场中的"间断"附近产生数值振荡.本文采用多种格式,通过求解一维Burgers方程来充分展示了这种振荡现象.为了消除这种数值振荡,一种半结点型的高阶紧致格式被用来求解粘性项.Fourier频谱分析表明,这种格式具有非常优越的保频谱性能.一维和二维的数值计算结果表明,通过"半结点-结点"交错采用这种半结点型格式的方法可以非常有效地避免粘性项高阶离散可能导致的数值振荡.     

考虑粘性时激波反射及壁面热流计算

本文采用显式ＴＶＤ格式，求解薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，计算了激波与楔面相互干扰的流场情况，得出了马赫反射的各种等强度级及固壁上详细的物理分布，根据计算得出的流场，计算了壁面摩擦阻力及热流分布。     

一类新的二阶精度非交错形式差分格工工的构造

本文对二维标量非线性双曲型守恒律，以非交错形式上Ｌａｘ－Ｆｒｉｅｄｒｉｃｈｓ格式为基本模块，构造了一类新的不须解Ｒｉｅｍａｎｎ问题的差分格式，并证明其具有时空二阶精度。数值实验表明，格式具有较高的分辨激波，稀疏波及接触不连续的能力。     

一种MWENO谱体积（SVMWEN05）格式

构造一种新的谱体积（spectral volume）格式来求解双曲型守恒律，记为SVMWEN05（Spectral Volume method by Multi—Weighted Essentially Non-Oscillatory）格式。其主要思想是：第一步，将空间计算区域划分为一系列单元，称为谱体积，等分每个谱体积为一些子单元，称为控制体积（controlvolume）。第二步，在谱体积内部采用类似MWEN05（Multi-Weighted Essentially Non—Oscillatory）的格式进行重构，而在谱体积的边界处采用传统WEN05格式进行处理。第三步，利用Runge—KutmTVD离散方法对半离散格式进行时间离散，得到时空一致高精度全离散格式。最后．在文中给出几个经典数值算例用以验证本格式的计算能力。