Computational fluid dynamics simulation of flow over a three-dimensional 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 large-eddy simulations 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 HRL models can significantly improve predictions over RANS models, but only when the development of turbulent velocity fluctuations in the separated shear layer and recirculation region are well resolved.