Error bounds of finite difference schemes for multi-dimensional scalar conservation laws with source terms

This paper studies explicit and semi=implicit finite differenceschemes approximating non-homogeneous multi-dimensional scalarconservation laws with both stiff and non-stiff source terms.Error bounds of order O(t) are presented for both cases.     

Diffusion-dispersion limits for multidimensional scalar conservation laws with source terms

In this paper we consider conservation laws with diffusion and dispersion terms. We study the convergence for approximation applied to conservation laws with source terms. The proof is based on the Hwang and Tzavaras's new approach [Seok Hwang, Athanasios E. Tzavaras, Kinetic decomposition of approximate solutions to conservation laws: Application to relaxation and diffusion-dispersion approximations, Comm. Partial Differential Equations 27 (5-6) (2002) 1229-1254] and the kinetic formulation developed by Lions, Perthame, and Tadmor [P.-L. Lions, B. Perthame, E. Tadmor, A kinetic formulation of multidimensional scalar conservation laws and related equations, J. Amer. Math. Soc. 7 (1) (1994) 169-191].     

Convergence of relaxation schemes for hyperbolic conservation laws with stiff source terms

We focus in this study on the convergence of a class of relaxation numerical schemes for hyperbolic scalar conservation laws including stiff source terms. Following Jin and Xin, we use as approximation of the scalar conservation law, a semi-linear hyperbolic system with a second stiff source term. This allows us to avoid the use of a Riemann solver in the construction of the numerical schemes. The convergence of the approximate solution toward a weak solution is established in the cases of first and second order accurate MUSCL relaxed methods.

     

On stability of difference schemes. Central schemes for hyperbolic conservation laws with source terms

The stability of nonlinear explicit difference schemes with not, in general, open domains of the scheme operators are studied. For the case of path-connected, bounded, and Lipschitz domains, we establish the notion that a multi-level nonlinear explicit scheme is stable iff (if and only if) the corresponding scheme in variations is stable. A new modification of the central Lax–Friedrichs (LxF) scheme is developed to be of the second-order accuracy. The modified scheme is based on nonstaggered grids. A monotone piecewise cubic interpolation is used in the central scheme to give an accurate approximation for the model in question. The stability of the modified scheme is investigated. Some versions of the modified scheme are tested on several conservation laws, and the scheme is found to be accurate and robust. As applied to hyperbolic conservation laws with, in general, stiff source terms, it is constructed a second-order nonstaggered central scheme based on operator-splitting techniques.     

Implicit kinetic schemes for scalar conservation laws

Based on kinetic formulation for scalar conservation laws, we present implicit kinetic schemes. For time stepping these schemes require resolution of linear systems of algebraic equations. The scheme is conservative at steady states. We prove that if time marching procedure converges to some steady state solution, then the implicit kinetic scheme converges to some entropy steady state solution. We give sufficient condition of the convergence of time marching procedure. For scalar conservation laws with a stiff source term we construct a stiff numerical scheme with discontinuous artificial viscosity coefficients that ensure the scheme to be equilibrium conserving. We couple the developed implicit approach with the stiff space discretization, thus providing improved stability and equilibrium conservation property in the resulting scheme. Numerical results demonstrate high computational capabilities (stability for large CFL numbers, fast convergence, accuracy) of the developed implicit approach. © 2002 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 18: 26–43, 2002     

Weakly nonoscillatory schemes for scalar conservation laws

A new class of Godunov-type numerical methods (called here weakly nonoscillatory or WNO) for solving nonlinear scalar conservation laws in one space dimension is introduced. This new class generalizes the classical nonoscillatory schemes. In particular, it contains modified versions of Min-Mod and UNO. Under certain conditions, convergence and error estimates for WNO methods are proved.

     

On convergence of numerical schemes for hyperbolic conservation laws with stiff source terms

We deal in this study with the convergence of a class of numerical schemes for scalar conservation laws including stiff source terms. We suppose that the source term is dissipative but it is not necessarily a Lipschitzian function. The convergence of the approximate solution towards the entropy solution is established for first and second order accurate MUSCL and for splitting semi-implicit methods.

     

Equilibrium schemes for scalar conservation laws with stiff sources

We consider a simple model case of stiff source terms in hyperbolic conservation laws, namely, the case of scalar conservation laws with a zeroth order source with low regularity. It is well known that a direct treatment of the source term by finite volume schemes gives unsatisfactory results for both the reduced CFL condition and refined meshes required because of the lack of accuracy on equilibrium states. The source term should be taken into account in the upwinding and discretized at the nodes of the grid. In order to solve numerically the problem, we introduce a so-called equilibrium schemes with the properties that (i) the maximum principle holds true; (ii) discrete entropy inequalities are satisfied; (iii) steady state solutions of the problem are maintained. One of the difficulties in studying the convergence is that there are no estimates for this problem. We therefore introduce a kinetic interpretation of upwinding taking into account the source terms. Based on the kinetic formulation we give a new convergence proof that only uses property (ii) in order to ensure desired compactness framework for a family of approximate solutions and that relies on minimal assumptions. The computational efficiency of our equilibrium schemes is demonstrated by numerical tests that show that, in comparison with an usual upwind scheme, the corresponding equilibrium version is far more accurate. Furthermore, numerical computations show that equilibrium schemes enable us to treat efficiently the sources with singularities and oscillating coefficients.

     

Monotone first‐order weighted schemes for scalar conservation laws

In this work, we present a monotone first‐order weighted (FORWE) method for scalar conservation laws using a variational formulation. We prove theoretical properties as consistency, monotonicity, and convergence of the proposed scheme for the one‐dimensional (1D) Cauchy problem. These convergence results are extended to multidimensional scalar conservation laws by a dimensional splitting technique. For the validation of the FORWE method, we consider some standard bench‐mark tests of bidimensional and 1D conservation law equations. Finally, we analyze the accuracy of the method with L¹ and L^∞ error estimates. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Finite volume schemes for nonhomogeneous scalar conservation laws: error estimate

Summary. In this paper, we study finite volume schemes for the nonhomogeneous scalar conservation law with initial condition . The source term may be either stiff or nonstiff. In both cases, we prove error estimates between the approximate solution given by a finite volume scheme (the scheme is totally explicit in the nonstiff case, semi-implicit in the stiff case) and the entropy solution. The order of these estimates is in space-time -norm (h denotes the size of the mesh). Furthermore, the error estimate does not depend on the stiffness of the source term in the stiff case. Received October 21, 1999 / Published online February 5, 2001     

Finite-difference methods for one-dimensional hyperbolic conservation laws

This article contains a survey of some important finite-difference methods for one-dimensional hyperbolic conservation laws. Weak solutions of hyperbolic conservation laws are introduced and the concept of entropy stability is discussed. Furthermore, the Riemann problem for hyperbolic conservation laws is solved. An introduction to finite-difference methods is given for which important concepts such as, e.g., conservativity, stability, and consistency are introduced. Godunov-type methods are elaborated for general systems of hyperbolic conservation laws. Finally, flux limiter methods are developed for the scalar nonlinear conservation law. © 1994 John Wiley & Sons, Inc.     

Convergence of fluctuation-splitting schemes for two dimensional scalar conservation laws with a kinetic solver

Summary. The “fluctuation-splitting schemes” (FSS in short) have been introduced by Roe and Sildikover to solve advection equations on rectangular grids and then extended to triangular grids by Roe, Deconinck, Struij... For a two dimensional nonlinear scalar conservation law, we consider the case of a triangular grid and of a kinetic approach to reduce the discretization of the nonlinear equation to a linear equation and apply a particular FSS called N-scheme. We show that the resulting scheme converges strongly in in a finite volume sense. Received February 25, 1997 / Revised version received November 8, 1999 / Published online August 24, 2000     

A study of splitting scheme for hyperbolic conservation laws with source terms

This study deals with the convergence of a numerical scheme for conservation laws including source terms. A splitting method for source term integration is presented. More precisely, the convergence of the numerical solution towards the entropy solution is proved in the scalar case. Because of the effect of source term, the constructed scheme is total variation bounded. Numerical experiments for one-dimensional shallow water equation are presented to demonstrate the performance of the scheme.     

Convergence of finite volume monotone schemes for scalar conservation laws on bounded domains

Summary. This paper is devoted to the study of the finite volume methods used in the discretization of conservation laws defined on bounded domains. General assumptions are made on the data: the initial condition and the boundary condition are supposed to be measurable bounded functions. Using a generalized notion of solution to the continuous problem (namely the notion of entropy process solution, see [9]) and a uniqueness result on this solution, we prove that the numerical solution converges to the entropy weak solution of the continuous problem in for every . This also yields a new proof of the existence of an entropy weak solution. Received May 18, 2000 / Revised version received November 21, 2000 / Published online June 7, 2001     

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     

Stochastic scalar conservation laws

We introduce a notion of stochastic entropic solution à la Kruzkov, but with Ito's calculus replacing deterministic calculus. This results in a rich family of stochastic inequalities defining what we mean by a solution. A uniqueness theory is then developed following a stochastic generalization of L¹ contraction estimate. An existence theory is also developed by adapting compensated compactness arguments to stochastic setting. We use approximating models of vanishing viscosity solution type for the construction. While the uniqueness result applies to any spatial dimensions, the existence result, in the absence of special structural assumptions, is restricted to one spatial dimension only.     

Multiresolution schemes for conservation laws

Summary. In recent years a variety of high–order schemes for the numerical solution of conservation laws has been developed. In general, these numerical methods involve expensive flux evaluations in order to resolve discontinuities accurately. But in large parts of the flow domain the solution is smooth. Hence in these regions an unexpensive finite difference scheme suffices. In order to reduce the number of expensive flux evaluations we employ a multiresolution strategy which is similar in spirit to an approach that has been proposed by A. Harten several years ago. Concrete ingredients of this methodology have been described so far essentially for problems in a single space dimension. In order to realize such concepts for problems with several spatial dimensions and boundary fitted meshes essential deviations from previous investigations appear to be necessary though. This concerns handling the more complex interrelations of fluxes across cell interfaces, the derivation of appropriate evolution equations for multiscale representations of cell averages, stability and convergence, quantifying the compression effects by suitable adapted multiscale transformations and last but not least laying grounds for ultimately avoiding the storage of data corresponding to a full global mesh for the highest level of resolution. The objective of this paper is to develop such ingredients for any spatial dimension and block structured meshes obtained as parametric images of Cartesian grids. We conclude with some numerical results for the two–dimensional Euler equations modeling hypersonic flow around a blunt body. Received June 24, 1998 / Revised version received February 21, 2000 / Published online November 8, 2000     

Localization effects and measure source terms in numerical schemes for balance laws

This paper investigates the behavior of numerical schemes for nonlinear conservation laws with source terms. We concentrate on two significant examples: relaxation approximations and genuinely nonhomogeneous scalar laws. The main tool in our analysis is the extensive use of weak limits and nonconservative products which allow us to describe accurately the operations achieved in practice when using Riemann-based numerical schemes. Some illustrative and relevant computational results are provided.

     

Third-order scalar evolution equations with conservation laws

Associated to the family of third-order quasilinear scalar evolution equations is the geometry of point transformations. This geometry provides a framework from which to study the structure of conservation laws of the equation, and to study the special nature of the geometry of those equations which do possess conservation laws. There is an easy and obvious normal form for equations which possess at least one conservation law. The geometric structure of the equation gives rise to a simple yet much less obvious normal form for equations which possess at least two conservation laws.