A new symmetric linearly implicit exponential integrator preserving polynomial invariants or Lyapunov functions for conservative or dissipative systems

Abstract

We present a new linearly implicit exponential integrator that preserves the polynomial first integrals or Lyapunov functions for the conservative and dissipative stiff equations, respectively. The method is tested by both oscillated ordinary differential equations and partial differential equations, e.g., an averaged system in wind-induced oscillation, the Fermi-Pasta-Ulam systems, and the polynomial pendulum oscillators. The numerical simulations confirm the conservative properties of the proposed method and demonstrate its good behavior in superior running speed when compared with fully implicit schemes for long-time simulations.

References in corpus (5)

• On energy stable, maximum-principle preserving, second order BDF scheme with variable steps for the Allen-Cahn equation
• Exponential integrators preserving first integrals or Lyapunov functions for conservative or dissipative systems
• A unified approach to Hamiltonian systems, Poisson systems, gradient systems, and systems with Lyapunov functions and/or first integrals
• Analysis of the second order BDF scheme with variable steps for the molecular beam epitaxial model without slope selection
• A linearly implicit energy-preserving exponential integrator for the nonlinear Klein-Gordon equation