Computational fluid dynamics (CFD) simulation of flow over a three-dimensional (3D) axisymmetric hill presents a unique set of challenges for turbulence modeling. The flow past the crest of the hill is characterized by boundary layer separation, complex vortical structures, and unsteady wake flow. As a result, traditional eddy-viscosity Reynolds-averaged Navier–Stokes (RANS) models have been found to perform poorly for this benchmark test case. Recent studies have focused on the use of large eddy simulation (LES) and hybrid RANS–LES (HRL) methods to improve accuracy. In this study, several different HRL models are investigated, and results from the different models are evaluated relative to each other, to an eddy-viscosity RANS model, and to previously documented high-fidelity LESs and experimental data. Results obtained from the simulations in terms of mean flow statistics, surface pressure distribution, and turbulence characteristics are presented and discussed in detail. Results indicate that eddy-viscosity based HRL models can improve predictions over comparable eddy-viscosity based RANS models such as the shear stress transport (SST) k–ω model used in this study, but only when the development of turbulent velocity fluctuations in the separated shear layer and recirculation region are well resolved.