The MATLAB program (psychro.m) (shown below) is used to plot the Simplified Psychrometric Chart shown above. Data is read from 2 datafiles containing water saturation temperature/pressure data, file t_pg over 1°C intervals for plotting the saturation and constant relative humidity curves, and file t_pg1 over 5°C intervals for plotting the wet-bulb temperature and enthalpy lines. The source of the data: NIST Chemistry WebBook - accessed Feb 2008
The following four equations (refer back to Part 1 and Part 2) are evaluated in the program:
specific humidity: (sometimes known as: humidity ratio or absolute humidity)::
Notice in the psychrometric chart that the oblique
enthalpy axis is given by
= T+5, for T varying between 0°C to 25°C. Substituting this value
of
in the enthalpy equation (above) leads to the intercept of the
enthalpy line with the enthalpy axis as follows:
h = T + 2.5 (T + 5)
or finally: T = (h - 12.5)/3.5
Thus the enthalpy lines can be plotted parallel to, but independant of the wet bulb temperature lines.
The complete MATLAB program follows. Note that it is
the convention in programming that all variable names begin with
lower case letters. Thus t represents temperature [°C], p represents
pressure [kPa], w (in place of
)
represents specific humidity [grams / kg dry air] (aka: absolute
humidity, humidity ratio), h represents enthalpy (kJ/kg dry air) and
the suffix g represents the saturated vapor state (following the
convention used in steam tables).
% Plotting the Simplified Psychrometric Chart
% Izzi Urieli 2/27/2008
tpg = dlmread('t_pg','\t'); % saturation temp/pressure
t = tpg(:,1); % temperature (C)
pg = tpg(:,2); % saturation vapor pressure (kPa)
patm = 101.325; % standard atmosphere (kPa)
rair = 0.287; % gas constant of air (kJ/kg.K)
wg = 622*pg./(patm-pg); % saturation specific humidity
plot(t,wg,'r-')
hold
grid
for phi = 0.1:0.1:0.4, % phi = relative humidity 10% - 40%
w = 622*phi*pg./(patm-phi*pg);
plot(t,w)
end
for phi = 0.6:0.2:0.8, % phi = 60%, 80%
w = 622*phi*pg./(patm-phi*pg);
plot(t,w)
end
% specific volume and enthalpy/wet-bulb-temp
tpg1 = dlmread('t_pg1','\t');
t1 = tpg1(:,1); % saturation temperature (C)
pg1 = tpg1(:,2); % saturation pressure (kPa)
wg1 = 622*pg1./(patm-pg1); % saturation specific humidity
% specific volume of dry air (cubic m/kg dry air) (green)
vol = rair.*(t1+273)./(patm-pg1) % specific vol at saturation
tv0 = patm*vol/rair-273; % air temperature at zero humidity
for i = 1:7,
plot([t1(i),tv0(i)],[wg1(i),0],'g-')
end
% wet bulb temperature (also enthalpy) lines (red)
h = t1 + 2.5*wg1 % enthalpy (kJ/kg-dry-air) (displayed)
t0 = h; % temperature at zero humidity for enthalpy h
for i = 1:6,
plot([t1(i),t0(i)],[wg1(i),0],'r-')
end
% enthalpy axis and enthalpy lines (black)
for h = 10:10:110, % enthalpy (kJ/kg-dry-air)
t0 = h; % temperature at zero humidity
t1 = (h - 12.5)/3.5; % temperature on the enthalpy axis
w1 = t1 + 5; % specific humidity on the enthalpy axis
plot([t0,t1],[0,w1],'k-')
end
plot([0,25],[5,30],'k-') % the oblique enthalpy axis
axis([0,50,0,30]) % limit the range of the chart
title('Simplified Psychrometric Chart')
xlabel('Dry Bulb Temperature (deg C)')
ylabel('Specific Humidity (gm vap/kg dry air)')
|
Notice in the program that both the specific volume (vol) and the enthalpy values (h) are displayed (no semicolon) thus they can be subsequently added to the plot together with the relevant saturation/wet bulb temperatures. relative humidity values.and enthalpy values on the enthalpy axis.
