psychrometrics

Install and load psychrometrics Toolbox for Scilab

1. Download v0.1.3.zip file from Releases.

2. Set Scilab current directory to the download folder.

3. To install and load the toolbox, type:

atomsInstall('v0.1.3.zip'),..
atomsLoad('psychrometrics')

Citation of psychrometrics Toolbox for Scilab

You can cite all versions (both released and pre-released), by using DOI 105281/zenodo.6927994.

This DOI represents all versions, and will always resolve to the latest one.

To cite the last released version, please check https://zenodo.org/account/settings/github/repository/aumpierre-unb/Psychrometrics-for-Scilab.

---

The following is a very short introduction to the steady internal flow of an incompressible and inviscid fluid and to the Psychrometrics Toolbox for Scilab.

Internal flow is a pretty extensive topic in fluid mechanics and there are a lot of important and interesting observations related to it that are not taken into account in this text, because they have no direct impact the computation performed by the functions in this toolbox. Our focus here is a small set of equations that described the phenomenon and are required to solve problems on Psychrometrics.

This text is divided in two main sections: The Theory and The Psychrometrics Toolbox for Scilab.

The Theory

Psychrometrics

Psychrometrics is the field of physics concerned with the thermodynamic properties of a mixture of a condensable vapor and a non condensable gas. The most usual situation in engineering is the mixture of water vapor in the atmospheric air at atmospheric pressure.

A psychrometric chart is a two dimensional diagram used for engineers to predict the thermodynamical state of humid air. In most cases, engineers are concerned with the properties of mixtures of water vapor and atmospheric air at atmospheric pressure. The prediction of the properties of such systems are mandatory for a series of industrial processes related to humidification and dehumidification as well as with air-conditioning processes.

The basic readings from a psychrometric chart are the thermodynamic temperature and the humidity. Humidity is the mass ratio of water vapor to dry air.

Alongside with thermochemical properties of water vapor and dry air, the material and energy balances allow to calculate the specific enthalpy and the specific volume of the system using some equation of state. The ideal gas equations of state is usually used since water vapor and dry air present negligible deviation from ideality at room temperature.

Also, alongside with mass and energy transfer coefficients of water vapor, the combined mass and energy transfer phenomena allow to calculate the temperature of a thin layer of liquid water providing water vapor to the gaseous mixture. This temperature is refereed to as wet bulb temperature, while the thermodynamic temperature is refereed to as dry bulb temperature.

Most psychrometric charts show sets of lines of constant specific volume, constant specific enthalpy and constant wet bulb temperature. Set all together, they produce a fairly complete chart of the thermodynamic state of humid air.

Saturation & Dew Point

Consider a constant pressure control volume filled with humid air only. As temperature decreases, all particles in the system lose energy. Eventually, some particles of water will lose energy ate the point they condensate. At this point, the gaseous mixture contains the maximum possible amount of water particles, it is said to be saturated. In psychrometrics, saturation is the condition where the maximum amount of water vapor is in the gaseous phase.

Once pressure is an indirect measure of the number of particles in the system, the amount of particles in the gaseous phase is indirectly refereed to as saturation pressure.

The thermodynamic state where the smallest amount of energy removed from the gaseous phase produces an incipient condensed phase is called dew point. Dew point is characterized by the dew point temperature and the saturation pressure.

Humidity & Relative Humidity

Humidity is the mass ratio of water vapor and dry air. Since both water vapor and dry air are taken as ideal gases, the masses can be replaced by their the partial pressures,

$$ W = 0.621945 {p_{vap} \over {p - p_{vap}}} $$

or

$$ {1 \over {p - p_{vap}}} = {\displaystyle{1+{W \over 0.621945}} \over p} $$

where $p$ is the total pressure and $p_{vap}$ is the partial pressure of water vapor. Here, the total pressure is the atmospheric pressure at sea level, 101325 Pa. Analogously, the humidity of saturated air is the saturation humidity,

$$ W_{sat} = 0.621945 {p_{sat} \over {p-p_{sat}}} $$

Relative humidity is the material ratio of water vapor to the water vapor at saturation,

$$ \phi = {p_{vap} \over p_{sat}} $$

Note that relative humidity is not the ratio of humidity to saturation humidity. This is so because humidities are not fractions.

Specific Enthalpy & Specific Volume

Consider the adiabatic saturation of humid air with water. The amount of water required is the difference of humidity between the saturation and the humid air. The amount of dry gas is unchanged in the process. That is all about material balances. Taking water at the saturation temperature as reference for enthalpy, the enthalpies per mass of dry air, or specific enthalpy, of the inlet humid air, the inlet water and the outlet saturated air are given by

$$ h = c_{dry}\ (T - T_{sat}) + W\ (c_{vap}\ (T - T_{sat}) + \lambda_{sat}) $$

$$ h_{liq} = c_{liq}\ (W_{sat} - w)\ (T_{liq} - T_{sat}) $$

$$ h_{sat} = \lambda_{sat}\ W_{sat} $$

where $c_{dry}$ and $c_{vap}$ are the heat capacities of the dry gas and of water vapor and $\lambda_{sat}$ is the vaporization latent heat at the saturation temperature. All together, the energy balance gives

$$ {{W - W_{sat}} \over {T - T_{sat}}} = {c \over {-\lambda_{sat} + c_{vap}\ (T_{liq} - T_{sat})}} $$

where $c$ is the heat capacity of the inlet humid air,

$$ c = c_{dry} + c_{vap}\ (W_{sat} - W) $$

As the vaporization latent heat is usually much higher than the sensible heat, the variation of humidity in the gaseous phase is closely proportional to its variation in temperature for constant specific enthalpy, producing fairly straight lines in the psychrometric chart.

The volume of the gaseous mixture per mass of dry air, or specific volume, is given by

$$ v = {\displaystyle{nRT_{dry} \over p} \over {m_{dry}}} = R_{air}T_{dry} {\displaystyle{1+{W \over 0.621945}} \over p} $$

At room temperature at atmospheric pressure, humidity is closely proportional to dry bulb temperature for constant specific volume, producing fairly straight lines in the psychrometric chart.

Dry Bulb Temperature & Wet Bulb Temperature

If the gaseous phase in contact with water is not saturated with water vapor, the system is not at thermodynamic equilibrium. By removing sensible heat from its surroundings, some molecules overcome the vaporization heat and escape from the condensed to the gaseous phase. It happens spontaneously increasing the amount of water vapor in the gaseous phase and decreasing the temperature of the system.

The temperature in the surroundings of the evaporating molecules is the wet bulb temperature. This is so because of the construction of the simplest apparatus to indirectly read the air humidity, composed of two bulb thermometers, one in direct contact with the gaseous phase and one in contact with a thin layer of water in contact with the gaseous phase. The temperature of the wet bulb is affected by the evaporation of water from the thin layer to the gaseous phase. The temperature of the gaseous phase is the dry bulb temperature, read at the dry bulb.

The spontaneous heat and mass transfer phenomena are given by

$$ q = h\ (T - T_{wet}) $$

$$ N = k\ (p_{wet} - p_{vap}) $$

where $h$ and $k$ are the heat and mass transfer coefficients, $p_{wet}$ is the saturation pressure at the wet bulb temperature and $p_{vap}$ is water vapor pressure of the gaseous phase. In most cases, both pressures are much smaller than the total pressure, so the mass flux can be approximate to

$$ N = k^*\ (W_{wet} - W) $$

where $k^*$ is a mass transfer coefficient. The energy removed from water adjacent to the wet bulb is used to evaporate part of that water

$$ q = N\ \lambda_{wet} $$

where $\lambda_{wet}$ is the heat of vaporization at the wet bulb temperature. Combining mass and energy transfer phenomena, one has

$$ {{W - W_{wet}} \over {T - T_{wet}}} = -{h \over {\lambda_{wet}\ k^*}} $$

Therefore, the variation of humidity in the gaseous phase is closely proportional to its variation in temperature, producing fairly straight lines in the psychrometric chart for constant wet bulb temperature.

Empirical Equations

Equations used in psychrometrics toolbox come from the first chapter of the 2017 ASHRAE Handbook Fundamentals Systems - International Metric System, published by the American Society of Heating, Refrigerating and Air-Conditioning Engineers.

For ice in the range -100 °C to 0 °C, the water vapor pressure in equilibrium with pure ice is given by

$$ \ln p_{sat}^* = {C_1 \over T} + C_2 + C_3\ T + C_4\ T^2 + C_5\ T^3 + C_6\ T^4 + C_7 \ln T $$

and for water in the range 0 °C to 200 °C, the water vapor pressure in equilibrium with pure water is given by

$$ \ln p_{sat}^* = {C_8 \over T} + C_9 + C_{10}\ T + C_{11}\ T^2 + C_{12}\ T^3 + C_{13} \ln T $$

where $p^*_{sat}$ is given in Pa and dry bulb temperature $T$, in K.

As the saturation of air in ice and water is negligible, the water vapor pressure over pure ice or water is almost the same as the saturation pressure over ice or water when there is air in the gaseous phase.

The specific volume and the specific enthalpy (volume and enthalpy of the gaseous phase per unit of mass of dry air) are given by

$$ v = {{0.287042\ (t_{dry} + 273.15)\ (1 + 1.607858\ W)} \over p} $$

$$ h = 1.006\ t_{dry} + W\ (2501 + 1.86\ t_{dry}) $$

with $v$ given in m³/kg of dry air, $h$ in kJ/kg, dry bulb temperature $t_{dry}$ in °C, $W$ in kg/kg_dry, and total pressure $p$ in kPa.

Dew point is given by

$$ t_{dew} = C_{14} + C_{15}\ \alpha + C_{16}\ \alpha^2 + C_{17}\ \alpha^3 + C_{18}\ p_{vap}^{0.1984} $$

where $\alpha = ln\ p_{vap}$, with water vapor pressure $p_{vap}$ given in kPa and dew point temperature $t_{dew}$ in °C. The parameters for the equations used in psychrometrics toolbox are

$$ C= \begin{bmatrix} \begin{align*} -5&.6745359 \times 10^3 \\ 6&.3925247 \\ -9&.6778430 \times 10^{-3} \\ 6&.2215701 \times 10^{-7} \\ 2&.0747825 \times 10^{-9} \\ -9&.4840240 \times 10^{-13} \\ 4&.1635019 \\ -5&.8002206 \times 10^3 \\ 1&.3914993 \\ -4&.8640239 \times 10^{-2} \\ 4&.1764768 \times 10^{-5} \\ -1&.4452093 \times 10^{-8} \\ 6&.5459673 \\ 6&.54 \\ 14&.526 \\ 0&.7389 \\ 0&.09486 \\ 0&.4569 \\ \end{align*} \end{bmatrix} $$

The psychrometrics Toolbox for Scilab

psychrometrics provides the following functions:

• psychro
• humidity
• satPress
• enthalpy
• volume
• adiabSat

psychro

psychro computes the dry bulb temperature, the wet bulb temperature, the dew point temperature, the adiabatic saturation temperature, the humidit, the saturation humidity, the saturation humidity at the wet bulb temperature, the adiabatic saturation humidity, the specific enthalpy, the specific volume, the the relative humidity, the water vapor pressure, the water saturation pressure, the saturation pressure at the wet bulb temperature and the density, given any two input arguments, except for the combination of water vapor pressure and dew point temperature, which are not mutually independent. A schematic psychrometric chart can the presented as a graphical representation of the solution.

Syntax:

[Tdry,Twet,Tdew,Tadiab,W,Wsat,Wsatwet,Wadiab,h,v,phi,pw,psat,psatwet,rho]=..
psychro(Tdry,Twet,Tdew,W,h,v,phi,fig)

Example:

Compute the dry bulb temperature, the wet bulb temperature, the dew point temperature, the adiabatic saturation temperature, the humidity, the saturation humidity, the saturation humidity at wet bulb temperature, the adiabatic saturation humidity, the specific enthalpy, the specific volume, the relative humidity, the water vapor pressure, the saturation pressure, the saturation pressure at wet bulb temperature and the density given the specific enthalpy is 79.5 kJ/kg and the relative humidity is 0.29 and plot a graphical representation of the answer in a schematic psychrometric chart:

[Tdry,Twet,Tdew,Tadiab,W,Wsat,Wsatwet,Wadiab,h,v,phi,pw,psat,psatwet,rho]=..
psychro(h=79.5e3,phi=.29,fig=%t) // inputs and outputs in SI units

humidity

humidity computes the humidity of humid air given the water vapor pressure and the total pressure. By default, total pressure is assumed to be the atmospheric pressure at sea level (101325 Pa).

Syntax:

[W]=humidity(pw,p)

Example:

Compute the humidity of humid air at atmospheric pressure given water vapor pressure is 1 kPa at 1 atm total pressure:

pw=1e3; // water vapor pressure in Pa
W=humidity(pw) // saturation pressure in kg/kg of dry air

satPress

satPress computes the saturation pressure of humid air given the dry bulb temperature.

Syntax:

[psat]=satPress(Tdry)

Example:

Compute the saturation pressure given the dry bulb temperature is 25 °C:

Tdry=25+273.15; // dry bulb temperature in K
psat=satPress(Tdry) // saturation pressure in Pa

enthalpy

enthalpy computes the specific enthalpy of humid air given the dry bulb temperature and the humidity.

Syntax:

[h]=enthalpy(Tdry,W)

Example:

CCompute the specific enthalpy given the dry bulb temperature is 25 °C and the humidity is 7 g/kg of dry air:

Tdry=25+273.15; // dry bulb temperature in K
W=7e-3; // humidity in kg/kg of dry air
h=enthalpy(Tdry,W) // specific enthalpy in J/kg of dry air

volume

volume computes the specific volume of humid air given the dry bulb temperature the humidity and the total pressure. By default, total pressure is assumed to be the atmospheric pressure at sea level (101325 Pa).

Syntax:

[v]=volume(Tdry,W,p)

Example:

CCompute the specific volume given the dry bulb temperature is 25 °C and the humidity is 7 g/kg of dry air at 1 atm total pressure:

Tdry=25+273.15; // dry bulb temperature in K
W=7e-3; // humidity in kg/kg of dry air
v=volume(Tdry,W) // specific volume in cu. m/kg of dry air

adiabSat

adiabSat Compute the temperature and the humidity at saturation condition given the specific enthalpy. A schematic psychrometric chart can the presented as a graphical representation of the solution.

Syntax:

[Tadiab,Wadiab]=adiabSat(h,fig)

Example:

Compute the temperature and the humidity at saturation condition given the specific enthalpy 82.4 kJ/kg of dry air and plot a graphical representation of the answer in a schematic psychrometric chart:

h=82.4e3; // enthalpy in J/kg dry air
[T,W]=adiabSat(h,%t) // inputs and outputs in SI units

See Also

Internal-Fluid-Flow-for-Scilab.

Copyright © 2023 Alexandre Umpierre

email: aumpierre@gmail.com