Abstract:We present an approach for the efficient implementation of self-adjusting multi-rate Runge-Kutta methods and we extend the previously available stability analyses of these methods to the case of an arbitrary number of sub-steps for the active components. We propose a physically motivated model problem that can be used to assess the stability of different multi-rate versions of standard Runge-Kutta methods and the impact of different interpolation methods for the latent variables. Finally, we present the results of several numerical experiments, performed with implementations of the proposed methods in the framework of the \textit{OpenModelica} open-source modelling and simulation software, which demonstrate the efficiency gains deriving from the use of the proposed multi-rate approach for physical modelling problems with multiple time scales.

What problem does this paper attempt to address?

This paper aims to solve the problems of numerical solution efficiency and stability in multi - scale ordinary differential equations (ODEs) systems. Specifically, the paper proposes an adaptive multi - rate Runge - Kutta method and conducts an extended analysis of its stability to be applicable to the active components with an arbitrary number of sub - step lengths. The main contributions of the paper include: 1. **Improved implementation of the adaptive multi - rate Runge - Kutta method**: - The paper proposes a simplified and improved adaptive multi - rate method, which combines the standard time - step adaptive technique. By marking only a small part of the variables as fast variables, the computational cost can be reduced while ensuring that the global time - step meets the given error tolerance. - This method is similar to the strategy used in the adaptive finite - element technique, in which a certain proportion of degrees of freedom can be marked for refinement according to the error indicator value exceeding the given tolerance. 2. **Extension of stability analysis**: - The paper extends the previous stability analysis of multi - rate methods to be applicable to an arbitrary number of sub - step lengths in explicit Runge - Kutta methods (ERK), diagonally implicit Runge - Kutta methods (DIRK), and single - diagonally implicit Runge - Kutta methods (ESDIRK). - To verify the stability of these methods, the paper proposes a physically meaningful model problem to evaluate the influence of different interpolation methods on latent variables and the influence of different numbers of sub - step lengths. 3. **Numerical experiments**: - The paper implements the proposed multi - rate method using the OpenModelica open - source modeling and simulation software and demonstrates its performance through several numerical benchmark problems. - The experimental results show that for physical modeling problems with multiple time scales, using the multi - rate method can significantly improve computational efficiency. ### Main research contents 1. **Background and definition of multi - rate methods**: - The paper first reviews the development history of multi - rate methods, especially the adaptive multi - rate methods that do not require a priori partitioning. - It introduces the basic concepts of multi - rate Runge - Kutta methods, including how to identify fast and slow variables and how to dynamically adjust these variables according to the error indicator during runtime. 2. **Implementation of the adaptive multi - rate method**: - A new multi - rate algorithm is proposed, which combines the standard time - step adaptive technique. - Through error estimation and sorting, the fast variables that need to be subdivided into sub - step lengths are determined, thereby reducing the computational cost while ensuring accuracy. 3. **Stability analysis**: - A detailed analysis of the stability of multi - rate methods is carried out, considering different types of Runge - Kutta methods (ERK, DIRK, ESDIRK). - A specific physical model problem is proposed to evaluate the influence of different interpolation methods and the number of sub - step lengths on the stability of multi - rate methods. 4. **Numerical experiments and results**: - The proposed multi - rate method is implemented using OpenModelica software, and its effectiveness is verified through multiple numerical experiments. - The experimental results show that the multi - rate method can significantly improve computational efficiency when dealing with physical modeling problems with multiple time scales. ### Conclusion Through the improved adaptive multi - rate Runge - Kutta method and its stability analysis, the paper provides an effective tool for solving the numerical solution problems of multi - scale ordinary differential equations systems. The experimental results show that this method can significantly reduce the computational cost while maintaining accuracy, especially when dealing with complex physical models with multiple time scales. Future research directions may include further optimizing the implementation of multi - rate methods and applying them to a wider range of physical and engineering problems.

Multi-rate Runge-Kutta methods: stability analysis and applications

Many-Stage Optimal Stabilized Runge–Kutta Methods for Hyperbolic Partial Differential Equations

Many-Stage Optimal Stabilized Runge-Kutta Methods for Hyperbolic Partial Differential Equations

Stability of multistep Runge-Kutta methods for systems of functional-differential and functional equations

Performance analysis of relaxation Runge–Kutta methods

Explicit stabilized multirate method for stiff differential equations

A class of generalized additive Runge-Kutta methods

Construction and mean-square stability analysis of a new family of stochastic Runge-Kutta methods

Runge-Kutta methods: some historical notes

Entropy-Preserving and Entropy-Stable Relaxation IMEX and Multirate Time-Stepping Methods

High order strong stability preserving multi-derivative implicit and IMEX Runge--Kutta methods with asymptotic preserving properties

Algebraically Stable SDIRK Methods with Controllable Numerical Dissipation for First/second-Order Time-Dependent Problems

An introduction to “Almost Runge-Kutta” methods

Development of a new Runge‐Kutta method and its economical implementation

Runge–Kutta methods for a semi-analytical prediction of milling stability

Diagonally drift-implicit Runge-Kutta methods of weak order one and two for Itô SDEs and stability analysis

Three-stage Semi-implicit Stochastic Runge-Kutta Methods for Stochastic Differential Equations

Strong Stability Preserving Integrating Factor Runge-Kutta Methods

A novel high-order linearly implicit and energy-stable additive Runge-Kutta methods for gradient flow models

Some new Runge-Kutta methods with interpolation properties and their application to the Magnetic-Binary Problem

Stability Analysis and Error Estimate of the Explicit Single-Step Time-Marching Discontinuous Galerkin Methods with Stage-Dependent Numerical Flux Parameters for a Linear Hyperbolic Equation in One Dimension