Data–driven model reduction for the Bayesian solution of inverse problems
Journal Article
·
· International Journal for Numerical Methods in Engineering
- Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States); Office of Scientific and Technical Information (OSTI)
- Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States)
One of the prime challenges in the Bayesian solution of inverse problems governed by partial differential equations (PDEs) is the computational cost of repeatedly evaluating numerical PDE models, as required by Markov chain Monte Carlo (MCMC) methods for posterior sampling. This paper presents a data-driven projection-based model reduction technique to reduce this computational cost. The proposed technique has two distinctive features. To begin, the model reduction strategy is tailored to inverse problems: the snapshots used to construct the reduced-order model are computed adaptively from the posterior distribution. Posterior exploration and model reduction are thus pursued simultaneously. Second, to avoid repeated evaluations of the full-scale numerical model as in a standard MCMC method, we couple the full-scale model and the reduced-order model together in the MCMC algorithm. This maintains accurate inference while reducing its overall computational cost. In numerical experiments considering steady-state flow in a porous medium, the data-driven reduced-order model achieves better accuracy than a reduced-order model constructed using the classical approach. It also improves posterior sampling efficiency by several orders of magnitude compared with a standard MCMC method.
- Research Organization:
- Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR) (SC-21)
- Grant/Contract Number:
- SC0009297
- OSTI ID:
- 1557833
- Journal Information:
- International Journal for Numerical Methods in Engineering, Journal Name: International Journal for Numerical Methods in Engineering Journal Issue: 5 Vol. 102; ISSN 0029-5981
- Publisher:
- WileyCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Scalable posterior approximations for large-scale Bayesian inverse problems via likelihood-informed parameter and state reduction
An adaptive Hessian approximated stochastic gradient MCMC method
Scalable Algorithms for Inverse Problems With High-Dimensional Parameter Spaces
Journal Article
·
Mon Mar 28 20:00:00 EDT 2016
· Journal of Computational Physics
·
OSTI ID:1548326
An adaptive Hessian approximated stochastic gradient MCMC method
Journal Article
·
Wed Feb 03 19:00:00 EST 2021
· Journal of Computational Physics
·
OSTI ID:1853727
Scalable Algorithms for Inverse Problems With High-Dimensional Parameter Spaces
Technical Report
·
Mon Aug 25 00:00:00 EDT 2025
·
OSTI ID:2584729