A Well-Balanced HLL Scheme for Hyperbolic Conservation Systems With Source Terms北大核心CSCD

针对含源项的双曲守恒方程给出了一种新的有限体积格式.经典的有限体积格式不能正确地模拟对流通量项和外力之间的平衡所产生的动力学问题.为解决这个问题,仿照经典的HLL近似Riemann求解器设计思路设计了含源项的近似Riemann求解器.针对含重力源项的一维流体Euler方程和理想磁流体方程,通过对通量计算格式的修正得到了保平衡HLL格式（WB-HLL）,并给出了保平衡的证明.针对一维Euler方程和理想磁流体给出了两个算例,比较了传统HLL格式和提出的WB-HLL格式的计算精度.计算结果表明,WB-HLL格式精度更高,收敛更快.     

Fine structures for the solutions of the two-dimensional Riemann problems by high-order WENO schemes

The two-dimensional Riemann problem with polytropic gas is considered. By a restriction on the constant states of each quadrant of the computational domain such that there is only one planar centered wave connecting two adjacent quadrants, there are nineteen genuinely different initial configurations of the problem. The configurations are numerically simulated on a fine grid and compared by the 5th-order WENO-Z5, 6th-order WENO-??6, and 7th-order WENO-Z7 schemes. The solutions are very well approximated with high resolution of waves interactions phenomena and different types of Mach shock reflections. Kelvin-Helmholtz instability-like secondary-scaled vortices along contact continuities are well resolved and visualized. Numerical solutions show that WENO-??6 outperforms the comparing WENO-Z5 and WENO-Z7 in terms of shock capturing and small-scaled vortices resolution. A catalog of the numerical solutions of all nineteen configurations obtained from the WENO-??6 scheme is listed. Thanks to their excellent resolution and sharp shock capturing, the numerical solutions presented in this work can be served as reference solutions for both future numerical and theoretical analyses of the 2D Riemann problem.     

Application of a multi-component Eulerian approach to capture interface evolution and vorticity generation in compressible multiphase flow

The hyperbolic Eularian model is used as a mathematical framework for compressible multiphase flows. The formulation was obtained after an averaging process of the single phase Navier-Stokes equations. The closure of multi-component system leads to the volume fraction equation containing a non-conservative term and a pressure equilibrium condition. As a result the model equations cannot be written in a conservative form. To solve the equations a finite volume Godunov type computational approach is developed which uses an approximate Riemann solver combined with a numerical scheme to tackle the non-conservative terms. The approach accounts for pressure non-equilibrium. It enables resolving interfaces separating compressible fluids and captures the baroclinic source of vorticity generation. The computations are performed for various initial conditions and compared with theoretical and experimental data for a shock-bubble interaction problem. The investigated cases include acoustic wave transmission through isolated bubbles of helium and krypton. The numerical results illustrate the characteristic features of the evolving interfaces. The impulsively generated flow perturbations are dominated by the reflection and refraction of the shock and by the vorticity generation within the media. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A roe-type Riemann solver for hyperbolic systems with relaxation based on time-dependent wave decomposition

Summary. The paper is devoted to the construction of a higher order Roe-type numerical scheme for the solution of hyperbolic systems with relaxation source terms. It is important for applications that the numerical scheme handles both stiff and non stiff source terms with the same accuracy and computational cost and that the relaxation variables are computed accurately in the stiff case. The method is based on the solution of a Riemann problem for a linear system with constant coefficients: a study of the behavior of the solutions of both the nonlinear and linearized problems as the relaxation time tends to zero enables to choose a convenient linearization such that the numerical scheme is consistent with both the hyperbolic system when the source terms are absent and the correct relaxation system when the relaxation time tends to zero. The method is applied to the study of the propagation of sound waves in a two-phase medium. The comparison between our numerical scheme, usual fractional step methods, and numerical simulation of the relaxation system shows the necessity of using the solutions of a fully coupled hyperbolic system with relaxation terms as the basis of a numerical scheme to obtain accurate solutions regardless of the stiffness. Received October 7, 1994 / Revised version received September 27, 1995     

Nonconservative scheme with the isentropic condition in rarefaction waves as applied to the compressible Euler equations

A previously developed second-order accurate quasi-monotone scheme is tested using the Riemann problem with high initial pressure and density ratios. For shock waves, the scheme is conservative, while, in rarefaction waves, the isentropic condition along the trajectory of a Lagrangian particle is used instead of conservativeness in energy. It is shown that the shock front position produced by the scheme has no considerable errors typical of a representative set of conservative quasi-monotone schemes of various orders of accuracy. The numerical accuracy is significantly improved in the case of moving grids with a contact discontinuity explicitly introduced in the form of a grid node. It is shown how the method can be extended to cover the multidimensional case and the presence of additional terms in the original equations.     

Numerical solution of the two-dimensional shallow water equations by the application of relaxation methods

A generalization and extension of a finite difference method for calculating numerical solutions of the two dimensional shallow water system of equations is investigated. A previously developed non-oscillatory relaxation scheme is generalized as to included problems with source terms in two dimensions, with emphasis given to the bed topography, resulting to a class of methods of first- and second-order in space and time. The methods are based on classical relaxation models combined with TVD Runge–Kutta time stepping mechanisms where neither Riemann solvers nor characteristic decompositions are needed. Numerical results are presented for several test problems with or without the source term present. The wetting and drying process is also considered. The presented schemes are verified by comparing the results with documented ones.     

Well-balanced central finite volume methods for the Ripa system

We propose a new well-balanced central finite volume scheme for the Ripa system both in one and two space dimensions. The Ripa system is a nonhomogeneous hyperbolic system with a non-zero source term that is obtained from the shallow water equations system by incorporating horizontal temperature gradients. The proposed numerical scheme is a second-order accurate finite volume method that evolves a non-oscillatory numerical solution on a single grid, avoids the process of solving Riemann problems arising at the cell interfaces, and follows a well-balanced discretization that ensures the steady state requirement by discretizing the geometrical source term according to the discretization of the flux terms. Furthermore the proposed scheme mimics the surface gradient method and discretizes the water height according to the discretization of the water level. The proposed scheme is then applied and classical one and two-dimensional Ripa problems with flat or variable bottom topographies are successfully solved. The obtained numerical results are in good agreement with corresponding ones appearing in the recent literature, thus confirming the potential and efficiency of the proposed method.     

Analysis of Nonlinear Resonance in Conservation Laws with Point Sources and Well-Balanced Scheme

The purpose of this paper is to study the wave behavior of the hyperbolic conservation law with concatenation of point sources: for i ∈ I some finite index set, and where δ( ) is the Dirac measure. Special features of this problem are the discontinuities that appear along the t -axis at the point sources x = x _{i +1/2}. Resonance occurs when the speed of the nonlinear wave is close to zero. In addition to classical shock waves, the equation exhibits overcompressive and undercompressive waves. The Riemann problem "with resonance" is solved, and we show global existence via the Glimm scheme. Analytical understanding is used to design a well-balanced numerical scheme, of the Godunov type, which preserves the balance between the sources terms and the fluxes terms. Some numerical tests are reported.     

A new composite scheme for two-layer shallow water flows with shocks

This paper is devoted to solve the system of partial differential equations governing the flow of two superposed immiscible layers of shallow water flows. The system contains source terms due to bottom topography, wind stresses, and nonconservative products describing momentum exchange between the layers. The presence of these terms in the flow model forms a nonconservative system which is only conditionally hyperbolic. In addition, two-layer shallow water flows are often accompanied with moving discontinuities and shocks. Developing stable numerical methods for this class of problems presents a challenge in the field of computational hydraulics. To overcome these difficulties, a new composite scheme is proposed. The scheme consists of a time-splitting operator where in the first step the homogeneous system of the governing equations is solved using an approximate Riemann solver. In the second step a finite volume method is used to update the solution. To remove the non-physical oscillations in the vicinity of shocks a nonlinear filter is applied. The method is well-balanced, non-oscillatory and it is suitable for both low and high values of the density ratio between the two layers. Several standard test examples for two-layer shallow water flows are used to verify high accuracy and good resolution properties for smooth and discontinuous solutions.     

Riemann problem and elementary wave interactions in isentropic magnetogasdynamics

In this paper, we consider the Riemann problem for a quasilinear hyperbolic system of equations governing the one dimensional unsteady simple wave flow of an isentropic, inviscid and perfectly conducting compressible fluid, subjected to a transverse magnetic field. This class of equations includes, as a special case, the equations of isentropic gasdynamics. We study the shock and rarefaction waves and their properties, and discuss the geometry of shock curves using the Riemann invariant coordinates. Under certain conditions, we show the existence and uniqueness of the solution to the Riemann problem for arbitrary initial data, and then discuss the vacuum state in isentropic magnetogasdynamics. Finally, we discuss numerical results for different initial data, and discuss all possible interactions of elementary waves. It is noticed that although the magnetogasdynamic system is more complex than the corresponding gasdynamic system, all the parallel results remain identical. However, unlike the ordinary gasdynamic case, the solution inside rarefaction waves in magnetogasdynamics cannot be obtained directly and explicitly; indeed, it requires an extra iteration procedure. It is also observed that the presence of a magnetic field makes both the shock and rarefaction stronger compared to what they would have been in the absence of a magnetic field.     

An Efficient Method for Computing Hyperbolic Systems with Geometrical Source Terms Having Concentrations

We propose a simple numerical method for calculating both unsteady and steady state solution of hyperbolic system with geometrical source terms having concentrations. Physical problems under consideration include the shallow water equations with topography,and the quasi one-dimensional nozzle flows. We use the interface value, rather than the cell-averages, for the source terms, which results in a well-balanced scheme that can capture the steady state solution with a remarkable accuracy. This method approximates the source terms via the numerical fluxes produced by an (approximate) Riemann solver for the homogeneous hyperbolic systems with slight additional computation complexity using Newton‘s iterations and numerical integrations. This method solves well the subor super-critical flows, and with a transonic fix, also handles well the transonic flows over the concentration. Numerical examples provide strong evidence on the effectiveness of this new method for both unsteady and steady state calculations.     

Well-balanced unstaggered central schemes for one and two-dimensional shallow water equation systems

We propose a new well-balanced unstaggered central finite volume scheme for hyperbolic balance laws with geometrical source terms. In particular we construct a new one and two-dimensional finite volume method for the numerical solution of shallow water equations on flat/variable bottom topographies. The proposed scheme evolves a non-oscillatory numerical solution on a single grid, avoids the time consuming process of solving Riemann problems arising at the cell interfaces, and is second-order accurate both in space and time. Furthermore, the numerical scheme follows a well-balanced discretization that first discretizes the geometrical source term according to the discretization of the flux terms, and then mimics the surface gradient method and discretizes the water height according to the discretization of the water level. The resulting scheme exactly satisfies the C-property at the discrete level. The proposed scheme is then applied and classical one and two-dimensional shallow water equation problems with flat or variable bottom topographies are successfully solved. The obtained numerical results are in good agreement with corresponding ones appearing in the recent literature, thus confirming the potential and efficiency of the proposed method.     

Global entropy solutions to a class of quasi-linear wave equations with large time-oscillating sources

This research explores the Cauchy problem for a class of quasi-linear wave equations with time dependent sources. It can be transformed into the Cauchy problem of hyperbolic integro-differential systems of nonlinear balance laws. We introduce the generalized Glimm scheme in new version and study its stability which is proved by Glimm-type interaction estimates in a dissipativity assumption. The generalized solutions to the perturbed Riemann problems, the building blocks of generalized Glimm scheme, are constructed by Riemann problem method modeled on the source free equations. The global existence for the Lipschitz continuous solutions and weak solutions to the systems is established by the consistency of scheme and the weak convergence of source. Finally, the weak solutions are also the entropy solutions which satisfy the entropy inequality.     

An asymptotic preserving scheme for hydrodynamics radiative transfer models

In this paper, we present a numerical scheme for a hydrodynamics radiative transfer model consisting of two steps: the first one is based on a relaxation method and the second one on the well balanced scheme. The derivation of the scheme relies on the resolution of a stationary Riemann problem with source terms. The obtained scheme preserves the limited flux property and it is compatible with the diffusive regime of hydrodynamics radiative transfer models. These properties are illustrated by numerical tests, one of them involves a radiative transfer model coupled with an equation for the temperature of the material.     

A hybrid numerical method for the shallow water equations with source term

In coastal oceanography there is interest in problems modeled by the shallow water equations, where variations in channel depth are accounted for by the presence of source terms. A numerical treatment for the solution of such problems is presented here, in terms of a hybrid approach, which combines a second-order TVD scheme for conservation law equations (assuming no source terms) with an eigenvector projection scheme that incorporates the effects of nonzero source terms (in regions where the bottom is not flat). For the case where an initially sharp wave profile is assumed, the progress of a wave as it traverses an estuary whose channel depth varies is calculated. Excellent numerical results are obtained. © 1996 John Wiley & Sons, Inc.     

一维气液两相漂移模型全隐式AUSMV算法研究

气液两相漂移模型显式AUSMV（advection upstream splitting method combined with flux vector splitting method）算法的时间步长受限于CFL（Courant-Friedrichs-Lewy）条件，为了提高计算效率，建立了一种全隐式AUSMV算法求解气液两相漂移模型.采用AUSM格式结合FVS（flux vector splitting）格式构造连续方程和运动方程的对流项数值通量，AUSM格式构造压力项数值通量.离散控制方程是非线性方程组，采用六阶Newton(牛顿)法结合数值Jacobi矩阵求解.计算经典算例Zuber-Findlay激波管问题和复杂漂移关系变质量流动问题，结果分析表明:全隐式AUSMV算法，色散效应小，无数值震荡，计算精度高.在压力波波速高的条件下，可以显著提高计算效率，耗散效应小.     

Riemann Problems for 5×5 Systems of Fully Non-linear Equations Related to Hypoplasticity

The equations of motion for two-dimensional deformations of an incompressible elastoplastic material involve five equations, two equations expressing conservation of momentum, and three constitutive laws, which we take in the rate form, i.e. relating the stress rate to the strain rate. In hypoplasticity, the constitutive laws are homogeneous of degree one in the stress and strain rates. This property has the consequence that although the equations are not in conservation form, there is nonetheless a natural way to characterize planar shock waves. The Riemann problem is the initial value problem for plane waves, in which the initial data for stress and velocity consist of two constant vectors separated by a single discontinuity. The main result is that, under appropriate assumptions, the Riemann problem has a scale invariant piecewise constant solution. The issue of uniqueness is left unresolved. Indeed, we give an example satisfying the conditions for existence, for which there are many solutions. Using asymptotics, we show how solutions of the Riemann problem are approximated by smooth solutions of a system regularized by the addition of viscous terms that preserve the property of scale invariance.     

A Global Solution to a Two-dimensional Riemann Problem Involving Shocks as Free Boundaries

We present a global solution to a Riemann problem for the pressure gradient system of equations.The Riemann problem has initially two shock waves and two contact discontinuities. The angle between the two shock waves is set initially to be close to 180 degrees. The solution has a shock wave that is usually regarded as a free boundary in the self-similar variable plane. Our main contribution in methodology is handling the tangential oblique derivative boundary values.     

A Finite Volume Scheme for Savage-Hutter Equations on Unstructured Grids

A Godunov-type finite volume scheme on unstructured grids is proposed to numerically solve the Savage-Hutter equations in curvilinear coordinate. We show the direct observation that the model isn't a Galilean invariant system. At the cell boundary, the modified Harten-Lax-van Leer (HLL) approximate Riemann solver is adopted to calculate the numerical flux. The modified HLL flux is not troubled by the lack of Galilean invariance of the model and it is helpful to handle discontinuities at free interface. Rigidly the system is not always a hyperbolic system due to the dependence of flux on the velocity gradient. Even so, our numerical results still show quite good agreements with reference solutions. The simulations for granular avalanche flows with shock waves indicate that the scheme is applicable.     

Numerical modelling of free-surface shallow flows over irregular topography with complex geometry

A well-balanced Godunov-type finite volume algorithm is developed for modelling free-surface shallow flows over irregular topography with complex geometry. The algorithm is based on a new formulation of the classical shallow water equations in hyperbolic conservation form. Unstructured triangular grids are used to achieve the adaptability of the grid to the geometry of the problem and to facilitate localised refinement. The numerical fluxes are calculated using HLLC approximate Riemann solver, and the MUSCL-Hancock predictor–corrector scheme is adopted to achieve the second-order accuracy both in space and in time where the solutions are continuous, and to achieve high-resolution results where the solutions are discontinuous. The novelties of the algorithm include preserving well-balanced property without any additional correction terms and the wet/dry front treatments. The good performance of the algorithm is demonstrated by comparing numerical and theoretical results of several benchmark problems, including the preservation of still water over a two-dimensional hump, the idealised dam-break flow over a frictionless flat rectangular channel, the circular dam-break, and the shock wave from oblique wall. Besides, two laboratory dam-break cases are used for model validation. Furthermore, a practical application related to dam-break flood wave propagation over highly irregular topography with complex geometry is presented. The results show that the algorithm can correctly account for free-surface shallow flows with respect to its effectiveness and robustness thus has bright application prospects.