1. Introduction
In the pursuit of enhancing the performance and efficiency of gas turbine engines operating with a Joule–Brayton cycle, the next generation of engines demands higher overall pressure ratios.
Figure 1 shows the evolution of engine pressure ratio over time showcasing this trend. Current flying engines from major manufacturers are indicated with crosses in the figures, illustrating the ongoing efforts in this direction.
This drive toward increased higher overall pressure ratios has led to the design of compressors with reduced blade heights, resulting in a more compact core. Increasing the overall pressure ratio at the compressor outlet, in turn, requires a reduction in the channel core dimensions due to the need to maintain a fixed Mach number at the combustion chamber inlet. This reduction leads to a subsequent decrease in compressor blade height. This compactness poses challenges for traditional measurement devices, which struggle to match the miniaturization process, causing increased blockage and variations in the local flow field.
The quest for enhancing the performance and efficiency of gas turbine engines demands innovative measurement techniques to overcome challenges posed by increased pressure ratios and compact core designs.
Efforts to address these challenges through machine learning data-driven approaches, particularly in Computational Fluid Dynamics (CFD) modeling, have shown promise [
1]. However, their application to real, noisy experimental measurements remains limited.
Traditional optical techniques and high-speed particle image velocimetry, while promising, face challenges in industrial and engine testing environments. Data-driven methods offer an alternative by requiring fewer measurements and less instrumentation, minimizing the instrumentation blockage effect without compromising accuracy.
Venturi et al. [
2] used DNS data; Bertream et al. [
3] used a fusion of CFD with experimental data to reconstruct flow fields with proper orthogonal decomposition, and Lou et al. [
4]; and Seshadri et al. [
5] demonstrated successful flow reconstruction and optimal sensor placement, in the turbomachinery field with optimization and Bayesian methodologies, respectively.
In response to current machine learning trends [
6] and recognizing the need to mitigate instrumentation impact and minimize data acquisition times, Gaussian Processes (GPs), or Kriging, has emerged as a key algorithm to deal with real data, due to this Bayesian formulation naturally dealing with uncertainty. GPs are a non-parametric Bayesian method that allows for the prediction of a mean unknown response and covariance, providing straightforward uncertainty propagation. GPs are typically used in the fluid mechanics literature as surrogate models, for example, to replace expensive CFD codes in blade design studies [
7].
This paper focuses on evaluating and comparing various Gaussian Process (GP) modeling approaches within a proposed hybrid measurement technique, aiming to mitigate instrumentation impact and minimize testing duration. The schematic representation of this hybrid methodology is illustrated in
Figure 2, depicting a process that seeks comprehensive information through localized experimental measurements, requiring significantly fewer samples.
The methodology begins by selecting an operating point (OP) and defining the essential boundary conditions for a CFD RANS simulation. The resulting flow field quantities across the turbomachinery domain guide a Design of Experiments (DoEs) framework to select optimal experimental acquisition locations. Experimental measurements acquired at these locations provide the basis for the subsequent GP modeling approaches, which will be systematically compared.
This paper seeks to evaluate the efficacy of different GP modeling approaches at achieving accurate flow assessments with a reduced number of measurement samples. The trained GP model allows for the inference of the flow field across the entire spatial domain, providing accurate assessments comparable to physical measurements, with the Bayesian framework facilitating uncertainty quantification.
The paper is outlined as it follows:
Section 2 presents the H25 axial compressor test case used in work, with both experimental and numerical tools being presented. In
Section 3, an overview of the different GP modeling approaches we propose is given.
Section 4 shows the mean flow results and a comparison study between the different modeling approaches. Finally, the work concludes and speculates on future implementations to the proposed measurement technique.
2. Test Case—H25 Axial Compressor
This work focuses on the H25 compressor stage, which is a low aspect ratio research stage developed under the LEMCOTEC project. The stage is intended to evaluate the impact of core size on the performance and flow field of the stage of a high-speed axial compressor for aircraft gas turbine applications. The stage is highly three-dimensional and features a reduced blade height, from which it derives its name, with the blade height measuring a mere 25 mm.
The configuration of the H25 axial compressor, characterized by its reduced blade height, makes it an ideal candidate for the methodology’s assessment since it aligns with the research goals of improving experimental flow assessment methods in next-generation compressor technologies that are more sensitive to instrumentation since their size is decreasing and probe minimization cannot follow.
More importantly, data availability played a major role in the selection of this test case. The H25 axial compressor has been extensively and comprehensibly experimentally tested by [
8], providing access and in-depth knowledge of the flow. These data include a wide array of measurements related to the flow characteristics of the compressor that are described below.
2.1. Experimental Setup
The experimental data were acquired in the R4 test rig of the von Karman Institute for Fluid Dynamics (VKI). The layout of this facility is depicted in
Figure 3. This facility incorporates an inlet plenum equipped with a heat exchanger and honeycomb (1), which regulate the total temperature and flow homogeneity. The flow is then channeled to the H25 test section (2), via a smooth convergent bell-mouth. After passing through the test section, the compressed flow is released into a large collector (3). A return duct (4) connects this collector back to the inlet plenum, forming a closed-loop system. The facility includes a precision throttle valve (5) that enables the accurate control of mass flow.
The R4 facility operates in a closed-loop configuration, providing the flexibility to independently adjust key parameters, including the Reynolds number (Re), rotational speed, and throttling transients. All of the experimental data used for this work were obtained under atmospheric conditions.
A meridional view of the H25 test section, labeled as point (2) in
Figure 3, is presented in
Figure 4. In this view, the rotating components are highlighted in red, while the stator row is marked in green. The instrumentation is axially distributed along four measurement planes, which are referenced as MP (Measurement Plane). MP0: Located after the convergent inlet bell-mouth, seven rotor chords upstream of the rotor. MP1: One chord upstream of the rotor. MP2: In the inter-row region between rotor and stator. MP4: At the stage outlet, one chord downstream of the stator blades.
The overall performance, particularly the pressure ratio and efficiency, is assessed using combined pressure-total temperature rakes deployed at MP0 and MP4. Additionally, the radial traversing of probes is possible in all measurement planes.
At MP4, a motor was installed to allow for the pitch-wise traversing of the probes over one and a half stator pitches. This setup enabled the sampling of a radial-azimuthal two-dimensional map of the flow field. The acquisition process for this outlet flow cartography, while essential for evaluating the flow inside the compressor, was the most time-consuming part of the experimental campaign, taking approximately 3 h to collect data points from 950 locations, distributed in a matrix of 25 radial points over 38 pitch-wise points, featuring refinement near the walls in the radial direction for improved gradient evaluation and uniform discretization in the tangential direction. For this reason, and given the critical role of MP4 in assessing the compressor’s performance, this experimental acquisition of the total pressure ratio flow field (P) at the machine’s design point serves as the test case for applying the proposed GP modeling approaches.
For this test, a custom-designed miniaturized three-hole probe was utilized to measure flow total pressure, angle, and Mach number. The systematic uncertainty of the miniaturized pressure probe is determined by the pressure scanner used. In this case, it amounts to ±0.08% of the full span. At nominal speed, the total pressure ratio uncertainty budget was estimated to be around 1%.
The probes were installed with a standard mechanical safety margin of 0.02 mm to ensure a gap between the probe and hub/casing walls. Additionally, the radial size of the probe head was added to this value when determining the accessible span for each probe. This led to a 2% span end-wall proximity.
2.2. Numerical Simulations—CFD RANS
A CFD numerical model is used to integrate highly discretized flow data with the undersampled experiment in the modeling approaches described above, and thus the domain matching between CFD and experiments is critical. The simulation domain is extended from experimental MP0 to a location two rotor chords downstream of MP4. This avoids issues related to outlet boundary conditions at MP4.
The numerical domain retains technological features like blade fillets in rotor and stator sections but excludes experimental cavities. Assuming periodicity, only one blade passage is simulated, utilizing a mixing plane approach for conservative pitchwise coupling.
Mesh generation employs a structured grid with Autogrid5, adopting an O-4H topology for blade rows and OH topology in the rotor tip gap. The non-dimensional wall distance () is maintained below 2.
Figure 5 illustrates the meridional and blade-to-blade views of the computational grid. A multi-block structured grid approach results in three mesh levels, with grid convergence analysis on 0.6, 4.55, and 19.45 million cells, selecting the 4.55 million cells mesh for the Reynolds-Averaged Navier–Stokes (RANS) data.
NUMECA FINE Turbo 12.1 solves fully turbulent compressible steady RANS equations using the k-
SST turbulence model [
9]. Operating at the compressor’s design point, the inlet imposes experimentally measured pitchwise averaged total pressure and total temperature profiles at MP0. At the outlet, absolute mass flow is imposed, with an adiabatic assumption for all domain walls.
3. Gaussian Process Modeling
The objective of this work is to retrieve missing or unavailable data using a limited set of experimental measurements through a Bayesian mathematical framework. This can be formulated as a regression problem:
where
represents the vector of observed output variables with noise-associated
.
is the matrix of input variables and
denotes the underlying regression function that maps inputs to outputs.
In the context of this work, the available experimental flow field measurements at measurement locations are used to build a regression model , enabling the estimation of unmeasured flow field data points at locations where probe measurements are unavailable . This repression model is based on Gaussian Processes.
Gaussian processes (GPs) are a non-parametric method that incorporates probabilistic distributions directly into the regression model
, constructing “Gaussian distributions over functions” [
10] rather than in a fixed number of parameters. This unique characteristic makes GPs a powerful probabilistic tool, especially for regression problems, as formulated in Equation (
1), for predicting unknown data
from observed data
, which include some level of measurement error or uncertainty.
A GP defines a multivariate Gaussian prior distribution over functions for every point in the input space
:
where
is a multivariate random variable
with a mean function
and a covariance matrix
that is defined with a positive definite kernel function. Due to the flexibility of GPs in modeling, the mean is assumed to be zero
and the GP distribution over
is fully defined by
. Therefore, the choice of the kernel function
is a key modeling decision in the GP framework.
Multiple kernel function families exist in the literature for different function modeling behaviors (e.g., smoothness, stationary), making this choice a prior input of knowledge of the problem. In the present work, the base kernel function selected is a stationary anisotropic Matérn-5/2 kernel, which is given by:
where
ℓ is the length scales for a single input dimension and controls the smoothness of the function modeling and
acts as a scaling factor. These are usually called hyperparameters
, and ensuring their correct values plays a pivotal role in GPs. The length scales
ℓ control the smoothness of the function by directly quantifying how much a single data point influences the space around it. The scaling factor
controls the vertical scale of the function around its mean, influencing the uncertainty of the model.
Given the above assumptions and available noisy observations
, where the noise
follows an independent, identically distributed Gaussian distribution with zero mean and fixed variance
, the log marginal likelihood can be derived as:
where the first term is a model complexity term, defined by the selection of the kernel function. The second term is a likelihood data-fit term showing the Bayesian weighting of the prior with the observed data. The third term is a constant term. Based on the available data, the log marginal likelihood can be maximized to obtain the optimal hyperparameters mentioned above.
The inherent noise of the available data is taken as a hyperparameter to infer along with the kernel parameters. This choice allows the model to determine the optimal level of noise during hyperparameter training. It enhances model flexibility by enabling adaptation to varying noise levels, and it provides valuable uncertainty information of what the model thinks of the data.
Predictions of the unknown states
can be obtained by conditioning their input locations
on the prior distribution of Equation (
2), giving a predictive posterior distribution:
where
represents the predicted mean values for new input locations
, and where
is the covariance matrix associated with the predicted values, computed, respectively, as:
3.1. Baseline GP Approach
The conventional single-fidelity approach to Gaussian Process modeling, here taken as a baseline, involves using a set of domain inputs,
, where
is a single location where a noisy experimental measurement
is taken. The unknown values,
, are the equivalent quantity of interest evaluated at non-sampled locations
.
This method provides a good approximation of unknown values close to input locations where data are available but tends to become less accurate as the inference moves away from these locations, for example, in regions where few measurements are available or none exist, like in the boundary layer region of the flow.
3.2. CFD to Experiments Approach
The approach proposed in [
11] involves coupling numerical data with experimental measurements in GPs. Instead of mapping the physical input domain
,
J RANS simulations flow fields (at various operating conditions of the machine) are used as inputs
. Using the case of total pressure at the compressor outlet, one input
is a j-dimensional vector that contains all the
J RANS total pressure results evaluated at a grid point
, where one experimental measurement
is available. The GP functional mapping of this work is represented by the following matrix format in Equation (
9) for the total pressure at the outlet of the compressor.
where the matrix
represents the inputs for the GP model when incorporating CFD data. Each column of this matrix corresponds to a location
in the physical domain, and the rows encapsulate the results of various RANS simulations at that location. Specifically,
denotes the total pressure obtained from the j-th RANS simulation at the i-th grid point
.
The use of an anisotropic kernel in the GP model plays a crucial role in this context. An anisotropic kernel allows the model to adapt its lengthscales independently for each input dimension. In the case of multiple CFD simulations, this flexibility is essential for balancing lengthscales between different simulations. It enables the GP model to account for variations in the physical domain along the radial and pitch directions separately, ensuring a more accurate representation of the complex flow field. This adaptability is particularly beneficial when dealing with simulations under different operating conditions as it allows the GP model to capture variations in the flow characteristics across the entire parameter space.
3.3. Multi-Fidelity Approach
The last approach, previously applied to fluid mechanics by multiple authors in [
12,
13], leverages the auto-regressive formulation proposed by O’Hagan [
14] for Multi-Fidelity GP modeling. The formulation assumes a linear dependency between each fidelity model. In a general formulation, it assigns a GP prior to each fidelity model
t, where the higher fidelity model
is expressed as a function of the lower fidelity
, multiplied by a scaling factor
plus a bias function
, which is itself a GP
.
The scaling factor weights the cross-correlation between fidelity levels and it is considered an hyperparameter to be learned in the inference process. The above formulation ensures that at each fidelity level t, the conditional distribution of the GP is influenced by thourgh the Markov property, , which means that, assuming is known, no more can be learned about from any other lower fidelity model output , for .
In the present case, only two levels of fidelity are available, the high fidelity experiments and the low fidelity CFD, which is the default case in this work, and assuming a constant scaling factor
, the formulation of Equation (
10) simplifies to:
where
and
represent the high- low-fidelity datasets, respectively, and
is a GP
. The prior GP model of the formulation of Equation (
11) can be expressed as:
where
and
are the mean functions.
is the covariance matrix associated with the low-fidelity data,
is the cross correlation matrix between fidelities, and
is the covariance matrix associated with the high-fidelity data.
The recursive approach of Le Gratiet [
15] decouples the inference of the fidelities into standard GP inference. The key step of the formulation is to infer the posterior of the low-fidelity model
independently and replace it with its prior in Equation (
11). With this formulation, the multi-fidelity GP posterior distribution for the high-fidelity
is fully defined, and predictions are made with following mean and covariance:
The proposed approach offers several benefits over the traditional GP regression approach. By using CFD as low fidelity data in the model, its high data density in the full domain enables the prediction of the unknown flow values at non-sampled domain locations where the experimental probe cannot reach. This makes it possible to infer the flow at regions where few or no measurements are available, such as the hub boundary layer in the case of the H25 test case.
In summary, three distinct, increasing in complexity, modeling approaches—single-fidelity experiments, single-fidelity CFD to experiments, and multi-fidelity modeling—were presented for flow field reconstruction from undersampled measurements. The subsequent sections will present a detailed analysis of the performance in accurately reconstructing flow fields, especially in challenging scenarios with limited measurements, analyzing each approach’s strengths, limitations, and practical implications, shedding light on their effectiveness in real-world applications.
4. Results
This section presents a comprehensive analysis of various GP modeling approaches, with the goal of evaluating their performance within a data-driven hybrid measurement technique presented in
Figure 2. An initial investigation involves a systematic comparison of all GP modeling approaches, utilizing a consistent measurement grid. This ensures a consistent evaluation of mean flow inference and its associated uncertainty across different approaches.
The data-driven hybrid measurement technique integrates a Design of Experiments (DoEs) phase, currently in development and not presented in this work. A preliminary application of this DoEs phase revealed that the optimal number of measurements is around one-third (33%) of the complete experimental reference test dataset.
The complete 950 experiment points dataset is divided into two distinct subsets: the training subset and the validation subset. The training subset corresponds to an optimal 33% measurement grid and is utilized in the optimization process, maximizing the marginal log-likelihood of Equation (
4), enabling hyper-parameter learning and flow reconstruction. The validation subset, corresponding to the remaining 67% measurements, serves as the ground truth and is compared with the predicted flow reconstruction to evaluate the performance of the GP modeling approach.
To assess the different GP modeling approaches’ performance in predicting flow characteristics, the following fundamental metrics are assessed:
Root Mean Squared Error (RMSE): RMSE is a standard error metric that quantifies prediction accuracy. It calculates the square root of the average squared differences between predicted and actual values.
where
represents a predicted value and
represents the corresponding ground truth validation measurement.
Maximum Absolute Error (MaxAE): MaxAE identifies the most significant prediction error, enabling a direct comparison with experimental measurement uncertainty. A notable difference between MaxAE and RMSE indicates problematic model performance in specific flow regions.
4.1. Baseline GP
This section presents a comprehensive analysis of the Baseline GP model within the context of the data-driven hybrid measurement technique, serving as a reference for subsequent comparisons with more advanced GP models. The baseline approach adopts a standard single-fidelity GP, associating each outlet pressure measurement with its domain location .
The systematic application of the baseline GP model on a consistent measurement grid provides insights into mean flow inference and its associated uncertainty. Trained on a carefully selected 33% subset of the complete experimental dataset, the model demonstrates its capability to reconstruct the flow field across the entire spatial domain.
Figure 6 illustrates the comparison between the predicted total pressure flow values (
), on the right, and the ground truth validation data (
). The selected experimental grid is represented with white crosses.
The Baseline GP model shows a smooth flow field prediction, revealing a generally consistent match between predicted and actual flows. It accurately captures high-pressure regions and secondary flow structures in terms of both shape and location. Notably, the model successfully captures the radial contouring of the strong gradient wake and other secondary flows, such as the boundary layer hub region (A) and losses emanating from the rotor tip vortex in the casing (B). These are correctly predicted in terms of size and magnitude.
However, there are notable differences, particularly in predicting the minimum pressure value in the wake. The GP model struggles to infer the correct wake lower pressure region core (C), where higher losses are expected. Another observation is the transition between the hub secondary flow structure and the blade wake region, which appears more pointwise in the prediction but occupies a larger pitch size in the experimental reference.
The Bayesian modeling enables direct uncertainty propagation, as illustrated in
Figure 7. This serves a dual purpose—showcasing the GP framework capability to predict the mean flow and its associated uncertainty, and providing insights into areas with potentially less certain predictions.
The model predicted uncertainty tends to increase rapidly outside the studied domain, which is particularly evident around the pitch close to zero. Furthermore, the model exhibits higher uncertainty in regions where no experimental measurements are available, as expected. The highest uncertainty within the domain of interest is observed in the strong gradient wake region, primarily around a radial span of 0.4, corresponding to the location of the wake core highlighted above. This increased uncertainty in regions of the wake and transition zones indicates that the model does not overstate its confidence in inherently challenging-to-predict flow regions. Importantly, even in regions with “high uncertainty”, the predicted uncertainty remains within the bounds of the full experimental test’s uncertainty. In other parts of the predicted flow, the uncertainty is comparatively lower.
The performance metrics presented in
Table 1 provide a quantitative assessment of the Baseline GP model. Despite specific challenges in certain flow regions, the model demonstrates strong predictive capabilities, with both RMSE and MaxAE values at magnitudes of 1 ×
or lower. Notably, the MaxAE discrepancy, observed in the wake core, is still within a reasonable range. Importantly, both errors are smaller than the experimental uncertainty, suggesting that the Baseline GP model provides predictions aligning well with the expected variability in experimental measurements. This is a positive indicator of the model’s robustness in practical applications.
In summary, the Baseline GP model, while encountering specific challenges in certain flow regions, serves as a strong foundation for the subsequent comparison with other modeling approaches. The performance metrics validate its capability to predict the flow field with high accuracy and provide a realistic representation of uncertainty.
4.2. CFD to Experiments
This subsection presents the CFD to Experiments GP modeling approach, as proposed in [
11]. Unlike the Baseline GP model, this approach introduces an innovative coupling of numerical data with experimental measurements within the GP framework.
In this approach, a total of J RANS simulations are executed, each corresponding to various operating conditions of the turbomachinery. These simulations provide comprehensive flow field data across the entire domain, effectively capturing the system’s behavior under different scenarios. Instead of directly mapping the physical input domain, each input becomes a j-dimensional vector containing all J RANS total pressure results evaluated at a specific grid point . This strategic choice aligns with the case of total pressure at the compressor outlet, where each experimental measurement is available.
Figure 8 visually compares the predicted total pressure flow values (
) using the CFD to Experiments GP approach (on the right) with the ground truth validation data (
) on the left.
The qualitative analysis of the CFD to Experiments GP model reveals distinctive features in the flow patterns. The model successfully captures key aspects, such as high-pressure regions and the wake. However, notable differences compared to the Baseline GP model are observed. Unlike the Baseline approach, the flow in this case exhibits sharper pointwise gradients in the radial contouring of the wake. Additionally, flow discrepancies appear larger, with the wake core not being correctly represented and its shape not aligning with the experimental reference (C).
In the tip region of the flow (B), the GP prediction displays a different flow feature shape, with the tip flow occupying a higher span and pitch area in the blade wake region and its having direction changed. This results in an overprediction of the pressure losses associated with this secondary flow.
Another effect of using CFD is the presence of a passage vortex signature in the hub region around pitch 0.5 in the predicted flow. This passage vortex is clearly visible in all CFD simulations at this span high, but it is not evident in the experimental reference (D). Despite extracting information from the CFD, this information might be incompatible with the experimental reference, hindering the GP modeling approach flow assessment capabilities.
Figure 9 extends the qualitative examination by illustrating the uncertainty in the predicted pressure ratio (
). This provides insights into the areas where the CFD to Experiments GP approach exhibits higher uncertainty.
Overall, the uncertainty levels are equivalent or higher than in the Baseline GP case, with the lowest and highest predicted uncertainty values being double. However, the regions of higher uncertainty differ. Unlike the Baseline GP, where the highest uncertainty was in the domain limits, in the CFD to Experiments approach, the regions of higher uncertainty are in the region of the blade wake and in the boundary layer regions, both at the hub and tip, where the mean flow differences were highlighted above. Notably, the wake radial gradient contouring is clearly detected as a region of high uncertainty. This is believed to be due to the fact that the different CFD simulations, at different operating points, have different wakes. Similarly, the passage vortex region, not present in the experiments but present in the CFD, shows high uncertainty as well.
Table 2 presents the performance metrics of the CFD to Experiments Gaussian Process (GP) modeling approach and compares them with the Baseline GP using 33% of the measurements.
The increase in both RMSE and MaxAE for the CFD to Experiments modeling highlights a less accurate representation of the flow field compared to the Baseline GP. Specifically, the CFD to Experiments model exhibits larger discrepancies in predicting total pressure values across the spatial domain, particularly in some important flow regions, with the MaxAE being one order of magnitude higher. While the Baseline GP demonstrated robust performance, the integration of CFD data in this current framework seems to introduce complexities leading to larger errors.
Differences observed in the CFD to Experiments model, such as sharper gradients and discrepancies in specific flow features, are reflected in the quantitative metrics. Utilizing CFD data, despite capturing certain flow characteristics, introduces uncertainties and discrepancies contributing to higher error metrics. Challenges in integrating CFD data with experimental measurements are evident in the notable increase in errors.
While being a novel approach, the present integration of CFD data introduces challenges, resulting in higher prediction errors. The comparison with the Baseline GP emphasizes the importance of carefully evaluating the impact of additional data sources on overall model performance. Looking ahead, the subsequent subsection explores the Multi-Fidelity GP model as an alternative, aiming to address the challenges highlighted in the CFD to Experiments modeling approach.
4.3. Multi-Fidelity GP
The section presents the Multi-Fidelity GP modeling approach, building upon the framework proposed in [
12]. The MFGP model integrates CFD data as low-fidelity information and experimental measurements as high-fidelity data to map the input domain space
to pressure measurements
. This approach introduces CFD information and complexity over the Baseline GP, aiming to enhance the accuracy of flow field predictions.
Figure 10 visually compares the Multi-Fidelity GP model predicted total pressure flow values (
) with the ground truth validation data (
). The model demonstrates a visually seamless prediction, accurately capturing high-pressure regions and secondary flow structures (A) and (B). A distinctive enhancement is observed in the prediction of the wake lower pressure region (C), signifying improved accuracy in areas with higher losses. Additionally, the transition between the hub secondary flow structure and the blade wake region is more accurately represented, addressing a previous pointwise discrepancy in the Baseline GP.
Another characteristic of the Multi-Fidelity GP model is its ability to smooth the flow, providing a super-resolution effect. This inherent smoothing facilitates flow assessment with reduced experimental measurements, indicating potential benefits in terms of testing time and instrumentation effort.
Figure 11 shows the uncertainty of the Multi-Fidelity GP model, demonstrating a close resemblance to the Baseline GP. The model uncertainty pattern aligns with expectations, showing a rapid increase beyond the studied domain and higher uncertainty in regions without experimental measurements.
Contrary to initial expectations, the contribution of CFD uncertainty to the overall uncertainty is minimal. The scaling factor hyperparameter , reflecting the relationship between CFD and experiments, is below 0.2, indicating a low contribution. This result is justified by a substantial discrepancy between RANS simulations and experiments, emphasizing the MFGP model’s adaptability when incorporating and propagating uncertainties in the presence of challenging data sources.
Table 3 provides a quantitative assessment of the Multi-Fidelity GP modeling approach, comparing its performance metrics with the Baseline GP.
The RMSE for the Multi-Fidelity GP is 4.5 × , showcasing a slight improvement of −4.3% compared to the Baseline GP. This reduction in RMSE indicates that the MFGP model provides a more accurate overall prediction of total pressure values, aligning closely with the experimental reference. The MaxAE, the Multi-Fidelity GP model, achieves a value of 1.6 × , representing a −5.9% decrease concerning the Baseline GP. The negative percentage change implies a reduction in the maximum prediction error.
The improvements observed in both RMSE and MaxAE highlight that the Multi-Fidelity GP model improved the capability to capture flow field characteristics compared to the Baseline GP. The model’s ability to better represent the wake core lower pressure region and refine the transition between flow structures contributes to these improvements.
Overall, the Multi-Fidelity GP approach demonstrates promising results, providing a more accurate and detailed representation of the flow field. The negative percentage changes in error metrics emphasize the advancement achieved by integrating CFD data as a supplementary information source in the present modeling framework. The Multi-Fidelity GP model’s effectiveness at mitigating the limitations observed in the Baseline GP positions is based in it being a valuable alternative for comprehensive flow field predictions in practical applications, which makes it an optimal data fusion approach for the proposed data-driven hybrid measurement technique in development.
4.4. Sensitivity Analysis
In this final segment of the results section, a comprehensive sensitivity analysis is conducted to evaluate the robustness and performance of the three GP modeling approaches—Baseline GP, CFD to Experiments, and Multi-Fidelity GP—under varying conditions of measurements subsampling and different measurement grids.
The sensitivity analysis evaluates subsampling percentages of 10%, 20%, 33%, and 50%, corresponding to absolute measurement points of 95, 190, 313, and 475, respectively. This systematic exploration aims to assess the models’ sensitivity to the quantity and distribution of experimental measurements, providing insights into their adaptability and performance across different measurement scenarios and grids.
Figure 12 presents the sensitivity of each modeling approach depicted through a series of RMSE box plots for different measurement subsampling percentages (10%, 20%, 33%, and 50%) and across different randomly generated measurement grids. The box plots provide a visual representation of the distribution of RMSE values, allowing for a comparative analysis of model performance under diverse measurement conditions. In these box plots, the whiskers span from the lowest to the highest values. The box itself represents the first and third quartiles, with the median shown as a line.
For all the different modeling approaches, an increase in the number of measurements leads to a reduction in RMSE. This reduction is particularly pronounced for the Baseline and the Multi-Fidelity GP when increasing the number of measurements. This behavior is due to the convergence of the hyperparameters, as discussed earlier. Once a sufficient number of measurements are used, the models effectively learn the correct data hyperparameters, and additional data do not significantly improve the modeling.
The results clearly indicate that the CFD to Experiments GP consistently performs worse than both the Baseline GP and Multi-Fidelity GP across all subsampling percentages. The box plots reveal higher median RMSE values for the CFD to Experiments GP, indicating a less accurate flow prediction compared to the other two models.
In contrast, the Baseline GP and Multi-Fidelity GP exhibit comparable performance. The box plots for both models show similar medians, suggesting a consistent level of accuracy across different subsampling percentages. However, the Multi-Fidelity GP stands out as a more robust option, as evidenced by its lower box plot spread and fewer outliers. This indicates that the Multi-Fidelity GP is less sensitive to variations in the experimental grid and maintains a stable performance even with reduced measurement data since the spread between the box plot whiskers can be primarily attributed to the influence of the randomly selected experimental measurement grids. The outliers correspond to random generated grids where the measurement points are clustered in regions of low-gradients. In these cases, the Multi-Fidelity GP can still rely on the gradients from the CFD, while the Baseline GP model cannot due to lack of information.
These findings highlight the advantages of the Multi-Fidelity GP model, showcasing its robustness and ability to provide accurate predictions with a reduced number of experimental measurements for the proposed data-driven Bayesian hybrid measurement technique in development.
5. Conclusions
In this work, three Gaussian Process (GP) modeling approaches—Baseline GP, CFD to Experiments GP, and Multi-Fidelity GP—are compared in the context of a data-driven hybrid measurement technique for turbomachinery applications. The goal was to evaluate their effectiveness in predicting mean flow characteristics and associated uncertainties, with a particular focus on informing the development of a hybrid measurement approach.
The Baseline GP served as a robust foundation, showcasing accurate predictions of the flow field and providing realistic uncertainty estimates. However, the novel integration of CFD data in the CFD into Experiments GP introduced complexities, resulting in more prediction errors and uncertainties, highlighting the need for careful consideration of additional data sources.
The Multi-Fidelity GP approach integrates CFD data as low-fidelity data and experimental measurements as high-fidelity data. It demonstrated enhanced accuracy in predicting critical flow features, such as the wake core and hub and casing regions.
The sensitivity analysis showed the Multi-Fidelity GP robustness, even with a reduced number of experimental measurements that might not be optimal. This is relevant in practical applications where resource constraints may limit the number of measurements. The Multi-Fidelity GP is, therefore, a reliable data fusion method for the proposed hybrid measurement technique, offering the potential to streamline instrumentation efforts and greatly reduce testing time.
This study lays the groundwork for further investigation into the scalability and practical implementation of these GP modeling approaches in real-world turbomachinery testing environments. Ongoing research is focusing on predicting the optimal number and location of measurement points based on CFD information. This method will then be experimentally validated in an industrial test case, to evaluate its feasibility and impact on reducing instrumentation efforts and testing time.