Una Planta Requiere Que 15 KG
Una Planta Requiere Que 15 KG
Una Planta Requiere Que 15 KG
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/mL. s y una rapidez del
gas de 2.0 kg/m* . 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.
SOLUCION:
Los cálculos de 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 y
una grafica de las curvas de operación son calculadas con el siguiente
algoritmo programado en scilab: (archivo adjunto)
//CALCULO DE UNA TORRE ENFRAIDORA DE AGUA CON AIRE
clear
clc
L2=input('DIGITE EL FLUJO DE LIQUIDO A ENFRIAR EN Kg POR SEGUNDO');
TL2=input('DIGITE LA TEMPERATURA DE ENTRADA DEL LIQUIDO');
TL1=input('DIGITE LA TEMPERATURA A LA QUE DEBE SALIR LA CORRIENTE DE AGUA DEL
ENFRIADOR');
e=0.1;
Qs=input('DIGITE EL FLUX MINIMO DE GAS EN UNIDADES INTERNACIONALES');
P=input('DIGITE LA PRESION DE LA ENTRADA DE AIRE EN PASCALES');
TG1=input('DIGITE LA TEMPERATURA DE BULBO SECO DEL AIRE ENFRIADOR DISPONIBLE');
TH=input('DIGITE LA TEMPERATURA HUMEDA DEL AIRE ENTRANTE');
hla=150;
kya=input('DIGITE kya EN UNIDADES INTERNACIONALES');
B=input('DIGITE 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;
LAMBDAHUMEDO=0;
if TH<57
PSH1=exp(A1-(B1/(C1+TH)));
else PSH1=exp(A2-(B2/(C2+TH)));
end
YSH1=0.62069*(PSH1/(P-PSH1));
//CALCULO DE CALOR LATENTE DE VAPORIZACI?N
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;
LAMBDAHUMEDO=lambda(T2+1)-((lambda(T1+1)-lambda(T2+1))*(TH-(T2)));
W=(TG1-TH)/LAMBDAHUMEDO;
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('ESTA BENDITA TORRE 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=1:1:(N-1)
POTEN(n)=(1/(Hi(n)-H(n)));
end
s=1;
SUM1=0;
for s=1:1:((N-2)/2)
SUM1=POTEN(2*s)+SUM1;
end
SUM2=0;
t=1;
for t=1:1:((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=1:1:(N-1)
TLI(i)=TL(i);
end
plot2d(TLI,Hi);
H1 =
H2 =
37.
INTEGR =
Es necesario que sea de 5.50m 2 con tal de que se mantenga la velocidad de flujo de gas
mínima para que Kya=0.9
Ejemplo 7.12
En el enfriador del ejemplo 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?
Cpl=4.184;
f=5;
LAMBDA=0;
dHdT=0;
Tc=0;
Gsmin=0;
Gs=11.55
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)+2504.1*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)+2504.1*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)+LAMBDATij*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)+LAMBDATik*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('ESTA BENDITA TORRE 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)+LAMBDATGkl);
k=k+1;
end
end
POTEN=zeros(N);
for n=1:1:(N-1)
POTEN(n)=(1/(Hi(n)-H(n)));
end
s=1;
SUM1=0;
for s=1:1:((N-2)/2)
SUM1=POTEN(2*s)+SUM1;
end
SUM2=0;
t=1;
for t=1:1:((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=1:1:(N-1)
TLI(i)=TL(i);
end
plot2d(TLI,Hi);
Ejemplo 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 de la situación:
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
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)
101325.
Digite el flujo másico de aire seco que se va a soplar (kg/s)-->3.562
Digite la longitud de la cámara enfriadora (m) -->2
Digite la temperatura de salida de la cámara sencilla (ºC)-->42
Favor digite TG1 de la mezcla en °C-->65
Favor digite Y1 de la mezcla en Kg agua/Kg gas seco-->0.017
Digite cuantas cámaras es serie va a usar-->2
TSAcalculada =
32.612237ºC
YSA =
0.0307765
Y2 =
0.0447688 Kg/s
Operaciones unitarias 1