A grad-div stabilized projection finite element method for a double-diffusive natural convection model