DeepAI AI Chat
Log In Sign Up

A full multigrid multilevel Monte Carlo method for the single phase subsurface flow with random coefficients

by   Yang Liu, et al.
King Abdullah University of Science and Technology

The subsurface flow is usually subject to uncertain porous media structures. In most cases, however, we only have partial knowledge about the porous media properties. A common approach is to model the uncertain parameters of porous media as random fields, then the statistical moments (e.g. expectation) of the Quantity of Interest(QoI) can be evaluated by the Monte Carlo method. In this study, we develop a full multigrid-multilevel Monte Carlo (FMG-MLMC) method to speed up the evaluation of random parameters effects on single-phase porous flows. In general, MLMC method applies a series of discretization with increasing resolution and computes the QoI on each of them. The effective variance reduction is the success of the method. We exploit the similar hierarchies of MLMC and multigrid methods and obtain the solution on coarse mesh Q^c_l as a byproduct of the full multigrid solution on fine mesh Q^f_l on each level l. In the cases considered in this work, the computational saving due to the coarse mesh samples saving is 20% asymptotically. Besides, a comparison of Monte Carlo and Quasi-Monte Carlo (QMC) methods reveals a smaller estimator variance and a faster convergence rate of the latter approach in this study.


page 1

page 2

page 3

page 4


Quasi continuous level Monte Carlo

Continuous level Monte Carlo is an unbiased, continuous version of the c...

Multilevel Monte Carlo Simulations of Composite Structures with Uncertain Manufacturing Defects

By adopting a Multilevel Monte Carlo (MLMC) framework, we show that only...

Solution of the 3D density-driven groundwater flow problem with uncertain porosity and permeability

As groundwater is an essential nutrition and irrigation resource, its po...

Constrained Ensemble Langevin Monte Carlo

The classical Langevin Monte Carlo method looks for i.i.d. samples from ...

Frozen Gaussian Sampling: A Mesh-free Monte Carlo Method For Approximating Semiclassical Schrödinger Equations

In this paper, we develop a Monte Carlo algorithm named the Frozen Gauss...

Learning Patterns in Sample Distributions for Monte Carlo Variance Reduction

This paper investigates a novel a-posteriori variance reduction approach...