Numerical solution of fractional elliptic stochastic PDEs with spatial white noise

David Bolin, Kristin Kirchner*, Mihály Kovács

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

35 Scopus citations


The numerical approximation of solutions to stochastic partial differential equations with additive spatial white noise on bounded domains in Rd is considered. The differential operator is given by the fractional power Lβ, β ∈ (0,1) of an integer-order elliptic differential operator L and is therefore nonlocal. Its inverse L-β is represented by a Bochner integral from the Dunford-Taylor functional calculus. By applying a quadrature formula to this integral representation the inverse fractional-order operator L-β is approximated by a weighted sum of nonfractional resolvents (I + exp(2 yℓ) L)-1 at certain quadrature nodes tj > 0. The resolvents are then discretized in space by a standard finite element method. This approach is combined with an approximation of the white noise, which is based only on the mass matrix of the finite element discretization. In this way an efficient numerical algorithm for computing samples of the approximate solution is obtained. For the resulting approximation the strong mean-square error is analyzed and an explicit rate of convergence is derived. Numerical experiments for L=κ2-Δ, κ> 0 with homogeneous Dirichlet boundary conditions on the unit cube (0,1)d in d=1,2,3 spatial dimensions for varying β ∈ (0,1) attest to the theoretical results.

Original languageEnglish (US)
Pages (from-to)1051-1073
Number of pages23
JournalIMA Journal of Numerical Analysis
Issue number2
StatePublished - Apr 24 2020


  • Gaussian white noise
  • Matérn covariances
  • finite element methods
  • fractional operators
  • spatial statistics
  • stochastic partial differential equations

ASJC Scopus subject areas

  • General Mathematics
  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'Numerical solution of fractional elliptic stochastic PDEs with spatial white noise'. Together they form a unique fingerprint.

Cite this