OPERACIONES UNITARIAS I
Jennyfer Cortés Peña Cód: 2033526 Grupo: O 2
Ej. 7-11:
Una planta requiere que 15 kg/s (1 984 Ib/min) de agua refrigerante fluyan a través
del equipo de condensación para la destilación; se eliminan así 270 W (55 270
Btu/min) de los condensadores. El agua sale de los condensadores a 45 ºC. Para
volver a utilizar el agua, se planea enfriarla mediante el contacto con aire en una
torre de enfriamiento de tiro inducido. Las condiciones del diseño son: aire
entrante a 30 ºC temperatura de bulbo seco, 24 ºC temperatura de bulbo húmedo;
el agua se va a enfriar a 5 ºC de la temperatura de bulbo húmedo del aire
entrante, o sea, a 29 ºC; se usará una relación de aire/vapor de agua de 1.5 veces
el mínimo. El agua de compensación va a entrar de una represa a 10 ºC, dureza
de 500 ppm (partes/millón) de sólidos disueltos. El agua circulante no debe tener
una dureza superior a 2 000 ppm. Respecto aI empaque que se va a utilizar, se
espera que Kya sea de 0.90kg/m3 *s (Y)(202 Ibm/h ft2 Y).
Para una rapidez del líquido como mínimo, de 2.7 kg/m 2. s y una rapidez del gas
de 2.0 kg/m2.s (1 991 y 1 474 Ib,/h ft2, respectivamente). Calcular las dimensiones
de la sección empacada y las necesidades de agua de compensación.
Solución:
Con los datos proporcionados y a partir del siguiente algoritmo se puede calcular:
el Gs mínimo, las condiciones de entalpía de entrada y salida de la corriente
gaseosa, la altura de la torre, la cantidad de agua de reposición.
// TORRE ENFRIADORA DE AGUA CON AIRE
clear
clc
L2=input('INGRESE EL VALOR DEL FLUJO DE LIQUIDO A ENFRIAR EN Kg /S');
TL2=input('INGRESE EL VALOR DE LA TEMPERATURA DE ENTRADA DEL LIQUIDO');
TL1=input(' INGRESE EL VALOR DE LA TEMPERATURA A LA QUE DEBE SALIR LA CORRIENTE DE AGUA
DEL ENFRIADOR');
//PARA UN ERROR FIJO DE;
e=0.1;
Qs=input('INGRESE EL FLUJO MINIMO DE GAS(Kg/m2*s )');
P=input(' INGRESE LA PRESION DE LA ENTRADA DE AIRE EN PASCALES');
TG1=input('INGRESE LA TEMPERATURA DE BULBO SECO DEL AIRE ENFRIADOR DISPONIBLE');
TH=input('INGRESE EL VALOR DE LA TEMPERATURA HUMEDA DEL AIRE ENTRANTE');
hla=150;
kya=input('INGRESE EL VALOR DE kya (S.I.)');
B=input('CUANTAS VECES QUIERE USAR EL MINIMO FLUJO DE GAS?');
A1=23.7093;
B1=4111.1;
C1=237.7;
A2=23.863;
B2=3809.4;
C2=226.7;
PSH1=0;
YSH1=0;
Y1=0;
W=0;
LAMBDATG1=0;
//LAMBDA HUMEDO=0
LW=0;
if TH<57
PSH1=exp(A1-(B1/(C1+TH)));
else PSH1=exp(A2-(B2/(C2+TH)));
end
YSH1=0.62069*(PSH1/(P-PSH1));
//CALOR LATENTE DE VAPORIZACION
lambda=[2501.4 2499.0 2496.7 2494.3 2491.9 2489.6 2487.1 2484.8 2482.5 2480.1 2477.7 2475.4 2473.0
2470.7 2468.3 2465.9 2463.6 2461.2 2458.8 2456.5 2454.1 2451.8 2449.4 2447.0 2444.7 2442.3 2439.9
2437.6 2435.2 2432.8 2430.5 2428.1 2425.7 2423.4 2421.0 2418.6 2416.2 2413.9 2411.5 2409.1 2406.7
2404.3 2401.9 2399.5 2397.2 2394.8 2392.4 2390.0 2387.6 2385.2 2382.7 2380.3 2377.9 2375.5 2373.1
2370.7 2368.2 2365.8 2363.4 2360.9 2358.5 2356.0 2353.6 2351.1 2348.7 2346.2 2343.7 2341.3 2338.8
2336.3 2333.8 2331.4 2338.9 2326.4 2323.9 2321.4 2318.9 2316.3 2312.8 2311.3 2308.8 2306.2 2303.7
2301.1 2298.6 2296.0 2293.5 2290.9 2288.3 2285.8 2283.2 2280.6 2278.0 2275.4 2272.8 2270.2 2267.6
2264.9 2262.3 2259.7 2257.0];
T1=int(TH);
T2=T1+1;
LW=lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(TH-(T2)));
W=(TG1-TH)/LW;
Y1=((YSH1-(W*1.0035))/(W*1.8723+1));
T1=int(TG1);
T2=T1+1;
LAMBDATG1=lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(TG1-(T2)));
H1=(1.0035+1.8723*Y1)*TG1+LAMBDATG1*Y1
//CALCULO DEL FLUJO MINIMO DE GAS PARA UNA TORRE EMPACADA PARA EL ENFRIAMIENTO DE
AGUA
T=TL1+1;
A1=23.7093;
B1=4111.1;
C1=237.7;
A2=23.1863;
B2=3809.4;
C2=226.7;
Cpl=4.184;
f=5;
LAMBDA=0;
dHdT=0;
Tc=0;
Gsmin=0;
while abs(f)>=e
if T<57
Ps=exp(A1-B1/(C1+T));
else
Ps=exp(A2-B2/(C2+T));
end
Ysa=0.62069*(Ps/(P-Ps));
//CALCULO DE CALOR LATENTE DE VAPORIZACIÓN
T1=int(T);
T2=T1+1;
//LAMBDA=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(T-(T2))));
H=(1.0035+1.8723*Ysa)*T+2501.4*Ysa;
if T<57
dHdT=(B1*P*Ysa)/(((C1+T)^2)*(P-Ps))*(1.8723*T+2501.4)+1.8723*Ysa+1.0035;
else
dHdT=(B2*P*Ysa)/(((C2+T)^2)*(P-Ps))*(1.8723*T+2501.4)+1.8723*Ysa+1.0035;
end
Tc=TL1+(H-H1)/(dHdT);
f=T-Tc;
T=Tc;
end
if Tc>=TL2
T=TL2;
else T=Tc;
end
Gsmin=(L2*Cpl*(T-TL1))/(H-H1)
Gs=Gsmin*B
H2=((L2*Cpl)/Gs)*(TL2-TL1)+H1
N=11;
deltaH=5
while abs(deltaH)>=2.5
deltaH=((H2-H1)/(N-1));
N=N+2;
end
deltaTL=(Gs/(L2*Cpl))*deltaH
N=N-2
TL=zeros(N);
TG=zeros(N);
H=zeros(N);
Y=zeros(N);
deltaTG=zeros(N);
k=1;
TL(k)=TL1;
TG(k)=TG1;
H(k)=H1;
Y(k)=Y1;
Ti=zeros(N);
Hi=zeros(N);
while k<N
if kya==0
Ti(k)=TL(k);
T1=int(Ti(k));
T2=T1+1;
//LAMBDATik=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(Ti(k)-(T2))));
if Ti(k)<57
psi=exp(A1-(B1/(C1+Ti(k))));
else psi=exp(A2-(B2/(C2+Ti(k))));
end
Ysi=0.6209*(psi/(P-psi));
Hsi=(1.0035+1.8723*Ysi)*Ti(k)+2501.4*Ysi;
Hi(k)=Hsi;
else
j=1;
Ti(j)=TL(k)-0.5*deltaTL;
TC=zeros(N);
Hic=zeros(N);
T1=int(Ti(j));
T2=T1+1;
//LAMBDATij=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(Ti(j)-(T2))));
if Ti(j)<57
psj=exp(A1-B1/(C1+Ti(j)));
else psj=exp(A2-B2/(C2+Ti(j)));
end
Ysj=0.62069*(psj/(P-psj));
Hsj=(1.0035+1.8723*Ysj)*Ti(j)+2501.4*Ysj;
Hic(j)=Hsj;
TC(j)=TL(k)+(Hic(j)-H(k))/(-1*(hla/kya));
AUX1=8;
while AUX1>e
j=j+1;
Ti(j)=TC(j-1);
T1=int(Ti(j));
T2=T1+1;
//LAMBDATij=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(Ti(j)-(T2))));
if Ti(j)<57
psj=exp(A1-B1/(C1+Ti(j)));
else psj=exp(A2-B2/(C2+Ti(j)));
end
Ysj=0.62069*(psj/(P-psj));
Hsj=(1.0035+1.8723*Ysj)*Ti(j)+2501.4*Ysj;
Hic(j)=Hsj;
TC(j)=TL(k)+(Hic(j)-H(k))/(-1*(hla/kya));
AUX1=abs(TC(j)-Ti(j));
end
Ti(k)=TC(j);
T1=int(Ti(k));
T2=T1+1;
//LAMBDATik=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(Ti(k)-(T2))));
if Ti(k)<57
psk=exp(A1-B1/(C1+Ti(k)));
else psk=exp(A2-B2/(C2+Ti(k)));
end
Ysk=0.62069*(psk/(P-psk));
Hsk=(1.0035+1.8723*Ysk)*Ti(k)+2501.4*Ysk;
Hi(k)=Hsk;
if TG(k)<57
Ps=exp(A1-B1/(C1+TG(k)));
else Ps=exp(A2-B2/(C2+TG(k)));
end
PA=(Y(k)*P)/(1+Y(k));
if PA >Ps
disp('LATORRE ES INOPERABLE: CONSIDERE UN FLUJO MAYOR DE GAS')
break
else
deltaTG(k)=deltaH*(Ti(k)-TG(k))/(Hi(k)-H(k));
end
H(k+1)=H(k)+deltaH;
TL(k+1)=TL(k)+deltaTL;
TG(k+1)=TG(k)+deltaTG(k);
T1=int(TG(k+1));
T2=T1+1;
//LAMBDATGkl=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(TG(k+1)-(T2))));
Y(k+1)=(H(k+1)-1.0035*TG(k+1))/(1.8723*TG(k+1)+2501.4);
k=k+1;
end
end
POTEN=zeros(N);
for n=[Link](N-1)
POTEN(n)=(1/(Hi(n)-H(n)));
end
s=1;
SUM1=0;
for s=[Link]((N-2)/2)
SUM1=POTEN(2*s)+SUM1;
end
SUM2=0;
t=1;
for t=[Link]((N-4)/2)
SUM2=SUM2+POTEN(2*t+1);
end
INTEGR=(deltaH*(POTEN(1)+POTEN(N-1)+4*SUM1+2*SUM2))/3
Z=(INTEGR*Qs)/kya
AG=Gs*(Y(N)-Y(1))
plot2d(TL,H);
for i=[Link](N-1)
TLI(i)=TL(i);
end
plot2d(TLI,Hi);
Resultados obtenidos
INGRESE EL VALOR DEL FLUJO DE LIQUIDO A ENFRIAR EN Kg /S-->15
INGRESE EL VALOR DE LA TEMPERATURA DE ENTRADA DEL LIQUIDO-->45
INGRESE EL VALOR DE LA TEMPERATURA A LA QUE DEBE SALIR LA CORRIENTE DE AGUA
DEL ENFRIADOR-->29
INGRESE EL FLUJO MINIMO DE GAS(Kg/m2*s )-->2
INGRESE LA PRESION DE LA ENTRADA DE AIRE EN PASCALES-->101325
INGRESE LA TEMPERATURA DE BULBO SECO DEL AIRE ENFRIADOR DISPONIBLE-->30
INGRESE EL VALOR DE LA TEMPERATURA HUMEDA DEL AIRE ENTRANTE-->24
INGRESE EL VALOR DE kya (S.I.)-->0.9
CUANTAS VECES QUIERE USAR EL MINIMO FLUJO DE GAS?-->1.5
Entalpía de entrada de la corriente gaseosa, H1 = 70.583359 KJ/Kg de aire seco
Gsmin = 7.32 Kg aire seco/s
Gs = 10.98 Kg aire seco/s
Entalpía de salida de la corriente gaseosa H2 =162.03 KJ/Kg de aire seco
deltaTL = 0.4210526
N = 39
INTEGR = 3.3023159 valor de la integral
Altura de la torre, Z= 7.338 m
AG = 0.3439605
Para el flujo líquido de por lo menos 2.7kg/m2.s
Área transversal= 15 Kg/s.m2.s = 5.56 m2
2.7 Kg
Para un flujo de gas de por lo menos 2 kg/m2.s
Área transversal= 10.98 Kg aire seco/s.m2.s = 5.49 m2
2 Kg
Es necesario que sea de 5.49 m 2 con tal de que se mantenga la velocidad de flujo de gas
mínima para que Kya=0.9
Fig1. Gráfico temperatura del líquido [Link] de aire-vapor de agua
Ej. 7-12:
En el enfriador del Ej. 7.11, ¿a que temperatura se enfriara el agua si, después de
que la torre se construyó a los valores de diseño L ’ y G’, el aire entrante entrase a
una temperatura de bulbo seco TG= 32 ºC, temperatura de bulbo húmedo Tw = 28
ºC?
Solución
Para calcular la temperatura de salida del líquido se realiza una prueba y error en
la siguiente subrutina, programada en scilab. Puesto que el flujo tanto de gas
como de líquido permanecen constantes, el valor de la integral sera el mismo que
en el caso anterior y por lo tanto el valor de Z; es decir el valor de TL 1 será
finalmente el correspondiente a una Z de 7.37m, igual a la diseñada en el anterior
ejercicio.
// TORRE ENFRIADORA DE AGUA CON AIRE
clear
clc
L2=input('INGRESE EL VALOR DEL FLUJO DE LIQUIDO A ENFRIAR EN Kg /S');
TL2=input('INGRESE EL VALOR DE LA TEMPERATURA DE ENTRADA DEL LIQUIDO');
TL1=input(' INGRESE EL VALOR DE LA TEMPERATURA A LA QUE DEBE SALIR LA CORRIENTE DE AGUA
DEL ENFRIADOR');
//PARA UN ERROR FIJO DE;
e=0.1;
Qs=input('INGRESE EL FLUJO MINIMO DE GAS(Kg/m2*s )');
P=input(' INGRESE LA PRESION DE LA ENTRADA DE AIRE EN PASCALES');
TG1=input('INGRESE LA TEMPERATURA DE BULBO SECO DEL AIRE ENFRIADOR DISPONIBLE');
TH=input('INGRESE EL VALOR DE LA TEMPERATURA HUMEDA DEL AIRE ENTRANTE');
hla=150;
kya=input('INGRESE EL VALOR DE kya (S.I.)');
B=input('CUANTAS VECES QUIERE USAR EL MINIMO FLUJO DE GAS?');
A1=23.7093;
B1=4111.1;
C1=237.7;
A2=23.863;
B2=3809.4;
C2=226.7;
PSH1=0;
YSH1=0;
Y1=0;
W=0;
LAMBDATG1=0;
//LAMBDA HUMEDO=0
LW=0;
if TH<57
PSH1=exp(A1-(B1/(C1+TH)));
else PSH1=exp(A2-(B2/(C2+TH)));
end
YSH1=0.62069*(PSH1/(P-PSH1));
//CALOR LATENTE DE VAPORIZACION
lambda=[2501.4 2499.0 2496.7 2494.3 2491.9 2489.6 2487.1 2484.8 2482.5 2480.1 2477.7 2475.4 2473.0
2470.7 2468.3 2465.9 2463.6 2461.2 2458.8 2456.5 2454.1 2451.8 2449.4 2447.0 2444.7 2442.3 2439.9
2437.6 2435.2 2432.8 2430.5 2428.1 2425.7 2423.4 2421.0 2418.6 2416.2 2413.9 2411.5 2409.1 2406.7
2404.3 2401.9 2399.5 2397.2 2394.8 2392.4 2390.0 2387.6 2385.2 2382.7 2380.3 2377.9 2375.5 2373.1
2370.7 2368.2 2365.8 2363.4 2360.9 2358.5 2356.0 2353.6 2351.1 2348.7 2346.2 2343.7 2341.3 2338.8
2336.3 2333.8 2331.4 2338.9 2326.4 2323.9 2321.4 2318.9 2316.3 2312.8 2311.3 2308.8 2306.2 2303.7
2301.1 2298.6 2296.0 2293.5 2290.9 2288.3 2285.8 2283.2 2280.6 2278.0 2275.4 2272.8 2270.2 2267.6
2264.9 2262.3 2259.7 2257.0];
T1=int(TH);
T2=T1+1;
LW=lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(TH-(T2)));
W=(TG1-TH)/LW;
Y1=((YSH1-(W*1.0035))/(W*1.8723+1));
T1=int(TG1);
T2=T1+1;
LAMBDATG1=lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(TG1-(T2)));
H1=(1.0035+1.8723*Y1)*TG1+LAMBDATG1*Y1
//CALCULO DEL FLUJO MINIMO DE GAS PARA UNA TORRE EMPACADA PARA EL ENFRIAMIENTO DE
AGUA
T=TL1+1;
A1=23.7093;
B1=4111.1;
C1=237.7;
A2=23.1863;
B2=3809.4;
C2=226.7;
Cpl=4.184;
f=5;
LAMBDA=0;
dHdT=0;
Tc=0;
Gsmin=0;
while abs(f)>=e
if T<57
Ps=exp(A1-B1/(C1+T));
else
Ps=exp(A2-B2/(C2+T));
end
Ysa=0.62069*(Ps/(P-Ps));
//CALCULO DE CALOR LATENTE DE VAPORIZACIÓN
T1=int(T);
T2=T1+1;
//LAMBDA=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(T-(T2))));
H=(1.0035+1.8723*Ysa)*T+2501.4*Ysa;
if T<57
dHdT=(B1*P*Ysa)/(((C1+T)^2)*(P-Ps))*(1.8723*T+2501.4)+1.8723*Ysa+1.0035;
else
dHdT=(B2*P*Ysa)/(((C2+T)^2)*(P-Ps))*(1.8723*T+2501.4)+1.8723*Ysa+1.0035;
end
Tc=TL1+(H-H1)/(dHdT);
f=T-Tc;
T=Tc;
end
if Tc>=TL2
T=TL2;
else T=Tc;
end
Gsmin=(L2*Cpl*(T-TL1))/(H-H1)
Gs=Gsmin*B
H2=((L2*Cpl)/Gs)*(TL2-TL1)+H1
N=11;
deltaH=5
while abs(deltaH)>=2.5
deltaH=((H2-H1)/(N-1));
N=N+2;
end
deltaTL=(Gs/(L2*Cpl))*deltaH
N=N-2
TL=zeros(N);
TG=zeros(N);
H=zeros(N);
Y=zeros(N);
deltaTG=zeros(N);
k=1;
TL(k)=TL1;
TG(k)=TG1;
H(k)=H1;
Y(k)=Y1;
Ti=zeros(N);
Hi=zeros(N);
while k<N
if kya==0
Ti(k)=TL(k);
T1=int(Ti(k));
T2=T1+1;
//LAMBDATik=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(Ti(k)-(T2))));
if Ti(k)<57
psi=exp(A1-(B1/(C1+Ti(k))));
else psi=exp(A2-(B2/(C2+Ti(k))));
end
Ysi=0.6209*(psi/(P-psi));
Hsi=(1.0035+1.8723*Ysi)*Ti(k)+2501.4*Ysi;
Hi(k)=Hsi;
else
j=1;
Ti(j)=TL(k)-0.5*deltaTL;
TC=zeros(N);
Hic=zeros(N);
T1=int(Ti(j));
T2=T1+1;
//LAMBDATij=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(Ti(j)-(T2))));
if Ti(j)<57
psj=exp(A1-B1/(C1+Ti(j)));
else psj=exp(A2-B2/(C2+Ti(j)));
end
Ysj=0.62069*(psj/(P-psj));
Hsj=(1.0035+1.8723*Ysj)*Ti(j)+2501.4*Ysj;
Hic(j)=Hsj;
TC(j)=TL(k)+(Hic(j)-H(k))/(-1*(hla/kya));
AUX1=8;
while AUX1>e
j=j+1;
Ti(j)=TC(j-1);
T1=int(Ti(j));
T2=T1+1;
//LAMBDATij=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(Ti(j)-(T2))));
if Ti(j)<57
psj=exp(A1-B1/(C1+Ti(j)));
else psj=exp(A2-B2/(C2+Ti(j)));
end
Ysj=0.62069*(psj/(P-psj));
Hsj=(1.0035+1.8723*Ysj)*Ti(j)+2501.4*Ysj;
Hic(j)=Hsj;
TC(j)=TL(k)+(Hic(j)-H(k))/(-1*(hla/kya));
AUX1=abs(TC(j)-Ti(j));
end
Ti(k)=TC(j);
T1=int(Ti(k));
T2=T1+1;
//LAMBDATik=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(Ti(k)-(T2))));
if Ti(k)<57
psk=exp(A1-B1/(C1+Ti(k)));
else psk=exp(A2-B2/(C2+Ti(k)));
end
Ysk=0.62069*(psk/(P-psk));
Hsk=(1.0035+1.8723*Ysk)*Ti(k)+2501.4*Ysk;
Hi(k)=Hsk;
if TG(k)<57
Ps=exp(A1-B1/(C1+TG(k)));
else Ps=exp(A2-B2/(C2+TG(k)));
end
PA=(Y(k)*P)/(1+Y(k));
if PA >Ps
disp('LATORRE ES INOPERABLE: CONSIDERE UN FLUJO MAYOR DE GAS')
break
else
deltaTG(k)=deltaH*(Ti(k)-TG(k))/(Hi(k)-H(k));
end
H(k+1)=H(k)+deltaH;
TL(k+1)=TL(k)+deltaTL;
TG(k+1)=TG(k)+deltaTG(k);
T1=int(TG(k+1));
T2=T1+1;
//LAMBDATGkl=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(TG(k+1)-(T2))));
Y(k+1)=(H(k+1)-1.0035*TG(k+1))/(1.8723*TG(k+1)+2501.4);
k=k+1;
end
end
POTEN=zeros(N);
for n=[Link](N-1)
POTEN(n)=(1/(Hi(n)-H(n)));
end
s=1;
SUM1=0;
for s=[Link]((N-2)/2)
SUM1=POTEN(2*s)+SUM1;
end
SUM2=0;
t=1;
for t=[Link]((N-4)/2)
SUM2=SUM2+POTEN(2*t+1);
end
INTEGR=(deltaH*(POTEN(1)+POTEN(N-1)+4*SUM1+2*SUM2))/3
Z=(INTEGR*Qs)/kya
AG=Gs*(Y(N)-Y(1))
plot2d(TL,H);
for i=[Link](N-1)
TLI(i)=TL(i);
end
plot2d(TLI,Hi);
Para llegar a un TL1 que cumpliera con las condiciones descritas anteriormente se
realizo la prueba suponiendo esta temperatura, hasta llegar a un valor
aproximado de Z= 7.34 m
Los resultados obtenidos son:
DIGITE EL FLUJO DE LÍQUIDO A ENFRIAR EN Kg POR SEGUNDO-->15
DIGITE LA TEMPERATURA SUPUESTA A LA QUE DEBE SALIR LA CORRIENTE DE AGUA DEL
ENFRIADOR-->31.58
DIGITE LA TEMPERATURA DEL AIRE ENFRIADOR DISPONIBLE-->32
DIGITE LA TEMPERATURA HUMEDA DEL AIRE ENTRANTE-->28
DIGITE CUANTAS VECES QUIERE USAR EL MINIMO FLUJO DE GAS-->1.5
Temperatura de entrada del liquido TL2 = 47.57 en ºC
Entalpía de entrada del gas, H1 = 87.57en KJ/Kg.
Gs = 11.55 en kg aire seco/s valor de diseño
Entalpía del gas a la salida de la torre, H2 =174.47 KJ/Kg.
El valor de Z sigue siendo el mismo:
Z = 7.34 m valor de diseño.
Estos resultados se obtuvieron con una temperatura de salida de la corriente de agua del
enfriador de 31.58 ºC=TL1
Ej. 7-13:
Se va a utilizar una cámara de aspersión horizontal con recirculación de agua para
la humidificación adiabática y el enfriamiento de aire. La parte activa de
la cámara tiene 2 m de longitud y una sección transversal de 2 m 2. El flujo del aire
es de 3.5 m3/s a la temperatura de bulbo seco de 65.0 ºC, Y’ = 0.0170 kg agua/kg
aire seco; el aire se enfría y humidifica a la temperatura de bulbo húmedo 42.0 ºC.
Si se va a colocar una cámara de aspersión igual, manejada en la misma forma,
en serie con la cámara existente. ¿Que condiciones de salida se preverían para el
aire?
Solución
Datos:
L=2m
A=2m2
Q=3.5m/s
Y1=0.0170kg agua/Kg. Aire seco Gs =3.562kg aire seco/s
T1=65ºC
T2=42ºC
Las condiciones de salida de la cámara sencilla y de varias cámaras colocadas en
serie según se desee, son calculadas con el siguiente rutina programada en
scilab.
//CALCULO DE LAS CONDICIONES DE SALIDA DE UNA CAMARA ADIABATICA PARA HUMIDIFICACIÓN Y
ENFRIAMIENTO DE AIRE CON AGUA CON UNA O VARIAS CAMARAS CONECTADAS EN SERIE
P=input('DIGITE LA PRESION TOTAL DEL SISTEMA EN PASCALES')
ERROR=0.3;
Gs=input('digite el flujo masico de aire seco que se va a soplar');
Z=input('digite la longitud de la camara enfriadora');
TG2=input('digite la temperatura de salida de la camara sencilla');
TG1=input('favor digite TG1 de la mezcla en °C');
Y1=input('favor digite Y1 de la mezcla en Kg agua/Kg gas seco');
B=input(' Numero camaras en serie ');
ch=1.0035+1.8723*Y1;
TSA=TG1;
e=10;
while e>=ERROR
T=TSA;
lambda=[2501.4 2499.0 2496.7 2494.3 2491.9 2489.6 2487.1 2484.8 2482.5 2480.1 2477.7 2475.4 2473.0
2470.7 2468.3 2465.9 2463.6 2461.2 2458.8 2456.5 2454.1 2451.8 2449.4 2447.0 2444.7 2442.3 2439.9
2437.6 2435.2 2432.8 2430.5 2428.1 2425.7 2423.4 2421.0 2418.6 2416.2 2413.9 2411.5 2409.1 2406.7
2404.3 2401.9 2399.5 2397.2 2394.8 2392.4 2390.0 2387.6 2385.2 2382.7 2380.3 2377.9 2375.5 2373.1
2370.7 2368.2 2365.8 2363.4 2360.9 2358.5 2356.0 2353.6 2351.1 2348.7 2346.2 2343.7 2341.3 2338.8
2336.3 2333.8 2331.4 2338.9 2326.4 2323.9 2321.4 2318.9 2316.3 2312.8 2311.3 2308.8 2306.2 2303.7
2301.1 2298.6 2296.0 2293.5 2290.9 2288.3 2285.8 2283.2 2280.6 2278.0 2275.4 2272.8 2270.2 2267.6
2264.9 2262.3 2259.7 2257.0];
T1=int(T);
T2=T1+1;
LAMBDA=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(T-T2)));
YSA=(Y1+(ch*(TG1-TSA)/(LAMBDA)));
Ps=(P*((YSA*1.65806)/(1+(YSA*1.65806))));
A1=23.7093;
B1=4111;
C1=237.7;
A2=23.1863;
B2=3809.4;
C2=226.7;
T=((B1/(A1-log(Ps)))-C1);
if T<57
TSAcalculada=T;
else
TSAcalculada=(B2/(A2-log(Ps)))-C2;
end
e=abs(TSA-TSAcalculada);
TSA=TSAcalculada;
end
TSAcalculada
YSA
T1=int(TSAcalculada);
T2=T1+1;
LAMBDAs=(lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(TSAcalculada-T2)));
Y2=((YSA-(1.0035/LAMBDAs)*(TG2-TSAcalculada))/(1+(1.8723/LAMBDAs)*(TG2-TSAcalculada)))
INVHTG=((log((YSA-Y1)/(YSA-Y2))/Z))
Yfinal=YSA-((YSA-Y1)/(exp(INVHTG*Z*B)))
chn=1.0035+1.8723*Yfinal;
Tfinal=TSA+(LAMBDAs*(YSA-Yfinal))/(chn)
AGUA=Gs*(Yfinal-Y1)
-->DIGITE LA PRESION TOTAL DEL SISTEMA EN PASCALES-->101325
P =
101325.
Digite el flujo masico de aire seco que se va a soplar-->3.562
Digite la longitud de la camara enfriadora-->2
Digite la temperatura de salida de la camara sencilla-->42
Digite TG1 de la mezcla en °C-->65
Digite Y1 de la mezcla en Kg agua/Kg gas seco-->0.017
Numero camaras en serie -->2
Resultados
TSAcalculada = 32.61 ºC
YSA = 0.0308
Y2 = 0.0267 salida de la primera camara
INVHTG = 0.6084950
Las condiciones a la salida de la segunda camara son:
Yfinal = 0.0296
Tfinal = 35.38 ºC
