Abstract:
This study presents an extensive numerical investigation on how the addition of polymer molecules into
a Newtonian solvent would tend to influence the natural convection heat transfer phenomena in a differentially heated square cavity. The rheological behaviour of the resulting polymer solution is modelled
based on the FENE-P viscoelastic constitutive equation, which along with other governing equations,
namely, mass, momentum, and energy equations, have been solved using the open source CFD code
OpenFOAM over the following ranges of conditions: Rayleigh number, 103 ≤ Ra ≤ 106; Weissenberg number, 1 ≤ Wi ≤ 100; polymer extensibility parameter, 10 ≤ L2 ≤ 500, polymer viscosity ratio, 0.5 ≤ β ≤ 0.9
and for a fixed value of the Prandtl number of Pr = 7. At low values of the Rayleigh number, the average Nusselt number gradually increases with the Weissenberg number. However, at high values of it, the
average Nusselt number first increases steeply up to a critical value of the Weissenberg number, and beyond that, it remains almost constant as the Weissenberg number further increases. The average Nusselt
number is seen to decrease with the increasing values of both the polymer extensibility paramter and
polymer viscosity ratio. The viscous dissipation tends to deteriorate the heat transfer rate in comparison
to that seen in the absence of it; however, the extent of this deterioration is found to be independent on
the values of Wi, L2 and β. Furthermore, a detailed discussion of the results in terms of the streamline
profiles, isotherm contours, distribution of local Nusselt number, variation of velocity components, etc.,
is also presented. Finally, from an application standpoint, a simple correlation for the average Nusselt
number is presented, which can be used for the interpolation of the present results for the intermediate
values of the governing parameters in a new application.
