Modeling the turbulent heat fluxes in low Prandtl number shear flows