Universidad Autnoma de
Yucatn
Facultad de Ingeniera Qumica
Ingeniera de reactores
IQI 6 Semestre
Grupo 1
Ejercicios de reactores
Integrantes:
Czarez Escalante Karla Mariana
Chan Catzn Liliana Jazmn
Lpez Falcn Alfonso
Maldonado Pinzn Cinthya Guadalupe
Noh Lpez Enrique Aldair
Ortiz Salazar Jos ngel
Quintal Flores Alan Israel
Rivera Arjona David Alejandro
Dr. ngel G. Zitlalpopoca Soriano
Jueves 26 de febrero de
2015
1- Considere la reaccin A+B P con k= 1 m3mol-1s-1. Suponga bin =10 mol/m3. Se
desea alcanzar una concentracin bout= 0.01 mol/ m3.
(a) Encuentre el tiempo de residencia necesario para alcanzar el objetivo
asumiendo un flujo pistn y con ain = bin.
%Ecuacin de diseo del reactor
dCi/dtresidencia=ri
function fx = ejercicio1 (t,x)
global k1
ca=x(1); cb=x(2); cp=x(3);
%
% reaction rates
%
ra=-k1*ca*cb;
rb=-k1*ca*cb;
rp=k1*ca*cb;
fx(1)=ra;
fx(2)=rb;
fx(3)=rp;
fx=fx';
%End main_ejercicio1.m
clear all; clc
%
%Multiple reactions in a tubular reactor
%
global k1
%Suponiendo los siguientes datos
cao=10; cbo=cao; cpo=0;
k1=1;
x0=[cao cbo cpo];
%Suposicion de un tiempo de residencia
time =linspace(0,99); %Es tiempo de residencia
[t,x]=ode15s('ejercicio1',time,x0);
ca=x(:,1); cb=x(:,2); cp=x(:,3);
%Presentar en una figura todos los resultados
figure (1)
plot(t,cb,'blue')
xlabel('Tiempo (seg)'), ylabel('Concentracion (mol/lt)')
legend('C_{B}')
t=99 s
%End main_ejercicio1.m
(b) Repita el inciso anterior asumiendo la reaccin ocurre en un RCTA.
clc;clear
%Ecuacin de diseo del reactor
%dCi=(q/v)*(cbo-cb)+ra
%Parametros
cbo=10;
cao=cbo;
cb=.01;
k=1;
%Relaciones auxiliares
xb=(1-cb/cbo);
ca=cao-cbo*xb;
ra=k*ca*cb;
rb=ra;
%Sabiendo que theta=v/q
theta=(cbo-cb)/(rb);
theta= 99900.0000000021 s
%ejercicio1b.m
(c) Repita a y b asumiendo ain = 10bin
%Repitiendo los cdigos de los incisos a y b pero cao=10*cbo obtenemos
que:
Para flujo pistn: t= 0.0756 s
Para RCTA: theta =11.0987668036885 s
2- La reaccin de trimetilamina con n-bromuro de propilo da una sal cuaternaria de
amonio:
N(CH3)3 + C3H7Br (CH3)3(C3H7)NBr
Suponga que los resultados de laboratorio a 110C usando tolueno como un solvente
muestra que la reaccin es de segundo orden con una constante de velocidad k=
5.610E-7m3mol-1s-1. Suponga [N(CH3)3]0 = [C3H7Br]0= 80 mol m-3.
(a) Estime el tiempo requerido para alcanzar 99% de conversin en un reactor
batch.
%Ecuacin de diseo del reactor
dt/dx=cao/ra
function dtdx=ejercicio2a(xa,t)% el tiempo depende de la conversin
porque integro la conversin.
%Parmetros
k=5.610e-7;
cao=80;
cbo=cao;
%Relaciones auxiliares
ca=cao*(1-xa);
cb=cbo*(1-xa);
cc=cao*(xa);
ra=k*ca*cb
%Ecuacin diferencial
dtdx=cao/ra;%expresin de diseo donde haremos la integral para
determinar el tiempo requerido xa
t=2222041.71549720 s
%End ejercicio2a.m
clear, clc
%Programa principal para obtener la simulacinde un reactor batch
t=0;%tiempo de simulacin
xa0=[0,0.99];%conversin deseada
[ys,ts]=ode15s('ejercicio2a',xa0,t);
%Graficacin
figure(1)
plot(ys,ts)
ylabel('Tiempo [seg]')
xlabel('Conversin')
%End main_ejercicio2a.m
(b) Estime el volumen requerido en un RCTA para alcanzar 99% de conversin si se
desea una produccin de sal de 100 kg h-1
%Ecuacin de diseo del reactor
dcadt=q*(cao-ca)/v-k*ca
clear all; clc
%Ecuacion de diseo
%dca/dt=q*(cao-ca)/v-k*ca*cb
%Parmetros
sal=100*1000/3600;%g/s
pmsal=182.1019;%g/mol
fa=sal/pmsal;%flujo alimentado mol/s
cbo=80;%mol/m3
cao=cbo;%mol/m3
k=5.610e-7;%m3/(mol*s)
xa=.99
%Relaciones auxiliares
ca=cao*(1-xa);%Concentracin mol/m3
q=fa/cao;%Flujo volumtrico alimentado m3/s
cb=cbo-(cao-ca);
%Clculo del volumen
V=q*(cao-ca)/(k*ca*cb);
V= 420605.920393484 m3
%End ejercicio2b.m
3- Las siguientes reacciones ocurren en un reactor batch isotrmico a volumen
constante:
A + BC k1
B + CD k2
Los parmetros para las reacciones son a0 = b0 = 10 mol/ m3, c0 = d0 = 0,
k1 = k2 = 0.01 m3mol-1h-1, y tbatch = 4 h.
(a) Encuentre la concentracin de C al final del ciclo batch.
%Ecuacin de diseo del reactor
dca/dt=ra
function fx = ejercicio3 (time,x)
global k1 k2
ca=x(1); cb=x(2); cc=x(3); cd=x(4);
%
% reaction rates
%
ra = -k1*ca*cb;
rb = -k1*ca*cb-k2*cb*cc;
rc = k1*ca*cb-k2*cb*cc;
rd = k2*cb*cc;
fx (1) = ra;
fx (2) = rb;
fx (3) = rc;
fx (4) = rd;
fx = fx';
%End ejercicio3.m
clear all; clc
%
% Multiple reactions in an isothermal batch reactor
%
global k1 k2
cao=10; cbo=10; cco=0; cdo=0;
k1=.01; k2=.01;
x0 = [cao cbo cco cdo];
time = linspace(0,4);
[t,x] = ode15s('ejercicio3',time,x0);
ca=x(:,1); cb=x(:,2); cc=x(:,3); cd=x(:,4);
%Presentar en una figura todos los resultados
figure (1)
subplot(221)
plot(t,ca,'red')
xlabel('Tiempo (hr)'), ylabel('Concentracion (mol/lt)')
legend('C_{A}')
subplot(222)
plot(t,cb,'blue')
xlabel('Tiempo (hr)'), ylabel('Concentracion (mol/lt)')
legend('C_{B}')
subplot(223)
plot(t,cc,'green')
xlabel('Tiempo (hr)'), ylabel('Concentracion (mol/lt)')
legend('C_{C}')
subplot(224)
plot(t,cd,'magenta')
xlabel('Tiempo (hr)'), ylabel('Concentracion (mol/lt)')
legend('C_{D}')
[C]=2.37355079327372 mol/lt
%End main_ejercicio3.m
4. Calcular la longitud del reactor tubular isotrmico que se requiere para consumir
90% del reactivo A si ocurre la reaccin elemental A +2B C. Se alimenta al reactor
una corriente de 0.2 m3/min de una mezcla de los reactivos A y B cuyas
concentraciones son 0.1 y 0.4 mol/m3, respectivamente. El rea transversal del reactor
es 0.1 m2 y la constante de velocidad de reaccin k es igual a 30 m6/(mol2-min).
%Ecuacin de diseo del reactor
dx/dv=-ra/fao
function dvdx = ejercicio4 (xa,v)
%Parmetros
k1=30;%constante de rapidez
cao=.1;%mol/m3
cbo=.4;%mol/m3
cco=0;%mol/m3
q=.2;%Flujo volumtrico alimentado m3/min
%Relaciones auxiliares
fao=q*cao;
ca=cao*(1-xa);
cb=cbo-2*cao*xa;
ra=k1*ca*cb*cb;
dvdx=fao/ra;
%End ejercicio4.m
clear, clc
%Programa principal para obtener la simulacin
%de un reactor tubular
v=0;%volumen de simulacin
xa0=[0,0.9];%conversin deseada
[ys,vs]=ode15s('ejercicio4',xa0,v);
%Graficacin
figure(1)
plot(ys,vs)
ylabel('Voluemn [m3]')
xlabel('Conversin')
V= 0.217084014188891 m3
%Sabiendo que V=Atransversal*l
L= 2.17084014188891 m
%End main_ejercicio4.m
a) Qu porcentaje del reactivo B alimentado se consume?
%Sabiendo que cb=cbo-2*cao*xa
cb= 0.22 mol/m3
xb =1-cb/cbo=45%
Por qu es necesario que CBo y CAo?
Para poder encontrar las concentraciones final de Cb y Ca mediante estequiometria.
Qu longitud se requerir si CBo= 0.3 mol/m3?, Qu porcentaje del reactivo B se
consume en este caso?
%Utilizando el mismo cdigo que para CBo=0.4 mol/m3 pero haciendo ahora
CBo=0.3 mol/m3
V=0.594176593571712 m3
L= 5.94176593571712 m
cb= 0.12 mol/m3
xb =1-cb/cbo=60%
Existir algn valor de CBo para el cual las conversiones de los reactivos A y B sean
altas simultneamente (p.e. 90% o ms)? Determinar este valor.
Cbo= 0.2 mol/m3
cb= 0.02 mol/m3
xb =1-cb/cbo=90%
