An Accelerated Convergence Scheme for Solving Stochastic Fractional Diffusion Equation

Abstract

An accelerated convergence scheme for temporal approximation of stochastic partial differential equation is presented. First, the regularity of the mild solution is provided. Combining the Itô formula and the remainder term of the exponential Euler scheme, this paper proposes a high accuracy time discretization method. Based on regularity results, a strong convergence rate for the discretization error

O({\tau }^{\frac{3}{2}-ϵ})

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$O(\tau ^{\frac{3}{2}-\epsilon })$$\end{document}

is proved for arbitrarily small

ϵ>0

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\epsilon >0$$\end{document}

. Here

\tau

\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\tau $$\end{document}

is the uniform time step size. Finally, the theoretical results are verified by several numerical experiments.

Keywords

Accelerated convergence scheme / Temporal approximation / Itô formula / Remainder term

References

[1] BréhierCE. Approximation of the invariant measure with an Euler scheme for stochastic PDEs driven by space-time white noise. Potential Anal., 2014, 40: 1-40.

[2] BréhierCE, CuiJ, HongJ. Strong convergence rates of semidiscrete splitting approximations for the stochastic Allen–Cahn equation. IMA J. Numer. Anal., 2019, 39: 2096-2134.

[3] BréhierCE, KopecM. Approximation of the invariant law of SPDEs: error analysis using a Poisson equation for a full-discretization scheme. IMA J. Numer. Anal., 2017, 37: 1375-1410.

[4] BréhierCE, VilmartG. High order integrator for sampling the invariant distribution of a class of parabolic stochastic PDEs with additive space-time noise. SIAM J. Sci. Comput., 2016, 38: 2283-2306.

[5] CuiJ, HongJ. Strong and weak convergence rates of a spatial approximation for stochastic partial differential equation with one-sided Lipschitz coefficient. SIAM J. Numer. Anal., 2019, 57: 1815-1841.

[6] CuiJ, HongJ, SunL. Weak convergence and invariant measure of a full discretization for parabolic SPDEs with non-globally Lipschitz coefficients. Stoch. Process. Appl., 2021, 134: 55-93.

[7] CusimanoN, del TesoF, Gerardo-GiordaL, PagniniG. Discretizations of the spectral fractional Laplacian on general domains with Dirichlet, Neumann, and Robin boundary conditions. SIAM J. Numer. Anal., 2018, 56: 1243-1272.

[8] GeissertM, KovácsM, LarssonS. Rate of weak convergence of the finite element method for the stochastic heat equation with additive noise. BIT, 2009, 49: 343-356.

[9] GrebenkovDS, NguyenBT. Geometrical structure of Laplacian eigenfunctions. SIAM Rev., 2013, 55: 601-667.

[10] JentzenA. Higher order pathwise numerical approximations of SPDEs with additive noise. SIAM J. Numer. Anal., 2011, 49: 642-667.

[11] KruseR, WuY. A randomized and fully discrete Galerkin finite element method for semilinear stochastic evolution equations. Math. Comp., 2019, 88: 2793-2825.

[12] LiuX. High-accuracy time discretization of stochastic fractional diffusion equation. J. Sci. Comput., 2022, 90: 1-24.

[13] LiuX. Strong approximation for fractional wave equation forced by fractional Brownian motion with Hurst parameter H∈(0,12)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$${H}\in (0,\frac{1}{2})$$\end{document}. J. Comput. Appl. Math., 2023, 432. 115285

[14] LiuX, DengWH. Numerical approximation for fractional diffusion equation forced by a tempered fractional Gaussian noise. J. Sci. Comput., 2020, 84: 1-28.

[15] LiuX, DengWH. Higher order approximation for stochastic space fractional wave equation forced by an additive space-time Gaussian noise. J. Sci. Comput., 2021, 87: 1-29.

[16] LiuZ, QiaoZ. Strong approximation of monotone stochastic partial differential equations driven by white noise. IAM J. Numer. Anal., 2020, 40: 1074-1093.

[17] NochettoRH, OtárolaE, SalgadoAJ. A PDE approach to fractional diffusion in general domains: a priori error analysis. Comput. Math., 2015, 15: 733-791.

[18] SongR, VondracekZ. Potential theory of subordinate killed Brownian motion in a domain. Probab. Theory Relat. Fields., 2003, 125: 578-592.

[19] TambueA, NgnotchouyeJMT. Weak convergence for a stochastic exponential integrator and finite element discretization of stochastic partial differential equation with multiplicative & additive noise. Appl. Numer. Math., 2016, 108: 57-86.

[20] WangX. Strong convergence rates of the linear implicit Euler method for the finite element discretization of SPDEs with additive noise. IAM J. Numer. Anal., 2017, 37: 965-984.

[21] WangX. An efficient explicit full-discrete scheme for strong approximation of stochastic Allen–Cahn equation. Stoch. Process. Appl., 2020, 130: 6271-6299.

[22] WangX, GanS. Weak convergence analysis of the linear implicit Euler method for semilinear stochastic partial differential equations with additive noise. J. Math. Anal. Appl., 2013, 398: 151-169.

RIGHTS & PERMISSIONS

Shanghai University