Andreas Rößler

Kristian Debrabant, Andreas Rößler

In the present paper, a class of stochastic Runge-Kutta methods containing the second order stochastic Runge-Kutta scheme due to E. Platen for the weak approximation of Itô stochastic differential equation systems with a multi-dimensional Wiener process is considered. Order one and order two conditions for the coefficients of explicit stochastic Runge-Kutta methods are solved and the solution space of the possible coefficients is analyzed. A full classification of the coefficients for such stochastic Runge-Kutta schemes of order one and two with minimal stage numbers is calculated. Further, within the considered class of stochastic Runge-Kutta schemes coefficients for optimal schemes in the sense that additionally some higher order conditions are fulfilled are presented.

Claudine Leonhard, Andreas Rößler

Higher order numerical schemes for stochastic partial differential equations that do not possess commutative noise require the simulation of iterated stochastic integrals. In this work, we extend the algorithms derived by Kloeden, Platen, and Wright (1992) and by Wiktorsson (2001) for the approximation of two-times iterated stochastic integrals involved in numerical schemes for finite dimensional stochastic ordinary differential equations to an infinite dimensional setting. These methods clear the way for new types of approximation schemes for SPDEs without commutative noise. Precisely, we analyze two algorithms to approximate two-times iterated integrals with respect to an infinite dimensional $Q$-Wiener process in case of a trace class operator $Q$ given the increments of the $Q$-Wiener process. Error estimates in the mean-square sense are derived and discussed for both methods. In contrast to the finite dimensional setting, which is contained as a special case, the optimal approximation algorithm cannot be uniquely determined but is dependent on the covariance operator $Q$. This difference arises as the stochastic process is of infinite dimension.

Kristian Debrabant, Andreas Rößler

Recently, a new class of second order Runge-Kutta methods for Itô stochastic differential equations with a multidimensional Wiener process was introduced by Rößler. In contrast to second order methods earlier proposed by other authors, this class has the advantage that the number of function evaluations depends only linearly on the number of Wiener processes and not quadratically. In this paper, we give a full classification of the coefficients of all explicit methods with minimal stage number. Based on this classification, we calculate the coefficients of an extension with minimized error constant of the well-known RK32 method to the stochastic case. For three examples, this method is compared numerically with known order two methods and yields very promising results.

Kristian Debrabant, Andreas Rößler

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.

Claudine Leonhard, Andreas Rößler

We consider a higher-order Milstein scheme for stochastic partial differential equations with trace class noise which fulfill a certain commutativity condition. A novel technique to generally improve the order of convergence of Taylor schemes for stochastic partial differential equations is introduced. The key tool is an efficient approximation of the Milstein term by particularly tailored nested derivative-free terms. For the resulting derivative-free Milstein scheme the computational cost is, in general, considerably reduced by some power. Further, a rigorous computational cost model is considered and the so called effective order of convergence is introduced which allows to directly compare various numerical schemes in terms of their efficiency. As the main result, we prove for a broad class of stochastic partial differential equations, including equations with operators that do not need to be pointwise multiplicative, that the effective order of convergence of the proposed derivative-free Milstein scheme is significantly higher than for the original Milstein scheme. In this case, the derivative-free Milstein scheme outperforms the Euler scheme as well as the original Milstein scheme due to the reduction of the computational cost. Finally, we present some numerical examples that confirm the theoretical results.

Kristian Debrabant, Andreas Rößler

A convergence theorem for the continuous weak approximation of the solution of stochastic differential equations by general one step methods is proved, which is an extension of a theorem due to Milstein. As an application, uniform second order conditions for a class of continuous stochastic Runge-Kutta methods containing the continuous extension of the second order stochastic Runge-Kutta scheme due to Platen are derived. Further, some coefficients for optimal continuous schemes applicable to Itô stochastic differential equations with respect to a multi-dimensional Wiener process are presented.