An Empirical Interpolation and Model-Variance Reduction Method for Computing Statistical Outputs of Parametrized Stochastic Partial Differential Equations
Name
15m1016783.pdf
Size
916.78 KB
Format
Adobe PDF
Checksum (MD5)
2c84a986b89f38468c5d1f4dcf396b76
Author(s)
Giles, M. B.
Vidal-Codina, Ferran
Nguyen, Ngoc Cuong
Peraire, Jaime
Date Issued
January 2016
Journal
SIAM/ASA Journal on Uncertainty Quantification
Publisher
Society of Industrial and Applied Mathematics
Citation
Vidal-Codina, F. et al. “An Empirical Interpolation and Model-Variance Reduction Method for Computing Statistical Outputs of Parametrized Stochastic Partial Differential Equations.” SIAM/ASA Journal on Uncertainty Quantification 4, 1 (January 2016): 244–265 © 2016 Society for Industrial and Applied Mathematics
Version
Final published version
Abstract
We present an empirical interpolation and model-variance reduction method for the fast and reliable computation of statistical outputs of parametrized stochastic elliptic partial differential equations. Our method consists of three main ingredients: (1) the real-time computation of reduced basis (RB) outputs approximating high-fidelity outputs computed with the hybridizable discontinuous Galerkin (HDG) discretization; (2) the empirical interpolation for an efficient offline-online decoupling of the parametric and stochastic inuence; and (3) a multilevel variance reduction method that exploits the statistical correlation between the low-fidelity approximations and the high-fidelity HDG dis- cretization to accelerate the convergence of the Monte Carlo simulations. The multilevel variance reduction method provides efficient computation of the statistical outputs by shifting most of the computational burden from the high-fidelity HDG approximation to the RB approximations. Fur- thermore, we develop a posteriori error estimates for our approximations of the statistical outputs. Based on these error estimates, we propose an algorithm for optimally choosing both the dimensions of the RB approximations and the size of Monte Carlo samples to achieve a given error tolerance. In addition, we extend the method to compute estimates for the gradients of the statistical out- puts. The proposed method is particularly useful for stochastic optimization problems where many evaluations of the objective function and its gradient are required.
MIT Department
Massachusetts Institute of Technology. Department of Aeronautics and Astronautics
Terms of Use
Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use.
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1137/15M1016783
