Implicit-Explicit Variational Integration of Highly Oscillatory Problems

Abstract

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örmer/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.

Hamiltonian Interpolation of Splitting Approximations for Nonlinear PDEs

“…The "mid-split" cases correspond to the approximation of the system ξ a = −iω a ξ a , a ∈ N , η a = iω a η a , a∈ N , by the midpoint rule [1,28]. Starting from a given point (ξ 0 a , η 0 a ) we have by definition for the first equation…”

A Brief Introduction to Variational Integrators

“…Other applications of variational integration can be reviewed, e.g., in (Kharevych et al 2006;Kraus 2013;Ober-Blöbauma et al 2013;Stern and Grinspun 2009). …”

Variational time integrators for finite‐dimensional thermo‐elasto‐dynamics without heat conduction

“…This example also showcases the conservation of linear momentum. Each spring has the constitutive behavior defined by (31), with C 0 = 5, = 0.2, l 0 = 1 and 0 = 300. Initially, the springs form an equilateral triangle with side length 1, with the masses located at the vertices.…”