Taller - reactor CSTR
Enviado por Bluemoon1018 • 6 de Noviembre de 2015 • Práctica o problema • 776 Palabras (4 Páginas) • 278 Visitas
Nombre: Silvia Alejandra Díaz Suarez. Cód.: 2102470. Fecha: 18/06/15
Profesor: Gustavo Ramírez.
Materia: Diseño de Reactores.
TALLER #3
- Teniendo en cuenta la siguiente información, para un reactor CSTR donde ocurre una
reacción irreversible con un solo reactivo : grafique el efecto de la temperatura de la corriente de los productos (entre y ) y las conversiones , sobre el volumen del reactor , el diámetro del reactor , el área de transferencia de calor , el calor transferido a la chaqueta , la temperatura de salida del fluido de enfriamiento , y el flujo del fluido de enfriamiento . Analice los resultados.[pic 2][pic 3][pic 4][pic 5][pic 6][pic 7][pic 8][pic 9][pic 10][pic 11][pic 12]
Factor preexponencial [pic 13] | [pic 14] | [pic 15] |
Energía de activación [pic 16] | [pic 17] | [pic 18] |
Peso molecular | 100 | |
Densidades y [pic 19][pic 20] | [pic 21] | [pic 22] |
Densidad del enfriador [pic 23] | [pic 24] | [pic 25] |
Capacidad calorífica y [pic 26][pic 27] | [pic 28] | [pic 29] |
Capacidad calorífica enfriador [pic 30] | [pic 31] | [pic 32] |
Calor de reacción [pic 33] | [pic 34] | [pic 35] |
Temperatura de entrada [pic 36] | [pic 37] | [pic 38] |
Caudal de entrada [pic 39] | [pic 40] | [pic 41] |
Composición alimento [pic 42] | [pic 43] | [pic 44] |
Temperatura entrada enfriador [pic 45] | [pic 46] | [pic 47] |
Coeficiente de transferencia de calor [pic 48] | [pic 49] | [pic 50] |
Relación [pic 51] | [pic 52] | |
[pic 53] | [pic 54] | [pic 55] |
Codigo en scilab
//TALLER #3: BALANCE REACTOR CSTR
clc
clear
k0=20.75*10^6
E=69.71*10^6
PM=100
d=801 //d0=d
dj=1000
cp=3137 //cp0=cp
cj=4183
l=-69.71*10^6
T0=294
q0=4.377*10^(-3)
ca0=8.01
Tcin=294
U=851
REL=2 // (L/D)
R=8314
Fa0=3.506
x=[0.5 0.8 0.9]
j=1
for x=x
i=1
for Tr=320:360
ca(j)=ca0-x*ca0
k(i)=k0*exp(-E/(R*Tr))
Vr(i,j)=(q0*ca0-q0*ca(j))/(k(i)*ca(j))
D(i,j)=(2*Vr(i,j)/3.1415)^(1/3)
Aj(i,j)=2*3.1415*D(i,j)^2
Q(i,j)=q0*ca0*PM*(cp*T0-cp*Tr)-l*k(i)*ca(j)*Vr(i,j)
Tj(i,j)=Tr-Q(i,j)/(U*Aj(i,j))
qj(i,j)=(-Q(i,j))/(cp*(Tcin-Tj(i,j)))
T(i)=Tr
i=i+1
end
j=j+1
end
plot(T,Vr)
xtitle("Volumen del reactor vs Temperatura del reactor (x=0.5 ; x=0.8 x=0.9)",..
"Tr","Vr")
legend("x=0.5","x=0.8","x=0.9")
plot(T,D)
xtitle("Diametro del reactor vs Temperatura del reactor (x=0.5 ; x=0.8 x=0.9)",..
"Tr","D")
legend("x=0.5","x=0.8","x=0.9")
plot(T,Aj)
...