This study presents a novel approach for enhancing Reynolds-averaged Navier-Stokes (RANS) turbulence modeling through the application of a Relative Importance Term Analysis (RITA) methodology to develop a new zonally-augmented k−ω SST model. Traditional Linear Eddy Viscosity Models often struggle with separated flows. Our approach introduces a physics-based binary classifier that systematically identifies separated shear layers requiring correction by analyzing the relative magnitudes of terms in the turbulent kinetic energy equation. Using symbolic regression, we develop compact correction terms for Reynolds stress anisotropy and turbulent kinetic energy production. Trained on 2D configurations, our model demonstrates significant improvements in predicting separation dynamics while maintaining baseline performance in fully attached flows. Generalization tests on Ahmed body and Faith hill 3D configurations confirm robust transferability, establishing an effective methodology for targeted enhancement of RANS predictions in separated flows.

Data-driven techniques have improved the accuracy of Reynolds-averaged Navier-Stokes (RANS) models in fluid dynamics. However, modeling separated flows remains challenging due to their complex physics and sensitivity to local conditions. Existing approaches often struggle with generalization beyond training cases and lack robust uncertainty quantification frameworks, limiting their utility in complex flow regimes. We propose a Bayesian neural network (BNN)-based framework specifically designed for two-dimensional separated flows. By focusing on flow zones near separated regions, we ensure targeted training and enhance predictive reliability. The BNN framework incorporates physics-guided, invariant inputs to maintain consistency with turbulence physics. Correction terms predicted by the BNN are selectively applied to specific regions of the flow domain using a novel classifier, improving accuracy. A key feature of this approach is propagating BNN-derived corrections to flow solutions, enabling uncertainty quantification in unseen test cases. This probabilistic characterization of modeling errors offers insights into the reliability of RANS predictions across geometries with similar topologies. Preliminary results demonstrate that this method accurately predicts correction terms for Reynolds stress anisotropy and turbulent kinetic energy production in separated flow regions, effectively addressing dominant modeling errors and advancing turbulence modeling through uncertainty quantification.

We present a data-driven approach to Reynolds-averaged Navier-Stokes (RANS) turbulence closure modelling in magnetohydrodynamic (MHD) flows. In these flows the magnetic field interacting with the conductive fluid induces unconventional turbulence states such as quasi two-dimensional (2D) turbulence, and turbulence suppression, which are poorly represented by standard Boussinesq models. Our data-driven approach uses time-averaged Large Eddy Simulation (LES) data of annular pipe flows, at different Hartmann numbers, to derive corrections for the - SST model. Correction fields are obtained by injecting time averaged LES fields into the MHD RANS equations, and examining the remaining residuals. The correction to the Reynolds-stress anisotropy is approximated with a modified Tensor Basis Neural Network (TBNN). We extend the generalised eddy hypothesis with a traceless antisymmetric tensor representation of the Lorentz force to obtain MHD flow features, thus keeping Galilean and frame invariance while including MHD effects in the turbulence model. The resulting data-driven models are shown to reduce errors in the mean flow, and to generalise to annular flow cases with different Hartmann numbers from those of the training cases.