Numerical Study of Vortical Separation From a Three-Dimensional Hill Using Eddy-Viscosity Models

Turbulent flow over an axisymmetric hill is highly three-dimensional (3D) due to the presence of both streamwise and spanwise pressure gradients. Complex vortical separations and reattachments of the turbulent boundary layer are observed on the lee side, accurate prediction of which presents a demanding task for linear eddy-viscosity models (EVMs) when compared to attached boundary layer flows. In this study, an axisymmetric hill is investigated using three Reynolds-averaged Navier-Stokes (RANS) models — fully turbulent model (SST k-ω), transition-sensitive model (k-k_L-ω), and a new four-equation model (k-k_L-ω-v²). The new model is designed to exhibit physically correct responses to flow transition, streamline curvature, and system rotation effects. The test case includes a hill mounted in a channel with hill height H = 2δ, where δ is the approach turbulent boundary layer thickness. The flow Reynolds number (Re) based on the hill height is Re_H = 1.3 × 10⁵. Computational fluid dynamics (CFD) simulation results obtained using the new model are compared with the other two RANS models and with experimental data. Improved mean flow statistics are obtained using the new model that match well with the experiments. The results from this study highlight the need for a model that is able to resolve both flow transition and streamline curvature effects over blunt/curved bodies with reasonable engineering accuracy and computational cost.