Author
CROCI, M
GILES, M
Rognes, M
Farrell, P
Journal title
SIAM/ASA Journal on Uncertainty Quantification
DOI
10.1137/18M1175239
Issue
4
Volume
6
Last updated
2024-03-26T20:03:35.93+00:00
Page
1630-1655
Abstract
When solving stochastic partial differential equations (SPDEs) driven by
additive spatial white noise, the efficient sampling of white noise
realizations can be challenging. Here, we present a new sampling technique that can be used to efficiently compute white noise samples in a finite element method and multilevel Monte Carlo (MLMC) setting. The key idea is to exploit the finite element matrix assembly procedure and factorize each local mass matrix independently, hence avoiding the factorization of a large matrix. Moreover, in a MLMC framework, the white noise samples must be coupled between subsequent levels. We show how our technique can be used to enforce this coupling even in the case of non-nested mesh hierarchies. We demonstrate the efficacy of our method with numerical experiments. We observe optimal convergence rates for the finite element solution of the elliptic SPDEs of interest in 2D and 3D and we show convergence of the sampled field covariances. In a MLMC setting, a good coupling is enforced and the telescoping sum is respected.
Symplectic ID
865752
Favourite
On
Publication type
Journal Article
Publication date
20 Nov 2018
Please contact us with feedback and comments about this page. Created on 10 Jul 2018 - 14:40.