Modern Monte Carlo methods for efficient uncertainty quantification and propagation: A survey

Abstract

Uncertainty quantification (UQ) includes the characterization, integration, and propagation of uncertainties that result from stochastic variations and a lack of knowledge or data in the natural world. Monte Carlo (MC) method is a sampling‐based approach that has widely used for quantification and propagation of uncertainties. However, the standard MC method is often time‐consuming if the simulation‐based model is computationally intensive. This article gives an overview of modern MC methods to address the existing challenges of the standard MC in the context of UQ. Specifically, multilevel Monte Carlo (MLMC) extending the concept of control variates achieves a significant reduction of the computational cost by performing most evaluations with low accuracy and corresponding low cost, and relatively few evaluations at high accuracy and corresponding high cost. Multifidelity Monte Carlo (MFMC) accelerates the convergence of standard Monte Carlo by generalizing the control variates with different models having varying fidelities and varying computational costs. Multimodel Monte Carlo method (MMMC), having a different setting of MLMC and MFMC, aims to address the issue of UQ and propagation when data for characterizing probability distributions are limited. Multimodel inference combined with importance sampling is proposed for quantifying and efficiently propagating the uncertainties resulting from small data sets. All of these three modern MC methods achieve a significant improvement of computational efficiency for probabilistic UQ, particularly uncertainty propagation. An algorithm summary and the corresponding code implementation are provided for each of the modern MC methods. The extension and application of these methods are discussed in detail. This article is categorized under: Statistical and Graphical Methods of Data Analysis > Monte Carlo Methods Statistical and Graphical Methods of Data Analysis > Sampling

B1inhomogeneity correction of RARE MRI at low SNR: Quantitative in vivo19F MRI of mouse neuroinflammation with a cryogenically‐cooled transceive surface radiofrequency probe

“…Given the sparse nature of 19 F images and the spatially varying B 1 fields of the 19 F‐CRP, we computed concentration uncertainty maps after B 1 correction as follows (Figure 2): Step 1.Monte Carlo SNR simulations 34,35 (1000 iterations) were performed using measured (T 1 values) and synthetic data (SI computed using the simulated RARE SI model). Simulation parameters (Table 2) were defined to mimic realistic excitation FAs, ‐values, and SNRs within the sample.…”

“…Step 1. Monte Carlo SNR simulations 34,35 (1000 iterations) were performed using measured (T 1 values) and synthetic data (SI computed using the simulated RARE SI model). Simulation parameters (Table 2) were defined to mimic realistic excitation FAs, B − 1 -values, and SNRs within the sample.…”

Uncertainty Propagation via Optimal Transport Ambiguity Sets

“…Since exact propagation is feasible only for a very limited class of processes (e.g., Gaussian distributions and linear processes), numerous approaches have been investigated in the literature. Monte Carlo methods (and more generally simulation-based methods) rely on the propagation of a finite number of realizations, which can then be used to infer the statistics of the propagated probability distribution [47,48].…”

Simulation Intelligence: Towards a New Generation of Scientific Methods

“…In some cases we may be able to build a surrogate model of the ABM simulation (i.e. statistical emulation) for efficient execution and uncertainty quantification [266,267].…”