Subcooled flow boiling and microbubble emission boiling phenomena in a partially heated microchannel

A simultaneous visualization and measurement study has been carried out to investigate subcooled flow boiling and microbubble emission boiling (MEB) phenomena of deionized water in a partially heated Pyrex glass microchannel, having a hydraulic diameter of 155 μm, which was integrated with a Platinum microheater. Effects of mass flux, inlet water subcooling and surface condition of the microheater on subcooled flow boiling in microchannels are investigated. It is found that MEB occurred at high inlet subcoolings and at high heat fluxes, where vapor bubbles collapsed into microbubbles after contacting with the surrounding highly subcooled liquid. In the fully-developed MEB regime where the entire microheater was covered by MEB, the mass flux, the inlet water subcooling and the heater surface condition have only small effects on the boiling curves. The occurrence of MEB in microchannel can remove a large amount of heat flux, as high as 14.41 MW/m² at a mass flux of 883.8 kg/m² s, with only a moderate rise in wall temperature. Therefore, MEB is a very promising method for cooling of microelectronic chips. Heat transfer in the fully-developed MEB in the microchannel is presented, which is compared with existing subcooled flow boiling heat transfer correlations for macrochannels.     

Subcooled flow boiling heat transfer and associated bubble characteristics of FC-72 on a heated micro-pin-finned silicon chip

Experiments are conducted here to investigate subcooled flow boiling heat transfer and associated bubble characteristics of FC-72 on a heated micro-pin-finned silicon chip flush-mounted on the bottom of a horizontal rectangular channel. In the experiments the mass flux is varied from 287 to 431 kg/m² s, coolant inlet subcooling from 2.3 to 4.3 °C, and imposed heat flux from 1 to 10 W/cm². Besides, the silicon chips contain three different geometries of micro-structures, namely, the smooth, pin-finned 200 and pin-finned 100 surfaces. The pin-finned 200 and 100 surfaces, respectively, contain micro-pin-fins of size 200 μm × 200 μm × 70 μm (width × length × height) and 100 μm × 100 μm × 70 μm. The measured data show that the subcooled flow boiling heat transfer coefficient is reduced at increasing inlet liquid subcooling but is little affected by the coolant mass flux. Besides, adding the micro-pin-fin structures to the chip surface can effectively raise the single-phase convection and flow boiling heat transfer coefficients. Moreover, the mean bubble departure diameter and active nucleation site density are reduced for rises in the FC-72 mass flux and inlet liquid subcooling. Increasing coolant mass flux or reducing inlet liquid subcooling results in a higher mean bubble departure frequency. Furthermore, larger bubble departure diameter, higher bubble departure frequency, and higher active nucleation site density are observed as the imposed heat flux is increased. Finally, empirical correlations for the present data for the heat transfer and bubble characteristics in the FC-72 subcooled flow boiling are proposed.     

Heat transfer characteristics during subcooled flow boiling of a kerosene kind hydrocarbon fuel in a 1 mm diameter channel

The subcooled flow boiling heat transfer characteristics of a kerosene kind hydrocarbon fuel were investigated in an electrically heated horizontal tube with an inner diameter of 1.0 mm, in the range of heat flux: 20–1500 kW/m², fluid temperature: 25–400 °C, mass flux: 1260–2160 kg/m² s, and pressure: 0.25–2.5 MPa. It was proposed that nucleate boiling heat transfer mechanism is dominant, as the heat transfer performance is dependent on heat flux imposed on the channel, rather than the fuel flow rate. It was found that the wall temperatures along the test section kept constant during the fully developed subcooled boiling (FDSB) of the non-azeotropic hydrocarbon fuel. After the onset of nucleate boiling, the temperature differences between inner wall and bulk fluid begin to decrease with the increase of heat flux. Experimental results show that the complicated boiling heat transfer behavior of hydrocarbon fuel is profoundly affected by the pressure and heat flux, especially by fuel subcooling. A correlation of heat transfer coefficients varying with heat fluxes and fuel subcooling was curve fitted. Excellent agreement is obtained between the predicted values and the experimental data.     

Correlation of critical heat flux and two-phase friction factor for subcooled convective boiling in structured surface microchannels

Critical heat flux (CHF) and pressure drop of subcooled flow boiling are measured for a microchannel heat sink containing 75 parallel 100 μm × 200 μm structured surface channels. The heated surface is made of a Cu metal sheet with/without 2 μm thickness diamond film. Tests and measurements are conducted with de-ionized water, de-ionized water +1 vol.% MCNT additive solution, and FC-72 fluids over a mass velocity range of 820–1600 kg/m² s, with inlet temperatures of 15(8.6)°C, 25(13.6)°C, 44(24.6)°C, and 64(36.6)°C for DI water (FC-72), and heat fluxes up to 600 W/cm². The CHF of subcooled flow boiling of the test fluids in the microchannels is measured parametrically. The two-phase pressure drop is also measured. Both CHF and the two-phase friction factor correlation for one-side heating with two other side-structured surface microchannels are proposed and developed in terms of the relevant parameters.     

Subcooled and saturated boiling of ammonia–lithium nitrate solution in a plate-type generator for absorption machines

This paper presents the experimental heat transfer evaluation during subcooled and saturated boiling of ammonia–lithium nitrate solution in a fusion plate heat exchanger, acting as a vapor generator under operating conditions representative of single-effect absorption machines. The solution flow rate and outlet temperature were modified in the ranges of 0.041–0.083 kg/s and 78–97 °C, respectively. The region where vapor bubbles begin to arise is estimated using a correlation for the wall superheat required for the onset of nucleate boiling. Results show that subcooled boiling is present in the generator. The initial boiling temperature is about 3.1 °C lower than the saturation temperature. The influence of the heat and mass fluxes on the boiling heat transfer coefficient is analyzed. The paper offers a correlation for the Nusselt number, including the subcooled and saturated boiling regions.     

Unstable and stable flow boiling in parallel microchannels and in a single microchannel

A simultaneous visualization and measurement study have been carried out to investigate flow boiling instabilities of water in microchannels at various heat fluxes and mass fluxes. Two separate flow boiling experiments were conducted in eight parallel silicon microchannels (with flow interaction from neighboring channels at headers) and in a single microchannel (without flow interaction), respectively. These microchannels, at a length of 30 mm, had an identical trapezoidal cross-section with a hydraulic diameter of 186 μm. At a given heat flux and inlet water temperature, it was found that stable and unstable flow boiling regimes existed, depending on the mass flux. A flow boiling map, in terms of heat flux vs mass flux, showing stable flow boiling regime and unstable flow boiling regime is presented for parallel microchannels as well as for a single microchannel, respectively, at an inlet water temperature of 35 °C. In the stable flow boiling regime, isolated bubbles were generated and were pushed away by the incoming subcooled liquid. Two unstable flow boiling regimes, with long-period oscillation (more than 1 s) and short-period oscillation (less than 0.1 s) in temperature and pressure, were identified. The former was due to the expansion of vapor bubble from downstream while the latter was owing to the flow pattern transition from annular to mist flow. A comparison of results of flow boiling in parallel microchannels and in a single microchannel shows that flow interaction effects from neighboring channels at the headers are significant.     

Enhanced boiling of HFE-7100 dielectric liquid on porous graphite

Enhanced boiling of HFE-7100 dielectric liquid on porous graphite measuring 10 mm × 10 mm is investigated, and results are compared with those for smooth copper (Cu) of the same dimensions. Although liquid is out-gassed for hours before performing the pool boiling experiments, air entrapped in re-entrant type cavities, ranging in size from tens to hundreds of microns, not only enhanced the nucleate boiling heat transfer and the critical heat flux (CHF), but also, the mixing by the released tiny air bubbles from the porous graphite prior to boiling incipience enhanced the natural convection heat transfer by ∼19%. No temperature excursion is associated with the nucleate boiling on porous graphite, which ensues at very low surface superheat of 0.5–0.8 K. Conversely, the temperature overshoot at incipient boiling on Cu is as much as 39.2, 36.6, 34.1 and 32.8 K in 0 (saturation), 10, 20 and 30 K subcooled boiling, respectively. Nucleate boiling ensues on Cu at a surface superheat of 11.9, 10.9, 9.5 and 7.5 K in 0 (saturation), 10, 20 and 30 K subcooled boiling, respectively. The saturation nucleate boiling heat flux on porous graphite is 1700% higher than that on Cu at a surface superheat of ∼10 K and decreases exponentially with increased superheat to ∼60% higher near CHF. The CHF values of HFE-7100 on porous graphite of 31.8, 45.1, 55.9 and 66.4 W/cm² in 0 (saturation), 10, 20 and 30 K subcooled boiling, are 60% higher and the corresponding superheats are 25% lower than those on Cu. In addition, the rate of increase in CHF with increased liquid subcooling is 50% higher than that on Cu.     

Microchannel size effects on local flow boiling heat transfer to a dielectric fluid

Heat transfer with liquid–vapor phase change in microchannels can support very high heat fluxes for use in applications such as the thermal management of high-performance electronics. However, the effects of channel cross-sectional dimensions on the two-phase heat transfer coefficient and pressure drop have not been investigated extensively. In the present work, experiments are conducted to investigate the local flow boiling heat transfer of a dielectric fluid, Fluorinert FC-77, in microchannel heat sinks. Experiments are performed for mass fluxes ranging from 250 to 1600 kg/m² s. Seven different test pieces made from silicon and consisting of parallel microchannels with nominal widths ranging from 100 to 5850 μm, all with a nominal depth of 400 μm, are considered. An array of temperature sensors on the substrate allows for resolution of local temperatures and heat transfer coefficients. The results of this study show that for microchannels of width 400 μm and greater, the heat transfer coefficients corresponding to a fixed wall heat flux as well as the boiling curves are independent of channel size. Also, heat transfer coefficients and boiling curves are independent of mass flux in the nucleate boiling region for a fixed channel size, but are affected by mass flux as convective boiling dominates. A strong dependence of pressure drop on both channel size and mass flux is observed. The experimental results are compared to predictions from a number of existing correlations for both pool boiling and flow boiling heat transfer.     

Experimental investigation and modelling of heat transfer during convective boiling in a minichannel

Flow boiling heat transfer characteristics of water are experimentally studied in a circular minichannel with an inner diameter of 1500 μm. The fluid flows upwards and the test section, made of the nickel alloy Inconel 600, is directly electrically heated. Thus, the evaporation takes place under the defined boundary condition of constant heat flux. Mass fluxes between 50 and 100 kg/(m² s) and heat fluxes from 10 to 115 kW/m² at an inlet pressure of 3 bar are examined.Infrared thermography is applied to measure the outer wall temperatures of the minichannel. This experimental method permits the identification of different boiling regions, boiling mechanisms and the determination of local heat transfer coefficients. Measurements are carried out in single-phase flow, subcooled and saturated boiling regions. The experimental heat transfer coefficients in the region of saturated boiling are compared with correlations available in literature and with a physically founded model developed for convective boiling.     

On the rise paths of single vapor bubbles after the departure from nucleation sites in subcooled upflow boiling

A photographic study was carried out for the subcooled flow boiling of water to elucidate the rise characteristics of single vapor bubbles after the departure from nucleation sites. The test section was a transparent glass tube of 20 mm in inside diameter and the flow direction was vertical upward; liquid subcooling was parametrically changed within 0–16 K keeping system pressure and liquid velocity at 120 kPa and 1 m/s, respectively. The bubble rise paths were analyzed from the video images that were obtained at the heat flux slightly higher than the minimum heat flux for the onset of nucleate boiling. In the present experiments, all the bubbles departed from their nucleation sites immediately after the inception. In low subcooling experiments, bubbles slid upward and consequently were not detached from the vertical heated wall; the bubble size was increased monotonously with time in this case. In moderate and high subcooling experiments, bubbles were detached from the wall after sliding for several millimeters and migrated towards the subcooled bulk liquid. The bubbles then reversed the direction of lateral migration and were reattached to the wall at moderate subcooling while they collapsed due to the condensation at high subcooling. It was hence considered that the mechanisms of the heat transfer from heated wall and the axial growth of vapor volume were influenced by the difference in bubble rise path. It was observed after the inception that bubbles were varied from flattened to more rounded shape. This observation suggested that the bubble detachment is mainly caused by the change in bubble shape due to the surface tension force.     

Measurement of surface dryout near heating surface at high heat fluxes in subcooled pool boiling

The authors have conducted measurements of liquid–vapor behavior in the vicinity of a heating surface for saturated and subcooled pool boiling on an upward-facing copper surface by using a conductance probe method. A previous paper [A. Ono, H. Sakashita, Liquid–vapor structure near heating surface at high heat flux in subcooled pool boiling, Int. J. Heat Mass Transfer 50 (2007) 3481–3489] reported that thicknesses of a liquid rich layer (a so-called macrolayer) forming in subcooled boiling are comparable to or thicker than those formed near the critical heat flux (CHF) in saturated boiling. This paper examines the dryout behavior of the heating surface by utilizing the feature that a thin conductance probe placed very close to the heating surface can detect the formation and dryout of the macrolayer. It was found that the dryout of the macrolayer formed beneath a vapor mass occurs in the latter half of the hovering period of the vapor mass. Two-dimensional measurements conducted at 121 grid points in a 1-mm × 1-mm area at the center of the heating surface showed that the dryout commences at specific areas and spreads over the heating surface as the heat flux approaches the CHF. Furthermore, transient measurements of wall void fractions from nucleate boiling to transition boiling were conducted under the transient heating mode, showing that the wall void fraction has small values (<10%) in the nucleate boiling region, and then steeply increases in the transition boiling region. These findings strongly suggest that the macrolayer dryout model is the most appropriate model of the CHF for saturated and subcooled pool boiling of water on upward facing copper surfaces.     

Correlation of critical heat flux for flow boiling of water in mini-channels

In view of practical significance of a correlation of critical heat flux (CHF) in the aspects of engineering design and prediction, this study is aiming at evaluation of existing CHF correlations for flow boiling of water with available databases taken from small-diameter tubes, and then development of a new, simple CHF correlation. Available CHF databases in the literature for flow boiling of water in small-diameter tubes (0.33 < D_h < 6.22 mm) are collected, covering wide parametric ranges. Three correlations by Bowring, Katto and Shah are evaluated with the CHF data for saturated flow boiling, and three correlations by Inasaka–Nariai, Celata et al. and Hall–Mudawar evaluated with the CHF data for subcooled flow boiling. The Hall–Mudawar correlation and the Shah correlation seem to be the most reliable tools for CHF prediction in the subcooled and saturated flow boiling regions, respectively. In order to avoid the defect of predictive discontinuities often encountered when applying previous correlations, a simple, nondimensional, inlet conditions dependent CHF correlation for saturated flow boiling has been formulated. Its functional form is determined by the application of the artificial neural network and parametric trend analyses to the collected database. Superiority of this correlation has been verified by the database. The new correlation has a mean deviation of 16.8% for this collected databank, smallest among all tested correlations. Compared to many inordinately complex correlations, this new correlation consists only of a single equation.     

Non-linear dynamical analyses of transient surface temperature fluctuations during subcooled pool boiling on a horizontal disk

The class of dynamics in pool boiling on a large-size heater is assessed under subcooled pool boiling conditions. Transient surface temperature measurements are obtained using surface micro-machined K-type thin film thermocouples (TFT) in 10 °C subcooled pool boiling experiments on a 62.23 mm diameter silicon wafer using PF-5060 as the test liquid. Surface temperature data is obtained at each steady state condition to generate the boiling curve. The fraction of false-nearest neighbors, recurrence plots and space–time separation plots are obtained using the TISEAN package. The correlation dimension is then estimated from the re-constructed phase space data using a naïve algorithm. The correlation dimension varies from ～11.2 to 11.5 near onset of nucleate boiling (ONB), to ～7–10 in fully developed nucleate boiling (FDNB) ～7–9 near critical heat flux (CHF) condition, and from ～6.6 to 7.7 in film boiling. False-nearest neighbor estimates and recurrence plots show that nucleate boiling may be dominated by statistical processes near ONB and in partial nucleate boiling (PNB). In contrast, FDNB, CHF and film boiling seem chaotic and governed by deterministic processes.     

Effects of water subcooling on heat transfer in vertical annuli

To obtain effects of major geometric parameters on subcooled heat transfer in vertical annuli three gap sizes (7.05, 18.15, and 28.2 mm) and two bottom conditions (open or closed) have been investigated experimentally in water at atmospheric pressure. Up to 50 °C of liquid subcooling has been tested and 429 data points were obtained at both single phase and boiling regions. The increase in pool subcooling results in much change in heat transfer coefficients. The governing mechanisms are suggested as single-phase natural convection and liquid agitation for the annuli with open bottoms while liquid agitation and bubble coalescence are the major factors at the bottom-closed annuli. Four empirical correlations for the heat transfer coefficient have been suggested in terms of the gap size, the degree of subcooling, and the heat flux. The correlations predict the heat transfer data of single phase and of boiling within ±5% and ±20%, respectively.     

Boiling flow characteristics in microchannels with very hydrophobic surface to super-hydrophilic surface

Boiling flow process plays a very important role to affect the heat transfer in a microchannel. Different boiling flow modes have been found in the past which leads to different oscillations in temperatures and pressures. However, a very important issue, i.e. the surface wettability effects on the boiling flow modes, has never been discussed. The current experiments fabricated three different microchannels with identical sizes at 105 × 1000 × 30000 μm but at different wettability. The microchannels were made by plasma etching a trench on a silicon wafer. The surface made by the plasma etch process is hydrophilic and has a contact angle of 36° when measured by dipping a water droplet on the surface. The surface can be made hydrophobic by coating a thin layer of low surface energy material and has a contact angle of 103° after the coating. In addition, a vapor–liquid–solid growth process was adopted to grow nanowire arrays on the wafer so that the surface becomes super-hydrophilic with a contact angle close to 0°. Different boiling flow patterns on a surface with different wettability were found, which leads to large difference in temperature oscillations. Periodic oscillation in temperatures was not found in both the hydrophobic and the super-hydrophilic surface. During the experiments, the heat flux imposed on the wall varies from 230 to 354.9 kW/m² and the flow of mass flux into the channel from 50 to 583 kg/m²s. Detailed flow regimes in terms of heat flux versus mass flux are also obtained.     

Refrigerant flow boiling heat transfer in parallel microchannels as a function of local vapor quality

Flow boiling of refrigerant HFC-134a in a multi-microchannel copper cold plate evaporator is investigated. The heat transfer coefficient is measured locally for the entire range of vapor qualities starting from subcooled liquid to superheated vapor. The test piece contains 17 parallel, rectangular microchannels (0.762 mm wide) of hydraulic diameter 1.09 mm and aspect ratio 2.5. The design of the test facility is validated by a robust energy balance as well as a comparison of single-phase heat transfer coefficients with results from the literature. Results are presented for four different mass fluxes of 20.3, 40.5, 60.8, and 81.0 kg m^?2 s^?1, which correspond to refrigerant mass flow rates of 0.5–2.0 g s^?1, and at three different pressures 400, 550 and 750 kPa corresponding to saturation temperatures of 8.9, 18.7, and 29 °C. The wall heat flux varies from 0 to 20 W/cm² in the experiments. The heat transfer coefficient is found to vary significantly with refrigerant inlet quality and mass flow rate, but only slightly with saturation pressure for the range of values investigated. The peak heat transfer coefficient is observed for a vapor quality of approximately 20%.     

Drift-flux model in a sub-channel of rod bundle geometry

In view of the practical interest of the drift-flux model for two-phase flow analysis, the distribution parameter and drift velocity constitutive equations have been obtained for subcooled boiling flow in a sub-channel of rod bundle geometry. The constitutive equation of the distribution parameter for subcooled boiling flow in a sub-channel is obtained from the bubble-layer thickness model. In this derivation an existing constitutive equation for subcooled boiling flow in a round pipe is modified by taking account of the difference in the flow channel geometry between the sub-channel and round pipe. The constitutive equation of the drift velocity is proposed based on an existing correlation and considering the rod wall and sub-channel geometry effects. The prediction accuracy of the newly developed correlations has been checked against experimental data in a 3 × 3 rod bundle sub-channel, obtaining better predicting errors than the existing correlations most used in literature.     

Local heat transfer distribution and effect of instabilities during flow boiling in a silicon microchannel heat sink

Flow boiling of the perfluorinated dielectric fluid FC-77 in a silicon microchannel heat sink is investigated. The heat sink contains 60 parallel microchannels each of 100 μm width and 389 μm depth. Twenty-five evenly distributed temperature sensors in the substrate yield local heat transfer coefficients. The pressure drop across the channels is also measured. Experiments are conducted at five flow rates through the heat sink in the range of 20–80 ml/min with the inlet subcooling held at 26 K in all the tests. At each flow rate, the uniform heat input to the substrate is increased in steps so that the fluid experiences flow regimes from single-phase liquid flow to the occurrence of critical heat flux (CHF). In the upstream region of the channels, the flow develops from single-phase liquid flow at low heat fluxes to pulsating two-phase flow at high heat fluxes during flow instability that commences at a threshold heat flux in the range of 30.5–62.3 W/cm² depending on the flow rate. In the downstream region, progressive flow patterns from bubbly flow, slug flow, elongated bubbles or annular flow, alternating wispy-annular and churn flow, and wall dryout at highest heat fluxes are observed. As a result, the heat transfer coefficients in the downstream region experience substantial variations over the entire heat flux range, based on which five distinct boiling regimes are identified. In contrast, the heat transfer coefficient midway along the channels remains relatively constant over the heat flux range tested. Due to changes in flow patterns during flow instability, the heat transfer is enhanced both in the downstream region (prior to extended wall dryout) and in the upstream region. A previous study by the authors found no effect of instabilities during flow boiling in a heat sink with larger microchannels (each 300 μm wide and 389 μm deep); it appears therefore that the effect of instabilities on heat transfer is amplified in smaller-sized channels. While CHF increases with increasing flow rate, the pressure drop across the channels has only a minimal dependence on flow rate once boiling is initiated in the microchannels, and varies almost linearly with increasing heat flux.     

New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside horizontal tubes

A new flow boiling heat transfer model and a new flow pattern map based on the flow boiling heat transfer mechanisms for horizontal tubes have been developed specifically for CO₂. Firstly, a nucleate boiling heat transfer correlation incorporating the effects of reduced pressure and heat flux at low vapor qualities has been proposed for CO₂. Secondly, a nucleate boiling heat transfer suppression factor correlation incorporating liquid film thickness and tube diameters has been proposed based on the flow boiling heat transfer mechanisms so as to capture the trends in the flow boiling heat transfer data. In addition, a dryout inception correlation has been developed. Accordingly, the heat transfer correlation in the dryout region has been modified. In the new flow pattern map, an intermittent flow to annular flow transition criterion and an annular flow to dryout region transition criterion have been proposed based on the changes in the flow boiling heat transfer trends. The flow boiling heat transfer model predicts 75.5% of all the CO₂ database within ±30%. The flow boiling heat transfer model and the flow pattern map are applicable to a wide range of conditions: tube diameters (equivalent diameters for non-circular channels) from 0.8 to 10 mm, mass velocities from 170 to 570 kg/m² s, heat fluxes from 5 to 32 kW/m² and saturation temperatures from −28 to 25 °C (reduced pressures from 0.21 to 0.87).     

A novel time strip flow visualisation technique for investigation of intermittent dewetting and dryout in elongated bubble flow in a microchannel evaporator

A flow visualisation study of flow boiling of R245fa in silicon multi-microchannels at low mass flux and moderate heat flux has been carried out with a high speed digital camera. The micro-evaporator had 67 channels of length 20 mm, width 223 μm, and height 680 μm while the fin width between adjacent channels was 80 μm. The base heat flux ranged from 2 to 26 W cm^?2 for a mass velocity of 100 kg s^?1 m^?2, resulting in exit vapour qualities ranging from 10% to 70%. In particular, a novel time strip technique was developed to analyse the recorded image sequences and significantly highlight the various phenomena occurring along given channels. Notably, this technique was able to reveal profound details regarding the intermittent dryout mechanism of liquid films trapped between the elongated bubbles and the heated channel walls. The results show that the intermittent dryout of the evaporating liquid film is comprised of four stages with distinct time scales and dynamics: (i) the growth of liquid film thinning perturbations to a critical amplitude causing the rupture of the metastable liquid film, (ii) a dewetting stage involving expanding dry spots leading to a rivulet flow regime, (iii) evaporation of the rivulets leading to full dryout, and (iv) a rewetting stage. This intermittent dryout mechanism appears to explain the many seemingly contradictory heat transfer coefficient trends observed with changes in vapour quality in microchannels, thus resolving an important heat transfer dilemma. Furthermore, since dryout is an undesirable event during the practical application of a microchannel evaporator, it is important to delay or even suppress the initial rupture of the liquid film that leads to dryout. This can be achieved by manufacturing or treating the channel surfaces to be highly wettable with the chosen refrigerant.