Enhanced Multi-Index Monte Carlo by means of Multiple Semi-Coarsened Multigrid for Anisotropic Diffusion Problems

by   Pieterjan Robbe, et al.

In many models used in engineering and science, material properties are uncertain or spatially varying. For example, in geophysics, and porous media flow in particular, the uncertain permeability of the material is modelled as a random field. These random fields can be highly anisotropic. Efficient solvers, such as the Multiple Semi-Coarsened Multigrid (MSG) method, see [15, 16, 17], are required to compute solutions for various realisations of the uncertain material. The MSG method is an extension of the classic Multigrid method, that uses additional coarse grids that are coarsened in only a single coordinate direction. In this sense, it closely resembles the extension of Multilevel Monte Carlo (MLMC) [6] to Multi-Index Monte Carlo (MIMC) [9]. We present an unbiased MIMC method that reuses the MSG coarse solutions, similar to the work in [11]. Our formulation of the estimator can be interpreted as the problem of learning the unknown distribution of the number of samples across all indices, and unifies the previous work on adaptive MIMC [19] and unbiased estimation [18]. We analyse the cost of this new estimator theoretically and present numerical experiments with various anisotropic random fields, where the unknown coefficients in the covariance model are considered as hyperparameters. We illustrate its robustness and superiority over unbiased MIMC without sample reuse.


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

The subsurface flow is usually subject to uncertain porous media structu...

Smoothed Circulant Embedding with Applications to Multilevel Monte Carlo Methods for PDEs with Random Coefficients

We consider the computational efficiency of Monte Carlo (MC) and Multile...

Multilevel Monte Carlo estimators for elliptic PDEs with Lévy-type diffusion coefficient

General elliptic equations with spatially discontinuous diffusion coeffi...

On the Selection of Random Field Evaluation Points in the p-MLQMC Method

Engineering problems are often characterized by significant uncertainty ...

An Efficient Method for Uncertainty Propagation in Robust Software Performance Estimation

Software engineers often have to estimate the performance of a software ...

Robust Optimization of PDEs with Random Coefficients Using a Multilevel Monte Carlo Method

This paper addresses optimization problems constrained by partial differ...

A randomized Multi-index sequential Monte Carlo method

We consider the problem of estimating expectations with respect to a tar...

Please sign up or login with your details

Forgot password? Click here to reset