Multilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise sampling

11/27/2019
by   M. Croci, et al.
0

When solving partial differential equations with random fields as coefficients the efficient sampling of random field realisations can be challenging. In this paper we focus on the fast sampling of Gaussian fields using quasi-random points in a finite element and multilevel quasi Monte Carlo (MLQMC) setting. Our method uses the SPDE approach combined with a new fast (ML)QMC algorithm for white noise sampling. We express white noise as a wavelet series expansion that we divide in two parts. The first part is sampled using quasi-random points and contains a finite number of terms in order of decaying importance to ensure good QMC convergence. The second part is a correction term which is sampled using standard pseudo-random numbers. We show how the sampling of both terms can be performed in linear time and memory complexity in the number of mesh cells via a supermesh construction, yielding an overall linear cost. Furthermore, our technique can be used to enforce the MLQMC coupling even in the case of non-nested mesh hierarchies. We demonstrate the efficacy of our method with numerical experiments.

READ FULL TEXT
research
06/23/2023

Smoothed Circulant Embedding with Applications to Multilevel Monte Carlo Methods for PDEs with Random Coefficients

We consider the computational efficiency of Monte Carlo (MC) and Multile...
research
12/15/2020

On the Selection of Random Field Evaluation Points in the p-MLQMC Method

Engineering problems are often characterized by significant uncertainty ...
research
07/28/2020

Multilevel Hierarchical Decomposition of Finite Element White Noise with Application to Multilevel Markov Chain Monte Carlo

In this work we develop a new hierarchical multilevel approach to genera...
research
08/01/2019

On the Numerical Approximation of the Karhunen-Loève Expansion for Lognormal Random Fields

The Karhunen-Loève (KL) expansion is a popular method for approximating ...
research
07/15/2022

Quasi-Monte Carlo and discontinuous Galerkin

In this study, we design and develop Quasi-Monte Carlo (QMC) cubatures f...
research
02/05/2022

Analyticity of Parametric and Stochastic Elliptic Eigenvalue Problems with an Application to Quasi-Monte Carlo Methods

In the present paper, we study the analyticity of the leftmost eigenvalu...
research
10/27/2022

Quasi-Monte Carlo finite element approximation of the Navier-Stokes equations with initial data modeled by log-normal random fields

In this paper, we analyze the numerical approximation of the Navier-Stok...

Please sign up or login with your details

Forgot password? Click here to reset