高精度ENO格式在射流数值模拟中的应用

本文分析了ENO格式的特点，并应用于全NS方程的迁移项和压力顶，模拟了轴对称和三维射流。首先对轴对称自由射流为算例比较了高阶TVD格式和ENO格式，证实了ENO格式捕获复杂激波波系的能力；其次计算了四喷管湍射流干扰流场，给出了相应的计算结果；最后计算了四股射流流入圆管内的流场，给出了物理量沿壁面和截面的分布。     

高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

CPR方法与高分辨率二阶格式的混合算法

重构修正方法（correction procedure via reconstruction,CPR）具有紧致高效的优点,但对较强激波的捕捉能力还相对较弱,而加权紧致非线性格式（weighted compact nonlinear scheme,WCNS）具有很强的激波捕捉能力。将基于高阶WCNS插值的二阶格式引入到高阶CPR方法中,构造了一种高效高分辨率的混合激波捕捉格式。首先,基于非线性权偏离线性权的程度的激波侦测器侦测出问题单元,并在问题单元附近引入缓冲单元,其余单元则标记为光滑单元。然后,针对问题单元和缓冲单元采用二阶格式计算,光滑单元采用CPR方法计算,构造混合格式。通过对等熵涡问题、含激波的问题以及激波旋涡干扰问题的数值模拟,测试了混合格式的精度、激波捕捉能力和计算效率。数值模拟结果充分说明了该混合格式具有很强的激波捕捉能力,同时在光滑区具有高分辨特性,可以应用于高超声速流动问题的高效数值模拟中。相比于基于高阶WCNS插值的二阶格式,此格式具有更高的计算效率和更高的分辨率。     

MLP高阶重构格式应用

在有限体积框架下,利用MLP(multi-dimensional limiting process)系列重构格式结合HLLHLLC(Harten-Lax-van Leer with contact)近似黎曼求解器,同时引入激波探测函数进一步降低MLP在光滑流动区域的数值耗散,数值模拟了超声速前台阶流动、结冰翼型的非定常流动、高超声速双楔流动和DLR F6-WB跨声速流场,研究了MLP系列格式在可压缩复杂流场中的表现.结果表明:在多维空间中,MLP格式能够在如强斜激波与网格线不重合等复杂流场数值模拟中保持严格的流场单调性;具有和传统MUSCL(monotone upstream-centered schemes for conservation laws)格式类似的计算效率,可以实现5阶,甚至更高阶重构;数值耗散更低,捕获更准确的激波位置,对航空工程数值模拟具有重要意义.     

用隐式WENO格式计算悬停旋翼跨声速流场

发展了一种基于高阶迎风格式的计算悬停旋翼跨声速流场的隐式有限体积法。对流项采用Roe通量差分分裂格式,使用五阶WENO格式进行左右状态重构,并与MUSCL插值进行比较;粘性项采用中心有限体积法。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。数值模拟采用了一种能够有效传递网格间流场信息的重叠网格,其中使用了三层内边界和贡献边界的方法以便插值的直接进行。用该方法对一跨声速悬停旋翼粘性流场进行了数值计算,数值结果表明:与MUSCL格式相比,WENO格式对激波位置捕捉得更准确,具有更强的涡捕捉能力,同时还表明了WENO格式在很大程度上能够克服涡耗散问题。     

激波绕流与激波碰撞的数值模拟

本文应用基于Poe平均构造的最小耗散的扩散与反扩散守恒差分格式，数值模拟单激波绕流三我形障碍物，及等强度激波碰撞和不等强度激波碰撞的二维非定常复杂流场。计算结果表明，该格式具有很高的单波精度及较好的激波与接触间断分辨度，可直接用于复杂流场的数值模拟。     

基于无网格算法的化学非平衡流数值模拟

为了验证AUSMPW+格式在无网格算法的精度等以及无网格算法能否用于复杂化学非平衡流数值模拟中,首先将AUSMPW+格式推广到无网格算法中,给出了AUSMPW+在无网格算法中的具体形式。其次,将无网格算法应用于带化学反应的多组分气体Euler方程的求解,给出无网格条件下求解化学非平衡流的具体过程。最后,对NACA0012跨声速绕流场、高超声速弹丸绕流场、氢气/氧气诱导燃烧流场、球锥高超声速绕流场、楔体斜爆轰流场进行了数值模拟。数值模拟结果与相关文献计算结果和实验结果吻合较好,表明了AUSMPW+格式在无网格算法中表现较高精度。同时,采用无网格算法能较好地模拟复杂化学非平衡流场,正确分辨复杂的物理现象,为化学非平衡流数值模拟提供一种新的算法,拓宽了无网格算法的应用范围。     

三维侧压高超声速进气道试验与计算研究

针对三维侧压高超声速进气道开展了试验与计算研究。试验在中国航天空气动力技术研究院的高超声速脉冲风洞中进行,在试验中,来流马赫数为6,文中采用三维有限体积并行计算程序对风洞试验条件下的流动进行了数值模拟,采用k-g湍流模型,时间迭代采用LU隐式格式,对流通量采用了AUSM+、AUSMPW两种不同的空间离散格式进行了对比计算研究,结果表明两个不同的格式都能较好地模拟进气道流场。顶板中心线和离开中心线0.02m位置处的压力分布计算值和试验结果基本一致,表明AUSM+、AUSMPW格式的激波分辨率较高,粘性耗散比较低。计算结果还显示在侧压板壁面出现激波诱导边界层分离现象。     

高阶格式及其在内外流场计算中的应用

高阶格式这里包含两种不同的含义,一种是人们常说的高精度格式;另一种是高分辨率格式,因为它常与高精度有密切联系,但两者毕竟有严格的区别.在格式的构建上,相比之下高分辨率格式较高精度格式难得多.本文从4个方面讨论了高精度格式与高分辨率格式在相关领域中的应用并给出了大量算例,这对澄清人们对高阶格式的认识,弄清它们各自的适用范围是十分必要的.大量的数值计算与研究表明:对于复杂流场,尤其是既考虑激波与边界层相互作用,又要考虑壁面传热的流场,用低阶数值格式计算是无法得到与实验较接近的数值结果的,因此发展高阶格式对处理高温涡轮的气膜冷却问题以及高速飞行器的气动热计算是至关重要的.      

非平衡流计算方法及其模拟激波诱导振荡燃烧

给出一种数值模拟高超声速化学非平衡流动的计算方法，采用流动项与化学反应生成源项解耦处理，即可以有效地解决非平衡流动计算中遇到的刚性问题，还具有公式推导简单和应用方便的优点。从薄层近似N－S方程出发，采用ENO差分格式，数值模拟了超音速H2／Air预混合气体钝体绕流场，较好分辨出流场中激波诱导振荡燃烧非定常流动现象，表明计算方法具有较高的精度。对计算方法的效率和网格数的影响也进行研究。     

LU-AUSMPW混合格式在可压缩无粘和粘性流动数值分析中的应用

本文首先将 AUSMPW格式与三阶 MUSCL格式融合, 给出了其在任意曲线坐标下的三维形式, 并将其与 LU-SGS格式结合, 应用于可压缩 Euler和 Navier-Stokes方程的求解。其次, 构造了一种新的通量限制器。最后, 为了验证 LU -AUSMPW混合格式的性能, 将平面叶栅跨音速无粘流动以及喷管超音速粘性流动作为算例。本文计算结果与文献计算结果和实验数据相符很好, 表明采用 LU-AUSMPW混合格式数值模拟可压缩流场, 具有较高的计算精度、较快的收敛速度和良好的稳定性。      

构建非结构网格高精度有限体积方法的新途径

综述了笔者所在研究团队在发展非结构网格紧致模板高精度有限体积方法方面的研究进展。非结构网格二阶精度有限体积方法在各类商用和自研计算流体力学（CFD）软件中得到了广泛应用。当进一步提高精度时,遇到的主要困难是高阶有限体积方法重构模板过大的问题。这已成为发展非结构网格高精度有限体积方法的主要技术瓶颈之一。近年来,为解决此问题开展了系统研究。基于首先提出的操作紧致性概念,先后提出了3种紧致模板高精度重构方法,包括紧致最小二乘重构、变分重构和多步重构。这些重构方法的共同特点是可在只包含面相邻单元的紧致模板上实施,并达到任意高阶精度。综述了这3种方法,对这些方法的构造思路、实施策略和进一步发展做了概要的阐述。其中变分重构方法将作为非结构网格高精度有限体积方法的方案之一,在国家数值风洞（NNW）工程中得到发展及应用。     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

高速粘性内流的高分辨率高精度迎风型杂交格式

本文在有限体积离散和ＬＵ分解的基础上，构造出一个新的隐式迎风型杂交格式，并用于求解定常流动的稳态解。     

非结构网格紧致高精度气体动理学格式

简要介绍了非结构网格上基于内自由度高效紧致重构方法的几种高精度气体动理学格式(GKS),包括结合三阶GKS通量求解器与子单元有限体积法(SCFV)及通量重构方法(CPR)构造的两种单步时空三阶格式SCFV-GKS和CPR-GKS,以及结合二阶GKS通量和两步四阶方法构造的时空四阶CPR-GKS.进而在CPR-GKS中混...     

贴体曲线坐标系中恢复函数的一种构造方法及应用

本文在结构网格下，从三点迎风紧致逼近出发提出了一种适合于有限体积离散的恢复函数生成办法，在光滑区它具有三阶精度，并且在捕捉激波时有较高的激波分辨率。典型的几个内流算例表明：用该方法得到的恢复函数去计算数值能量时能得到与实验数据较接近的数值解。     

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

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

非结构网格下Euler方程的高分辨率高精度解

提出了一种非结构网格下求解Ｅｕｌｅｒ方程的高分辨率高精度迎风格式。以Ｒｏｅ的矢通量差分分裂为基础,吸收了ＮＮＤ格式的优点,使其具有捕捉激波和滑移线的良好性能;在时间方向上采用Ｊａｍｅｓｏｎ的Ｒｕｎｇｅ－Ｋｕｔａ积分,并结合局部最大时间步长和残差光滑技术加速收敛。最后成功地完成了二维平板激波反射、跨音速Ｌａｖａｌ喷管内的流动和ＧＡＭＭ超音速前台阶绕流等算例,显示了该方法的有效性