Discussion: “Heat Transfer and Wall Heat Flux Partitioning During Subcooled Flow Nucleate Boiling–A Review” (, and , 2006, Journal of Heat Transfer, 128, pp. 1243–1256) PUBLIC ACCESS

[+] Author and Article Information
Satish G. Kandlikar

 Rochester Institute of Technology, Rochester, NYsgkeme@rit.edu

J. Heat Transfer 129(9), 1300-1301 (Mar 26, 2007) (2 pages) doi:10.1115/1.2745830 History: Received February 05, 2007; Revised March 26, 2007

In their paper (1), the authors state on p. 1245 that Boyd and Meng (2) in 1996 suggested an interpolation method for calculating the heat transfer characteristics in the partial boiling region, along with Eqs. (13)-(14). The authors of (1) further state that Kandlikar (3) in 1998 proposed a similar scheme and that the constants a, b, and m were assumed to be constant.

In this discussion, we would like to point out a few errors/omissions in the above referenced paper (1).

1. In their technical note, Boyd and Meng (2) refer to the paper by Kandlikar (4), published in 1990, for the Eqs. (13)-(14) referred to in (1). Unfortunately, Boyd and Meng quoted the wrong reference for the Kandlikar work in which the original model and equations were reported. The 1990 paper by Kandlikar (4), erroneously referred by Boyd and Meng, presents a correlation in the saturated flow boiling region and makes no reference to subcooled flow boiling. The correct reference in Boyd and Meng’s paper should have been Kandlikar (5), which was published in 1991.

2. In his 1991 paper, Kandlikar (5) presented Eqs. (21)-(29) with an accompanying Fig. 4 explaining the construction in the subcooled partial boiling region. Immediately following Eq. (26), a note appears regarding the exponent m, which matches the slope at the two ends. Kandlikar (5), however, contains typographical errors in Eqs. (25) and (26) for coefficients a and b.

3. Those typographical errors were later corrected by Kandlikar (3) in 1998.

4. Another error that appears in (1) is the incorrect year of publication of the paper by Boyd and Meng (2), cited as Ref. [12]; the correct year of publication is 1995, not 1996.

5. Discussion related to Fig. 6 appearing in (1) is correct, except that the 1998 paper by Kandlikar (3) referred therein simply corrects the typographical error and provides a more detailed comparison with the available experimental data.

6. Equation (21) appearing in Warrier and Dhir (1) is based on the paper by Kandlikar (3) as correctly reported by the authors of (1). However, the equation itself is incorrectly reproduced. The correct form of the equation (appearing as Eq. 14 in Ref (3)) is as follows:Display Formula

where ṁ is the mass flux (kgm2s). The negative sign in the exponent 0.7 is missing in Eq. (21) in (1), and the wall superheat ΔTsat is incorrectly replaced by (ΔTw+ΔTsub).

Subcooled Flow Boiling Model Description

For convenience, the correct form of the subcooled flow boiling model and coefficients are presented below.

Figure 1 shows a schematic representation of the subcooled flow boiling curve extending from the single-phase region at point C to the fully developed boiling at point E. In the single-phase region, to the left of C in Fig. 1, heat flux q̇ is given by the following equation:Display Formula

where hLO is the single-phase heat transfer coefficient with all flow as liquid, the local wall superheat is ΔTsat(=TwTsat) and the local liquid subcooling is ΔTsub(=TsatTf), the local saturation temperature is Tsat, the wall temperature is Tw and the liquid temperature is Tf.

In the fully developed boiling region, the heat transfer rate is related to the local wall superheat by the following equation:Display Formula

where Bo is the boiling number (q̇(Ghfg))G is the total mass flux (kgm2s), and hfg is the latent heat of vaporization (J/kg), and FFl is the fluid-dependent parameter in the Kandlikar correlation (3). The value of FFl is 1 for water and all other fluids flowing in stainless steel tubes. For specific fluids in different tube materials, refer to (4) or other more recent publications.

The equation for the q̇ΔTsat plot in the partial boiling region, the main focus of the current discussion, is given by the following equation:Display Formula

The constants a, b, and m are functions of heat flux q̇. The slope of the heat flux versus wall superheat in the partial boiling region is matched with the two limiting values, i.e., m=1 in the single-phase region at the beginning of the partial boiling region, identified by point C, and m=10.3 at the beginning of the fully developed boiling region identified by point E. Thus, the values of a and b are obtained in terms of the heat fluxes and wall superheats at C and E, and the value of m is obtained in terms of the heat fluxes at C, E and at the desired location, where heat flux is q̇ and wall superheat is ΔTsat.Display Formula
andDisplay Formula
Note that there were typographical errors in Kandlikar (5) that erroneously omitted the exponent m in Eqs. 5,6.

The value of m depends on the heat flux, and is allowed to vary linearly from m=1 at C to m=10.3 at D. Thus,Display Formula

and the values of n and p are obtained as follows:Display Formula
andDisplay Formula
Note that n and p are constants for a system (for a given geometry and operating conditions), whereas the values of m, a, and b depend of the local value of q̇. Thus, the value of ΔTsat can be obtained directly from a known value of q̇, while an iterative scheme is needed to calculate q̇ for a given value of ΔTsat in this region.

Further details on calculating q̇C and q̇E are given in Kandlikar (3).

Copyright © 2007 by American Society of Mechanical Engineers
View article in PDF format.



Grahic Jump Location
Figure 1

A representative log-log plot of q̇ versus ΔTsat for describing the Kandlikar model in the partial boiling region, adapted from Kandlikar (5)




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