Commun. Comput. Phys., 16 (2014), pp. 1-34. |
Adaptive Bayesian Inference for Discontinuous Inverse Problems, Application to Hyperbolic Conservation Laws Alexandre Birolleau ^{1}, Gael Poette ^{2}, Didier Lucor ^{3*} 1 CEA, DAM, DIF, F-91297 Arpajon, France; UPMC Univ Paris 06, UMR 7190, Institut Jean Le Rond d'Alembert, F-75005 Paris, France.2 CEA, DAM, DIF, F-91297 Arpajon, France. 3 UPMC Univ Paris 06, UMR 7190, Institut Jean Le Rond d'Alembert, F-75005 Paris, France; CNRS, UMR 7190, Institut Jean Le Rond d'Alembert, F-75005 Paris, France. Received 24 January 2013; Accepted (in revised version) 7 November 2013 Available online 28 March 2014 doi:10.4208/cicp.240113.071113a 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 {\it 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. AMS subject classifications: 52B10, 65D18, 68U05, 68U07Notice: Undefined variable: pac in /var/www/html/readabs.php on line 165 Key words: Stochastic inverse problems, Bayesian inference, iterative generalized Polynomial Chaos, compressible gas dynamics. *Corresponding author. Email: alexandre.birolleau@cea.fr (A. Birolleau), gael.poette@cea.fr (G. Poette), didier.lucor@upmc.fr (D. Lucor) |