## Services on Demand

## Journal

## Article

## Indicators

- Cited by SciELO
- Access statistics

## Related links

- Cited by Google
- Similars in SciELO
- Similars in Google

## Share

## Journal of the Brazilian Society of Mechanical Sciences and Engineering

##
*Print version* ISSN 1678-5878*On-line version* ISSN 1806-3691

### J. Braz. Soc. Mech. Sci. & Eng. vol.26 no.2 Rio de Janeiro Apr./June 2004

#### http://dx.doi.org/10.1590/S1678-58782004000200015

**TECHNICAL PAPERS**

**Evaluation of the rohsenow correlation through experimental pool boiling of halocarbon refrigerants on cylindrical surfaces**

**J. M. Saiz Jabardo; E. Fockink da Silva; G. Ribatski; S. F. de Barros**

Laboratório de Refrigeração Departamento de Engenharia Mecânica Escola de Engenharia de São Carlos Universidade de São Paulo Avenida Trabalhador Sãocarlense, 400 Centro 13566-590 São Carlos, SP. Brazil mjabardo@sc.usp.br

**ABSTRACT**

This paper reports results from an investigation on the performance of the Rohsenow's type of correlation when applied to the nucleate boiling of halocarbon refrigerants over cylindrical surfaces of different material. Experimental data for refrigerants R-11, R-123, R-12, and R-134a have been raised and fitted according to two different procedures. It has been determined that exponents m and n are weakly affected by the refrigerant, and surface material and finishing. Liquid/surface combinations and surface roughness affect the coefficient C_{sf}. A correlation has been developed for C_{sf} in terms of the liquid/surface combination and surface roughness. The resulting correlation has been evaluated through the experimental data used in the fitting process and results obtained elsewhere. Deviations of correlation with respect to the experimental heat transfer coefficient are within acceptable ranges.

**Keywords:** Nucleate boiling, halocarbon, refrigerant, cylinder

**Introduction**

Nucleate boiling has been under intense scrutiny since the early thirties, when pioneering research studies came to light. Ribatski (2002) rightly points out that the number of publications has grown steadily since that time. This interest is justified by the heat transfer enhancement that this mechanism promotes on heated surfaces. Generally speaking, nucleate boiling research has been traditionally focused into the investigation of the physical mechanism and the development of general correlations for the heat transfer coefficient. Physical models tend to reproduce the particular mechanism that prevails under specific operational conditions of the boiling surface such as pressure, liquid subcooling, liquid/surface combination, and geometry. As a result, most of the semi-empirical correlations are strongly dependent upon operational conditions not being amenable to generalizations. The Rohsenow's correlation is one of them, having endure up to present days. Based upon straightforward arguments, Rohsenow (1952) suggested that the heat transfer enhancement under boiling conditions is the result of local liquid circulation in the region close to the heating surface promoted by successive bubble detachments. Given this physical scenario, a locally defined Stanton number could be written in terms of an equally local Reynolds number and the liquid Prandtl number. The local Stanton number was referred to the mass velocity of vapor leaving the surface whereas the Reynolds number was written in terms of the same mass velocity and a characteristic dimension proportional to the bubble detachment diameter. After some manipulation the following equation results:

which can be written in a reduced form as:

The left hand side of Eq. (1) corresponds to the inverse of Stanton number, *St**, whereas the first term of the right hand side could be considered a Reynolds number referred to the bubble departure diameter. Values of exponents *m* and *n* were found to be equal to 0.33 and 1.7 by curve fitting experimental results obtained elsewhere for liquids other than water. Rohsenow suggested to change the exponent of the Prandtl number to 1.0 for water. Furthermore, according to Rohsenow, whereas *m* is not affected, *n* might depend upon the surface finishing, attaining values in the range between 0.8 and 2.0. The numerical coefficient, *C _{sf}*, was related to the effect of liquid/surface combination through the contact angle, q . Thus for water boiling on a smooth copper surface, Rohsenow determined the value of

*C*to be equal to 0.013. Later, Vachon et al. (1968) performed extensive experiments involving several liquid/surface combinations to obtain the associated values of

_{sf}*C*and

_{sf}*n*. Neither Rohsenow nor Vachon and coworkers performed experiments involving halocarbon refrigerants boiling on cylindrical surfaces which are the main objective of present study. Noteworthy are studies by Sauer et al. (1975), Saiz Jabardo and Silva (1991) and more recently Pioro (1999) who fitted data from experiments involving a limited number of halocarbon refrigerants. Their data shed some light over the behavior of the Rohsenow correlation with respect to the aforementioned refrigerants. Table 1 presents a summary of studies carried out in relation with the Rohsenow correlation. The table includes information regarding the fluid(s), surface material, operational conditions such as the ranges of pool pressures and heat fluxes, and the values of

*C*,

_{sf}*m*, and

*n*each study has come up with by fitting experimental results. Interesting to note in this table is that the few studies that have dealt with halocarbon refrigerants did it with the low pressure ones, such as R-11 and R-113, both of them of the banned family of CFCs.

The present study has been set forth to experimentally investigate nucleate boiling on cylindrical copper, brass and stainless steel surfaces of both low pressure R-11 (CFC) and R-123 (HCFC), and medium pressure halocarbon refrigerants, R-12 (CFC), and R-134a (HFC). Constant parameters of the Rohsenow correlation are evaluated and analyzed based upon data involving this set of refrigerants and over an extensive range of reduced pressures and heat fluxes along with different surface conditions.

**Nomenclature**

C= surface/liquid parameter of the Rohsenow correlation_{sf}

c= specific heat of the liquid [J/Kg.K]_{pl}

g= gravitational acceleration [m/s^{2}]

h= heat transfer coefficient [W/m^{2}.K]

h= correlation heat transfer coefficient [W/m_{corr}^{2}.K]

h= experimental heat transfer coefficient [W/m_{exp.}^{2}.K]

h= latent heat of evaporation [J/kg]_{lg}

k= thermal conductivity of liquid [W/m.K]_{l}

m,n= exponents of Rohsenow correlation

p= reduced pressure_{r}

Pr= Prandtl number of the liquid_{l}

q"= specific heat flux [W/m^{2}]

Ra= roughness arithmetic average [µm]

Re= Reynolds number referred to the bubble diameter_{b}

Rq= roughness root mean square [µm]

St*= inverse of the Stanton number

T= saturation temperature [K]_{sat}

T= wall temperature [K]_{w}

Greek Symbols

D T= surface superheat (T) [K]_{w}-T_{sat}

m= dynamic viscosity of the liquid [Pa.s]_{l}

q= contact angle [rad]

r= vapor density [kg/m_{g}^{3}]

r= liquid density [kg/m_{l}^{3}]

s= surface tension [N/m]

**Experimental Set Up**

The experimental set up comprises the refrigerant and cooling circuits, as shown in Fig. 1. The charge of refrigerant is basically contained in the boiler in which the liquid is kept at a reasonable level above the test surface (tube) so that the column head does not affect significantly the equilibrium saturation temperature. The cooling circuit is intended to control the equilibrium pressure in the boiler by condensing the refrigerant boiled in the heating surface. The condensing effect is obtained by a 60% solution of ethylene glycol/water that operates as intermediate fluid between the condenser and the cooling system not shown in Fig. 1. The ethylene glycol/water solution is cooled by either a refrigeration circuit or water from a cooling tower, depending upon the operating pressure. This solution is intended to operate in the range between -26ºC and 90ºC.

The boiler is a 40 liters carbon steel container with two lateral circular windows for visualization. It contains the boiling surface in addition to a 1500W/220V electrical heater, installed at the bottom, and two sheathed type T thermocouples. The boiler is also fitted with openings for connections to a pressure transducer, a safety valve, and vapor and liquid return copper lines, as shown in Fig. 1. The sheathed thermocouples are installed in such a way to measure and monitor the temperature of the liquid pool and the vapor in equilibrium with it. Under normal operating conditions these thermocouples indicate temperatures which are very close to each other and to the saturation temperature at the boiler internal pressure measured by the pressure transducer.

The test (boiling) surface is placed in the middle of the boiler so that the boiling mechanism can easily be visualized through the glass windows. It is made up of a 19.0 mm diameter and 3.1 mm thick copper tube, a cut way view of it is shown in Fig 2. The test tube is supported by a brass piece which is thread attached to the flanged cover of the boiler. The boiling surface is heated by a 12.6 mm diameter and 210 mm long cartridge electrical heater. The electrical power to the boiling surface is controlled by a manually operated voltage converter and measured by a power transducer. Surface temperature is measured through eight 30 AWG type T thermocouples installed in grooves carved by an electro-erosion process in locations indicated in Fig. 2. Thermocouples are kept in place by a thermally conductive epoxy resin. Electrical signals from the transducers are processed by a data acquisition system which includes two 12 bit A/D converter boards with 16 channels each, and three connection panels. Two of these panels are dedicated to thermocouple connections.

**Experimental Procedure**

The boiling surface used to be treated prior to the beginning of the tests. Sandpaper scales with mesh size varying in the range from 220 to 1,200 were used to obtain the final surface roughness and applied through a regular lathe machine run at 1,200 rpm. Experiments were also conducted with a polished surface, which required a special treatment, and a sand blast surface. After treatment, the boiling surface used to be thoroughly cleaned with a solvent (normally refrigerant R-11) and the roughness measured at 10 randomly selected regions before attaching it to the boiler. After testing the boiling surface, 10 randomly selected regions were again taken for roughness measurement so that conditions of the surface before and after the tests could be compared. The roughness was measured in terms of the CLA arithmetic average, *Ra.* The treatment suggested above allowed experiments to be run in the range of surface roughness, *Ra*, between 0.02± 0.01 µm and 3.3± 0.4 µm.

The internal surface of the boiler used to be cleaned and kept under a vacuum of less than 2 kPa during a period of 12 hours before the attachment of the boiling surface and the introduction of the test refrigerant. Tests were conducted under saturated conditions of the refrigerant. This condition was continuously monitored and adjusted as needed. The datum point would only be logged if the readings of the sheathed thermocouples were close enough (within 0.2 K) to each other and to the saturation temperature inside the boiler obtained from the pressure transducer reading. For analysis purposes, the saturation temperature of the pool was determined as the average of the readings of the sheathed thermocouples. Tests were conducted by gradually increasing the heat flux up to its predicted maximum. Once the maximum was attained, the heat flux was gradually reduced down to zero. Only downward heat flux data were considered for analysis purposes. Several procedures were tried to check for possible effects on the results. Two of such procedures consisted in keeping the boiling surface active for some time before logging data and starting directly from the maximum heat flux.

In measuring the surface temperature care was exercised in evaluating the thermal resistance of the copper wall between the couple location and the actual boiling surface. In addition, axially located thermocouples helped in evaluating axial heat conduction. It has been determined that in the location corresponding to section 2 of the test tube, Fig. 2, the axial heat flux was negligibly small. A thorough discussion of surface temperature measurement can be found in Ribatski (2002). The temperature considered for analysis purposes was the one from the thermocouple located midway between those at top and bottom of the heating surface at section 2 (Fig. 2). The temperature indicated by this thermocouple is equivalent to the average of the readings of the three section 2 thermocouples.

Instruments were calibrated and the uncertainty of measured parameters evaluated according to the procedure suggested by Abernethy and Thompson (1973) with results summarized in Table 2.

**A Summary of Experimental Results**

Tests conducted under the investigation reported herein have raised information regarding the effects on boiling heat transfer caused by the change in refrigerant, pressure and surface material and finishing condition. Though the main objective of the paper is to correlate results according to the form proposed by Rohsenow, a summarized analysis of the obtained data will be presented in such a way to clear up some important effects related the main physical parameters investigated in present research. These effects can clearly be seen in the boiling curves of Figures 3 to 5. Close examination of these figures allows one to draw the following general conclusions:

- The slopes of boiling curves vary from values close to zero at the incipience of boiling up to an asymptotic level in the high heat flux range. The effect of pressure on the asymptotic slope seems to be of limited extent, especially in the high pressure range (Fig. 3).
- The region of the boiling curve where significant variations of the slope occur corresponds to the one known as Partial Nucleate Boiling, Collier and Thome (1994). The extent of this region depends upon the boiling liquid, pressure, and surface roughness. For conditions considered in present study, heat fluxes corresponding to the Partial Nucleate Boiling region are limited to a maximum of 30 kW/m
^{2}. - Pressure affects positively the heat transfer, and so does roughness as the displacement of the boiling curves to the left suggests. This should be an expected result considering effects of both the surface tension and the range of cavity size, which are respectively associated to boiling incipience and density of active cavities. Refrigerants are affected differently by the pressure and especially by the surface roughness. In fact, as shown in Fig. 4 (a), the displacement of the boiling curve to the left with
*Ra*is more significant for R-11 than for R-123. For the case of medium pressure refrigerants, Fig. 4 (b), R-12 is affected more by surface roughness than R-134a. - Effects of pressure and roughness on the boiling curve are more effective at lower pressures as can clearly be seen in Figs. 3 and 4. This trend confirms a suggestion previously made by Cooper (1984).
- Another important nucleate boiling parameter is the heating surface material. As previously suggested, present research involved tests with three materials: copper, brass and stainless steel. Figure 5 shows boiling curves of refrigerants R-123 and R-134a on surfaces made up of the aforementioned materials. The effect of the heating surface material over the boiling curve is strikingly significant, the extent depending upon the particular fluid. It can be noted in Fig. 5 (a) for refrigerant R-123 that differences between boiling curves for copper and brass surfaces are minimal as compared to that for stainless steel. On the other hand, in the case of refrigerant R-134a, Fig. 5 (b), boiling curves corresponding to different materials present comparable differences between each other, with copper presenting better performance and stainless the worst.
- Results from part 5 suggest the influence of a combined surface/liquid effect instead of an isolated surface material effect, as suggested by some authors. In fact, if this were the case, boiling curves from Fig. 5 (a) for copper and brass would present comparable differences between each other as those from Fig. 5 (b), opposite to what actually is observed in these figures. This conclusion is in accordance with the proposed Rohsenow surface/liquid effect, characterized through the
*C*numerical coefficient._{sf}

**Correlation of Results**

Boiling curves as the ones shown in Figs. 3 to 5 include the following three regions with increasing heat flux: (1) natural convection; (2) partial nucleate boiling; and (3) fully developed nucleate boiling. Data considered for correlation purposes should be only those corresponding to fully developed nucleate boiling. However, in some previous studies experimental data for correlation purposes have been chosen by arbitrarily selecting a minimum heat flux. This heat flux would be one that would guarantee the occurrence of nucleate boiling and, as a result, the elimination of natural convection data. Thus, by adopting this procedure, data for both partial and fully developed nucleate boiling would be included in the set chosen for the fitting procedure. A more rigorous procedure would involve the determination of the onset of fully developed nucleate boiling in such a way that only experimental conditions corresponding to this boiling regime would be considered for fitting purposes. However, it must be recognized that the onset of this regime is not very well established since it depends upon parameters such as surface material and roughness as well as the particular boiling liquid, as clearly displayed in Figs. 3 to 5. Both procedures have been followed in this paper in order to compare results from each one regarding the level of correlation of experimental data. Thus, in the first case, a minimum heat flux was chosen in an "ad hoc" manner so that natural convection related data were not included in the set used for fitting purposes. A preliminary analysis of data allowed us to choose the aforementioned minimum heat flux as being equal to 5 kW/m^{2}. The second procedure involved a preliminary examination of experimental data along with the proposition of a criterion for the onset of fully developed nucleate boiling, related to the maximum value of the derivative of the boiling curve (*q"* vs. *D**T*). Initially, the variation of the derivative, *dq"/d(DT)*, of the boiling curves was raised followed by the determination of the onset heat flux and the selection of data points corresponding to fully developed nucleate boiling. The heat flux at the onset of fully developed nucleate boiling varied between 3 and 30 kW/m^{2}, depending upon the refrigerant, pressure and surface material and roughness.

Curve fitting of experimental results was performed according to the following procedures: (1) by keeping the original values of *m* and *n* proposed by Rohsenow, i. e. 0.33 and 1.7, and determining the coefficient *C _{sf}* for each refrigerant/surface material and roughness combination; (2) by letting

*C*,

_{sf}*m*, and

*n*free.

*C _{sf} Evaluation for m=0.33 and n=1.7*

In this case, fitting consisted in determining the coefficient *C _{sf}* for the refrigerant/surface material and roughness combinations considered in present study. Results are summarized in Table 3 for both

*(q")*5 kW/m

_{min=}^{2}and fully developed nucleate boiling data sets. The following general conclusions can be drawn from this table:

*C*values are not significantly affected by the data set considered in the fitting procedure, except in the case R-11 boiling on the stainless steel surface. However the absolute average deviation is generally higher for the_{sf}*(q")*5 kW/m_{min=}^{2}correlation with respect to its corresponding data set than that for the fully developed nucleate boiling one with respect to its associated data set.- The effect of the liquid/surface material combination over
*C*, previously suggested by Rohsenow and others, can clearly be noted. Unfortunately data from other sources are scarce. Thus a detailed comparison with present results is rather questionable since surface material and roughness, and operational conditions are generally not compatible._{sf}

A significant effect of the surface roughness *(Ra)* over *C _{sf}* can be noted in the table, a trend which had not been previously observed. As a general rule,

*C*diminishes with the roughness of the surface, a result which is consistent with the observed trends in the boiling curves. It must be noted that the originally proposed dependency upon the surface/liquid combination by the part of

_{sf}*C*is related to the wetting level of the liquid with respect to the heating surface. Thus

_{sf}*C*must keep a closed relationship with the contact angle, which in turn depends upon the surface roughness, as suggested elsewhere, Bikerman (1970), Kandlikar (2001).

_{sf}**General Procedure**

A second and more general fitting procedure has been pursued consisting in the determination of three parameters of the Rohsenow's correlation: *m*, *n* and *C _{sf}*. The adopted procedure has been performed according to the following steps.

**· m Determination**

Initially all data points were plotted in terms of curves *St** vs. *Re _{b}*. Values of

*m*for each refrigerant/surface combination have been determined with most of the results varying in a relatively narrow range. This is a clear indication that the refrigerant, surface material and roughness, and pressure do not affect significantly the value of

*m*. The plot of Fig. 6 of all fully developed nucleate boiling data points in terms of

*St**vs.

*Re*seems to confirm this trend. It can be noted that the cloud of data points is characterized by a common general slope, related to the

_{b}*m*value of Eq. (2). It must be stressed that, by doing so, individual trends involving different refrigerants, pressures, and surface material and condition are overviewed. According to Fig. 6, the resultant value of

*m*is 0.21, whereas for data corresponding to

*(q")*5 kW/m

_{min=}^{2},

*m*is equal to 0.18. These results allow one to conclude that the inclusion of partial nucleate boiling data do not affect significantly the value of

*m*.

**· n and C_{sf} Determination**

The characteristic value of *n* has been determined by plotting data in terms of vs *Pr _{l}*, as shown in Fig. 7 for the fully developed data set. It can be noted that, despite significant dispersion, plotted data follow a common trend characterized by a single slope, illustrated by the superposed straight line in the plot. Given that

*n*is related to this slope, one can immediately conclude that, as a first approximation, its value is not significantly affected either by the refrigerant, surface condition and material or pressure. The resulting value of

*n*is 1.03, for fully developed nucleate boiling data, and 1.15, for data corresponding to

*(q")*5 kW/m

_{min=}^{2}. Similarly to

*m*, differences in the value of

*n*associated to the data sets considered in this paper are relatively small.

It is interesting to note at this point that the liquid Prandtl number for refrigerant R-12 presents an odd behavior compared to that of the other three refrigerants in the range of pressures considered in present study, as shown in Fig. 8. It can be noted that whereas *Pr _{l}* diminishes with pressure for refrigerants R-11, R-123, and R134a, for refrigerant R-12,

*Pr*presents a minimum value at a reduced pressure of the order of 0.1. In addition, the relative variation of liquid refrigerant R-12 Prandtl number in the experimental pressure range (~ 10%) is much smaller than that for the other refrigerants (higher than 25%). Thus, the inclusion of R-12 experimental data in Fig. 7 affects the value of n due to the aforementioned shift in the slope with the reduced pressure though the relative variation of

_{l}*Pr*for this refrigerant is of limited extent. As an example of this effect, the value of

_{l}*n*changes from 1.15 to 1.55 when R-12 data are removed from a plot similar to that of Fig. 7 for the

*(q")*5 kW/m

_{min=}^{2}data set. In present study the value of

*n*was obtained including R-12 data for the sake of generality.

The coefficient *C _{sf}* of the Rohsenow's correlation has been obtained as the value where the

*(St*/Re*vs

_{b}^{m})*Pr*curve crosses the vertical axis. The dispersion of data in the plot of Fig. 7 is a clear indication that

_{l}*C*is affected not only by the liquid/surface combination, as previously suggested by Rohsenow and succeeding researchers, but also by the surface roughness,

_{sf}*Ra*, as pointed out in Section 5.1. As a result, data points of the plot of Fig. 7 must be sorted so that the effects of the aforementioned parameters over

*C*can quantitatively be determined. An example of the kind of results obtained by this procedure can be seen in Fig. 9, where the variation of C

_{sf}_{sf}with the reduced pressure for different surface conditions and refrigerants R-123 and R-11 boiling over copper surfaces is shown. It can be noted that, as a rule,

*C*diminishes with both the reduced pressure and surface roughness. In addition, the higher the roughness of the surface the lesser the effect of pressure.

_{sf}Individual *C _{sf}* curves for each refrigerant and surface material similar to those of Fig. 9 have been obtained. These results have been fitted by a correlation with the following general form:

(3)

valid for the refrigerants and refrigerant/surface combinations considered in this paper. Individual values of coefficients *a* to *d* have been obtained, requiring the use of one equation for each refrigerant. The procedure could be simplified by fitting data extensive to all the refrigerants with a single equation of the same general form as Eq. (3). Results form this general equation correlate data with adequate precision though not being as accurate as the particular ones as should be expected. The general equation has been adopted in this paper, according to it the coefficients *a* to *d* are constant and the refrigerant/surface material combinations are taken by the coefficient *C*. Table 4 presents values of the obtained coefficients from the data sets corresponding to fully developed nucleate boiling and *(q") _{min=}* 5 kW/m

^{2}.

Table 4 indicates that differences between the values of *C _{sf}* resulting from both data sets are minimal. Considering that the differences between exponents

*m*and

*n*are also small, as previously noted, one can conclude that the proposed correlations are very close to each other as the plots in Figs. 10 (a), (b) clearly display. The correlations are represented in terms of the heat transfer coefficient,

*h*, versus the heat flux,

*q"*, in these figures. The heat transfer coefficient is determined as

(4)

The maximum deviation between results from both correlations is of the order of 6.0% for case of Fig. 10 (a), for refrigerant R-11 and low reduced pressure and roughness, *p _{r}*=0.023 and 0.1 µm. The case of Fig. 10 (b), for the same refrigerant as in case (a), and

*p*=0.1 and

_{r}*Ra*=2.5 µm, presents a maximum deviation of the order of 11.0%. As expected, differences are minimal between the proposed correlations, a result that could explain why some of the literature correlations developed from data sets involving partial nucleate boiling data produce reasonable results.

**Evaluation of the Proposed Correlation**

In the previous section the procedure for fitting data was explained in detail leading to a general Rohsenow's correlation. The next step is to evaluate the proposed correlation with respect not only to data obtained as part of the research reported in this paper but also to data raised elsewhere.

In order to evaluate the performance of the proposed correlations with respect to experimental data, the relative deviation between the correlation and the experimental heat transfer coefficients is plotted against the corresponding heat flux. This plot is shown in Figs. 11 (a), (b), respectively for the fully developed nucleate boiling and *(q") _{min=}* 5 kW/m

^{2}data sets. It can be noted that most of the plotted data points fall within the ± 20% range, especially for the fully developed nucleate boiling data. Most of the data outside this range though of limited extent correspond to refrigerants R-123 and R-134a. It can also be noted that, as general rule, in the range of reduced heat fluxes the correlation tends to under predict the heat transfer coefficient, showing the opposite trend for higher heat fluxes. These trends are also discernible in the plots of Figs.12 (a), (b), where the correlation heat transfer coefficient,

*h*, is plotted against the experimental,

_{corr}*h*. The absolute average deviation of results from the fully developed nucleate boiling correlation with respect to their experimental counterparts is of the order of 10.1% whereas that relative to the

_{exp}*(q")*5 kW/m

_{min=}^{2}data set is of the order of 11.2%. These results are clearly noticeable in both figures, Figs. 11 and 12. Results displayed in these figures demonstrate the adequacy of Rohsenow's type of correlation in fitting nucleate boiling of halocarbon refrigerants over cylindrical surfaces.

Data from independent sources have been used in the evaluation of the proposed correlation. For refrigerant R-11 boiling on a brass surface, data from Silva (1989) have been used. Data from Webb and Pais (1992) for the same refrigerants considered in present paper boiling over copper surfaces have also been used for comparison purposes with the correlation for *(q") _{min=}* 5 kW/m

^{2}. An average roughness,

*Ra*, of 0.6 m m has been assumed for the surfaces since the authors only informed that the working surface was a commercial one. The comparison can be seen in Figs. 13 (a) and (b) in

*h*vs

_{corr}*h*plots. It can be noted that the proposed correlation fits very well data from both studies, with deviations being well within the ± 20% range, a result that can be considered adequate since the only adjustment made in the application of the proposed correlation was the assumption of the surface roughness.

_{exp}Data from Silva (1989), for refrigerant R-113 and a brass surface, and from Jensen (1985), for the same refrigerant and a stainless steel surface, have been fitted using the correlations proposed in the present investigation. The resulting values of the coefficient *C*, Eq. (3), are the following: 0.9 for brass, and 1.5, for stainless steel. The average roughness, *Ra*, of both surfaces has been assumed as being equal to 0.6 m m. Given the uncertainty of the experimental data regarding the roughness of the surfaces, the plots of Figs. 14 (a), (b), in terms of *h _{corr}* vs

*h*, display reasonable results from the correlation. Results from Silva (1989) for a brass surface, Fig,. 14 (a), are not as well correlated as the ones from Jensen (1985), Fig. 14 (b), for a stainless steel surface. In fact, the absolute average deviation of correlation heat transfer coefficient with respect to the experimental ones from Silva (1989) is of the order of 27% whereas with respect to those from Jensen (1985) is 6.2%.

_{exp}

**Conclusions**

The research reported herein has been set forth in order to investigate the performance of the Rohsenow type of correlation when applied to nucleate boiling of halocarbon refrigerants over cylindrical surfaces of different material. Two fitting procedures have been tried, including one consisting in keeping the exponents originally proposed by Rohsenow and adjusting the coefficient, *C _{sf}*. The other procedure, considered as general, fitted the exponents as well as the coefficient. It has been determined that exponents

*m*and

*n*are weakly affected by the refrigerant, and surface material and finishing. Liquid/surface combinations and surface finishing affect the coefficient

*C*. A correlation has been developed for

_{sf}*C*in terms of the liquid/surface combination and surface roughness. For the purpose of curve fitting, two data sets have been considered: one arbitrarily chosen to include data corresponding to heat fluxes higher than 5 kW/m

_{sf}^{2}whereas the other involved only fully developed nucleate boiling conditions. Though slightly different values of

*C*,

_{sf}*m*and

*n*have resulted from both data sets, differences in results could be considered minimal for all practical purposes. The correlations have been evaluated through the experimental results used in the fitting process and results obtained elsewhere. Deviations between the correlation and the experimental heat transfer coefficient are within acceptable ranges.

**Acknowledgements**

The authors gratefully acknowledge the support given to the reported research by the Fundação de Amparo à Pesquisa do Estado de São Paulo, FAPESP, and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, CAPES, Brazil. The technical support given to this investigation by Mr. José Roberto Bogni is also appreciated and recognized.

**References**

Abernethy, R. B., and Thompson, J. W., 1973, *Handbook, uncertainty in gas turbine measurements,* Arnold Engineering Development Center, Arnold Air Force Station, Tennessee

Bikerman, J. J., 1970, *Physical Surfaces,* Academic Press, Inc*,* New York.

Collier, J. G., and Thome, J. R., 1994, *Convective Boiling and Condensation*, 3rd Edition, Oxford University Press, Oxford.

Cooper, M. G., 1984, "Heat flow rates in saturated nucleate pool boiling - A wide ranging examination using reduced properties", *Advances in Heat Transfer,* Vol. 16, pp.157-238.

Jensen, M. K., 1985, "Physical Properties Data Tables and Experimental Boiling Data", Personal Communication.

Kandlikar, S. G., 2001, "A Theoretical Model to Predict Pool Boiling CHF Incorporating Effects of Contact Angle And Orientation", *Journal of Heat Transfer*, Vol. 123, pp. 1071-1079.

Kartsounes, G. T., 1975, "A Study of Treatment on Pool Boiling Heat Transfer in Refrigerant 12", *ASHRAE Transactions*, Vol. 18, pp. 320-326.

Liaw, S., and Dhir, V. K., 1993, "Void Fraction Measurements During Saturated Pool Boiling of Water on Partially Wetted Vertical Surfaces", *Journal of Heat Transfer*, Vol. 111, pp.731-738.

Pioro, I. L., 1999, "Experimental Evaluation of Constants for the Rohsenow Pool Boiling Correlation", *Int. J. Heat Mass Transfer*, Vol. 40, pp. 2379-2392.

Ribatski, G., 2002, "Theoretical and Experimental Analysis of Pool Boiling of Halocarbon Refrigerants", Doctoral thesis, Escola de Engenharia de São Carlos, University of São Paulo, Brazil.

Rohsenow, W. M., 1952, "A Method of Correlating Heat Transfer Data for Surface Boiling Liquids", *Transactions of ASME*, Vol.74, pp. 969-976.

Saiz Jabardo, J. M., and Silva, C. L., 1991, "Modeling of the nucleate boiling of refrigerant-oil mixtures", in: *18 ^{th} International Congress of Refrigeration*, Montreal, Canada.

Sauer, Jr. H. J., Medrow, R. A., and Sinnarwalla, A. M., 1975, "Effects of Surface Condition on Nucleate Boiling of Refrigerant-11", *ASHRAE Transactions*, Vol. 81, part. 2, pp. 274-28.

Silva, C. L., 1989, "*Experimental Investigation of the Nucleate Boiling of Refrigerant-Oil Mixtures"*, Doctoral thesis, Escola Politécnica, University of São Paulo, Brazil.

Tewari, P. K., Verma, R. K., and Ramani, M. P. S, 1986, "Effect of Surface Roughness and Polymeric Additive on Nucleate Pool Boiling at Sub-Atmospheric Pressures", *Int. Comm. Heat mass transfer*, Vol. 13, pp. 503-514.

Vachon, R. I., Nix, G. H., and Tanger, G. E., 1968, "Evaluation of Constants for the Rohsenow Pool-Boiling Correlation", *Journal of Heat Transfer*, Vol. 90, pp. 239-247.

Webb, R. L. and Pais, C., 1992, "Nucleate Pool Boiling Data for Five Refrigerants on Plain, Integral-Fin and Enhanced Tube Geometries", *Int. J. Heat Mass Transfer*, Vol. 35, pp. 1893-1904.

Paper accepted March, 2004. Technical Editor: Clóvis Raimundo Maliska