Local Discontinuous Galerkin Methods for the Two-Dimensional Camassa–Holm Equation

Abstract

In this paper, the local discontinuous Galerkin method is developed to solve the two-dimensional Camassa–Holm equation in rectangular meshes. The idea of LDG methods is to suitably rewrite a higher-order partial differential equations into a first-order system, then apply the discontinuous Galerkin method to the system. A key ingredient for the success of such methods is the correct design of interface numerical fluxes. The energy stability for general solutions of the method is proved. Comparing with the Camassa–Holm equation in one-dimensional case, there are more auxiliary variables which are introduced to handle high-order derivative terms. The proof of the stability is more complicated. The resulting scheme is high-order accuracy and flexible for arbitrary h and p adaptivity. Different types of numerical simulations are provided to illustrate the accuracy and stability of the method.

Keywords

Local discontinuous Galerkin method / Two-dimensional Camassa–Holm equation / Stability

Funding

National Natural Science Foundation of China(11722112)