Motivated by improving the performance of particle-based Monte-Carlo simulations in the transitional regime, Fokker–Planck kinetic models have been devised and studied as approximations of the Boltzmann collision operator. By generalizing the linear drift model, the cubic Fokker–Planck (cubic-FP) and ellipsoidal Fokker–Planck (ES-FP) have been proposed, in order to obtain the correct Prandtl number of 2/3 for a dilute monatomic gas. This study provides a close comparison between both models in low Mach and supersonic settings. While direct simulation Monte-Carlo (DSMC) here serves as the benchmark, overall close performance between cubic-FP, ES-FP, and DSMC is observed. Furthermore, while the ES-FP outperforms the cubic-FP model in the shock region of the supersonic flow around a cylinder, the latter shows a better accuracy in the near continuum regime. It is argued that the reason behind these discrepancies lies in the entropy law besides the transport properties.