Abstract:The local size of computational grids used in partial differential equation (PDE)-based probabilistic inverse problems can have a tremendous impact on the numerical results. As a consequence, numerical model identification procedures used in structural or material engineering may yield erroneous, mesh-dependent result. In this work, we attempt to connect the field of adaptive methods for deterministic and forward probabilistic finite-element (FE) simulations and the field of FE-based Bayesian inference. In particular, our target setting is that of exact inference, whereby complex posterior distributions are to be sampled using advanced Markov Chain Monte Carlo (MCMC) algorithms. Our proposal is for the mesh refinement to be performed in a goal-oriented manner. We assume that we are interested in a finite subset of quantities of interest (QoI) such as a combination of latent uncertain parameters and/or quantities to be drawn from the posterior predictive distribution. Next, we evaluate the quality of an approximate inversion with respect to these quantities. This is done by running two chains in parallel: (i) the approximate chain and (ii) an enhanced chain whereby the approximate likelihood function is corrected using an efficient deterministic error estimate of the error introduced by the spatial discretisation of the PDE of interest. One particularly interesting feature of the proposed approach is that no user-defined tolerance is required for the quality of the QoIs, as opposed to the deterministic error estimation setting. This is because our trust in the model, and therefore a good measure for our requirement in terms of accuracy, is fully encoded in the prior. We merely need to ensure that the finite element approximation does not impact the posterior distributions of QoIs by a prohibitively large amount. We will also propose a technique to control the error introduced by the MCMC sampler, and demonstrate the validity of the combined mesh and algorithmic quality control strategy.

Numerical posterior distribution error control and expected Bayes Factors in the bayesian Uncertainty Quantification of Inverse Problems

Posterior distribution existence and error control in Banach spaces in the Bayesian approach to UQ in inverse problems

Probabilistic Numerical Methods for PDE-constrained Bayesian Inverse Problems

A-posteriori QMC-FEM error estimation for Bayesian inversion and optimal control with entropic risk measure

On the convergence of the Laplace approximation and noise-level-robustness of Laplace-based Monte Carlo methods for Bayesian inverse problems

Bayesian Numerical Methods for Nonlinear Partial Differential Equations

A Bayesian Approach to Modeling Finite Element Discretization Error

A-posteriori error estimates for inverse problems

Generalized Bayes Approach to Inverse Problems with Model Misspecification

Inverse Problems in a Bayesian Setting

Adaptivity in Bayesian Inverse Finite Element Problems: Learning and Simultaneous Control of Discretisation and Sampling Errors

Convergence of Probability Densities using Approximate Models for Forward and Inverse Problems in Uncertainty Quantification

Statistical guarantees for Bayesian uncertainty quantification in non-linear inverse problems with Gaussian process priors

Solving linear Bayesian inverse problems using a fractional total variation-Gaussian (FTG) prior and transport map

Nonlinear Expectation Inference for Direct Uncertainty Quantification of Nonlinear Inverse Problems

Bayesian Approach to Inverse Problems for Functions with A Variable-Index Besov Prior

Multilevel Sequential${}^2$ Monte Carlo for Bayesian Inverse Problems

Posterior-Variance–Based Error Quantification for Inverse Problems in Imaging

Maximum-a-posteriori estimation with Bayesian confidence regions

Towards a variational principle for motivated vehicle motion

Bayesian probabilistic propagation of hybrid uncertainties: Estimation of response expectation function, its variable importance and bounds