In this paper, we consider a natural convection flow of an incompressible viscous fluid subject to Newtonian heating and constant mass diffusion. The proposed model has been described by the Caputo fractional operator. The used derivative is compatible with physical initial and boundaries conditions. The exact analytical solutions of the proposed model have been provided using the Laplace transform method. The obtained solutions are expressed using some special functions as the Gaussian error function, Mittag–Leffler function, Wright function, and
G
-function. The influences of the order of the fractional operator, parameters used in modeling the considered fluid, Nusselt number, and Sherwood number have been analyzed and discussed. The physical interpretations of the influences of the parameters of our fluid model have been presented and analyzed as well. We use the graphical representations of the exact solutions of the model to support the findings of the paper.