Extension of the applicable range of the implicit curve-fitting method for refrigerant thermodynamic properties to critical pressure

The method of implicit curve-fitting and explicit-calculation has been used for fast and stable calculations of thermodynamic properties of subcritical refrigerants. In order to extend that method to the critical pressure, a method of sectional implicit curve-fitting and explicit-calculation for refrigerant thermodynamic properties is introduced in this paper. The whole data range is divided into several subsections. The requirements on the continuity of thermodynamic properties and the first order derivative of thermodynamic properties in the intersection points of subsections are indicated, and the methods to meet the requirements are presented. Quadric equations are constructed instead of curve-fitting when no data can be given. With the source data obtained from REFPROP 7.1, explicit fast calculation formulae for thermodynamic properties of R410A, covering the saturated temperature of 213.15–344.51 K and superheat of 0–65 K, are given as an example. The calculation speeds of the formulae of R410A are more than 7000 times faster than those of REFPROP 7.1 while the total mean relative deviation of the fast calculation formulae from REFPROP 7.1 is only 0.04%.     

Comparison of R-290 and two HFC blends for walk-in refrigeration systems

To help provide a clear understanding of the relative performance potential of HFCs (R-404A and R-410A) as compared to R-290 for walk-in refrigeration systems representing direct expansion commercial refrigeration systems with small charge, an experimental evaluation of the three refrigerants was investigated. To compare the environmental impact of refrigerants over the entire life cycle of fluid and equipment, including power consumption, the life cycle climate performance (LCCP) of the three refrigerants were evaluated based on measured data. The estimated LCCPs at various emission rates indicate that the LCCP of R-290 is always lower than that of R-404A. The LCCP of R-410A is lower than that of R-290 as long as the annual emission is kept below 10%. It was concluded that R-410A has less or equivalent environmental impact as compared to R-290 when safety (toxicity and flammability), environmental impact (climate change), cost and performance (capacity and COP) are considered.     

Modeling and performance analyses of evaporators in frozen-food supermarket display cabinets at low temperatures

This paper presents modeling and experimental analyses of evaporators in “in situ” frozen-food display cabinets at low temperatures in the supermarket industry. Extensive experiments were conducted to measure store and display cabinet relative humidities and temperatures, and pressures, temperatures and mass flow rates of the refrigerant. The mathematical model adopts various empirical correlations of heat transfer coefficients and frost properties in a fin-tube heat exchanger in order to investigate the influence of indoor conditions on the performance of the display cabinets. The model is validated with the experimental data of “in situ” cabinets. The model would be a good guide tool to the design engineers to evaluate the performance of supermarket display cabinet heat exchangers under various store conditions.     

Liquid side heat transfer and pressure drop in finned-tube cooling-coils operated with secondary refrigerants

In this study full-scale experiments with two different conventional cooling-coils aimed for display cabinets were performed. Heat transfer and pressure drop on the liquid side for three different single phase secondary refrigerants were studied and compared to predictions by existing correlations. Predominantly, the laminar flow regime was studied. The results show that when predicting the heat transfer performance on the liquid side of a cooling-coil the Gnielinski correlation for thermally developing flow and uniform wall temperature boundary conditions (T) leads to good agreement for 0.0014 < x^* < 0.017 if 50 < Re < 1700, assuming a new entrance length is formed after each U-bend. In addition, these entrance lengths must also be accounted for, when predicting the pressure drop on the liquid side of the cooling-coil. The uncertainty of measurement can be a problem in this type of investigations and this has been taken into consideration when analysing the results.     

Measurements of the surface tension for R290, R600a and R290/R600a mixture

The surface tensions of R290, R600a and R290/R600a mixture have been measured by the modified differential capillary-rise method. Twenty-two data points for R290 and 21 data points for R600a were obtained in the temperature range between 273 K and 354 K, and 43 data points for R290/R600a mixture on three isotherms of 278 K, 300 K and 320 K were obtained. The experimental uncertainties of temperature and surface tension are estimated to be within 20 mK and 0.2 mN m⁻¹, respectively. Surface tension correlations as a function of temperature for pure R290 and R600a were formulated in the temperature range between 253 K and critical temperature, and the correlation as a function of the composition for R290/R600a mixture was discussed at 278 K, 300 K and 320 K. It is found that the surface tension for R290/R600a mixture can be reproduced by the simple mixing rule by mole fraction with the correlations of both pure components.     

Calculation methods for comparing the performance of pure and mixed working fluids in heat pump applications

Three methods for comparing cycle performance of working fluids, pure as well as non-azeotropic mixtures, are investigated for two applications and for two mixture pairs, HCFC22-CFC114 and HCFC22-HCFC142b, and their pure components. The methods differ in the way of calculating the heat exchange processes. They assume, respectively, equal minimum approach temperatures, equal mean temperature differences and equal heat transfer areas. Changes of coefficient of performance (COP) with composition are explained for all methods. It is shown that transport properties must be taken into account when making rigorous comparisons between working fluids. To predict the relations between fluids with high accuracy, one must use the method with equal heat transfer areas. By the method with equal mean temperature differences, the COP can be estimated with the same accuracy for mixtures as for pure fluids, and can be used for rough estimations of the COP level with different fluids. The method of equal minimum approach temperatures should be avoided for non-azeotropic mixtures.     

Prediction of evaporation heat transfer coefficient and pressure drop of refrigerant mixtures

A study on the prediction of heat transfer coefficient (HTC) and pressure drop of refrigerant mixtures is reported. HTCs and pressure drops of prospective mixtures to replace R12 and R22 are predicted on the same cooling capacity basis. Results indicate that nucleate boiling is suppressed at qualities greater than 20.0% for all mixtures and evaporation becomes the main heat transfer mechanism. For the same capacity, some mixtures containing R32 and R152a show 8.0–10.0% increase in HTCs. Some mixtures with large volatility difference exhibit as much as 55.0% reduction compared with R12 and R22, caused by mass transfer resistance and property degradation due to mixing (32.0%) and reduced mass flow rates (23.0%). Other mixtures with moderate volatility difference exhibit 20.0–30.0% 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 HTC 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.0% in transport properties causes a change of less than 6% in heat transfer prediction.     

An equation of state for 1,1-difluroethane (HFC 152a)

An equation of state for 1,1-difluoroethane (HFC 152a, CH₃CHF₂) has been developed on the basis of reliable experimental data including PVT, liquid C_p, and saturated-liquid-density data measured by our group. It is a non-dimensionalized virial equation whose functional form is the same as that originally developed for 1,1,1,2-tetrafluoroethane (HFC 134a) in our group. The effective range is for pressures up to 15 MPa, temperatures from 230 to 450 K, and densities to 1000 kg m⁻³. The equation represents reliable PVT measurements within ± 1% in pressure for the superheated vapour and supercritical fluid, while within ±0.5% in density for the compressed liquid. In addition, it should be noted that the equation represents the other essential thermodynamic properties including vapour pressures, saturated-liquid/ vapour densities, isobaric/isochoric specific heats and sound velocity in both the liquid and gaseous phase of HFC 152a.     

Thermodynamic properties of HFC-32 (difluoromethane)

An 18-coefficient modified Benedict–Webb–Rubin equation of state of HFC-32 (difluoromethane) has been developed, based on the updated available PVT measurements, heat capacity measurements and speed of sound measurements. Correlations of vapor pressure and saturated liquid density are also presented. The correlations have been developed based on the reported experimental saturation properties data. This equation of state is effective both in the superheated gaseous phase and compressed liquid phase at pressures up to 70 MPa, densities to 1450 kg/m³, and temperatures from 150 to 475 K, respectively.     

A generalized correlation for refrigerant mass flow rate through adiabatic capillary tubes

A capillary tube is a common expansion device used in small sized refrigeration and air-conditioning systems. A generalized correlation for refrigerant flow rate in adiabatic capillary tubes is developed by implementing dimensionless parameters based on extensive experimental data for R-22, R-290, and R-407C measured in this study. Dimensionless parameters are derived from the Buckingham Pi theorem, considering the effects of tube inlet conditions, capillary tube geometry, and refrigerant properties on mass flow rate. The generalized correlation yields good agreement with the present data for R-22, R-290, and R-407C with average and standard deviations of 0.9 and 5.0%, respectively. Approximately 97% of the present data are correlated within a relative deviation of ±10%. Further assessments of the correlation are made by comparing the predictions with measured data for R-12, R-134a, R-152a, R-410A, and R-600a in the open literature. The correlation predicts the data for those five refrigerants with average and standard deviations of −0.73 and 6.16%, respectively.