We study a class of semi-Lagrangian schemes which can be interpreted as a discrete version of the Hopf-Lax-Oleinik representation formula for the exact viscosity solution of first order evolutive Hamilton-Jacobi equations. That interpretation shows that the scheme is potentially accurate to any pres
New High-Resolution Semi-discrete Central Schemes for Hamilton–Jacobi Equations
✍ Scribed by Alexander Kurganov; Eitan Tadmor
- Publisher
- Elsevier Science
- Year
- 2000
- Tongue
- English
- Weight
- 397 KB
- Volume
- 160
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
✦ Synopsis
We introduce a new high-resolution central scheme for multidimensional Hamilton-Jacobi equations. The scheme retains the simplicity of the non-oscillatory central schemes developed by C.-T. Lin and E. Tadmor (in press, SIAM J. Sci. Comput.), yet it enjoys a smaller amount of numerical viscosity, independent of 1/ t. By letting t ↓ 0 we obtain a new second-order central scheme in the particularly simple semi-discrete form, along the lines of the new semi-discrete central schemes recently introduced by the authors in the context of hyperbolic conservation laws. Fully discrete versions are obtained with appropriate Runge-Kutta solvers. The smaller amount of dissipation enables efficient integration of convection-diffusion equations, where the accumulated error is independent of a small time step dictated by the CFL limitation. The scheme is non-oscillatory thanks to the use of nonlinear limiters. Here we advocate the use of such limiters on second discrete derivatives, which is shown to yield an improved high resolution when compared to the usual limitation of first derivatives. Numerical experiments demonstrate the remarkable resolution obtained by the proposed new central scheme.
📜 SIMILAR VOLUMES
Central schemes may serve as universal finite-difference methods for solving nonlinear convection-diffusion equations in the sense that they are not tied to the specific eigenstructure of the problem, and hence can be implemented in a straightforward manner as black-box solvers for general conservat