Abstract:Strong stability preserving (SSP) Runge-Kutta methods are often desired when evolving in time problems that have two components that have very different time scales. Where the SSP property is needed, it has been shown that implicit and implicit-explicit methods have very restrictive time-steps and are therefore not efficient. For this reason, SSP integrating factor methods may offer an attractive alternative to traditional time-stepping methods for problems with a linear component that is stiff and a nonlinear component that is not. However, the strong stability properties of integrating factor Runge-Kutta methods have not been established. In this work we show that it is possible to define explicit integrating factor Runge-Kutta methods that preserve the desired strong stability properties satisfied by each of the two components when coupled with forward Euler time-stepping, or even given weaker conditions. We define sufficient conditions for an explicit integrating factor Runge--Kutta method to be SSP, namely that they are based on explicit SSP Runge--Kutta methods with non-decreasing abscissas. We find such methods of up to fourth order and up to ten stages, analyze their SSP coefficients, and prove their optimality in a few cases. We test these methods to demonstrate their convergence and to show that the SSP time-step predicted by the theory is generally sharp, and that the non-decreasing abscissa condition is needed in our test cases. Finally, we show that on typical total variation diminishing linear and nonlinear test-cases our new explicit SSP integrating factor Runge-Kutta methods out-perform the corresponding explicit SSP Runge-Kutta methods, implicit-explicit SSP Runge--Kutta methods, and some well-known exponential time differencing methods.

New embedded pairs of explicit Runge–Kutta methods with FSAL properties adapted to the numerical integration of oscillatory problems

A New Pair of Explicit ARKN Methods for the Numerical Integration of General Perturbed Oscillators

Extended RKN Methods with FSAL Property for Oscillatory Systems.

A trigonometrically fitted explicit Numerov-type method for second-order initial value problems with oscillating solutions

Extended RKN-type Methods for Numerical Integration of Perturbed Oscillators

A trigonometrically fitted explicit hybrid method for the numerical integration of orbital problems

Embedded Exponential Runge–Kutta–Nyström Methods for Highly Oscillatory Hamiltonian Systems

Strong Stability Preserving Integrating Factor Runge-Kutta Methods

Functionally-fitted energy-preserving methods for solving oscillatory nonlinear Hamiltonian systems

Two-step Extended RKN Methods for Oscillatory Systems.

A Phase-Fitted Runge-Kutta-Nyström method for the Numerical Solution of Initial Value Problems with Oscillating Solutions

Cost-reduction implicit exponential Runge-Kutta methods for highly oscillatory systems

Some numerical algorithms for solving the highly oscillatory second-order initial value problems.

Two new classes of exponential Runge-Kutta integrators for efficiently solving stiff systems or highly oscillatory problems

Trigonometrically Fitted Explicit Numerov-type Method for Periodic IVPs with Two Frequencies

A new numerical algorithm based on least squares method for solving stochastic Itô-Volterra integral equations

A Multistep Method for Integration of Perturbed and Damped Second-Order ODE Systems

A Robust Trigonometrically Fitted Embedded Pair for Perturbed Oscillators

Explicit Multi-Symplectic Extended Leap-Frog Methods for Hamiltonian Wave Equations

Extended RKN-type methods with minimal dispersion error for perturbed oscillators

Semi-analytical Exponential RKN Integrators for Efficiently Solving High-Dimensional Nonlinear Wave Equations Based on FFT Techniques.