A fast and scalable computational framework for large-scale and high-dimensional Bayesian optimal experimental design

10/28/2020
by   Keyi Wu, et al.
0

We develop a fast and scalable computational framework to solve large-scale and high-dimensional Bayesian optimal experimental design problems. In particular, we consider the problem of optimal observation sensor placement for Bayesian inference of high-dimensional parameters governed by partial differential equations (PDEs), which is formulated as an optimization problem that seeks to maximize an expected information gain (EIG). Such optimization problems are particularly challenging due to the curse of dimensionality for high-dimensional parameters and the expensive solution of large-scale PDEs. To address these challenges, we exploit two essential properties of such problems: the low-rank structure of the Jacobian of the parameter-to-observable map to extract the intrinsically low-dimensional data-informed subspace, and the high correlation of the approximate EIGs by a series of approximations to reduce the number of PDE solves. We propose an efficient offline-online decomposition for the optimization problem: an offline stage of computing all the quantities that require a limited number of PDE solves independent of parameter and data dimensions, and an online stage of optimizing sensor placement that does not require any PDE solve. For the online optimization, we propose a swapping greedy algorithm that first construct an initial set of sensors using leverage scores and then swap the chosen sensors with other candidates until certain convergence criteria are met. We demonstrate the efficiency and scalability of the proposed computational framework by a linear inverse problem of inferring the initial condition for an advection-diffusion equation, and a nonlinear inverse problem of inferring the diffusion coefficient of a log-normal diffusion equation, with both the parameter and data dimensions ranging from a few tens to a few thousands.

READ FULL TEXT

page 17

page 18

page 20

page 21

page 22

research
01/20/2022

Derivative-informed projected neural network for large-scale Bayesian optimal experimental design

We address the solution of large-scale Bayesian optimal experimental des...
research
01/27/2023

A Greedy Sensor Selection Algorithm for Hyperparameterized Linear Bayesian Inverse Problems

We consider optimal sensor placement for a family of linear Bayesian inv...
research
11/30/2020

Derivative-Informed Projected Neural Networks for High-Dimensional Parametric Maps Governed by PDEs

Many-query problems, arising from uncertainty quantification, Bayesian i...
research
02/19/2018

Goal-Oriented Optimal Design of Experiments for Large-Scale Bayesian Linear Inverse Problems

We develop a framework for goal oriented optimal design of experiments (...
research
12/18/2019

Optimal experimental design under irreducible uncertainty for inverse problems governed by PDEs

We present a method for computing A-optimal sensor placements for infini...

Please sign up or login with your details

Forgot password? Click here to reset