Implicit-Explicit Variational Integration of Highly Oscillatory Problems
- 1 January 2009
- journal article
- Published by Society for Industrial & Applied Mathematics (SIAM) in Multiscale Modeling & Simulation
- Vol. 7 (4), 1779-1794
- https://doi.org/10.1137/080732936
Abstract
In this paper, we derive a variational integrator for certain highly oscillatory problems in mechanics. To do this, we take a new approach to the splitting of fast and slow potential forces: rather than splitting these forces at the level of the differential equations or the Hamiltonian, we split the two potentials with respect to the Lagrangian action integral. By using a different quadrature rule to approximate the contribution of each potential to the action, we arrive at a geometric integrator that is implicit in the fast force and explicit in the slow force. This can allow for significantly longer time steps to be taken (compared to standard explicit methods, such as St\"ormer/Verlet) at the cost of only a linear solve rather than a full nonlinear solve. We also analyze the stability of this method, in particular proving that it eliminates the linear resonance instabilities that can arise with explicit multiple-time-stepping methods. Next, we perform some numerical experiments, studying the behavior of this integrator for two test problems: a system of coupled linear oscillators, for which we compare against the resonance behavior of the r-RESPA method; and slow energy exchange in the Fermi--Pasta--Ulam problem, which couples fast linear oscillators with slow nonlinear oscillators. Finally, we prove that this integrator accurately preserves the slow energy exchange between the fast oscillatory components, which explains the numerical behavior observed for the Fermi--Pasta--Ulam problem.Comment: 17 pages, 3 figures, to appear in SIAM Multiscale Modeling and Simulation. v2: added new section, new figure, minor edits; v3: minor revision and clarificatioKeywords
Other Versions
This publication has 14 references indexed in Scilit:
- Variational Integrators for the GravitationalN‐Body ProblemThe Astrophysical Journal, 2007
- Geometric numerical integration illustrated by the Störmer–Verlet methodActa Numerica, 2003
- Asynchronous Variational IntegratorsArchive for Rational Mechanics and Analysis, 2003
- Discrete mechanics and variational integratorsActa Numerica, 2001
- The Midpoint Scheme and Variants for Hamiltonian Systems: Advantages and PitfallsSIAM Journal on Scientific Computing, 1999
- Implicit-explicit Runge-Kutta methods for time-dependent partial differential equationsApplied Numerical Mathematics, 1997
- Dangers of Multiple Time Step MethodsJournal of Computational Physics, 1993
- On the problem of energy equipartition for large systems of the Fermi-Pasta-Ulam type: analytical and numerical estimatesPhysica D: Nonlinear Phenomena, 1992
- Une méthode multipas implicite-explicite pour l'approximation des équations d'évolution paraboliquesNumerische Mathematik, 1980
- A study of extrapolation methods based on multistep schemes without parasitic solutionsZeitschrift für angewandte Mathematik und Physik, 1979