Gaussian Process Modeling of a Functional Output with Information from Boundary and Initial Conditions and Analytical Approximations

<p>A partial differential equation (PDE) models a physical quantity as a function of space and time. These models are often solved numerically with the finite element (FE) method and the computer output consists of values of the solution on a fine grid over the spatial and temporal domain. When the simulations are time-consuming, Gaussian process (GP) models can be used to approximate the relationship between the functional output and the computer inputs, which consists of boundary and initial conditions. The Dirichlet boundary and initial conditions give the functional output values on parts of the space-time domain boundary. Although this information can help improve prediction of the output, it has not been used to construct GP models. In addition, analytical solutions of the PDE derived by simplifying the PDE can often be obtained, which can help further improve performance of the GP model. This article proposes a Karhunen–Loève (KL) expansion-based GP model that satisfies the Dirichlet boundary and initial conditions almost surely, and effectively uses information from analytical approximations to the PDE solution. Real examples demonstrate the improved prediction performance achieved by using these sources of prior information. Supplementary materials for this article are available online.