We consider the new boundary value problem for the generalized Boussinesq model of heat transfer under the inhomogeneous Dirichlet boundary condition for the velocity and under mixed boundary conditions for the temperature. It is assumed that the viscosity, thermal conductivity and buoyancy force in the model equations, as well as the heat exchange boundary coefficient, depend on the temperature. The mathematical apparatus for studying the inhomogeneous boundary value problem under study based on the variational method is being developed. Using this apparatus, we prove the main theorem on the global existence of a weak solution of the mentioned boundary value problem and establish sufficient conditions for the problem data ensuring the local uniqueness of the weak solution that has the additional property of smoothness with respect to temperature. This article is part of the theme issue 'Non-smooth variational problems with applications in mechanics'.