Effect of throat diameters of the ejector on the performance of the refrigeration cycle using a two-phase ejector as an expansion device

This paper is a part in a series that reports on the experimental study of the performance of the two-phase ejector expansion refrigeration cycle. In the present study, three two-phase ejectors are used as an expansion device in the refrigeration cycle. The effects of throat diameter of the motive nozzle, on the coefficient of performance, primary mass flow rate of the refrigerant, secondary mass flow rate of the refrigerant, recirculation ratio, average evaporator pressure, compressor pressure ratio, discharge temperature and cooling capacity, which have never before appeared in open literature, are presented. The effects of the heat sink and heat source temperatures on the system performance are also discussed.     

A comparison of critical flow models for estimating two-phase flow of HCFC22 and HFC134a through short tube orifices

An experimental study to investigate the critical flow of refrigerants through short tube orifices has been performed by measuring the mass flowrates and pressure profiles along the short tube orifice. Eight critical flow models have been examined and their results compared with the experimental data for HCFC22 and HFC134a. These models include four homogeneous equilibrium models, two homogeneous frozen models, and two non-homogeneous equilibrium models. The data indicate that the flow was choked when downstream pressures were lower than the saturation pressure corresponding to the upstream temperature. The observed flows through short tube orifices included a lack of equilibrium due to short time of expansion and homogeneous mist flow at the exit plane. These flow trends would be more consistent with the basic assumptions of the homogeneous frozen models. Based on the comparison of the existing critical flow models and experimental data, the homogeneous frozen models showed the best agreement with the measured data except for exit qualities below 0.06.     

基于延迟均衡的CO_2两相引射器模型研究

引射器对跨临界CO_2引射制冷系统性能有极大的影响。本文考虑CO_2两相引射器中存在的非平衡相变、超音速和壅塞等复杂流动现象,构建了CO_2两相引射器的1D分布模型,并采用延迟均衡理论分析喷嘴中的非平衡相变过程。与实验结果比较显示,所建立的延迟均衡模型能够很好的预测引射器的性能。此外,通过与均衡模型的相比显示,在本文所选工况下,延迟均衡模型计算所得的主动流流量比均衡模型预测值低12.39%~25.30%,同时非平衡现象将延缓喷嘴中的膨胀过程,使得喷嘴出口压力比均衡模型预测值高。本文采用所建模型进一步分析了引射器的结构对性能的影响,结果显示在一定的工况下存在最优的混合室直径使得引射系数和升压比都较高;而当混合室直径一定时,较长的混合室有利于提高引射器的升压比。     

Underpressure of vaporization of refrigerant R-134a through a diabatic capillary tube

Experimental test results with R-134a have demonstrated the phenomenon of metastable flow through a diabatic capillary tube, and showed that the metastable flow of the refrigerant exists when the heat transfer is weak between the capillary tube and the suction line flows. Using homogeneous nucleation theory, a correlation predicting the underpressure of vaporization (a characteristic quantity for the metastable flow) for a diabatic capillary tube has been developed based on the experimental results.     

Distributed and non-steady-state modelling of an air cooler

The refrigerant flow inside the coils of a dry expansion plate-finned air cooler can be distinguished into two completely different types: two-phase flow and single-phase flow. The most difficult part of non-steady-state modelling of an air cooler is to describe the liquid and vapour mass transport phenomena occurring in the two-phase flow region, as this determines the boundary position between the two regions and then the superheat temperature, which is in turn the feedback signal of the thermostatic expansion valve. In fact, the mass transport is mainly governed by the momentum exchange between refrigerant liquid and vapour, which is usually called slip-effect. Because the momentum or force equilibrium is so fast compared to the thermal equilibrium, the slip-effect can be considered as a steady-state phenomenon. With this assumption, the mass transport in an air cooler can be described by using a simple propagation equation. The steady-state slip-effect, however, is found by solving the momentum equations for one-dimensional two-phase flow using advanced computer packages such as . This paper presents the derivation of the equations in non-steady-state modelling of an air cooler as well as the results obtained from the model. Because the model is purely distributed, it is applicable to various kinds of tube circuit arrangements of air coolers. The purpose of the model is studying and optimization of non-steady-state behaviour of refrigerating systems with capacity control.     

Flow visualization of supersonic two-phase transcritical flow of CO2 in an ejector of a refrigeration system

Supersonic two-phase flow of CO₂ in an ejector was investigated in flow visualization experiments. The ejector was installed in a transcritical CO₂ ejector-expansion refrigeration system with a convergent primary nozzle and a rectangular mixing chamber. The flow fields in the suction chamber and the mixing chamber of the ejector were visualized by direct photography for various operating conditions. The results showed that the liquid in the primary flow after the primary nozzle exit increased with increasing primary flow and secondary flow pressures. The expansion angle of the primary flow at the nozzle exit decreased with increasing secondary flow pressures. The primary flow was blocked in some cases at the mixing chamber entrance due to the large expansion angle which reduced the entrainment performance. The entrainment ratio was inversely related to the expansion angle. The primary and the secondary flows had a short mixing region in the mixing chamber with the mixed flow quickly becoming uniform.     

Numerical simulation of non-adiabatic capillary tubes considering metastable region. Part I: Mathematical formulation and numerical model

A detailed one-dimensional steady and transient numerical simulation of the thermal and fluid-dynamic behavior of capillary tube–suction line heat exchangers has been carried out. The governing equations (continuity, momentum, energy and entropy) for fluid flows, together with the energy equation in solids, are solved iteratively in a segregated manner. The discretized governing equations in the zones with fluid flow are coupled using a fully implicit step-by-step method. An implicit central difference numerical scheme and a line-by-line solver were used in solids. A special treatment has been implemented in order to consider transitions (subcooled liquid region, metastable liquid region, metastable two-phase region and equilibrium two-phase region). All the flow variables (enthalpies, temperatures, pressures, mass fractions, heat fluxes, etc.) together with the thermophysical and transport properties are evaluated at each point of the grid in which the domain is discretized. The numerical model allows analysis of aspects such as geometry, type of fluid, critical or non-critical flow conditions, metastable regions and transient cases. Comparison of the numerical simulation with experimental data presented in the technical literature will be shown in Part II of the present paper.     

Numerical analysis of evaporation performance in a finned-tube heat exchanger

This study discusses the effects of the heat exchanger type, refrigerant, inner tube configuration, and fin geometry on evaporator performance by adopting updated correlations of EVSIM, a numerical analysis model based on the tube-by-tube method developed by Domanski. The heat exchanger types considered are the cross-counter flow type and cross-parallel flow type. The refrigerants considered for the numerical test as a working fluid are R-134a, R-410A and R-22. For inner tube configuration, enhanced tube and smooth tube cases are considered. For the air side evaporation performance, heat exchangers using plate fins, wavy fins and slit fins are analyzed. Results show that the heat transfer rate of the cross-counter flow type heat exchanger is 3% higher than that of the cross-parallel flow type with R-22. The total heat transfer rate of the evaporator using R-410A is higher than those using R-22 and R-134a, while the total pressure drop of R-410A is lower than those of R-22 and R-134a. The heat transfer rate of the evaporator using enhanced tubes is two times higher than that using smooth tubes, but the pressure drop of the enhanced tube is 45–50% higher than that of the smooth tubes. The evaporation performance of slit fins is superior to that of plate fins by 54%.     

Flow pattern, mass flow rate, pressure distribution, and temperature distribution of two-phase flow of HFC-134a inside short-tube orifices

New experimental data on the influence of short-tube orifice configuration, including diameter, length, length-to-diameter ratio (L/D), and orientation on the flow pattern, mass flow rate, and pressure distribution of HFC-134a inside the short-tube orifice are presented. Short-tube orifice diameters ranging between 0.605 and 1.2 mm with L/D ranging between 1.87 and 33 are used in the experiments. Three different forms of the metastable liquid flow, which are metastable liquid core flow, conical metastable liquid core flow, and full metastable liquid flow are visually observed. The short-tube orifice diameter has a significant effect on the increase in the flow rate. Conversely, the change in the orientation of the test section has no significant effect on the flow rate. The choke flow phenomenon disappears inside the short-tube orifice when L/D is less than 2.91. Based on the present data, a correlation for predicting the mass flow rate through short-tube orifices is proposed.     

Two-phase split of refrigerants at a T-junction

The present study experimentally investigated the two-phase flow split of refrigerants at a T-junction. As geometric parameters, the direction of the inlet or branch tube and the tube diameter ratio of branch to inlet tube were chosen. As inlet flow parameters, the inlet mass flux and quality were varied from 100 to 700 kg m⁻² s⁻¹ and from 0.1 to 0.9, respectively, for the condition of distribution header of a multi-pass evaporator in the general refrigeration system. All experiments were performed for R-22, R-134a, and R-410A. The measured data were compared with the values predicted by the models developed for air–water or steam–water mixture in the literature. We propose a modified model for application to the reduced T-junction and vertical orientation of tubes. Among the geometric parameters, the branch tube direction showed the largest sensitivity to the mass flow rate ratio for the gas phase, while the inlet quality showed the largest sensitivity to the mass flow rate ratio among the inlet flow parameters.     

Numerical analysis and experimental investigation into the performance of a wire-on-tube condenser of a retrofitted refrigerator

The present paper discusses (a) the analysis of a wire-on-tube condenser under varying operating conditions of free convection using FEM, and (b) experimental verification of the performance of two wire-on-tube condensers in a retrofitted domestic refrigerator using refrigerant R-134a. The study is motivated by the desire to investigate if the wire-on-tube condensers used in R-12 based refrigerators could be used in a modified refrigerator using R-134a refrigerant. Experiments were conducted in a climate chamber under controlled and varying ambient temperatures and mass flow rates to determine the locations where phase change occurs and the degree of subcooling achieved. In terms of initial and final phase change point locations the predicted results agree with the experimental results to within ±10%. The analysis and the experiments also lead to the information about the adequacy of the number of tubes for complete condensation of the refrigerant vapour under given operating conditions. The methodology can be used as a design tool for the design of wire-on-tube condenser of a small refrigerator as well as the suitability of specific decommissioned condensers for use in a retrofitted refrigerator. It also indicates that R-12 based refrigerators using wire-on-tube condensers retrofitted with R-134a compressor and refrigerant deserve and warrant further studies for adoption.     

Effects of accumulator heat exchangers on the performance of a refrigeration system

An accumulator heat exchanger (AHX) consists of an accumulator and an inner heat exchanger (IHX) contained in a shell. The AHX has been used in multi-air-conditioners to obtain system reliability and high performance by providing liquid refrigerant into expansion devices and preventing wet-compression. Energy is exchanged between the evaporator exit and the condenser exit in the AHX. In this study, the heat transfer characteristics of the AHX were investigated experimentally, and the effects of the AHX on the performance of a refrigeration system using R22 were measured. The operating characteristics of the refrigeration system with the AHX were considerably different from those without the AHX. The AHX system showed higher refrigerant flow rate than the non-AHX system at a constant EEV (electronic expansion valve) opening because of higher subcooling, resulting in better performance and reliability of the refrigeration system. At 50% EEV opening, the cooling capacity and COP of the AHX system were higher than those of the non-AHX system by 7.5% and 3.2%, respectively.     

Prediction of two-phase pressure gradients of refrigerants in horizontal tubes

Two-phase pressure drop data were obtained for evaporation in two horizontal test sections of 10.92 and 12.00 mm diameter for five refrigerants (R-134a, R-123, R-402A, R-404A and R-502) over mass velocities from 100 to 500 kg/m² s and vapor qualities from 0.04 to 1.0. These data have then been compared against seven two-phase frictional pressure drop prediction methods. Overall, the method by Müller-Steinhagen and Heck (Müller-Steinhagen H, Heck K. A simple friction pressure drop correlation for two-phase flow in pipes. Chem. Eng. Process 1986;20:297–308) and that by Grönnerud (Grönnerud R. Investigation of liquid hold-up, flow-resistance and heat transfer in circulation type evaporators, part IV: two-phase flow resistance in boiling refrigerants. Annexe 1972-1, Bull. de l'Inst. du Froid, 1979) were found to provide the most accurate predictions while the widely quoted method of Friedel (Friedel L. Improved friction drop correlations for horizontal and vertical two-phase pipe flow. European Two-phase Flow Group Meeting, paper E2; June 1979; Ispra, Italy) gave the third best results. The data were also classified by two-phase flow pattern using the Kattan-Thome-Favrat (Kattan N, Thome JR, Favrat D. Flow boiling in horizontal tubes. Part 1: development of a diabatic two-phase flow pattern map. J. Heat Transfer 1998;120:140–7; Kattan N, Thome JR, Favrat D. Flow boiling in horizontal tubes. Part 2; new heat transfer data for five refrigerants. J Heat Transfer 1998;120:148–55; Kattan N, Thome JR, Favrat D. Flow boiling in horizontal tubes. Part 3: development of a new heat transfer model based on flow patterns. J. Heat Transfer 1998;120:156–65) flow pattern map. The best available method for annular flow was that of Müller-Steinhagen and Heck. For intermittent flow and stratified-wavy flow, the best method in both cases was that of Grönnerud. It was observed that the peak in the two-phase frictional pressure gradient at high vapor qualities coincided with the onset of dryout in the annular flow regime.     

Mathematical modeling of falling liquid film evaporation process

A mathematical model of evaporation process from a laminar falling liquid film on a vertical plate of constant temperature is presented. The model is developed with and without interfacial shear stress due to the vapor flow at the liquid film surface. The vapor pressure drop, vapor exit velocity and cooling rate are calculated for different liquid mass flow values. It is shown that lower liquid mass flow produces higher cooling rate. The results also show that the interfacial shear stress has a considerable negative effect on the cooling rate. It is proved that there exists an optimum distance between the plates, which gives the maximum volumetric cooling rate.     

Prediction of condensation in steam ejector for a refrigeration system

An experimental refrigeration system based on a two-stage steam ejector was set-up in the Thermodynamics and Heat Transfer Laboratory of our Department. The system optimization and realization have been described elsewhere ( [Grazzini and Mariani, 1998] and [Grazzini and Rocchetti, 2008] ). In both stages, primary flows are highly supersonic and reach low pressure and temperature levels. As usual in the literature, an ideal gas model was used during the design process. This paper is intended to check the validity of this assumption. In order to understand the actual working condition of our system, several models have been compared. The presence of high flow speed suggests the existence of metastable conditions. To set the border for the metastable region, the spinodal curve has been drawn. Isentropic expansion of vapour through the nozzle, modelled as ideal gas, seems well within the metastable zone. However, the Classic Nucleation Theory shows that the Wilson line is crossed at the nozzle throat. Condensation produces a marked difference in the conditions at the nozzle exit. Results coming from the present analysis will be used in further optimization of the experimental ejector design.     

Evaluation of heat and mass transfer coefficients in a gauze-type structured packing air dehumidifier operating with liquid desiccant

In this paper, high performance packing, namely, structured packing that has good heat and mass transfer characteristics, is proposed for dehumidification of air using liquid desiccants and for regeneration of liquid desiccants. In order to design a structured packing tower for liquid desiccant — air contacting operations, heat and mass transfer coefficients for each phase are required. This paper is concerned with the interface transfer of heat and mass when air is brought into contact with the liquid desiccant solution. A theoretical study of evaluating heat and mass coefficients in an air-desiccant contact system employing three liquid desiccants, namely calcium chloride, lithium chloride, and a mixture of 50% calcium chloride and 50% lithium chloride (called cost effective liquid desiccant, CELD) is investigated. Moreover, air phase transfer coefficients are correlated with flow rates of air and liquid and the temperature of air, whereas liquid phase coefficients are correlated with rates of air and liquid flow, and the temperature and concentration of the liquid. The findings for the three liquid desiccants are compared and discussed.     

Forced convective boiling heat transfer of R-410A in horizontal minichannels

Convective boiling heat transfer experiments were performed in horizontal minichannels with binary mixture refrigerant, R-410A. The test section is made of stainless steel tubes with inner diameters of 1.5 mm and 3.0 mm and with lengths of 1500 mm and 3000 mm, respectively, and is uniformly heated by applying electric current directly to the tubes. Local heat transfer coefficients were obtained for a heat flux range of 10–30 kW m⁻², a mass flux range of 300–600 kg m⁻² s⁻¹, and quality ranges of up to 1.0. The experimental results were mapped on Wang et al.'s (C.C. Wang, C.S. Chiang, D.C. Lu, Visual observation of two-phase flow pattern of R-22, R-134a, and R-407C in a 6.5-mm smooth tube, Experimental, Thermal and Fluid Science 15 (1997) 395–405) and Wojtan et al.'s (L. Wojtan, T. Ursenbacher, J.R. Thome, Investigation of flow boiling in horizontal tubes: part I – a new diabatic two-phase flow pattern map, International Journal of Heat and Mass Transfer 48 (2005) 2955–2969) flow pattern maps to observe the flow regimes. Laminar flow appears in flow with minichannels. A new boiling heat transfer coefficient correlation based on the superposition model for R-410A was developed with 11.20% mean deviation; it showed a good agreement between the measured data and the calculated heat transfer coefficients.     

Prediction of evaporation heat transfer coefficient and pressure drop of refrigerant mixtures in horizontal tubes

A study on the prediction of heat transfer coefficient and pressure drop of refrigerant mixtures is reported. Heat transfer coefficients and pressure drops of prospective mixtures to replace R12 and R22 are predicted on the same cooling capacity basis assuming evaporation in horizontal tubes. Results indicate that nucleate boiling is suppressed at qualities greater than 20% for all mixtures, and evaporation becomes the main heat transfer mechanism. For the same capacity, some mixtures containing R32 and R152a show 8–10% increase in heat transfer coefficients. Some mixtures with large volatility difference exhibit as much as 55% reduction compared to R12 and R22, caused by mass transfer resistance and property degradation due to mixing (32%) and reduced mass flow rates (23%). Other mixtures with moderate volatility difference exhibit 20–30% degradation due mainly to reduced mass flow rates. The overall impact of heat transfer degradation, however, is insignificant if major heat transfer resistance exists in the heat transfer fluid side (air system). If the resistance in the heat transfer fluid side is of the same order of magnitude as that on the refrigerant side (water system), considerable reduction in overall heat transfer coefficient of up to 20% is expected. A study of the effect of uncertainties in transport properties on heat transfer shows that transport properties of liquid affect heat transfer more than other properties. Uncertainty of 10% in transport properties causes a change of less than 6% in heat transfer prediction.     

Assessment of boiling and condensation heat transfer correlations in the modelling of plate heat exchangers

This paper studies refrigeration cycles in which plate heat exchangers are used as either evaporators or condensers. The performance of the cycle is studied by means of a method introduced in previous papers which consists of assessing the goodness of a calculation method by looking at representative variables such as the evaporation or the condensation temperature depending on the case evaluated. This procedure is also used to compare several heat transfer coefficients in the refrigerant side. As in previous works the models of all the cycle components are considered together with the heat exchanger models in such a way that the system of equations they provide is solved by means of a Newton–Raphson algorithm. Calculated and measured values of the evaporation and the condensation temperatures are also compared. The experimental results correspond to the same air-to-water heat pump studied in other papers and they have been obtained by using refrigerants R-22 and R-290.     

Maximum mass flow ratio due to secondary flow choking in an ejector refrigeration system

The occurrence of flow choking in an ejector of an ejector refrigeration system (ERS) was analysed and a model for predicting the maximum flow ratio of the ejector was developed. The multi-parameter equation to calculate the mass flow ratio takes into account the performance of the primary nozzle, the flow entrainment and mixing relating to ejector geometry and operating conditions. We validated the model using the reported experimental data of refrigerant R113, R141b and steam ERS. The present model was shown to provide better accuracy compared with results obtained by applying the existing 1-D ejector theory. We discussed the application of the model and highlighted the significance of the parameters for future work.