Physics-Information-Aided Kriging: Constructing Covariance Functions
  using Stochastic Simulation Models
  In this work, we propose a new Gaussian process regression (GPR) method:
physics information aided Kriging (PhIK). In the standard data-driven Kriging,
the unknown function of interest is usually treated as a Gaussian process with
assumed stationary covariance with hyperparameters estimated from data. In
PhIK, we compute the mean and covariance function from realizations of
available stochastic models, e.g., from realizations of governing stochastic
partial differential equations solutions. Such constructed Gaussian process
generally is non-stationary, and does not assume a specific form of the
covariance function. Our approach avoids the optimization step in data-driven
GPR methods to identify the hyperparameters. More importantly, we prove that
the physical constraints in the form of a deterministic linear operator are
guaranteed in the resulting prediction. We also provide an error estimate in
preserving the physical constraints when errors are included in the stochastic
model realizations. To reduce the computational cost of obtaining stochastic
model realizations, we propose a multilevel Monte Carlo estimate of the mean
and covariance functions. Further, we present an active learning algorithm that
guides the selection of additional observation locations. The efficiency and
accuracy of PhIK are demonstrated for reconstructing a partially known modified
Branin function, studying a three-dimensional heat transfer problem and
learning a conservative tracer distribution from sparse concentration
measurements.
