Optimal energy-conserving local discontinuous Galerkin method for the one-dimensional sine-Gordon equation

ABSTRACT

The nonlinear sine-Gordon equation arises in various problems in science and engineering. In this paper, we propose and analyze a high-order and energy-conserving local discontinuous Galerkin (LDG) method for the sine-Gordon nonlinear hyperbolic equation in one space dimension. We prove the energy-conserving property and the $L^2$ stability for the semi-discrete LDG method. Optimal a priori error estimates for the solution and for the auxiliary variable that approximates the first-order derivative are derived in the $L^2$-norm for the semi-discrete formulation. In particular, we identify a special numerical flux and a particular projection of the initial conditions for the LDG scheme for which the $L^2$-norm of the solution and its spatial derivative are of order $p+1$, when piecewise polynomials of degree at most p are used. Our numerical experiments demonstrate optimal order of convergence. Several numerical results are presented to validate the theoretical analyze of the proposed algorithm. It appears that similar conclusions are valid for the two-dimensional case.

KEYWORDS:

• Sine-Gordon equation
• energy conservation
• $L^2$ stability
• local discontinuous Galerkin method
• a priori error estimates

2010 AMS SUBJECT CLASSIFICATIONS:

• 65M12
• 65M15
• 65M60
• 35Q51
• 35Q53

Acknowledgments

The author would also like to thank the anonymous referees for their constructive comments and remarks which helped improve the quality and readability of the paper.

Disclosure statement

No potential conflict of interest was reported by the only one author.

Additional information

Funding

This research was supported by the University Committee on Research and Creative Activity (UCRCA Proposal 2015-01-F) at the University of Nebraska at Omaha.