A multi-resolution, non-parametric, Bayesian framework for identification of spatially-varying model parameters

Abstract

This paper proposes a hierarchical, multi-resolution framework for the identification of model parameters and their spatially variability from noisy measurements of the response or output. Such parameters are frequently encountered in PDE-based models and correspond to quantities such as density or pressure fields, elasto-plastic moduli and internal variables in solid mechanics, conductivity fields in heat diffusion problems, permeability fields in fluid flow through porous media etc. The proposed model has all the advantages of traditional Bayesian formulations such as the ability to produce measures of confidence for the inferences made and providing not only predictive estimates but also quantitative measures of the predictive uncertainty. In contrast to existing approaches it utilizes a parsimonious, non-parametric formulation that favors sparse representations and whose complexity can be determined from the data. The proposed framework in non-intrusive and makes use of a sequence of forward solvers operating at various resolutions. As a result, inexpensive, coarse solvers are used to identify the most salient features of the unknown field(s) which are subsequently enriched by invoking solvers operating at finer resolutions. This leads to significant computational savings particularly in problems involving computationally demanding forward models but also improvements in accuracy. It is based on a novel, adaptive scheme based on Sequential Monte Carlo sampling which is embarrassingly parallelizable and circumvents issues with slow mixing encountered in Markov Chain Monte Carlo schemes.

References in corpus (2)

• Central limit theorem for sequential Monte Carlo methods and its application to Bayesian inference
• Comment on "Sequential Monte Carlo for Bayesian Computation" (P. Del Moral, A. Doucet, A. Jasra)

Cited by in corpus (13)

• Multi-fidelity classification using Gaussian processes: accelerating the prediction of large-scale computational models
• Multilevel Sequential${}^2$ Monte Carlo for Bayesian Inverse Problems
• SINDy-BVP: Sparse Identification of Nonlinear Dynamics for Boundary Value Problems
• Gaussian process regression + deep neural network autoencoder for probabilistic surrogate modeling in nonlinear mechanics of solids
• Bayesian multiscale deep generative model for the solution of high-dimensional inverse problems
• A variational Bayesian method for inverse problems with impulsive noise
• Multilevel Sequential Importance Sampling for Rare Event Estimation
• Parameter Identification in Viscoplasticity using Transitional Markov Chain Monte Carlo Method
• Calibration of multi-physics computational models using Bayesian networks
• Bayesian calibration of coupled computational mechanics models under uncertainty based on interface deformation
• Sparse Variational Bayesian Approximations for Nonlinear Inverse Problems: applications in nonlinear elastography
• Nonlinear elastodynamic material identification of heterogeneous isogeometric Bernoulli-Euler beams
• Multimodal, high-dimensional, model-based, Bayesian inverse problems with applications in biomechanics