Abstract Natural convection processes are pivotal in various engineering applications, necessitating accurate and reliable computational fluid dynamics (CFD) models for their simulation. This study evaluates the efficacy of different turbulence models implemented in OpenFOAM for a natural convection scenario, aiming to identify the most suitable model for capturing complex thermal and fluid dynamic behaviors. We compared several turbulence models, including the k − ε, k − ω, and k − ω SST, using a benchmark thermal convection case. Our methodology involved setting up the simulations to reflect realistic thermal gradients and boundary conditions, followed by a rigorous analysis of temperature distribution, vertical velocity profiles, and computational efficiency. The results indicate that the k − ω SST achieved the lowest averaged root mean square error (RMSE) values for temperature (0.076) and vertical velocity (0.03) while also requiring the fewest convergence iterations (314) compared to k − ε (411) and k − ω (1474). These findings demonstrate that the k − ω SST model is a suitable compromise between accuracy and computational cost for engineering applications. This study underscores the importance of selecting an appropriate turbulence model in OpenFOAM to enhance the accuracy of natural convection simulations, which can significantly influence design and safety considerations in engineering systems.