Modified trigonometric integrators

DOI10.1137/130921118MaRDI QIDQ5254142zbMATH OpenOpenAlexFDO

Authors Ari Stern, Robert I. McLachlan

Publication date 9 June 2015

Published in SIAM Journal on Numerical Analysis (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1305.3216

zbMATH Keywords

Fermi-Pasta-Ulamgeometric integratorshighly oscillatory problemsIMEXtrigonometric integrators

Mathematics Subject Classification ID

Discretization methods and integrators (symplectic, variational, geometric, etc.) for dynamical systems (37M15) Numerical methods for Hamiltonian systems including symplectic integrators (65P10) Systems with slow and fast motions for nonlinear problems in mechanics (70K70)

Abstract: We study modified trigonometric integrators, which generalize the popular class of trigonometric integrators for highly oscillatory Hamiltonian systems by allowing the fast frequencies to be modified. Among all methods of this class, we show that the IMEX (implicit-explicit) method, which is equivalent to applying the midpoint rule to the fast, linear part of the system and the leapfrog (St"ormer/Verlet) method to the slow, nonlinear part, is distinguished by the following properties: (i) it is symplectic; (ii) it is free of artificial resonances; (iii) it is the unique method that correctly captures slow energy exchange to leading order; (iv) it conserves the total energy and a modified oscillatory energy up to to second order; (v) it is uniformly second-order accurate in the slow components; and (vi) it has the correct magnitude of deviations of the fast oscillatory energy, which is an adiabatic invariant. These theoretical results are supported by numerical experiments on the Fermi-Pasta-Ulam problem and indicate that the IMEX method, for these six properties, dominates the class of modified trigonometric integrators.

Recommendations

• Reconsidering trigonometric integrators
• Implicit-explicit variational integration of highly oscillatory problems
• Long-term analysis of numerical integrators for oscillatory Hamiltonian systems under minimal non-resonance conditions
• Conservation properties of numerical integrators for highly oscillatory Hamiltonian systems
• Diagonally implicit trigonometrically fitted symplectic Runge-Kutta methods

Cited in

(16)

• A long-term numerical energy-preserving analysis of symmetric and/or symplectic extended RKN integrators for efficiently solving highly oscillatory Hamiltonian systems
• Symplectic P-stable additive Runge-Kutta methods
• Long-time analysis of an extended RKN integrator for Hamiltonian systems with a solution-dependent high frequency
• Numerical integrators based on the Magnus expansion for nonlinear dynamical systems
• Long-time momentum and actions behaviour of energy-preserving methods for semi-linear wave equations via spatial spectral semi-discretisations
• The averaged Lagrangian method
• Long-term analysis of the Störmer-Verlet method for Hamiltonian systems with a solution-dependent high frequency
• Reconsidering trigonometric integrators
• Implicit-explicit variational integration of highly oscillatory problems
• Dispersive fractalisation in linear and nonlinear Fermi-Pasta-Ulam-Tsingou lattices
• Long-time oscillatory energy conservation of total energy-preserving methods for highly oscillatory Hamiltonian systems
• Long time energy and kinetic energy conservations of exponential integrators for highly oscillatory conservative systems
• Structure-preserving numerical integrators for Hodgkin-Huxley-type systems
• High order numerical methods for highly oscillatory problems
• Numerical conservations of energy, momentum and actions in the full discretisation for nonlinear wave equations
• Long-term analysis of numerical integrators for oscillatory Hamiltonian systems under minimal non-resonance conditions

This page was built for publication: Modified trigonometric integrators

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5254142)