A Class of Multirate Infinitesimal GARK Methods

Abstract

Differential equations arising in many practical applications are characterized by multiple time scales. Multirate time integration seeks to solve them efficiently by discretizing each scale with a different, appropriate time step, while ensuring the overall accuracy and stability of the numerical solution. In a seminal paper Knoth and Wolke (APNUM, 1998) proposed a hybrid solution approach: discretize the slow component with an explicit Runge-Kutta method, and advance the fast component via a modified fast differential equation. The idea led to the development of multirate infinitesimal step (MIS) methods by Wensch et al. (BIT, 2009.)Günther and Sandu (BIT, 2016) explained MIS schemes as a particular case of multirate General-structure Additive Runge-Kutta (MR-GARK) methods. The hybrid approach offers extreme flexibility in the choice of the numerical solution process for the fast component. This work constructs a family of multirate infinitesimal GARK schemes (MRI-GARK) that extends the hybrid dynamics approachin multiple ways. Order conditions theory and stability analyses are developed, and practical explicit and implicit methods of up to order four are constructed. Numerical results confirm the theoretical findings. We expect the new MRI-GARK family to be most useful for systems of equations with widely disparate time scales, where the fast process is dispersive, and where the influence of the fast component on the slow dynamics is weak.

Cited by in corpus (15)

• Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers
• ARKODE: a flexible IVP solver infrastructure for one-step methods
• Implicit multirate GARK methods
• A DG-IMEX method for two-moment neutrino transport: Nonlinear solvers for neutrino-matter coupling
• Implicit-explicit multirate infinitesimal GARK methods
• Implicit-Explicit Multirate Infinitesimal Stage-Restart Methods
• A fast time-stepping strategy for dynamical systems equipped with a surrogate model
• Multirate Exponential Rosenbrock Methods
• Advancing Simulations of Coupled Electron and Phonon Nonequilibrium Dynamics Using Adaptive and Multirate Time Integration
• Mass-Conserving Implicit-Explicit Methods for Coupled Compressible Navier-Stokes Equations
• Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks
• Moment-based adaptive time integration for thermal radiation transport
• Coupled Multirate Infinitesimal GARK Schemes for Stiff Systems with Multiple Time Scales
• A unified formulation of splitting-based implicit time integration schemes
• New Time Integrators and Capabilities in SUNDIALS Versions 6.2.0-7.4.0