Research Papers: Natural and Mixed Convection

Natural Convection From Horizontal Cylinders at Near-Critical Pressures—Part II: Numerical Simulations

[+] Author and Article Information
Gopinath R. Warrier

e-mail: gwarrier@ucla.edu

Vijay K. Dhir

Henry Samueli School of
Engineering and Applied Science,
Mechanical and Aerospace
Engineering Department,
University of California, Los Angeles,
Los Angeles, CA 90095

1Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the Journal of Heat Transfer. Manuscript received September 9, 2011; final manuscript received June 17, 2012; published online December 28, 2012. Assoc. Editor: Ali Ebadian.

J. Heat Transfer 135(2), 022502 (Dec 28, 2012) (10 pages) Paper No: HT-11-1438; doi: 10.1115/1.4007673 History: Received September 09, 2011; Revised June 17, 2012

A numerical investigation of laminar natural convection heat transfer from small horizontal cylinders at near-critical pressures has been carried out. Carbon dioxide is the test fluid. The parameters varied are: pressure (P), (ii) bulk fluid temperature (Tb), (iii) wall temperature (Tw), and (iv) wire diameter (D). The results of the numerical simulations agree reasonably well with available experimental data. The results obtained are as follows: (i) At both subcritical and supercritical pressures, h is strongly dependent on Tb and Tw. (ii) For Tw < Tsat (for P < Pc) and Tw < Tpc (for P > Pc), the behavior of h as a function of Tw is similar; h increases with increase in Tw. (iii) For P > Pc and large Tw (Tw > Tpc), natural convection heat transfer occurring on the cylinder is similar that observed during film boiling on a cylinder. The heat transfer coefficient decreases as Tw increases. (iv) For subcritical pressures, the dependence of h on D is h ∝ D0.5 in the range 25.4 ≤ D ≤ 100 μm. For larger values of D (500–5000 μm), h ∝ D0.24. (v) For supercritical pressures, the dependence of h on D is h ∝ D0.47 in the range 25.4 ≤ D ≤ 100 μm. For larger values of D (500–5000 μm), h ∝ D−0.27. (vi) For a given P, the maximum heat transfer coefficient is obtained for conditions where Tb < Tpc and Tw ≥ Tpc. Analysis of the temperature and flow field shows that this peak in h occurs when k, Cp, and Pr in the fluid peak close to the heated surface.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Hahne, E., 1985, “Natural Convection in the Near Critical Region and Its Application in Heat Pipes,” Natural Convection: Fundamentals and Applications, S. Kakac, W. Aung, and R. Viskanta, eds., Springer-Verlag, Berlin, pp. 774–826.
Simon, H. A., and Eckert, E. R. G., 1963, “Laminar Free Convection in Carbon Dioxide Near Its Critical Point,” Int. J. Heat Mass Transfer, 6, pp. 681–690. [CrossRef]
Knapp, K. K., and Sabersky, R. H., 1966, “Free Convection Heat Transfer to Carbon Dioxide Near the Critical Point,” Int. J. Heat Mass Transfer, 9, pp. 41–45. [CrossRef]
Abadzic, E., and Goldstein, R. J., 1970, “Film Boiling and Free Convection Heat Transfer to Carbon Dioxide Near the Critical State,” Int. J. Heat Mass Transfer, 13, pp. 1163–1175. [CrossRef]
Neumann, R. J., and Hahne, E. W. P., 1980, “Free Convective Heat Transfer to Supercritical Carbon Dioxide,” Int. J. Heat Mass Transfer, 23, pp. 1643–1652. [CrossRef]
Nishikawa, K., and Ito, T., 1969, “An Analysis of Free Convective Heat Transfer From an Isothermal Vertical Plate to Supercritical Fluids,” Int. J. Heat Mass Transfer, 12, pp. 1449–1463. [CrossRef]
Hilal, M. A., 1978, “Analytical Study of Laminar Free Convection Heat Transfer to Supercritical Helium,” Cryogenics, 18, pp. 545–551. [CrossRef]
Seetharam, T. R., and Sharma, G. K., 1979, “Free Convective Heat Transfer to Fluids in the Near-Critical Region From Vertical Surfaces With Uniform Heat Flux,” Int. J. Heat Mass Transfer, 22, pp. 13–20. [CrossRef]
Rousselet, Y., Warrier, G. R., and Dhir, V. K., 2013, “Natural Convection From Horizontal Cylinders at Near-Critical Pressures—Part I: Experimental Study,” ASME J. Heat Transfer, 135(2), p. 022501. [CrossRef]
Paolucci, S., 1982, “On the Filtering of Sound From the Navier-Stokes Equations,” Scandia National Laboratories, Technical Report No. SAND82-8257.
Zappoli, B., Amiroudine, S., Carlès, P., and Ouazzani, J., 1996, “Thermoacoustic and Buoyancy-Driven Transport in a Square Side-Heated Cavity Filled With a Near-Critical Fluid,” J. Fluid Mech., 316, pp. 53–72. [CrossRef]
Muller, B., 1998, “Low Mach Number Asymptotics of the Navier-Stokes Equations,” J. Eng. Math., 34, pp. 97–109. [CrossRef]
Le Quéré, P., Weisman, C., Paillère, H., Vierendeels, J., Dick, E., Becker, R., Braack, M., and Locke, J., 2005, “Modelling of Natural Convection Flows With Large Temperature Differences: A Benchmark Problem for Low Mach Number Solvers. Part 1. Reference Solutions,” Math. Model. Numer. Anal., 39, pp. 609–616. [CrossRef]
Accary, G., and Raspo, I., 2006, “A 3D Finite Volume Method for the Prediction of a Supercritical Fluid Buoyant Flow in a Differentially Heated Cavity,” Comput. Fluids, 35, pp. 1316–1331. [CrossRef]
Mazumder, S., 2007, “On the Use of the Fully Compressible Navier-Stokes Equations for the Steady-State Solution of Natural Convection Problems in Closed Cavities,” ASME J. Heat Transfer, 129, pp. 387–390. [CrossRef]
Hauke, G., and Hughes, T. J. R., 1998, “A Comparative Study of Different Sets of Variables for Solving Compressible and Incompressible Flows,” Comput. Methods Appl. Mech. Eng., 153, pp. 1–44. [CrossRef]
NIST Standard Reference Database 23, REFPROP, Reference Fluid Thermodynamic and Transport Properties, Version 8.0.
Ferziger, J. H., and Peric, M., 1999, Computational Methods for Fluid Dynamics, Springer-Verlag, Berlin, Germany.
Khosla, P. K., and Rubin, S. G., 1974, “A Diagonally Dominant Second-Order Accurate Implicit Scheme,” Comput. Fluids, 2, pp. 207–209. [CrossRef]
Waterson, N. P., and Deconinck, H., 2007, “Design Principles for Bounded Higher-Order Convection Schemes – A Unified Approach,” J. Comput. Phys., 224, pp. 182–207. [CrossRef]
Van Doormaal, J. P., and Raithby, G. D., 1984, “Enhancements of the SIMPLE Method for Predicting Incompressible Fluid Flows,” Numer. Heat Transfer, 7, pp. 147–163. [CrossRef]
Stone, H. L., 1968, “Iterative Solution of Implicit Approximations of Multidimensional Partial Differential Equations,” SIAM Journal of Numerical Analysis, 5, pp. 530–558. [CrossRef]
Saad, Y., and Schultz, M. H., 1986, “GMRES: A Generalized Residual Algorithm for Solving Non-Symmetric Linear Systems,” SIAM J. Sci. Stat. Comput., 7, pp. 856–869. [CrossRef]
Kuehn, T. H., and Goldstein, R. J., 1980, “Numerical Solution to the Navier-Stokes Equations for Laminar Natural Convection About a Horizontal Isothermal Circular Cylinder,” Int. J. Heat Mass Transfer, 23, pp. 971–979. [CrossRef]
Saitoh, T., Sajiki, T., and Maruhara, K., 1993, “Bench Mark Solutions to Natural Convection Heat Transfer Problem Around a Horizontal Circular Cylinder,” Int. J. Heat Mass Transfer, 36, pp. 1251–1259. [CrossRef]
Churchill, S. W., and Chu, H. H. S., 1975, “Correlating Equations for Laminar and Turbulent Free Convection From a Horizontal Cylinder,” Int. J. Heat Mass Transfer, 18, pp. 1049–1053. [CrossRef]
Mills, A. F., 1999, Heat Transfer, 2nd ed., Prentice-Hall, Englewood Cliffs, NJ.
Kato, H., Nishiwaki, N., and Hirata, M., 1968, “Studies on the Heat Transfer of Fluids at a Supercritical Pressure,” Bull. Jpn. Soc. Mech. Eng., 11, pp. 654–663. [CrossRef]


Grahic Jump Location
Fig. 1

Property variation as a function of temperature for CO2 [8] (a) subcritical pressure (P = 6.99 MPa) and (b) supercritical pressure (P = 8.10 MPa)

Grahic Jump Location
Fig. 2

Computational domain

Grahic Jump Location
Fig. 5

Effect of Tb on h for subcritical pressures (P = 7.34 MPa, D = 100 μm)

Grahic Jump Location
Fig. 3

Effect of number of grid points (a) varying Nη and (b) varying Nθ

Grahic Jump Location
Fig. 4

Comparison of numerical results with experimental data for subcritical pressures, P = 7.34 MPa, 10 °C ≤ Tb ≤ 35 °C, D = 100 μm

Grahic Jump Location
Fig. 6

Effect of pressure for D = 25.4 μm, Tb = 10 °C (P < Pc)

Grahic Jump Location
Fig. 7

Variation of h with D for subcritical pressures

Grahic Jump Location
Fig. 10

Effect of pressure for D = 100 μm, Tb = 10, and 31 °C (P > Pc)

Grahic Jump Location
Fig. 11

Effect of diameter for supercritical pressures, P = 7.50 MPa, Tb = 25, 31, and 35 °C

Grahic Jump Location
Fig. 13

Radial profiles of T, ∂T/∂y and relevant fluid properties

Grahic Jump Location
Fig. 8

Comparison of numerical results with experimental data for supercritical pressures, P = 8.10 MPa, 10 °C ≤ Tb ≤ 33.3 °C, and D = 100 μm

Grahic Jump Location
Fig. 9

Effect of bulk liquid temperature for supercritical pressures

Grahic Jump Location
Fig. 12

Variation of Tw,peak as a function of Tb and P (P > Pc)



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In