Numerical solution of fractional elliptic stochastic PDEs with spatial white noise

Abstract

The numerical approximation of solutions to stochastic partial differential equations with additive spatial white noise on bounded domains in ${\mathbb{R}}^d$ is considered. The differential operator is given by the fractional power $L^{\beta }$, $\beta \in (0,1)$, of an integer order elliptic differential operator $L$ and is therefore non-local. Its inverse $L^{-\beta }$ is represented by a Bochner integral from the Dunford-Taylor functional calculus. By applying a quadrature formula to this integral representation, the inverse fractional power operator $L^{-\beta }$ is approximated by a weighted sum of non-fractional resolvents $(I+t_j^{2}L{)}^{-1}$ at certain quadrature nodes $t_j>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={\kappa }^2-\Delta$, $\kappa >0$, with homogeneous Dirichlet boundary conditions on the unit cube $(0,1{)}^d$ in $d=1,2,3$ spatial dimensions for varying $\beta \in (0,1)$ attest the theoretical results.

21 pages, 1 figure

References in corpus (1)

• Spatial models generated by nested stochastic partial differential equations, with an application to global ozone mapping

Cited by in corpus (27)

• The rational SPDE approach for Gaussian random fields with general smoothness
• Efficient white noise sampling and coupling for multilevel Monte Carlo with non-nested meshes
• Matérn Gaussian processes on Riemannian manifolds
• Fast sampling of parameterised Gaussian random fields
• Weak convergence of Galerkin approximations for fractional elliptic stochastic PDEs with spatial white noise
• Deformed SPDE models with an application to spatial modeling of significant wave height
• Regularity and convergence analysis in Sobolev and Hölder spaces for generalized Whittle-Matérn fields
• A fractional PDE model for turbulent velocity fields near solid walls
• Covariance-based rational approximations of fractional SPDEs for computationally efficient Bayesian inference
• Equivalence of measures and asymptotically optimal linear prediction for Gaussian random fields with fractional-order covariance operators
• Graph-based Prior and Forward Models for Inverse Problems on Manifolds with Boundaries
• Surface finite element approximation of spherical Whittle--Matérn Gaussian random fields
• Multilevel approximation of Gaussian random fields: Covariance compression, estimation and spatial prediction
• Galerkin--Chebyshev approximation of Gaussian random fields on compact Riemannian manifolds
• Hybrid simulation scheme for volatility modulated moving average fields
• Hilbert--Schmidt regularity of symmetric integral operators on bounded domains with applications to SPDE approximations
• The SPDE approach for Gaussian and non-Gaussian fields: 10 years and still running
• An a posteriori error estimator for the spectral fractional power of the Laplacian
• Regularity theory for a new class of fractional parabolic stochastic evolution equations
• Finite Element Representations of Gaussian Processes: Balancing Numerical and Statistical Accuracy
• A statistical framework for generating microstructures of two-phase random materials: application to fatigue analysis
• Multiple and weak Markov properties in Hilbert spaces with applications to fractional stochastic evolution equations
• Multilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise sampling
• Numerical analysis for stochastic time-space fractional diffusion equation driven by fractional Gaussion noise
• Approximation of SPDE covariance operators by finite elements: A semigroup approach
• Spectral approximation of a new class of stochastic fractional evolution equations
• Weak convergence of fully discrete finite element approximations of semilinear hyperbolic SPDE with additive noise