Bayesian learning with Gaussian processes for low-dimensional representations of time-dependent nonlinear systems

Shane A. McQuarrie,Anirban Chaudhuri,Karen E. Willcox,Mengwu Guo
2024-08-07
Abstract:This work presents a data-driven method for learning low-dimensional time-dependent physics-based surrogate models whose predictions are endowed with uncertainty estimates. We use the operator inference approach to model reduction that poses the problem of learning low-dimensional model terms as a regression of state space data and corresponding time derivatives by minimizing the residual of reduced system equations. Standard operator inference models perform well with accurate training data that are dense in time, but producing stable and accurate models when the state data are noisy and/or sparse in time remains a challenge. Another challenge is the lack of uncertainty estimation for the predictions from the operator inference models. Our approach addresses these challenges by incorporating Gaussian process surrogates into the operator inference framework to (1) probabilistically describe uncertainties in the state predictions and (2) procure analytical time derivative estimates with quantified uncertainties. The formulation leads to a generalized least-squares regression and, ultimately, reduced-order models that are described probabilistically with a closed-form expression for the posterior distribution of the operators. The resulting probabilistic surrogate model propagates uncertainties from the observed state data to reduced-order predictions. We demonstrate the method is effective for constructing low-dimensional models of two nonlinear partial differential equations representing a compressible flow and a nonlinear diffusion-reaction process, as well as for estimating the parameters of a low-dimensional system of nonlinear ordinary differential equations representing compartmental models in epidemiology.
Numerical Analysis
What problem does this paper attempt to address?
The paper aims to address the following issues: 1. **Prediction Uncertainty**: Traditional Operator Inference-based methods cannot provide uncertainty estimates for prediction results when dealing with nonlinear dynamical systems. 2. **Stability under Sparse or Noisy Data**: When state data is sparse or noisy, traditional operator inference methods struggle to produce stable and accurate models. To tackle these challenges, the authors propose a new method that combines Gaussian Processes (GP) with the operator inference framework, called GP-BayesOpInf. This method introduces Gaussian Processes to describe the uncertainty in state predictions and provides analytical time derivative estimates and their uncertainty quantification. Specifically, this method can learn physics-based surrogate models in low-dimensional representations and endow these predictions with uncertainty estimates. Additionally, the method can be effectively applied to parameter estimation in two nonlinear partial differential equations and a class of nonlinear ordinary differential equation systems in epidemiology.