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.

Convergence of adaptive stochastic collocation with finite elements

Convergence analysis of the adaptive stochastic collocation finite element method

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

An adaptive hierarchical sparse grid collocation algorithm for the solution of stochastic differential equations

On the convergence of adaptive Galerkin FEM for parametric PDEs with lognormal coefficients

Finite Element Methods for the Stochastic Allen-Cahn Equation with Gradient-type Multiplicative Noises

Strong convergence of an adaptive time-stepping Milstein method for SDEs with monotone coefficients

Convergence of Adaptive Finite Element Approximations for Nonlinear Eigenvalue Problems

A finite element approximation for the stochastic Landau--Lifshitz--Gilbert equation with multi-dimensional noise

Strong and weak convergence rates of finite element method for stochastic partial differential equation with non-sided Lipschitz coefficient

Strong convergence of a class of adaptive numerical methods for SDEs with jumps

Convergence and Quasi-Optimality of Adaptive FEM with Inhomogeneous Dirichlet Data

On the convergence of adaptive approximations for stochastic differential equations

A spectral collocation method for nonlocal diffusion equations

Superconvergence and Postprocessing of Collocation Methods for Fractional Differential Equations

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

A two-level stochastic collocation method for semilinear elliptic equations with random coefficients

Convergence of Adaptive Mixed Finite Element Methods for the Hodge Laplacian Equation: without harmonic forms

Convergence of adaptive finite element methods for eigenvalue problems

On a finite-volume approximation of a diffusion-convection equation with a multiplicative stochastic force

Approximation of Stochastic Partial Differential Equations by a Kernel-based Collocation Method