An Efficient Bayesian Calibration Approach Using Dynamically Biorthogonal Field Equations

Author(s):  
Piyush M. Tagade ◽  
Han-Lim Choi

Present paper proposes new dynamic-biorthogonality based Bayesian formulation for calibration of computer simulators with parametric uncertainty. The formulation uses decomposition of solution field into mean and random field. The random field is represented as a convolution of separable Hilbert spaces in stochastic and spacial dimensions. Both the dimensions are spectrally represented using respective orthogonal bases. In particular, present paper investigates polynomial chaos basis for stochastic dimension and eigenfunction basis for spacial dimension. Dynamic evolution equations are derived such that basis in stochastic dimension is retained while basis in spacial dimension is changed such that dynamic orthogonality is maintained. Resultant evolution equations are used to propagate prior uncertainty in input parameters to the solution output. Whenever new information is available through experimental observations or expert opinion, Bayes theorem is used to update the basis in stochastic dimension. Efficacy of the proposed methodology is demonstrated for calibration of 2D transient diffusion equation with uncertainty in source location. Computational efficiency of the method is demonstrated against Generalized Polynomial Chaos and Monte Carlo method.

2017 ◽  
Vol 27 (2) ◽  
pp. 229-243
Author(s):  
Piyush M. Tagade ◽  
Han-Lim Choi

AbstractThis paper proposes a novel computationally efficient stochastic spectral projection based approach to Bayesian inversion of a computer simulator with high dimensional parametric and model structure uncertainty. The proposed method is based on the decomposition of the solution into its mean and a random field using a generic Karhunen-Loève expansion. The random field is represented as a convolution of separable Hilbert spaces in stochastic and spatial dimensions that are spectrally represented using respective orthogonal bases. In particular, the present paper investigates generalized polynomial chaos bases for the stochastic dimension and eigenfunction bases for the spatial dimension. Dynamic orthogonality is used to derive closed-form equations for the time evolution of mean, spatial and the stochastic fields. The resultant system of equations consists of a partial differential equation (PDE) that defines the dynamic evolution of the mean, a set of PDEs to define the time evolution of eigenfunction bases, while a set of ordinary differential equations (ODEs) define dynamics of the stochastic field. This system of dynamic evolution equations efficiently propagates the prior parametric uncertainty to the system response. The resulting bi-orthogonal expansion of the system response is used to reformulate the Bayesian inference for efficient exploration of the posterior distribution. The efficacy of the proposed method is investigated for calibration of a 2D transient diffusion simulator with an uncertain source location and diffusivity. The computational efficiency of the method is demonstrated against a Monte Carlo method and a generalized polynomial chaos approach.


2019 ◽  
Vol 19 (1) ◽  
pp. 39-53 ◽  
Author(s):  
Martin Eigel ◽  
Johannes Neumann ◽  
Reinhold Schneider ◽  
Sebastian Wolf

AbstractThis paper examines a completely non-intrusive, sample-based method for the computation of functional low-rank solutions of high-dimensional parametric random PDEs, which have become an area of intensive research in Uncertainty Quantification (UQ). In order to obtain a generalized polynomial chaos representation of the approximate stochastic solution, a novel black-box rank-adapted tensor reconstruction procedure is proposed. The performance of the described approach is illustrated with several numerical examples and compared to (Quasi-)Monte Carlo sampling.


Sign in / Sign up

Export Citation Format

Share Document