A full-discrete exponential Euler approximation of the invariant measure for parabolic stochastic partial differential equations

Abstract We discrete the ergodic semilinear stochastic partial differential equations in space dimension d ≤ 3 with additive noise, spatially by a spectral Galerkin method and temporally by an exponential Euler scheme. It is shown that both the spatial semi-discretization and the spatio-temporal full discretization are ergodic. Further, convergence orders of the numerical invariant measures, depending on the regularity of noise, are recovered based on an easy time-independent weak error analysis without relying on Malliavin calculus. To be precise, the convergence order is 1 − ϵ in space and 1 2 − ϵ in time for the space-time white noise case and 2 − ϵ in space and 1 − ϵ in time for the trace class noise case in space dimension d = 1 , with arbitrarily small ϵ > 0 . Numerical results are finally reported to confirm these theoretical findings.

[1]  Jialin Hong,et al.  Approximation of Invariant Measures for Stochastic Differential Equations with Piecewise Continuous Arguments via Backward Euler Method , 2019, ArXiv.

[2]  P. Kloeden,et al.  Overcoming the order barrier in the numerical approximation of stochastic partial differential equations with additive space–time noise , 2009, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[3]  M. Röckner,et al.  A Concise Course on Stochastic Partial Differential Equations , 2007 .

[4]  Xiaojie Wang,et al.  Full discretisation of semi-linear stochastic wave equations driven by multiplicative noise , 2015, 1503.00073.

[5]  Andrew M. Stuart,et al.  Convergence of Numerical Time-Averaging and Stationary Measures via Poisson Equations , 2009, SIAM J. Numer. Anal..

[6]  A. Borovkov.  Ergodicity and stability of stochastic processes , 1998 .

[7]  D. Talay.  Second-order discretization schemes of stochastic differential systems for the computation of the invariant law , 1990 .

[8]  J. T. Ngnotchouye,et al.  Weak convergence for a stochastic exponential integrator and finite element discretization of stochastic partial differential equation with multiplicative & additive noise , 2016 .

[9]  Stig Larsson,et al.  Weak convergence for a spatial approximation of the nonlinear stochastic heat equation , 2012, Math. Comput..

[10]  D. Talay.  Stochastic Hamiltonian Systems : Exponential Convergence to the Invariant Measure , and Discretization by the Implicit Euler Scheme , 2002 .

[11]  Xiaojie Wang,et al.  An Exponential Integrator Scheme for Time Discretization of Nonlinear Stochastic Wave Equation , 2013, J. Sci. Comput..

[12]  Stig Larsson,et al.  Finite Element Approximation of the Linear Stochastic Wave Equation with Additive Noise , 2010, SIAM J. Numer. Anal..

[13]  S. Cerrai.  Second Order Pde's in Finite and Infinite Dimension: A Probabilistic Approach , 2001 .

[14]  Stig Larsson,et al.  A Trigonometric Method for the Linear Stochastic Wave Equation , 2012, SIAM J. Numer. Anal..

[15]  R. Kruse,et al.  Duality in refined Sobolev–Malliavin spaces and weak approximation of SPDE , 2013, 1312.5893.

[16]  Arnulf Jentzen,et al.  Weak Convergence Rates for Euler-Type Approximations of Semilinear Stochastic Evolution Equations with Nonlinear Diffusion Coefficients , 2015, Found. Comput. Math..

[17]  Xicheng Zhang,et al.  Strong Feller properties for degenerate SDEs with jumps , 2013, 1312.7380.

[18]  Arnulf Jentzen,et al.  Weak convergence rates of spectral Galerkin approximations for SPDEs with nonlinear diffusion coefficients , 2014, The Annals of Applied Probability.

[19]  Xiaojie Wang,et al.  Weak convergence analysis of the linear implicit Euler method for semilinear stochastic partial differential equations with additive noise , 2013 .

[20]  A. Stuart,et al.  Weak Error Estimates for Trajectories of SPDEs Under Spectral Galerkin Discretization , 2016, Journal of Computational Mathematics.

[21]  Assyr Abdulle,et al.  High Order Numerical Approximation of the Invariant Measure of Ergodic SDEs , 2014, SIAM J. Numer. Anal..

[22]  Charles-Edouard Bréhier,et al.  Approximation of the Invariant Measure with an Euler Scheme for Stochastic PDEs Driven by Space-Time White Noise , 2012, 1202.2707.

[23]  Peter E. Kloeden,et al.  The exponential integrator scheme for stochastic partial differential equations: Pathwise error bounds , 2011, J. Comput. Appl. Math..

[24]  G. Prato.  An Introduction to Infinite-Dimensional Analysis , 2006 .

[25]  P. Kloeden,et al.  The Numerical Approximation of Stochastic Partial Differential Equations , 2009 .

[26]  Jialin Hong,et al.  Approximation of Invariant Measure for Damped Stochastic Nonlinear Schrödinger Equation via an Ergodic Numerical Scheme , 2015, 1509.09148.

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

[28]  Di Liu,et al.  Ergodic Approximation to Chemical Reaction System with Delay , 2019, SIAM J. Numer. Anal..

[29]  G. Lord,et al.  Stochastic exponential integrators for finite element discretization of SPDEs for multiplicative and additive noise , 2011, 1103.1986.

[30]  David Cohen,et al.  A fully discrete approximation of the one-dimensional stochastic heat equation , 2015, IMA Journal of Numerical Analysis.

[31]  G. Sell,et al.  Dynamics of Evolutionary Equations , 2002 .

[32]  J. Zabczyk,et al.  Ergodicity for Infinite Dimensional Systems: Invariant measures for stochastic evolution equations , 1996 .

[33]  Jonathan C. Mattingly,et al.  Ergodicity for SDEs and approximations: locally Lipschitz vector fields and degenerate noise , 2002 .

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

[35]  J. Zabczyk,et al.  Stochastic Equations in Infinite Dimensions , 2008 .

[36]  Charles-Edouard Br'ehier,et al.  Kolmogorov equations and weak order analysis for SPDEs with nonlinear diffusion coefficient , 2017, Journal de Mathématiques Pures et Appliquées.

[37]  R. Kruse.  Strong and Weak Approximation of Semilinear Stochastic Evolution Equations , 2013 .

[38]  Jerzy Zabczyk,et al.  Ergodicity for Infinite Dimensional Systems: Stochastic delay systems , 1996 .

[39]  Arnaud Debussche,et al.  Weak approximation of stochastic partial differential equations: the nonlinear case , 2008, Math. Comput..

[40]  Xiaojie Wang,et al.  Weak error estimates of the exponential Euler scheme for semi-linear SPDEs without Malliavin calculus , 2014, 1408.0713.

[41]  M. V. Tretyakov,et al.  Computing ergodic limits for Langevin equations , 2007 .

[42]  Charles-Edouard Br'ehier,et al.  Approximation of the invariant law of SPDEs: error analysis using a Poisson equation for a full-discretization scheme , 2013, 1311.7030.

[43]  Arnulf Jentzen,et al.  Efficient simulation of nonlinear parabolic SPDEs with additive noise , 2011, 1210.8320.

[44]  Liying Zhang,et al.  Numerical Analysis on Ergodic Limit of Approximations for Stochastic NLS Equation via Multi-symplectic Scheme , 2017, SIAM J. Numer. Anal..

[45]  Charles-Edouard Bréhier,et al.  High Order Integrator for Sampling the Invariant Distribution of a Class of Parabolic Stochastic PDEs with Additive Space-Time Noise , 2015, SIAM J. Sci. Comput..

[46]  Arnulf Jentzen,et al.  Weak convergence rates for temporal numerical approximations of stochastic wave equations with multiplicative noise , 2019, 1901.05535.