Flow boiling heat transfer of R134a,R236fa and R245fa in a horizontal 1.030 mm circular channel

This research focuses on acquiring accurate flow boiling heat transfer data and flow pattern visualization for three refrigerants, R134a, R236fa and R245fa in a 1.030 mm channel. We investigate trends in the data, and their possible mechanisms, for mass fluxes from 200 to 1600 kg/m²s, heat fluxes from 2.3 kW/m² to 250 kW/m² at T_sat = 31 °C and ΔT_sub from 2 to 9 K. The local saturated flow boiling heat transfer coefficients display a heat flux and a mass flux dependency but no residual subcooling influence. The changes in heat transfer trends correspond well with flow regime transitions. These were segregated into the isolated bubble (IB) regime, the coalescing bubble (CB) regime, and the annular (A) regime for the three fluids. The importance of nucleate boiling and forced convection in these small channels is still relatively unclear and requires further research.     

Flow boiling of halogenated refrigerants at high saturation temperature in a horizontal smooth tube

Several correlations are available in the open literature for computing the heat transfer coefficient during flow boiling inside plain channels. With respect to halogenated refrigerants, these correlations are usually compared to data taken in a limited range of evaporation temperature and reduced pressure. More recently, the adoption of new refrigerants, such as high pressure HFCs and carbon dioxide, requires to largely extend the pressure range of application of such correlations. Besides, the design of evaporators for some heat pumping applications, where temperatures are set at higher values as compared to usual evaporating temperatures in air-conditioning equipment, requires proper validation of the computing methods.The present paper aims at comparing four well-known predicting models to a new database collected during flow boiling of HCFC (R22) and HFC refrigerants (R134a, R125 and R410A) in a horizontal 8 mm internal diameter tube. This database is characterized by saturation temperature ranging between 25 °C and 45 °C, reduced pressure spanning between 0.19 and 0.53. Mass velocity ranges between 200 and 600 kg m^?2 s^?1 and heat flux between 9 and 53 kW m^?2.Evaporating heat transfer coefficients of halogenated refrigerants at such high temperatures have not been reported in the open literature so far. The discussion of the results will enlighten some similarities with experimental trends presented in the literature for evaporating carbon dioxide.Two models tested here show good prediction capabilities of the present experimental data, but not for all the data sets in the same way. For the purpose of practical use, a simple modification of the correlation by Gungor and Winterton [1] is proposed, showing that this is able to catch the experimental trends of the present database with good agreement.     

Experimental study on condensation heat transfer in vertical minichannels for new refrigerant R1234ze(E) versus R134a and R236fa

Experimental condensation heat transfer data for the new refrigerant R1234ze(E), trans-1,3,3,3-tetrafluoropropene, are presented and compared with refrigerants R134a and R236fa for a vertically aligned, aluminum multi-port tube. Local condensation heat transfer measurements with such a multi-microchannel test section are very challenging due to the large uncertainties related to the heat flux estimation. Presently, a new experimental test facility was designed with a test section to directly measure the wall temperature along a vertically aligned aluminum multi-port tube with rectangular channels of 1.45 mm hydraulic diameter. Then, a new data reduction process was developed to compute the local condensation heat transfer coefficients accounting for the non-uniform distribution of the local heat flux along the channels. The condensation heat transfer coefficients showed the expected decrease as the vapor quality decreased (1.0-0.0) during the condensation process, as the mass velocity decreased (260-50 kg m⁻² s⁻¹) and as the saturation temperature increased (25-70 °C). However, the heat transfer coefficients were not affected by the condensing heat flux (1-62 kW m⁻²) or by the entrance conditions within the tested range. It was found that the heat transfer performance of R1234ze(E) was about 15-25% lower than for R134a but relatively similar to R236fa. The experimental data were then compared with leading prediction methods from the literature for horizontal channels. In general, the agreement was poor, over-predicting the high Nusselt number data and under-predicting the low Nusselt number data, but capturing the mid-range quite well. A modified correlation was developed and yielded a good agreement with the current database for all three fluids over a wide range of operating conditions.     

Two-phase flow operational maps for multi-microchannel evaporators

The current paper presents new operational maps for several different multi-microchannel evaporators, with and without any inlet restrictions (micro-orifices), for the two-phase flow of refrigerants R245fa, R236fa, and R1234ze(E). The test fluids flowed in 67 parallel channels, each having a cross-sectional area of 100 × 100 μm². In order to emulate the power dissipated by active components in a 3D CMOS CPU chip, two aluminium microheaters were sputtered onto the back-side of the test section providing a 0.5 cm² each. Without any inlet restrictions in the micro-evaporator, significant parallel channel flow instabilities, vapor back flow, and flow maldistribution led to high-amplitude and high-frequency temperature and pressure oscillations. Such undesired phenomena were then prevented by placing restrictions at the inlet of each channel. High-speed flow visualization distinguished eight different operating regimes of the two-phase flow depending on the tested operating conditions. Therefore, the preferred operating regimes can be easily traced. In particular, flashing two-phase flow without back flow appeared to be the best operating regime without any flow and temperature instabilities.     

Pressure drop prediction in annular two-phase flow in macroscale tubes and channels

A new prediction method for the frictional pressure drop in annular two-phase flow is presented. This new prediction method focuses on the aerodynamic interaction between the liquid film and the gas core in annular flows, and explicitly takes into account the asymmetric liquid film distribution in the tube cross section induced by the action of gravity in horizontal tubes operated at low mass fluxes. The underlying experimental database contains 6291 data points from the literature with 13 fluid combinations (both single-component saturated fluids such as water, carbon dioxide and refrigerants R12, R22, R134a, R245fa, R410a, R1234ze, and two-component fluids such as water-argon, water-nitrogen, alcohol-argon, water plus alcohol-argon and water-air), vertical and horizontal tubes and annuli with diameters from 3 mm to 25 mm, and both adiabatic and evaporating flow conditions. The new prediction method is very simple to implement and use, is physically based and outperforms existing pressure drop correlations (mean absolute error of 12.9%, and 7 points out of 10 captured to within ±15%).     

Boiling heat transfer,pressure drop,and flow pattern in a horizontal square minichannel

This study experimentally investigated the flow boiling heat transfer, pressure drop, and flow pattern in a horizontal square minichannel with a hydraulic diameter of 2.0 mm, and the effects of mass flux, vapor quality, heat flux, and refrigerant properties on the flow boiling characteristics were clarified. The heat transfer coefficient and pressure drop of R32 and R1234yf were measured in a mass flux range of 50–400 kgm⁻²s⁻¹ at a saturation temperature of 15 °C. The flow pattern of the square minichannel outlet was observed and was classified as plug, wavy, churn, and annular flows. The heat transfer coefficients in the square minichannel were larger than those in the circular minichannel with a similar hydraulic diameter at low mass flux conditions. The heat transfer coefficients of R32 indicated higher values compared with those of R1234yf at same mass flux and qualities. An empirical heat transfer model taking into account the forced convection, nucleate boiling, and thin liquid film evaporation was developed for horizontal square and circular minichannels. The frictional pressure drop of R32 was 1.5–2 times higher than that of R1234yf at same mass flux and vapor quality condition, and the effect of channel shape on the frictional pressure drop was small unlike the boiling heat transfer.     

Flow regime visualization and pressure drops of HFO-1234yf,R-134a and R-410A during downward two-phase flow in vertical return bends

This paper provides a qualitative visual observation of the two-phase flow patterns for HFO-1234yf and R-134a during downward flow in a vertical 6.7 mm inner diameter glass return bend. The different flow regimes observed are: slug, intermittent and annular flows. Bubble and vapor slug dynamical behaviors in downward slug flow are reported for HFO-1234yf. In addition, to determine the perturbation lengths up- and downstream of the return bend, the total pressure drop has been measured at different pressure tap location up- and downstream of the singularity. Furthermore, 285 pressure drop data points measured for two-phase flow of HFO-1234yf, R-134a and R-410A in vertical downward flow return bends are presented. The flow behavior in the return bend, which is subjected to the complex combined actions of gravity and centrifugal force was expressed in terms of the vapor Froude number. This experimental pressure drop database, which is included in the appendix, is compared to four well-known prediction methods available in the literature.     

Flow boiling heat transfer study of R-134a/R-290/R-600a mixture in 9.52 and 12.7 mm smooth horizontal tubes: Experimental investigation

A detailed experimental investigation is carried out to study the flow boiling heat transfer behavior of R-134a/R-290/R-600a (91%/4.068%/4.932% by mass) refrigerant mixture in smooth horizontal tubes of diameter 9.52 and 12.7 mm. The heat transfer coefficients of the mixture are experimentally measured under varied heat flux conditions for stratified flow patterns using a coaxial counter-current heat exchanger test section. The tests are conducted for refrigerant inlet temperatures between ?9 and 5 °C and mass flow rates ranging from 3 to 5 g s^?1. Kattan–Thome–Favrat maps are used to confirm the flow patterns for the tested conditions. The magnitude of the heat transfer coefficient with respect to flow patterns and different mechanisms of boiling are discussed. The heat transfer coefficient of the refrigerant mixture is also compared with that of R-134a for selected working conditions. The significance of nucleate boiling in the overall heat transfer process under these testing conditions is highlighted.     

Experimental and numerical investigation on heat transfer characteristics of supercritical R1234ze(E) cooled in horizontal tubes

The heat transfer characteristics of supercritical fluids in tubes have been considered indispensable for the design and optimization of the heat exchanger and the energy conversion system. Specifically the cooling heat transfer of supercritical R1234ze(E) in horizontal tubes is a promising heat-power conversion technology; however, there is a scarcity of conducted research in available literature. The present work, the first-ever study in this direction, aims to thoroughly investigate the heat transfer characteristics of supercritical R1234ze(E) which is cooled in horizontal tubes. Experimental work was performed to thoroughly explore and inspect the heat transfer characteristics of supercritical R1234ze(E) passing through the tube of 4.12 mm diameter at 4–5 MPa pressure and 240–400 kg/m² s mass flux. Furthermore, the simulation study, supporting the experimental investigation under the same conditions of pressure and mass flux, extended the range of tube diameter up to 9.44 mm. The effects of pressure, mass flux and tube diameter on the heat transfer coefficient were carefully analyzed in the present research work. Based on the simulation results and experimental results, heat transfer correlations were newly developed by separating the region above and below the pseudo-critical temperature. The average absolute deviation between the calculated Nusselt numbers by the numerical correlation and the simulation results was found 2.87%; the average absolute deviation between the calculated Nusselt numbers by the experimental correlation and the experimental results was found 5.3%.     

An investigation of a model of the flow pattern transition mechanism in relation to the identification of annular flow of R134a in a vertical tube using various void fraction models and flow regime maps

In the present study, new experimental data are presented for literature on the prediction of film thickness and identification of flow regime during the co-current downward condensation in a vertical smooth copper tube having an inner diameter of 8.1 mm and a length of 500 mm. R134a and water are used as working fluids in the tube side and annular side of a double tube heat exchanger, respectively. Condensation experiments are done at mass fluxes of 300 and 515 kg m^?2 s^?1. The condensing temperatures are between 40 and 50 °C; heat fluxes are between 12.65 and 66.61 kW m^?2. The average experimental heat transfer coefficient of the refrigerant HFC-134a is calculated by applying an energy balance based on the energy transferred from the test section. A mathematical model by Barnea et al. based on the momentum balance of liquid and vapor phases is used to determine the condensation film thickness of R134a. The comparative film thickness values are determined indirectly using relevant measured data together with various void fraction models and correlations reported in the open literature. The effects of heat flux, mass flux, and condensation temperature on the film thickness and condensation heat transfer coefficient are also discussed for the laminar and turbulent flow conditions. There is a good agreement between the film thickness results obtained from the theoretical model and those obtained from six of 35 void fraction models in the high mass flux region of R134a. In spite of their different valid conditions, six well-known flow regime maps from the literature are found to be predictive for the annular flow conditions in the test tube in spite of their different operating conditions.     

Roughness and surface material effects on nucleate boiling heat transfer from cylindrical surfaces to refrigerants R-134a and R-123

This paper presents results of an experimental investigation carried out to determine the effects of the surface roughness of different materials on nucleate boiling heat transfer of refrigerants R-134a and R-123. Experiments have been performed over cylindrical surfaces of copper, brass and stainless steel. Surfaces have been treated by different methods in order to obtain an average roughness, Ra, varying from 0.03 μm to 10.5 μm. Boiling curves at different reduced pressures have been raised as part of the investigation. The obtained results have shown significant effects of the surface material, with brass being the best performing and stainless steel the worst. Polished surfaces seem to present slightly better performance than the sand paper roughened. Boiling on very rough surfaces presents a peculiar behavior characterized by good thermal performance at low heat fluxes, the performance deteriorating at high heat fluxes with respect to smoother surfaces.     

Flow regimes and two-phase pressure gradient in horizontal straight tubes: Experimental results for HFO-1234yf, R-134a and R-410A

Two-phase flow regime visualizations of HFO-1234yf and R-134a in a 6.70 mm inner diameter glass straight tube have been simultaneous investigated by top and side views with a high speed high resolution camera. No major difference was observed between both refrigerants. HFO-1234yf flow regimes were satisfactorily predicted by the Wojtan et al. [1] flow pattern map. In addition, 819 pressure drop data points measured during two-phase flow of refrigerants HFO-1234yf, R-134a and R-410A in horizontal straight tubes are presented. The tube diameter (D) varies from 7.90 to 10.85 mm. The mass velocity ranges from 187 to 1702 kg m⁻² s⁻¹ and the saturation temperatures from 4.8 °C to 20.7 °C. The results are compared against 10 well-known two-phase frictional pressure drop prediction methods. For the entire database, the best accuracy is given by the method of Müller-Steinhagen and Heck [2] with around 90% of the data predicted within a ±30% error band. An analysis was carried out on the maximum pressure gradient and on the corresponding vapor quality. A statistical analysis for each flow regime was also carried out.     

Saturated critical heat flux in a multi-microchannel heat sink fed by a split flow system

An extensive experimental campaign has been carried out for the measurement of saturated critical heat flux in a multi-microchannel copper heat sink. The heat sink was formed by 29 parallel channels that were 199 μm wide and 756 μm deep. In order to increase the critical heat flux and reduce the two-phase pressure drop, a split flow system was implemented with one central inlet at the middle of the channels and two outlets at either end. The base critical heat flux was measured using three HFC Refrigerants (R134a, R236fa and R245fa) for mass fluxes ranging from 250 to 1500 kg/m² s, inlet subcoolings from ?25 to ?5 K and saturation temperatures from 20 to 50 °C. The parametric effects of mass velocity, saturation temperature and inlet subcooling were investigated. The analysis showed that significantly higher CHF was obtainable with the split flow system (one inlet–two outlets) compared to the single inlet–single outlet system, providing also a much lower pressure drop. Notably several existing predictive methods matched the experimental data quite well and quantitatively predicted the benefit of higher CHF of the split flow.     

Periodic boiling in parallel micro-channels at low vapor quality

Based on experimental investigations the present study evaluates instability and heat transfer phenomenon under condition of periodic flow boiling of water and ethanol in parallel triangular micro-channels. Tests were performed in the range of hydraulic diameter 100–220 μm, mass flux 32–200 kg/m² s, heat flux 120–270 kW/m², vapor quality x = 0.01–0.08. The period between successive events depends on the boiling number and decreases with an increase in the boiling number. The initial film thickness decreases with increasing heat flux. When the liquid film reached the minimum initial film thickness CHF regime occurred. Temporal variations of pressure drop, fluid and heater temperatures were periodic. Oscillation frequency is the same for the pressure drop, for the fluid temperature at the outlet manifold, and for the mean and maximum heater temperature fluctuations. All these fluctuations are in phase. The CHF phenomenon is different from that observed in a single channel of conventional size. A key difference between micro-channel heat sink and single conventional channel is amplification of parallel-channel instability prior to CHF. The dimensionless experimental values of the heat transfer coefficient are presented as the Nusselt number dependence on the Eotvos number and the boiling number.     

Stratified flow model for convective condensation in an inclined tube

Experimental data are reported for condensation of R134a in an 8.38 mm inner diameter smooth tube in inclined orientations with a mass flux of 200 kg/m² s. Under these conditions, the flow is stratified and there is an optimum inclination angle, which leads to the highest heat transfer coefficient. There is a need for a model to better understand and predict the flow behaviour. In this paper, the state of the art of existing models of stratified two-phase flows in inclined tubes is presented, whereafter a new mechanistic model is proposed. The liquid–vapour distribution in the tube is determined by taking into account the gravitational and the capillary forces. The comparison between the experimental data and the model prediction showed a good agreement in terms of heat transfer coefficients and pressure drops. The effect of the interface curvature on the heat transfer coefficient has been quantified and has been found to be significant. The optimum inclination angle is due to a balance between an increase of the void fraction and an increase in the falling liquid film thickness when the tube is inclined downwards. The effect of the mass flux and the vapour quality on the optimum inclination angle has also been studied.     

Circumferential temperature distribution during nucleate pool boiling outside smooth and modified horizontal tubes

In the work an approach to avoid a circumferential temperature distribution existing during nucleate pool boiling on a horizontal cylinder within low heat flux densities is presented. The idea of the approach is local heat transfer enhancement by a porous layer application on a part of the heating surface. An experiment on nucleate pool boiling heat transfer from horizontal cylinders to saturated R141b and water under atmospheric pressure is reported. Experiments have been conducted using stainless steel tubes with the outside diameter between 8 mm and 23 mm with the active length of 250 mm. The outside surface of the tubes was smooth or partially coated with a porous metallic layer. In particular, measurements of inside circumferential temperature distribution have been performed.     

Flow boiling heat transfer characteristics of R123 and R134a in a micro-channel

Flow boiling heat transfer in a single circular micro-channel of 0.19 mm ID has been experimentally investigated with R123 and R134a for various experimental conditions: heat fluxes (10, 15, 20 kW/m²), mass velocities (314, 392, 470 kg/m² s), vapor qualities (0.2–0.85) and different saturation pressures (158, 208 kPa for R123; 900, 1100 kPa for R134a). The heat transfer trends between R123 and R134a are clearly distinguished. Whether nucleate boiling is suppressed at low vapor quality or not determines the heat transfer trend and mechanism in the flow boiling of micro-channels. High convective heat transfer coefficients in the two-phase flow of micro-channels enables nucleate boiling to be suppressed even at low vapor quality, depending on the wall superheat requirement for nucleate boiling. In the case of early suppression of nucleate boiling, specifically R123, heat transfer is dominated by evaporation of thin liquid films around elongated bubbles. In the contrary case, namely R134a, nucleate boiling is dominant heat transfer mechanism until its suppression at high vapor quality and then two-phase forced convection heat transfer becomes dominant. It is similar to the heat transfer characteristic of macro-channels except the enhancement of nucleate boiling and the short forced convection region.     

A self-consistent,physics-based boiling heat transfer modeling framework for use in computational fluid dynamics

Computational Fluid Dynamics (CFD) offers the opportunity to investigate physically and geometrically complex systems with high fidelity. Its applicability to multiphase flow, and particularly boiling heat transfer, is currently limited by the lack of appropriate closure models to describe all relevant phenomena. In this paper, we present an original subcooled flow boiling modeling framework for CFD, which aims at consistently and accurately characterizing the key physics that affect heat transfer at the boiling surface. The new framework introduces a fully mechanistic representation of heterogeneous boiling that improves numerical robustness and reduces sensitivity to closure coefficients. The proposed formulation is inspired by new experimental insight, and significantly extends the existing boiling models by capturing the effects of (i) the microlayer on surface evaporation, (ii) the boiling surface, and (iii) bubbles sliding along the boiling surface. A new statistical treatment of the location and mutual interactions of bubbles on the surface allows for mechanistic prediction of the dry surface area, an important quantity that affects the boiling heat transfer coefficient. This approach lends itself naturally to extension to very high heat fluxes, potentially up to the critical heat flux. An assessment and sensitivity study of the model is presented for a range of mass fluxes (500–1250 kg/m²/s), heat fluxes (100–1600 kW/m²), inlet subcoolings (5, 10, 15 K), and pressures (1, 1.5, 2 bars), demonstrating improved robustness and predictive accuracy at all tested conditions in comparison to traditional heat partitioning approaches, including high heat fluxes, where classic models often fail to converge. Lastly, the framework proposed here should not be viewed as another heat partitioning model, but rather as a general platform that allows incorporation of advanced models for each physical phenomenon considered, leveraging the growing insight generated by modern experimental diagnostics for boiling heat transfer.     

Transient heating effects in high pressure Diesel injector nozzles

The tendency of today’s fuel injection systems to reach injection pressures up to 3000 bar in order to meet forthcoming emission regulations may significantly increase liquid temperatures due to friction heating; this paper identifies numerically the importance of fuel pressurization, phase-change due to cavitation, wall heat transfer and needle valve motion on the fluid heating induced in high pressure Diesel fuel injectors. These parameters affect the nozzle discharge coefficient (C_d), fuel exit temperature, cavitation volume fraction and temperature distribution within the nozzle. Variable fuel properties, being a function of the local pressure and temperature are found necessary in order to simulate accurately the effects of depressurization and heating induced by friction forces. Comparison of CFD predictions against a 0-D thermodynamic model, indicates that although the mean exit temperature increase relative to the initial fuel temperature is proportional to (1 − C_d²) at fixed needle positions, it can significantly deviate from this value when the motion of the needle valve, controlling the opening and closing of the injection process, is taken into consideration. Increasing the inlet pressure from 2000 bar, which is the pressure utilized in today’s fuel systems to 3000 bar, results to significantly increased fluid temperatures above the boiling point of the Diesel fuel components and therefore regions of potential heterogeneous fuel boiling are identified.     

Two-phase pressure drop and flow visualization of FC-72 in a silicon microchannel heat sink

The rapid development of two-phase microfluidic devices has triggered the demand for a detailed understanding of the flow characteristics inside microchannel heat sinks to advance the cooling process of micro-electronics. The present study focuses on the experimental investigation of pressure drop characteristics and flow visualization of a two-phase flow in a silicon microchannel heat sink. The microchannel heat sink consists of a rectangular silicon chip in which 45 rectangular microchannels were chemically etched with a depth of 276 μm, width of 225 μm, and a length of 16 mm. Experiments are carried out for mass fluxes ranging from 341 to 531 kg/m² s and heat fluxes from 60.4 to 130.6 kW/m² using FC-72 as the working fluid. Bubble growth and flow regimes are observed using high speed visualization. Three major flow regimes are identified: bubbly, slug, and annular. The frictional two-phase pressure drop increases with exit quality for a constant mass flux. An assessment of various pressure drop correlations reported in the literature is conducted for validation. A new general correlation is developed to predict the two-phase pressure drop in microchannel heat sinks for five different refrigerants. The experimental pressure drops for laminar-liquid laminar-vapor and laminar-liquid turbulent-vapor flow conditions are predicted by the new correlation with mean absolute errors of 10.4% and 14.5%, respectively.