An extension of Osher's Riemann solver for chemical and vibrational non-equilibrium gas flows

In this paper we study an extension of Osher's Riemann solver to mixtures of perfect gases whose equation of state is of the form encountered in hypersonic applications. As classically, one needs to compute the Riemann invariants of the system to evaluate Osher's numerical flux. For the case of interest here it is impossible in general to derive simple enough expressions which can lead to an efficient calculation of fluxes. The key point here is the definition of approximate Riemann invariants to alleviate this difficulty. Some of the properties of this new numerical flux are discussed. We give 1D and 2D applications to illustrate the robustness and capability of this new solver. We show by numerical examples that the main properties of Osher's solver are preserved; in particular, no entropy fix is needed even for hypersonic applications.     

Strategies to cure numerical shock instability in the HLLEM Riemann solver

The HLLEM scheme is a popular contact and shear preserving approximate Riemann solver that is known to be plagued by various forms of numerical shock instability. In this paper, we clarify that the shock instability exhibited by this scheme is primarily triggered by the spurious activation of the antidiffusive terms present in the first and third Riemann flux components on the transverse interfaces adjoining the shock front due to numerical perturbations. These erroneously activated terms are shown to counteract the favorable damping mechanism provided by its inherent HLL-type diffusive terms, causing an unphysical variation of the conserved quantity ρu both along and across the numerical shock. To prevent this, two distinct strategies are proposed termed as S elective W ave M odification and A nti D iffusion C ontrol. The former focuses on enhancing the quantity of the favorable HLL-type dissipation available on these critical flux components by carefully increasing the magnitudes of certain nonlinear wave speed estimates, while the latter focuses on directly controlling the magnitude of these critical antidiffusive terms. A linear perturbation analysis is performed to gauge the effectiveness of these cures and to estimate a von Neumann–type stability bounds on the CFL number associated with their use. Results from a variety of classic shock instability test cases show that the proposed strategies are able to provide excellent shock stable solutions even on grids that are highly elongated across the shock front without compromising the accuracy on inviscid contact or shear dominated viscous flows.     

A contribution to the great Riemann solver debate

The aims of this paper are threefold: to increase the level of awareness within the shock-capturing community of the fact that many Godunov-type methods contain subtle flaws that can cause spurious solutions to be computed; to identify one mechanism that might thwart attempts to produce very-high-resolution simulations; and to proffer a simple strategy for overcoming the specific failings of individual Riemann solvers.     

A linearized Riemann solver for the steady supersonic Euler equations

A very simple linearization of the solution to the Riemann problem for the steady supersonic Euler equations is presented. When used locally in conjunction with the Godunov method, computing savings by a factor of about four relative to the use of exact Riemann solvers can be achieved. For severe flow regimes, however, the linearization loses accuracy and robustness. We then propose the use of a Riemann solver adaptation procedure. This retains the accuracy and robustness of the exact Riemann solver and the computational efficiency of the cheap linearized Riemann solver. Numerical results for two- and three-dimensional test problems are presented.     

Calculations of 3D compressible flows using an efficient low diffusion upwind scheme

A newly suggested E‐CUSP upwind scheme is employed for the first time to calculate 3D flows of propulsion systems. The E‐CUSP scheme contains the total energy in the convective vector and is fully consistent with the characteristic directions. The scheme is proved to have low diffusion and high CPU efficiency. The computed cases in this paper include a transonic nozzle with circular‐to‐rectangular cross‐section, a transonic duct with shock wave/turbulent boundary layer interaction, and a subsonic 3D compressor cascade. The computed results agree well with the experiments. The new scheme is proved to be accurate, efficient and robust for the 3D calculations of the flows in this paper. Copyright © 2005 John Wiley & Sons, Ltd.     

Harten-Lax-van Leer-contact (HLLC) approximation Riemann solver with elastic waves for one-dimensional elastic-plastic problems

A Harten-Lax-van Leer-contact (HLLC) approximate Riemann solver is built with elastic waves (HLLCE) for one-dimensional elastic-plastic flows with a hypoelastic constitutive model and the von Mises’ yielding criterion. Based on the HLLCE, a third-order cell-centered Lagrangian scheme is built for one-dimensional elastic-plastic problems. A number of numerical experiments are carried out. The numerical results show that the proposed third-order scheme achieves the desired order of accuracy. The third-order scheme is used to the numerical solution of the problems with elastic shock waves and elastic rarefaction waves. The numerical results are compared with a reference solution and the results obtained by other authors. The comparison shows that the presented high-order scheme is convergent, stable, and essentially non-oscillatory. Moreover, the HLLCE is more efficient than the two-rarefaction Riemann solver with elastic waves (TRRSE).     

基于WENO重构的真正多维黎曼求解器

具有良好守恒性与网格适应性的有限体积格式在流体力学的数值计算中占有重要地位。其中，求解数值流通量是实施有限体积法的关键步骤。一维情形下，通过求解局部黎曼问题来获得数值流通量的相关理论已经比较成熟。但是在计算多维问题时，传统的维度分裂方法仅考虑沿界面法向传播的信息，这不仅影响格式的精度，还可能会造成数值不稳定性从而诱发非物理现象。本文基于对流-压力通量分裂方法来构造真正多维的黎曼求解器，通过求解网格顶点处的多维黎曼问题来实现格式的多维特性。采用五阶WENO重构方法来获得空间的高阶精度，时间离散采用三阶TVD龙格-库塔格式。一系列数值实验的结果表明，真正多维的黎曼求解器不仅具有更高的分辨率还能有效克服多维强激波模拟中的数值不稳定性。     

A finite volume scheme for two-phase compressible flows

In gas-particle two-phase flows, when the concentration of the disperesed phase is low, certain assumptions may be made which simplify considerably the equations one has to solve. The gas and particle flows are then linked only via the interaction terms. One may therefore uncouple the full system of equations into two subsystems: one for the gas phase, whose homogeneous part reduces to the Euler equations; and a second system for the particle motion, whose homogeneous part is a degenerate hyperbolic system. The equations governing the gas phase flow may be solved using a high-resolution scheme, while the equations describing the motion of the dispersed phase are treated by a donor-cell method using the solution of a particular Riemann problem. Coupling is then achieved via the right-hand-side terms. To illustrate the capabilities of the proposed method, results are presented for a case specially chosen from among the most difficult to handle, since it involves certain geometrical difficulties, the treatment of regions in which particles are absent and the capturing of particle fronts.     

一种简单的精确捕捉接触间断的黎曼求解器

基于Godunov型数值格式的有限体积法是求解双曲型守恒律系统的主流方法,其中用来计算界面数值通量的黎曼求解器在很大程度上决定了数值格式在计算中的表现。单波的Rusanov求解器和双波的HLL求解器具有简单、高效和鲁棒性好等优点,但是在捕捉接触间断时耗散太大。全波的HLLC格式能够精确捕捉接触间断,但是在计算中出现的激波不稳定现象限制了其在高马赫数流动问题中的应用。本文利用双曲正切函数和五阶WENO格式来重构界面两侧的密度值,并且结合边界变差下降算法来减小Rusanov格式耗散项中的密度差,从而提高格式对于接触间断的分辨率。研究表明,相比于全波的HLLC求解器,本文构造的黎曼求解器不仅具有更高的接触分辨率,而且还具有更好的激波稳定性。     

A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

An approximate‐state Riemann solver for the two‐dimensional shallow water equations with porosity

PorAS, a new approximate‐state Riemann solver, is proposed for hyperbolic systems of conservation laws with source terms and porosity. The use of porosity enables a simple representation of urban floodplains by taking into account the global reduction in the exchange sections and storage. The introduction of the porosity coefficient induces modified expressions for the fluxes and source terms in the continuity and momentum equations. The solution is considered to be made of rarefaction waves and is determined using the Riemann invariants. To allow a direct computation of the flux through the computational cells interfaces, the Riemann invariants are expressed as functions of the flux vector. The application of the PorAS solver to the shallow water equations is presented and several computational examples are given for a comparison with the HLLC solver. Copyright © 2009 John Wiley & Sons, Ltd.     

The Osher scheme for non-equilibrium reacting flows

An extension of the Osher upwind scheme to non-equilibrium reacting flows is presented, Owing to the presence of source terms, the Riemann problem is no longer self-similar and therefore its approximate solution becomes tedious. With simplicity in mind, a linearized approach which avoids an iterative solution is used to define the intermediate states and sonic points. The source terms are treated explicitly. Numerical computations are presented to demonstrate the feasibility, efficiency and accuracy of the proposed method. The test problems include a ZND (Zeldovich-Neumann-Doring) detonation problem for which spurious numerical solutions which propagate at mesh speed have been observed on coarse grids. With the present method, a change of limiter causes the solution to change from the physically correct CJ detonation solution to the spurious weak detonation solution.     

用鲁棒Riemann求解器和运动重叠网格计算旋翼粘性绕流

发展了一种基于鲁棒Riemann求解器和运动重叠网格技术计算直升机悬停旋翼流场的方法。基于惯性坐标系,悬停旋翼流场是非定常流场,控制方程为可压缩Reynolds平均Navier-Stoke方程,其对流项采用Roe近似Reimann求解器离散,使用改进的五阶加权基本无振荡格式进行高阶重构,非定常时间推进采用含牛顿型LUSGS子迭代的全隐式双时间步方法。为实施旋转运动和便于捕捉尾迹,计算采用运动重叠网格技术。计算得到的桨叶表面压力分布及桨尖涡涡核位置都与实验结果吻合较好。数值结果表明:所发展方法对桨尖涡具有较高的分辨率,对激波具有较好的捕捉能力,该方法可进一步推广到前飞旋翼粘性绕流的计算。     

管系内一维可压缩流动数值建模与特性分析

针对复杂管系内可压缩流体,基于有限体积法,采用HLLC(Harten-Lax-van Leer Contact)格式和黎曼求解器构建了有限控制体数值离散方法,引入虚拟节点用于连接有限控制体,借助虚拟节点给出控制体之间数值通量的计算格式,发展了一种管道内一维流动数值建模方法.针对含有分支管路的管系,在管道连接部位构建了分...     

An approximate solution of the Riemann problem for a realisable second-moment turbulent closure

Abstract. An approximate solution of the Riemann problem associated with a realisable and objective turbulent second-moment closure, which is valid for compressible flows, is examined. The main features of the continuous model are first recalled. An entropy inequality is exhibited, and the structure of waves associated with the non-conservative hyperbolic convective system is briefly described. Using a linear path to connect states through shocks, approximate jump conditions are derived, and the existence and uniqueness of the one-dimensional Riemann problem solution is then proven. This result enables to construct exact or approximate Riemann-type solvers. An approximate Riemann solver, which is based on Gallou?t's recent proposal is eventually presented. Some computations of shock tube problems are then discussed. Received 2 March 1999 / Accepted 24 August 2000     

A finite volume shock‐capturing solver of the fully coupled shallow water‐sediment equations

This paper describes a numerical solver of well‐balanced, 2D depth‐averaged shallow water‐sediment equations. The equations permit variable horizontal fluid density and are designed to model water‐sediment flow over a mobile bed. A Godunov‐type, Harten–Lax–van Leer contact (HLLC) finite volume scheme is used to solve the fully coupled system of hyperbolic conservation laws that describe flow hydrodynamics, suspended sediment transport, bedload transport and bed morphological change. Dependent variables are specially selected to handle the presence of the variable density property in the mathematical formulation. The model is verified against analytical and semi‐analytical solutions for bedload transport and suspended sediment transport, respectively. The well‐balanced property of the equations is verified for a variable‐density dam break flow over discontinuous bathymetry. Simulations of an idealised dam‐break flow over an erodible bed are in excellent agreement with previously published results, validating the ability of the model to capture the complex interaction between rapidly varying flow and an erodible bed and validating the eigenstructure of the system of variable‐density governing equations. Flow hydrodynamics and final bed topography of a laboratory‐based 2D partial dam breach over a mobile bed are satisfactorily reproduced by the numerical model. Comparison of the final bed topographies, computed for two distinct sediment transport methods, highlights the sensitivity of shallow water‐sediment models to the choice of closure relationships. Copyright © 2016 John Wiley & Sons, Ltd.     

A problem for diffusion with relaxation in polymers

A nonlinear problem for penetrant diffusion with relaxation in polymers is considered. A numerical approach to solving this type of problems is developed. The proposed numerical scheme based on a finite element domain approximation and a time difference method can be used for numerical simulation of the considered penetrant diffusion in 2-D and 3-D domains. A numerical procedure and a corresponding computer code are created and tested for some examples in 1-D and 2-D domains.     

A Riemann solver for the two-dimensional MHD equations

This paper presents how the equations of magnetohydrodynamics (MHD) in primitive form should be written in conservative form with the inclusion of a divergence source along with a divergence wave and how a physically correct sonic fix can be embedded directly in the fluxes. The numerical scheme was applied to a blast wave problem in which a circular energetic plasma is released in a free and magnetized medium with a reflected wall. The results show that the method with the new sonic fix can handle the divergence condition on the magnetic field and produces an almost uniform shock compression in all directions, resolving the shocks and discontinuities rather sharply. © 1997 John Wiley & Sons, Ltd.