Explicit multi-symplectic extended leap-frog methods for Hamiltonian wave equations

  • Authors:
  • Wei Shi;Xinyuan Wu;Jianlin Xia

  • Affiliations:
  • Department of Mathematics, Nanjing University, State Key Laboratory for Novel Software Technology at Nanjing University, Nanjing 210093, PR China;Department of Mathematics, Nanjing University, State Key Laboratory for Novel Software Technology at Nanjing University, Nanjing 210093, PR China;Department of Mathematics, Purdue University, West Lafayette, IN 47907, USA

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2012

Quantified Score

Hi-index 31.45

Visualization

Abstract

In this paper, we study the integration of Hamiltonian wave equations whose solutions have oscillatory behaviors in time and/or space. We are mainly concerned with the research for multi-symplectic extended Runge-Kutta-Nystrom (ERKN) discretizations and the corresponding discrete conservation laws. We first show that the discretizations to the Hamiltonian wave equations using two symplectic ERKN methods in space and time respectively lead to an explicit multi-symplectic integrator (Eleap-frogI). Then we derive another multi-symplectic discretization using a symplectic ERKN method in time and a symplectic partitioned Runge-Kutta method, which is equivalent to the well-known Stormer-Verlet method in space (Eleap-frogII). These two new multi-symplectic schemes are extensions of the leap-frog method. The numerical stability and dispersive properties of the new schemes are analyzed. Numerical experiments with comparisons are presented, where the two new explicit multi-symplectic methods and the leap-frog method are applied to the linear wave equation and the Sine-Gordon equation. The numerical results confirm the superior performance and some significant advantages of our new integrators in the sense of structure preservation.