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 language | English (US) |
---|---|
Pages (from-to) | 5761-5787 |
Number of pages | 27 |
Journal | Journal of Computational Physics |
Volume | 231 |
Issue number | 17 |
DOIs | |
State | Published - Jul 1 2012 |
Externally published | Yes |
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