Abstract:In this study, we present an adaptive multilevel Monte Carlo (AMLMC) algorithm for approximating deterministic, real-valued, bounded linear functionals that depend on the solution of a linear elliptic PDE with a lognormal diffusivity coefficient and geometric singularities in bounded domains of Rd . Our AMLMC algorithm is built on the results of the weak convergence rates in the work (Moon et al., 2006) for an adaptive algorithm using isoparametric d-linear quadrilateral finite element approximations and the dual weighted residual error representation in a deterministic setting. Designed to suit the geometric nature of the singularities in the solution, our AMLMC algorithm uses a sequence of deterministic, non-uniform auxiliary meshes as a building block . The above-mentioned deterministic adaptive algorithm generates these meshes, corresponding to a geometrically decreasing sequence of tolerances. In particular, for a given realization of the diffusivity coefficient and accuracy level, AMLMC constructs its approximate sample using the first mesh in the hierarchy that satisfies the corresponding bias accuracy constraint. This adaptive approach is particularly useful for the lognormal case treated here, which lacks uniform coercivity and thus produces functional outputs that vary over orders of magnitude when sampled. Furthermore, we discuss iterative solvers and compare their efficiency with direct ones. To reduce computational work, we propose a stopping criterion for the iterative solver with respect to the quantity of interest, the realization of the diffusivity coefficient, and the desired level of AMLMC approximation. From the numerical experiments, based on a Fourier expansion of the diffusivity coefficient field, we observe improvements in efficiency compared with both standard Monte Carlo (MC) and standard MLMC (SMLMC) for a problem with a singularity similar to that at the tip of a slit modeling a crack.

Finite Element Error Analysis of Elliptic PDEs with Random Coefficients and Its Application to Multilevel Monte Carlo Methods

Engineering Analysis Error Estimation when Removing Finite-Sized Features in Nonlinear Elliptic Problems.

Application of Randomized Quadrature Formulas to the Finite Element Method for Elliptic Equations

Quasi-Monte Carlo sparse grid Galerkin finite element methods for linear elasticity equations with uncertainties

Corrector Analysis of a Heterogeneous Multi-scale Scheme for Elliptic Equations with Random Potential

Galerkin Finite Element Approximations of Stochastic Elliptic Partial Differential Equations

A convergent adaptive finite element stochastic Galerkin method based on multilevel expansions of random fields

Goal-oriented adaptive finite element multilevel Monte Carlo with convergence rates

Nonasymptotic Convergence Rate of Quasi-Monte Carlo: Applications to Linear Elliptic PDEs with Lognormal Coefficients and Importance Samplings

Reliable Error Estimates for Optimal Control of Linear Elliptic PDEs with Random Inputs

Quasi continuous level Monte Carlo for random elliptic PDEs

Multilevel Monte Carlo methods for stochastic convection-diffusion eigenvalue problems

Local decay rates of best-approximation errors using vector-valued finite elements for fields with low regularity and integrable curl or divergence

Projected Multilevel Monte Carlo Method for PDE with random input data

A row-sampling-based randomised finite element method for elliptic partial differential equations

Quasi-Monte Carlo and discontinuous Galerkin

A New Error Analysis for Finite Element Methods for Elliptic Neumann Boundary Control Problems with Pointwise Control Constraints

Monte Carlo convergence rates for kth moments in Banach spaces

Forward and Inverse Uncertainty Quantification using Multilevel Monte Carlo Algorithms for an Elliptic Nonlocal Equation

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

Circulant embedding with QMC -- analysis for elliptic PDE with lognormal coefficients