Sampling-free linear Bayesian update of polynomial chaos representations

Bojana V. Rosić*, Alexander Litvinenko, Oliver Pajonk, Hermann G. Matthies

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

56 Scopus citations

Abstract

We present a fully deterministic approach to a probabilistic interpretation of inverse problems in which unknown quantities are represented by random fields or processes, described by possibly non-Gaussian distributions. The description of the introduced random fields is given in a " white noise" framework, which enables us to solve the stochastic forward problem through Galerkin projection onto polynomial chaos. With the help of such a representation the probabilistic identification problem is cast in a polynomial chaos expansion setting and the Baye's linear form of updating. By introducing the Hermite algebra this becomes a direct, purely algebraic way of computing the posterior, which is comparatively inexpensive to evaluate. In addition, we show that the well-known Kalman filter is the low order part of this update. The proposed method is here tested on a stationary diffusion equation with prescribed source terms, characterised by an uncertain conductivity parameter which is then identified from limited and noisy data obtained by a measurement of the diffusing quantity.

Original languageEnglish (US)
Pages (from-to)5761-5787
Number of pages27
JournalJournal of Computational Physics
Volume231
Issue number17
DOIs
StatePublished - Jul 1 2012
Externally publishedYes

Keywords

  • Kalman filter
  • Linear Bayesian update
  • Minimum squared error estimate
  • Minimum variance estimate
  • Polynomial chaos expansion

ASJC Scopus subject areas

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • General Physics and Astronomy
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Sampling-free linear Bayesian update of polynomial chaos representations'. Together they form a unique fingerprint.

Cite this