Multilevel quasi-Monte Carlo for random elliptic eigenvalue problems I: Regularity and error analysis

10/02/2020
by   Alexander D. Gilbert, et al.
0

Random eigenvalue problems are useful models for quantifying the uncertainty in several applications from the physical sciences and engineering, e.g., structural vibration analysis, the criticality of a nuclear reactor or photonic crystal structures. In this paper we present a simple multilevel quasi-Monte Carlo (MLQMC) method for approximating the expectation of the minimal eigenvalue of an elliptic eigenvalue problem with coefficients that are given as a series expansion of countably-many stochastic parameters. The MLQMC algorithm is based on a hierarchy of discretisations of the spatial domain and truncations of the dimension of the stochastic parameter domain. To approximate the expectations, randomly shifted lattice rules are employed. This paper is primarily dedicated to giving a rigorous analysis of the error of this algorithm. A key step in the error analysis requires bounds on the mixed derivatives of the eigenfunction with respect to both the stochastic and spatial variables simultaneously. An accompanying paper [Gilbert and Scheichl, 2020], focusses on practical extensions of the MLQMC algorithm to improve efficiency, and presents numerical results.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
03/05/2021

Multilevel quasi-Monte Carlo for random elliptic eigenvalue problems II: Efficient algorithms and numerical results

Stochastic PDE eigenvalue problems often arise in the field of uncertain...
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
03/07/2023

Multilevel Monte Carlo methods for stochastic convection-diffusion eigenvalue problems

We develop new multilevel Monte Carlo (MLMC) methods to estimate the exp...
research
08/06/2023

Semilinear elliptic eigenvalue problem: Parametric analyticity and the uncertainty quantification

In this paper, to the best of our knowledge, we make the first attempt a...
research
02/01/2023

Multilevel Monte Carlo FEM for Elliptic PDEs with Besov Random Tree Priors

We develop a multilevel Monte Carlo (MLMC)-FEM algorithm for linear, ell...
research
02/01/2023

Multilevel Markov Chain Monte Carlo for Bayesian Elliptic Inverse Problems with Besov Random Tree Priors

We propose a multilevel Monte Carlo-FEM algorithm to solve elliptic Baye...
research
04/26/2020

Quasi-Monte Carlo finite element analysis for wave propagation in heterogeneous random media

We propose and analyze a quasi-Monte Carlo (QMC) algorithm for efficient...

Please sign up or login with your details

Forgot password? Click here to reset