Explicit solution to the exact Riemann problem and application in nonlinear shallow‐water equations

The Riemann solver is the fundamental building block in the Godunov‐type formulation of many nonlinear fluid‐flow problems involving discontinuities. While existing solvers are obtained either iteratively or through approximations of the Riemann problem, this paper reports an explicit analytical solution to the exact Riemann problem. The present approach uses the homotopy analysis method to solve the nonlinear algebraic equations resulting from the Riemann problem. A deformation equation defines a continuous variation from an initial approximation to the exact solution through an embedding parameter. A Taylor series expansion of the exact solution about the embedding parameter provides a series solution in recursive form with the initial approximation as the zeroth‐order term. For the nonlinear shallow‐water equations, a sensitivity analysis shows fast convergence of the series solution and the first three terms provide highly accurate results. The proposed Riemann solver is implemented in an existing finite‐volume model with a Godunov‐type scheme. The model correctly describes the formation of shocks and rarefaction fans for both one and two‐dimensional dam‐break problems, thereby verifying the proposed Riemann solver for general implementation. Copyright © 2007 John Wiley & Sons, Ltd.     

Riemann solution for ideal isentropic magnetogasdynamics

In the present paper, we study the Riemann problem for quasilinear hyperbolic system of partial differential equations governing the one dimensional ideal isentropic magnetogasdynamics with transverse magnetic field. We discuss the properties of rarefaction waves, shocks and contact discontinuities. Differently from single equation methods rooted in the ideal gasdynamics, the new approach is based on the system of two nonlinear equations imposing the equality of total pressure and velocity, assuming as unknowns the two values of densities, on both sides of the contact discontinuity. Newton iterative method is used to obtain densities. The resulting exact solver is implemented with the examples of general applicability of the proposed approach. For comparisons with exact solution we also shown numerical results obtained by the total variation diminishing slope limiter centre scheme. It is shown that both analytical and numerical results demonstrate the broad applicability and robustness of the new Riemann solver.     

The Riemann problem for a hyperbolic model of two-phase flow in conservative form

This article is to continue the present author's work (International Journal of Computational Fluid Dynamics (2009) 23 (9), 623–641) on studying the structure of solutions of the Riemann problem for a system of three conservation laws governing two-phase flows. While existing solutions are limited and found quite recently for the Baer and Nunziato equations, this article presents the first instance of an exact solution of the Riemann problem for two-phase flow in gas–liquid mixture. To demonstrate the structure of the solution, we use a hyperbolic conservative model with mechanical equilibrium and without velocity equilibrium. The Riemann problem solution for the model equations comprises a set of elementary waves, rarefaction and discontinuous waves of various types. In particular, such a solution treats both the wave structure and the intermediate states of the two-phase gas–liquid mixture. The resulting exact Riemann solver is fully non-linear, direct and complete. On this basis then, we use locally the exact Riemann solver for the two-phase flow in gas–liquid mixture within the framework of finite volume upwind Godunov methods. In order to demonstrate the effectiveness and accuracy of the proposed solver, we consider a series of test problems selected from the open literature and compare the exact and numerical results by using upwind Godunov methods, showing excellent oscillation-free results in two-phase fluid flow problems.     

A compact high‐order HLL solver for nonlinear hyperbolic systems

We present a new finite‐volume method for calculating complex flows on non‐uniform meshes. This method is designed to be highly compact and to accurately capture all discontinuities that may arise within the solution of a nonlinear hyperbolic system. In the first step, we devise a fourth‐degree Hermite polynomial to interpolate the solution. The coefficients defining this polynomial are calculated by using a least‐square method. To introduce monotonicity conditions within the procedure, two constraints are added into the least‐square system. Those constraints are derived by locally matching the high‐order Hermite polynomial with a low‐order TVD polynomial. To emulate these constraints only in regions of discontinuities, data‐depending weights are defined; these weights are based upon normalized indicators of smoothness of the solution and are parameterized by an O(1) quantity. The reconstruction so generated is highly compact and is fifth‐order accurate when the solution is smooth; this reconstruction becomes first order in regions of discontinuities. In the second step, this reconstruction is inserted in an HLL approximate Riemann solver. This solver is designed to correctly capture all discontinuities that may arise into the solution. To this aim, we introduce the contribution of a possible contact discontinuity into the HLL Riemann solver. Thus, a spatially fifth‐order non‐oscillatory method is generated. This method evolves in time the solution and its first derivative. In a one‐dimensional context, a linear spectral analysis and extensive numerical experiments make it possible to assess the robustness and the advantages of the method in computing multi‐scale problems with embedded discontinuities. Copyright © 2008 John Wiley & Sons, Ltd.     

Restoration of the contact surface in the HLL-Riemann solver

The missing contact surface in the approximate Riemann solver of Harten, Lax, and van Leer is restored. This is achieved following the same principles as in the original solver. We also present new ways of obtaining wave-speed estimates. The resulting solver is as accurate and robust as the exact Riemann solver, but it is simpler and computationally more efficient than the latter, particulaly for non-ideal gases. The improved Riemann solver is implemented in the second-order WAF method and tested for one-dimensional problems with exact solutions and for a two-dimensional problem with experimental results.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

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

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

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

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

弱不连续问题高阶有限元离散系统的GAMG法

弱不连续问题(如含夹杂问题)是固体力学计算中的一类重要问题。高阶有限元方法由于其具有更好的逼近效果,是确保数值解在界面保持较高精度的计算方法之一。但与线性元相比,高阶单元需要更多的计算机存储单元,具有更高的计算复杂性。本文利用两水平算法的思想,将高阶有限元离散系统化归于线性元离散系统的求解,为弱不连续问题高阶有限元离散系统设计了一种新的基于几何与分析信息的代数多重网格(GAMG)法,并应用于圆形求解域含单夹杂问题的高阶有限元离散系统的求解。数值试验结果表明,相比于常用GAMG法,新方法的迭代次数基本不依赖于问题规模、单元阶次以及杨氏模量的间断性,CPU计算时间得到明显改善,具有更好的计算效率和鲁棒性,可大大提高弱不连续问题有限元分析的整体效率。     

Derivative artificial compression method for conservation laws with multi-dimensional extension

A new resolution-enhancing technique called derivative artificial compression method is developed with multi-dimensional extension. The method is constructed via applying high-resolution difference schemes on derivative equations of conservation laws. In this way, one could overcome the defect of accuracy decay at extreme points that has plagued almost all high-resolution schemes. The new method has high resolution, low dissipation and low diffusion properties, and could enhance the resolution (of numerical solution) both at discontinuities and at extreme points. Numerical experiments are implemented using initial value problems of single conservation law, one-dimensional shock-tube problem, two-dimensional Riemann problems, double Mach reflection problem, and a shock reflection from a wedge. Resolutions of discontinuities, extremes and fine structures are compared between the original TVD scheme, TVD scheme with artificial compression method and TVD scheme with derivative artificial compression method.     

Towards numerical simulations of supersonic liquid jets using ghost fluid method

A computational tool based on the ghost fluid method (GFM) is developed to study supersonic liquid jets involving strong shocks and contact discontinuities with high density ratios. The solver utilizes constrained reinitialization method and is capable of switching between the exact and approximate Riemann solvers to increase the robustness. The numerical methodology is validated through several benchmark test problems; these include one-dimensional multiphase shock tube problem, shock–bubble interaction, air cavity collapse in water, and underwater-explosion. A comparison between our results and numerical and experimental observations indicate that the developed solver performs well investigating these problems. The code is then used to simulate the emergence of a supersonic liquid jet into a quiescent gaseous medium, which is the very first time to be studied by a ghost fluid method. The results of simulations are in good agreement with the experimental investigations. Also some of the famous flow characteristics, like the propagation of pressure-waves from the liquid jet interface and dependence of the Mach cone structure on the inlet Mach number, are reproduced numerically. The numerical simulations conducted here suggest that the ghost fluid method is an affordable and reliable scheme to study complicated interfacial evolutions in complex multiphase systems such as supersonic liquid jets.     

A hybrid pressure‐based solver for nonideal single‐phase fluid flows at all speeds

This paper describes the implementation of a numerical solver that is capable of simulating compressible flows of nonideal single‐phase fluids. The proposed method can be applied to arbitrary equations of state and is suitable for all Mach numbers. The pressure‐based solver uses the operator‐splitting technique and is based on the PISO/SIMPLE algorithm: the density, velocity, and temperature fields are predicted by solving the linearized versions of the balance equations using the convective fluxes from the previous iteration or time step. The overall mass continuity is ensured by solving the pressure equation derived from the continuity equation, the momentum equation, and the equation of state. Nonphysical oscillations of the numerical solution near discontinuities are damped using the Kurganov‐Tadmor/Kurganov‐Noelle‐Petrova (KT/KNP) scheme for convective fluxes. The solver was validated using different test cases, where analytical and/or numerical solutions are present or can be derived: (1) A convergent‐divergent nozzle with three different operating conditions; (2) the Riemann problem for the Peng‐Robinson equation of state; (3) the Riemann problem for the covolume equation of state; (4) the development of a laminar velocity profile in a circular pipe (also known as Poiseuille flow); (5) a laminar flow over a circular cylinder; (6) a subsonic flow over a backward‐facing step at low Reynolds numbers; (7) a transonic flow over the RAE 2822 airfoil; and (8) a supersonic flow around a blunt cylinder‐flare model. The spatial approximation order of the scheme is second order. The mesh convergence of the numerical solution was achieved for all cases. The accuracy order for highly compressible flows with discontinuities is close to first order and, for incompressible viscous flows, it is close to second order. The proposed solver is named rhoPimpleCentralFoam and is implemented in the open‐source CFD library OpenFOAM^®. For high speed flows, it shows a similar behavior as the KT/KNP schemes (implemented as rhoCentralFoam‐solver, Int. J. Numer. Meth. Fluids 2010), and for flows with small Mach numbers, it behaves like solvers that are based on the PISO/SIMPLE algorithm.     

A note on reservoir simulation for heterogeneous porous media

We present a front-tracking method for solving two-phase reservoir simulation problems arising in reservoirs with varying geological properties. The method preserves exactly the characteristic features of saturation fronts crossing media discontinuities. Two test examples are presented.     

A fast Godunov method for the water‐hammer problem

An efficient Godunov‐type numerical method with second‐order accuracy was developed to simulate the water‐hammer problem in piping. The exact solutions of the Riemann problem were analysed and illustrated on the intriguing solution diagram by properly introducing dimensionless variables within reasonably practical ranges. Based on the solution diagram, an efficient fast Riemann solver was also developed. Moreover, small perturbation analysis was performed to demonstrate the relations between the primitive variables, velocity and pressure, for the Riemann problem. The typical shock‐tube problem and the water‐hammer problem were implemented as sample ones to test the numerical method. It was shown that the present numerical method incorporated with Van Leer's flux limiter is a promising one to simulate fluid transient problem for piping in the present study. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

Riemann solvers with evolved initial conditions

The scope of this paper is three fold. We first formulate upwind and symmetric schemes for hyperbolic equations with non‐conservative terms. Then we propose upwind numerical schemes for conservative and non‐conservative systems, based on a Riemann solver, the initial conditions of which are evolved non‐linearly in time, prior to a simple linearization that leads to closed‐form solutions. The Riemann solver is easily applied to complicated hyperbolic systems. Finally, as an example, we formulate conservative schemes for the three‐dimensional Euler equations for general compressible materials and give numerical results for a variety of test problems for ideal gases in one and two space dimensions. Copyright © 2006 John Wiley & Sons, Ltd.     

An exact Riemann solver for detonation products

Numerical methods based upon the Riemann Problem are considered for solving the general initial-value problem for the Euler equations applied to real gases. Most of such methods use an approximate solution of the Riemann problem when real gases are involved. These approximate Riemann solvers do not yield always a good resolution of the flow field, especially for contact surfaces and expansion waves. Moreover, approximate Riemann solvers cannot produce exact solutions for the boundary points. In order to overcome these shortcomings, an exact solution of the Riemann problem is developed, valid for real gases. The method is applied to detonation products obeying a 5th order virial equation of state, in the shock-tube test case. Comparisons between our solver, as implemented in Random Choice Method, and finite difference methods, which do not employ a Riemann solver, are given.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Construction of approximate solutions of boundary value impact strain problems

The method for constructing approximate solutions of boundary value problems of impact strain dynamics in the form of ray expansions behind the strain discontinuity fronts is generalized to the case of curvilinear and diverging rays. This proposed generalization is illustrated by an example of dynamics of an antiplane motion of an elastic medium. The ray method is one of the methods for constructing approximate solutions of nonstationary boundary value problems of strain dynamics. It was proposed in [1, 2] and then widely used in nonstationary problems of mathematical physics involving surfaces on which the desired function or its derivatives have discontinuities [3–7]. A complete, qualified survey of papers in this direction can be found in [8]. This method is based on the expansion of the solution in a Taylor-type series behind the moving discontinuity surface rather than in a neighborhood of a stationary point. The coefficients of this series are the jumps of the derivatives of the unknown functions, for which, as a consequence of the compatibility conditions, one can obtain ordinary differential equations, i.e., discontinuity damping equations. In the case where the problem with velocity discontinuity surfaces is considered in a nonlinear medium, this method cannot be used directly, because one cannot obtain the damping equation. A modification of this method for the purpose of using it to solve problems of that type was proposed in [9–11], where, as an example, the solutions of several one-dimensional problems were considered. In the present paper, we show how this method can be transferred to the case of multidimensional impact strain problems in which the geometry of the ray is not known in advance and the rays become curvilinear and diverging. By way of example, we consider a simple problem on the antiplane motion of a nonlinearly elastic incompressible medium.     

An Approximate Roe-Type Riemann Solver for a Class of Realizable Second Order Closures

A realizable, objective second-moment turbulence closure, allowing for an entropy characterisation, is analyzed with respect to its convective subset. The distinct characteristic wave system of these equations in non-conservation form is exposed. An approximate solution to Ihe associated one-dimensional Riemann problem is constructed making use of approximate jump conditions obtained by assuming a linear path across shock waves. A numerical integration method based on a new approximate Riemann solver (flux-difference-splitting) is proposed for use in conjunction with either unstructured or structured grids. Test calculations of quasi one-dimensional flow cases demonstrate the feasibility of the current technique even where Euler-based approaches fail.