Abstract:We note a fact that stiff systems or differential equations that have highly oscillatory solutions cannot be solved efficiently using conventional methods. In this paper, we study two new classes of exponential Runge-Kutta (ERK) integrators for efficiently solving stiff systems or highly oscillatory problems. We first present a novel class of explicit modified version of exponential Runge-Kutta (MVERK) methods based on the order conditions. Furthermore, we consider a class of explicit simplified version of exponential Runge-Kutta (SVERK) methods. Numerical results demonstrate the high efficiency of the explicit MVERK integrators and SVERK methods derived in this paper compared with the well-known explicit ERK integrators for stiff systems or highly oscillatory problems in the literature.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: traditional methods cannot efficiently solve rigid systems or highly oscillatory problems. Specifically, for rigid differential equations or differential equations with highly oscillatory solutions, when using traditional methods (such as standard explicit methods) for numerical solution, a very small time step is required to achieve acceptable accuracy, which leads to a long running time. Therefore, these traditional methods are inefficient in dealing with such problems. To solve this problem, the author proposes two new types of exponential Runge - Kutta (ERK) integrators: 1. **Explicit modified exponential Runge - Kutta (MVERK) method**: - By introducing a matrix - valued function \( w_s(hM) \), the update step is modified while keeping the internal stages unchanged. - The formula is expressed as: \[ \begin{cases} Y_i = y_0 + h \sum_{j = 1}^s \bar{a}_{ij}(MY_j + f(Y_j)), & i = 1,\dots,s, \\ y_1 = e^{hM}y_0 + h \sum_{i = 1}^s \bar{b}_i f(Y_i)+w_s(hM), \end{cases} \] where \( \bar{a}_{ij} \) and \( \bar{b}_i \) are real constants, \( w_s(hM) \) is an appropriate matrix - valued function, and when \( M\rightarrow0 \), \( w_s(hM) = 0 \). 2. **Explicit simplified exponential Runge - Kutta (SVERK) method**: - The matrix exponential of \( hM \) is used in the internal stages, but all coefficients are independent of the matrix exponential. - The formula is expressed as: \[ \begin{cases} Y_i = e^{\bar{c}_i hM}y_0 + h \sum_{j = 1}^s \bar{a}_{ij}f(Y_j), & i = 1,\dots,s, \\ y_1 = e^{hM}y_0 + h \sum_{i = 1}^s \bar{b}_i f(Y_i)+\tilde{w}_s(hM), \end{cases} \] where \( \bar{a}_{ij} \) and \( \bar{b}_i \) are real constants, \( \bar{c}_i=\sum_{j = 1}^s \bar{a}_{ij} \), \( \tilde{w}_s(hM) \) is an appropriate matrix - valued function, and when \( M\rightarrow0 \), \( \tilde{w}_s(hM) = 0 \). By introducing these two new types of ERK integrators, the author aims to reduce the computational cost and improve the efficiency of solving rigid systems or highly oscillatory problems. Numerical experiments show that the newly proposed explicit MVERK and SVERK methods have higher accuracy and efficiency compared with the explicit ERK integrators in the existing literature.

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

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

Two new families of fourth-order explicit exponential Runge--Kutta methods with four stages for first-order differential systems

Error bounds for explicit ERKN integrators for systems of multi-frequency oscillatory second-order differential equations

Explicit Multi-Frequency Symmetric Extended RKN Integrators for Solving Multi-Frequency and Multidimensional Oscillatory Reversible Systems

Symplectic Exponential Runge-Kutta Methods for Solving Nonlinear Hamiltonian Systems

ENERGY-PRESERVING CONTINUOUS-STAGE EXPONENTIAL RUNGE--KUTTA INTEGRATORS FOR EFFICIENTLY SOLVING HAMILTONIAN SYSTEMS

Effective Integrators for Nonlinear Second-Order Oscillatory Systems with a Time-Dependent Frequency Matrix

A new approach to constructing efficient stiffly accurate exponential propagation iterative methods of Runge-Kutta type (EPIRK)

Nonlinear Stability and Convergence of ERKN Integrators for Solving Nonlinear Multi-Frequency Highly Oscillatory Second-Order ODEs with Applications to Semi-Linear Wave Equations

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

Adaptive rational Krylov methods for exponential Runge--Kutta integrators

Adaptive Rational Krylov Methods for Exponential Runge–Kutta Integrators

ERKN Integrators Solving Multi-Frequency Highly Oscillatory Systems with Applications

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

A highly accurate explicit symplectic ERKN method for multi-frequency and multidimensional oscillatory Hamiltonian systems

ERKN integrators for systems of oscillatory second-order differential equations

AnEssential Extension of the Finite-Energy Condition for ExtendedRunge-Kutta-Nyström Integratorswhen Applied to Nonlinear Wave Equations

Efficient Implementation of the ARKN and ERKN Integrators for Multi-Frequency Oscillatory Systems with Multiple Time Scales

ERKN methods for long-term integration of multidimensional orbital problems

Explicit stabilized multirate method for stiff differential equations