One-sided stability and convergence of the Nessyahu–Tadmor scheme

Non-oscillatory schemes are widely used in numerical approximations of nonlinear conservation laws. The Nessyahu–Tadmor (NT) scheme is an example of a second order scheme that is both robust and simple. In this paper, we prove a new stability property of the NT scheme based on the standard minmod reconstruction in the case of a scalar strictly convex conservation law. This property is similar to the One-sided Lipschitz condition for first order schemes. Using this new stability, we derive the convergence of the NT scheme to the exact entropy solution without imposing any nonhomogeneous limitations on the method. We also derive an error estimate for monotone initial data.     

A non-parameterized entropy correction for Roe's approximate Riemann solver

Summary. The main drawback with Roe's approximate Riemann solver is that non-physical expansion shocks can occur in the vicinity of sonic points. Previous work aimed at enforcing the entropy condition is based on the representation of sonic rarefaction waves. We propose a new non-parameterized approach which is based on a nonlinear Hermite interpolation of an approximate flux function and the exact resolution of non convex scalar Riemann problems. Convergence and consistency with the entropy condition are proved for scalar convex conservation laws with arbitrarily large initial data. When considering strictly hyperbolic systems of conservation laws, consistency of the resulting scheme with the entropy condition is also proved for initial data sufficiently close to a constant. Numerical results on a one-dimensional shock-tube and a two-dimensional supersonic forward facing step confirm our theoretical results. Received March 1, 1993 / Revised version received August 26, 1994     

A priori error estimates for approximate solutions to convex conservation laws

Summary. We introduce a new technique for proving a priori error estimates between the entropy weak solution of a scalar conservation law and a finite–difference approximation calculated with the scheme of Engquist-Osher, Lax-Friedrichs, or Godunov. This technique is a discrete counterpart of the duality technique introduced by Tadmor [SIAM J. Numer. Anal. 1991]. The error is related to the consistency error of cell averages of the entropy weak solution. This consistency error can be estimated by exploiting a regularity structure of the entropy weak solution. One ends up with optimal error estimates. Received December 21, 2001 / Revised version received February 18, 2002 / Published online June 17, 2002     

Error estimate for the upwind finite volume method for the nonlinear scalar conservation law

In this paper we estimate the error of upwind first order finite volume schemes applied to scalar conservation laws. As a first step, we consider standard upwind and flux finite volume scheme discretization of a linear equation with space variable coefficients in conservation form. We prove that, in spite of their lack of consistency, both schemes lead to a first order error estimate. As a final step, we prove a similar estimate for the nonlinear case. Our proofs rely on the notion of geometric corrector, introduced in our previous paper by Bouche et al. (2005) [24] in the context of constant coefficient linear advection equations.     

Entropy production in second-order three-point schemes

Summary We discuss semi-discrete three-point finite difference methods for the numerical solution of system of conservation laws which are second order accurate in space in the sense of truncation error. Particular discretizations of the numerical entropy flux associated with such schemes are studied clarifying the importance of this discretization with regard to the production of numerical entropy. Using a numerical entropy flux constructed in a canonical way we prove that a wide class of finite difference methods cannot satisfy a discrete entropy inequality. Together with a well known result of Schonbek concerning Lax-Wendroff type schemes our result indicates a strong relationship between entropy production and oscillations in numerical solutions.The research reported here was supported by a grant from the Stiftung Volkswagenwerk, Federal Republic of Germany. It is a part of the doctoral thesis of the above author, Universität Stuttgart, 1991.     

An entropy satisfying MUSCL scheme for systems of conservation laws

Summary. For the high-order numerical approximation of hyperbolic systems of conservation laws, we propose to use as a building principle an entropy diminishing criterion instead of the familiar total variation diminishing criterion introduced by Harten for scalar equations. Based on this new criterion, we derive entropy diminishing projections that ensure, both, the second order of accuracy and all of the classical discrete entropy inequalities. The resulting scheme is a nonlinear version of the classical Van Leer's MUSCL scheme. Strong convergence of this second order, entropy satisfying scheme is proved for systems of two equations. Numerical tests demonstrate the interest of our theory. Received March 28, 1995 / Revised version received June 17, 1995     

Convergence of a staggered Lax-Friedrichs scheme for nonlinear conservation laws on unstructured two-dimensional grids

Summary. Based on Nessyahu and Tadmor's nonoscillatory central difference schemes for one-dimensional hyperbolic conservation laws [16], for higher dimensions several finite volume extensions and numerical results on structured and unstructured grids have been presented. The experiments show the wide applicability of these multidimensional schemes. The theoretical arguments which support this are some maximum-principles and a convergence proof in the scalar linear case. A general proof of convergence, as obtained for the original one-dimensional NT-schemes, does not exist for any of the extensions to multidimensional nonlinear problems. For the finite volume extension on two-dimensional unstructured grids introduced by Arminjon and Viallon [3,4] we present a proof of convergence for the first order scheme in case of a nonlinear scalar hyperbolic conservation law. Received April 8, 2000 / Published online December 19, 2000     

Second order scheme for scalar conservation laws with discontinuous flux

Burger, Karlsen, Torres and Towers in [9] proposed a flux TVD (FTVD) second order scheme with Engquist–Osher flux, by using a new nonlocal limiter algorithm for scalar conservation laws with discontinuous flux modeling clarifier thickener units. In this work we show that their idea can be used to construct FTVD second order scheme for general fluxes like Godunov, Engquist–Osher, Lax–Friedrich, … satisfying (A, B)-interface entropy condition for a scalar conservation law with discontinuous flux with proper modification at the interface. Also corresponding convergence analysis is shown. We show further from numerical experiments that solutions obtained from these schemes are comparable with the second order schemes obtained from the minimod limiter.     

Convergence of MUSCL and filtered schemes for scalar conservation laws and Hamilton-Jacobi equations

Summary. This paper considers the questions of convergence of: (i) MUSCL type (i.e. second-order, TVD) finite-difference approximations towards the entropic weak solution of scalar, one-dimensional conservation laws with strictly convex flux and (ii) higher-order schemes (filtered to ``preserve' an upper-bound on some weak second-order finite differences) towards the viscosity solution of scalar, multi-dimensional Hamilton-Jacobi equations with convex Hamiltonians. Received May 16, 1994     

A MUSCL method satisfying all the numerical entropy inequalities

We consider here second-order finite volume methods for one-dimensional scalar conservation laws. We give a method to determine a slope reconstruction satisfying all the exact numerical entropy inequalities. It avoids inhomogeneous slope limitations and, at least, gives a convergence rate of . It is obtained by a theory of second-order entropic projections involving values at the nodes of the grid and a variant of error estimates, which also gives new results for the first-order Engquist-Osher scheme.

     

A finite volume method for scalar conservation laws with stochastic time–space dependent flux functions

We propose a new finite volume method for scalar conservation laws with stochastic time–space dependent flux functions. The stochastic effects appear in the flux function and can be interpreted as a random manner to localize the discontinuity in the time–space dependent flux function. The location of the interface between the fluxes can be obtained by solving a system of stochastic differential equations for the velocity fluctuation and displacement variable. In this paper we develop a modified Rusanov method for the reconstruction of numerical fluxes in the finite volume discretization. To solve the system of stochastic differential equations for the interface we apply a second-order Runge–Kutta scheme. Numerical results are presented for stochastic problems in traffic flow and two-phase flow applications. It is found that the proposed finite volume method offers a robust and accurate approach for solving scalar conservation laws with stochastic time–space dependent flux functions.     

Third order nonoscillatory central scheme for hyperbolic conservation laws

Summary. A third-order accurate Godunov-type scheme for the approximate solution of hyperbolic systems of conservation laws is presented. Its two main ingredients include: 1. A non-oscillatory piecewise-quadratic reconstruction of pointvalues from their given cell averages; and 2. A central differencing based on staggered evolution of the reconstructed cell averages. This results in a third-order central scheme, an extension along the lines of the second-order central scheme of Nessyahu and Tadmor \cite{NT}. The scalar scheme is non-oscillatory (and hence – convergent), in the sense that it does not increase the number of initial extrema (– as does the exact entropy solution operator). Extension to systems is carried out by componentwise application of the scalar framework. In particular, we have the advantage that, unlike upwind schemes, no (approximate) Riemann solvers, field-by-field characteristic decompositions, etc., are required. Numerical experiments confirm the high-resolution content of the proposed scheme. Thus, a considerable amount of simplicity and robustness is gained while retaining the expected third-order resolution. Received April 10, 1996 / Revised version received January 20, 1997     

间断流函数守恒律方程的高精度有限差分格式

本文研究了具有间断流函数的守恒律方程,借助本质无振荡（ENO）的思想,利用Rankine—Hugoniot关系和全局熵条件设计出一种高精度计算格式;并利用此格式计算出相关情形的Riemann问题,显示了满意的数值解果．     

A Large Time Step explicit scheme (CFL>1) on unstructured grids for 2D conservation laws: Application to the homogeneous shallow water equations

In this work, a Large Time Step (LTS) explicit finite volume scheme designed to allow CFL > 1 is applied to the numerical resolution of 2D scalar and systems of conservation laws on triangular grids. Based on the flux difference splitting formulation, a special concern is put on finding the way of packing the information to compute the numerical solution when working on unstructured grids. Not only the cell areas but also the length of the interfaces and their orientation are questions of interest to send the information from each edge or interface. The information to update the cell variables is computed according to the local average discrete velocity and the orientation of the edges of the cells involved. The performance of these ideas is tested and compared with the conventional explicit first order and second order schemes in academic configurations for the 2D linear scalar equation and for 2D systems of conservation laws (in particular the shallow water equations) without source terms. The LTS scheme is demonstrated to preserve or even gain accuracy and save computational time with respect to the first order scheme.     

On the convergence of a second order approximation of conservation laws with discontinuous flux

A second order scheme is constructed for the scalar conservation laws with flux function allowed to be discontinuous in the space variable. The corresponding numerical solutionis shown to converge to the (A,B) entropysolution. Numerical results are presented to validate the convergence analysis. In addition, it is inferred that these numerical solutions are comparable with those obtained using the standard minmod limiter.     

A parallel well-balanced finite volume method for shallow water equations with topography on the cubed-sphere

A finite volume scheme for the global shallow water model on the cubed-sphere mesh is proposed and studied in this paper. The new cell-centered scheme is based on Osher’s Riemann solver together with a high-order spatial reconstruction. On each patch interface of the cubed-sphere only one layer of ghost cells is needed in the scheme and the numerical flux is calculated symmetrically across the interface to ensure the numerical conservation of total mass. The discretization of the topographic term in the equation is properly modified in a well-balanced manner to suppress spurious oscillations when the bottom topography is non-smooth. Numerical results for several test cases including a steady-state nonlinear geostrophic flow and a zonal flow over an isolated mountain are provided to show the flexibility of the scheme. Some parallel implementation details as well as some performance results on a parallel supercomputer with more than one thousand processor cores are also provided.     

Semidiscrete central-upwind scheme for conservation laws with a discontinuous flux function in space

In this paper, a modified semidiscrete central-upwind scheme is derived for the scalar conservation laws with a discontinuous flux function in space. The new scheme is based on dealing with the phase transition at the stationary discontinuity, where the unknown variable function is not continuous, but the flux function is continuous. The main advantages of the new scheme are the same as them of the original semidiscrete central-upwind scheme. Numerical results are displayed to illustrate the efficiency of the methods.     

The method of fractional steps for conservation laws

Summary The stability, accuracy, and convergence of the basic fractional step algorithms are analyzed when these algorithms are used to compute discontinuous solutions of scalar conservation laws. In particular, it is proved that both first order splitting and Strang splitting algorithms always converge to the unique weak solution satisfying the entropy condition. Examples of discontinuous solutions are presented where both Strang-type splitting algorithms are only first order accurate but one of the standard first order algorithms is infinite order accurate. Various aspects of the accuracy, convergence, and correct entropy production are also studied when each split step is discretized via monotone schemes, Lax-Wendroff schemes, and the Glimm scheme.Partially supported by an Alfred Sloan Foundation fellowship and N.S.F. grant MCS-76-10227Sponsored by US Army under contract No. DAA 629-75-0-0024     

A third-order semi-discrete genuinely multidimensional central scheme for hyperbolic conservation laws and related problems

Summary. We construct a new third-order semi-discrete genuinely multidimensional central scheme for systems of conservation laws and related convection-diffusion equations. This construction is based on a multidimensional extension of the idea, introduced in [17] – the use of more precise information about the local speeds of propagation, and integration over nonuniform control volumes, which contain Riemann fans. As in the one-dimensional case, the small numerical dissipation, which is independent of , allows us to pass to a limit as . This results in a particularly simple genuinely multidimensional semi-discrete scheme. The high resolution of the proposed scheme is ensured by the new two-dimensional piecewise quadratic non-oscillatory reconstruction. First, we introduce a less dissipative modification of the reconstruction, proposed in [29]. Then, we generalize it for the computation of the two-dimensional numerical fluxes. Our scheme enjoys the main advantage of the Godunov-type central schemes –simplicity, namely it does not employ Riemann solvers and characteristic decomposition. This makes it a universal method, which can be easily implemented to a wide variety of problems. In this paper, the developed scheme is applied to the Euler equations of gas dynamics, a convection-diffusion equation with strongly degenerate diffusion, the incompressible Euler and Navier-Stokes equations. These numerical experiments demonstrate the desired accuracy and high resolution of our scheme. Received February 7, 2000 / Published online December 19, 2000