Accurate prediction of unsteady aerodynamic loads remains a major challenge in turbomachinery design. High-fidelity Computational Fluid Dynamics (CFD) simulations are expensive, while aeroelastic Quantities of Interest (QoI) depend sensitively on the temporal evolution of the pressure field. This work introduces periodic Fourier Neural Mapping (p-FNM), a neural-operator framework for predicting unsteady pressure distributions on turbine rotor blades simulated using the chorochronic numerical hypothesis. The architecture embeds temporal periodicity into the model and learns a continuous mapping from operating conditions and time to pressure fields. Unlike sequential latent-space approaches, p-FNM predicts pressure fields independently at any time, avoiding error accumulation while preserving temporal continuity. The model is evaluated on a database of unsteady rotor-blade simulations and compared with a reduced-order baseline based on a variational autoencoder and recurrent neural network, refered as the Temporal Prediction Model (TPM). Performance is assessed for pressure fields and Generalized Aerodynamic Forces (GAFs), the primary aeroelastic QoI. Across all training datasets, p-FNM consistently outperforms TPM. On the largest dataset, p-FNM achieves a pressure-field mean absolute percentage error of 0.46% and a GAF-magnitude prediction error of 4.42%, corresponding to improvements of 60.7% and 77.6%, respectively. The minimum weighted phase error reaches 0.060 rad, demonstrating accurate preservation of the temporal characteristics of the aerodynamic response. The results show that GAF prediction is more challenging than pressure-field prediction and that temporal coherence is critical for accurately predicting spectral aerodynamic quantities. These findings demonstrate the potential of periodic neural operators for reduced-order modeling and aeroelastic analysis in turbomachinery.
Figures & tables
Figure 1: Mid-span ( h/H=0.5 ) schlieren-like density gradient and static pressure field on the reference blade, showing the correlation between flow structures and aerodynamic loading
Table 1: Number of simulations for the different and independent datasets.
Figure 2: FNO layer architecture
Figure 3: p-FNM architecture (a) and V2F block (b)
Figure 4: Comparison of the TPM, TPM-SOAP and p-FNM models with respect to the pressure prediction metrics across the different training datasets. The metrics are evaluated on the test dataset. The black error bars represent the metric uncertainty quantified by the ±1.96σ confidence interval computed over 30 independent training runs.
Table 2: Comparison between predicted and reference projected-force signals for the VAE, TPM, and p-FNM models across the different training datasets for two representative test-cases from the test dataset. The black curve denotes the reference solution, whereas the orange, red, and green curves correspond to the VAE, TPM, and p-FNM predictions, respectively. For visual reference, each grid cell corresponds to 0.0005 units of projected force along the vertical axis.
Figure 5: Comparison of the TPM, and p-FNM models with respect to the temporal spectrum of the projected-force signal evaluated on the test dataset. All the Fourier harmonics of the blade-passing frequency used in the chorochronic boundary conditions are displayed.
Figure 6: Weighted Mean Absolute Percentage Error (wMAPE) associated with the prediction of the first Fourier mode of the projected-force for the TPM and p-FNM models across the different training datasets. The wMAPE is evaluated on the test dataset. The black error bars represent the metric uncertainty quantified by the ±1.96σ confidence interval computed over 30 independent training runs.
Figure 7: Comparison of the TPM, TPM-SOAP and p-FNM models with respect to the GAF prediction metrics across the different training datasets. The metrics are evaluated on the test dataset. The black error bars represent the metric uncertainty quantified by the ±1.96σ confidence interval computed over 30 independent training runs.
Table 3: Comparison between predicted and reference GAF norm values for the TPM and p-FNM models across the different training datasets. The models are evaluated on the test dataset. The black dashed line indicates perfect agreement, while the green and red dashed lines correspond to relative errors of ±5% and ±10% , respectively. The markers represent the mean predicted GAF values, and the gray error bars denote the associated prediction uncertainty quantified by the ±1.96σ confidence interval. Each grid cell corresponds to 0.0005 GAF units along both axes.
Figure 8: Comparison of the GAF phase prediction error distributions (in radians) for the TPM and p-FNM models across the different training datasets. The models are evaluated on the test dataset. The phase error is defined within the interval [−π,π] . Optimal performance corresponds to error distributions centered around zero with minimal dispersion.
Figure 9: Comparison of the TPM, TPM-SOAP, and p-FNM models in terms of training time (GPU ⋅ h) across the different training datasets. The black error bars indicate the training-time uncertainty quantified by the ±1.96σ confidence interval computed over all independent training runs.
Figure 10: Comparison of the p-FNM, GAF-MLP and projected-force-MLP models on GAF metrics across the different training datasets. The metrics are evaluated on the test dataset. The black error bars indicate the training-time uncertainty quantified by the ±1.96σ confidence interval computed over all independent training runs.
Appendix figures & tables8 assets
Supplementary material from the paper’s appendix.
Appendix
Model
Parameters count
VAE
1,083,649
TPM (RNN module only)
75,104
TPM (full architecture)
1,158,753
p-FNM
1,156,094
GAF-MLP
1,169,427
Projected-force-MLP
1,199,425
Appendix
Table 4: Number of trainable parameters for the deep learning architectures considered in this study. The TPM parameter count is reported both for the recurrent prediction module alone and for the complete model, including the VAE component.
Dataset
Metric
VAE
VAE with SOAP
TPM
TPM-SOAP
Small
Pressure MAE [kPa]
0.97 ± 0.08
0.87 ± 0.03
2.05 ± 0.34
1.80 ± 0.31
Pressure MAPE [%]
1.59 ± 0.14
1.42 ± 0.047
3.35 ± 0.54
2.91 ± 0.48
Pressure MME [kPa]
18.74 ± 1.63
16.98 ± 1.10
28.53 ± 3.58
25.30 ± 3.04
Pressure R2
0.9904 ± 0.0014
0.9914 ± 0.0013
0.9501 ± 0.018
0.9610 ± 0.015
GAF norm MAPEρ [%]
27.29 ± 11.28
25.56 ± 6.59
67.11 ± 32.24
38.52 ± 13.75
GAF norm Rρ2
0.70 ± 0.16
0.71 ± 0.089
0.018 ± 0.43
0.45 ± 0.31
Appendix
Table 5: Comparison of the VAE and TPM models [ 12 ] with their variants trained using the SOAP optimizer. Metrics are reported as the mean ± 1.96 std for the test dataset across 30 training runs.
Dataset
VAE
VAE with SOAP
TPM alone
TPM alone (SOAP variant)
TPM total
TPM-SOAP total
Small
0.16 ± 0.08
0.20 ± 0.08
0.02 ± 0.003
0.03 ± 0.01
0.18 ± 0.08
0.23 ± 0.09
Medium
0.29 ± 0.12
0.32 ± 0.11
0.04 ± 0.01
0.05 ± 0.02
0.33 ± 0.13
0.37 ± 0.13
Large
0.46 ± 0.18
0.64 ± 0.28
0.11 ± 0.03
0.14 ± 0.07
0.57 ± 0.21
0.78 ± 0.35
Extra Large
0.83 ± 0.25
1.38 ± 0.56
0.29 ± 0.03
0.27 ± 0.1
1.12 ± 0.28
1.65 ± 0.66
Appendix
Table 6: Comparison of models training time [GPU.h] for the VAE and TPM models [ 12 ] with their variants trained using the SOAP optimizer. Metrics are reported as the mean ± 1.96 std for the test set across 30 training runs.
Table 8: Detailed metrics for the p-FNM model across the training datasets. Metrics are reported as the mean ± 1.96 std for the test dataset across 30 training runs.
Table 9: Example of predicted and reference pressure fields for a test-dataset simulation at a fixed time step. The red vertical line corresponds to the leading edge of the blade.
Table 10: Example of predicted and reference pressure fields for a test-dataset simulation at a fixed time step. The red vertical line corresponds to the leading edge of the blade.
Table 11: t-SNE projection of latent means and log-variances obtained by encoding all time steps of different test-dataset simulations using a VAE trained on the large dataset.
Institute of Turbomachinery, Xi’an Jiaotong University, Xi’an, 710049, China · ENN Science and Technology Development China Co., Ltd., Langfang, 065001, China · Hebei Key Laboratory of Compact Fusion, Langfang, 065001, China