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

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

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

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

上风格式的若干性能分析

采用理论分析和数值试验相结合的方法，对上风格式的三个典型代表：FDS、FVS及AUSM格式的性能进行了分析比较：FDS格式具有高的间断分辨率但可能出现非物理懈；FVS格式捕捉强非线性波的能力很强，计算效率高，但耗散比较大；AUSM格式兼有高间断分辨率和高计算效率且克服了两者的弱点。以一维无粘激波管流场为算例考察了这些格式性能与间断分辨率。     

捕捉间断的高精度数值方法

为发展适用于捕捉超声速流场中各种间断的高精度算法,将通量限制的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的通量限制型差分格式.通过在时间方向上利用一阶精度格式计算的一维定常激波,以及在时间方向采用多步Runge-Kutta方法计算的一维非定常激波管问题上的数值试验与二阶精度的TVD格式所计算的结果比较,表明新方法比二阶精度方法在间断的捕捉上具有明显的优势.通过新方法的计算结果与精确解的比较,表明新方法的准度也是非常令人满意的.     

基于边界变差最小的高精度有限差分格式构造

从加权紧致非线性（WCNS）高精度有限差分格式出发,在其重构形式的基础上,根据边界变差最小（BVD）原理,遵循单元边界两侧重构物理量值之差最小的准则,在每个单元内通过两步空间重构,构造了一种新的高精度有限差分格式。一般对WCNS等加权非线性格式的改进都是基于改善色散耗散特性、优化非线性权、提高分辨率等单一途径,将它们作为重构候补函数进行结合,既保持了各自优势所在,又能控制格式整体黏性,所得格式具有丰富的应用场景。通过数值实验,将结果与单一格式进行对比,新格式既能在流场光滑区保证设计的精度,对激波等间断附近的振荡也有很好的抑制作用,提高了对高波数区的分辨率,而且在长时间计算后也有较为精确的结果。面向广泛发展的数值格式,还可以构造出其他新方法,对包含强间断和多尺度的流动问题可以获得更好的结果。     

高阶精度非线性格式WCNS-E-5在二维流动中的应用研究

本文采用具有5阶精度的加权紧致非线性显式格式(WCNS-E-5)对定常与非定常二维流动进行数值模拟,研究表明该格式对各类间断有很好的分辨捕捉能力,而且对强间断如激波的计算,即使在高马赫数与高雷诺数条件下它仍具有很好的收敛性与可靠的计算结果.此外,WCNS-E-5在粗网格条件下也体现出优越性.类如WCNS-E-5的高精度激波捕捉方法将为以后开展湍流数值模拟工作提供坚实的技术保证.     

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

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

适用于超声速流场计算的有限体积紧致算法

结合紧致格式与有限体积法两者的优点,提出了一种基于有限体积的紧致算法.采用特征处理的通量限制方法抑制紧致格式在间断附近的非物理振荡,将该算法应用于航天领域超声速流场的数值模拟.通过计算典型算例,并与其他格式计算结果进行对比,证明了该算法能够无振荡捕捉各种间断,且对波系结构的刻画最为锐利,并能够精细分辨小尺度流场结构,计算结果与精确解最为接近.      

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

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

高阶精度格式WCNS-E-5在亚跨声速流动中的应用研究

采用高阶精度非线性紧致加权格式WCNS-E-5和Baldwin-Lomax模型,求解雷诺平均Navier-Stokes方程,开展了典型翼型与机翼的湍流流动数值模拟研究.对方程中粘性项采用的四阶精度差分近似以及网格导数求解与边界格式的四阶精度,保证了高精度算法的实现.计算结果表明:本文算法能够准确地模拟这些翼型与机翼的亚跨声速流场,得到与实验测量十分吻合的壁面压力分布,计算结果对网格的依赖性小.     

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

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

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

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

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

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

NND格式在多维理想磁流体方程组中的应用

采用修正的四步Runge-Kutta方法求解三维一般曲线坐标系下的理想磁流体方程组,为克服数值振荡,加特征型NND格式进行后处理.特征型NND格式推广到求解三维磁流体(MHD)问题需要知道雅可比通量的左右特征矩阵,在具体计算时需要克服矩阵的奇性.本文用三维程序采用推广的特征NND格式计算了一维MHD激波管和二维(MHD)喷管流动,计算结果表明,特征NND格式保持了TVD格式高精度的优点,又具有计算简单的特点,在包括强弱间断等复杂波系的定常和非定常MHD流场数值模拟中是成功的.     

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

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

A hybrid multidimensional Riemann solver to couple self-similar method with MULTV method for complex flows

Since proposed, the self-similarity variables based genuinely multidimensional Riemann solver is attracting more attentions due to its high resolution in multidimensional complex flows. However, it needs numerous logical operations in supersonic cases, which limit the method’s applicability in engineering problems greatly. In order to overcome this defect, a hybrid multidimensional Riemann solver, called HMTHS (Hybrid of MulTv and multidimensional HLL scheme based on Self-similar structures), is proposed. It simulates the strongly interacting zone by adopting the MHLLES (Multidimensional Harten-Lax-van Leer-Eifeldt scheme based on Self-similar structures) scheme at subsonic speeds, which is with a high resolution by considering the second moment in the similarity variables. Also, it adopts the MULTV (Multidimensional Toro and Vasquez) scheme, which is with a high resolution in capturing discontinuities, to simulate the flux at supersonic speeds. Systematic numerical experiments, including both one-dimensional cases and two-dimensional cases, are conducted. One-dimensional moving contact discontinuity case and sod shock tube case suggest that HMTHS can accurately capture one-dimensional expansion waves, shock waves, and linear contact discontinuities. Two-dimensional cases, such as the double Mach reflection case, the supersonic shock / boundary layer interaction case, the hypersonic flow over the cylinder case, and the hypersonic viscous flow over the double-ellipsoid case, indicate that the HMTHS scheme is with a high resolution in simulating multidimensional complex flows. Therefore, it is promising to be widely applied in both scholar and engineering areas.