Falling film evaporation on a tube bundle with plain and enhanced tubes

The complexities of two-phase flow and evaporation on a tube bundle present important problems in the design of heat exchangers and the understanding of the physical phenomena taking place. The development of structured surfaces to enhance boiling heat transfer and thus reduce the size of evaporators adds another level of complexity to the modeling of such heat exchangers. Horizontal falling film evaporators have the potential to be widely used in large refrigeration systems and heat pumps, in the petrochemical industry and for sea water desalination units, but there is a need to improve the understanding of falling film evaporation mechanisms to provide accurate thermal design methods. The characterization of the effect of enhanced surfaces on the boiling phenomena occurring in falling film evaporators is thus expected to increase and optimize the performance of a tube bundle. In this work, the existing LTCM falling film facility was modified and instrumented to perform falling film evaporation measurements on single tube row and a small tube bundle. Four types of tubes were tested including: a plain tube, an enhanced condensing tube (Gewa-C+LW) and two enhanced boiling tubes (Turbo-EDE2 and Gewa-B4) to extend the existing database. The current investigation includes results for two refrigerants, R134a and R236fa, at a saturation temperature of Tsat = 5 °C, liquid film Reynolds numbers ranging from 0 to 3000, at heat fluxes between 20 and 60 kW/m2 in pool boiling and falling film configurations. Measurements of the local heat transfer coefficient were obtained and utilized to improve the current prediction methods. Finally, the understanding of the physical phenomena governing the falling film evaporation of liquid refrigerants has been improved. Furthermore, a method for predicting the onset of dry patch formation has been developed and a local heat transfer prediction method for falling film evaporation based on a large experimental database has been proposed. These represent significant improvements for the design of falling film evaporators.

Evaporation du réfrigérant HFC 407C sur des tubes lisses ou à surface améliorée.

Daniel Favrat, John Richard Thome, Olivier Zürcher

The substitution of CFC refrigerants in refrigeration systems, heat pumps and organic Rankine cycles for heat recovery, requires good methods for predicting heat transfer of substitute fluids. The measurements in the LENI test facility (concentric tubes with water flowing in a countercurrent flow) with HFC 407C, HFC 134a, HCFC 123, HFC 404A and HFC/HCFC 402A provide a new data band for new refrigerants, and allow a coherent comparison with old refrigerants CFC 11, CFC 12, CFC/HCFC 502 and with existing correlations. The existing correlations were found to be inadequate. Because of this work, an improved flow pattern map and flow boiling model were developed, which resulted in a substantial progress in the accurate predict of heat transfer in plain, horizontal tubes for refigerants without oil. The Kattan et al. (1996, [6]) correlation was programmed to calculate and compare predicted heat transfer coefficients to the new HFC 407C test data. The flow pattern map proposed by Kattan et al. (1996) was also programmed and compared to flow regimes observed for HFC 407C. It predicted the HFC 407C flow pattern data quite accurately. The original objective of the HFC 407C measurements was the validation of the Kattan et al. (1996) correlation applied to a zeotropic refrigerant blend. Local flow boiling heat transfer coefficients were measured for HFC 407C evaporating inside a microfin and plain tube. In addition, microfin heat transfer augmentation relative to a plain tube was investigated. The presence of oil in the evaporator had a complex effect on heat transfer coefficients. Local flow boiling heat transfer coefficients were measured for refrigerant HFC 407C ester oil mixtures (Mobil EAL Arctic 68). A new thermodynamic approach for modeling mixtures of zeotropic refrigerant blends and lubricating oils was also developed. A very high accuracy, straight vibrating tube type of density flowmeter was used to measure oil concentrations of flowing HFC 407C oil mixtures. The test covered oil concentration from 0.5-5 wt.% oil. A method was proposed to predict thermodynamic and transport properties of the refrigerant-oil mixtures. In addition, a first empirical approach was proposed for predicting heat transfer for boiling of refrigerant-oil mixtures. Based on this work, the following phenomena important to the problem were identified that require further investigation: • oil retention in the evaporator tubes that was found and clearly proved in the present project. • the formation of foam in the HFC 407C/oil mixture (but not in the previous HFC 134a/oil tests) •the effect of oil on flow patterns, in particular stratified and stratified-wavy flows.

1996

Heat Transfer and Pressure Drop in the Dryout Region of Intube Evaporation with Refrigerant/Lubricant Mixtures.

Daniel Favrat, John Richard Thome, Olivier Zürcher

A summary of the objectives and conclusions reached in this project are: •A litterature search on post dryout heat transfer for pure refrigerants and for refrigerant-oil mixtures was performed and a comprehensive review written (Chapters 9, 10, and 11); •The correct thermodynamic definition of the heat transfer coefficient for evaporation of refrigerant-oil mixtures was presented and then used in the project, i.e. defining the heat transfer coefficient using the local bubble point temperature Tbub of the refrigerant-oil mixture rather than with the saturation temperature of the pure refrigerant Tsat incorrectly used in prior studies (Chapter 3); •A new generalized approach, that is simple to implement, was developed for predicting bubble point temperatures of refrigerant-oil mixtures for miscible oils that can be applied to any refrigerant-oil combination and is accurate over the local oil concentration range from 0-50 wt.% oil range confronted in direct- expansion evaporators with 0-5 wt.% oil in their refrigerant charge (Chapter 4); •A new thermodynamic method for preparation of temperature-enthalpy- vapor quality (T-h-x) curves for refrigerant-oil mixtures for general application was developed and validated against test data (Chapters 5 and 6), including a simple expression for calculating local oil concentrations as a function of vapor quality; •New, accurate thermodynamic methods for preparation of temperature-enthalpy-vapor (T-h-x) curves for the refrigerant blend R-407C and R-407C/oil mixtures were developed for thermal design of evaporators and for the present experimental test, and was validated against measured test data (Chapter 12); •Recommendations were made on how to include oil effects and refrigerant-oil T-h-x curves into evaporator design methods (Chapter 7); •The effects of oil on set-point parameters for control devices were analyzed and some recommendations made (Chapter 8); •An online oil concentration measurement system was perfected using a high accuracy density flowmeter that allowed oil concentrations to be measured accurately and continuously during experiments, such that oil holdup in the heat transfer test sections at high vapor quality could be identified and the correct inlet oil concentrations be measured, and thus utilized to reduce experimental data; a simpler industrial method for applying the measurement system was also developed and presented (Chapter 15); •Experimental results (local heat transfer coefficients, two-phase pressure drops and some two-phase flow patterns) were obtained for R-134a and R-134a/oil mixtures evaporating in plain and microfinned tubes over a wide range of test conditions (Chapters 16 and 17); •Experimental results (local heat transfer coefficients, two-phase pressure drops and some two-phase flow patterns ) were obtained for R-407C and R- 407C/oil mixtures evaporating in plain and microfinned tubes over a wide range of test conditions (Chapters 18 and 19); •The recently proposed flow boiling model and flow pattern map of Kattan-Thome-Favrat, that importantly models heat transfer coefficients based on local flow pattern and predicts the onset of dryout in horizontal tubes and local heat transfer coefficients in partial dryout regimes, was described in detail; the flow pattern map is also useful to designers wishing to design in specific flow pattern regimes and avoid inefficient heat transfer regimes, such as mist flow and stratified flow (Chapters 20 and 21); •The stratified-wavy flow regime model was modified based on the new test data at high vapor quality for pure R-134a and R-407C, which increased the accuracy of the method significantly in this flow regime, now referred to as the “modified” Kattan- Thome-Favrat flow boiling model (Chapter 22); •The modified Kattan- Thome-Favrat flow boiling model and flow pattern map were shown to accurately predict R-407C heat transfer and flow pattern data and hence is recommended for industrial use for designing direct-expansion evaporators with refrigerant blends [the original version was also verified for R-402 and R-404A blends, see Chapter 21] (Chapter 22); •The effect of local refrigerant-oil viscosity on liquid-phase convection in flow boiling heat transfer and on flow patterns was successfully incorporated into the modified Kattan-Thome-Favrat model, establishing the first flow boiling design method that predicts oil effects on heat transfer (with quite reasonable accuracy too), and included a simple method for calculating the values of µref-oil for use in the model (Chapter 23); •The oil effects on microfin heat transfer coefficients were analyzed and heat transfer augmentation ratios were determined for R-134a/oil and R-407C/oil mixtures (Chapter 25); •Plain tube, two-phase pressure drop gradients for R-134a and R-407C at three mass velocities were compared to the Friedel frictional two-phase pressure drop correlation, attaining accurate results and also microfin pressure drop augmentation ratios were determined (Chapter 25); •The ratio of [µoil/µref]mw, derived from the Arrhenius logrithmic viscosity mixing rule and local oil concentration w, was shown to accurately predict the effect of oil on two-phase pressure drop gradients at high vapor qualities, with an empirical exponent m = 0.18355 found for R-134a/oil mixtures (no foaming) and an empirical expression for m for R-407C/oil mixtures (with foaming). This ratio incorporated into the Friedel correlation accurately predicts two-phase pressure drops of evaporating refrigerant-oil flows (Chapter 25).

1996

Contribution to the heat transfer analysis of natural and substitute refrigerants evaporated in a horizontal tube

Olivier Zürcher

The substitution of CFC refrigerants in refrigeration systems, heat pumps and organic Rankine cycles, requests a good knowledge of the heat transfer and pressure drop properties of substitute fluids. A contribution to this international effort is proposed with the study of two hydrofluorocarbon refrigerants (HFC-134a and the zeotropic mixture HFC-407C) and the study of the natural refrigerant ammonia. The HFCs have been experimentally tested on the first test rig developed in the Laboratory of Industrial Energy Systems (LENI) which was substantially modified to cope with ammonia, both in terms of safety requirements and operational conditions. The experimental test section is composed of two concentric tubes, with in tube evaporation of the refrigerant inside the inner tube of counter-current Annular water heating. A new database of the local heat transfer coefficients for the refrigerants HFC-134a, HFC-407C and ammonia together with pressure drop measurements has been collected and used to define and calibrate a new and more general heat transfer model. The two HFC refrigerants have been tested on the test rig developed by Kattan [36]; due to the chemical and the thermophysical properties of ammonia, a new test rig has been designed for the ammonia tests and a new calculation procedure based on the mean temperature measurement on the tube wall has been used. The experimental database covers evaporation tests on plain and on microfin tubes, and with nominal oil concentrations varying from 0 to 5 [%] (by wt.). Detailed modeling is concentrated on the tests on plain tube without oil. Each of the three refrigerants has been evaporated at 4[°C]. The other experimental conditions for the HFCs refrigerants were an inside tube diameter of 10.92[mm], a heat flux range of 2 – 5 [kW/m2] and three mass velocities of G = 100, 200, 300 [kg/(m2s)] visualized as Stratified-Wavy flows (G = 100[kg/(m2s)]), and mainly Annular for the higher mass velocities (G = 200, 300[kg/(m2s)]). The other global conditions for ammonia were an internal tube diameter of 14.00 [mm], a heat flux range of 5 – 70[kW/m2] and eleven mass velocities of G = 10,20,30,40,45,50,55,60,80,120,140[kg/(m2s)], corresponding to Stratified, Stratified-Wavy, Intermittent and Annular flow patterns; complementary tests, varying the mass velocity at constant vapor qualities (x = 20,50,80[%]), have been made in particular to better characterize boundaries between flow patterns. The flow patterns are visualized through glass sections at both ends of the 3.1 [m] long test tube. An extensive review of void fraction models with sensitivity analysis on the actual available flow pattern map has been made. From the 17 correlations reviewed, the models of Taitel-Dukler [71] and Rouhani [79] have been used. Other improvements to the threshold criteria of the Stratified-Wavy to Annular transition, lead to an accurate diabatic map to predict the flow patterns and their transition for very different families of fluids like HFC's and ammonia. A new approach in the prediction of two-phase flow heat transfer has been proposed through the study of each flow pattern separately, according to a new criterion defining the onset of nucleate boiling as a function of the critical convective heat transfer coefficient representative of the location where nucleate boiling might occur. A function based on a pseudo Biot number allows, from the mean heat flux around the tube periphery, to base the model on two different mean heat fluxes applied respectively to the parts of the perimeter in contact with the liquid and the vapor in stratified types of flow. Considering pure convective heat transfer, or mixed convective and nucleate heat transfer, this division allows the use of a common criterion to be applied to each flow pattern. The two-phase flow heat transfer coefficient has finally been obtained as a weighted mean function of the vapor and the liquid heat transfer coefficient with respect to their contact surface with the heated tube. Based on the database of refrigerants HFC-134a and ammonia, the standard deviation of the new heat transfer model is of σ = 27.9[%]. Even if the database showed that the flow conditions were close to, or in the turbulent to laminar flow transition, and even if the major part of the experimental points were purposely obtained close to the various flow pattern transitions, the new model showed very good agreement with the experimental database. Due to the precision of the new flow pattern map and the effectiveness of the onset on nucleate boiling criterion, this new heat transfer model accurately predicts the heat transfer conditions during evaporation. Finally, a new method to model separated flows is proposed, based on partial hydraulic diameters and a mean interface velocity.

EPFL2000

Evaporation du réfrigérant HFC 407C sur des tubes lisses ou à surface améliorée.

Daniel Favrat, John Richard Thome, Olivier Zürcher

1996

Heat Transfer and Pressure Drop in the Dryout Region of Intube Evaporation with Refrigerant/Lubricant Mixtures.

Daniel Favrat, John Richard Thome, Olivier Zürcher

1996

Contribution to the heat transfer analysis of natural and substitute refrigerants evaporated in a horizontal tube

Olivier Zürcher

EPFL2000