Experimental and Analytical Study on the Influence of Saturation Pressure andSurfaceRoughness onPoolBoilingCHFofHFE-7100

Critical heat ux (CHF) determines the safety and application of pool boiling heat transfer in many industrial elds.e inuence of saturation pressure and surface roughness on saturated pool boiling CHF in hydrouoroether HFE-7100 was experimentally studied in this investigation. Visualization and heat transfer measurements were conducted at the critical and transition boiling state, and further, the accuracy of CHF prediction models and enhancement mechanism had been analyzed. e polished boiling surfaces had various surface roughness values ranging from 0.019 to 0.587 μm and their contact angles varied from 7° to 10°, while the experimental saturation pressure changed from 0.7 to 2.0 bar.e visual images showed that the pool boiling phenomenon at a critical state was composed of dierent-sized bubbles, vapor column, and large mushroom vapor, whereas the unsteady blanket of vapor continually injected bubbles at a transition state. e saturation pressure and surface roughness had an obvious improvement on pool boiling CHF, which might be ascribed to the eects of bubble momentum owing to evaporation, distribution and recovery period of a heat transfer boundary layer, capillary action of the working liquid, as well as ratio of vapor jets’ area. Compared with the well-known correlations reported in the literature, CHF correlation of Bailey et al. (2006) predicted the current results more accurately. To further improve the prediction accuracy, a new empirical correlation for CHF dimensionlessK considering the eects of saturation pressure and surface roughness was developed, and the predicted values were in better agreement with the experimental data.


Introduction
Pool boiling heat transfer is an e cient way of transferring high thermal loads in a variety of industries including power generation sectors, nuclear plants, refrigeration plants, and cooling of electronic equipment. e triggering mechanisms of critical heat ux (CHF) are complex and depend on a number of factors including properties of uid [1,2], characteristics of engineered surfaces [3,4], and operating conditions of experiments [5], which have signi cant inuence on the accumulation of bubbles and supplement of liquid phase on boiling surfaces. e system saturation pressure is an important parameter, which has an e ect on the thermophysical properties of working media and dynamic parameters of bubbles, and it plays a signi cant role on CHF. More attempts have been made to investigate the impact of pressure on pool boiling process at a critical state. Dahariya and Betz [6,7] conducted pool boiling experiments for water on horizontal smoothed and sintered-particle wick copper surfaces at pressures ranging from 0 to 413.7 kPa, and the results showed that higher pressure would result in the enhancement of CHF due to the formation of a thermal boundary layer, merging of vapor trains into neighboring trains and modulation of wavelength. Sakashita [8] measured CHF in saturated pool boiling for water and TiO 2 nano uid on a vertical copper surface at pressures of 0.1∼0.8 MPa, and they indicated that CHF of water and a nano uid increased with the increase of pressure, while CHF enhancement of a nano uid would decrease with the increasing pressure. Similar to the abovementioned experimental results, Mudawar and Anderson [9], Alvarino et al. [10], and Guan et al. [11] also reported that pool boiling heat transfer performance and CHF would be improved with the increasing system pressure. e conclusions on the effect of pressure are nearly the same, whereas the influence of surface roughness may be controversial. Alvarino et al. [10] pointed out that CHF variations with surface roughness revealed an increment with roughness up to a certain limit and diminishing for the surface with higher roughness. Kim et al. [12][13][14] investigated the effect of surface roughness on pool boiling of water on metallic surfaces with different wettability, and the experimental results showed that the positive impact of surface roughness for superhydrophilic and hydrophilic surfaces on CHF was more obvious than that for hydrophobic surfaces, and they considered that the enhancement of CHF could be attributed to capillary wicking. Walunj and Sathyabhama [15,16] performed pool boiling experiments on the samples with the surface roughness varying from 0.106 to 4.03 μm in saturated water at the pressures of 1, 5, and 10 bar, and their findings displayed that CHF increased with the increase of roughness as well as pressure, and steady-state CHF was found to be higher than transient CHF for all the samples. Ferjančič and Golobič [17] also demonstrated that surface roughness had a positive effect on pool boiling CHF for both FC-72 and water and the influence of boiling surface chemistry should also be considered. Shin et al. [18,19] prepared a rough micronano hybrid surface by the electrodeposition method and pointed out that the combined effects of nanowire and microcavity delayed bubble coalescence and maximized bubble density, resulting in the enhancement of pool boiling CHF. Chu et al. [20] investigated the effect of roughness-augmented wettability on CHF and indicated that roughness-amplified capillary forces were responsible for CHF enhancement on microstructured pillar surfaces. However, some experimental results suggested that surface roughness did not significantly affect boiling CHF. Berenson [21] measured a pool boiling curve for n-pentane at atmospheric pressure and considered that CHF was essentially independent of surface material, roughness, and cleanliness due to the difference between the average and local heat flux. O'Hanley et al. [22] investigated the effects of surface wettability, porosity, and roughness on pool boiling CHF of water, and they indicated that porosity enhanced CHF, whereas the wettability and roughness had little effects on CHF.
In order to provide the guidance for engineering applications, researchers also proposed large number of models and correlations to predict boiling CHF, including hydrodynamic instability model [23], macrolayer dryout model [24], dry spot model [25], interfacial lift-off model [26], forces balance model [27,28], and empirical correlations [29,30]. Most of the models and correlations were developed according to triggering mechanisms and/or experimental database of pool boiling CHF. e involved factors influencing CHF in the models and correlations have been reported and include thermophysical properties of working fluid, boiling surface characteristics, wettability, heat flux, surface orientation, and system pressure. Further and in-depth study on the effect of the above factors on pool boiling CHF has a significant theory and application values.
As mentioned above, saturation pressure and surface roughness are two important factors influencing pool boiling CHF. Currently, there are few studies on the combined influence of the two parameters on heat transfer. However, it is valuable to obtain the combination effects of saturation pressure and surface roughness, as well as analyze mechanisms for CHF enhancement in more detail. Meanwhile, the predictive CHF correlation is essential to further consider both effects. erefore, this study investigated the impacts of saturation pressure and surface roughness on pool boiling CHF of hydrofluoroether dielectric fluid HFE-7100 (methoxy-nonafluorobutane, C 4 F 9 OCH 3 ) and quantitatively analyzed the mechanisms for CHF enhancement with previous classical models. e selected working fluid HFE-7100 has wide applications in industry, such as heat transfer fluid and spray contact cleaner, owing to its low global warming potential (GWP 320), zero ozone depletion potential, chemical and thermal stability, high wettability, low toxicity, and nonflammability [31]. Visualization and heat transfer measurements of pool boiling have been conducted in this study. Experiments were carried out at four different pressures of 0.7, 1.0, 1.5, and 2.0 bar on four horizontal copper surfaces with the average roughness of 0.019, 0.205, 0.311, and 0.587 μm. e CHF values were also compared with the existing correlations and a new predictive correlation for CHF dimensionless K was established to make recommendations for design.

Experimental Apparatus and Method.
e experimental facility, as demonstrated in Figure 1, consisted of a copper block, a boiling chamber with an auxiliary heater, a vapor condenser, a water/R-134a heat exchanger, water and R-134a cooling system. Saturated vapor of HFE-7100 generated from boiling surface entered the condenser unit, with the condensate returning to the boiling chamber by gravity via a filter. Deionized water was used as a cooling medium in the condenser. e deionized water was cooled in the heat exchanger using a R-134a vapor compression refrigeration unit. e boiling chamber was a vertical stainless steel 304 cylinder with an inner diameter of 22 cm and equipped with two circular glass windows mounted on its sides for visualization. A pressure transducer and three thermocouples were applied to monitor the saturation state of the working fluid in the chamber. e chamber was thermally insulated and had a 1250 W auxiliary heater rubbed between the rubber insulation and the chamber to facilitate the maintenance of system saturated conditions. Heating was provided by six 250 W cartridge heaters with the diameter of 10 mm, which were installed in the lower part of the copper heater block, as seen in Figure 1. e heat input was regulated using a variable transformer and measured by a power meter. Six type-K thermocouples were inserted in the center of the copper block and keeping them 5 mm apart in the vertical direction helped calculate the heat flux to a boiling surface. e test surface was the upper surface of a disc, with 40 mm diameter and 5 mm height placed on the top of the 2 International Journal of Chemical Engineering copper block. A solder was used to maintain a good thermal contact between the bottom surface of the disc and top surface of the copper block. e experimental system was tested for leaks and degassed before use. All the experimental data were recorded after the system reached a steady state, which was con rmed when the system pressure and heat ux remained stable for 10 minutes (i.e., the uctuation of temperature and pressure were less than 1 K and 3 kPa, respectively). Tests were performed at the system pressures of 0.7, 1.0, 1.5, and 2.0 bar on smooth and rough copper surfaces. Saturated conditions were achieved by adjusting the ow and inlet temperature to the condenser and the electric input to the auxiliary heater. Experimental data were obtained for increasing heat ux, from nucleate boiling to transition boiling. In this study, the maximum heat ux at given operating conditions was considered as CHF, and the critical point indicated the boiling process would change from nucleate boiling to transition boiling regime.

Preparation and Characterization of Boiling Surfaces.
Four copper surfaces with di erent roughness were fabricated for the pool boiling experiment. e smooth copper surface was prepared using a diamond turning machine (250UPL, Nanotech), while the rough copper surface was polished by sandpapers with di erent grain sizes under the same weight (500 g). e surface characteristics were then studied using a scanning electron microscope (1450VP, LEO), a 3D optical surface pro lometer (ZeGage Plus, ZYGO) with the morphology repeatability less than 0.15 nm, and a contact angle meter (DSA100, KRüSS GmbH) with the uncertainty of ±0.1°.
Surface topography for boiling surfaces with di erent roughness values was detected by SEM, as shown in Figure 2. Compared with a smooth surface, the grooves of rough surfaces were more obvious, and the depth and width of grooves increased with the decrease of sandpaper grain size.
To quantitatively analyze the surface roughness of copper surfaces, a 3D optical surface pro lometer was used to measure average surface roughness (Ra), root mean square of surface feature heights (Rq), and average roughness peak distance (Rsm), and the results are listed in Table 1. e average surface roughness Ra was 0.019∼0.587 μm. e average value of the length of pro le element along the sampling length, Rsm, was 8.733∼14.334 μm, and the root mean square of surface feature heights, Rq, was 0.024∼0.767 μm.
A contact angle meter was applied to record the static contact angle θ of HFE-7100 on copper surfaces. Figure 3 displays the contact angle θ was around 7∼10°at 20°C, indicating that HFE-7100 had good wettability on copper surfaces. e contact angle would decrease with the increasing roughness, due to the fact that the liquid ran down the bigger grooves easily.

Data Reduction.
In this study, six thermocouples were employed to measure the temperature distribution in the copper block. e temperature gradient was then applied to calculate the boiling heat ux (q) as follows: where k cu is the thermal conductivity of copper block and dT/dy| y 0 is the vertical temperature gradient at the top surface of disc. CHF is considered as the maximum heat ux at the stable state under given experimental conditions. System calibration was rst performed to ensure the reliability of experiments. e thermocouples were calibrated using a platinum resistance thermometer (F250, OMEGA) with the uncertainty of ±0.1 K, and the uncertainty of test location was ±0.2 mm. e pressure transducers were also calibrated using a dead weight pressure gauge tester (Bryans Aeroquipment LTD) within 2 kPa. An error propagation analysis was also conducted according to the Mo at model [32]. e calculated uncertainty in the critical heat ux was between 1.5% and 2.7%.  International Journal of Chemical Engineering transition boiling, and lm boiling. e heat transfer performance at a critical boiling state is superior to the others with much more dramatic relative motion between vapor and liquid, indicating that it sets the upper limit of the nucleate boiling regime. In this study, the macroscopic visualization images for pool boiling at the critical state or transition state were collected, analyzed, and compared.

Results and Discussion
Pool boiling images were obtained for critical boiling (a-d) and transition boiling (e) as exhibited in Figure 4. At the critical boiling state, a two-phase interface uctuated violently due to high heat transfer ux, and the boiling twophase structures appeared as nucleated bubbles, large bubbles, a large vapor column, and a vapor mushroom cloud alternately.
e experimental macroscopic phenomenon under various operating conditions were similar; however, pool boiling process would achieve a critical state at lower heat ux under lower pressure for a boiling surface with smaller roughness. Bubbles with various sizes always covered the boiling surface. Small nucleated bubbles on the boiling surface would rapidly grow and coalesce to form large irregular bubbles, and the growth rate of bubbles was nearly 1.13 mm/ms (see a 1 -a 5 ); large bubbles would further merge together to form a large vapor column on the boiling surface (see b 1 -b 5 ), which constructed the vapor ow passage; the center line of the vapor column would pulsate along the horizontal direction (see c 1 -c 5 ) and inject a vapor mushroom cloud (see d 1 -d 5 ); meanwhile, the boiling surface was always covered with the bubbles of di erent sizes. As the boiling surface temperature further increased, the boiling process would terminate the critical state and enter the boiling transition regime. e boiling surface was covered with an uneven blanket of vapor and the heat transfer    Large vapor column

CHF Enhancement Mechanism
Analysis. e dynamic parameters of vapor bubbles and supplement of the liquid phase were the key factors influencing boiling CHF. Due to the limitation of viewing angle and camera distance, the microscale boiling process in a sublayer cannot be recorded, so the pool boiling microscale two-phase structure was constructed according to related physical models [11,12], as illustrated in Figure 6. e irregular big bubble, microscale vapor column, nucleated bubbles and liquid sublayer covered on the boiling surface, and the vapor of HFE-7100 moved upward constantly while the liquid rewetted the dry spots continuously to maintain the steady state of boiling. When the supplement and evaporation of working fluid were in equilibrium and the boiling heat flux reached the upper limit, the boiling process was regarded as a critical point. When nucleated bubbles and vapor column expanded continuously in the horizontal direction and covered most of the boiling surface, and the liquid could not supply effectively, a shift from the critical boiling state to transition boiling state would occur, resulting in the decrease of the heat transfer rate and the rapid increase of surface temperature correspondingly. erefore, pool boiling CHF was determined by spreading and rewetting of liquid on a boiling surface, liquid boundary layer thickness, horizonal momentum of bubbles, and vapor area ratio of the heated surface synergistically. e mechanism of saturation pressure and surface roughness on the above factors would be discussed as follows.
Saturation pressure has a significant effect on thermal properties of the working fluid and further influences the dynamic parameters of vapor phase and boundary layer of the liquid phase. Meanwhile, surface roughness has an impact on the spreading wettability of the liquid working medium on a boiling surface and then determines the effective replenishment of liquid.
Kandlikar [27] established a theoretical model to analyze the hydrodynamic behavior of a bubble on CHF conditions and indicated that the change in momentum owing to evaporation caused a two-phase interface to move rapidly along the boiling surface, leading to the initiation of CHF condition. e force due to the change in momentum could be expressed by the evaporation mass flow rate and vapor velocity and be calculated as follows: In formula (2), the momentum force F M along the horizontal direction was influenced by interface heat flux of bubbles q I , characteristic scale of bubbles D b and H b , and properties of working fluid ρ V , ρ L , h LV , σ, and static contact angle θ. Figure 7 shows the value of momentum force F M for smooth and rough copper surfaces under different saturation pressures according to the above force balance model. e results demonstrated that the saturation pressure had significant effects on momentum force F M , while the effect of surface roughness was less obvious, which owes to the effect that the fluid properties were influenced by pressure remarkably, whereas the contact angle was impacted by surface roughness slightly. e momentum force F M would increase with the increase of interface heat flux and decrease of saturation pressure. Larger momentum force F M caused the vapor bubbles to expand along the boiling surface and blanket the heater surface more easily, triggering CHF condition.
Haramura and Katto [24] developed a macrolayer dryout CHF model and described the two-phase structure at a boiling critical state that the vapor-liquid interface of columnar vapor stems distributed in a liquid layer wetting the heated surface, indicating that the growth of large bubbles was the result of consumption of macrolayer by evaporation, and the CHF condition was triggered by total evaporation of liquid on a boiling surface. Subsequently, Rajvanshi et al. [38] measured the initial macrolayer thickness for di erent liquids by the electrical resistance probe method, and their experimental results agreed with twice magnitude of the predictive values derived from Haramura and Katto's correlation [24], and they suggested the macrolayer thickness could be de ned as follows: In general, the typical nucleation pool boiling process is composed of bubble nucleation, bubble growth, bubble departure, and liquid rewetting. e liquid sublayer would change the periodicity during the process. Zhao and Williams [39] and Ding et al. [40] developed a prediction model for the recovery period of the liquid boundary layer, which was expressed as follows:   International Journal of Chemical Engineering Figures 8 and 9 display the in uence of system saturation pressure on the thickness and recovery period of the liquid sublayer, respectively. e results obtained from the abovementioned models showed that the liquid sublayer thickness increased with the increase in pressure, while the recovery period of liquid sublayer decreased with the increase of pressure. e thicker liquid sublayer and shorter recovery period were more bene cial to supply the liquid phase on a boiling surface and reduce the dramatic changes in temperature gradient of thermal boundary. erefore, it was easier to prevent the formation of dry spots and hinder coalesce, as well as the aggregation of vapor bubbles or jets on boiling surface under higher saturation pressure. e surface microscopic geometry characteristics can be in uenced by surface roughness evidently. Quan et al. [28] and Son and Kim [41] pointed out that the augmenting surface roughness improved the ability of liquid spreading on a surface and e ectively supplied liquid to the heated surface, thus delaying the occurrence of CHF.
In order to investigate the e ect of surface roughness on dynamic spreading process of working uid HFE-7100, an experiment in the impacting of droplets on copper surfaces was conducted in this study. e droplets with the same volume (20 μL) were dropped on the copper surfaces with di erent roughness values by a microsyringe, and the distance between the microsyringe and copper surface was 15 mm. e spreading of HFE-7100 on di erent surfaces is shown in Figure 10. Initially, the droplet would collapse on the copper surface and change into a liquid lm. After nearly 7 ms, it entered the spreading stage, and the liquid lm would extend continuously along three-phase contact line until it reached the maximum spreading length. Figure 11 exhibits the spreading length of the liquid lm on copper surfaces at di erent times. It was found that the spreading length for rougher surfaces was larger than that of the smooth surface, with the average spreading velocity of 0.178, 0.188, 0.191, and 0.213 mm/ms for Ra of 0.019, 0.205, 0.311, and 0.587 μm, respectively. is implied that the liquid working uid was easier to ll with the microstructure with larger scale and crush into larger gaps on the copper surface, which was conducive to provide e ective liquid supplement.
Similar to our nding, Son and Kim [41] also experimentally investigated the characteristics of receding capillary ow under adiabatic conditions and demonstrated that the precondition for capillary ow was due to the fact that the capillary pressure of liquid overcome recoil pressure of gas, and a rougher surface was bene cial to enhance capillarity and restrain expansion of dry area. Meanwhile, they theoretically developed a model to predict the advancing capillary ow velocity considering mass ux balance at the liquid-vapor interface under evaporation conditions. e equation was expressed as follows:  Figure 12 displays the values of capillary ow velocity under evaporation conditions from equation (5). e results showed that the in uence of surface roughness on capillary ow velocity was much more obvious than that of saturation pressure, and the advancing capillary ow velocity increased with the increase of roughness. Higher capillary ow velocity resulted in larger liquid supplement, which could compensate a higher evaporative potential and enhance the boiling CHF. e ratio of vapor jets area on boiling surface and area of heated surface A V /A W at a critical state is also an important parameter, determining the distribution of vapor and liquid on the boiling surface, which impacts the velocity di erence at a vapor-liquid interface, as described by Haramura and Katto [24]. e equation was expressed as follows:   Figure 13 presents the area ratio A V /A W for boiling surfaces with di erent roughness values under four saturation pressures. It could be seen that the area ratio A V / A W increased with the increasing saturation pressure and surface roughness. e greater region could be allowed to be occupied by vapor jets, indicating that both e ects e ectively controlled the liquid-vapor phase separation, owing to more liquid supplement. e area ratio A V /A W was in proportion to liquid supplement capacity, based on continuity momentum balance. As discussed above, thicker liquid sublayer and higher capillary ow velocity under the condition of higher saturation pressures and larger surface roughness would provide more e ective liquid replenishment and enhance the uid ow momentum.

Models and Correlations of Pool Boiling CHF.
Prediction accuracy of pool boiling CHF is vital to safety and reliability of cooling system. Numerous pool boiling CHF models and correlations have been developed by theoretical analysis and/or correlating experimental data, which are useful within the given range of experimental conditions. In the present study, the experimental CHF results of HFE-7100 had been compared with the values predicted by related models and correlations as summarized in Table 2.   Table 2 are shown in Figure 14. As depicted in this gure, the mean absolute deviation ranged from 5.38% to 42.9%. e predicted results of correlations derived from Guan et al. [11], Bailey et al. [30], Bailey et al. [30], and Kandlikar [27] were in good agreement with CHF values for boiling surfaces with the roughness of 0.019 μm, 0.205 μm, 0.311 μm, and 0.587 μm under four saturation pressures, and the mean absolute deviation was 4.58%, 3.31%, 2.60%, and 2.68%, respectively. In summary, the correlation developed by Bailey et al. [30] could predict the current data on CHF more accurately, with the absolute deviation of 5.38% for all experimental operating conditions.
Most models and correlations used parameter K de ned by Kutateladze [43] to predict CHF. In the proposed equation (7), CHF parameter K re ected the e ect of experimental conditions, such as wettability, surface roughness, saturation pressure, surface orientation, liquid properties, etc., while other items illustrated the impact of thermophysical properties of the working uid. As shown in Table 2, parameter K was usually expressed as dimensionless parameters or constants.
In this study, both saturation pressure and surface roughness had a positive e ect on pool boiling CHF, so the CHF correlation required to be correlated with the two e ects. erefore, the dimensionless parameter Pr (reduced pressure), Ra/Rsm (the ratio between average surface roughness and average roughness peak distance), and ρ v /ρ L (density ratio) could be introduced to establish a new correlation. Using the curve tting method based on experimental data, a new predicting empirical correlation for CHF was proposed as follows.
e new correlation would be applied for sandblasted or polished metal surfaces under the following conditions, 0.019 μm ≤ Ra ≤ 0.587 μm, 0.03 ≤ Pr ≤ 0.09, β < 10°. Figure 15 demonstrates the comparison between experimental data and predictive values from correlations for CHF parameter K. Compared with the experimental results, although some models were more accurate in predicting the values of CHF parameter K under individual conditions, the predicted values of K from the new correlation established in this study were in better agreement with the overall experimental values. e average absolute deviation range of CHF parameter K for the previous correlations was from 5.66% to 42.35%, whereas the average absolute deviation for the new correlation was 2.72%, indicating that the prediction accuracy was signi cantly improved.
In order to certify the reliability of the new correlation established in this study, the experimental data from literature and present study were used to compare with the predicted values. e selected working media with high wettability and pool boiling experimental conditions are listed in Table 3. e comparison results between the CHF experimental data (working uid: HFE-7100 [10], FC-72 [11,44], Pentane [11], PF-5060 [45,46], HFE-7000 [46]) and the predicted values from the new correlation are displayed in Figure 16. As seen in the gure, the new correlation could predict the CHF of working uids well and the absolute deviation of the most experimental data was within 10%.
e comparing results showed that the new developed correlation could provide more convenience in practical engineering applications for pool boiling CHF of highly wetting working media.
Interfacial lift-off model, considering the instability of a liquid-vapor interface and wetting front behavior on boiling surfaces, reflecting the influence of fluid thermophysical properties.
Bailey et al. [30] Empirical correlation, experimental data of pentane, methanol, and water， reflecting the influence of fluid thermophysical properties.
Guan et al. [11] q CHF � 0.2445 · ((ρ Macrolayer lift-off model, describing the slice of the liquid microlayer, considering the momentum, conservation, and Laplace condition, reflecting the influence of fluid thermophysical properties, especially the density of the working fluid.

Conclusions
Pool boiling of HFE-7100 on copper surfaces with di erent roughness values ranging from 0.019 to 0.587 μm under four di erent saturation pressures of 0.7, 1.0, 1.5, and 2.0 bar at a critical state was presented in this study. e conclusions drawn from the experimental results of visualization and heat transfer were summarized as follows: (1) Pool boiling macroscopic visualization images at the CHF point was composed of bubbles with di erent size, large vapor column, and vapor mushroom cloud. e movement of vapor and liquid was dramatic, and the patterns showed multiple variations. e non-at blanket of vapor at the transition state would cover the boiling surface and continually inject the dispersed bubbles and the dynamic process was relatively simple.
(2) Saturation pressure and surface roughness had a positive impact on pool boiling CHF, and the e ect of saturation pressure was more remarkable. On the basis of mechanism analysis, it indicated that saturation pressure had a noticeable e ect on the momentum force of vapor bubbles, thickness, and recovery time of the liquid layer, while surface roughness would improve the ability of liquid spreading and both of them a ected the area ratio between vapor jets and boiling surface. (3) e present CHF results compared well with the predictions by Bailey et al. (2006) for the entire pressure ranged covering sub-atmospheric pressure of 0.7 bar to 2.0 bar and the surface roughness ranging from 0.019 to 0.587 μm, and the mean absolute deviation was 5.38%. In order to consider the combined e ect of saturation pressure and surface roughness and improve the prediction accuracy, a new correlation for CHF dimensionless parameter K was developed, and the predicted values were in good agreement with the experimental data in literature and this study. CHF dimensionless parameter P: Pressure (Pa) q: Heat ux (kW m − 2 ) Rq: Root mean square of the surface feature heights (μm) y: Vertical distance (m) A W : Area of heated surface (m 2 ) D b : Bubble diameter (m) g: Gravitational acceleration (m s − 2 ) h LV : Latent heat of vaporization (kJ kg − 1 ) k: ermal conductivity (W m − 1 K − 1 ) Pr: Reduced pressure Ra: Average surface roughness (μm) Rsm: Average roughness peak distance (μm) Greek Symbols   [45,46] Square, 10 × 10 mm 2 Copper 0.039∼0.58 0.05 HFE-7000 [46] Square, 10 × 10 mm 2 Copper 0.039∼0.58 0.04 Experimental data of HFE-7100 [10] Experimental data of FC-72 [11,41] Experimental data of Pentane [11] Experimental data of PF-5060 [42,43] Experimental data of HFE-7000 [43] 450 Figure 16: Comparison of experimental pool boiling CHF data and predicted values from new correlation.