Efficient white noise sampling and coupling for multilevel Monte Carlo with non-nested meshes

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.

28 pages, 10 figures

References in corpus (2)

• Spatial models generated by nested stochastic partial differential equations, with an application to global ozone mapping
• Numerical solution of fractional elliptic stochastic PDEs with spatial white noise

Cited by in corpus (12)

• Optimal sensor and actuator placement for feedback control of vortex shedding
• A fractional PDE model for turbulent velocity fields near solid walls
• A probabilistic finite element method based on random meshes: Error estimators and Bayesian inverse problems
• Multi-output multilevel best linear unbiased estimators via semidefinite programming
• Fast uncertainty quantification of tracer distribution in the brain interstitial fluid with multilevel and quasi Monte Carlo
• hIPPYlib: An Extensible Software Framework for Large-Scale Inverse Problems Governed by PDEs; Part I: Deterministic Inversion and Linearized Bayesian Inference
• Resolvent-based approach for H2-optimal estimation and control: an application to the cylinder flow
• Goal-Oriented Adaptive Finite Element Multilevel Monte Carlo with Convergence Rates
• Multifidelity Approximate Bayesian Computation with Sequential Monte Carlo Parameter Sampling
• Complexity bounds on supermesh construction for quasi-uniform meshes
• A statistical framework for generating microstructures of two-phase random materials: application to fatigue analysis
• Multilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise sampling