Fully discrete stochastic maximal regularity and $H^{\infty }$-calculus for second-order elliptic operators

Abstract

This paper establishes the fully discrete stochastic maximal $L^p$-regularity and the accompanying sharp maximal estimate for numerical approximations of parabolic stochastic partial differential equations. We consider the spatial finite element discretization $A_h$ of a general second-order elliptic operator $A=-\nabla \cdot a\nabla +b\cdot \nabla +c$ with Dirichlet boundary conditions on a smooth, bounded, convex domain in ${\mathbb{R}}^3$, coupled with a broad class of temporal schemes, including rational approximations and the exponential Euler method. To obtain these optimal discrete regularity results, we establish a bounded $H^{\infty }$-calculus for the discrete spatial operator $A_h$, uniformly in the mesh size $h$. As a direct byproduct, we also establish the discrete-in-space stochastic maximal regularity for the corresponding spatial semi-discretizations.