PVT measurements of water-ammonia refrigerant mixture in the critical and supercritical regions

The PVTx properties of H₂O + NH₃ mixture (0.2607 mole fraction of ammonia) have been measured in the near- and supercritical regions. Measurements were made along 40 liquid and vapor isochores in the range from 120.03 to 727.75 kg m⁻³ and at temperatures from 301 to 634 K and at pressures up to 28 MPa. Temperatures and densities at the liquid–gas phase transition curve, dew- and bubble-pressure points, and the critical parameters for the 0.7393 H₂O + 0.2607 NH₃ mixture were obtained using the quasi-static thermograms and isochoric (P–T) break-point techniques. The expanded uncertainty of the density, pressure, temperature measurements at the 95% confidence level with a coverage factor of k = 2 is estimated to be 0.06%, 0.02–0.09%, and 15 mK, respectively.     

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.     

Experimental researches on characteristics of vapor–liquid equilibrium of NH3–H2O–LiBr system

An improved system of NH₃–H₂O–LiBr was proposed for overcoming the drawback of NH₃–H₂O absorption refrigeration system. The LiBr was added to NH₃–H₂O system anticipating a decrease in the content of water in the NH₃–H₂O–LiBr system. An equilibrium cell was used to measure thermal property of the ternary NH₃–H₂O–LiBr mixtures. The pressure–temperature data for their vapor–liquid equilibrium (VLE) data were measured at ten temperature points between 15–85 °C, and pressures up to 2 MPa. The LiBr concentration of the solution was chosen in the range of 5–60% of mass ratio of LiBr in pure water. The VLE for the NH₃–H₂O–LiBr ternary solution was measured statically. The experimental results show that the equilibrium pressures reduced by 30–50%, and the amount of component of water in the gas phase reduced greatly to 2.5% at T=70 °C. The experimental results predicted much better characteristics of the new ternary system than the NH₃–H₂O system for the applications.     

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.     

Vapor–liquid equilibria of the (hexafluoroethane + 1,1,1-trifluoroethane) binary system from 258 to 343 K up to 3.89 MPa

Development of modern refrigeration systems is critical for the success of new global environmental protection efforts. The binary system of refrigerants: Hexafluoroethane (R116) + 1,1,1-trifluoroethane (R143a), has been studied with the aim of providing PTxy data. Isothermal vapor–liquid equilibrium data have been generated using the “static–analytic” method from 258 to 328 K at pressures from 0.39 to 3.89 MPa. The model composed of the Peng–Robinson equation of state, the Mathias–Copeman alpha function, the Wong–Sandler mixing rules and the NRTL cell theory is applied herein to correlate the data and calculate the critical line.     

Density of organic binary mixtures from equilibrium measurements

A reduction model for equilibrium measurements data (pressure–temperature–weight data) was developed. The model allows indirect measurements of the refrigerant concentration at the liquid and the gas phases and the solution density. The reduction model was validated by both comparing the obtained pressure–temperature–concentration relation at equilibrium conditions and solution densities with experimentally direct measurements carried out by Kriebel and Loffler [M. Kriebel, J. Loffler, Therodynamische eigenschaften des binary systems difluormono-chlormethan (R22) – tetraäthylenglykoldimethyläther (E181), Kältetechnik 17 (1965) 266–271] (up to 80 °C) and by us (up to 120 °C) over a wide range of pressure and temperature.     

Condensation heat transfer of binary refrigerant mixtures of R22 and R114 inside a horizontal tube with internal spiral grooves

An experimental study of the condensation of pure and mixed refrigerants of R22 and R114 inside a spirally grooved horizontal copper tube has been carried out. A double-tube counterflow condenser in the pressure range 3–21 bar and at a mass flow-rate 26–70 kg h⁻¹ was used. The axial distributions of refrigerant, tube wall and cooling water temperatures, wall heat flux density and vapour quality are shown graphically. The variation of tube wall temperature around the circumference of the tube is also shown. The local Nusselt number depends on the molar fraction, whereas the average Nusselt number can be correlated by an equation which is modified from a previously established equation for pure refrigerants inside a horizontal smooth tube. The frictional pressure drop evaluated is correlated well by the Lockhart-Martinelli parameters and is independent of the concentration of the mixture.     

Analysis of the gas compressibility effects on the constant-entropy reversible processes of refrigerants and refrigerant mixtures

Thermally and calorically real gas modelling based on the Martin–Hou equation of state is assumed for pure and mixed refrigerants in the superheated vapour phase. It allows the constant-entropy reversible processes which take place within the work transfer components of ideal vapour compression cycles to be properly analysed. These processes, in fact, occur in a region of the Mollier diagram close to the saturated vapour curve where covolume and molecular forces alter the equation of state of an ideal gas. Thus, real gas effects are significant and cannot be ignored. They give a more accurate indication of the refrigerant end temperature associated with an isentropic compression as well as of the corresponding work exchanged and volumetric efficiency. In particular, it is shown that the gas compressibility effects play a ‘favourable’ role during the isentropic compression processes since they allow the work transferred to be reduced up to 10% for HFC-refrigerant 134a, and HFC-refrigerant mixtures 407C and 410A. But, at the same time, they play an ‘unfavourable’ role since they can reduce the compressor volumetric efficiency (i.e. refrigerant mass flow rate) and, consequently, the cooling (or heating) capacity of the vapour compression system up to 7%.     

Refrigerant R134a condensation heat transfer and pressure drop inside a small brazed plate heat exchanger

This paper presents the experimental tests on HFC-134a condensation inside a small brazed plate heat exchanger: the effects of refrigerant mass flux, saturation temperature and vapour super-heating are investigated.A transition point between gravity controlled and forced convection condensation has been found for a refrigerant mass flux around 20 kg/m² s. For refrigerant mass flux lower than 20 kg/m² s, the saturated vapour heat transfer coefficients are not dependent on mass flux and are well predicted by the Nusselt [Nusselt, W., 1916. Die oberflachenkondensation des wasserdampfes. Z. Ver. Dt. Ing. 60, 541–546, 569–575] analysis for vertical surface. For refrigerant mass flux higher than 20 kg/m² s, the saturated vapour heat transfer coefficients depend on mass flux and are well predicted by the Akers et al. [Akers, W.W., Deans, H.A., Crosser, O.K., 1959. Condensing heat transfer within horizontal tubes. Chem. Eng. Prog. Symp. Ser. 55, 171–176] equation. In the forced convection condensation region, the heat transfer coefficients show a 30% increase for a doubling of the refrigerant mass flux. The condensation heat transfer coefficients of super-heated vapour are 8–10% higher than those of saturated vapour and are well predicted by the Webb [Webb, R.L., 1998. Convective condensation of superheated vapour. ASME J. Heat Transfer 120, 418–421] model. The heat transfer coefficients show weak sensitivity to saturation temperature. The frictional pressure drop shows a linear dependence on the kinetic energy per unit volume of the refrigerant flow and therefore a quadratic dependence on the refrigerant mass flux.     

A theoretical study of a novel regenerative ejector refrigeration cycle

There has been a demand for developments of the ejector refrigeration systems using low grade thermal energy, such as solar energy and waste heat. In this paper, a novel regenerative ejector refrigeration cycle was described, which uses an auxiliary jet pump and a conventional regenerator to enhance the performance of the novel cycle. The theoretical analysis on the performance characteristics was carried out for the novel cycle with the refrigerant R141b. Compared with the conventional cycle, the simulation results show that the coefficient of performance (COP) of the novel cycle increases, respectively, by from 9.3 to 12.1% when generating temperature is in a range of 80–160 °C, the condensing temperature is in a range of 35–45 °C and the evaporating temperature is fixed at 10 °C. Especially due to the enhanced regeneration with increasing the pump outlet pressure, the improvement of COP of the novel cycle is approached to 17.8% compared with that in the conventional cycle under the operating condition that generating temperature is 100 °C, condensing temperature is 40 °C and evaporating temperature is 10 °C. Therefore, the characteristics of the novel cycle performance show its promise in using low grade thermal energy for the ejector refrigeration system.     

A simple set of equations of state for process calculations and its application to R134a and R152a

This paper reports a useful set of equations which enables the consistent and reliable calculation of thermodynamic properties. This set of equations consists of a vapour pressure equation, an equation for the gas phase p, v, T properties, an equation giving the saturated liquid densities and an equation for the specific heat capacity in the ideal gas domain. These equations are of a simple structure because the critical region is excluded. Therefore, for a preliminary investigation only few experimental data points are required for parameter regression, which makes this set of equations suitable for ‘new’ refrigerants. The relationships for enthalpy and entropy are derived from these equations and evaluation procedures are summarized. Examples are given for the refrigerants R134a and R152a.     

Solubility, density and viscosity of a mixture of R-600a and polyol ester oil

Experimental data on solubility, liquid phase density and viscosity of a mixture of R-600a and a POE ISO 7 lubricant oil are presented. A specially designed experimental facility for simultaneous measurements of the physical properties was used in the experiments at temperatures ranging from 10 to 60 °C. The VLE data were correlated with the Heil–Prausnitz and Flory–Huggins activity models and the Peng–Robinson equation of state (EoS). Liquid density was correlated with the Peng–Robinson EoS and with a first-order Redlich–Kister expansion for the excess molar volume. Liquid viscosity was correlated with an excess-property approach based on the classical Eyring liquid viscosity model. Satisfactory agreement was obtained between models and experiments; maximum root mean square (RMS) deviations of models used in the VLE, density and viscosity predictions were 1.1% (VLE EoS), 0.2% (Redlich–Kister) and 3.0% (Grunberg–Nissan), respectively.     

Experimental studies on the evaporative heat transfer and pressure drop of CO2 in smooth and micro-fin tubes of the diameters of 5 and 9.52 mm

Carbon dioxide among natural refrigerants has gained a considerable attention as an alternative refrigerant due to its excellent thermophysical properties. In-tube evaporation heat transfer characteristics of carbon dioxide were experimentally investigated and analyzed as a function of evaporating temperature, mass flux, heat flux and tube geometry. Heat transfer coefficient data during evaporation process of carbon dioxide were measured for 5 m long smooth and micro-fin tubes with outer diameters of 5 and 9.52 mm. The tests were conducted at mass fluxes of from 212 to 656 kg m⁻² s⁻¹, saturation temperatures of from 0 to 20 °C and heat fluxes of from 6 to 20 kW m⁻². The difference of heat transfer characteristics between smooth and micro-fin tubes and the effect of mass flux, heat flux, and evaporation temperature on enhancement factor (EF) and penalty factor (PF) were presented. Average evaporation heat transfer coefficients for a micro-fin tube were approximately 150–200% for 9.52 mm OD tube and 170–210% for 5 mm OD tube higher than those for the smooth tube at the same test conditions. The effect of pressure drop expressed by measured penalty factor of 1.2–1.35 was smaller than that of heat transfer enhancement.     

Supersonic two-phase flow of CO2 through converging–diverging nozzles for the ejector refrigeration cycle

CO₂ is environmentally friendly, safe and more suitable to ejector refrigeration cycle than to vapor compression cycle. Supersonic two-phase flow of CO₂ in the diverging sections of rectangular converging–diverging nozzles was investigated. The divergence angles with significant variation of decompression were 0.076°, 0.153°, 0.306° and 0.612°. This paper presents experimental decompression phenomena which can be used in designing nozzles and an assessment of Isentropic Homogeneous Equilibrium (IHE). Inlet conditions around 6–9 MPa, 20–37 °C were used to resemble ejector nozzles of coolers and heat pumps. For inlet temperature around 37 °C, throat decompression boiling from the saturated liquid line, supersonic decompression and IHE solution were obtained for the two large divergence angles. For divergence angles larger than 0.306°, decompression curves for inlet temperature above 35 °C approached IHE curves. For divergence angles smaller than 0.306° or for nozzles with inlet temperature below 35 °C, IHE had no solution.     

An experimental investigation and modelling of the thermodynamic properties of isobutane–compressor oil solutions: Some aspects of experimental methodology

This paper presents experimental data for the solubility, density and capillary constant for solutions of natural refrigerant isobutane with commercial mineral compressor oil Azmol over a wide range of temperatures and concentrations. Based on information for the capillary constant, the surface tension of the solutions isobutane/Azmol is determined. The experimental data were obtained in the temperature range from 303 K to 363 K and at pressures up to 1.7 MPa using static methods. The experimental data obtained for the solutions of the natural refrigerant isobutane with the commercial mineral compressor oil Azmol are sufficiently described with the help of correlations based on the theory of thermodynamic similarity. The paper reports variation of the vapor pressure, density, capillary constant and surface tension as a function of concentration for the isobutane/Azmol solutions. The enthalpy of liquid phase of the isobutane/Azmol solutions is calculated. The analysis of the behaviour of the excess thermodynamic functions is carried out. The paper examines experimental and methodical uncertainties in the investigation of thermodynamic properties of the refrigerant/oil solutions (ROS). The influence of the time taken to establish thermodynamic equilibrium in the experimental cell on the uncertainty of the experimental data for gas-saturated mixtures such as ROS is discussed. Information about the changing concentration of refrigerant in the liquid phase of the ROS and in the surface layer of the liquid phase of the ROS at increasing temperature is presented. In addition, the experimental data for the density, surface tension and refractive index of the mineral compressor oil Azmol are reported.     

Dynamic simulation of a thermostatically controlled counter-flow evaporatorSimulation dynamique d'un évaporateur à détendeur thermostatique avec débit à contre-courant

A numerical model is developed to simulate the transient behaviour of a counter-flow water cooling evaporator controlled by a thermostatic expansion valve (TEV) in a vapour compression refrigeration system. The liquid–vapour slip in the two-phase region of the evaporator is accounted for by a void fraction model (VFM). The thermal capacitance of the TEV is included in the analysis. For the purpose of comparison with predictions of the model, experimental data available are filtered to obtain the best estimate of the mean variation of the liquid–vapour transition plane. The predictions are in good trend-wise agreement with the filtered experimental data. The results of the transient simulations demonstrate the dependence of the stability of the evaporator–TEV system on the characteristics of the TEV, the thermal capacitance of the bulb, thermal conductance between the bulb and wall and on the nature of the input disturbance.     

Modélisation des données thermodynamiques du mélange ammoniac/eauModelling of the thermodnyamic properties of the ammonia/water mixture

The present work deals with the modeling of the thermodynamic properties of the ammonia/water mixture using the Gibbs free energy function. For the liquid phase, a three constant Margules model of the excess free enthalpy is formulated. The vapour phase is considered as a perfect mixture of real gases, each pure gas being described by a virial equation state in pressure truncated after the third term. The model developed describes with a good accuracy the mixture in the three states, subcooled liquid, superheated vapour and liquid-vapour saturation for temperatures from 200 to 500 K and pressures up to 100 bar.