Arc length based WENO scheme for Hamilton-Jacobi Equations

In this article, novel smoothness indicators are presented for calculating the nonlinear weights of weighted essentially non-oscillatory scheme to approximate the viscosity numerical solutions of Hamilton-Jacobi equations. These novel smoothness indicators are constructed from the derivatives of reconstructed polynomials over each sub-stencil. The constructed smoothness indicators measure the arc-length of the reconstructed polynomials so that the new nonlinear weights could get less absolute truncation error and gives a high-resolution numerical solution. Extensive numerical tests are conducted and presented to show the performance capability and the numerical accuracy of the proposed scheme with the comparison to the classical WENO scheme.Comment: 14 pages, 9 figure

A New Discontinuous Galerkin Finite Element Method for Directly Solving the Hamilton-Jacobi Equations

In this paper, we improve upon the discontinuous Galerkin (DG) method for Hamilton-Jacobi (HJ) equation with convex Hamiltonians in (Y. Cheng and C.-W. Shu, J. Comput. Phys. 223:398-415,2007) and develop a new DG method for directly solving the general HJ equations. The new method avoids the reconstruction of the solution across elements by utilizing the Roe speed at the cell interface. Besides, we propose an entropy fix by adding penalty terms proportional to the jump of the normal derivative of the numerical solution. The particular form of the entropy fix was inspired by the Harten and Hyman's entropy fix (A. Harten and J. M. Hyman. J. Comput. Phys. 50(2):235-269, 1983) for Roe scheme for the conservation laws. The resulting scheme is compact, simple to implement even on unstructured meshes, and is demonstrated to work for nonconvex Hamiltonians. Benchmark numerical experiments in one dimension and two dimensions are provided to validate the performance of the method

Local-Structure-Preserving Discontinuous Galerkin Methods with Lax-Wendroff Type Time Discretizations for Hamilton-Jacobi Equations

NSF [DMS-0652481, DMS-0847241]; Alfred P. Sloan Research Fellowship; NSFC [10931004, 1081112028, 10871093]In this paper, a family of high order numerical methods are designed to solve the Hamilton-Jacobi equation for the viscosity solution. In particular, the methods start with a hyperbolic conservation law system closely related to the Hamilton-Jacobi equation. The compact one-step one-stage Lax-Wendroff type time discretization is then applied together with the local-structure-preserving discontinuous Galerkin spatial discretization. The resulting methods have lower computational complexity and memory usage on both structured and unstructured meshes compared with some standard numerical methods, while they are capable of capturing the viscosity solutions of Hamilton-Jacobi equations accurately and reliably. A collection of numerical experiments is presented to illustrate the performance of the methods