Backward error analysis for multi-symplectic integration methods

Summary.A useful method for understanding discretization error in the numerical solution of ODEs is to compare the system of ODEs with the modified equations obtained through backward error analysis, and using symplectic integration for Hamiltonian ODEs provides more incite into the modified equations. In this paper, the ideas of symplectic integration are extended to Hamiltonian PDEs, and this paves the way for the development of a local modified equation analysis solely as a useful diagnostic tool for the study of these types of discretizations. In particular, local conservation laws of energy and momentum are not preserved exactly when symplectic integrators are used to discretize, but the modified equations are used to derive modified conservation laws that are preserved to higher order along the numerical solution. These results are also applied to the nonlinear wave equation.

[1]  L. Vu-Quoc,et al.  Finite difference calculus invariant structure of a class of algorithms for the nonlinear Klein-Gordon equation , 1995 .

[2]  R. F. Warming,et al.  The modified equation approach to the stability and accuracy analysis of finite-difference methods , 1974 .

[3]  Thomas J. Bridges,et al.  Unstable eigenvalues and the linearization about solitary waves and fronts with symmetry , 1999, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[4]  Brian E. Moore,et al.  Multi-symplectic integration methods for Hamiltonian PDEs , 2003, Future Gener. Comput. Syst..

[5]  Catherine Sulem,et al.  The nonlinear Schrödinger equation , 2012 .

[6]  Bengt Fornberg,et al.  A practical guide to pseudospectral methods: Introduction , 1996 .

[7]  A. Durán,et al.  The numerical integration of relative equilibrium solutions. Geometric theory , 1998 .

[8]  R. McLachlan.  Symplectic integration of Hamiltonian wave equations , 1993 .

[9]  G. Benettin,et al.  On the Hamiltonian interpolation of near-to-the identity symplectic mappings with application to symplectic integration algorithms , 1994 .

[10]  T. Bridges.  Multi-symplectic structures and wave propagation , 1997, Mathematical Proceedings of the Cambridge Philosophical Society.

[11]  A. Durán,et al.  The numerical integration of relative equilibrium solutions. The nonlinear Schrödinger equation , 2000 .

[12]  J. Marsden,et al.  Multisymplectic Geometry, Variational Integrators, and Nonlinear PDEs , 1998, math/9807080.

[13]  T. Bridges.  A geometric formulation of the conservation of wave action and its implications for signature and the classification of instabilities , 1997, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[14]  Salvador Jiménez,et al.  Derivation of the discrete conservation laws for a family of finite difference schemes , 1994 .

[15]  Benedict Leimkuhler,et al.  Computational Molecular Dynamics: Challenges, Methods, Ideas , 1999, Computational Molecular Dynamics.

[16]  Robert L. Pego,et al.  Solitary waves on FPU lattices: I. Qualitative properties, renormalization and continuum limit , 1999 .

[17]  Brian E. Moore.  A modified equations approach for multi-symplectic integration methods. , 2003 .

[18]  Ernst Hairer,et al.  The life-span of backward error analysis for numerical integrators , 1997 .

[19]  S. Reich.  Backward Error Analysis for Numerical Integrators , 1999 .

[20]  W. Strauss,et al.  Numerical solution of a nonlinear Klein-Gordon equation , 1978 .

[21]  S. Reich.  Multi-Symplectic Runge—Kutta Collocation Methods for Hamiltonian Wave Equations , 2000 .

[22]  J. C. Simo,et al.  On the stability of symplectic and energy-momentum algorithms for non-linear Hamiltonian systems with symmetry , 1996 .

[23]  J. M. Sanz-Serna,et al.  Numerical Hamiltonian Problems , 1994 .

[24]  L. Vu-Quoc,et al.  INVARIANT-CONSERVING FINITE DIFFERENCE ALGORITHMS FOR THE NONLINEAR KLEIN-GORDON EQUATION , 1993 .

[25]  S. Reich,et al.  Multi-symplectic integrators: numerical schemes for Hamiltonian PDEs that conserve symplecticity , 2001 .

[26]  O. Gonzalez.  Time integration and discrete Hamiltonian systems , 1996 .