arrow
Volume 16, Issue 1
Adaptive Bayesian Inference for Discontinuous Inverse Problems, Application to Hyperbolic Conservation Laws

Alexandre Birolleau, Gaël Poëtte & Didier Lucor

Commun. Comput. Phys., 16 (2014), pp. 1-34.

Published online: 2014-07

Export citation
  • Abstract

Various works from the literature aimed at accelerating Bayesian inference in inverse problems. Stochastic spectral methods have been recently proposed as surrogate approximations of the forward uncertainty propagation model over the support of the prior distribution. These representations are efficient because they allow affordable simulation of a large number of samples from the posterior distribution. Unfortunately, they do not perform well when the forward model exhibits strong nonlinear behavior with respect to its input.
In this work, we first relate the fast (exponential) $L^2$-convergence of the forward approximation to the fast (exponential) convergence (in terms of Kullback-Leibler divergence) of the approximate posterior. In particular, we prove that in case the prior distribution is uniform, the posterior is at least twice as fast as the convergence rate of the forward model in those norms. The Bayesian inference strategy is developed in the framework of a stochastic spectral projection method. The predicted convergence rates are then demonstrated for simple nonlinear inverse problems of varying smoothness.
We then propose an efficient numerical approach for the Bayesian solution of inverse problems presenting strongly nonlinear or discontinuous system responses. This comes with the improvement of the forward model that is adaptively approximated by an iterative generalized Polynomial Chaos-based representation. The numerical approximations and predicted convergence rates of the former approach are compared to the new iterative numerical method for nonlinear time-dependent test cases of varying dimension and complexity, which are relevant regarding our hydrodynamics motivations and therefore regarding hyperbolic conservation laws and the apparition of discontinuities in finite time.

  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-16-1, author = {}, title = {Adaptive Bayesian Inference for Discontinuous Inverse Problems, Application to Hyperbolic Conservation Laws}, journal = {Communications in Computational Physics}, year = {2014}, volume = {16}, number = {1}, pages = {1--34}, abstract = {

Various works from the literature aimed at accelerating Bayesian inference in inverse problems. Stochastic spectral methods have been recently proposed as surrogate approximations of the forward uncertainty propagation model over the support of the prior distribution. These representations are efficient because they allow affordable simulation of a large number of samples from the posterior distribution. Unfortunately, they do not perform well when the forward model exhibits strong nonlinear behavior with respect to its input.
In this work, we first relate the fast (exponential) $L^2$-convergence of the forward approximation to the fast (exponential) convergence (in terms of Kullback-Leibler divergence) of the approximate posterior. In particular, we prove that in case the prior distribution is uniform, the posterior is at least twice as fast as the convergence rate of the forward model in those norms. The Bayesian inference strategy is developed in the framework of a stochastic spectral projection method. The predicted convergence rates are then demonstrated for simple nonlinear inverse problems of varying smoothness.
We then propose an efficient numerical approach for the Bayesian solution of inverse problems presenting strongly nonlinear or discontinuous system responses. This comes with the improvement of the forward model that is adaptively approximated by an iterative generalized Polynomial Chaos-based representation. The numerical approximations and predicted convergence rates of the former approach are compared to the new iterative numerical method for nonlinear time-dependent test cases of varying dimension and complexity, which are relevant regarding our hydrodynamics motivations and therefore regarding hyperbolic conservation laws and the apparition of discontinuities in finite time.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.240113.071113a}, url = {http://global-sci.org/intro/article_detail/cicp/7031.html} }
TY - JOUR T1 - Adaptive Bayesian Inference for Discontinuous Inverse Problems, Application to Hyperbolic Conservation Laws JO - Communications in Computational Physics VL - 1 SP - 1 EP - 34 PY - 2014 DA - 2014/07 SN - 16 DO - http://doi.org/10.4208/cicp.240113.071113a UR - https://global-sci.org/intro/article_detail/cicp/7031.html KW - AB -

Various works from the literature aimed at accelerating Bayesian inference in inverse problems. Stochastic spectral methods have been recently proposed as surrogate approximations of the forward uncertainty propagation model over the support of the prior distribution. These representations are efficient because they allow affordable simulation of a large number of samples from the posterior distribution. Unfortunately, they do not perform well when the forward model exhibits strong nonlinear behavior with respect to its input.
In this work, we first relate the fast (exponential) $L^2$-convergence of the forward approximation to the fast (exponential) convergence (in terms of Kullback-Leibler divergence) of the approximate posterior. In particular, we prove that in case the prior distribution is uniform, the posterior is at least twice as fast as the convergence rate of the forward model in those norms. The Bayesian inference strategy is developed in the framework of a stochastic spectral projection method. The predicted convergence rates are then demonstrated for simple nonlinear inverse problems of varying smoothness.
We then propose an efficient numerical approach for the Bayesian solution of inverse problems presenting strongly nonlinear or discontinuous system responses. This comes with the improvement of the forward model that is adaptively approximated by an iterative generalized Polynomial Chaos-based representation. The numerical approximations and predicted convergence rates of the former approach are compared to the new iterative numerical method for nonlinear time-dependent test cases of varying dimension and complexity, which are relevant regarding our hydrodynamics motivations and therefore regarding hyperbolic conservation laws and the apparition of discontinuities in finite time.

Alexandre Birolleau, Gaël Poëtte & Didier Lucor. (2020). Adaptive Bayesian Inference for Discontinuous Inverse Problems, Application to Hyperbolic Conservation Laws. Communications in Computational Physics. 16 (1). 1-34. doi:10.4208/cicp.240113.071113a
Copy to clipboard
The citation has been copied to your clipboard