Computational fluid dynamics (CFD) methods and schemes have been evolving at a rate that significantly outpaces the equipment needed to readily utilize them at scale. This lack of computational resources has resulted in an increased reliance on turbulence models and the need to know where turbulence models do well, where they do poorly, and where/how they can be improved upon. The BeVERLI Hill experiments aim to address this issue by providing experimental data that achieves a completeness level of three, which has never been done for this type of project. The experimental data collected is studied along side computational results from CFD solvers in order to help address and answer these questions. This paper provides an overview of the current computational status of the BeVERLI Hill project at Virginia Tech. The computational grids used for the analyses are presented such that the reader can gain an appreciation for the modeling techniques and methods being implemented. An analysis of the numerical error associated with the computational results is presented to provide confidence in the results obtained. An in-depth analysis will be presented that shows the results for the various grid levels that are being utilized to determine any grid based effects that are occurring within the solutions. Then, an analysis of the influence of the Reynolds numbers being run is shown. An investigation into the differences between the two different solvers being utilized, SENSEI and Fluent, is shown. An analysis of the effects on the solutions due to numerical limiters is presented to assist in increasing the computational efficiency of the workflow while not adversely affecting the results. Finally, an analysis of the differences between the two turbulence models being utilized is presented. Computational results are compared to available experimentally obtained data to further motivate and identify flow features. / Master of Science / An analysis has been done with high-fidelity computational fluid dynamic solvers that are utilized in order to solve for the flow over a three-dimensional bump called BeVERLI. An analysis is provided that discusses the use of different computational meshes, solvers, turbulence models, and numerical limiters within the computational tools to characterize the flow over the bump. An analysis of the estimated amount of numerical error within the solutions is provided along with a comparison to experimentally obtained data.