Empirical Model to Calculate the Thermodynamic Wet-Bulb Temperature of Moist Air ()
1. Introduction
The condensation and vaporization processes are important in different conventional and non-conventional systems of energy, using fossil fuel or renewable energies, and in systems applied to refrigeration, evaporation, air conditioning, humidification, condensation, dehumidification and airing. For example, as in [1] it shows a complete abstract about refrigeration for evaporation under different climate characteristics, and efficiency of several evaporative cooling equipments is indicated, through wet and dry bulb, and its relationship with mass and heat transference processes. These processes can be carried out using solar energy, through the use of different technologies as in [2] . In the industry of electrical production, they have application in the cogeneration systems, throughout compressed air to feed a gas turbine, and in cooling systems of incoming air to the turbines which allows output power increase, and improves thermal efficiency; these efficiencies are obtained through local climate conditions. The calculations began with determination of temperature of wet bulb ambient as in [3] . Also, they have application in the different gas turbine cycles as water injection cycles, steam injection cycles and evaporative cycles with humidification [4] [5] . Design and analysis of all these systems require knowledge of properties corresponding to dry air and steam water mixture from local environment. Calculus methodology of those properties from dry bulb temperature (Tdb), relative humidity (ϕ), local altitude (Z) is known and is shown by ASHRAE [6] . However, to determine thermodynamic wet bulb temperature(Twb), according to the methodology, a numerical implementation is necessary [7] -[11] .
In the present work, an equation to determine wet bulb temperature for mentioned conditions, is presented. It was obtained by processing numerical data of a Meteorological Station located in the region Ciénega of Chapala in Michoacán México, with latitude 20˚0'52.24''N and length 102˚44.37'37.11''W. During the numerical analysis 91,519 data were processed, for each of the variables involved, and they correspond to a typical local year.
2. Methodology
The considered variables for the analysis were dry bulb temperature and relative humidity, in ranges from 3˚C to 35˚C, and 7% to 97% respectively, with local altitude Z = 1526 m [12] . Wet bulb temperature was determined considering data equation and according to methodology described by ASHRAE, throughout a numerical implementation. Data were grouped as a function of relative humidity; the ranges consisted of 5 units, as it can be seen in Figure 1.
For each humidity range Twb vs Tdb were plotted and a linear adjusted equation was obtained, with general form: Twb = ATdb + B. Table 1 shows adjustment equations for each humidity range and the coefficients of determination.
Then, coefficients for A and B were correlated, for each humidity range considered (average of the range) as shown in Figure 2. Thus, Equation (1) was obtained and polynomial for coefficients A and B.
(1)
In Equation (1), relative humidity is considered dimensionless, with Tdb and Twb given in ˚C. Adjusted coefficient values for each polynomial, are shown in Table2
3. Results and Discussion
Equation (1) has in its general form a linear behavior for dry bulb temperature. But coefficients Ai and Bi de-
Figure 1. Percentage distribution for humidity ranges.
Figure 2. Polynomial coefficients for A and B with coefficients of determination.
Table 1. Adjustment equations from Twb and Tdb for each interval humidity.
Table 2. Fit coefficients for proposal equation.
pend on relative humidity, and they are 3th and 4th degree polynomials.
Figure 3 shows the correlation of Tdb, ϕ, Twb for data station and it is compared with the obtained values from Equation (1), for mentioned ranges, with acceptable fit.
In order to evaluate Equation (1), two errors are estimated: normalized error (En) and real error (Er). These are obtained from base data and with calculated data, according to the proposed equation. Normalized error evaluation was obtained using the following equation:. And real error using equation:.
Figure 4 shows the error between the proposed model and wet bulb temperature data from the meteorological station. The maximum normalized error is 5.4% (secondary vertical axis, En) and is equivalent to maximum variation of 0.6˚C as shown in the main vertical axis (Er). Absolute average error is 0.065% (En_av) about meteorological station, representing an average variation of a hundredth degree Celsius from real error (Er_av = 0.011˚C).
Also, Figure 4 shows that normalized error is present in a range of −6 and 3, and to estimate how this error impacts in base data, absolute frequencies were settled, from normalized error.
Table 3 shows absolute frequencies, this is to estimate the data quality from normalized error. This table shows an En range (Xj) according to Figure 4. Absolute frequency is determined by data selection and it is clear that between 1.0 and −1.0 a percentage of 96.4 of the obtained data are present.
Figure 5 shows data from Table 1 plotted, and 88,265 data are between 1.0 and −1.0, which indicates that the
Figure 3. Comparison of Tdb, ϕ, Twb correlation of meteorological data (a); and from proposed model according to Table 1(b).
Figure 4. Comparison between real error and normalized error evaluated with data from meteorological station.
Table 3. Absolute frequency with its corresponding percentage to En data.
difference of real error and normalized error are in an acceptable value (only 3.6% of total data are out of 1.0 and −1.0).
4. Conclusion
A direct equation was obtained in order to predict wet bulb temperature, from relative humidity and dry bulb temperature. Methodology used to obtain the model, was in an empiric way, grouping relative humidity in ranges of 5 units, and for each range wet and dry bulb temperatures were correlated, obtaining a general linear adjust. Also, coefficients A and B were evaluated. Quality of generated data for the model was settled through error analysis. Normalized errors are less than 5.5%. It means a maximum variation of 0.55˚C from data. The normalized error applies only to 3.6% of the data considered in the present work. The model has an acceptable fit. It is valid for a range of temperature and relative humidity of 3˚C to 100˚C and 7% to 97% respectively. The pressure was calculated from the local altitude and the model behavior was not explored at different pressures. The equation obtained represents an important tool that facilitates the analysis and engineering calculations of various important energy processes.
Nomenclature
Fit Coefficients ASHRAE American Society of Heating, Refrigerating and Air-Conditioning Engineers En Normalized Error, % En_av Absolute Normalized Average Error Er Real Error, ˚C Er_av Real Average Error fj Absolute Frequency N North W West R2 Coefficient of Determination
Temperature, ˚C Xj Range, Corresponding to En Values.
Z Local Altitude, m
Φ Relative humidity, %
ϕ Relative Humidity, Dimensionless
Subscript
db Dry bulb DB Database i Component i-th of fit Coefficients j Data from En
n Normalized r Real wb Wet bulb