二维非结构网格上的高精度有限体积WENO格式

构造了非结构网格上二维双曲型守恒律的一类新的高精度有限体积WENO格式。其主要思想是:根据格式精度的要求,按照谱体积方法对三角形单元网格进行剖分,通过选取适当的子单元组成模板,利用WENO重构方法重构二阶和三阶多项式,利用有限体积公式和高阶Runge-Kutta TVD时间离散方法,构造了非结构网格上二维双曲型守恒律的一致二阶和三阶精度的有限体积WENO格式。然后,推广到二维Euler方程组。最后,给出几个数值算例,验证了格式的稳定性、高阶精度和高分辨捕捉激波等间断的能力。     

求解双曲守恒律方程的自适应人工黏性熵稳定格式

采用熵守恒格式求解双曲守恒律方程不满足熵稳定条件.为了达到熵稳定,需要在熵守恒格式的基础上增加黏性机制,使总熵耗散.鉴于自适应人工黏性在不同的计算区域上其黏性大小的变化特点,在熵守恒通量的基础上添加经过修正的自适应人工黏性通量构造出一种新的结构简单的熵稳定格式.根据自适应人工黏性的自适应性,并经过简单调节黏性比例系数C之后,该格式可以在不使用限制器的情况下达到整体的高分辨率,表现出:在间断区域有足够的数值耗散保证稳定;在光滑区域,黏性很小,不会影响格式在此的高精度.最后的几个数值算例以熵稳定ERoe格式和高分辨率熵稳定EYee格式的数值结果作为参照,分析说明这种格式的特性.     

基于二维结构化网格的间断Galerkin方法研究

在二维结构化网格上,对求解Euler方程的间断Galerkin方法进行了研究。应用二阶、三阶精度的问断Galerkin方法对绕NACA0012翼型的跨音速无粘流动进行了数值模拟。在整体逼近精度基本保持不变的前提下,采用moment限制器以消除数值解在激波附近的伪振荡。数值结果表明,间断Galerkin方法能够很好地模拟流场,准确捕捉激波;moment限制器成功抑制了激波附近的伪振荡;高阶格式具有比低阶格式更强的激波捕捉能力。     

适用于气动学问题的DRP类有限差分方法

研究了一类低耗散低色散的高阶精度有限差分方法,目的是直接计算非定常欧拉方程用于气动声学问题.采用的数值方法是空间四阶、时间三阶精度的色散关系保持(DRP:dispersion-relation-preserving)类有限差分格式,通过波动方程算例验证了格式模拟波动问题的能力;采用加变幅值的高波数人工耗散项来抑制高频数值振荡,使得该格式可推广用于含激波、初始间断或非线性波动的问题.通过含激波、初始间断条件的算例验证了所用变幅值人工耗散的有效性,最后作为初步应用计算了二维亚、超音速均匀流中简单点声源辐射问题,得到了很好的结果,表明此类方法在气动声学问题计算上深有潜力.     

混合动理学通量WENO方法拓展

对基于气体动理论的高精度WENO方法进行了拓展研究,将Int.J.Numer.Meth.Fluids 79(6), 290-305(2015)中提出的5阶混合动理学通量WENO方法进一步拓展到7阶和9阶;在5阶混合动理学WENO方法框架下,比较了采用不同激波探测技术的计算准确度和效率。在时间方向采用三阶Runge-Kutta方法进行推进,空间离散采用混合动理学通量WENO方法。通过一维和二维算例,验证了构造的7阶和9阶混合动理学方法比传统的矢通量分裂技术具有更高的分辨率和更小的数值耗散,并保持了较好的激波捕捉能力;对比发现Ohwada等提出的激波探测技术具备良好的激波探测能力和计算效率。     

非结构网格上的TV-van Leer混合迎风格式

为了改进迎风型对流通量格式的激波捕捉计算稳定性,对迎风格式产生数值激波不稳定现象的原因进行分析,并且在分析结果基础上构造新型迎风型格式。为此,首先设计了数值算例,分析非结构网格上迎风格式出现数值激波不稳定现象的主要影响因素。结果表明,在一定条件下,网格和流动状态共同作用导致产生不合理的数值结果,并影响计算的稳定。其中,网格间的动量通量对数值激波不稳定现象影响显著。因此,基于TV格式构造了一种自适应调节动量通量耗散项的混合格式,从而在激波附近引入耗散,抑制扰动的产生和传播。数值计算结果显示新型混合格式在非结构网格离散条件下能够稳定地捕捉激波,同时混合格式保持了TV格式原有的黏性流动计算分辨率。     

双曲守恒型方程的二阶摄动有限差分格式

对双曲守恒型方程，将其一阶迎风格式空间差商的常系数摄动展开为时间步长和空间步长的幂级数，通过确定幂级数系数而获得二阶精度的摄动有限差分(PFD)格式。进而从双曲守恒型方程的通量分裂型一阶迎风格式出发，通过类似的摄动展开方法，获得空间精度为二阶的通量分裂形式的摄动有限差分(FPFD)格式。这两类格式保留了一阶守恒迎风格式的简洁结构形式，使用三节点即可达到二阶精度，又避免了三点二阶格式的非物理数值振荡。并将这两类格式推广应用到双曲守恒型方程组，最后通过模型方程和一维激波管流动的数值算例验证了格式的高精度、高分辨率性质。     

一个基于通量分裂的高精度MmB差分格式

本文研究双曲型守恒律的高精度差分方法.在新的算法中,首先将计算区间划分为互不相交的小区间,再根据精度要求等分小区间;其次,根据流动方向进行通量分裂,重构小区间交界面上的正、负数值通量,并进行校正;然后,采用高阶Runge-Kutta TVD方法进行时间离散,构造了一维非线性双曲型守恒律方程的一个高精度、高分辨率的守恒型差分格式.推广到二维双曲型守恒律方程,证明了格式的MmB特性.进而推广到二维守恒型方程组情形.最后对二维Burgers方程及Euler方程进行了数值试验,数值结果令人满意.     

基于人工粘性的间断Galerkin有限元方法

发展了一种基于人工粘性的间断Galerkin有限元方法,作出改进以增强适应性,并向非均匀网格推广.选取了典型算例对方法进行验证.一维激波管算例表明,改进的方法在保持计算结果高分辨率的同时,能够更好地抑制非物理振荡.分析得知,当人工粘性方法用于二阶DG格式时,所得计算结果的数值耗散较大,而当格式精度大于二阶时,采用人工粘性方法所得的结果的分辨率较高.通过计算圆柱高超声速粘性绕流,将三阶DG格式与三阶MUSCL格式和五阶WENO格式的结果进行对比,结果表明,该人工粘性方法对于高超声速流动计算也具有一定的优势.     

高超声速湍流直接数值模拟技术

概述了近年来国内外高超声速湍流直接数值模拟(DNS)技术方面最新的研究进展,主要集中在高精度、高鲁棒性数值方法方面,同时也介绍了近年来典型的高超声速湍流DNS算例。在数值方法方面,主要介绍了高精度激波捕捉格式以及保持计算稳定的数值技术,重点是WENO格式及高阶保单调格式的最新进展。在高超声速湍流DNS算例方面,介绍了压缩性影响、壁温影响、真实气体效应以及高超声速转捩等方面的DNS研究。此外,还简要介绍了作者开发的可压缩高精度计算流体力学软件OpenCFD。     

一种修正的Osher通量在高阶WCNS中的特性

将一种基于多维修正的Osher通量运用于高阶加权紧致非线性格式(WCNS)中,该修正通量主要在垂直于激波面的界面上增加耗散,能够改善Osher通量的激波捕捉稳定性,同时对边界层和接触间断的分辨率影响非常小。对修正的Osher通量在高阶WCNS中的特性进行研究,通过数值模拟测试了基于Osher通量的WCNS的激波稳定性、热流预测精度、边界层模拟能力、激波边界层干扰模拟能力,并与Steger-Warming通量和Roe通量进行了对比。结果表明修正后的Osher通量比Harten修正的Roe通量具有更好的激波捕捉鲁棒性,而边界层、驻点热流值和激波边界层干扰的模拟则明显优于Steger-Warming通量。上述结果说明了基于修正的Osher通量的高阶WCNS具有较好的激波捕捉特性、热流预测精度和边界层计算能力。     

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

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

适用于超声速的一种通量限制型紧致格式

紧致格式因其结构简单、在相同的网格点上能达到比非紧致格式更高的精度以及与谱方法相近的分辨率等优点,日益受到人们的重视。用紧致格式模拟超声速流场的主要问题之一是如何保证高阶紧致格式能光滑地捕捉到流场的各种间断。本文借鉴NND格式的思想,构造出一种总体上具有三阶精度的通量限制型紧致(LFC)格式,并成功地应用于含有激波、滑移面等复杂流动现象的数值模拟。计算结果表明这种格式不仅具有较高的精度和分辨率,而且还保证了在间断附近基本无虚假波动。     

高精度差分求解气动方程的几个问题

本文探讨了发展高精度格式的必要性，研究了高精度格式与网格以及熵增条件的关系，并发展一种半离散化的空间为三阶精度的格式。从此半离散化的格式出发，可建立多步显式格式及隐式格式。模型问题的计算表明，该三阶精度的格式具有好的精度，且激波附近基本上没有波动。     

求解双曲型守恒律的高精度,无波动样条逼近有限体积方法

本文从三次样条逼近出发，提出了求解流体力学双曲型守恒律的一种高精度、无波动的数值方法。针对样条函数的特性，本文提出了一种新的通量限制技术，使该方法在光滑区可以达到四阶精度，在流动参数的空间分布出现拐点或极值点时分别退化为二阶或一阶精度的格式。数值实验表明，该方法对流场中的激波和接触间断有很高的分辨率，优于二阶精度的ＴＶＤ格式。     

HIGH ORDER ACCURACY SCHEME FOR 2-D TRANSONIC FLOWS

ＨＩＧＨＯＲＤＥＲＡＣＣＵＲＡＣＹＳＣＨＥＭＥＦＯＲ２－ＤＴＲＡＮＳＯＮＩＣＦＬＯＷＳ￥ＬｉＨａｉｄｏｎｇ；ＬｉｕＱｉｕｓｈｅｎｇ；ＳｈｅｎＭｅｎｇｙｕ（ＤｅｐａｒｔｍｅｎｔｏｆＥｎｇｉｎｅｅｒｉｎｇＭｅｃｈａｎｉｃｓ，ＴｓｉｎｇｈｕａＵｎｉｖｅｒｓ...     

关于建立高阶差分格式的问题

为了能在不太密的网格上捕捉到流场的细致结构。通常采用高阶精度的差分 式进行数值模拟。为能抑制计算中和在激波附近产生的虚假波动,本文从物理构思出发,提出了建立高阶格式的两个基本原则,作者称之为抑制波动的原则和稳定性原则。     

非结构网格上新型的NND有限元格式

张涵信发展的ＮＮＤ差分格式是由中心差分格式,二阶迎风格式和一阶迎风格式混合组成的杂交型格式。众所周知,和中心差分格式相对应的是Ｇａｌｅｒｋｉｎ有限元格式。通过对中心型有限元格式加修正项的方法本文成功地构造出二阶迎风型有限元格式和一阶迎风型有限元格式。     

高超声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析，研究了Steger-Warming通矢量分裂、Van Leer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。