Modified trigonometric integrators

DOI10.1137/130921118zbMATH Open1320.65194arXiv1305.3216OpenAlexW3102638609MaRDI QIDQ5254142FDOQ5254142

Authors: Ari Stern, Robert I. McLachlan 

Publication date: 9 June 2015

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

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.

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

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

zbMATH Keywords

IMEXhighly oscillatory problemsgeometric integratorsFermi-Pasta-Ulamtrigonometric integrators

Mathematics Subject Classification ID

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

Cited In (16)

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

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)