CSTR No Isotermico

REACTORES CSTR EN ESTADO ESTACIONARIO
1. INTRODUCCION
En el diseño de reactores continuos tipo tanque con agitación perfecta (CSTR) se requiere resolver las
ecuaciones del balance de materia y energía.
Se toma como referencia una reacción general con la siguiente estequiometria:
b c d
A + B→ C + D [1]
a a a
2. BALANCE DE MATERIA
La ecuación del balance de materia para un CSTR en estado estacionario puede tomar diversas formas:
0=F Ao −F A +r A V
Recordando que la diferencia del flujo molar que entra y el flujo molar que sale del tanque se debe a
los moles convertidos.
F Ao −F A=F Ao X
El balance de materia se puede escribir como:
0=F Ao X +r A V
De la cual se puede despejar el volumen del reactor:
F Ao X
V=
−r A [2]
A la cual se le denomina ecuación de diseño.
Para poder calcular el volumen es necesario especificar la cinética de la reacción, es decir determinar el
orden de reacción y su dependencia con la temperatura, usualmente con la ecuación de Arrhenius
r A =kϕ( X , T ) [3]
[4]
k =k o e−E /RT
3. BALANCE DE ENERGIA
La ecuación general de energía en estado estacionario para los reactores continuos (CSTR y PFR) es:
n T
0=−F A 0 ∑ θi ∫ Cpi dT −ΔH R(T ) F A 0 X + Q̇−Ẇ [4]
i=1 T io
FLUJO MOLAR RELATIVO: θ
Fi 0 Cio v 0 Cio Y io
θi = = = =
F Ao C Ao v 0 C Ao Y Ao
ENTALPIA DE REACCION:
nr T nr [5]
i=1
o
ΔH R (T )=∑ υ i H i( T )+ ∫
R [
T R i=1
]
∑ υ i Cpi dT
PARA CAPACIDADES CALORIFICAS CONSTANTES:
_
H R (T)  H oR (TR )   Cp  T  TR 
PARA CAPACIDADES CALORIFICAS VARIABLES:
Cpi=ai +bi T +cT 2i
H R (T)  H oR (TR )   (a  bT  cT 2 )dT

TR
4. CASO DE ESTUDIO: PRODUCCIÓN DE GLICOL PROPILENO
Se tomará como ejemplo el caso desarrollado por Fogler S. en su texto: Elements of Chemical Reaction
Engenieering. 2nd edition para la producción de glicol propileno a partir del óxido de propileno.
A temperatura ambiente y ácido sulfúrico como catalizador, la reacción es:

H2SO4
CH2 CH CH3 H2O CH2 CH CH3
O OH OH
Oxido de propileno agua propilenglicol
Se debe determinar la conversión y temperatura de reacción para un reactor CSTR de 40 lt que opera
adiabáticamente, para lo cual se mezcla agua con óxido de propileno según el diagrama siguiente:
A, B, I
To = 75°F
X, T
V = 300 gal
A, B, I, C
 Datos de la reacción:
Estequiometría:
A + B→C
A: oxido de propileno
B: agua
C: glicol propileno
M: Metanol (inerte)
Orden de Reacción: primer orden
−r A =kC A
C A=C Ao (1− X)
Como dato se tiene que: C Ao=¿0.132
Ecuación de Arrhenius:
k =k o e− E / RT
-E = -32400 ko=16.96 x 10 12 R=1.987
 Datos de la alimentación
Flujos molares de entrada:
F Ao=40.3
F Bo=802.8
F Io =71.87
Temperatura de entrada:
T o=535 ° R
 Capacidades Caloríficas y Entalpias Estándar
A B C I
Cp 35 18 46 19.5
H oi -66600 -123000 -226000
Tref = 528°R
BALANCE DE MATERIA
τ k o e−E / RT
x= − E /RT
1+τ k o e
Si se considera que:
k =k o e− E / RT
Entonces:
τk
x=
1+τk
Para usar fsolve es necesario que la ecuación se iguale a cero
x ( 1+ τk )−τk =0
−E −E
(
x 1+ τ k o e RT
)−τ k e
o
RT
=0
V
τ=
vo
V = 40.1
vo = 326.34
BALANCE DE ENERGIA
Considerando las capacidades caloríficas constantes:

n ¿
∑ θi Cpi (T −T io )+ ΔH R (T ) X =0
i=1
ARCHIVOS DE MATLAB
Script: Archivo GLICOL.m
%========================================
% REACTOR CSTR NO ADIABATICO
% Cp CONSTANTES.
% ESTADO ESTACIONARIO
% CASO: GLICOL PROPILENO
% NOMENCLATURA
% A: OXIDO DE PROPILENO
% B: AGUA
% C: GLICOL PROPILENO
% M: METANOL
%======================================
% DATOS DEL REACTOR
V = 40.1; % VOLUMEN DEL REACTOR
vo = 326.34; % FLUJO VOLUMETRICO
tau = V/vo
% DATOS DE LA CORRIENTE DE ALIMENTACION

FAo= 43.04;
FBo = 802.2;
FMo = 71.87;
CAo = 0.132;
To = 535;
% DATOS PARA LA ECUACION DE ARRHENIUS

A = 16.96e12;
E = 32400;
R = 1.987;
% DATOS DE CAPACIDADES CALORIFICAS (CONSTANTES)

CPA= 35;
CPB = 18;
CPC = 46;
CPM = 19.5;
% DATOS DE ENTALPIAS ESTANDAR

HEA = -66600;
HEB = -123000;
HEC = -226000;
TR = 528;
% CALCULO DE LOS COEFICIENTES TITA

TITA = 1;
TITB = FBo/FAo;
TITM = FMo/FAo;
% CALCULO DEL CALOR DE REACCION ESTANDAR (DHE) Y DELTA CP (DCP)

DHE = HEC - HEA - HEB;
DCP = CPC - CPA - CPB;
% CALCULO DE SUMATORIA DE TITAi*CPi

STCP = TITA*CPA + TITB*CPB + TITM*CPM;
% CALCULO DE X Y T DE MANERA SIMULTANEA
VI =[0.1,800];
S= fsolve(@(Z)BEM1(Z,A,R,CAo,E,To,DHE,DCP,STCP,TR,tau), VI);
S
%==========================
Archivo de función
%======== ARCHIVO DE FUNCIÓN ==================
function F = BEM1( Z,A,R,CAo,E,To,DHE,DCP,STCP,TR,tau)
X = Z(1);
T = Z(2);
K = A*exp(-E/(R*T));
RA = -K*CAo*(1-X);
DHR = DHE+ DCP*(T-TR);
F(1) = X*(1+tau*K)-tau*K;
F(2) = STCP*(T-To)+ DHR*X;
end
%======================================
RESULTADOS
Al ejecutar el script se obtienen los resultados para la fracción de conversión es para X = 0.8578 y la
temperatura T = 613.7499
Los resultados coinciden con los presentados por Fogler.
ACTIVIDADES
ANALISIS DEL EFECTO DE VARIABLES
Se podrían tomar como independientes a ciertas variables del reactor y de la alimentación. Por ejemplo
podría preguntarse qué pasaría si se cambia el tiempo espacial Tau ( τ = V/vo). También se podría
analizar el efecto de la temperatura de la corriente de alimentación y su composición.
a) Utilizando el script complete la tabla siguiente:
Volumen de Flujo Tau Conversión Temperatura

reactor volumétrico de (X) (T)
(V) Alimentación
(vo)
Grafique X y T versus tau. Interprete sus resultados y obtenga conclusiones.
b) De igual modo evalúe el efecto de la temperatura de la alimentación sobre la conversión y

temperatura de reacción. Complete la tabla, analice los resultados y obtenga conclusiones.
Temperatura de Conversión Temperatura

alimentación (To) (X) (T)
c) Los análisis anteriores se pueden realizar ejecutando sucesivas veces el script anterior, pero
también se le podría modificar para que realice varios cálculos de una sola vez. Modifique el script
presentado para que realice varios análisis de variables.

CSTR No Isotermico

Cargado por

Copyright:

Formatos disponibles

CSTR No Isotermico

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

CSTR No Isotermico

Cargado por

Copyright:

Formatos disponibles

REACTORES CSTR EN ESTADO ESTACIONARIO

Se toma como referencia una reacción general con la siguiente estequiometria:

A la cual se le denomina ecuación de diseño.

FLUJO MOLAR RELATIVO: θ

PARA CAPACIDADES CALORIFICAS CONSTANTES:

PARA CAPACIDADES CALORIFICAS VARIABLES:

Cpi=ai +bi T +cT 2i

H R (T)  H oR (TR )   (a  bT  cT 2 )dT

4. CASO DE ESTUDIO: PRODUCCIÓN DE GLICOL PROPILENO

A temperatura ambiente y ácido sulfúrico como catalizador, la reacción es:

Orden de Reacción: primer orden

Como dato se tiene que: C Ao=¿0.132

Flujos molares de entrada:

 Capacidades Caloríficas y Entalpias Estándar

Para usar fsolve es necesario que la ecuación se iguale a cero

Considerando las capacidades caloríficas constantes:

Script: Archivo GLICOL.m

% DATOS DE LA CORRIENTE DE ALIMENTACION

% DATOS PARA LA ECUACION DE ARRHENIUS

% DATOS DE CAPACIDADES CALORIFICAS (CONSTANTES)

% DATOS DE ENTALPIAS ESTANDAR

% CALCULO DE LOS COEFICIENTES TITA

% CALCULO DEL CALOR DE REACCION ESTANDAR (DHE) Y DELTA CP (DCP)

% CALCULO DE SUMATORIA DE TITAi*CPi

% CALCULO DE X Y T DE MANERA SIMULTANEA

%======== ARCHIVO DE FUNCIÓN ==================

function F = BEM1( Z,A,R,CAo,E,To,DHE,DCP,STCP,TR,tau)

DHR = DHE+ DCP*(T-TR);

Los resultados coinciden con los presentados por Fogler.

ANALISIS DEL EFECTO DE VARIABLES

a) Utilizando el script complete la tabla siguiente:

Volumen de Flujo Tau Conversión Temperatura

Grafique X y T versus tau. Interprete sus resultados y obtenga conclusiones.

b) De igual modo evalúe el efecto de la temperatura de la alimentación sobre la conversión y

Temperatura de Conversión Temperatura

También podría gustarte