Higher order numerical schemes for SPDEs with additive Noise

Abstract

We present high-order numerical schemes for linear stochastic heat and wave equations with Dirichlet boundary conditions, driven by additive noise. Standard Euler schemes for SPDEs are limited to an order convergence between 1/2 and 1 due to the low temporal regularity of noise. For the stochastic heat equation, a modified Crank-Nicolson scheme with proper numerical quadrature rule for the noise term in its reformulation as random PDE achieves a strong convergence rate of 3/2. For the stochastic wave equation with additive noise a corresponding approach leads to a scheme which is of order 2.

14 pages, 2 figures