Picture water droplets

Developing mathematical theories of the physical world: Open Access research on fluid dynamics from Strathclyde

Strathprints makes available Open Access scholarly outputs by Strathclyde's Department of Mathematics & Statistics, where continuum mechanics and industrial mathematics is a specialism. Such research seeks to understand fluid dynamics, among many other related areas such as liquid crystals and droplet evaporation.

The Department of Mathematics & Statistics also demonstrates expertise in population modelling & epidemiology, stochastic analysis, applied analysis and scientific computing. Access world leading mathematical and statistical Open Access research!

Explore all Strathclyde Open Access research...

Mean-square and asymptotic stability of numerical methods for stochastic ordinary differential equations

Higham, D.J. (2000) Mean-square and asymptotic stability of numerical methods for stochastic ordinary differential equations. SIAM Journal on Numerical Analysis, 38 (3). pp. 753-769. ISSN 0036-1429

Full text not available in this repository. Request a copy from the Strathclyde author

Abstract

Stability analysis of numerical methods for ordinary differential equations (ODEs) is motivated by the question 'for what choices of stepsize does the numerical method reproduce the characteristics of the test equation?' We study a linear test equation with a multiplicative noise term, and consider mean-square and asymptotic stability of a stochastic version of the theta method. We extend some mean-square stability results in [Saito and Mitsui, SIAM. J. Numer. Anal., 33 (1996), pp. 2254--2267]. In particular, we show that an extension of the deterministic A-stability property holds. We also plot mean-square stability regions for the case where the test equation has real parameters. For asymptotic stability, we show that the issue reduces to finding the expected value of a parametrized random variable. We combine analytical and numerical techniques to get insights into the stability properties. For a variant of the method that has been proposed in the literature we obtain precise analytic expressions for the asymptotic stability region. This allows us to prove a number of results. The technique introduced is widely applicable, and we use it to show that a fully implicit method suggested by [Kloeden and Platen, Numerical Solution of Stochastic Differential Equations, Springer-Verlag, 1992] has an asymptotic stability extension of the deterministic A-stability property. We also use the approach to explain some numerical results reported in [Milstein, Platen, and Schurz, SIAM J. Numer. Anal., 35 (1998), pp. 1010--1019.]