skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Nonlinear sparse Bayesian learning for physics-based models

Journal Article · · Journal of Computational Physics
 [1];  [2]; ORCiD logo [3];  [4];  [1]
  1. Carleton Univ., Ottawa, ON (Canada). Dept. of Civil & Environmental Engineering
  2. Sandia National Lab. (SNL-CA), Livermore, CA (United States). Quantitative Modeling & Analysis Dept.
  3. United States Naval Academy, Annapolis, MD (United States). Dept. of Aerospace Engineering
  4. Royal Military College of Canada, Kingston, ON (Canada). Dept. of Mechanical & Aerospace Engineering

This paper addresses the issue of overfitting while calibrating unknown parameters of over-parameterized physics-based models with noisy and incomplete observations. Here, a semi-analytical Bayesian framework of nonlinear sparse Bayesian learning (NSBL) is proposed to identify sparsity among model parameters during Bayesian inversion. NSBL offers significant advantages over machine learning algorithm of sparse Bayesian learning (SBL) for physics-based models, such as 1) the likelihood function or the posterior parameter distribution is not required to be Gaussian, and 2) prior parameter knowledge is incorporated into sparse learning (i.e. not all parameters are treated as questionable). NSBL employs the concept of automatic relevance determination (ARD) to facilitate sparsity among questionable parameters through parameterized prior distributions. The analytical tractability of NSBL is enabled by employing Gaussian ARD priors and by building a Gaussian mixture-model approximation of the posterior parameter distribution that excludes the contribution of ARD priors. Subsequently, type-II maximum likelihood is executed using Newton's method whereby the evidence and its gradient and Hessian information are computed in a semi-analytical fashion. We show numerically and analytically that SBL is a special case of NSBL for linear regression models. Subsequently, a linear regression example involving multimodality in both parameter posterior pdf and model evidence is considered to demonstrate the performance of NSBL in cases where SBL is inapplicable. Next, NSBL is applied to identify sparsity among the damping coefficients of a mass-spring-damper model of a shear building frame. These numerical studies demonstrate the robustness and efficiency of NSBL in alleviating overfitting during Bayesian inversion of nonlinear physics-based models.

Research Organization:
Sandia National Lab. (SNL-CA), Livermore, CA (United States)
Sponsoring Organization:
USDOE National Nuclear Security Administration (NNSA)
Grant/Contract Number:
AC04-94AL85000; NA0003525
OSTI ID:
1670717
Report Number(s):
SAND-2019-14718J; 682468; TRN: US2203877
Journal Information:
Journal of Computational Physics, Vol. 426; ISSN 0021-9991
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English

References (20)

Updating Models and Their Uncertainties. I: Bayesian Statistical Framework journal April 1998
Probabilistic parameter estimation of a fluttering aeroelastic system in the transitional Reynolds number regime journal July 2013
Global sensitivity analysis using polynomial chaos expansions journal July 2008
Tracking noisy limit cycle oscillation with nonlinear filters journal January 2010
Stochastic spectral methods for efficient Bayesian solution of inverse problems journal June 2007
Model Selection Using Response Measurements: Bayesian Probabilistic Approach journal February 2004
Bayesian model selection using automatic relevance determination for nonlinear dynamical systems journal June 2017
Weak convergence and optimal scaling of random walk Metropolis algorithms journal February 1997
Bayes factors: Prior sensitivity and model generalizability journal December 2008
Transitional Markov Chain Monte Carlo Method for Bayesian Model Updating, Model Class Selection, and Model Averaging journal July 2007
Bayesian calibration of computer models journal August 2001
Dimensionality reduction and polynomial chaos acceleration of Bayesian inference in inverse problems journal April 2009
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
Variational inference of cluster-weighted models for local and global sensitivity analysis journal January 2014
Calculation of Posterior Probabilities for Bayesian Model Class Assessment and Averaging from Posterior Samples Based on Dynamic System Data: Calculation of posterior probabilities journal February 2010
The Wiener--Askey Polynomial Chaos for Stochastic Differential Equations journal January 2002
The Homogeneous Chaos journal October 1938
Bayesian inference of nonlinear unsteady aerodynamics from aeroelastic limit cycle oscillations journal July 2016
Bayesian Compressive Sensing Using Laplace Priors journal January 2010
Bayesian system identification based on probability logic journal October 2010