On the Numerical Prediction of Transitional Flows With Reynolds-Averaged Navier-Stokes and Scale-Resolving Simulation Models

The present work investigates the transitional flow around a smooth circular cylinder at Reynolds number Re = 140,000. The flow is resolved using the viscous-flow solver ReFRESCO, and distinct mathematical models are applied to assess their ability to handle transitional flows. The selected mathematical models are the Reynolds-Averaged Navier-Stokes equations (RANS), Scale-Adaptive Simulation (SAS), Delayed Detached-Eddy Simulation (DDES), eXtra Large-Eddy Simulation (XLES) and Partially-Averaged Navier-Stokes (PANS) equations. The RANS equations are supplemented with the k–ω Shear-Stress Transport (SST) with and without the Local Correlation Transition Model (LCTM). The numerical simulations are carried out using structured grids ranging from 9.32 × 10⁴ to 2.24 × 10⁷ cells, and a dimensionless time-step of 1.50 × 10⁻³. As expected, the outcome demonstrates that transition from laminar to turbulent regime is incorrectly predicted by the k–ω SST model. Transition occurs upstream of the flow separation, which is typical of the supercritical regime and so the flow physics is incorrectly modelled. Naturally, all Scale-Resolving Simulation (SRS) models that rely on RANS to solve the boundary-layer, called hybrid models, will exhibit a similar trend. On the other hand, mathematical models capable to resolve part of the turbulence field in the boundary layer (PANS) lead to a better agreement with the experimental data. Furthermore, the k–ω SST LCTM is also able to improve the modelling accuracy when compared to the k–ω SST. Therefore, it might be a valuable engineering tool if its computational demands are considered (in the RANS context). Therefore, the results confirm that the choice of the most appropriate mathematical model for the simulation of turbulent flows is not straightforward and it may depend on the details of the flow physics.