this paper we introduce a new family of central schemes which retain the simplicity of being independent of the eigenstructure of the problem, yet which enjoy a much smaller numerical viscosity (of the corresponding order )).In particular, our new central schemes maintain their high-resolution independent of O(1/#t ), and letting #t 0, they admit a particularly simple semi-discrete formulation

. 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 thirdorder central scheme, an extension along the lines of the second-order central scheme of Nessyahu and Tadmor [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 highresolution content of...

We introduce new Godunov-type semi-discrete central schemes for hyperbolic systems of conservation laws and Hamilton-Jacobi equations. The schemes are based on the use of more precise information about the local speeds of propagation, and can be viewed as a generalization of the schemes from [26, 24, 25] and [27]. The main advantages of the proposed central schemes are the high resolution, due to the smaller amount of the numerical dissipation, and the simplicity. There are no Riemann solvers and characteristic decomposition involved, and this makes them a universal tool for a wide variety of applications. At the same time, the developed schemes have an upwind nature, since they respect the directions of wave propagation by measuring the one-sided local speeds. This is the reason why we call them central-upwind schemes. The constructed schemes are applied to various problems, such as the Euler equations of gas dynamics, the Hamilton-Jacobi equations with convex and nonconvex Hamiltoni...

In this work we propose a new set of partial differential equations (PDEs) which can be seen as a generalization of the classical eikonal and transport equations, to allow for solutions with multiple phases. The traditional geometrical optics pair of equations suffer from the fact that the class of physically relevant solutions is limited. In particular, it does not include solutions with multiple phases, corresponding to crossing waves. Our objective has been to generalize these equations to accommodate solutions containing more than one phase. The new equations are based on the same high frequency approximation of the scalar wave equation as the eikonal and the transport equations. However, they also incorporate a finite superposition principle. The maximum allowed number of intersecting waves in the solution can be chosen arbitrarily, but a higher number means that a larger system of PDEs must be solved. The PDEs form a hyperbolic system of conservation laws with source terms. Altho...

We present new third- and fifth-order Godunov-type central schemes for approximating solutions of the Hamilton-Jacobi (HJ) equation in an arbitrary number of space dimensions. These are the first central schemes for approximating solutions of the HJ equations with an order of accuracy that is greater than two. In two space dimensions we present two versions for the third-order scheme: one scheme that is based on a genuinely two-dimensional Central WENO reconstruction, and another scheme that is based on a simpler dimension-by-dimension reconstruction. The simpler dimension-by-dimension variant is then extended to a multi-dimensional fifth-order scheme. Our numerical examples in one, two and three space dimensions verify the expected order of accuracy of the schemes. Key words. Hamilton-Jacobi equations, central schemes, high order, WENO, CWENO.

During the recent decades there was an enormous amount of activity related to the construction and analysis of modern algorithms for the approximate solution of nonlinear hyperbolic conservation laws and related problems. To present some aspects of this successful activity, we discuss the analytical tools which are used in the development of convergence theories for these algorithms. These include classical compactness arguments (based on BV a priori estimates), the use of compensated compactness arguments (based on H^-1-compact entropy production), measure valued solutions (measured by their negative entropy production), and finally, we highlight the most recent addition to this bag of analytical tools -- the use of averaging lemmas which yield new compactness and regularity results for nonlinear conservation laws and related equations. We demonstrate how these analytical tools are used in the convergence analysis of approximate solutions for hyperbolic conservation laws and related equations. Our discussion includes examples of Total Variation Diminishing (TVD) finite-difference schemes; error estimates derived from the one-sided stability of Godunov-type methods for convex conservation laws (and their multidimensional analogue -- viscosity solutions of demi-concave Hamilton-Jacobi equations); we outline, in the one-dimensional case, the convergence proof of finite-element streamline-diffusion and spectral viscosity schemes based on the div-curl lemma; we also address the questions of convergence and error estimates for multidimensional finite-volume schemes on non-rectangular grids; and finally, we indicate the convergence of approximate solutions with underlying kinetic formulation, e.g., finite-volume and relaxation schemes, once their regularizing effect is quantified in terms of the averaging lemma.

We present a general procedure to convert schemes which are based on staggered spatial grids into nonstaggered schemes. This procedure is then used to construct a new family of nonstaggered, central schemes for hyperbolic conservation laws by converting the family of staggered central schemes recently introduced in [H. Nessyahu and E. Tadmor, J. Comput. Phys., 87 (1990), pp. 408--463; X. D. Liu and E. Tadmor, Numer. Math., 79 (1998), pp. 397--425; G. S. Jiang and E. Tadmor, SIAM J. Sci. Comput., 19 (1998), pp. 1892--1917]. These new nonstaggered central schemes retain the desirable properties of simplicity and high resolution, and in particular, they yield Riemann-solver-free recipes which avoid dimensional splitting. Most important, the new central schemes avoid staggered grids and hence are simpler to implement in frameworks which involve complex geometries and boundary conditions.

We present a family of high-order, essentially non-oscillatory, central schemes for approximating solutions of hyperbolic systems of conservation laws. These schemes are based on a new centered version of the Weighed Essentially NonOscillatory (WENO) reconstruction of point-values from cell-averages, which is then followed by an accurate approximation of the fluxes via a natural continuous extension of Runge-Kutta solvers. We explicitly construct the third and fourthorder scheme and demonstrate their high-resolution properties in several numerical tests.

In this paper we study the numerical transition from a Hamilton-Jacobi (H-J) equation to its associated system of conservation laws in arbitrary space dimensions. We first study how, in a very generic setting, one can recover the viscosity solutions of the H-J equation using the numerical solutions to the system of conservation laws. We then introduce a simple, second-order relaxation scheme to solve the underlying weakly hyperbolic system of conservation laws.

In this paper, we construct second-order central schemes for multidimensional Hamilton-Jacobi equations and we show that they are nonoscillatory in the sense of satisfying the maximum principle. Thus, these schemes provide the first examples of nonoscillatory second-order Godunov-type schemes based on global projection operators. Numerical experiments are performed; L 1 /L # -errors and convergence rates are calculated. For convex Hamiltonians, numerical evidence confirms that our central schemes converge with second-order rates, when measured in the L 1 -norm advocated in our recent paper [Numer Math, to appear]. The standard L # -norm, however, fails to detect this second-order rate.