Hamiltonian-Preserving Discontinuous Galerkin Methods for the Liouville Equation With Discontinuous Potential

Abstract

. Numerically solving the Liouville equation in classical mechanics with a discontinuous potential often leads to the 4 challenges of how to preserve the Hamiltonian across the potential barrier and a severe time step constraint according to the CFL 5 condition. Motivated by the Hamiltonian-preserving finite volume schemes by Jin and Wen [19], we introduce a Hamiltonian- 6 preserving discontinuous Galerkin (DG) scheme for the Liouville equation with discontinuous potential in this paper. The DG 7 method can be designed with arbitrary order of accuracy, and offers many advantages including easy adaptivity, compact stencils 8 and the ability of handling complicated boundary condition and interfaces. We propose to carefully design the numerical fluxes 9 of the DG methods to build the behavior of a classical particle at the potential barrier into the numerical scheme, which ensures 10 the continuity of the Hamiltonian across the potential barrier and the correct transmission and reflection condition. Our scheme 11 is proved to be positive and stable in L 1 norm if the positivity-preserving limiter is applied. Numerical examples are provided to 12 illustrate the accuracy and effectiveness of the proposed numerical scheme. 13 results show 2D2V test discontinuity of HPDG