English
Related papers

Related papers: Multilevel Monte Carlo estimators for elliptic PDE…

200 papers

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…

Numerical Analysis · Mathematics 2022-12-07 Joakim Beck , Yang Liu , Erik von Schwerin , Raúl Tempone

This paper provides a framework in which multilevel Monte Carlo and continuous level Monte Carlo can be compared. In continuous level Monte Carlo the level of refinement is determined by an exponentially distributed random variable, which…

Numerical Analysis · Mathematics 2023-10-13 Cedric Aaron Beschle , Andrea Barth

We propose a multilevel Markov chain Monte Carlo (MCMC) method for the Bayesian inference of random field parameters in PDEs using high-resolution data. Compared to existing multilevel MCMC methods, we additionally consider level-dependent…

Numerical Analysis · Mathematics 2025-08-19 Pieter Vanmechelen , Geert Lombaert , Giovanni Samaey

The Richards' equation is a model for flow of water in unsaturated soils. The coefficients of this (nonlinear) partial differential equation describe the permeability of the medium. Insufficient or uncertain measurements are commonly…

Numerical Analysis · Mathematics 2020-03-10 Andrea Barth , Andreas Stein

The efficient approximation of quantity of interest derived from PDEs with lognormal diffusivity is a central challenge in uncertainty quantification. In this study, we propose a multilevel quasi-Monte Carlo framework to approximate…

Numerical Analysis · Mathematics 2025-08-06 Joakim Beck , Yang Liu , Erik von Schwerin , Raúl Tempone

We develop new multilevel Monte Carlo (MLMC) methods to estimate the expectation of the smallest eigenvalue of a stochastic convection-diffusion operator with random coefficients. The MLMC method is based on a sequence of finite element…

Numerical Analysis · Mathematics 2024-02-13 Tiangang Cui , Hans De Sterck , Alexander D. Gilbert , Stanislav Polishchuk , Robert Scheichl

The Multilevel Monte Carlo (MLMC) method has proven to be an effective variance-reduction statistical method for Uncertainty Quantification (UQ) in Partial Differential Equation (PDE) models, combining model computations at different levels…

Mathematical Software · Computer Science 2023-05-24 Santiago Badia , Jerrad Hampton , Javier Principe

We present a multilevel Monte Carlo (MLMC) method for the uncertainty quantification of variably saturated porous media flow that are modeled using the Richards' equation. We propose a stochastic extension for the empirical models that are…

Numerical Analysis · Mathematics 2019-03-22 Prashant Kumar , Carmen Rodrigo , Francisco J. Gaspar , Cornelis W. Oosterlee

This paper addresses optimization problems constrained by partial differential equations with uncertain coefficients. In particular, the robust control problem and the average control problem are considered for a tracking type cost…

Optimization and Control · Mathematics 2017-11-08 Andreas Van Barel , Stefan Vandewalle

In a previous paper (J. Comp. Phys. 230 (2011), 3668--3694), the authors proposed a new practical method for computing expected values of functionals of solutions for certain classes of elliptic partial differential equations with random…

Numerical Analysis · Mathematics 2018-04-03 Ivan G. Graham , Frances Y. Kuo , Dirk Nuyens , Rob Scheichl , Ian H. Sloan

Stochastic PDE 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…

Numerical Analysis · Mathematics 2022-10-07 Alexander D. Gilbert , Robert Scheichl

We consider the numerical solution of scalar, nonlinear degenerate convection-diffusion problems with random diffusion coefficient and with random flux functions. Building on recent results on the existence, uniqueness and continuous…

Analysis of PDEs · Mathematics 2013-11-08 U. Koley , N. H. Risebro , Ch. Schwab , F. Weber

The multilevel Monte Carlo (MLMC) method has proven to be an effective variance-reduction statistical method for Uncertainty quantification in PDE models. It combines approximations at different levels of accuracy using a hierarchy of…

Numerical Analysis · Mathematics 2019-11-28 Santiago Badia , Jerrad Hampton , Javier Principe

In this article we consider static Bayesian parameter estimation for partially observed diffusions that are discretely observed. We work under the assumption that one must resort to discretizing the underlying diffusion process, for…

Computation · Statistics 2017-01-23 Ajay Jasra , Kengo Kamatani , Kody J. H. Law , Yan Zhou

This manuscript presents a framework for using multilevel quadrature formulae to compute the solution of optimal control problems constrained by random partial differential equations. Our approach consists in solving a sequence of optimal…

Numerical Analysis · Mathematics 2025-05-19 Fabio Nobile , Tommaso Vanzan

In this paper, we present a generalisation of the Multilevel Monte Carlo (MLMC) method to a setting where the level parameter is a continuous variable. This Continuous Level Monte Carlo (CLMC) estimator provides a natural framework in PDE…

Numerical Analysis · Mathematics 2018-02-22 Gianluca Detommaso , Tim Dodwell , Rob Scheichl

We consider elliptic diffusion problems with a random anisotropic diffusion coefficient, where, in a notable direction given by a random vector field, the diffusion strength differs from the diffusion strength perpendicular to this notable…

Numerical Analysis · Mathematics 2018-02-13 Helmut Harbrecht , Marc Schmidlin

We introduce the multivariate decomposition finite element method for elliptic PDEs with lognormal diffusion coefficient $a=\exp(Z)$ where $Z$ is a Gaussian random field defined by an infinite series expansion $Z(\boldsymbol{y}) =…

Numerical Analysis · Mathematics 2021-09-28 Dong T. P. Nguyen , Dirk Nuyens

In this article we consider a Bayesian inverse problem associated to elliptic partial differential equations (PDEs) in two and three dimensions. This class of inverse problems is important in applications such as hydrology, but the…

Computation · Statistics 2014-12-16 Alex Beskos , Ajay Jasra , Ege Muzaffer , Andrew Stuart

The Multilevel Monte Carlo (MLMC) approach usually works well when estimating the expected value of a quantity which is a Lipschitz function of intermediate quantities, but if it is a discontinuous function it can lead to a much slower…

Numerical Analysis · Mathematics 2023-09-06 Michael B Giles