A boundary value problem is formulated for a stationary model of mass transfer, which generalizes the Boussinesq approximation in the case when the coefficients in the model equations can depend on the concentration of a substance or on spatial variables. The global existence of a weak solution of this boundary value problem is proved. Some fundamental properties of its solutions are established. In particular, the validity of the maximum principle for the substance’s concentration has been proved. Sufficient conditions on the input data of the boundary value problem under consideration, which ensure the local existence of the strong solution from the space H2, and conditions that ensure the conditional uniqueness of the weak solution with additional property of smoothness for the substance’s concentration are established.