A spectral Galerkin exponential Euler time-stepping scheme for parabolic SPDEs on two-dimensional domains with a C2-boundary

TL;DR

This paper develops a spectral Galerkin exponential Euler scheme for numerically solving semi-linear parabolic SPDEs on complex 2D domains with C2 boundaries, combining boundary integral methods and eigenfunction expansions.

Contribution

It introduces a novel spectral Galerkin and boundary element approach for SPDEs on irregular 2D domains, including error analysis and numerical validation.

Findings

Effective eigenvalue and eigenfunction computation for irregular shapes

Error bounds for the proposed numerical scheme

Numerical results demonstrating the method on asymmetric domains

Abstract

We consider the numerical approximation of second-order semi-linear parabolic stochastic partial differential equations interpreted in the mild sense which we solve on general two-dimensional domains with a $\mathcal{C}^2$ boundary with homogeneous Dirichlet boundary conditions. The equations are driven by Gaussian additive noise, and several Lipschitz-like conditions are imposed on the nonlinear function. We discretize in space with a spectral Galerkin method and in time using an explicit Euler-like scheme. For irregular shapes, the necessary Dirichlet eigenvalues and eigenfunctions are obtained from a boundary integral equation method. This yields a nonlinear eigenvalue problem, which is discretized using a boundary element collocation method and is solved with the Beyn contour integral algorithm. We present an error analysis as well as numerical results on an exemplary asymmetric shape, and point out limitations of the approach.