Paper Info
Title
Enhanced Multi-Index Monte Carlo By Means Of Multiple Semicoarsened Multigrid For Anisotropic Diffusion Problems
Abstract
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 modeled as a random field. These random fields can be highly anisotropic. Efficient solvers, such as the multiple semicoarsened multigrid (MSG) method are required to compute solutions for various realizations of the uncertain material. The MSG method is an extension of the classic multigrid method, which 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 to multi-index Monte Carlo (MIMC). We present an unbiased MIMC method that reuses the MSG coarse solutions. 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 and unbiased estimation. We analyze 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.
Year
DOI
Venue
2021
10.1002/nla.2281
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS
Keywords
DocType
Volume
anisotropic diffusion problems, multi-index Monte Carlo, multiple semicoarsened multigrid
Journal
28
Issue
ISSN
Citations 
3
1070-5325
0
PageRank 
References 
Authors
0.34
0
3
Name
Order
Citations
PageRank
Robbe Pieterjan100.34
Dirk Nuyens216817.97
Stefan Vandewalle350162.63