This paper proposes an analytical model for natural convection in a closed rectangular enclosure filled by a fluid, with imposed heat fluxes at the vertical walls and adiabatic horizontal walls. The analytical model offers a simplified, but easy to handle, description of the temperature and velocity fields. The predicted temperature, velocity, and pressure fields are shown to be in agreement with those obtained from a reliable numerical model. The Nusselt numbers for both the analytical and numerical solutions are then calculated and compared, varying both the aspect ratio of the enclosure and the Rayleigh number. Based on the comparisons, it is possible to assess the dependence of the reliability of the analytical model on the aspect ratio of the enclosure, showing that the prediction error rapidly decreases with the increase of the enclosure slenderness.