Fast sampling of parameterised Gaussian random fields

04/30/2018
by   Jonas Latz, et al.
0

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-Loève expansions for sampling. To reduce the cost we construct a reduced basis surrogate built from snapshots of Karhunen-Loève eigenvectors. In particular, we consider Matérn-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 ∼ 10^4 degrees of freedom.

READ FULL TEXT

page 3

page 21

page 33

page 34

research
11/24/2022

Numerical Approximation of Gaussian random fields on Closed Surfaces

We consider the numerical approximation of Gaussian random fields on clo...
research
12/05/2021

On the Numerical Approximation of the Karhunen-Loève Expansion for Random Fields with Random Discrete Data

Many physical and mathematical models involve random fields in their inp...
research
08/29/2023

Robust topology optimisation of lattice structures with spatially correlated uncertainties

The uncertainties in material and other properties of structures are usu...
research
01/24/2020

Certified and fast computations with shallow covariance kernels

Many techniques for data science and uncertainty quantification demand e...
research
06/23/2023

Smoothed Circulant Embedding with Applications to Multilevel Monte Carlo Methods for PDEs with Random Coefficients

We consider the computational efficiency of Monte Carlo (MC) and Multile...
research
09/11/2017

A determinant-free method to simulate the parameters of large Gaussian fields

We propose a determinant-free approach for simulation-based Bayesian inf...
research
05/31/2019

Unified Analysis of Periodization-Based Sampling Methods for Matérn Covariances

The periodization of a stationary Gaussian random field on a sufficientl...

Please sign up or login with your details

Forgot password? Click here to reset