An exactly divergence-free finite element method for non-isothermal flow problems