Mesh-clustered Gaussian Process Emulator for Partial Differential Equation Boundary Value Problems

Chih-Li Sung,Wenjia Wang,Liang Ding,Xingjian Wang
DOI: https://doi.org/10.1080/00401706.2024.2320211
2024-01-01
Technometrics
Abstract:Partial differential equations (PDEs) have become an essential tool formodeling complex physical systems. Such equations are typically solvednumerically via mesh-based methods, such as finite element methods, withsolutions over the spatial domain. However, obtaining these solutions are oftenprohibitively costly, limiting the feasibility of exploring parameters in PDEs.In this paper, we propose an efficient emulator that simultaneously predictsthe solutions over the spatial domain, with theoretical justification of itsuncertainty quantification. The novelty of the proposed method lies in theincorporation of the mesh node coordinates into the statistical model. Inparticular, the proposed method segments the mesh nodes into multiple clustersvia a Dirichlet process prior and fits Gaussian process models with the samehyperparameters in each of them. Most importantly, by revealing the underlyingclustering structures, the proposed method can provide valuable insights intoqualitative features of the resulting dynamics that can be used to guidefurther investigations. Real examples are demonstrated to show that ourproposed method has smaller prediction errors than its main competitors, withcompetitive computation time, and identifies interesting clusters of mesh nodesthat possess physical significance, such as satisfying boundary conditions. AnR package for the proposed methodology is provided in an open repository.
What problem does this paper attempt to address?