Numerical results for two-dimensional axisymmetic flow of an incompressible fluid over ellipsoidal particles are presented for Reynolds numbers from 0.1 to 200. An bp-finite element method is employed to obtain values of drag around the ellipsoids, which range from spheres to discs. The results are compared to previous numerical analytical, and experimental work presented in the literature. Using this data, a new drag model is created that predicts the drag on oblate ellipsoids. The new model is accurate to within 1.5% over the desired range of parameters.