Journal
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING
Volume 348, Issue -, Pages 978-1012Publisher
ELSEVIER SCIENCE SA
DOI: 10.1016/j.cma.2019.02.003
Keywords
Uncertainty quantification; Bayesian inverse problem; Reduced basis methods; Spatial statistics; Karhunen-Loeve expansion
Funding
- Isaac Newton Institute for Mathematical Sciences, UK
- DFG, Germany through the International Graduate School of Science and Engineering at the Technical University of Munich [10.02 BAYES]
- EPSRC, UK [EP/K032208/1]
- EPSRC [EP/K032208/1] Funding Source: UKRI
Ask authors/readers for more resources
Gaussian random fields are popular models for spatially varying uncertainties, arising for instance in geotechnical engineering, hydrology or image processing. A Gaussian random field is fully characterised by its mean function and covariance operator. In more complex models these can also be partially unknown. In this case we need to handle a family of Gaussian random fields indexed with hyperparameters. Sampling for a fixed configuration of hyperparameters is already very expensive due to the nonlocal nature of many classical covariance operators. Sampling from multiple configurations increases the total computational cost severely. In this report we employ parameterised Karhunen-Loeve expansions for sampling. To reduce the cost we construct a reduced basis surrogate built from snapshots of Karhunen-Loeve eigenvectors. In particular, we consider Matern-type covariance operators with unknown correlation length and standard deviation. We suggest a linearisation of the covariance function and describe the associated online-offline decomposition. In numerical experiments we investigate the approximation error of the reduced eigenpairs. As an application we consider forward uncertainty propagation and Bayesian inversion with an elliptic partial differential equation where the logarithm of the diffusion coefficient is a parameterised Gaussian random field. In the Bayesian inverse problem we employ Markov chain Monte Carlo on the reduced space to generate samples from the posterior measure. All numerical experiments are conducted in 2D physical space, with non-separable covariance operators, and finite element grids with similar to 10(4 )degrees of freedom. (C) 2019 Elsevier B.V. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available