Multifidelity Monte Carlo estimation for large-scale uncertainty propagation
Name
MFMCAIAASciTech2018.pdf
Size
1.04 MB
Format
Adobe PDF
Checksum (MD5)
eb8fd354048426b964019d92a8e2430f
Author(s) • •
Peherstorfer, Benjamin
Beran, Philip S.
Willcox, Karen E
Date Issued
January 2018
Journal
2018 AIAA Non-Deterministic Approaches Conference
Publisher
American Institute of Aeronautics and Astronautics
Citation
Peherstorfer, Benjamin, Philip S. Beran, and Karen E. Willcox. “Multifidelity Monte Carlo Estimation for Large-Scale Uncertainty Propagation.” 2018 AIAA Non-Deterministic Approaches Conference, 2018 AIAA Non-Deterministic Approaches Conference, 7 January, 2018, Kissimmee, Florida, AIAA, 2018.
Version
Author's final manuscript
Abstract
One important task of uncertainty quantification is propagating input uncertainties through a system of interest to quantify the uncertainties’ effects on the system outputs; however, numerical methods for uncertainty propagation are often based on Monte Carlo estimation, which can require large numbers of numerical simulations of the numerical model describing the system response to obtain estimates with acceptable accuracies. Thus, if the model is computationally expensive to evaluate, then Monte-Carlo-based uncertainty propagation methods can quickly become computationally intractable. We demonstrate that multifidelity methods can significantly speedup uncertainty propagation by leveraging low-cost low-fidelity models and establish accuracy guarantees by using occasional recourse to the expensive high-fidelity model. We focus on the multifidelity Monte Carlo method, which is a multifidelity approach that optimally distributes work among the models such that the mean-squared error of the multifidelity estimator is minimized for a given computational budget. The multifidelity Monte Carlo method is applicable to general types of low-fidelity models, including projection-based reduced models, data-fit surrogates, response surfaces, and simplified-physics models. We apply the multifidelity Monte Carlo method to a coupled aero-structural analysis of a wing and a flutter problem with a high-aspect-ratio wing. The low-fidelity models are data-fit surrogate models derived with standard procedures that are built in common software environments such as Matlab and numpy/scipy. Our results demonstrate speedups of orders of magnitude compared to using the high-fidelity model alone.
MIT Department
Massachusetts Institute of Technology. Department of Aeronautics and Astronautics
Terms of Use
Creative Commons Attribution-Noncommercial-Share Alike
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.2514/6.2018-1660