The primary goal of the work presented in this thesis is to assess the effectiveness of Computational Fluid Dynamics (CFD) tools in predicting flow fields around golf discs. The quantities of interest involve aerodynamic force and moment coefficients, surface pressure distributions, and velocity and pressure distributions in the near field region. The veracity of the CFD predictions are investigated by comparing the CFD results against wind tunnel experiments. To date, direct comparisons between CFD and wind tunnel data for commercially available golf discs have not been published, and this will be the first work to do so. In doing so, this thesis also outlines the best practices for the CFD analysis of golf discs. The methods detailed in this thesis can be used to evaluate future disc designs. Three new disc design concepts are conceived and evaluated with the method. Steady-state Reynolds Average Navier-Stokes (RANS) simulations on highly resolved grids are performed using the k-ω Shear Stress Transport (SST) and Lag elliptic blending k-ε turbulence models. The latter one is a newer model and existing literature shows a very limited number of studies carried out using this turbulence model. The simulations presented in this work were carried out using a commercial finite volume CFD code, STAR-CCM+. It was observed that, compared to the SST model, the Lag elliptic blending k-ε turbulence model produces better correlation with experimental data. However, further improvement of experimental correlation requires that the turbulence model closure coefficients used in the Lag elliptic blending k-ε model be tuned to better correlate the RANS simulations to the large database of experimental data for a particular disc. Additionally, in order to understand the limitations of the steady-state solution of an inherently unsteady phenomenon, transient Detached Eddy Simulations (DES) are also performed and the results are compared to the steady-state RANS and experimental data.