Numerical investigation of buoyancy effect on heat transfer to carbon dioxide flow in a tube at supercritical pressures

• Published in: International journal of heat and mass transfer 2018-02, Vol.117, p.595-606
• Main Authors: Liu, Sheng-hui, Huang, Yan-ping, Wang, Jun-feng, Leung, Laurence K.H.
• Format: Article
• Language: English
• Subjects: Buoyancy; Buoyancy effect; Carbon dioxide; Computational fluid dynamics; Critical point; Data correlation; Deteriorated and post-deteriorated heat transfer; Diffusion layers; Eddy viscosity; Flow and turbulence fields; Fluid dynamics; Fluid flow; Heat recovery; Heat transfer; Numerical analysis; Numerical simulation; Software; Supercritical CO2 heat transfer; Supercritical pressures; Turbulence models; Turbulent flow; Velocity; Velocity gradient; Viscous sublayers
• ISSN: 0017-9310; 1879-2189
• DOI: 10.1016/j.ijheatmasstransfer.2017.10.037

•Deteriorated and post-deteriorated heat transfer was studied numerically.•The mechanism of deteriorated heat transfer was analyzed.•Enhancement heat transfer was found at the post-deteriorated heat transfer section.•Two special characteristics of radial velocity gradient near the wall were summarized.•Enhancement heat transfer at post-deteriorated heat transfer section was explained. A numerical analysis of the buoyancy effect on heat transfer to carbon dioxide flow at supercritical pressures in a vertical heated tube has been carried out using the FLUENT computational fluid dynamics (CFD) tool with the shear-stress-transport low Reynolds number eddy viscosity turbulence model. It covered the 1200-mm heated section and included 500-mm unheated section upstream for flow development. Wall temperatures calculated from the analysis were compared against those obtained from the experiments. Although the analysis overestimated the experimental wall temperatures, it captured appropriately the general trends of the experimental results. Based on the examinations of the calculated flow and turbulence fields, the effect of buoyancy on heat transfer has been quantified. Results showed that there are two obvious characteristics of the velocity gradient profile influenced by severe buoyancy effect, i.e., (1) zero-velocity-gradient region may come into being; (2) the slope of the velocity gradient profile near the wall may become very large. When the zero-velocity-gradient region comes into the rim of viscous sub-layer, turbulent production and diffusion will be damped, finally deteriorated heat transfer occurs. Enhancement heat transfer may happen at the post-deteriorated heat transfer section because of the large slope of the mean velocity gradient profile, and the possible mechanism was introduced. A striking phenomenon has been found, i.e., enhancement heat transfer might also happen when the buoyancy parameter is lower the value corresponding to heat transfer recovery. This phenomenon relates to the difficulty of heat transfer data correlation near the pseudo-critical point.