Abstract
The present article discusses the numerical simulation results of laminar mixed convection flow in vertical tubes. A comparative analysis of the thermal and hydrodynamic features of both buoyancy-assisted and opposed flows was performed for Reynolds number (), Grashof number (
) and Richardson number (
) with uniform heat flux boundary condition. 2-D axisymmetric steady state simulations were carried out for a length-to-diameter (
) ratio of
with water as the working fluid. Numerical simulations were performed by employing SIMPLE scheme for pressure–velocity coupling in momentum equations and second order UPWIND scheme for solving energy equations. In case of assisting flow (
), at fully developed state the centerline velocity decreases, and velocity is increased near the tube wall due to heat flux induced free convection. With increasing heat flux, the decrement in centerline velocity compared to the no-heat flux condition increases. Further, with increasing heat flux, the increase in friction factor and Nusselt number was observed. Therefore, at the same
the variation of heat flux led to unique velocity profiles and temperature gradients. The variation of centerline velocity and temperature in developing region was also studied. While centerline temperature was monotonically increasing with length, the centerline velocity increased up to a maximum in the developing region and then attained the steady state at a lower value in the fully developed state. Subsequently we studied the dependence of
on the hydrodynamic and thermal features. For constant
friction factor and Nusselt number was observed to increase with increase in
from 0.1 to 1.5 range. At fixed
variation of heat flux leads to variation of
or
Hence, the buoyancy effect has a significant role in the entrance length development. The hydrodynamic entry length after an initial increase, attained an almost constant value with increasing
On the other hand, the thermal entry length exhibited a decreasing trend with increasing
Similarly, at constant
Nusselt number increased with increasing
for a range of 100–1000. It was evident that for a given
and
the heat transfer in the developing flow is always higher as that of the fully developed flow. Contrasting observations were observed for buoyancy-opposed flow. Velocity is accelerated at the center as compared to the tube wall in case of opposing flow for same
and
For a fixed
the friction factor and Nusselt number decreased with increasing the
Both the hydrodynamic and thermal entry length increases for opposing flows. Finally, we have developed correlations for fully developed friction factor (
) with
as well as
and Nusselt number (
) with
for buoyancy-assisting and buoyancy-opposing flows. Two independent correlations are also produced for
with Graetz number (
) for developing and fully developed regimes in both assisting as well as opposing flows.
Acknowledgments
SG would like to thank Ministry of Education, Govt. of India. SKD and DS would like to thank IIT Ropar for providing the computational resources.