一种基于AUSM分裂的真正多维HLL格式

文章给出了一种真正多维的HLL Riemann解算器.采用AUSM分裂将通量分解成为对流通量和压力通量, 其中对流通量的计算采用迎风格式, 压力通量的计算采用HLL格式, 且将HLL格式的耗散项中的密度差用压力差代替, 从而使得格式能够分辨接触间断.为了实现数值格式真正多维的特性, 分别计算了网格界面中点和角点上的数值通量, 并且采用Simpson公式加权组合中点和角点上的数值通量得到网格界面的数值通量.为了减少重构角点处状态时的模板宽度, 计算中采用基于SDWLS梯度的线性重构获得2阶空间精度, 而时间离散采用2阶保强稳Runge-Kutta方法.数值实验表明, 相比于传统的一维HLL格式, 文章的真正多维HLL格式具有能够分辨接触间断, 以及更大的时间步长等优点.与其他能够分辨接触间断的格式(例如HLLC格式)不同, 真正多维的HLL格式在计算二维问题时不会出现激波不稳定现象.      

一种改进的HLLEM格式及其激波稳定性分析

使用低耗散激波捕捉格式对高超声速流动问题进行数值模拟时经常会遭受不同形式的激波不稳定性.本文基于二维无黏可压缩Euler方程,对低耗散HLLEM格式进行激波稳定性分析.结果表明:激波面横向通量中切向速度的扰动增长诱发了格式的不稳定性.通过增加耗散来治愈HLLEM格式的激波不稳定性.为了避免引入过多的耗散进而影响剪切层的分辨率,定义激波探测函数和亚声速区探测函数,使得只有在计算激波层亚声速区的横向数值通量时才增加耗散,其余地方的数值通量依然采用低耗散的HLLEM格式来计算.稳定性分析和数值模拟的结果表明,改进的HLLEM格式不仅保留了原格式高分辨率的优点,还大大提高了格式的鲁棒性,在计算强激波问题时能够有效地抑制不稳定现象的发生.     

高精度非定常激波装配法

为正确模拟高超声速绕流中,来流小扰动与弓形激波之间的干扰对流动特征的影响,将弓形激波作为动边界,利用非定常特征关系处理激波处的边界条件.应用五阶精度迎风紧致格式和六阶精度的对称格式与三阶精度的R-K方法相结合,建立高精度非定常激波装配方法.采用该方法数值模拟钝锥高超声速定常流场和二维抛物外形高超声速边界层流动的感受性问题,数值模拟来流小扰动与弓形激波干扰激波后非定常扰动流场,研究扰动波进入边界层产生边界层不稳定波的特征.     

浅水波方程的黏性正则化PINN算法

针对经典PINN(Physics-informed Neural Networks)在求解浅水波方程间断问题时的不足，提出一种黏性耗散机制的正则化PINN算法。该算法利用黏性正则化的浅水波方程作为网络构建中的物理约束，并在损失函数中作为惩罚项，训练网络用正则化方程的光滑解逼近原方程的间断解，采用网格加密熵稳定格式的数值解作为参考，学习得原方程在整个区域的解。对满足不同初始条件的一维、二维浅水问题进行数值模拟，并与经典PINN算法进行比较，数值结果表明新算法泛化能力强，可预测任意时刻的解，分辨率高，不会出现抹平和伪振荡现象。     

一类新型激波捕捉格式的耗散性与稳定性分析

高超声速流动是高复杂性的可压缩黏性流动, 其中存在激波、剪切层、激波/激波干扰、激波/边界层干扰、旋涡与分离流动等复杂流场结构. 对其进行准确模拟需要使用低耗散、强鲁棒性的激波捕捉方法. 本文基于一类新型的通量项分裂方法, 提出了一种耗散低且鲁棒性好的激波捕捉格式K-CUSP-X. 对该格式的耗散性和激波稳定性进行了详细的理论分析, 得到了格式激波稳定的数值条件. 推论认为, 迎风格式激波稳定的充分条件为速度扰动量具有衰减性, 数值实验验证了该推论. 研究表明, 该格式与Toro提出的通量分裂格式K-CUSP-T相比, 在保证精确捕捉接触间断的同时, 又具有更好的稳定性, 在激波处不会产生“红玉”现象.     

分辨率优化的混合WENO格式

为提高有限差分格式的分辨率,利用傅里叶分析对WENO格式进行色散及耗散优化,并给出优化的线性权重.用优化后的WENO格式与保单调格式（MP）进行加权混合,得到新的加权混合WENO格式（H-WENO）.通过一维激波管问题、Shu-Osher问题及二维双Mach反射问题及R-T不稳定性问题对格式进行数值测试.结果显示,新格式具有强健的激波捕捉能力和对小尺度波结构的高分辨率,与原WENO格式相比改进明显.     

基于波传播算法的火焰不稳定性

基于波传播算法构造了多组分反应流的数值格式,利用CH4空气基元反应动力学模型,并采用分离算法,对CH4空气混合物中,入射激波与火焰的相互作用,以及反射激波与火焰的二次作用过程进行了数值模拟.根据计算结果,讨论了激波诱导火焰失稳的发展过程及其特点.结果表明,Helmholtz不稳定、RichtmyerMeshkov不稳定以及反应放热速率对火焰失稳过程有重要影响.计算结果与实验结果进行了比较,对数值方法的有效性进行了验证.     

气相爆轰波在分叉管中传播现象的数值研究

数值研究气相爆轰波在分叉管中的传播现象.用二阶附加半隐龙格-库塔法和5阶WENO格式求解二维欧拉方程,用基元反应描述爆轰化学反应过程,得到了密度、压力、温度、典型组元质量分数场及数值胞格结构和爆轰波平均速度.结果表明:气相爆轰波在分叉管中传播,分叉口左尖点的稀疏波导致诱导激波后压力、温度急剧下降,诱导激波和化学反应区分离,爆轰波衰减为爆燃波(即爆轰熄灭).分离后的诱导激波在垂直支管右壁面反射,并导致二次起爆.畸变的诱导激波在水平和垂直支管中均发生马赫反射.分叉口上游均匀胞格区和分叉口附近大胞格区的边界不是直线,其起点通常位于分叉口左尖点上游或恰在左尖点.水平支管中马赫反射三波点迹线始于右尖点下游.分叉口左尖点附近的流场中出现了复杂的旋涡结构、未反应区及激波与旋涡作用.旋涡加速了未反应区的化学反应速率.反射激波与旋涡作用并使旋涡破碎.反射激波与未反应区作用,加速其反应消耗,并形成一个内嵌的射流.数值计算得到的波系演变和胞格结构与实验定性一致.     

非定常欠膨胀射流的正格式数值模拟

 采用二阶正格式方法对非定常欠膨胀射流进行了数值模拟。将二维守恒方程的正格式方法推广到轴对称Euler方程组的求解,并对不同马赫数下的燃气射流进行了数值计算。计算结果表明,该方法能够较好地捕捉到包含膨胀波、入射激波、反射激波、马赫盘、射流边界以及三波点等复杂射流流场的波系结构,与实验照片反映的流动特征以及已有的数值结果相吻合。表明该方法对间断解具有较强的捕捉能力,在激波阵面上不会出现数值振荡。     

激波计算不稳定现象机理研究

激波计算不稳定问题是激波捕获格式的一个重要缺陷,多年来受到学术界的广泛关注与研究。本文针对Roe格式,对这一问题进行了机理研究,提出了新的观点,认为造成这一问题的原因在于激波计算格式未能真正考虑流动的多维性,而保留了本应消失的类动量插值机制,即界面修正速度中的压力梯度项。通过这一机理认识,提出了简明的修正方案,改进了Roe格式。数值算例表明这一改进的Roe格式可以有效地消除激波不稳定现象。     

A class of hybrid DG/FV methods for conservation laws I: Basic formulation and one-dimensional systems

By comparing the discontinuous Galerkin (DG) and the finite volume (FV) methods, a concept of ‘static reconstruction’ and ‘dynamic reconstruction’ is introduced for high-order numerical methods. Based on the new concept, a class of hybrid DG/FV schemes is presented for one-dimensional conservation law using a ‘hybrid reconstruction’ approach. In the hybrid DG/FV schemes, the lower-order derivatives of a piecewise polynomial solution are computed locally in a cell by the DG method based on Taylor basis functions (called as ‘dynamic reconstruction’), while the higher-order derivatives are re-constructed by the ‘static reconstruction’ of the FV method, using the known lower-order derivatives in the cell itself and its adjacent neighboring cells. The hybrid DG/FV methods can greatly reduce CPU time and memory required by the traditional DG methods with the same order of accuracy on the same mesh, and they can be extended directly to unstructured and hybrid grids in two and three dimensions similar to the DG and/or FV methods. The hybrid DG/FV methods are applied to one-dimensional conservation law, including linear and non-linear scalar equation and Euler equations. In order to capture the strong shock waves without spurious oscillations, a simple shock detection approach is developed to mark ‘trouble cells’, and a moment limiter is adopted for higher-order schemes. The numerical results demonstrate the accuracy, and the super-convergence property is shown for the third-order hybrid DG/FV schemes. In addition, by analyzing the eigenvalues of the semi-discretized system in one dimension, we discuss the spectral properties of the hybrid DG/FV schemes to explain the super-convergence phenomenon.     

混合通量差分格式的比较

系统研究了几种混合通量差分格式的构造方法和耗散模型,分别对低速平板绕流、二维跨音速喷管流动和高超音速钝头体无粘绕流进行了数值模拟,结合先进的EASM湍流模型对格式的粘性分辨率和激波稳定性进行了细致的比较分析.结果表明混合通量差分格式兼顾了FDS和FVS格式的优点,具有较高的间断分辨率和数值稳定性.     

多介质五方程简化模型及界面捕捉的人工压缩算法

根据两介质五方程简化模型的基本假设,发展了适用于任意多种介质的体积分数方程。为了捕捉多介质界面,将HLLC-HLLCM混合型数值通量的计算格式推广应用于二维平面和柱几何的多介质复杂流动问题,在高阶精度的数据重构过程中采用斜率修正型人工压缩方法ACM。通过一维、二维多介质黎曼问题算例测试,结果表明：发展的计算格式能够较好地分辨接触间断和激波,间断附近物理量无振荡;对于添加了初始扰动的激波问题,能够有效抑制激波数值不稳定性;使用二维柱球SOD问题和接触间断型黎曼问题检验计算格式对多介质复杂流动问题的适应性。     

ADER schemes for the shallow water equations in channel with irregular bottom elevation

This paper deals with the construction of high-order ADER numerical schemes for solving the one-dimensional shallow water equations with variable bed elevation. The non-linear version of the schemes is based on ENO reconstructions. The governing equations are expressed in terms of total water height, instead of total water depth, and discharge. The ENO polynomial interpolation procedure is also applied to represent the variable bottom elevation. ADER schemes of up to fifth order of accuracy in space and time for the advection and source terms are implemented and systematically assessed, with particular attention to their convergence rates. Non-oscillatory results are obtained for discontinuous solutions both for the steady and unsteady cases. The resulting schemes can be applied to solve realistic problems characterized by non-uniform bottom geometries.     

A variable order approach to improve differential quadrature accuracy in dynamic analysis

Differential quadrature (DQ) is a numerical technique which can produce highly accurate results by using a considerably small number of grid points. When it is applied to dynamic equations, however, DQ may exhibit dynamic numerical instability. The present paper analyzed the sources of dynamic numerical instability through a simple example, and the main finding is that dynamic stability is dominated by the grid points near and on boundaries. Based on this, we propose a variable order approach which is characterized by applying different DQ schemes to the grid points near boundaries and grid points far away from boundaries. Numerical examples of both linear and non-linear dynamic equations show that the variable order approach presented in this paper may greatly improve dynamic stability, producing convincing results.     

Long wave run-up on random beaches

The estimation of the maximum wave run-up height is a problem of practical importance. Most of the analytical and numerical studies are limited to a constant slope plain shore and to the classical nonlinear shallow water equations. However, in nature the shore is characterized by some roughness. In order to take into account the effects of the bottom rugosity, various ad hoc friction terms are usually used. In this Letter, we study the effect of the roughness of the bottom on the maximum run-up height. A stochastic model is proposed to describe the bottom irregularity, and its effect is quantified by using Monte Carlo simulations. For the discretization of the nonlinear shallow water equations, we employ modern finite volume schemes. Moreover, the results of the random bottom model are compared with the more conventional approaches.     

Nonlinear spectral-like schemes for hybrid schemes

In spectral-like resolution-WENO hybrid schemes,if the switch function takes more grid points as discontinuity points,the WENO scheme is often turned on,and the numerical solutions may be too dissipative.Conversely,if the switch function takes less grid points as discontinuity points,the hybrid schemes usually are found to produce oscillatory solutions or just to be unstable.Even if the switch function takes less grid points as discontinuity points,the final hybrid scheme is inclined to be more stable,provided the spectral-like resolution scheme in the hybrid scheme has moderate shock-capturing capability.Following this idea,we propose nonlinear spectral-like schemes named weighted group velocity control(WGVC)schemes.These schemes show not only high-resolution for short waves but also moderate shock capturing capability.Then a new class of hybrid schemes is designed in which the WGVC scheme is used in smooth regions and the WENO scheme is used to capture discontinuities.These hybrid schemes show good resolution for small-scales structures and fine shock-capturing capabilities while the switch function takes less grid points as discontinuity points.The seven-order WGVC-WENO scheme has also been applied successfully to the direct numerical simulation of oblique shock wave-turbulent boundary layer interaction.     

A large time step 1D upwind explicit scheme (CFL > 1): Application to shallow water equations

It is possible to relax the Courant–Friedrichs–Lewy condition over the time step when using explicit schemes. This method, proposed by Leveque, provides accurate and correct solutions of non-sonic shocks. Rarefactions need some adjustments which are explored in the present work with scalar equation and systems of equations. The non-conservative terms that appear in systems of conservation laws introduce an extra difficulty in practical application. The way to deal with source terms is incorporated into the proposed procedure. The boundary treatment is analysed and a reflection wave technique is considered. In presence of strong discontinuities or important source terms, a strategy is proposed to control the stability of the method allowing the largest time step possible. The performance of the above scheme is evaluated to solve the homogeneous shallow water equations and the shallow water equations with source terms.     

Well-balanced and energy stable schemes for the shallow water equations with discontinuous topography

We consider the shallow water equations with non-flat bottom topography. The smooth solutions of these equations are energy conservative, whereas weak solutions are energy stable. The equations possess interesting steady states of lake at rest as well as moving equilibrium states. We design energy conservative finite volume schemes which preserve (i) the lake at rest steady state in both one and two space dimensions, and (ii) one-dimensional moving equilibrium states. Suitable energy stable numerical diffusion operators, based on energy and equilibrium variables, are designed to preserve these two types of steady states. Several numerical experiments illustrating the robustness of the energy preserving and energy stable well-balanced schemes are presented.