Abstract:Uncertainty propagation in complex engineering systems often poses significant computational challenges related to modeling and quantifying probability distributions of model outputs, as those emerge as the result of various sources of uncertainty that are inherent in the system under investigation. Gaussian Processes regression (GPs) is a robust meta-modeling technique that allows for fast model prediction and exploration of response surfaces. Multi-fidelity variations of GPs further leverage information from cheap and low fidelity model simulations in order to improve their predictive performance on the high fidelity model. In order to cope with the high volume of data required to train GPs in high dimensional design spaces, a common practice is to introduce latent design variables that are typically projections of the original input space to a lower dimensional subspace, and therefore substitute the problem of learning the initial high dimensional mapping, with that of training a GP on a low dimensional space. In this paper, we present a Bayesian approach to identify optimal transformations that map the input points to low dimensional latent variables. The “projection” mapping consists of an orthonormal matrix that is considered a priori unknown and needs to be inferred jointly with the GP parameters, conditioned on the available training data. The proposed Bayesian inference scheme relies on a two-step iterative algorithm that samples from the marginal posteriors of the GP parameters and the projection matrix respectively, both using Markov Chain Monte Carlo (MCMC) sampling. In order to take into account the orthogonality constraints imposed on the orthonormal projection matrix, a Geodesic Monte Carlo sampling algorithm is employed, that is suitable for exploiting probability measures on manifolds. We extend the proposed framework to multi-fidelity models using GPs including the scenarios of training multiple outputs together. We validate our framework on three synthetic problems with a known lower-dimensional subspace. The benefits of our proposed framework, are illustrated on the computationally challenging aerodynamic optimization of a last-stage blade for an industrial gas turbine, where we study the effect of an 85-dimensional shape parameterization of a three-dimensional airfoil on two output quantities of interest, specifically on the aerodynamic efficiency and the degree of reaction.

Low-Dimensional Input and High-Dimensional Output Modelling Using Gaussian Process

Industrial Data Modeling with Low-Dimensional Inputs and High-Dimensional Outputs

Spatial Batch Optimal Design Based on Self-Learning Gaussian Process Models for LPCVD Processes

Towards Efficient Modeling and Inference in Multi-Dimensional Gaussian Process State-Space Models

A modified collocation modeling framework for dynamic evolution of molecular weight distributions in general polymer kinetic systems

Generalized initialization for the dynamic simulation and optimization of grade transition processes using two-dimensional collocation

Data-Driven Design via Scalable Gaussian Processes for Multi-Response Big Data With Qualitative Factors

A Fast Optimal Latin Hypercube Design for Gaussian Process Regression Modeling

Scalable Gaussian Processes for Data-Driven Design Using Big Data With Categorical Factors

Thin and Deep Gaussian Processes

Gaussian Process Emulation For Big Data In Data-Driven Metamaterials Design

Large Linear Multi-output Gaussian Process Learning

A Variant of Gaussian Process Dynamical Systems.

Multi-output Gaussian process prediction for computationally expensive problems with multiple levels of fidelity

Fast multi-fidelity gaussian processes with derivatives for complex system modeling

Mixed-effects Gaussian Process Modeling Approach with Application in Injection Molding Processes

Gradient-Enhanced Multi-Output Gaussian Process Model for Simulation-Based Engineering Design

A Graphical Multi-Fidelity Gaussian Process Model, with Application to Emulation of Heavy-Ion Collisions

Deep Learning Gaussian Processes For Computer Models with Heteroskedastic and High-Dimensional Outputs

Dimensionality Reduction for Multi-Fidelity Gaussian Processes Using Bayesian Adaptation

Bayesian Learning of Orthogonal Embeddings for Multi-Fidelity Gaussian Processes