Scribd
Cargar
57 vistas5 páginas

ASIGNACION-ejemplo-4.9-matlab-1 1

El documento presenta un estudio sobre la producción de bromuro de metilo en un reactor semicontinuo isotérmico, describiendo la reacción química y las condiciones iniciales. Se incluyen ecuaciones diferenciales para modelar la variación de concentraciones de reactivos y productos a lo largo del tiempo, así como un código MATLAB para simular el sistema. Se busca calcular las concentraciones de los reactivos y productos, así como la velocidad de reacción en función del tiempo.

Cargado por

vanesa
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
57 vistas5 páginas

ASIGNACION-ejemplo-4.9-matlab-1 1

El documento presenta un estudio sobre la producción de bromuro de metilo en un reactor semicontinuo isotérmico, describiendo la reacción química y las condiciones iniciales. Se incluyen ecuaciones diferenciales para modelar la variación de concentraciones de reactivos y productos a lo largo del tiempo, así como un código MATLAB para simular el sistema. Se busca calcular las concentraciones de los reactivos y productos, así como la velocidad de reacción en función del tiempo.

Cargado por

vanesa
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

Integrantes: Jose Contreras,Mario Robledo.

Diseño de reactores.
Universidad del Atlántico. Facultad de Ingeniería. Ingeniería Química
Fecha: 9 de Nov del 2017
Grupo 2.
EJEMPLO 4,9. REACTOR SEMICONTINUO ISOTÉRMICO.
Elementos de ingeniería de las reacciones químicas. H. Scott Fogler
La producción de bromuro de metilo es una reacción irreversible en fase líquida que sigue
una ley de velocidad elemental. La reacción

CNBr +C H 3 N H 2 → C H 3 Br + NCN H 2

Se efectúa isotérmicamente en un reactor semicontinuo. Se alimenta una solución acuosa de


metilamina (B) con una concentración de 0.025 mol /dm3 con flujo volumétrico de
3
ϑ 0=0.025 dm /s a una solución acuosa de cianuro de bromo (A) contenida en un reactor
recubierto de vidrio. El volumen inicial del líquido en el recipiente debe ser de 5 dm3 con
una concentración de cianuro de bromo de 0.05 mol /dm 3. La constante de velocidad
específica de la reacción es:
3
k =2.2 dm /(s∗mol)
Encuentre las concentraciones de cianuro de bromo y bromuro de metilo, así como la
velocidad de reacción en función del tiempo

SOLUCIÓN
Datos
Constante de velocidad especifica de la reacción k =2.2 dm3 /(s∗mol)
Flujo volumétrico inicial ϑ 0=0.025 dm3 /¿
Concentración inicial de B C B=0.025 mol/dm3
Volumen inicial V 0=5 dm3
Concentración inicial de A C A=0.05 mol/dm3

La reacción a estudiar es:


CNBr +C H 3 N H 2 → C H 3 Br + NCN H 2

A+ B →C + D

Ley de velocidad (elemental): −r a =k C A C B (Ec.1)


dCA ϑ0
Balance molar de la especie A: =r a− C A (Ec. 2)
dt V

d CB ϑ0
Balance molar de la especie B: =r a + (C ¿ ¿ BO−C B )¿ (Ec. 3)
dt V

d CC ϑ0
Balance molar de la especie C: =−r a − C C (Ec. 4)
dt V

dCD ϑ0
Balance molar de la especie D: =−r a − C D (Ec. 5)
dt V

dV
Del balance global de masa, para un sistema a densidad constante =ϑ 0 (Ec. 6)
dt

Con la condición inicial V = V0 en t = 0, e integrando para el caso de flujo volumétrico


constante ϑ 0 se obtiene

V =V 0+ ϑ 0∗t (Ec. 7)

Sustituyendo la (Ec. 1) en (Ec. 2; 3; 4; 5) se obtienen las ecuaciones diferenciales a


ingresar en el programa MATLAB.

dCA ϑ0
=−k C A C B − C A (Ec. 8)
dt V

d CB ϑ0
=−k C A C B + (C ¿ ¿ BO−C B )¿ (Ec. 9)
dt V

d CC ϑ0
=k C A C B− C C (Ec. 10)
dt V
dCD ϑ0
=k C A C B − C D (Ec. 11)
dt V

También podríamos calcular la conversión de A.


N A 0−N A
X=
NA0
C A 0 V 0−C A V 0
X=
C A0V 0

Codigo
%Ejemplo 4,9 Reactor semicontinuo%
%Variacion de las concentraciones de reactivos y productos con respecto
al tiempo
%CA=x(1)
%CB=x(2)
%CC=x(3)
%CF=x(4)
Eq=inline('[-2.2*x(1)*x(2)-(0.05/5)*x(1);-2.2*x(1)*x(2)+(0.05/5)*(0.025-
x(2));2.2*x(1)*x(2)-(0.05/5)*x(3);2.2*x(1)*x(2)-
(0.055/5)*x(4)]','t','x');
[t,x]=ode15s(Eq,[0:10:500],[0.05 0 0 0]);
disp([t,x])
figure (1)

plot(t,x(:,1),'-b',t,x(:,2),'-r',t,x(:,3),'-
g',t,x(:,4),'y','LineWidth',2),grid on

title('Variación de concentraciones','Fontsize',15);
text(50,0.03,'CA');text(250,0.022,'CB');text(250,0.005,'CC');text(150,0.0
1,'CD');
xlabel('Tiempo');ylabel('Concentración');
figure(2)
r=2.2*x(:,1).*x(:,2);
plot (t,r)
title(' Velocidad de reaccion vs tiempo');
xlabel('tiempo');ylabel('velocidad de reaccion');

También podría gustarte