An adaptive local reduced basis method for solving PDEs with uncertain inputs and evaluating risk

Zou, Zilong; Kouri, Drew; Aquino, Wilkins

doi:10.1016/j.cma.2018.10.028

Title: An adaptive local reduced basis method for solving PDEs with uncertain inputs and evaluating risk

Journal Article · 25 October 2018 · Computer Methods in Applied Mechanics and Engineering

DOI: https://doi.org/10.1016/j.cma.2018.10.028 · OSTI ID:1485450

Zou, Zilong ^[1]; Kouri, Drew ^[2]; Aquino, Wilkins ^[1]

Duke Univ., Durham, NC (United States)
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

Many physical systems are modeled using partial differential equations (PDEs) with uncertain or random inputs. For such systems, naively propagating a fixed number of samples of the input probability law (or an approximation thereof) through the PDE is often inadequate to accurately quantify the "risk" associated with critical system responses. In this paper, we develop a goal-oriented, adaptive sampling and local reduced basis approximation for PDEs with random inputs. Our method determines a set of samples and an associated (implicit) Voronoi partition of the parameter domain on which we build local reduced basis approximations of the PDE solution. The samples are selected in an adaptive manner using an a posteriori error indicator. A notable advantage of the proposed approach is that the computational cost of the approximation during the adaptive process remains constant. We provide theoretical error bounds for our approximation and numerically demonstrate the performance of our method when compared to widely used adaptive sparse grid techniques. In addition, we tailor our approach to accurately quantify the risk of quantities of interest that depend on the PDE solution. We demonstrate our method on an advection-diffusion example and a Helmholtz example. Keywords: stochastic partial differential equation, adaptive sampling, reduced basis, a posteriori error estimation, coherent risk measure, conditional value-at-risk.

View Accepted Manuscript (DOE)

Cite

Export

Save

Research Organization:: Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

Sponsoring Organization:: DARPA

Grant/Contract Number:: AC04-94AL85000

OSTI ID:: 1485450

Report Number(s):: SAND2018-12270J; 669167

Journal Information:: Computer Methods in Applied Mechanics and Engineering, Journal Name: Computer Methods in Applied Mechanics and Engineering Journal Issue: C Vol. 345; ISSN 0045-7825

Publisher:: ElsevierCopyright Statement

Country of Publication:: United States

Language:: English

References (34)

Data-driven model reduction for the Bayesian solution of inverse problems: DATA-DRIVEN MODEL REDUCTION FOR INVERSE PROBLEMS Cui, Tiangang; Marzouk, Youssef M.; Willcox, Karen E. International Journal for Numerical Methods in Engineering, Vol. 102, Issue 5 https://doi.org/10.1002/nme.4748	journal	August 2014
Dimension?Adaptive Tensor?Product Quadrature Gerstner, T.; Griebel, M. Computing, Vol. 71, Issue 1 https://doi.org/10.1007/s00607-003-0015-5	journal	August 2003
An adaptive Monte Carlo algorithm for computing mixed logit estimators Bastin, Fabian; Cirillo, Cinzia; Toint, Philippe L. Computational Management Science, Vol. 3, Issue 1 https://doi.org/10.1007/s10287-005-0044-y	journal	January 2006
Comparison Between Reduced Basis and Stochastic Collocation Methods for Elliptic Problems Chen, Peng; Quarteroni, Alfio; Rozza, Gianluigi Journal of Scientific Computing, Vol. 59, Issue 1 https://doi.org/10.1007/s10915-013-9764-2	journal	August 2013
Reduced Basis Approximation and a Posteriori Error Estimation for Affinely Parametrized Elliptic Coercive Partial Differential Equations: Application to Transport and Continuum Mechanics Rozza, G.; Huynh, D. B. P.; Patera, A. T. Archives of Computational Methods in Engineering, Vol. 15, Issue 3 https://doi.org/10.1007/s11831-008-9019-9	journal	May 2008
Reduced Basis Techniques for Stochastic Problems Boyaval, S.; Le Bris, C.; Lelièvre, T. Archives of Computational Methods in Engineering, Vol. 17, Issue 4 https://doi.org/10.1007/s11831-010-9056-z	journal	October 2010
Solution of stochastic partial differential equations using Galerkin finite element techniques Deb, Manas K.; Babuška, Ivo M.; Oden, J. Tinsley Computer Methods in Applied Mechanics and Engineering, Vol. 190, Issue 48 https://doi.org/10.1016/S0045-7825(01)00237-7	journal	September 2001
A reduced basis approach for variational problems with stochastic parameters: Application to heat conduction with variable Robin coefficient Boyaval, Sébastien; Bris, Claude Le; Maday, Yvon Computer Methods in Applied Mechanics and Engineering, Vol. 198, Issue 41-44 https://doi.org/10.1016/j.cma.2009.05.019	journal	September 2009
Accurate and efficient evaluation of failure probability for partial different equations with random input data Chen, Peng; Quarteroni, Alfio Computer Methods in Applied Mechanics and Engineering, Vol. 267 https://doi.org/10.1016/j.cma.2013.08.016	journal	December 2013
Sparse-grid, reduced-basis Bayesian inversion Chen, Peng; Schwab, Christoph Computer Methods in Applied Mechanics and Engineering, Vol. 297 https://doi.org/10.1016/j.cma.2015.08.006	journal	December 2015
An ‘empirical interpolation’ method: application to efficient reduced-basis discretization of partial differential equations Barrault, Maxime; Maday, Yvon; Nguyen, Ngoc Cuong Comptes Rendus Mathematique, Vol. 339, Issue 9 https://doi.org/10.1016/j.crma.2004.08.006	journal	November 2004
Sparse grid collocation schemes for stochastic natural convection problems Ganapathysubramanian, Baskar; Zabaras, Nicholas Journal of Computational Physics, Vol. 225, Issue 1 https://doi.org/10.1016/j.jcp.2006.12.014	journal	July 2007
An adaptive hierarchical sparse grid collocation algorithm for the solution of stochastic differential equations Ma, Xiang; Zabaras, Nicholas Journal of Computational Physics, Vol. 228, Issue 8 https://doi.org/10.1016/j.jcp.2009.01.006	journal	May 2009
An adaptive high-dimensional stochastic model representation technique for the solution of stochastic partial differential equations Ma, Xiang; Zabaras, Nicholas Journal of Computational Physics, Vol. 229, Issue 10 https://doi.org/10.1016/j.jcp.2010.01.033	journal	May 2010
A method for solving stochastic equations by reduced order models and local approximations Grigoriu, M. Journal of Computational Physics, Vol. 231, Issue 19 https://doi.org/10.1016/j.jcp.2012.06.013	journal	August 2012
On the efficacy of stochastic collocation, stochastic Galerkin, and stochastic reduced order models for solving stochastic problems Field, R. V.; Grigoriu, M.; Emery, J. M. Probabilistic Engineering Mechanics, Vol. 41 https://doi.org/10.1016/j.probengmech.2015.05.002	journal	July 2015
Adaptive Monte Carlo analysis for strongly nonlinear stochastic systems Shields, Michael D. Reliability Engineering & System Safety, Vol. 175 https://doi.org/10.1016/j.ress.2018.03.018	journal	July 2018
Coherent Measures of Risk Artzner, Philippe; Delbaen, Freddy; Eber, Jean-Marc Mathematical Finance, Vol. 9, Issue 3 https://doi.org/10.1111/1467-9965.00068	journal	July 1999
High-Order Collocation Methods for Differential Equations with Random Inputs Xiu, Dongbin; Hesthaven, Jan S. SIAM Journal on Scientific Computing, Vol. 27, Issue 3 https://doi.org/10.1137/040615201	journal	January 2005
A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data Babuška, Ivo; Nobile, Fabio; Tempone, Raúl SIAM Journal on Numerical Analysis, Vol. 45, Issue 3 https://doi.org/10.1137/050645142	journal	January 2007
A Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data Nobile, F.; Tempone, R.; Webster, C. G. SIAM Journal on Numerical Analysis, Vol. 46, Issue 5 https://doi.org/10.1137/060663660	journal	January 2008
An Anisotropic Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data Nobile, F.; Tempone, R.; Webster, C. G. SIAM Journal on Numerical Analysis, Vol. 46, Issue 5 https://doi.org/10.1137/070680540	journal	January 2008
Model Reduction for Large-Scale Systems with High-Dimensional Parametric Input Space Bui-Thanh, T.; Willcox, K.; Ghattas, O. SIAM Journal on Scientific Computing, Vol. 30, Issue 6 https://doi.org/10.1137/070694855	journal	January 2008
A Trust-Region Algorithm with Adaptive Stochastic Collocation for PDE Optimization under Uncertainty Kouri, D. P.; Heinkenschloss, M.; Ridzal, D. SIAM Journal on Scientific Computing, Vol. 35, Issue 4 https://doi.org/10.1137/120892362	journal	January 2013
A Weighted Reduced Basis Method for Elliptic Partial Differential Equations with Random Input Data Chen, Peng; Quarteroni, Alfio; Rozza, Gianluigi SIAM Journal on Numerical Analysis, Vol. 51, Issue 6 https://doi.org/10.1137/130905253	journal	January 2013
Response Statistics for Random Heterogeneous Microstructures Grigoriu, M. SIAM/ASA Journal on Uncertainty Quantification, Vol. 2, Issue 1 https://doi.org/10.1137/130921490	journal	January 2014
Inexact Objective Function Evaluations in a Trust-Region Algorithm for PDE-Constrained Optimization under Uncertainty Kouri, D. P.; Heinkenschloss, M.; Ridzal, D. SIAM Journal on Scientific Computing, Vol. 36, Issue 6 https://doi.org/10.1137/140955665	journal	January 2014
The ROMES Method for Statistical Modeling of Reduced-Order-Model Error Drohmann, Martin; Carlberg, Kevin SIAM/ASA Journal on Uncertainty Quantification, Vol. 3, Issue 1 https://doi.org/10.1137/140969841	journal	January 2015
Reduced Basis Methods for Uncertainty Quantification Chen, Peng; Quarteroni, Alfio; Rozza, Gianluigi SIAM/ASA Journal on Uncertainty Quantification, Vol. 5, Issue 1 https://doi.org/10.1137/151004550	journal	January 2017
Galerkin Finite Element Approximations of Stochastic Elliptic Partial Differential Equations Babuska, Ivo; Tempone, Raúl; Zouraris, Georgios E. SIAM Journal on Numerical Analysis, Vol. 42, Issue 2 https://doi.org/10.1137/S0036142902418680	journal	January 2004
The Wiener--Askey Polynomial Chaos for Stochastic Differential Equations Xiu, Dongbin; Karniadakis, George Em SIAM Journal on Scientific Computing, Vol. 24, Issue 2 https://doi.org/10.1137/S1064827501387826	journal	January 2002
Optimization of conditional value-at-risk Rockafellar, R. Tyrrell; Uryasev, Stanislav The Journal of Risk, Vol. 2, Issue 3 https://doi.org/10.21314/JOR.2000.038	journal	January 2000
An Adaptive Sampling Approach for Solving PDEs with Uncertain Inputs and Evaluating Risk Zou, Zilong; Kouri, Drew P.; Aquino, Wilkins 19th AIAA Non-Deterministic Approaches Conference https://doi.org/10.2514/6.2017-1325	conference	January 2017
A Review of Unified A Posteriori Finite Element Error Control Carstensen, C. Numerical Mathematics: Theory, Methods and Applications, Vol. 5, Issue 4 https://doi.org/10.4208/nmtma.2011.m1032	journal	June 2012

Cited By (1)

Adaptive greedy algorithms based on parameter-domain decomposition and reconstruction for the reduced basis method Jiang, Jiahua; Chen, Yanlai arXiv https://doi.org/10.48550/arxiv.1911.02214	preprint	January 2019