Abstract
In this article, a grad-div stabilized projection finite element method is proposed for the double-diffusive natural convection model. Moreover, for the presented method, the stability analysis, and error estimates are deduced. Finally, numerical tests are provided that demonstrate the efficiency of the method. It is found that the presented method can improve the accuracy of the solution compared with the standard projection method, and compute approximations with better mass conservation.