American Institute of Mathematical Sciences

Advanced
September  2018, 23(7): 2763-2773. doi: 10.3934/dcdsb.2018105

Numerical preservation of long-term dynamics by stochastic two-step methods

Raffaele D'Ambrosio 1, , Martina Moccaldi 2, and  Beatrice Paternoster 2,

1. 

Department of Engineering and Computer Science and Mathematics, University of L'Aquila, L'Aquila (AQ), Italy
2. 

Department of Mathematics, University of Salerno, Fisciano (SA), Italy

Received  March 2017 Revised  January 2018 Published  September 2018 Early access  March 2018

Fund Project: The work is supported by GNCS-Indam project.

The paper aims to explore the long-term behaviour of stochastic two-step methods applied to a class of second order stochastic differential equations. In particular, the treatment focuses on preserving long-term statistics related to the dynamics of a linear stochastic damped oscillator whose velocity, in the stationary regime, is distributed as a Gaussian variable and uncorrelated with the position. By computing the solution of a very simple matrix equality, we a-priori determine the long-term statistics characterizing the numerical dynamics and analyze the behaviour of a selection of methods.

Keywords: Stochastic differential equations, stochastic damped oscillators, stochastic two-step methods, long-term numerical dynamics.
Mathematics Subject Classification: 65C30.
Citation: Raffaele D'Ambrosio, Martina Moccaldi, Beatrice Paternoster. Numerical preservation of long-term dynamics by stochastic two-step methods. Discrete and Continuous Dynamical Systems - B, 2018, 23 (7) : 2763-2773. doi: 10.3934/dcdsb.2018105
References:
[1]

E. Buckwar and R. D'Ambrosio, Exponential mean-square stability properties of stochastic multistep methods, submitted.

[2]

E. Buckwar, R. Horvath-Bokor and R. Winkler, Asymptotic mean-square stability of two-step methods for stochastic ordinary differential equations, BIT Numer. Math., 46 (2006), 261-282.  doi: 10.1007/s10543-006-0060-5.

[3]

P. M. Burrage and K. Burrage, Structure-preserving Runge-Kutta methods for stochastic Hamiltonian equations with additive noise, Numer. Algor., 65 (2014), 519-532.  doi: 10.1007/s11075-013-9796-6.

[4]

P. M. Burrage and K. Burrage, Low rank Runge-Kutta methods, symplecticity and stochastic Hamiltonian problems with additive noise, J. Comput. Appl. Math., 236 (2012), 3920-3930.  doi: 10.1016/j.cam.2012.03.007.

[5]

K. Burrage, I. Lenane and G. Lythe, Numerical methods for second-order stochastic differential equations, SIAM J. Sci. Comput., 29 (2007), 245-264.  doi: 10.1137/050646032.

[6]

K. Burrage and G. Lythe, Accurate stationary densities with partitioned numerical methods for stochastic differential equations, SIAM J. Numer. Anal., 47 (2009), 1601-1618.  doi: 10.1137/060677148.

[7]

D. Conte, R. D'Ambrosio and B. Paternoster, On the stability of ϑ-methods for stochastic Volterra integral equations, Discr. Cont. Dyn. Sys. - B, accepted for publication, (2017).

[8]

P. Henrici, Discrete Variable Methods in Ordinary Differential Equations, Wiley, New York, 1962.

[9]

D. J. Higham, An algorithmic introduction to numerical simulation of stochastic differential equations, SIAM Rev., 43 (2001), 525-546.  doi: 10.1137/S0036144500378302.

[10]

P. E. Kloeden and E. Platen, The Numerical Solution of Stochastic Differential Equations, Springer-Verlag, 1992.

[11]

A. H. Strömmen and D. J. Melbö Higham, Numerical simulation of a linear stochastic oscillator with additive noise, Appl. Numer. Math., 51 (2004), 89-99.  doi: 10.1016/j.apnum.2004.02.003.

[12]

G. Vilmart, Weak second order multi-revolution composition methods for highly oscillatory stochastic differential equations with additive or multiplicative noise, SIAM J. Sci. Comput., 36 (2014), 1770-1796.  doi: 10.1137/130935331.

show all references

References:
[1]

E. Buckwar and R. D'Ambrosio, Exponential mean-square stability properties of stochastic multistep methods, submitted.

[2]

E. Buckwar, R. Horvath-Bokor and R. Winkler, Asymptotic mean-square stability of two-step methods for stochastic ordinary differential equations, BIT Numer. Math., 46 (2006), 261-282.  doi: 10.1007/s10543-006-0060-5.

[3]

P. M. Burrage and K. Burrage, Structure-preserving Runge-Kutta methods for stochastic Hamiltonian equations with additive noise, Numer. Algor., 65 (2014), 519-532.  doi: 10.1007/s11075-013-9796-6.

[4]

P. M. Burrage and K. Burrage, Low rank Runge-Kutta methods, symplecticity and stochastic Hamiltonian problems with additive noise, J. Comput. Appl. Math., 236 (2012), 3920-3930.  doi: 10.1016/j.cam.2012.03.007.

[5]

K. Burrage, I. Lenane and G. Lythe, Numerical methods for second-order stochastic differential equations, SIAM J. Sci. Comput., 29 (2007), 245-264.  doi: 10.1137/050646032.

[6]

K. Burrage and G. Lythe, Accurate stationary densities with partitioned numerical methods for stochastic differential equations, SIAM J. Numer. Anal., 47 (2009), 1601-1618.  doi: 10.1137/060677148.

[7]

D. Conte, R. D'Ambrosio and B. Paternoster, On the stability of ϑ-methods for stochastic Volterra integral equations, Discr. Cont. Dyn. Sys. - B, accepted for publication, (2017).

[8]

P. Henrici, Discrete Variable Methods in Ordinary Differential Equations, Wiley, New York, 1962.

[9]

D. J. Higham, An algorithmic introduction to numerical simulation of stochastic differential equations, SIAM Rev., 43 (2001), 525-546.  doi: 10.1137/S0036144500378302.

[10]

P. E. Kloeden and E. Platen, The Numerical Solution of Stochastic Differential Equations, Springer-Verlag, 1992.

[11]

A. H. Strömmen and D. J. Melbö Higham, Numerical simulation of a linear stochastic oscillator with additive noise, Appl. Numer. Math., 51 (2004), 89-99.  doi: 10.1016/j.apnum.2004.02.003.

[12]

G. Vilmart, Weak second order multi-revolution composition methods for highly oscillatory stochastic differential equations with additive or multiplicative noise, SIAM J. Sci. Comput., 36 (2014), 1770-1796.  doi: 10.1137/130935331.

Figure 1.  Patterns over $\eta$ of $|\widetilde{\sigma}_{x,EM}^2-\sigma_x^2|$ (continuous line), $|\widetilde{\sigma}_{v,EM}^2-\sigma_v^2|$ (dashed line) and $|\widetilde{\mu}_{EM}-\mu|$ (dashed-dotted line), for $g = 1$, $\Delta t = 10^{-2}$, $\varepsilon = 1$ for the Eulero-Maruyama method (16) applied to the stochastic problem (4).
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 2.  Patterns over $\eta$ of $|\widetilde{\sigma}_{x,TRAP}^2-\sigma_x^2|$ (continuous line), $|\widetilde{\sigma}_{v,TRAP}^2-\sigma_v^2|$ (dashed line, almost overlapping the continuous line) and $|\widetilde{\mu}_{TRAP}-\mu|$ (dashed-dotted line), for $g = 1$, $\Delta t = 10^{-2}$, $\varepsilon = 1$ for the trapezoidal method (17) applied to the stochastic problem (4).
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 3.  Patterns over $\eta$ of $|\widetilde{\sigma}_{x,AM}^2-\sigma_x^2|$ (continuous line), $|\widetilde{\sigma}_{v,AM}^2-\sigma_v^2|$ (dashed line) and $|\widetilde{\mu}_{AM}-\mu|$ (dashed-dotted line), for $g = 1$, $\Delta t = 10^{-2}$, $\varepsilon = 1$ for the Adams-Moulton method (19) applied to the stocastic problem (4).
Figure Options

Download full-size image

Download as PowerPoint slide

Figure 4.  Patterns over $\eta$ of $|\widetilde{\sigma}_{x,BDF}^2-\sigma_x^2|$ (continuous line), $|\widetilde{\sigma}_{v,BDF}^2-\sigma_v^2|$ (dashed line) and $|\widetilde{\mu}_{BDF}-\mu|$ (dashed-dotted line), for $g = 1$, $\Delta t = 10^{-2}$, $\varepsilon = 1$ for the BDF method (21) applied to the stochastic problem (4).
Figure Options

Download full-size image

Download as PowerPoint slide

[1]

Angelamaria Cardone, Dajana Conte, Beatrice Paternoster. Two-step collocation methods for fractional differential equations. Discrete and Continuous Dynamical Systems - B, 2018, 23 (7) : 2709-2725. doi: 10.3934/dcdsb.2018088
[2]

Quan Zhou, Yabing Sun. High order one-step methods for backward stochastic differential equations via Itô-Taylor expansion. Discrete and Continuous Dynamical Systems - B, 2022, 27 (8) : 4387-4413. doi: 10.3934/dcdsb.2021233
[3]

William F. Thompson, Rachel Kuske, Yue-Xian Li. Stochastic phase dynamics of noise driven synchronization of two conditional coherent oscillators. Discrete and Continuous Dynamical Systems, 2012, 32 (8) : 2971-2995. doi: 10.3934/dcds.2012.32.2971
[4]

Igor Chueshov, Peter E. Kloeden, Meihua Yang. Long term dynamics of second order-in-time stochastic evolution equations with state-dependent delay. Discrete and Continuous Dynamical Systems - B, 2018, 23 (3) : 991-1009. doi: 10.3934/dcdsb.2018139
[5]

Tingting Wu, Yufei Yang, Huichao Jing. Two-step methods for image zooming using duality strategies. Numerical Algebra, Control and Optimization, 2014, 4 (3) : 209-225. doi: 10.3934/naco.2014.4.209
[6]

Jiaohui Xu, Tomás Caraballo. Long time behavior of fractional impulsive stochastic differential equations with infinite delay. Discrete and Continuous Dynamical Systems - B, 2019, 24 (6) : 2719-2743. doi: 10.3934/dcdsb.2018272
[7]

Nguyen Dinh Cong, Nguyen Thi Thuy Quynh. Coincidence of Lyapunov exponents and central exponents of linear Ito stochastic differential equations with nondegenerate stochastic term. Conference Publications, 2011, 2011 (Special) : 332-342. doi: 10.3934/proc.2011.2011.332
[8]

Nikolaos Halidias. Construction of positivity preserving numerical schemes for some multidimensional stochastic differential equations. Discrete and Continuous Dynamical Systems - B, 2015, 20 (1) : 153-160. doi: 10.3934/dcdsb.2015.20.153
[9]

Minoo Kamrani. Numerical solution of partial differential equations with stochastic Neumann boundary conditions. Discrete and Continuous Dynamical Systems - B, 2019, 24 (10) : 5337-5354. doi: 10.3934/dcdsb.2019061
[10]

Yulan Lu, Minghui Song, Mingzhu Liu. Convergence rate and stability of the split-step theta method for stochastic differential equations with piecewise continuous arguments. Discrete and Continuous Dynamical Systems - B, 2019, 24 (2) : 695-717. doi: 10.3934/dcdsb.2018203
[11]

Ying Liu, Yabing Sun, Weidong Zhao. Explicit multistep stochastic characteristic approximation methods for forward backward stochastic differential equations. Discrete and Continuous Dynamical Systems - S, 2022, 15 (4) : 773-795. doi: 10.3934/dcdss.2021044
[12]

Dingshi Li, Lin Shi, Xiaohu Wang. Long term behavior of stochastic discrete complex Ginzburg-Landau equations with time delays in weighted spaces. Discrete and Continuous Dynamical Systems - B, 2019, 24 (9) : 5121-5148. doi: 10.3934/dcdsb.2019046
[13]

Philippe Angot, Pierre Fabrie. Convergence results for the vector penalty-projection and two-step artificial compressibility methods. Discrete and Continuous Dynamical Systems - B, 2012, 17 (5) : 1383-1405. doi: 10.3934/dcdsb.2012.17.1383
[14]

Ludwig Arnold, Igor Chueshov. Cooperative random and stochastic differential equations. Discrete and Continuous Dynamical Systems, 2001, 7 (1) : 1-33. doi: 10.3934/dcds.2001.7.1
[15]

Raffaele D'Ambrosio, Nicola Guglielmi, Carmela Scalone. Destabilising nonnormal stochastic differential equations. Discrete and Continuous Dynamical Systems - B, 2022  doi: 10.3934/dcdsb.2022140
[16]

Jianhua Huang, Tianlong Shen, Yuhong Li. Dynamics of stochastic fractional Boussinesq equations. Discrete and Continuous Dynamical Systems - B, 2015, 20 (7) : 2051-2067. doi: 10.3934/dcdsb.2015.20.2051
[17]

Igor Chueshov, Irena Lasiecka, Daniel Toundykov. Long-term dynamics of semilinear wave equation with nonlinear localized interior damping and a source term of critical exponent. Discrete and Continuous Dynamical Systems, 2008, 20 (3) : 459-509. doi: 10.3934/dcds.2008.20.459
[18]

Lijin Wang, Pengjun Wang, Yanzhao Cao. Numerical methods preserving multiple Hamiltonians for stochastic Poisson systems. Discrete and Continuous Dynamical Systems - S, 2022, 15 (4) : 819-836. doi: 10.3934/dcdss.2021095
[19]

Andrew Vlasic. Long-run analysis of the stochastic replicator dynamics in the presence of random jumps. Journal of Dynamics and Games, 2018, 5 (4) : 283-309. doi: 10.3934/jdg.2018018
[20]

Christopher Rackauckas, Qing Nie. Adaptive methods for stochastic differential equations via natural embeddings and rejection sampling with memory. Discrete and Continuous Dynamical Systems - B, 2017, 22 (7) : 2731-2761. doi: 10.3934/dcdsb.2017133

2021 Impact Factor: 1.497

Tools

Metrics

  • PDF downloads (257)
  • HTML views (377)
  • Cited by (4)

Other articles
by authors

[Back to Top]

Copyright © 2022 American Institute of Mathematical Sciences | Privacy Policy