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

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

一种求解激波问题的中心差分-WENO混合方法研究

为提高求解包含激波问题的计算精度和效率,发展了一种低耗散、高效率的中心差分-WENO(Weighted Essentially Non-Oscillatory)混合格式,Navier-Stokes方程的无黏项在光滑流场区域采用六阶中心差分格式离散,而间断附近采用五阶WENO格式求解;基于密度设计了一种新型格式开关实现两种格式在光滑区域和间断之间自动切换,确保数值解在间断附近基本无振荡。针对一维激波熵波作用问题验证了混合格式的低耗散特性;对二维Riemann问题的研究表明发展的基于密度的格式开关更为合理,混合格式的计算效率较WENO格式明显提高;将其应用到激波诱导燃烧问题中,混合格式能很好地捕捉定常流场中激波和化学反应锋面的位置以及非定常流场的振荡频率。     

特定条件下高阶WENO格式计算结果误差

通过比较一阶迎风格式和五阶WENO格式模拟激波、接触间断、曲线坐标系下的均匀流和激波正规反射等4个简单流场得到的数值结果,发现WENO格式模拟的激波和接触间断在从初始间断变化成数值过渡区的过程中出现的非物理波动比一阶迎风格式的结果更加明显,流场结构也更加复杂;同时,由坐标变换而产生的几何诱导误差和边界近似模型误差也明显...     

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

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

自适应的高精度中心﹣WE NO混合格式

为对可压缩流动进行高精度高分辨率的数值模拟,基于模板的光滑度量函数,构造了新型的加权因子,在此基础上将7阶精度的WENO格式与8阶中心格式进行加权组合,构造了一种自适应低耗散的中心-WENO混合格式(H-WENO7-CD8),并采用Fourier方法对离散格式数学特性进行了理论分析。该自适应的高精度中心-WENO混合格式相比于七阶的WENO格式具有更小的耗散误差。通过对激波-密度干扰问题、Rayleigh-Taylor不稳定性问题和双Mach反射问题的数值模拟,并与WENO-JS格式的计算结果对比,结果表明:该格式结合了WENO格式和中心格式的特点,能更好地捕捉激波,对物理脉动也具有较高的分辨率,适合于复杂流体流动的数值模拟。     

WENO格式的一种新型光滑因子及其应用

光滑因子是构造WENO格式的关键，决定了WENO格式能否达到最优收敛精度以及在间断附近能否保持本质无振荡特性。针对广泛应用的五阶WENO格式，采用算子函数近似导数关系的方法，设计了一种新型光滑因子。在间断区域，与传统的光滑因子相比，新型光滑因子对间断的识别更准确。在光滑区域，新型光滑因子使得新型WENO格式的非线性权更接近线性权。理论分析和数值验证表明新型WENO格式即使在一阶临界点也能保持一致五阶精度。一维典型激波问题的数值结果表明，与现有WENO格式相比，新型WENO格式提高了短波的分辨率，降低了格式的耗散，同时能够准确识别间断。对于典型包含激波和剪切层的流动问题，新型WENO格式不仅能够准确地捕捉强激波，而且能够精细模拟剪切层和声波等流场结构。     

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

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

ENO-AUSMPW格式对激波涡干扰的数值模拟

高阶ENO格式在复杂流场计算中具有重要的地位,它理论上可以达到任何阶的精度。本文采用二维高阶ENO插值方法,构造了高阶AUSMPW格式,并对包含复杂激波结构的激波涡干扰流场进行了数值模拟。计算结果表明了高阶ENO插值下的AUSMPW格式具有较高的激波分辨率和较低的数值耗散。     

一种求解双曲守恒律方程的新方法

近年来,在计算流体力学中,WENO(WeightedEssentiallyNon-oscillatory)格式得到了很大的发展,这也促进了基于WENO格式基础上的CWENO(CentralWENO)格式的发展。WENO格式的主要思想是通过凸组合的方式进行插值重构,以得到高阶逼近,并且在间断附近防止伪震荡的产生。本文构造了一种新的五阶紧凑CWENO格式,该格式所选模板节点少,并且保留了原有中心格式的优点。最后通过一些数值算例来验证该格式的高精度,本性无震荡性质。     

基于高阶WENO格式的旋翼非定常涡流场数值模拟

为提高直升机旋翼黏性涡流场计算流体力学（CFD）的计算精度及降低数值耗散,发展了一套基于五阶加权本质无振荡（WENO）格式的旋翼非定常涡流场数值计算方法。采用运动嵌套网格方法生成围绕前飞状态旋翼的网格系统,以Navier-Stokes方程作为主控方程,湍流模型采用一方程Spalart-Allmaras模型。为了提高旋翼流场中对涡的形成、演化等发展过程的模拟精度,采用高间断分辨率的Roe-WENO格式计算对流通量。对状态变量的插值,选取适当的加权因子,通过多个重构模板的凸组合来构造五阶WENO格式,在交界面附近获得具有五阶精度的状态变量。为了提高计算效率,采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）推进方法和并行计算策略。最后,运用该方法分别对悬停状态下C-T（Caradonna-Tung）旋翼和UH-60A旋翼以及前飞状态下C-T旋翼和SA349/2旋翼的涡流场及气动特性进行了数值模拟,并将数值结果与传统二阶精度格式的计算值进行了对比。结果表明：在计算时间增长10%～20%的前提下,五阶WENO格式能够捕捉到更精确的涡流场特性,反映了五阶WENO格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。     

低耗散高可信度格式及其在气动力热计算中的应用

本文简要介绍了TVD类型格式中近似黎曼解的耗散机理,该耗散机制对于粘性问题的计算有较大的负面影响.通过引入低耗散的格式,所发展的数值方法对粘性流场的气动力热问题的计算具有较高的精度.以该数值方法为基础发展的软件系统,在亚跨超和高超声速领域中的气动力热问题上均得到了验证,计算结果与实验进行了对比验证,表明该计算软件系统具有较高的可信性.     

NND格式的推广及在粘流计算中的应用

本文把ＮＮＤ格式推广到了特征变量，守恒变量和原始变量形式。通过引入通量差分分裂方法，由守恒变量形式的ＮＮＤ格式得到了Ｒｏｅ的近拟Ｒｉｅｍａｎｎ解法，并讨论了限制器的作用。在隐式ＮＮＤ格式中引入了一种分解算法，提高了计算效率。本文的数值计算表明，ＮＮＤ格式对复杂波系干扰，边界层分离与再附等流动现象均有较高的分辨率。     

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阶,甚至更高阶重构;数值耗散更低,捕获更准确的激波位置,对航空工程数值模拟具有重要意义.     

TVD格式计算一维Euler方程间断解存在的几个理论问题

本文讨论了TVD格式按本文第二节A段意义推广到拟线性方程组的几个理论问题。指出:1.对于Riemann间断分解,Euler方程的解和Riemann不变量都没有TVD性质,因此不能将TVD格式自然推广;2.TVD格式是提高分辨率的一种途径,它不能保证得到物理解,能得到物理解的TVD格式一定是耗散守恒格式,因此将TVD格式推广时,必须保证格式是耗散的,否则仍会得到非物理解。并给出了一阶迎风格式,TVD2,UNO2得到的若干非物理解算例;3.非物理解现象只是发生在具有剧烈膨胀的单波区的计算上,而某些TVD格式若会得到非物理解,则一定发生在特征值变号的单波区的计算上。     

利用紧致格式捕捉间断的数值方法研究

将通量限制或加权的思想引入到紧致格式中,构造了两类传统方法与紧致格式混合组成的差分格式:通量限制与加权型差分格式。通过方波、组合波、定常激波、非定常Sod问题、Shu问题和Lax问题上的计算,以及与精确解的比较,结果表明这两种方法在间断的捕捉上具有高精度、高分辨率,而在计算方波、组合波或Shu问题和Lax问题上,加权格式具有更大的优势。     

求解双曲型守恒律的五阶松弛格式

给出了一种求解一维双曲型守恒律的五阶松弛格式.该格式以五阶WENO重构和显隐式Runge-Kutta方法为基础.本文格式保持了松弛格式简单的优点,即不用Riemann解算器和计算非线性通量函数的雅可比矩阵.用该格式对一维Euler方程进行了数值试验,并与三阶和四阶松弛格式的计算结果进行了比较,结果表明本文的格式具有更低的数值耗散和更高的分辨率.     

跨声速压气机数值模拟格式效应研究

针对跨声速叶轮机械复杂内部粘性流动,开发了全三维高精度RANS数值求解器。分别采用FVS格式与AUSM+格式结合多种限制器对Rotor 67跨声速压气机转子进行数值模拟,分析和比较了各格式的计算效果并与试验结果对比。结果表明:与FVS格式相比,AUSM+格式的数值粘性更小,边界层的模拟精度更高;Hemker限制器的综合表现最优,Van Albada限制器的粘性分辨率略低于Hemker限制器,Minmod限制器对流动分离现象的捕捉能力较差,Van Leer限制器容易引入色散误差。     

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

简要介绍了非结构网格上基于内自由度高效紧致重构方法的几种高精度气体动理学格式(GKS),包括结合三阶GKS通量求解器与子单元有限体积法(SCFV)及通量重构方法(CPR)构造的两种单步时空三阶格式SCFV-GKS和CPR-GKS,以及结合二阶GKS通量和两步四阶方法构造的时空四阶CPR-GKS。进而在CPR-GKS中混合SCFV增强对流场间断的分辨能力,并拓展到三维六面体网格。通过几种典型数值算例,对比分析了这几种格式在可压缩流动问题中具有的高精度、高效率和良好的间断捕捉能力。GKS通量的多维时空演化特性与紧致重构方法的高效性为发展高精度格式提供了有力的支撑。     

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

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

Relativistic MHD Simulations of Jets with Toroidal Magnetic Fields

This paper presents an application of the recent relativistic HLLC approximate Riemann solver by Mignone & Bodo to magnetized flows with vanishing normal component of the magnetic field. The numerical scheme is validated in two dimensions by investigating the propagation of axisymmetric jets with toroidal magnetic fields. The selected jet models show that the HLLC solver yields sharper resolution of contact and shear waves and better convergence properties over the traditional HLL approach.