The flow around different prolate (needle-like) and oblate (disc-like) spheroids is studied using a multi-relaxation-time lattice Boltzmann method. We compute the mean drag coefficient C D,φ at different incident angles φ for a wide range of Reynolds numbers (Re). We show that the sine-squared drag law2 φ holds up to large Reynolds numbers Re = 2000. Further, we explore the physical origin behind the sine-squared law, and reveal that surprisingly, this does not occur due to linearity of flow fields. Instead, it occurs due to an interesting pattern of pressure distribution contributing to the drag at higher Re for different incident angles. The present results demonstrate that it is possible to perform just two simulations at φ = 0• and φ = 90• for a given Re and obtain particle shape specific C D at arbitrary incident angles. However, the model has limited applicability to flatter oblate spheroids, which do not exhibit the sine-squared interpolation, even for Re = 100, due to stronger wakeinduced drag. Regarding lift coefficients, we find that the equivalent theoretical equation can provide a decent approximation, even at high Re, for prolate spheroids.