双重加权实质无波动激波捕捉格式的改进和应用

本文重新研究了双重加权实质无波动(DWENO)激波捕捉格式的构造过程，对模板光滑系数的定义，以及两重权值的计算方法进行了改进，并通过一维Riemann问题测试了这些改进，确定了参数的选择范围。最后，本文将改进后的DWENO(r=3)格式进一步应用于一维、二维典型问题的数值模拟，并比较了DWENO格式和WENO格式的精度和分辨率。     

高超声速飞行器外形热流密度分布计算的高精度方法研究

为提高高超声速条件下飞行器表面热环境的计算精度,发展了高精度算法WCNS-E-5,并开展了热流密度分布的数值计算研究.该算法采用的是五阶精度加权显式非线性格式,结合四阶精度的二阶偏导数差分近似、四阶精度的边界格式,同时还对网格偏导数及温度梯度也离散为四阶精度.通过求解三维Navier-Stokes方程,数值研究了高超声速飞行器外形如钝锥和双椭球体的物面热流密度分布.模拟结果表明,这种保证全流场高精度的WCNS-E-5得到的流场图像清晰、真实.与通常的二阶算法MUSCL对比,WCNS-E-5具有更高的流线分辨率,获得的热流密度更接近于实验测得分布.     

高阶精度线性耗散紧致格式的渐近稳定性

文献〔２〕提出的高阶精度线性耗散紧致格式,（ＤＣＳ）能较好地抑制高波数振荡。本文进一步研究表明：内点耗散参数α能够控制耗散紧致格式ＤＣＳ３和ＤＣＳ５的渐近稳定性,边界格式参数α１和α对以上格式的渐近稳定性的影响是次要的,并给出了 些参数的取值范围。通过平面可压缩Ｃｏｕｅｔｔｅ流的稳定性数值模拟表明,当格式ＤＣＳ５满足渐近稳定性要求且具有耗散时,计算结果与线性稳定性分析的结果十分吻合,否则,会出现较     

DCS5在计算气动声学中的应用研究

对流项采用五阶线性耗散紧致格式(DCS5)离散,粘性项采用6阶中心紧致格式离散,时间项采用5级5阶龙格-库塔时间推进格式求解流体动力学方程组,并由此建立近声场的直接模拟方法。为了应用结构网格在复杂几何形状条件下对流场进行计算,发展了高精度的多块特征对接边界技术。完成了低速无粘振荡圆柱和粘性圆柱绕流等标准计算声学算例,获得了令人满意的结果。     

两种湍流模型及其可压缩修正在高超压缩拐角中的评估(英文)

采用在高阶精度差分格式对SA湍流模型和SST湍流模型及相应的可压缩修正模型在高超拐角流中进行了评估。可压缩修正方式考虑了密度梯度、压力膨胀和湍流马赫数等方法。为了减小数值误差与模型误差之间的混淆,控制方程的对流项采用了5阶精度的加权紧致非线性格式(WCNS-E-5),粘性项采用了一种半结点/结点交错的4阶中心格式。通过对马赫数为9.22的15度拐角和34度拐角湍流的模拟,考察了原始湍流模型及其修正模型的效果。计算表明:原始SST模型对高超拐角湍流的预测比原始SA模型准确,这种准确主要体现在对分离区预测、再附点附近压力和热流峰值预测上。通过混合采用Catris和Shur等的方法对 SA模型进行可压缩修正可以大大改进模拟效果。在SST模型的可压缩修正方法中,Catris的修正方法最好;考虑压力膨胀修正后得到的分离区远远偏大,比本原始SST模型更差;考虑湍流马赫数的修正方法得到的分离区偏小。本文还给了部分低阶格式的计算结果,高阶格式与低阶格式相比,对分离区大小、再附点附近的压力和热流峰值等的预测准度有所改进。     

辉光放电等离子体对边界层流动控制的机理研究

本文求解电位势方程得到电场分布,假定电场强度大于击穿阈值的区域为等离子体区,给定电荷密度,得到了作用于流体上的电场力。通过求解带源项的NS方程,研究了一个大气压下的均匀辉光放电等离子体对边界层流动的影响,考察了电场力做功对流动的影响。本文研究结果同文献[6]一致,即电场力总体上使边界层流动加速。另外,电场力做功对流动参数的影响可以忽略。文献[6]给出的线化电场模型和本文得到的电场相比,具有较大的差别,该差别引起了流动参数显著差别。