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

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

可压缩流动的高分辨率低耗散混合格式研究(英文)

针对可压缩流场的计算,提出了一种混合型激波捕捉格式来获得对于光滑区域和流场间断都具有高精度的解。该混合方法采用有限体积形式的五阶WENO格式来计算流场的光滑区域,而采用MUSCL格式来捕捉流场中的间断。该混合格式结合WENO格式的低耗散特性和MUSCL格式的高分辨率特性。两种格式之间通过一个间断感受因子进行切换。文中选择了三种典型的间断感受因子进行对比研究。WENO格式和MUSCL格式均为激波捕捉格式,因此混合开关中的常数的选择范围相对较宽。选取了若干典型算例对于混合格式的性能进行了研究,并对不同格式的计算精度和效率进行了对比。数值结果表明,与WENO和MUSCL相比,本文提出的混合方法在较广的范围内均表现出良好的特性。     

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

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

求解多维欧拉方程的二阶旋转输运格式

提出了一个基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法.不同于"网格相关"(Gridaligned)的有限体积方法或者维数分裂的有限差分方法,本格式在求解Riemann问题时不依赖网格的方向,而是沿具有一定的物理意义的两个方向(称为迎风方向).我们发现当迎风方向取为控制体界面两侧速度差矢量方向(及与之正交的方向)时,该格式能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或者所谓"红斑"(carbuncle)现象.为了提高格式的时间和空间精度,我们通过在控制体界面处求解线化的广义Riemann问题的方法体现输运过程对通量计算的影响.这种方案,使得我们有可能以此为基础,构造真正多维的有限体积型迎风格式.     

一种改进的紧致WENO混合格式

给出了一种用于激波捕捉计算的守恒型混合紧致WENO格式.混合格式可视为五阶WENO格式和五阶守恒紧致格式的加权平均.在格式构造中,本文采用了间断分辨率更好的WENO权因子计算方法以及新的子格式权因子计算方法.对于Euler方程,仅对混合格式中WENO格式部分做特征投影处理,不仅获得了特征型WENO格式良好的间断分辨率,同时又保证了计算效率.数值实验验证了本文混合格式的高时间效率和高流场分辨率.     

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

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

一种非结构网格上基于径向基函数重构的ENO格式

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

应用间断伽辽金方法求解二维欧拉方程

在二维非结构网格上,构造了规范正交基作为测试函数,应用不同精度的间断伽辽金方法(DGM)数值求解二维Euler方程,其中数值通量项采用Lax-Friedrichs格式计算.通过对绕NACA0012翼型的跨音速流场的数值模拟,分析了该方法在捕捉间断解方面的特性;同时借鉴有限体积法(FVM)中的一些技术,对求解方案进行了优化.把间断伽辽金方法和有限体积法计算的结果进行比较,认为前者能够更好地处理间断解问题.     

空间推进方法求解抛物化Navier-Stokes方程及其验证

采用伪时间空间推进方法求解完全气体PNS方程。PNS方程采用有限体积方法离散,在推进方向上采用一阶精度迎风格式,在推进面的展向和法向采用三阶精度的MUSCL插值和AUSM类通量构造格式;在推进面上的伪时间推进采用二维LU-SGS迭代方法。通过算例证明,空间推进方法能得到正确的压力、摩阻和热流分布,且计算时间比时间迭代方法快一个数量级以上。     

梯度加权局部分析解差分格式及其应用

提出了一种新的针对Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程符合含有梯度项的对流扩散方程的梯度加权局部分析解的有限差分格式，且给出了通用的系数表达式。ＧＷ－ＬＡＤ格式是基于积分型的有限体积法并且通过附加源项的方法来实现。能够适用于粘性和无粘性及不可压缩和可压缩的压力驱动流动。通过对叶栅流道中二维湍流流的计算表明ＧＷＬＡＤ格式与ＥＤ格式相比，具有较高的精度和很快的收敛速度并且得到更好的压力场和速度场。