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.

Enhanced Multilevel Monte Carlo Method Applied to FDTD for Probability Distribution Estimation

Probability Density Function Estimation in Multilevel Monte Carlo FDTD Method

Multilevel Monte Carlo FDTD Method for Uncertainty Quantification

Monte-Carlo Simulation of a Stochastic Differential Equation

Multilevel Monte Carlo methods for stochastic convection-diffusion eigenvalue problems

A Massively Parallel Implementation of Multilevel Monte Carlo for Finite Element Models

Multilevel Monte Carlo with Numerical Smoothing for Robust and Efficient Computation of Probabilities and Densities

Quantify uncertainty by estimating the probability density function of the output of interest using MLMC based Bayes method

Multilevel Monte Carlo methods for the Dean-Kawasaki equation from Fluctuating Hydrodynamics

A filtered multilevel Monte Carlo method for estimating the expectation of cell-centered discretized random fields

Efficient risk estimation via nested multilevel quasi-Monte Carlo simulation

Multilevel and Multi-index Monte Carlo methods for the McKean-Vlasov equation

An Efficient Stochastic HIE-FDTD Algorithm for Uncertainty Quantification of Electromagnetic Fields

A Stochastic FDTD Algorithm for Uncertainty Quantification of Electromagnetic-Thermal Simulation

Multilevel Markov Chain Monte Carlo Method for High-Contrast Single-Phase Flow Problems

Multilevel Monte Carlo Simulations of Composite Structures with Uncertain Manufacturing Defects

Multilevel Sequential Monte Carlo Samplers

Analysis and Application of Single Level, Multi-Level Monte Carlo and Quasi-Monte Carlo Finite Element Methods for Time-Dependent Maxwell's Equations with Random Inputs

Multilevel Asymptotic-Preserving Monte Carlo for Particle Simulations

Enhanced Feynman-Kac Formula for Nonlinear Uncertainty Quantification Problems

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