Abstract
Flow boiling heat transfer around the critical heat flux (CHF) condition at high reduced pressures of carbon dioxide in a 296-μm hydraulic diameter microchannel was experimentally studied. The CHF conditions for developing flow and fully developed flow were measured and compared to established correlations. The post-CHF heat transfer coefficient was obtained for l/d of 3.2, 7.4, and 11.6 for inlet Reynolds numbers, based on the homogeneous two-phase flow model, ranging from 6622 to 32,248. The critical heat flux conditions seemed to peak around a reduced pressure of about 0.5 and gradually decreased with reduced pressure. However, the typical rapid increase in the surface temperature following the CHF condition decreased with increasing pressure, and the post-CHF heat transfer coefficient was appreciably high (up to about 50 kW/m2K) at high reduced pressures. The enhancement in the heat transfer coefficient and CHF condition near the inlet were quantified. The experimental results were compared to established CHF correlations and heat transfer coefficient correlations with some limited success. Thus, the Katto CHF correlation (Katto and Ohno, 1984, “An Improved Version of the Generalized Correlation of Critical Heat Flux for the Forced Convective Boiling in Uniformly Heated Vertical Tubes,” Int. J. Heat Mass Transfer, 27(9), pp. 1641–1648) and the Bishop correlation (Bishop et al., 1964, “Forced-Convection Heat Transfer to Water at Near-Critical Temperatures and Supercritical Pressures,” Westinghouse Electric Corp, Atomic Power Division, Pittsburgh, PA.) for the post-CHF heat transfer coefficient were adjusted to better predict the experimental results. Additionally, an enhancement factor was derived to predict the increase in the heat transfer coefficient in the developing region.
1 Introduction
Flow boiling in microchannels has been a topic of much fundamental and practical interest over the better part of the 21st century. However, there are several major roadblocks hindering its wide use in cooling systems for high-power electronics. Chief among them are flow boiling instabilities, the critical heat flux (CHF) condition, high-pressure drop, temperature overshoot, high superheat temperature at the onset of nucleate boiling (ONB), and the prevalence of laminar flow, which is associated with a low-grade heat transfer coefficient. As a plethora of systems require low pumping power (also low-pressure drop) and low fluid inventory, it is desired to leverage the latent heat of vaporization to its fullest extent by enabling high mass quality flows (i.e., reduced flowrate, and thus, reduced pressure drop and low fluid inventory). However, the CHF conditions impose a stringent requirement on the ability of these systems to operate at desired high mass qualities. As a result, more recent research efforts have started to explore means to operate microscale flow boiling systems at mass qualities as high as 0.95 [1]. As part of this endeavor, several studies have explored the flow boiling of carbon dioxide (CO2) in microchannels at a reduced pressure of 0.5 and higher, as it offers several distinct advantages. Unlike many potential coolants, the critical temperature of CO2 of 31.1 °C is well aligned with the temperature domain permissible by electronic systems, and high reduced pressure has been shown to mitigate flow boiling instabilities [2–7]. Carbon dioxide has a low viscosity, allowing microscale cooling systems to operate with the enhancement offered by turbulent flow. It also has a low surface tension, which minimizes the ONB, eliminates temperature overshot, and helps mitigate the rapid bubble growth instability. In addition, at high reduced pressures, the thermal conductivity of the vapor can be comparable to that of liquid, albeit still lower, helping mitigate some of the harsh penalties imposed by the critical heat flux condition (i.e., the transition to post-CHF might not be so harmful since the vapor can still support heating loads at acceptable surface temperatures). (It should be noted that CO2 has a global warming potential of 1, which is an order of magnitude smaller than hydrofluorocarbons, such as R134a with a global warming potential of 1430.)
Equipped with the realization of these potential benefits, Parahovnik et al. [8–10] and Parahovnik and Peles [11–13] embarked on a journey to reveal the nature of flow boiling heat transfer of CO2 in microchannels at high reduced pressures. Initially, subcooled flow boiling was examined, and the ONB [8], subcooled heat transfer coefficient [8,14–17], saturated flow boiling patterns [12,18,19], bubble dynamics [11], the piston effect [9], and CHF conditions were examined [12]. Subsequently, saturated flow boiling heat transfer was studied and compared with previous models and correlations [12]. In parallel, Adeoye et al. [20,21] and Adeoye and Peles [22] studied flow boiling heat transfer with micro-impingement jets and elucidated the heat transfer coefficient and critical heat flux. Additionally, Asadzadeh et al. [23] investigated heat transfer characteristics in a silicon microchannel with a NACA pin fin configuration.
The current study extends these previous efforts and provides an in-depth analysis of the transition through the CHF condition. It delineates and quantifies the transition from pre-CHF condition to post-CHF condition and the associated heat transfer toll related to the transition. Specifically, it provides a measure of the criticality of the critical heat flux condition (i.e., the severity of the temperature increase, or the drop in the heat transfer coefficient, during the transition). The coefficients and exponents of the Katto CHF correlation [24] and the Bishop et al. heat transfer coefficient correlation [25] were examined and subsequently modified to better predict the experimental results.
2 Experimental Setup and Methods
2.1 Microfluidic Device and Experimental Setup.
The microfluidic device was made from two parts: a substrate with resistive temperature detectors (RTDs) and heaters, termed part 1, and a second custom-made substrate, part 2, fabricated by Mindrum Precision® with a 296-μm hydraulic diameter microchannel and an embedded 81-μm hydraulic diameter orifice with a length of 1.8 mm (i.e., l/d ratio of about 22). The width and height of the orifice were 500 μm and 45 μm, respectively, and the width and height of the microchannel were 1800 μm and 160 μm (Fig. 1(a)). (Note that the orifice had slightly rounded edges due to tooling and computer-aided design (cad) software and a measurement of the area and the perimeter out of there resulted in a hydraulic diameter of 81 μm rather than 82.6 μm.)

(a) An isometric sketch of the microfluidic device with the microchannel's inlet, the microchannel (Dh = 0.3 mm), the orifice, and the microchannel's outlet. (b) The package, pressure transducer, thermocouple, and mass flux. (c) The layout of the heater and the RTDs on part 1 (d) a layout of the orifice, heater, and RTDs positions.

(a) An isometric sketch of the microfluidic device with the microchannel's inlet, the microchannel (Dh = 0.3 mm), the orifice, and the microchannel's outlet. (b) The package, pressure transducer, thermocouple, and mass flux. (c) The layout of the heater and the RTDs on part 1 (d) a layout of the orifice, heater, and RTDs positions.
Part 1 was made of a 0.5-mm thick fused silica plate, which included the electrical circuit (i.e., the RTDs, heaters, electrical vias, and electrical pads). These were made by the deposition of metallic layers (i.e., titanium, platinum, and aluminum) and multiple silicon oxide layers for electrical insulation. The RTDs were custom-made from vapor-deposition layers. The heater and RTDs were made using similar microfabrication processes but with much different geometrical shapes to render their functionality. The rectangular heater had dimensions of 5.82 mm by 0.89 mm (Fig. 1(b)) with a low electrical resistance of about 70 Ω, generating heat when the DC power was turned on. The serpentine RTDs had footprints of 50 μm by 60 μm, with a resistance of about 2.5k Ω each. The RTDs were placed on top of the heater at locations specified in Fig. 1(b). The linear dependency of the metal resistivity of the RTDs on the temperature was leveraged to measure the local, temporal temperature of the heated wall.
Figure 1(c) depicts a top view of part 1, where the microchannel, inlet, outlet, and orifice were marked with a dashed white line. The flow is shown inside the dashed area, and the restriction in the middle represents the orifice. The heater is shown as a darker area inside the dashed region, split into two independently controlled zones (upstream and downstream of the orifice). (Only the downstream heater was used in this study.) The heater and the RTDs were made from metallic layers of titanium (7 nm) and platinum (30 nm) deposited on a fused silica wafer and later sliced into rectangular dies. These components were placed inside the microchannel and had intimate contact with the flow. The electrical path included electrical vias, electric pads, spring pins, custom-made printed circuit boards (PCBs) mounted on the package, and standard nine-pin connectors to sense and control the in-channel elements (i.e., the heater and the RTDs). The vias had an additional aluminum layer (∼1 μm) to significantly reduce their electrical resistance and render it negligible compared to the RTDs and heaters. For insulation, the metallic layers were covered with a silicon oxide layer with a thickness of ∼1.1 μm that was later selectively etched to expose the electrical pads that interfaced with the push pins. The power was supplied from a DC power supply (Kysight®, E3645A®) through electrical wiring, push pins, contact pads, and the vias.
The microfluidic device was enclosed in a special housing that kept its structural integrity and enabled sealing between the two pieces while supplying the flow at desired experimental conditions (Fig. 1(d)). Additionally, the package enabled to read the temperatures for three RTDs located at length to hydraulic diameter ratios of 3.2, 7.4, and 11.6, and to apply the desired thermal heat load. The two-phase flow conditions were achieved using a throttling process through the embedded micro-orifice and an external preheater (HeatBrisk®). The orifice allowed to study of a wide range of mass qualities and pressures in the microchannel, and its operational design and thermodynamic states were described in detail by Parahovnik et al. [10]. Because of the low mass flow rates through the microfluidic device, the carbon dioxide was allowed to flow through the open loop experimental setup and was released to the ambient. The experimental flow setup was intended to provide maximum control and flexibility to study a range of conditions, such as pressure, inlet temperature, and mass flux. The experimental setup had pressure transducers (Omega®, Omega Engineering Inc., Norwalk, CT), a mass flowmeter (Alicat®), and a high-precision T-type Thermocouple that enabled to measure of pressure, mass flowrate, and fluid temperature at the inlet of the microdevice, respectively. Additional details about the pressure and mass flow control are provided by Parahovnik et al. [8,11].
2.2 Experimental Procedure.
2.3 Uncertainty Analysis.
where ΔTsat = TRTD–Tsat and q″loss = (TRTD–Tamb)/Rt.
The mass of the microfluidic device and the microdeposited components were small such that the RTDs and the heater's thermal time response were practically instantaneous. Nevertheless, data was recorded only 30 s after a new heat load was set, allowing ample time for the temperature to stabilize. Subsequently, the data were sampled for an additional 20 s with a sampling frequency of 1 kHz and were averaged to obtain the measured point. The standard deviation was also calculated and added to the zero-shift accuracy to estimate the uncertainty of the reading.
3 Results and Discussion
The experimental data points, given in Table 1, cover pressures from 1.7 MPa to 6.9 MPa (i.e., reduced pressures of 0.23–0.935) and mass fluxes from 1405 kg/m2s to 2776 kg/m2s (based on the microchannel's cross sections of A = 0.284 × 10−6 m2), corresponding to a Reynolds number based on the homogeneous two-phase flow model ranging from 6622 to 32,248. Three length-to-diameter ratios were studied, l/d, of 3.2, 7.4, and 11.6. As the flow exited the orifice, it gradually expanded while passing through the first two RTDs (i.e., l/d = 3.2 and 7.4) before becoming fully developed at the third RTD at l/d = 11.6. (Since the flow was turbulent, it is believed that around l/d∼10 the flow was fully developed [27].)
List of experimental conditions used
Critical heat flux (W/cm2) | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
Exp. # | P (MPa) | Pr | Tin (K) | G (kg/m2s) | x | Re | @RTD1 | @RTD2 | @RTD3 | q″max (W/cm2) |
1 | 1.7 | 0.23 | 250.0 | 1405 | 0.66 | 12,098 | 146.7 | 155.1 | ||
2 | 3.5 | 0.47 | 274.4 | 2200 | 0.27 | 9876 | 214.7 | 157 | 274.4 | |
3 | 3.6 | 0.49 | 275.4 | 1693 | 0.94 | 12,856 | 191.6 | 161.5 | 225.8 | |
4 | 3.8 | 0.51 | 277.5 | 1867 | 0.8 | 12,452 | 243.3 | 203.4 | 146.9 | 206.4 |
5 | 4.0 | 0.54 | 279.4 | 2256 | 0.36 | 13,496 | 138.1 | 277.5 | ||
6 | 4.8 | 0.65 | 286.6 | 2141 | 0.14 | 12,797 | 186.5 | 154.7 | 219 | |
7 | 4.9 | 0.67 | 287.4 | 2207 | 0.15 | 15,112 | 195.5 | 175.9 | 140.2 | 216.9 |
8 | 5.1 | 0.70 | 289.1 | 1990 | 0.13 | 14,884 | 137.2 | 176.7 | ||
9 | 5.4 | 0.73 | 291.4 | 2164 | 0.08 | 11,828 | 199.6 | 160.4 | 133.1 | 210.3 |
10 | 5.5 | 0.75 | 292.2 | 2366 | 0.052 | 10,902 | 143.1 | 133.4 | 211.8 | |
11 | 5.6 | 0.76 | 292.9 | 1851.2 | 0.04 | 12,498 | 134.8 | 117.6 | 134.0 | |
12 | 5.6 | 0.76 | 292.9 | 2068 | 0.05 | 6622 | 186.8 | 140.1 | 131.2 | 209.8 |
13 | 5.8 | 0.78 | 293.2 | 2118.6 | –0.02 | 11,972 | 145.4 | 135.4 | 126.3 | 135.0 |
14 | 5.8 | 0.78 | 295 | 2467 | 0 | 13,159 | 167.8 | 133.2 | 101.3 | 211.4 |
15 | 6.0 | 0.81 | 295 | 2558 | –0.06 | 13,843 | 162 | 126 | 102.5 | 236.5 |
16 | 6.0 | 0.81 | 295 | 2041 | –0.03 | 12,092 | 120.5 | 97.4 | 164.8 | |
17 | 6.0 | 0.81 | 295.8 | 2427 | 0 | 32,248 | 167.8 | 134 | 109.7 | 222.3 |
18 | 6.0 | 0.82 | 293.2 | 2776 | –0.08 | 17,678 | 140 | 118.6 | 170.8 | |
19 | 6.2 | 0.85 | 296.4 | 2521 | –0.08 | 20,860 | 148.1 | 115.8 | 86.72 | 213.4 |
20 | 6.5 | 0.89 | 296.6 | 2333 | –0.21 | 24,230 | 129.5 | 90.45 | 75.94 | 235.9 |
21 | 6.7 | 0.91 | 296.3 | 2654 | –0.27 | 18,057 | 92.86 | 74.15 | 243 | |
22 | 6.9 | 0.94 | 297.3 | 2400 | –0.5 | 31,728 | 75.03 | 59.1 | 268.1 |
Critical heat flux (W/cm2) | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
Exp. # | P (MPa) | Pr | Tin (K) | G (kg/m2s) | x | Re | @RTD1 | @RTD2 | @RTD3 | q″max (W/cm2) |
1 | 1.7 | 0.23 | 250.0 | 1405 | 0.66 | 12,098 | 146.7 | 155.1 | ||
2 | 3.5 | 0.47 | 274.4 | 2200 | 0.27 | 9876 | 214.7 | 157 | 274.4 | |
3 | 3.6 | 0.49 | 275.4 | 1693 | 0.94 | 12,856 | 191.6 | 161.5 | 225.8 | |
4 | 3.8 | 0.51 | 277.5 | 1867 | 0.8 | 12,452 | 243.3 | 203.4 | 146.9 | 206.4 |
5 | 4.0 | 0.54 | 279.4 | 2256 | 0.36 | 13,496 | 138.1 | 277.5 | ||
6 | 4.8 | 0.65 | 286.6 | 2141 | 0.14 | 12,797 | 186.5 | 154.7 | 219 | |
7 | 4.9 | 0.67 | 287.4 | 2207 | 0.15 | 15,112 | 195.5 | 175.9 | 140.2 | 216.9 |
8 | 5.1 | 0.70 | 289.1 | 1990 | 0.13 | 14,884 | 137.2 | 176.7 | ||
9 | 5.4 | 0.73 | 291.4 | 2164 | 0.08 | 11,828 | 199.6 | 160.4 | 133.1 | 210.3 |
10 | 5.5 | 0.75 | 292.2 | 2366 | 0.052 | 10,902 | 143.1 | 133.4 | 211.8 | |
11 | 5.6 | 0.76 | 292.9 | 1851.2 | 0.04 | 12,498 | 134.8 | 117.6 | 134.0 | |
12 | 5.6 | 0.76 | 292.9 | 2068 | 0.05 | 6622 | 186.8 | 140.1 | 131.2 | 209.8 |
13 | 5.8 | 0.78 | 293.2 | 2118.6 | –0.02 | 11,972 | 145.4 | 135.4 | 126.3 | 135.0 |
14 | 5.8 | 0.78 | 295 | 2467 | 0 | 13,159 | 167.8 | 133.2 | 101.3 | 211.4 |
15 | 6.0 | 0.81 | 295 | 2558 | –0.06 | 13,843 | 162 | 126 | 102.5 | 236.5 |
16 | 6.0 | 0.81 | 295 | 2041 | –0.03 | 12,092 | 120.5 | 97.4 | 164.8 | |
17 | 6.0 | 0.81 | 295.8 | 2427 | 0 | 32,248 | 167.8 | 134 | 109.7 | 222.3 |
18 | 6.0 | 0.82 | 293.2 | 2776 | –0.08 | 17,678 | 140 | 118.6 | 170.8 | |
19 | 6.2 | 0.85 | 296.4 | 2521 | –0.08 | 20,860 | 148.1 | 115.8 | 86.72 | 213.4 |
20 | 6.5 | 0.89 | 296.6 | 2333 | –0.21 | 24,230 | 129.5 | 90.45 | 75.94 | 235.9 |
21 | 6.7 | 0.91 | 296.3 | 2654 | –0.27 | 18,057 | 92.86 | 74.15 | 243 | |
22 | 6.9 | 0.94 | 297.3 | 2400 | –0.5 | 31,728 | 75.03 | 59.1 | 268.1 |
The experiments were sorted according to their reduced pressure. Empty cells correspond to conditions in which the critical heat flux wasn't reached. For saturated flow boiling, Tin was calculated based on the saturation temperature, and for subcooled flow boiling, the liquid inlet temperature was used.
3.1 Transition Through the Critical Heat Flux and Post-Critical Heat Flux Heat Transfer Coefficient.
Figure 2 depicts the boiling curves for reduced pressures of 0.67 and 0.91 at three different length-to-diameter ratios. As the flow was still developing, the region closest to the orifice exhibited the highest CHF value, a trend that is also consistent with the model of Katto [24] (although outside of the correlation boundaries), and the highest heat transfer coefficient. (Regardless of the developing flow versus developed flow, CHF typically starts downstream at the trailing edge of the heater and then propagates upward, but since the mass quality didn't change much along the heater, it is believed that the higher CHF near the inlet was primarily because of the effect of the developing flow.) At low reduced pressures, the transition across the CHF condition can be rapid and accompanied by a large increase in surface temperature (Fig. 2(a)). As a matter of proper engineering practice, the CHF condition should be typically avoided. However, it was found that as the reduced pressure increased, the transition to post-CHF conditions (Fig. 2(b)) was gradually accompanied by only a moderate increase in the surface temperature. At the same time, the absolute CHF value decreased with pressure. The vapor-to-liquid ratios of the thermal diffusivity, dynamic viscosity, and density approach unity at high reduced pressure smoothen the transition from pre-CHF to post-CHF conditions.
![Boiling curves for reduced pressure of 0.67 (experiment #7) (a) and reduced pressure of 0.91 (experiment #21) (b). The temperature increased at RTD3 when the CHF was reached for a reduced pressure of 0.67 and for a reduced pressure of 0.91 was 71.6 K and 9 K, respectively. Predictions from the Katto [24] correlation and the current correlation are also shown.](https://asmedc.silverchair-cdn.com/asmedc/content_public/journal/heattransfer/145/4/10.1115_1.4056821/1/m_ht_145_04_041606_f002.png?Expires=1681270095&Signature=qUrYd8vNFYrAFKs8iRMUb861CWORmfL3MKA71RpMdj4KlHIBbxn7FxUsMeBu-M2Ws5PZPRJ4nF6BEV66iwimIgYHCbhm~4jHw16kpZ9N9VM160YbcLFrTTpGxFre6Fz42FRWtqQS0~LdOl3Av5XtbgaPLdMa4uT4X-pAczqRBm3fjSeOrjAPLokdoPj3oxkCQL7Q5DY8GyFbkrZj-CLMswrT~rjmmgXzIjYlMO4UvVHkEjyKPWC1g0IPDRSiXAH2jkezIuhWOauuK1u9P2oviPiwRKaemZ5mCnfncsHfd-2XIVovT6aMHCwiN4lxkNs2mROA-Aj~jALrtxDQI-xalA__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
Boiling curves for reduced pressure of 0.67 (experiment #7) (a) and reduced pressure of 0.91 (experiment #21) (b). The temperature increased at RTD3 when the CHF was reached for a reduced pressure of 0.67 and for a reduced pressure of 0.91 was 71.6 K and 9 K, respectively. Predictions from the Katto [24] correlation and the current correlation are also shown.
![Boiling curves for reduced pressure of 0.67 (experiment #7) (a) and reduced pressure of 0.91 (experiment #21) (b). The temperature increased at RTD3 when the CHF was reached for a reduced pressure of 0.67 and for a reduced pressure of 0.91 was 71.6 K and 9 K, respectively. Predictions from the Katto [24] correlation and the current correlation are also shown.](https://asmedc.silverchair-cdn.com/asmedc/content_public/journal/heattransfer/145/4/10.1115_1.4056821/1/m_ht_145_04_041606_f002.png?Expires=1681270095&Signature=qUrYd8vNFYrAFKs8iRMUb861CWORmfL3MKA71RpMdj4KlHIBbxn7FxUsMeBu-M2Ws5PZPRJ4nF6BEV66iwimIgYHCbhm~4jHw16kpZ9N9VM160YbcLFrTTpGxFre6Fz42FRWtqQS0~LdOl3Av5XtbgaPLdMa4uT4X-pAczqRBm3fjSeOrjAPLokdoPj3oxkCQL7Q5DY8GyFbkrZj-CLMswrT~rjmmgXzIjYlMO4UvVHkEjyKPWC1g0IPDRSiXAH2jkezIuhWOauuK1u9P2oviPiwRKaemZ5mCnfncsHfd-2XIVovT6aMHCwiN4lxkNs2mROA-Aj~jALrtxDQI-xalA__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
Boiling curves for reduced pressure of 0.67 (experiment #7) (a) and reduced pressure of 0.91 (experiment #21) (b). The temperature increased at RTD3 when the CHF was reached for a reduced pressure of 0.67 and for a reduced pressure of 0.91 was 71.6 K and 9 K, respectively. Predictions from the Katto [24] correlation and the current correlation are also shown.
The temperature increase following the transition to post-CHF condition (i.e., ΔTCHF) was inversely proportional to the reduced pressure (Fig. 3(a)) with a linearity Pearson coefficient [28] of –0.92, –0.9, and –0.74 for l/d of 3.2, 7.4, and 11.6, respectively. The coefficient indicates a strong linear dependency between the parameters. The CHF above a reduced pressure of 0.5 gradually declined (Fig. 3(b)). As the difference between the thermophysical properties of the liquid and the vapor diminished with pressure, ΔTCHF approached 0 K near the critical condition [29]. Although the ΔTCHF asymptotically approaches zero at the critical pressure, the CHF did not, and its lower limit was about 50 W/cm2. This is somewhat surprising since the surface tension and the differences between the properties of liquid and vapor diminish with pressure, which should eliminate any difference between the pre-CHF to post-CHF conditions at the critical state. This lower limit is hypothesized to be related to the piston effect, a unique thermalization mode of near-critical fluids [9].
Since the velocity distribution varied significantly from RTD1 to RTD3, the heat transfer process was considerably different. As such, the heat transfer coefficient and the CHF varied between RTD1 and RTD3 (Fig. 4). Since the velocity near the wall at RTD1 was much higher than downstream, the heat transfer coefficient was considerably higher.

Post-CHF heat transfer coefficients for RTD1, RTD2, and RTD3 located at l/d of 3.2, 7.4, and 11.6, respectively. The average heat transfer coefficients for RTD1, RTD2, and RTD3 were ∼50 kW/m2K, ∼30 kW/m2K, and ∼16 kW/m2K, respectively.
3.2 Critical Heat Flux Correlation.
As shown in Fig. 5, for reduced pressures above 0.6 (i.e., density ratio above 0.15), the correlation underpredicted the critical heat flux condition for all RTDs. The correlation overpredicted the experimental results for reduced pressures below 0.6 for RTD3. The MAE values for RTD1, RTD2, and RTD3 were 20.7%, 16.9%, and 12.4%, respectively.
Comparison of the fitted coefficients and exponents for the new CHF correlation
Coefficient or exponent . | Katto [24] . | New value . |
---|---|---|
a | 0.234 | 0.35 |
b | 0.27 | 0.581 |
e | 0.0031 | 0.155 |
f | 0.513 | 0.750 |
g | –0.433 | −0.405 |
Coefficient or exponent . | Katto [24] . | New value . |
---|---|---|
a | 0.234 | 0.35 |
b | 0.27 | 0.581 |
e | 0.0031 | 0.155 |
f | 0.513 | 0.750 |
g | –0.433 | −0.405 |
Since the exponents had the strongest impact on the accuracy of the correlation, a sensitivity test of the exponents f and g was performed and is shown in Fig. 6. Each exponent was independently varied while the other was kept constant at its optimum value, and the average MAE was obtain based on results for all three RTDs. The MAEs were sensitive to the coefficients, and a shift of 10% from the optimal values doubled the MAEs. Figure 7 compares the modified correlation (Eq. (3)) with the experimental results for the three length-to-diameter ratios. The CHF peaked around Pr = 0.5 and gradually decreased with pressure. (For comparison, Boyd [32] reported a peak value at a reduced pressure of 0.75 for water.) However, as stated earlier, the lower CHF at high reduced pressures resulted in a moderate decline in the post-CHF heat transfer coefficient.

Sensitivity study for the fitted exponents f and g. The MAEs of the three RTDs were averaged to calculate the overall correlation accuracy with adjusted exponents (i.e., MAEf and MAEg). A deviation of 10% from the optimal coefficient value doubled the MAE.
![Comparison of the new correlation (Eq. (3)) with the experimental results. The MAEs values for RTD1, RTD2, and RTD3 were 5.86%, 4.84%, and 6.69%, respectively, compared to 20.7%, 16.9%, and 12.4% for the original Katto [23] correlation. The correlation better predicted the HTC for all reduced pressures and for the three RTDs.](https://asmedc.silverchair-cdn.com/asmedc/content_public/journal/heattransfer/145/4/10.1115_1.4056821/1/m_ht_145_04_041606_f007.png?Expires=1681270095&Signature=oGyKyyZQp-~tg2Hf~A9rZoB2JRhO3cYLdYQMi4Rjsv-~ay9QOOlF-kMHiyyYdAxPp2fgFGWS7Wm4Yfna7zvKUobYLK61o1YYwYoECx0sLNxduFz3PvfU9PYAOb7C~h0CLWdjJ0VHzF6oH4yb-OBjNwQB8~lupVRQzRi~IaJCyIIZDd6KuxMMQ8rjh5NHvyo5q4bquY1IMGQ5-Dk0sLE07ikpaoUx3sOD1oWfVLFRmsDop14ECx2qg1t0oSb8M6PwdQrWjPkPqnj2qS0vWid~rSYFvUBeVDcZMx-~SVVTB4-aasP8wPxcfN1P2P~svIdqkBz8S2PmX1XE0sJjSnqRDg__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
Comparison of the new correlation (Eq. (3)) with the experimental results. The MAEs values for RTD1, RTD2, and RTD3 were 5.86%, 4.84%, and 6.69%, respectively, compared to 20.7%, 16.9%, and 12.4% for the original Katto [23] correlation. The correlation better predicted the HTC for all reduced pressures and for the three RTDs.
![Comparison of the new correlation (Eq. (3)) with the experimental results. The MAEs values for RTD1, RTD2, and RTD3 were 5.86%, 4.84%, and 6.69%, respectively, compared to 20.7%, 16.9%, and 12.4% for the original Katto [23] correlation. The correlation better predicted the HTC for all reduced pressures and for the three RTDs.](https://asmedc.silverchair-cdn.com/asmedc/content_public/journal/heattransfer/145/4/10.1115_1.4056821/1/m_ht_145_04_041606_f007.png?Expires=1681270095&Signature=oGyKyyZQp-~tg2Hf~A9rZoB2JRhO3cYLdYQMi4Rjsv-~ay9QOOlF-kMHiyyYdAxPp2fgFGWS7Wm4Yfna7zvKUobYLK61o1YYwYoECx0sLNxduFz3PvfU9PYAOb7C~h0CLWdjJ0VHzF6oH4yb-OBjNwQB8~lupVRQzRi~IaJCyIIZDd6KuxMMQ8rjh5NHvyo5q4bquY1IMGQ5-Dk0sLE07ikpaoUx3sOD1oWfVLFRmsDop14ECx2qg1t0oSb8M6PwdQrWjPkPqnj2qS0vWid~rSYFvUBeVDcZMx-~SVVTB4-aasP8wPxcfN1P2P~svIdqkBz8S2PmX1XE0sJjSnqRDg__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
Comparison of the new correlation (Eq. (3)) with the experimental results. The MAEs values for RTD1, RTD2, and RTD3 were 5.86%, 4.84%, and 6.69%, respectively, compared to 20.7%, 16.9%, and 12.4% for the original Katto [23] correlation. The correlation better predicted the HTC for all reduced pressures and for the three RTDs.
The new correlation also better captured the experimental trends and predicted the experimental data at reduced pressures above 0.6. The MAEs values for RTD1, RTD2, and RTD3 were 5.86%, 4.84%, and 6.69%, respectively, compared to 20.7%, 16.9%, and 12.4% for the original Katto [24] correlation.
3.3 Post-Critical Heat Transfer Coefficient Correlation.
For pre-CHF conditions, the liquid was the dominant phase that directly carried the heat away from the wall through vaporization or convection [34]. For post-CHF conditions, convection to the vapor phase was an essential part of the heat transfer process, and the liquid phase served to a large extent as a thermal capacitor [1,34]. As the reduced pressure increased, the liquid-to-vapor thermodynamic property ratios, such as the thermal conductivity, thermal diffusivity, asymptotically approach unity, and the surface tension diminished [36,37]. Thus, the post-CHF heat transfer coefficient was expected to be appreciably higher than for low pressures, which could potentially support moderate heat fluxes at acceptable surface temperatures.
The experimental results were compared to five correlations detailed in Song et al. [35] to obtain quantitative measures of the post-CHF heat transfer coefficient under the high reduced pressures examined in this study, given in Table 2. Through the maximum root-mean-square (RMS) defined as √((hexp.–hcalc.)/hexp.)2). The Miropolskiy [38] and the Groeneveld [39] correlations had similar terms (i.e., vapor-to-liquid density ratio, mass quality, vapor-based thermal conductivity and Reynolds number, and wall-based Prandtl number) but different functional dependencies (i.e., coefficient and exponents). The governing terms in the Bishop [25] correlation were similar (i.e., density ratios, Prandtl and Reynolds numbers, and thermal conductivity), but the terms were mainly evaluated based on the wall and vapor conditions. The Quinn [40] correlation accounted for a viscosity ratio (i.e., μv/μl), which wasn't considered in the other correlations. The Song et al. [35] correlation was formulated through coefficients, which depend on the Reynolds number, Prandtl number, and mass quality that are used as multipliers of the Nusselt number calculated through the Dittus–Boelter correlation [27] (where the Reynolds number is based on the homogeneous two-phase flow model [34] and the Prandtl number is evaluated based on the wall temperature).
The correlations were compared to the heat transfer coefficients at RTD3 (i.e., the location corresponding to a fully developed flow). To verify that the correlations were compared with stable post-CHF conditions, the heat transfer coefficients were calculated for a wall temperature of at least 40 K above the saturation temperature. The correlations were compared to 258 data points collected from 22 experiments. The calculated MAEs were 75.2%, 91.1%, 36.2%, 74.7%, and 95.2% for the Miropolskiy [38], Groeneveld [39], Bishop et al. [25], Quinn [40], and Song et al. [35], respectively. The Bishop et al. [25] correlation best predicted the experiments, which can be because it evaluates the thermophysical properties based on the wall temperature, suggesting that the superheated vapor near the wall was the dominant heat transfer medium. All the other correlations were mainly based on the saturated properties, which failed to predict the experimental results.
![The C ratio between the experimental data and the Bishop et al. [25] correlation for RTD3. The ratio was a linear function of the quality.](https://asmedc.silverchair-cdn.com/asmedc/content_public/journal/heattransfer/145/4/10.1115_1.4056821/1/m_ht_145_04_041606_f008.png?Expires=1681270095&Signature=n79fvJgMDaBRB-qmzlzB4ceLHN3UJUij~hGLYTxo7uW4bWUwHskvZt4G4vr6RpFX2MNaaJNt0puC4zJdqqRe0cYScQumkF60g0sCgtt0tADMLrzWb17V40~dQ4Icd3tEWWFmlsJI-RG523YnlWrEEvf4FOsuZ3tPww5svh9SQjyMjdRvkj-aYth9b27pZDBnehgC1V2k5EKPBabFK~HX9jReIkX0RtY0uvpFwN1ryfvtjDEWkaNRKajgoW3OwnWMuiejeCd6NruQUFIC17t7IOgqKUXXJLURnD4e0QnxKTwjxdTpg6IaEaJkiA-bLnNKW7Zaihhjtjb5i8HT1m6h8g__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
The C ratio between the experimental data and the Bishop et al. [25] correlation for RTD3. The ratio was a linear function of the quality.
![The C ratio between the experimental data and the Bishop et al. [25] correlation for RTD3. The ratio was a linear function of the quality.](https://asmedc.silverchair-cdn.com/asmedc/content_public/journal/heattransfer/145/4/10.1115_1.4056821/1/m_ht_145_04_041606_f008.png?Expires=1681270095&Signature=n79fvJgMDaBRB-qmzlzB4ceLHN3UJUij~hGLYTxo7uW4bWUwHskvZt4G4vr6RpFX2MNaaJNt0puC4zJdqqRe0cYScQumkF60g0sCgtt0tADMLrzWb17V40~dQ4Icd3tEWWFmlsJI-RG523YnlWrEEvf4FOsuZ3tPww5svh9SQjyMjdRvkj-aYth9b27pZDBnehgC1V2k5EKPBabFK~HX9jReIkX0RtY0uvpFwN1ryfvtjDEWkaNRKajgoW3OwnWMuiejeCd6NruQUFIC17t7IOgqKUXXJLURnD4e0QnxKTwjxdTpg6IaEaJkiA-bLnNKW7Zaihhjtjb5i8HT1m6h8g__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
The C ratio between the experimental data and the Bishop et al. [25] correlation for RTD3. The ratio was a linear function of the quality.
A similar procedure to the one described in Sec. 3.2 was performed to examine the validity of the exponents in the Bishop et al. correlation. The exponents, Table 4, were examined for their MAEs while each exponent was varied independently.
Post-critical heat flux correlations
Correlation | Correlation | Maximum RMS compared to Song et al. [35] database. | Current database MAE RTD3 |
---|---|---|---|
Miropolskiy [38] | h = 0.023 · kv/Dh · (Rev · (xm + (ρv/ρl) · (1–xm))0.8 · Prw0.8 (1-0.1(ρv/ρl–1)0.4 · (1–xm)0.4) | 63.6 % | 75.2 % |
w-property evaluated on wall temperature Rev = G · Dh/μv | |||
xm, as defined in this table under Song et al. correlation. | |||
Groeneveld. [39] | h = 0.00109 · kv/Dh · (Rev · (xm+(ρv/ρl) · (1–xm))0.989 · Prw1.41 · (1–0.1(ρv/ρ–l)0.4 · (1–xm)0.4)–1.15 | 70.4% | 91.1% |
Rev= G · Dh/μv | |||
Bishop et al. [25] | h = 0.098 · kw/Dh · (Rew · (ρw/ρv) · (xm+(ρv/ρl) · (1–xm)))0.8 · Prw0.83 · (ρv/ρ–l)0.5 | 57.7% | 36.2% |
Rew = G · Dh/μw | |||
xm, as defined in this table under Song et al. [35] correlation. | |||
Quinn [40] | h = 0.023 · kv/Dh · (Rev · (xm+(ρv/ρl) · (1–xm))0.8 · Prw0.4 · (μv/μl)0.14 | 58.3% | 74.7% |
Rev = G · Dh/μv | |||
xm, as defined in this table under Song et al. correlation. | |||
Song et al. [35] | For x < 0 xm = 0 For 0 ≤ x ≤ 1 xm = x For x > 1 xm = 1 where: x = (Hb–Hl )/ (Hv–Hl) | 20.2% | 95.2% |
μ = (xm/μv+(1–xm)/μl)–1 - Viscosity for homogeneous two-phase flow model [34] ReTP = G · Dh/μ | |||
F1 = 8.346 × 10–3 · ReTP0.319 F2 = max {(0.795 · Prw–1.752, 0.809 · Prw1.287)} w -property evaluated on wall temperature | |||
For x < 0.4 F3 = 0.715 · (1–x)–0.678 For x > = 0.4 F3 = 1 | |||
Nu = F1 · F2 · F3 · (0.023 · ReTP0.8 · Prw0.4) h = Nu · kv/Dh |
Correlation | Correlation | Maximum RMS compared to Song et al. [35] database. | Current database MAE RTD3 |
---|---|---|---|
Miropolskiy [38] | h = 0.023 · kv/Dh · (Rev · (xm + (ρv/ρl) · (1–xm))0.8 · Prw0.8 (1-0.1(ρv/ρl–1)0.4 · (1–xm)0.4) | 63.6 % | 75.2 % |
w-property evaluated on wall temperature Rev = G · Dh/μv | |||
xm, as defined in this table under Song et al. correlation. | |||
Groeneveld. [39] | h = 0.00109 · kv/Dh · (Rev · (xm+(ρv/ρl) · (1–xm))0.989 · Prw1.41 · (1–0.1(ρv/ρ–l)0.4 · (1–xm)0.4)–1.15 | 70.4% | 91.1% |
Rev= G · Dh/μv | |||
Bishop et al. [25] | h = 0.098 · kw/Dh · (Rew · (ρw/ρv) · (xm+(ρv/ρl) · (1–xm)))0.8 · Prw0.83 · (ρv/ρ–l)0.5 | 57.7% | 36.2% |
Rew = G · Dh/μw | |||
xm, as defined in this table under Song et al. [35] correlation. | |||
Quinn [40] | h = 0.023 · kv/Dh · (Rev · (xm+(ρv/ρl) · (1–xm))0.8 · Prw0.4 · (μv/μl)0.14 | 58.3% | 74.7% |
Rev = G · Dh/μv | |||
xm, as defined in this table under Song et al. correlation. | |||
Song et al. [35] | For x < 0 xm = 0 For 0 ≤ x ≤ 1 xm = x For x > 1 xm = 1 where: x = (Hb–Hl )/ (Hv–Hl) | 20.2% | 95.2% |
μ = (xm/μv+(1–xm)/μl)–1 - Viscosity for homogeneous two-phase flow model [34] ReTP = G · Dh/μ | |||
F1 = 8.346 × 10–3 · ReTP0.319 F2 = max {(0.795 · Prw–1.752, 0.809 · Prw1.287)} w -property evaluated on wall temperature | |||
For x < 0.4 F3 = 0.715 · (1–x)–0.678 For x > = 0.4 F3 = 1 | |||
Nu = F1 · F2 · F3 · (0.023 · ReTP0.8 · Prw0.4) h = Nu · kv/Dh |
The database used to develop these correlations includes water and R134a with uniform or nonuniform heating. The experimental parameters for the Song et al. database were reduced pressures from 0.045 to 0.975, mass fluxes from 2.5 kg/m2s to 3117.9 kg/m2s, hydraulic diameters from 10 mm to 1923 mm, heat fluxes from 1 × 10−4 W/cm2 to 158.7 W/cm2, mass qualities from 0.129 to 499.68 kg/m2s, and wall superheats from 2.98 K to 540.82 K.
Exponents m, n, and s varied between 0.6–0.9, 0.7–1.5, and 0.3–0.8, respectively. The MAE was invariant to the exponent of the Prw (i.e., n), with an MAE variation lower than 1.25%. Variations of the m and s exponents had a significant influence on the MAEs. However, contrary to the CHF expression, the exponents, in this case, were consistent with the original Bishop et al. [25] correlation and, therefore, were not changed, Fig. 9.
Equation 5 predicted the experiments well for a reduced pressure of 0.49 and a reduced pressure of 0.78, Figs. 10(a) and 10(b), while it underpredicted the subcooled flow boiling experiment of reduced pressures of 0.85, 0.89, 0.91, and 0.94. Regardless, Eq. (5) captured well the experimental trend. For experiment #3, the thermophysical properties of the vapor were much inferior for post-CHF conditions resulting in a lower post-CHF heat transfer coefficient, while for experiment #22, with a reduced pressure of 0.94, the transition to post-CHF was relatively smooth and resulted in comparably high post-CHF heat transfer coefficient (i.e., 30 kW/m2K at ΔTw = 40 K). A possible explanation for the underprediction of the heat transfer coefficient for experiment #22 can be related to the presence of the piston effect, which for a reduced pressure of 0.94, is predicted to account for 25% of the heat transfer rate [9]. By assuming an average calculated heat transfer coefficient of 20 kW/m2K for experiment #22. By adding the poison effect into consideration (i.e., 25% of the heat was dissipated into the bulk fluid through the piston effect), the actual heat transfer coefficient should be around 25 kW/m2K, which is approximately the value obtained from the experiments, Fig. 10(c). Since the Bishop et al. [25] correlation doesn't account for near-critical condition effects, it underpredicted the near-critical heat transfer coefficients.

Three of the 22 experiments were used to evaluate Eq. (5). (a) Pr = 0.49 and x = 0.97 (experiment #3). (b)Experiment #14 corresponded to saturated flow boiling with Pr = 0.78 and x = 0. (c) experiment #22 for subcooled flow boiling with Pr = 0.94, x = –0.5, and ΔTsub = –5.8 K. In all cases, the heat transfer coefficients were within a 10 kW/m2K to 30 kW/m2K band.

Three of the 22 experiments were used to evaluate Eq. (5). (a) Pr = 0.49 and x = 0.97 (experiment #3). (b)Experiment #14 corresponded to saturated flow boiling with Pr = 0.78 and x = 0. (c) experiment #22 for subcooled flow boiling with Pr = 0.94, x = –0.5, and ΔTsub = –5.8 K. In all cases, the heat transfer coefficients were within a 10 kW/m2K to 30 kW/m2K band.
3.4 Post-Critical Heat Transfer Coefficient for Developing Flow.
The center velocity of a flow exiting an orifice is the highest and narrowest at the orifice exit. As the flow propagates downstream, the center velocity declines, but the nonzero velocity region becomes wider, as shown in the cone sketched in Fig. 1(c). The enhancement factor was derived by minimizing the MAEs for the RTDs at l/d of 3.2 and 7.4, Fig. 11. The enhancement factor for RTD1 (i.e., t1) and RTD2 (i.e., t2) was found to be 2.3 and 0.8 with MAEs of 38.5% and 27.5%, respectively.

The heat transfer enhancement factor (t) was used to predict the heat transfer coefficient in the developing region
4 Conclusions
The critical heat flux condition and post-CHF heat transfer coefficient of carbon dioxide at high reduced pressures in a microchannel were experimentally studied and compared to available correlations. Based on the experimental results, two new correlations were proposed for the post-CHF heat transfer coefficient and the CHF condition. It was found that as the reduced pressure increased, the CHF decreased, and the post-CHF heat transfer coefficient increased. As a result, the typical rapid increase in the surface temperature following the transition to film boiling was moderate, unlike lower reduced pressures where the burnout phenomenon is common. Moreover, a micro-orifice and its associated developing flow were demonstrated to enhance the heat transfer coefficient and the CHF condition significantly. Additionally, it was found that the post-CHF heat transfer coefficient prediction was best when the thermophysical properties were evaluated at the wall temperature. Consequently, a new correlation for CHF condition, based on Katto [24] formulation, and a new correlation for fully developed post-CHF heat transfer coefficient were proposed, resulting in reduced MAEs. An enhancement factor was derived to predict the increase in the heat transfer coefficient in the developing region.
Acknowledgment
This work was funded by the Office of Naval Research (ONR) under award number N00014-21-1-2653. The authors would like to acknowledge the Cornell Nano-Scale Science & Technology Facility (CNF) for staff support, help, and use of facilities to fabricate the microdevices. The CNF is a member of the National Nanotechnology Coordinated Infrastructure (NNCI), supported by the National Science Foundation.
Funding Data
Office of Naval Research (Award No. N00014-21-1-2653l Funder ID: 10.13039/100000006).
Data Availability Statement
The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.
Nomenclature
- a =
CHF correlation coefficient, Eq. (2)
- A =
channel's cross section (0.284 × 10−6 m2)
- b =
CHF correlation exponent, Eq. (2)
- C =
correction coefficient for the developing flow, Eq. (4)
- CHF =
critical heat flux (W/cm2)
- Dh =
hydraulic diameter (m)
- e =
CHF correlation coefficient, Eq. (2)
- Exp # =
experiment #, Table 1
- f =
CHF correlation exponent, Eq. (2)
- F1, F2, F3 =
Song et al. [31] coefficients
- g =
CHF correlation exponent, Eq. (2)
- G =
mass flux (kg/m2s)
- h =
heat transfer coefficient (kW/m2K)
- H =
enthalpy (J)
- hfg =
latent heat of vaporization (J/kg)
- k =
thermal conductivity (W/mK)
- l/d =
length-to-diameter ratio
- MAE =
mean averaged value, Eq. (1)
- N =
number of samples, Eq. (1)
- Nu =
Nusselt number
- ONB =
onset of nucleate boiling
- P =
pressure (MPa)
- Pc =
critical pressure (i.e., =7.337 MPa [33])
- Pr =
Prandtl number
- q″ =
heat flux (W/cm2)
- q″loss =
heat losses to the environment from the package (W/cm2)
- Rt =
thermal resistance (K·cm2/W)
- R2 =
coefficient of determination
- Re =
Reynolds number
- RMS =
root-mean-square, Table 3
- t =
enhancement factor, Eq. (6)
- T =
temperature (K)
- x =
mass quality
- z =
dimensionless distance from the inlet (e.g., z = l/d)
- @RTD =
property measured at the location of the RTD (K)
Greek Symbols
- ΔT =
temperature difference (K)
- μ =
viscosity (Pa·s)
- ρ =
density (kg/m3)
- σ =
surface tension (N·m)
Subscripts
- amb =
ambient
- b =
bulk fluid
- c =
calculated
- exp =
experimental\measured
- fd =
fully developed
- i =
summation index, Eq. (1)
- in =
at the inlet of the microchannel after the orifice
- Katto =
based on the Katto correlation [24]
- l =
liquid phase
- max =
Maximum
- r =
reduced (e.g., reduced pressure)
- RTD =
property measured at RTD's location
- sat =
saturation
- TP =
two-phase
- v =
vapor phase
- w =
wall
- 1, 2, 3 =
numbers of the RTDs