Implicit-explicit variational integration of highly oscillatory problems

DOI10.1137/080732936MaRDI QIDQ3653216zbMATH OpenOpenAlexFDO

Authors Ari Stern, E. Grinspun

Publication date 21 December 2009

Published in Multiscale Modeling & Simulation (Search for Journal in Brave)

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

zbMATH Keywords

numerical experimentsHamiltonian systemsmechanical systemsimplicit-explicit methodssymplectic integratorsvariational integratorsFermi-Pasta-Ulam problemsplitting methodshighly oscillatory problemsdifferent time scalesnumerical resonance instabilityStörmer/Verlet/leapfrog method

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)

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.

Recommendations

• Modified trigonometric integrators
• Symplectic schemes for highly oscillatory Hamiltonian systems with varying fast frequencies
• Instabilities and inaccuracies in the integration of highly oscillatory problems
• Integrators for highly oscillatory Hamiltonian systems: an homogenization approach
• Long-Time-Step Methods for Oscillatory Differential Equations

Cited in

(35)

• Computational dynamics of a 3D elastic string pendulum attached to a rigid body and an inertially fixed reel mechanism
• High order variational integrators in the optimal control of mechanical systems
• Stable IMEX schemes for a Nitsche-based approximation of elastodynamic contact problems. Selective mass scaling interpretation
• Error analysis through energy minimization and stability properties of exponential integrators
• Construction and analysis of higher order Galerkin variational integrators
• General techniques for constructing variational integrators
• The averaged Lagrangian method
• Variational integrators for orbital problems using frequency estimation
• A long-term numerical energy-preserving analysis of symmetric and/or symplectic extended RKN integrators for efficiently solving highly oscillatory Hamiltonian systems
• Analysis of higher order phase fitted variational integrators
• High order numerical methods for highly oscillatory problems
• Long-time oscillatory energy conservation of total energy-preserving methods for highly oscillatory Hamiltonian systems
• Structure-preserving numerical integrators for Hodgkin-Huxley-type systems
• A variational approach to multirate integration for constrained systems
• Long-term analysis of numerical integrators for oscillatory Hamiltonian systems under minimal non-resonance conditions
• Long time energy and kinetic energy conservations of exponential integrators for highly oscillatory conservative systems
• H^1-parametrizations of complex planar physical domains in isogeometric analysis
• Hamiltonian interpolation of splitting approximations for nonlinear PDEs
• 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
• Variational integrators for the dynamics of thermo-elastic solids with finite speed thermal waves
• Numerical integrators for motion under a strong constraining force
• Instabilities and inaccuracies in the integration of highly oscillatory problems
• Modified trigonometric integrators
• Implicit-explicit-compact methods for advection diffusion reaction equations
• Efficient exponential time integration for simulating nonlinear coupled oscillators
• Nearly periodic maps and geometric integration of noncanonical Hamiltonian systems
• Lagrangian and Hamiltonian Taylor variational integrators
• Exponential Rosenbrock methods and their application in visual computing
• Variational time integrators for finite-dimensional thermo-elasto-dynamics without heat conduction
• Exponential integrators for stiff elastodynamic problems
• An implicit asynchronous variational integrator for flexible multibody dynamics
• Variational and linearly implicit integrators, with applications

This page was built for publication: Implicit-explicit variational integration of highly oscillatory problems

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