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

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

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

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

一种MWENO格式的构造和应用

在结构网格上提出了一种新的在混合模板集上守恒量插值的非线性权系数计算方法,构造了一种多重加权实质无波动激波捕捉格式,即所谓的MWENO(Multi-Weighted Essentially Non-Oscillatory)格式.与DWENO(Double Weighted Essentially Non-Oscillatory)格式相比非线性权系数的计算不用人为调节.随后文中给出了几个典型算例.通过MWENO格式和WENO格式的计算结果比较,表明MWENO格式具有较高的间断分辨率.     

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

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

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

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

三阶HWCNS的构造及其在高超声速流动中的应用

对网格质量要求高、计算稳定性差和计算效率低是制约高阶精度格式应用于高超声速复杂流动模拟的重要因素。针对这些问题,发展了三阶精度的混合节点半节点加权紧致非线性格式(HWCNS3),改进其光滑测试因子和非线性权得到了HWCNS3-OP,并给出了它们的频谱特性。利用Lax和Osher-Shu算例测试了格式对间断和高频波的捕捉能力;通过钝锥和航天飞机的高超声速绕流算例,考察了HWCNS3-OP在真实流动模拟中热流和气动力的预测精度及其计算效率。研究结果表明:HWCNS3-OP具有较高的分辨率和良好的间断捕捉能力,高频波捕捉能力相对HWCNS3提高了约3倍,相对守恒律的单调迎风中心格式(MUSCL)提高了约4倍;HWCNS3-OP计算稳定性较好,计算效率相对五阶HWCNS提高了2~3倍,HWCNS3-OP是一种较适合高超声速复杂流动模拟的高阶精度格式。     

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

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

基于高阶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格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。     

一种混合格式的构建及其在多尺度流动中的应用

面向非定常湍流数值模拟方法精度提升问题,构建了实用的矢通量分裂(FVS)形式高阶精度混合格式。该格式平滑区具有五阶精度线性紧致格式的高分辨率特证,而强间断和强梯度区呈二阶耗散型格式特征,进而保证全流场计算的稳定性。基于所建立的混合格式开展有限体积大涡模拟(LES)计算,将其应用于三维高Reynolds数非定常湍流预测中。结果显示,所建立的数值格式能够可靠预测湍流多尺度涡结构特征,且计算量小,具有较高分辨率和鲁棒性。     

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

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

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

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

Godunov型显式大时间步长格式研究进展

综述了Godunov型显式大时间步长格式的研究进展。首先介绍了显式大时间步长格式的概念、分类和优势。然后重点阐述了Godunov型显式大时间步长格式的构造方法、高阶精度推广方法、多维问题推广方法和收敛特性、分辨率及计算效率等性能，展示了其在典型问题中的应用和验证。最后给出了Godunov型显式大时间步长格式研究进一步可能的发展方向。     

消除粘性项高阶离散数值振荡的半结点-结点交错方法

通常的粘性项处理方法是连续两次采用结点型中心差分格式求得.但是,通过两次采用高阶结点型中心格式求得的粘性项易在流场中的"间断"附近产生数值振荡.本文采用多种格式,通过求解一维Burgers方程来充分展示了这种振荡现象.为了消除这种数值振荡,一种半结点型的高阶紧致格式被用来求解粘性项.Fourier频谱分析表明,这种格式具有非常优越的保频谱性能.一维和二维的数值计算结果表明,通过"半结点-结点"交错采用这种半结点型格式的方法可以非常有效地避免粘性项高阶离散可能导致的数值振荡.     

计算气动声学中的高精度紧致差分格式研究

优化了五对角紧致差分格式.通过Fourier分析,将优化目标转化为一个求带有约束的多元非线性函数的最小值问题,利用序列二次规划(SQP)方法获得最佳系数.通过三种措施保证优化的格式具有高精度和分辨率:①直接对量化波数误差积分求其最小值;②采用绝对误差准则,使各种波长的波具有同一误差限;③优化的波数空间和精确求解区间一致.通过调整内部和边界格式的Taylor精度及误差限,保证整个格式的稳定性,并从理论上证明了优化格式具有渐进稳定性.一维和二维基准算例体现了优化格式的性能改进.      

中国CFD史

本文简要回顾了中国的CFD史,早在20世纪70年代,在钱学森教授的建议下,一些中国学者对计算流体动力学(CFD)这一新型学科产生了浓厚兴趣并开始研究.从那时起,中国的CFD逐渐发展繁荣起来,特别是从世纪之交,中国CFD非常活跃,为航空航天和其他民用领域做出了巨大贡献.在中国,特别注重CFD的研究与物理分析的结合.本文首先介绍了M5A的概念,它概括了中国CFD的主要研究领域;随后,介绍了设计格式的原则,以此原则为指导发展了一系列格式,比如NND格式、ENN格式和紧致格式.最近,发展了一系列高阶精度数值格式并应用于实际外形的复杂流动计算.再后,讨论了求解气动方程的网格尺寸准则.最后,同顾了如流动分离、流动拓扑结构和结构稳定性理论、旋涡沿其轴向演化和Hopf分叉理论、动导数和飞行稳定性等代表性工作.     

火箭燃气喷流的数值和实验研究

通过数值计算和实验,研究了２１和９火箭发动机燃气喷流流场。数值计算采用激波高分辨率、二阶精度ＴＶＤ格式,得到包括射流边界、桶形激波、反射激波、Ｍａｃｈ盘和滑移线的清晰的流谱。实验测试应用Ｆ－Ｐ激光干涉仪和自动测压系统,得到流场Ｆ－Ｐ干涉照片和压强分布。计算结果与实验结果符合较好     

二维Euler方程SCB格式的并行实现

本文讨论二维双曲型守恒律方程组的高分辨ＳＣＢ格式，并应用这些格式求解二维Ｅｕｌｅｒ方程组。计算过程中在江南－Ⅱ型并行计算机上进行的，计算过程运用了许多并行计算技术，计算结果令人满意。     

Burgers方程高精度差分格式分析

以非线性扩散波动方程Burgers方程为模型,通过调节方程右端粘性项的大小,在分析了方程的性质的基础上.探讨方程左端项和右端项差分离散格式匹配问题,以及高精度格式稳定性影响和计算格式精度与粘性系数和计算网格之间的关系等问题.研究结果将对方程左右端项离散时应如何采取合适的差分格式有一定的指导意义.     

二维Euler方程组的高分辨率隐式差分格式的并行计算

本文首先对标量双曲型守恒律构造了一类新的二阶隐式ＴＶＤ差分格式，再利用Ｒｏｅ方法和隐式近似因子分解法推广格式到一维、二维守恒率方程组。为减少计算工作量，推导了二维Ｅｕｌｅｒ方程组隐格式的对角形式。最后，针对江南－Ⅱ并行机设计了适合于多处理机结构特征的并行计算程序，计算翼型的跨声速流场，结果令人满意。     

TVD格式数值模拟斜腔内旋涡运动

 采用ＴＶＤ有限体积格式,通过求解非定常可压缩Ｎ－Ｓ方程,对恒定出入口压力比驱动下的斜腔内部流动进行了数值模拟。计算后得到清晰的流场涡系结构及其发展过程；显示了空腔流动中的自维持振荡现象；并分析了剪切层不稳定对周期性振荡的作用机制。