Adaptive timestepping strategies for nonlinear stochastic systems

We introduce a class of adaptive timestepping strategies for stochastic differential equations with non-Lipschitz drift coefficients. These strategies work by controlling potential unbounded growth in solutions of a numerical scheme due to the drift. We prove that the Euler-Maruyama scheme with an adaptive timestepping strategy in this class is strongly convergent. Specific strategies falling into this class are presented and demonstrated on a selection of numerical test problems. We observe that this approach is broadly applicable, can provide more dynamically accurate solutions than a drift-tamed scheme with fixed stepsize, and can improve MLMC simulations.

[1]  Andrew M. Stuart,et al.  Strong Convergence of Euler-Type Methods for Nonlinear Stochastic Differential Equations , 2002, SIAM J. Numer. Anal..

[2]  Vassilios Sotiropoulos,et al.  An adaptive time step scheme for a system of stochastic differential equations with multiple multiplicative noise: chemical Langevin equation, a proof of concept. , 2008, The Journal of chemical physics.

[3]  Michael B. Giles Multilevel Monte Carlo methods , 2015, Acta Numerica.

[4]  K. Burrage,et al.  Adaptive stepsize based on control theory for stochastic differential equations , 2004 .

[5]  S. Sabanis A note on tamed Euler approximations , 2013, 1303.5504.

[6]  Wolf-Jürgen Beyn,et al.  Stochastic C-Stability and B-Consistency of Explicit and Implicit Euler-Type Schemes , 2014, J. Sci. Comput..

[7]  M. Hutzenthaler,et al.  Numerical Approximations of Stochastic Differential Equations With Non-globally Lipschitz Continuous Coefficients , 2012, 1203.5809.

[8]  Jonathan C. Mattingly,et al.  An adaptive Euler–Maruyama scheme for SDEs: convergence and stability , 2006, math/0601029.

[9]  X. Mao,et al.  Stochastic Differential Equations and Applications , 1998 .

[10]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[11]  István Gyöngy,et al.  Convergence of tamed Euler schemes for a class of stochastic evolution equations , 2014, Stochastics and Partial Differential Equations: Analysis and Computations.

[12]  Zhongqiang Zhang,et al.  A Fundamental Mean-Square Convergence Theorem for SDEs with Locally Lipschitz Coefficients and Its Applications , 2012, SIAM J. Numer. Anal..

[13]  Michael B. Giles,et al.  Adaptive Euler-Maruyama Method for SDEs with Non-globally Lipschitz Drift: Part I, Finite Time Interval , 2016, 1703.06743.

[14]  Raúl Tempone,et al.  Implementation and analysis of an adaptive multilevel Monte Carlo algorithm , 2014, Monte Carlo Methods Appl..

[15]  P. Kloeden,et al.  Strong and weak divergence in finite time of Euler's method for stochastic differential equations with non-globally Lipschitz continuous coefficients , 2009, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[16]  Catherine E. Powell,et al.  An Introduction to Computational Stochastic PDEs , 2014 .

[17]  Xuerong Mao,et al.  On the local dynamics of polynomial difference equations with fading stochastic perturbations , 2010 .

[18]  Wayne H. Enright,et al.  Adaptive time-stepping for the strong numerical solution of stochastic differential equations , 2015, Numerical Algorithms.

[19]  Xuerong Mao,et al.  Convergence rates of the truncated Euler-Maruyama method for stochastic differential equations , 2016, J. Comput. Appl. Math..

[20]  A. R. Humphries,et al.  Dynamical Systems And Numerical Analysis , 1996 .

[21]  T. Shardlow,et al.  On the pathwise approximation of stochastic differential equations , 2014, 1409.2362.

[22]  Konstantinos Dareiotis,et al.  On Tamed Euler Approximations of SDEs Driven by Lévy Noise with Applications to Delay Equations , 2014, SIAM J. Numer. Anal..

[23]  P. Kloeden,et al.  Strong convergence of an explicit numerical method for SDEs with nonglobally Lipschitz continuous coefficients , 2010, 1010.3756.

[24]  J. M. Sanz-Serna,et al.  Optimal tuning of the hybrid Monte Carlo algorithm , 2010, 1001.4460.

[25]  Alexandra Rodkina,et al.  On the Use of Adaptive Meshes to Counter Overshoot in Solutions of Discretised Nonlinear Stochastic Differential Equations , 2010 .

[26]  Xiaojie Wang,et al.  The tamed Milstein method for commutative stochastic differential equations with non-globally Lipschitz continuous coefficients , 2011, 1102.0662.