Dispersive properties of multisymplectic integrators

  • Authors:

  • C. M. Schober;T. H. Wlodarczyk

  • Affiliations:

  • Department of Mathematics, University of Central Florida, P.O. Box 161364, Orlando 32816-1364, United States;Department of Mathematics, University of Central Florida, P.O. Box 161364, Orlando 32816-1364, United States

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

Quantified Score

Hi-index 31.47

Visualization

Abstract

Multisymplectic (MS) integrators, i.e. numerical schemes which exactly preserve a discrete space-time symplectic structure, are a new class of structure preserving algorithms for solving Hamiltonian PDEs. In this paper we examine the dispersive properties of MS integrators for the linear wave and sine-Gordon equations. In particular a leapfrog in space and time scheme (a member of the Lobatto Runge-Kutta family of methods) and the Preissman box scheme are considered. We find the numerical dispersion relations are monotonic and that the sign of the group velocity is preserved. The group velocity dispersion (GVD) is found to provide significant information and succinctly explain the qualitative differences in the numerical solutions obtained with the different schemes. Further, the numerical dispersion relations for the linearized sine-Gordon equation provides information on the ability of the MS integrators to capture the sine-Gordon dynamics. We are able to link the numerical dispersion relations to the total energy of the various methods, thus providing information on the coarse grid behavior of MS integrators in the nonlinear regime.