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.

An MLMCE-HDG Method for the Convection Diffusion Equation with Random Diffusivity

An EMC-HDG Scheme for the Convection-Diffusion Equation with Random Diffusivity

A Multilevel Monte Carlo Ensemble and Hybridizable Discontinuous Galerkin Method for a Stochastic Parabolic Problem

Monte-Carlo Simulation of a Stochastic Differential Equation

Convergence Analysis and Cost Estimate of an MLMC-HDG Method for Elliptic PDEs with Random Coefficients

Multilevel Monte Carlo methods for stochastic convection-diffusion eigenvalue problems

Monte Carlo Finite Volume Element Methods for the Convection-Diffusion Equation with A Random Diffusion Coefficient

An Ensemble Monte Carlo HDG Method for Parabolic PDEs with Random Coefficients

Multi-Monte Carlo Method for Coagulation and Condensation/evaporation in Dispersed Systems.

An efficient multi-modes Monte Carlo homogenization method for random materials

Ensemble Time-Stepping Algorithm for the Convection-Diffusion Equation with Random Diffusivity

Stability Analysis and Error Estimates of Local Discontinuous Galerkin Methods for Convection–diffusion Equations on Overlapping Meshes

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

A Deep Learning Based Discontinuous Galerkin Method for Hyperbolic Equations with Discontinuous Solutions and Random Uncertainties

A discrete Hermite moments based multiple-relaxation-time lattice Boltzmann model for convection-diffusion equations

Multi-level Monte Carlo Weak Galerkin Method with Nested Meshes for Stochastic Brinkman Problem

Antithetic Multilevel Methods for Elliptic and Hypo-Elliptic Diffusions with Applications

Numerical Solutions of Convection-Diffusion Equations by Hybrid Discontinuous Galerkin Methods

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

An efficient Monte Carlo interior penalty discontinuous Galerkin method for the time-harmonic Maxwell's equations with random coefficients

An Adaptive Staggered Discontinuous Galerkin Method for the Steady State Convection–Diffusion Equation