Ejemplo 8-3 isomerizacin en fase lquida de n-butano

Abraham Castillo Espinosa

Pgina 1

function dxdv=isomerizacion(v,x); Ca0=9.3; Fa0=0.9*163; T=330+43.3*(x(1,:)); Kc=3.03*exp(-830.3*((T-333)/(T*333))); k=31.1*exp(7906*(T-360)/(T*360)); Xe=Kc/(1+Kc); ra=-k*Ca0*(1-(1+1/Kc)*(x(1,:))); rate=-ra y=Fa0/rate; dxdv(1,:)=-ra/Fa0;

clear clc %Concidiones iniciales v=[0:0.1:3.6]; x0=[0]; ra0=[40]; [v,x]=ode45(@isomerizacion,v,x0); plot(v,x(:,1)) xlabel('V') ylabel('X') plot(v,T) xlabel('V') ylabel('T')

plot(v,rate) xlabel('V') ylabel('ra')

Abraham Castillo Espinosa

Pgina 2

Ejemplo 8-5 produccin de anhdrido actico.

function dxdv=acetico(v,x); Fa0=0.0376; Cpa=163; delCp=-9; Ca0=18.8; T0=1035; deltaH=80770+delCp*(x(2,:)-298); ra=-Ca0*3.58*exp(34222*(1/T0-1/(x(2,:))))*(1-x(1,:))*(T0/(x(2,:)))/(1x(1,:)); %conversion(X) dxdv(1,:)=-ra/Fa0; %Temperatura dxdv(2,:)=-ra*(-deltaH)/(Fa0*(Cpa+x(1,:)*delCp));

clear clc v=[0:0.01:4]; x0=[0 1035]; [v,x]=ode45(@acetico,v,x0); plot(v,x(:,1)) xlabel('V') ylabel('X') plot(v,x(:,2)) xlabel('V') ylabel('T')

Abraham Castillo Espinosa

Pgina 3

function dxdv=acetico2(v,x); Fa0=0.0376; Cpa=163; delCp=-9; Ca0=18.8; T0=1035; deltaH=80770+delCp*(x(2,:)-298); ra=-Ca0*3.58*exp(34222*(1/T0-1/(x(2,:))))*(1x(1,:))*(T0/(x(2,:)))/(1+x(1,:)); Ua=16500; Ta=1150;

Abraham Castillo Espinosa

Pgina 4

%conversion(X) dxdv(1,:)=-ra/Fa0; %Temperatura dxdv(2,:)=(Ua*(Ta-x(2,:))+ra*(-deltaH))/(Fa0*(Cpa+x(1,:)*delCp));

clear clc v=[0:0.00001:0.001]; x0=[0 1035]; [v,x]=ode45(@acetico2,v,x0); plot(v,x(:,1)) xlabel('V') ylabel('X') plot(v,x(:,2)) xlabel('V') ylabel('T')

Abraham Castillo Espinosa

Pgina 5

function dxdv=acetico3(v,x); Fa0=0.0376; Cpa=163; delCp=-9; Ca0=18.8; T0=1035; deltaH=80770+delCp*(x(2,:)-298); ra=-Ca0*3.58*exp(34222*(1/T0-1/(x(2,:))))*(1x(1,:))*(T0/(x(2,:)))/(1+x(1,:)); Ua=16500; mc=0.111; Cpc=34.5; %conversion(X) dxdv(1,:)=-ra/Fa0; %Temperatura dxdv(2,:)=(Ua*(x(3,:)-x(2,:))+ra*(-deltaH))/(Fa0*(Cpa+x(1,:)*delCp)); %Ta dxdv(3,:)=Ua*(x(2,:)-x(3,:))/mc/Cpc; clear clc v=[0:0.00001:0.001]; x0=[0 1035 1150]; [v,x]=ode45(@acetico3,v,x0); plot(v,x(:,1)) xlabel('V') ylabel('X') plot(v,x(:,2),v,x(:,3))

Abraham Castillo Espinosa

Pgina 6

xlabel('V') ylabel('T')

Abraham Castillo Espinosa

Pgina 7

Ejemplo 8-10 Reacciones en paralelo en un PFR con efectos de calor

function dxdv=paralelo(v,x); k1a=10*exp(4000*(1/300-1/x(4,:))); k2a=0.09*exp(9000*(1/300-1/x(4,:))); Cto=0.1; Ft=x(1,:)+x(2,:)+x(3,:); T0=423; Ca=Cto*(x(1,:)/Ft)*(T0/x(4,:)); Cb=Cto*(x(2,:)/Ft)*(T0/x(4,:)); Cc=Cto*(x(3,:)/Ft)*(T0/x(4,:)); r1a=-k1a*Ca; r2a=-k2a*Ca^2; %Fa dxdv(1,:)=r1a+r2a; %Fb dxdv(2,:)=-r1a; %Fc dxdv(3,:)=-r2a/2; %T dxdv(4,:)=(4000*(373-x(4,:))+(-r1a)*20000+(r2a)*60000)/(90*x(1,:)+90*x(2,:)+180*x(3,:));

v=[0:0.01:1]; x0=[100 0 0 423]; [v,x]=ode45(@paralelo,v,x0); plot(v,x(:,4)) xlabel('V') ylabel('T') plot(v,x(:,1),v,x(:,2),v,x(:,3)) xlabel('V') ylabel('Fi')

Abraham Castillo Espinosa

Pgina 8

Ejemplo 8-11 Reacciones mltiples en un CSTR

function dTdt=multiples(t,T); Cp=200; Ca0=0.3; T0=283; tau=0.01; DH1=-55000; DH2=-71500; v0=1000; E2=27000; E1=9900;

Abraham Castillo Espinosa

Pgina 9

UA=40000; Ta=330; k2=4.58*exp((E2/1.987)*(1/500-1/T(1,:))); k1=3.3*exp((E1/1.987)*(1/300-1/T(1,:))); Ca=Ca0/(1+tau*k1); kappa=UA/(v0*Ca0)/Cp; G=-tau*k1/(1+k1*tau)*DH1-k1*tau*k2*tau*DH2/((1+tau*k1)*(1+tau*k2)); Tc=(T0+kappa*Ta)/(1+kappa); Cb=tau*k1*Ca/(1+k2*tau); R=Cp*(1+kappa)*(T(1,:)-Tc); Cc=Ca0-Ca-Cb; F=G-R; %T dTdt(1,:)=2;

t=[0:100]; T0=[283]; [t,T]=ode45(@multiples,t,T0); plot(G,T,R,T) xlabel('T') ylabel('G o R')

Problema 8.8 La reaccin elemental reversible en fase gaseosa A= B+C Se efectuara adiabticamente en un PFR empacado con catalizador. Entra A puro al reactor con un flujo volumtrico de 20 dm 3/s y presin de 10 atm y temperatura de 450 K a) Grafique la conversin y la temperatura a lo largo del PFR hasta que se logre una conversin del 80% (de ser posible).(El peso mximo del catalizador que puede empacarse en el PFR es de 50 kg). Asuma que P =0.0 Abraham Castillo Espinosa Pgina 10

A+B = C Balance de masa dX/dW= - rA/FA0 v0=20 dm3/s P0=10 atm Estequeometria: CA= CA0 (-X/1+X) T/T0 con =1 CA=CA0 (1-X/1+X) T/T0 Ley de velocidad -rA=kCA E= 31400 H=-20,000 J/mol Balance de Energia: T=T0 + X[-H*T0]/Cp +XCp Cp=15 + 24 -40 =0 T=450 + 20,000X/40 = 450 + 500X k= 0.133 exp[E/R (1/450 -1/T)]

Abraham Castillo Espinosa

Pgina 11

b) Qu peso de catalizador es necesario para lograr una conversin del 80% en un CSTR? WCSTR = FA0X/-rA T=450+500X = 450+500(0.8) = 850K k= 0.133 exp[31400/8.314 (1/450 -1/850)]= 6.9 Abraham Castillo Espinosa Pgina 12

WCSTR =39.42 kg

c) Ahora tome en cuenta la cada de presin en un PFR dP/dW = -/2 *(T/T0)* P0/(P/P0) *(1+X) El reactor puede empacarse con partculas de dos tamaos distintos. Elija uno de ellos. =0.019 /kg cat para el dimetro de partcula D1 =0.0075/kg cat para el dimetro de partcula D2 Grafique la temperatura , la conversin y la presin a lo largo del reactor. Vare los parmetros a y P0 para obtener el rango de valores en los cuales se afecta dramticamente la conversin. dP/dW=-/2 (T/T0)P0/(P/P0) (1+X) CA=CA0(1-X/1+X)(T/T0)(P/P0)

Abraham Castillo Espinosa

Pgina 13

Problema 8-9 Use los datos los datos del problema 8-8 para el caso en el caso en que se retira calor mediante un intercambiador de calor en forma de chaqueta en el reactor. El flujo del enfriador por la chaqueta es suficientemente alto como para que la temperatura ambiente del intercambiador sea constante a Ta= 50 C a) 1. Grafique los perfiles de temperatura y conversin para un reactor empacado con Ua/pb =0.08 J/s*kg cat*K Donde : pb=densidad volumetrica granel del catalizador (kg/m3) a=rea del intercambio de calor por unidad de volumen del reactor (m2/m3) U=coeficiente total de transferencia de calor (J/s*m2*K) 2. Como varian los perfiles si Ua/pb fue incrementada por un factor de 3000? 3. Si hubiera cada de presin con =0.019 kg-1

Abraham Castillo Espinosa

Pgina 14

dT/dW = Ua/p (Ta-T)+(-rA)(-H)/FA0CpA -H=20,000 J/mol T=323 K CpA=40 J/mol/K

Abraham Castillo Espinosa

Pgina 15

Si UA/p se incrementa en un factor de 3000, las graficas se modifican de la siguiente manera

Abraham Castillo Espinosa

Pgina 16

b) Repita el inciso a) para flujo a co-corriente y a contra corriente con mc= 0.2 kg/s Cp=5000 J/kg K y una temperatura de entrada del enfriador de 50 C Para una temperatura no constante en la chaqueta, se incorpora una ecuacin de flujo a cocorriente: Ta0= 50C Introducimos una nueva ecuacin diferencial dTa/dW= Ua/p (T-Ta) /mCp

Abraham Castillo Espinosa

Pgina 17

Abraham Castillo Espinosa

Pgina 18

Para flujo contra corriente Taf= 50C dTa/dW= - Ua/p (T-Ta) /mCp con una Taf=323 K a W=50 Hacemos una manipulaciones en el valor inicial, hasta encontrar un valor de Ti=438.8163K

Abraham Castillo Espinosa

Pgina 19

c) Encuentre X y T para un CSTR fluidizado con 80 kg de catalizador UA = 500 J/s K pb=1 Kg/m3

Abraham Castillo Espinosa

Pgina 20

Balance de las especies XMB= k/1+ k , =(W/p) CA0 /FA0 Balance de Energa XMB=(UA/p)(T-Ta)/FA0 + CpA (T-T0) XEB= XMB Resolviendo tenemos X= 0.95, T= 323 d) Repita los incisos a) y b) para W= 80 kg asumiendo reaccin reversible con una velocidad especifica de reaccin de :

Abraham Castillo Espinosa

Pgina 21

Problema 9-10 El reactor con calormetro de bomba adiabtico ARSST tambin puede emplearse para determinar los rdenes de reaccin. La hidrlisis del anhdrido actico para formar cido actico se efectua adiabticamente A+B =2C Se postula una ley de velocidad de la forma -rA= kCAa CBb Se obtienen los siguientes datos de temperatura contra tiempo para dos concentraciones crticas distintas del anhdrido actico en operacin adiabtica. La velocidad de calentamiento era de 2C/min CA0 = 6.7M CB0= 0.2 M

Abraham Castillo Espinosa

Pgina 22

a) Asuma que Cp=0 y muestre para una conversin de X=1, la diferencia entre la teperatura final Tf y la temperatura inicial Ti

b) Muestre que la concentracin de A puede ser escrita como:

Abraham Castillo Espinosa

Pgina 23

Con X=1 , Tf= T

c) Muestre que en estado no estacionario el balance de energa puede ser escrito como:

Para operacin adiabtica Q=0 asumiendo Cp=0 , despreciando W

Abraham Castillo Espinosa

Pgina 24

d) Asuma de primer orden en A y en B y que B=3, entonces demuestre

Con a=1, b=1, B=3

e) Reescriba la Ecuacin P9-10.6 en la forma

f)

Obtenga la energa de activacin y la velocidad de reaccin especifica k.

y=a0+a1*T=-8944.2*T+ 17.534 -E/R E=-8944.2*8.314=-74.36kJ/mol

a0=17.534 k1=12100/6.7 exp(17.534)=1.298x1011

g) Encuentre el calor de reaccin X=1 T=Tf

Abraham Castillo Espinosa

Pgina 25

Abraham Castillo Espinosa

Pgina 26