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

Kristian Debrabant,Andreas Rößler
DOI: https://doi.org/10.1016/j.apnum.2008.03.011
2016-05-09
Abstract:The class of stochastic Runge-Kutta methods for stochastic differential equations due to Rößler is considered. Coefficient families of diagonally drift-implicit stochastic Runge-Kutta (DDISRK) methods of weak order one and two are calculated. Their asymptotic stability as well as mean-square stability (MS-stability) properties are studied for a linear stochastic test equation with multiplicative noise. The stability functions for the DDISRK methods are determined and their domains of stability are compared to the corresponding domain of stability of the considered test equation. Stability regions are presented for various coefficients of the families of DDISRK methods in order to determine step size restrictions such that the numerical approximation reproduces the characteristics of the solution process.
Numerical Analysis
What problem does this paper attempt to address?
The problems that this paper attempts to solve are: **Develop and analyze the diagonally - drifted implicit Runge - Kutta (DDISRK) methods for weak first - order and second - order approximations of Ito stochastic differential equations (SDEs), and study the asymptotic stability and mean - square stability of these methods**. Specifically, the authors Kristian Debrabant and Andreas Rößler aim to: 1. **Compute coefficient families**: Compute the coefficient families of weak first - order and second - order for the diagonally - drifted implicit stochastic Runge - Kutta (DDISRK) method. 2. **Stability analysis**: Study the asymptotic stability and mean - square stability (MS - stability) of these methods when dealing with linear stochastic test equations with multiplicative noise. 3. **Stability function and region comparison**: Determine the stability function of the DDISRK method and compare its stability region with the corresponding stability region of the test equation to determine the step - size limitations for the numerical approximation to reproduce the characteristics of the solution process. ### Background and motivation Stochastic differential equations (SDEs) are widely used in the modeling of finance, physics, biology, etc. However, many SDEs do not have analytical solutions, so numerical methods are required to solve them. For the numerical solutions of SDEs, stability is a very important property, which determines the reliability of the numerical method in long - time simulations. Traditional explicit methods may be unstable when dealing with stiff problems, so implicit methods (such as the diagonally - drifted implicit method) are proposed to improve stability. This paper focuses on the diagonally - drifted implicit Runge - Kutta method. This type of method can better handle SDEs with complex noise structures while maintaining high precision. ### Main contributions 1. **Develop new numerical methods**: Propose the diagonally - drifted implicit Runge - Kutta methods of weak first - order and second - order, which are suitable for Ito - type stochastic differential equations. 2. **Stability theory**: Analyze in detail the asymptotic stability and mean - square stability of these methods and give specific stability conditions. 3. **Numerical experiment verification**: Verify the validity of the theoretical results through numerical experiments and show the advantages of the new method when dealing with SDEs of high - dimensional driving Wiener processes. In summary, the main objective of this paper is to develop and analyze efficient numerical methods to ensure both high precision and good stability when solving stochastic differential equations.