Reactor PFR isotérmico

REACTOR PFR ISOTÉRMICO

%Sistema Isotermico para un PFR

clear,clc

global h

%Datos

CA0 = 0.03;      %mol/L

v0 = 100;        %L/min

T0 = 288.15;     %K

Cpa = 3.8;       %KJ/Kg-K

Cpb = 4.184;     %KJ/Kg-K

FA0 = CA0*v0;

%Datos cineticos para k

A = 21400000;    %min^-1

E = 46.5;        %KJ/mol

dHRX = -209;     %KJ/mol

RSC = 8.3144;    %J/mol-k

R = RSC/1000;    %KJ/mol

Te = 0;

%Rango de integracion

%Condiciones iniciales de la variable independiente

Xi = 0;

Xf = 0.8;

%Tamaño de integracion de Metodo de Euler

h = 0.0001;

%Generacion de vector tiempo

X= Xi:h:Xf;

%Numero de iteraciones

n = length(X);

%Condiciones iniciales

V=0;

Conv = Xi;

%Inicio de Metodo de Euler (Desde la Iteracion 2)

for i =1:n-1

    %Funciones adicionales

    rA = -(A*exp(-E/(R*T0)))*(CA0*(1-X));

    %Ecuaciones de diseño

    dV_dX = FA0./-rA;

    %Llamar a la funcion INT_EULER para ocupar el metodo iterativo

    V = INT_EULER(dV_dX, V);

end

%Reeordenamiento de los valores

XX=X';

VV=V';

%Variables de interes

k = (A*exp(-(E ./(R*T0))));

rA = -(k.*(CA0*(1-X)));

%Tabla de resultados

resultados = [VV,XX]

%Graficas de conversión vs variables de interes

subplot(2,1,1),plot(X, V,'r'),title('Volumen vs Conversion')

xlabel('Conversion'),ylabel('Volumen')

legend({'Vol'},'Location','northeast')

subplot(2,1,2),plot(X, -rA,'g'),title('Conversion vs -rA')

xlabel('Conversion'),ylabel('-rA')

legend({'-rA'},'Location','northeast')

[pic 1]

REACTOR PFR ADIABÁTICO

#Datos

CA0 = 0.03

v0=100

T0 = 288.15

Cpa=3.8

Cpb = 4.184

Fa0=CA0*v0

#Operación de la temperatura

dHRX = -209

Te = 0

Op = Cpa+(Te*Cpb)

T = ((-dHRX*X)/Op)+T0

#Datos cinéticos

A = 21400000

E = 46.5

R=8.3144/1000

k = A*exp(-E/(R*T))

rA= -k*(CA0*(1-X))

#Ecuación de volumen vs conversión

d(V) / d(X) = Fa0/(-rA) 

V(0) = 0

#Rango de Conversión a analizar

X(0)=0

X(f)=0.8

[pic 2]

[pic 3]

[pic 4]

REACTOR CSTR ISOTÉRMICO

%Sistema Isotérmico en un CSTR

clear,clc

global h

%Datos

CA0 = 0.03;      %mol/L

v0 = 100;        %L/min

T0 = 288.15;     %K

Cpa = 3.8;       %KJ/Kg-K

Cpb = 4.184;     %KJ/Kg-K

...