Comparison of two low Mach number algorithms for thermobuoyant flows.
Both algorithms perform similarly in the Boussinesq limit.
Discrete mass conservation errors lead to numerical instabilities for low Prandtl non-Boussinesq convection.
Discrete mass conservation errors affect flow physics for unsteady non-Boussinesq convection.