A divergence-free low-order stabilized finite element method for a generalized steady state Boussinesq problem