In this paper, we consider the stationary double-diffusive natural convection model, which can model heat and mass transfer phenomena. Based on the fixed point theorem, the existence and uniqueness of the considered model are proved. Moreover, we design three finite element iterative methods for the considered problem. Under the uniqueness condition of a weak solution, iterative method I is stable. Compared with iterative method I, iterative method II is stable with a stronger condition. Moreover, iterative method III is stable with the strongest condition. From the perspective of viscosity, iterative method I displays well in the case of a low viscosity number, iterative method II runs well with slightly low viscosity, and iterative method III can deal with high viscosity. Finally, some numerical experiments are presented for testing the correctness of the theoretic analysis.