Quantify uncertainty by estimating the probability density function of the output of interest using MLMC based Bayes method
- 1 January 2023
- journal article
- research article
- Published by American Institute of Mathematical Sciences (AIMS) in Discrete and Continuous Dynamical Systems-Series B
- Vol. 28 (1) , 707-721
Abstract
In uncertainty quantification, the quantity of interest is usually the statistics of the space and/or time integration of system solution. In order to reduce the computational cost, a Bayes estimator based on multilevel Monte Carlo (MLMC) is introduced in this paper. The cumulative distribution function of the output of interest, that is, the expectation of the indicator function, is estimated by MLMC method instead of the classic Monte Carlo simulation. Then, combined with the corresponding probability density function, the quantity of interest is obtained by using some specific quadrature rules. In addition, the smoothing of indicator function and Latin hypercube sampling are used to accelerate the reduction of variance. An elliptic stochastic partial differential equation is used to provide a research context for this model. Numerical experiments are performed to verify the advantage of computational reduction and accuracy improvement of our MLMC-Bayes method.Keywords
This publication has 41 references indexed in Scilit:
- Inverse problems: A Bayesian perspectiveActa Numerica, 2010
- Importance sampling: a reviewWIREs Computational Statistics, 2009
- A Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input DataSIAM Journal on Numerical Analysis, 2008
- Stochastic Methods for Flow in Porous Media: Coping with UncertaintiesVadose Zone Journal, 2006
- REPRESENTING MODEL UNCERTAINTY IN WEATHER AND CLIMATE PREDICTIONAnnual Review of Earth and Planetary Sciences, 2005
- The Bayesian brain: the role of uncertainty in neural coding and computationTrends in Neurosciences, 2004
- The case for probabilistic forecasting in hydrologyJournal of Hydrology, 2001
- A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output From a Computer CodeTechnometrics, 2000
- QUANTIFICATION OF UNCERTAINTY IN COMPUTATIONAL FLUID DYNAMICSAnnual Review of Fluid Mechanics, 1997
- Random field finite elementsInternational Journal for Numerical Methods in Engineering, 1986