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.



page 3


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...

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 ...

Fast sampling of parameterised Gaussian random fields

Gaussian random fields are popular models for spatially varying uncertai...

Multilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise sampling

When solving partial differential equations with random fields as coeffi...

Stochastic phase-field modeling of brittle fracture: computing multiple crack patterns and their probabilities

In variational phase-field modeling of brittle fracture, the functional ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.