Saltar al contenido
Control Automático Educación

Modelado de un Reactor CSTR de Van de Vusse

En esta entrada aprenderás a Modelar y Linealizar un Reactor CSTR (Perfectamente agitado) usando sus ecuaciones diferenciales para describir el comportamiento dinámico y estacionario del sistema.

Adicionalmente aprenderás a simular el proceso del reactor CSTR en MATLAB y SIMULINK. Cuyos códigos y videos explicativos se encuentran al final de este post.

Antes de comenzar te invito a que veas nuestro CURSO GRATUITO DE MODELADO DE SISTEMAS.

Y que te suscribas al canal de YouTube para que sigas aprenduendo más sobre estos temas:

Reactor CSTR con Reacción de Van de Vusse

Vamos a estudiar el modelo fenomenológico de un reactor CSTR con reacción de Van de Vusse, el cual presenta una componente no lineal en su respuesta dinámica.

Aprovecharemos esta característica para obtener un modelo lineal usando las series de Taylor (matriz Jacobiana) que represente el proceso en diferentes puntos de operación. Este modelo lineal del reactor CSTR será escrito por medio de variables de estado y también por medio de funciones de transferencia.

Finalmente, realizaremos una comparación entre la respuesta dinámica del modelo del reactor CSTR no lineal vs la respuesta dinámica del reactor representado de forma lineal.

Reacción de Van de Vusse

El modelo CSTR (reactor perfectamente agitado, del ingles continuous stirred-tank reactor) de van de Vusse representa el sistema reacional utilizado para la síntesis de ciclopentenol (representado por B) a partir del ciclopentadieno (representado por A).

Este reactor envuelve entonces reacciones en serie y en paralelo entre la entrada del sistema (Reacción A) y el producto obtenido (Reacción B).

Sin embargo dado a la fuerte reatividade de los reagentes y produtos, se obtienen dos productos adicionales NO deseados, el diciclopentadieno (representado por D) es produzido como subproduto por la reacción de Diels-Alder y ciclopentenodiol (representado por C) es produzido como um produto consecutivo por la adición de otra molécula de água.

Ejemplo de reactor CSTR de Van de Vusse

Reactor de Van de Vusse
Reactor CSTR de Van de Vusse

A\overset{k_1}{\rightarrow}B\overset{k_2}{\rightarrow}C

2A\overset{k_3}{\rightarrow}D

donde los parámetros de la tasa de reacción, vienen dados por:

k_1=5/6 l/min

k_2=5/3 l/min

k_3=1/6 l/min

Balance de Masas del Reactor CSTR

Estudiemos las ecuaciones de diseño del reactor CSTR. Si se asume que la densidad y el volumen son constantes en el interior del reactor CSTR, podemos formular las siguientes dos hipotesis:

\dfrac{dV}{dt}=0

F\approx F_i

Con base a esto, procedemos a obtener el balance de cada uno de los componentes del reactor CSTR.

Balance del componente A:

Interior del reactor = Flujo de entrada – Flujo de Salida – sale reaccion 1 – sale reacción 2

\dfrac{d(VC_A)}{dt}=F(C_{Af}-C_A)-Vk_1C_A-VK_3C_A^2

Y debido a que consideramos el volumen constante, el balance de masa del componente A se resume a:

\dfrac{d(C_A)}{dt}=\dfrac{F}{V}(C_{Af}-C_A)-k_1C_A-K_3C_A^2

Balance del componente B:

De la misma forma podemos obtener el balance de masa para el componente B

\dfrac{d(C_B)}{dt}=-\dfrac{F}{V}C_B+k_1C_A-K_2C_B

Balance del componente C:

\dfrac{d(C_C)}{dt}=-\dfrac{F}{V}C_C+k_2C_B

Balance del componente D:

\dfrac{d(C_D)}{dt}=-\dfrac{F}{V}C_D+\dfrac{1}{2}k_3C_A^2

Como deseamos obtener el modelo del Componente B con relación al producto de entrada Componente A, observemos que estos dos balances de energia NO dependen de los balances C y D, por lo tanto, solo nos van a interesar estas dos ecuaciones diferenciales:

\dfrac{d(C_A)}{dt}=\dfrac{F}{V}(C_{Af}-C_A)-k_1C_A-K_3C_A^2

\dfrac{d(C_B)}{dt}=-\dfrac{F}{V}C_B+k_1C_A-K_2C_B

Estas dos ecuaciones representan el comportamiento fenomenológico no lineal de nuestro reactor de Van de Vusse.

Linealización del Reactor CSTR de Van de Vusse

Ahora como deseamos linealizar este sistema en torno a un punto de equilibrio o punto de operación el próximo paso sera encontrar dicho punto. Para lograr esto hacemos todas las derivadas (cambios) iguales a cero, lo que significa que el sistema no está cambiando y se encuentra en un estado estacionario.

Ecuaciones estacionarias del Reactor CSTR

A partir de la ecuación de C_A obtenemos lo siguiente:

0=\dfrac{F_s}{V}(C_{Afs}-C_{As})-k_1C_{As}-K_3C_{As}^2

-K_3C_{As}^2+\left(-k_1-\dfrac{F_s}{V}\right)C_{As} +\dfrac{F_s}{V}C_{Afs}=0

K_3C_{As}^2+\left(k_1+\dfrac{F_s}{V}\right)C_{As} -\dfrac{F_s}{V}C_{Afs}=0

Resolviendo la ecuación cuadrática anterior utilizando la formula general y escogiendo la raiz positiva (dado que una concentración negativa no tiene un sentido físico), tenemos que:

C_{As}=\dfrac{-\left(k_1+\dfrac{F_s}{V}\right)}{2k_3}+\dfrac{\sqrt{\left(k_1+\dfrac{F_s}{V}\right)^2+4k_3\dfrac{F_s}{V}C_{Afs}}}{2k_3}

Y para el componente B tenemos:

0=-\dfrac{F_s}{V}C_{Bs}+k_1C_{As}-K_2C_{Bs}

C_{Bs}=\dfrac{k_1C_{As}}{\dfrac{F_s}{V}+K_2}

La entrada de nuestro sistema es la tasa de disolución F/V l/min.

Si hacemos un barrido de nuestra entrada desde 1 l/min hasta 10 l/min. Obtendremos las curvas que relacionan la entrada y la salida del sistema en su estado estacionario:

Van de Vusse Estacionario


Si observamos el comportamiento del Componente B, vemos que tiene una inversión de fase (provocado por un cero de fase no mínima), que con una entrada de 1.292 l/min se obtiene la máxima producción del componente B (Lo más deseado) y que con una entrada superior, la producción de B comienza a caer.

A partir del gráfico anterior, podremos seleccionar un punto de operación, donde nosotros como ingenieros queremos que nuestro proceso opere.

Con ese punto de operación podemos proceder a realizar la linealización del sistema de las ecuaciones no lineales del reactor CSTR representado en Espacio de Estados utilizando el Jacobiano.

Representación Lineal del Reactor CSTR

Vamos a obtener una representación en espacio de estados de nuestro reactor.

\dot{x}=Ax+Bu

y=Cx

Definimos cuales serán nuestras variables de estado, en este caso serán las concentraciones de A y B:

x_1=C_A

\dot{x}_1=\dfrac{dC_A}{dt}

x_2=C_B

\dot{x}_2=\dfrac{dC_B}{dt}

u= F/V

Si renombramos las dos ecuaciones de balance de masa como función 1 (f_1) y función 2 (f_2).

\dfrac{d(C_A)}{dt}=f_1(C_A,C_B,F/V)=\dfrac{F}{V}(C_{Af}-C_A)-k_1C_A-K_3C_A^2

\dfrac{d(C_B)}{dt}=f_2(C_A,C_B,F/V)=-\dfrac{F}{V}C_B+k_1C_A-K_2C_B

Podremos obtener el jacobiano del sistema.

Jacobiano:

A=\begin{bmatrix}\dfrac{\partial f_1}{\partial x_1} & \dfrac{\partial f_1}{\partial x_2}\\ \\ \dfrac{\partial f_2}{\partial x_1} & \dfrac{\partial f_2}{\partial x_2}\end{bmatrix},B=\begin{bmatrix}\dfrac{\partial f_1}{\partial u} \\ \\\dfrac{\partial f_2}{\partial u}\end{bmatrix}

Matriz A:

A_{11}=\dfrac{\partial f_1}{\partial x_1}=-\dfrac{F_s}{V}-k_1-2k_3C_{As}

A_{12}=\dfrac{\partial f_1}{\partial x_2}=0

A_{21}=\dfrac{\partial f_2}{\partial x_1}=k_1

A_{22}=\dfrac{\partial f_2}{\partial x_2}=-\dfrac{F_s}{V}-k_2

Matriz B:

B_{11}=\dfrac{\partial f_1}{\partial u}=C_{Afs}-C_{As}

B_{21}=\dfrac{\partial f_2}{\partial u}=-C_{Bs}

El modelo en espacio de estado viene dado por:

\dot{x}=\begin{bmatrix}-\dfrac{F_s}{V}-k_1-2k_3C_{As} &0\\ \\k_1 & -\dfrac{F_s}{V}-k_2\end{bmatrix}x+\begin{bmatrix}C_{Afs}-C_{As} \\ \\-C_{Bs}\end{bmatrix}u

y=\begin{bmatrix}0&1\end{bmatrix}x

Teniendo la representación en el espacio de estado de nuestro sistema, vamos a seleccionar diferentes puntos de operación para analizar el comportamiento de este proceso.

Veremos tres casos diferentes. En el caso 1 veremos como aparece una respuesta inversa o respuesta de fase no mínima en la concentración de B mientras que en el caso 2 esta respuesta inversa no aparece. Y en el caso 3 se ilustra el punto de operación donde se alcanza la máxima producción del componente B.

Lazo Abierto y Lazo Cerrado

Lazo Abierto y Lazo Cerrado

Analisis de datos en MATLAB

Análisis de Datos en MATLAB

Reactor CSTR de Van de Vusse – CASO 1 (Fase no mínima)

Para este punto de operación vamos a considerar la tasa de flujo \dfrac{F_s}{V}=\dfrac{4}{7}min^{-1}. Notemos que esta tasa de flujo se ubica en el lado izquierdo del máximo pico del gráfico de estado estacionario que relaciona la entrada salida.

Inicialmente encontramos el estado estable del componente A (C_{As}) y el componente B(C_{Bs}).

C_{As}=\dfrac{-\left(k_1+\dfrac{F_s}{V}\right)}{2k_3}+\dfrac{\sqrt{\left(k_1+\dfrac{F_s}{V}\right)^2+4k_3\dfrac{F_s}{V}C_{Afs}}}{2k_3}

C_{As}=\dfrac{-\left(5/6+4/7\right)}{2(1/6)}+\dfrac{\sqrt{\left(5/6+4/7\right)^2+4(1/6)(4/7)(10)}}{2(1/6)}

C_{As}=3\dfrac{mol}{litro}

C_{Bs}=\dfrac{k_1C_{As}}{\dfrac{F_s}{V}+K_2}

C_{Bs}=\dfrac{3(5/6)}{4/7+5/3}

C_{Bs}=1.117\dfrac{mol}{litro}

Con esto obtenemos la siguiente representación en el espacio de estado

\dot{x}=\begin{bmatrix}-2.4048 &0\\ \\0.8333 & -2.2381\end{bmatrix}x+\begin{bmatrix}7 \\ \\-1.117\end{bmatrix}u

y=\begin{bmatrix}0&1\end{bmatrix}x

Y podemos obtener la función de transferencia de la misma forma que lo hemos visto en el curso de sistemas dinámicos lineales usando G(s)=C(sI-A)^{-1}B.

G(s)=\dfrac{-1.117s+3.1472}{s^2+4.6429s+5.3821}

Si notamos, vemos que esta función de transferencia tiene un cero positivo, lo que provoca la respuesta inversa o respuesta de fase no mínima.

Cero=\dfrac{3.1472}{1.117}=2.8175

Al final del post, puedes encontrar el código de implementación en Matlab para que lo copies y lo pegues en tu computador.

Veamos una comparación de los dos sistemas en el punto de operación linealizado

Van de Vusse caso 1

Ahora una comparación alejándonos más del punto de operación linealizado. Veamos como el sistema lineal comienza a dejar de valer.

Van de Vusse caso 1

Reactor CSTR de Van de Vusse – CASO 2 (Respuesta directa)

Para este punto de operación vamos a considerar la tasa de flujo \dfrac{F_s}{V}=2.8744min^{-1}. Notemos que esta tasa de flujo se ubica en el lado derecho del máximo pico del grafico de estado estacionario que relaciona la entrada salida.

Procedemos a encontrar el punto estable de ambos componentes de la misma forma que lo hicimos en el caso 1:

C_{As}=6.0870\dfrac{mol}{litro}

C_{Bs}=1.117\dfrac{mol}{litro}

Con esto obtenemos la siguiente representación en el espacio de estado

\dot{x}=\begin{bmatrix}-5.7367 &0\\ \\0.8333 & -4.5411\end{bmatrix}x+\begin{bmatrix}3.9130 \\ \\-1.117\end{bmatrix}u

y=\begin{bmatrix}0&1\end{bmatrix}x

Y podemos obtener la función de transferencia de la misma forma que lo hemos visto en el curso de sistemas dinámicos lineales usando G(s)=C(sI-A)^{-1}B.

G(s)=\dfrac{-1.117s-3.1472}{s^2+10.2778s+26.0508}

En este caso ya se tiene el cero en el semiplano negativo, por lo tanto la respuesta inversa desaparece.

Cero=-\dfrac{3.1472}{1.117}=-2.8175

Van de Vusse caso 2

Reactor CSTR de Van de Vusse – CASO 3 (Punto óptimo)

El punto de operación óptimo para la producción del componente B se logra utilizando la siguiente tasa de disolución:

\dfrac{F_s}{V}=1.2921min^{-1}

C_{As}=4.4949\dfrac{mol}{litro}

C_{Bs}=1.266\dfrac{mol}{litro}

Con esto obtenemos la siguiente representación en el espacio de estado

\dot{x}=\begin{bmatrix}-3.6337 &0\\ \\0.8333 & -2.9588\end{bmatrix}x+\begin{bmatrix}5.5051 \\ \\-1.2660\end{bmatrix}u

y=\begin{bmatrix}0&1\end{bmatrix}x

Y podemos obtener la función de transferencia de la misma forma que lo hemos visto en el curso de sistemas dinámicos lineales usando G(s)=C(sI-A)^{-1}B.

G(s)=\dfrac{-1.2660s}{s^2+6.5825s+10.7217}

El cero de esta función de transferencia se ubica en el origen. En este caso la condición de estado estable no cambia en la concentración de B para pequeños cambios en la entrada. Esto es debido a la ganancia cero de la función de transferencia.

Van de Vusse caso 3

Suscríbete a este sitio WEB para estar enterado de las nuevas entradas!

Simulación de un Reactor CSTR

Aprende como simular y modelar un Reactor CSTR utilizando el MATLAB y SIMULINK detalladamente.

A continuación te dejo dos videos explicando en detalle el procedimiento de simulación de estos dos procesos y al final podrás descargar los códigos implementados en los videos.

MATLAB:

SIMULINK:

Si estás interesado en aprender a utilizar perfectamente el SIMULINK para el modelado, simulación y control de procesos industriales, te invito a que te inscribas al curso de SIMULINK desde CERO con más de 10 horas de videos explicando en detalle este potente software de ingeniería:

A continuación se presenta el código de implementación para que lo puedas copiar y pegar en tu computador. Recuerda compartir el contenido de esta entrada en tus redes sociales, asi ayudas a que más personas conozcan la página y se beneficien con esta información.

%% Sergio Castaño
% Rio de Janeiro - 2017
% https://controlautomaticoeducacion.com/
%%
clc
close all
clear all
%Constantes
k1 = 5/6;
k2 = 5/3;
k3 = 1/6;
cafs = 10; %Concentracion de A en la entrada del reactor
fsov = 0:0.1:10;
%Vector de Parametros
par.k1=k1; par.k2=k2; par.k3=k3; par.cafs=cafs; par.fsov=fsov;
nm=100;
Ts=0.1;
t = 0:Ts:(nm-1)*Ts; %Tiempo de Simulación
%% PUNTOS DE EQULIBRIO
% Los puntos de equilibrio pueden ser encontrados utilizando cualquiera de
% las dos siguientes posibilidades:
%% Punto de Equilibrio (Posibilidad 1)
cas1 = -(k1 + fsov)/(2*k3);
cas2 = sqrt((k1+fsov).^2 + 4*k3*fsov.*cafs)/(2*k3);
cas  = cas1 + cas2;
cbs = k1*cas./(fsov + k2);
%% Punto de Equilibrio (Posibilidad 2)
x0=[1,1]; %Valores iniciales
for i=1:length(fsov)
%Encuentra el estado Estacionario
par.fsov=fsov(i);
X(i,:) = fsolve(@(x)vvmodel(t,x,par),x0); %X=[Ca Cb]
end
% Grafica de los Puntos de Equilibrio
figure
subplot(2,1,1)
plot(fsov,X(:,1),'Linewidth',2)
title('Relacion Concentraciones vs F/V');
ylabel('Concentracion C_{As} mol/litro');
legend('C_{As}','Location','SouthEast')
subplot(2,1,2)
plot(fsov,X(:,2),'Linewidth',2)
ylabel('Concentracion C_{Bs} mol/litro');
legend('C_{Bs}','Location','SouthEast')
xlabel('F_s/V (min^{-1})');
%% Linealizacion del sistema por Jacobiana
syms x1 x2 u
% Ecuaciones diferenciales
% x1 = ca
% x2 = cb
% u  = fsov
dx1 = u*(cafs - x1) -k1*x1 - k3*x1^2;
dx2 = -u*x2 + k1*x1 - k2*x2;
% Funciones para realizar la jacobiana
fx1 = u*(cafs - x1) -k1*x1 - k3*x1^2;
fx2 = -u*x2 + k1*x1 - k2*x2;
%Conforma vectores
fx = [fx1;fx2];
x = [x1;x2];
%Linealiza por medio de la Jacobiana
A = jacobian(fx,x);
B = jacobian(fx,u);
C = [0 1]; %Matriz de salida (Cb)
%***********************************************************************
%***********************************************************************
%% Caso 1 (Respuesta inversa)
% Punto de equilibrio
par.fsov=4/7;
Eq1 = fsolve(@(x)vvmodel(t,x,par),x0); %Eq1=[Ca Cb]
%Reemplazo los puntos de equilibrio en la Jacobiana
display('Representación en Espacio de Estado Caso 1')
A1 = double(subs(A,{x1,u},{Eq1(1),par.fsov}))
B1 = double(subs(B,{x1,x2},{Eq1}))
%Determino la Funcion de Transferencia
display('Representación en Función de Transferencia Caso 1')
[num1,den1]=ss2tf(A1,B1,C,0);
ft1=tf(num1,den1)
% Comparación del Modelo Lineal vs Modelo No Lineal
du=0.01;   %Pequeña Variación en la entrada del sistema
par.fsov=4/7+du;
% Simula el Modelo No Lineal
[ts,X1] = ode45(@(t,x)vvmodel(t,x,par), t , Eq1, par.fsov);
%tsm=0:0.001:0.1;
u2=zeros(length(ts),1); %Creo el vector de entrada para la FT
u2(:,1)=du*ones(length(ts),1); %Lo lleno con el pequeño incremento
[ylin]=lsim(ft1,u2,ts);
%Grafica
figure
plot(ts,X1(:,2),'-b',ts,ylin+Eq1(2),'--r','LineWidth',2);
title('Comparación modelo Lineal vs No Lineal (caso 1)');
xlabel('tiempo (min)');
ylabel('C_b mol/litro');
legend('No Lineal','Lineal','Location','SouthEast')
%***********************************************************************
%***********************************************************************
%% Caso 2 (Respuesta NO inversa)
% Punto de equilibrio
par.fsov=2.8744;
Eq2 = fsolve(@(x)vvmodel(t,x,par),x0); %Eq2=[Ca Cb]
%Reemplazo los puntos de equilibrio en la Jacobiana
display('Representación en Espacio de Estado Caso 2')
A2 = double(subs(A,{x1,u},{Eq2(1),par.fsov}))
B2 = double(subs(B,{x1,x2},{Eq2}))
%Determino la Funcion de Transferencia
display('Representación en Función de Transferencia Caso 2')
[num2,den2]=ss2tf(A2,B2,C,0);
ft2=tf(num2,den2)
% Comparación del Modelo Lineal vs Modelo No Lineal
du=0.1;   %Pequeña Variación en la entrada del sistema
par.fsov=2.8744+du;
% Simula el Modelo No Lineal
[ts,X2] = ode23t(@(t,x)vvmodel(t,x,par), t , Eq2, par.fsov);
%tsm=0:0.001:0.1;
u2=zeros(length(ts),1); %Creo el vector de entrada para la FT
u2(:,1)=du*ones(length(ts),1); %Lo lleno con el pequeño incremento
[ylin]=lsim(ft2,u2,ts);
%Grafica
figure
plot(ts,X2(:,2),'-b',ts,ylin+Eq2(2),'--r','LineWidth',2);
title('Comparación modelo Lineal vs No Lineal (caso 2)');
xlabel('tiempo (min)');
ylabel('C_b mol/litro');
legend('No Lineal','Lineal','Location','SouthEast')
%***********************************************************************
%***********************************************************************
%% Caso 3 (Punto Optmino)
% Punto de equilibrio
par.fsov=1.2921;
Eq3 = fsolve(@(x)vvmodel(t,x,par),x0); %Eq3=[Ca Cb]
%Reemplazo los puntos de equilibrio en la Jacobiana
display('Representación en Espacio de Estado Caso 3')
A3 = double(subs(A,{x1,u},{Eq3(1),par.fsov}))
B3 = double(subs(B,{x1,x2},{Eq3}))
%Determino la Funcion de Transferencia
display('Representación en Función de Transferencia Caso 3')
[num3,den3]=ss2tf(A3,B3,C,0);
ft3=tf(num3,den3)
% Comparación del Modelo Lineal vs Modelo No Lineal
du=0.01;   %Pequeña Variación en la entrada del sistema
par.fsov=1.2921+du;
% Simula el Modelo No Lineal
[ts,X3] = ode23t(@(t,x)vvmodel(t,x,par), t , Eq3, par.fsov);
%tsm=0:0.001:0.1;
u2=zeros(length(ts),1); %Creo el vector de entrada para la FT
u2(:,1)=du*ones(length(ts),1); %Lo lleno con el pequeño incremento
[ylin]=lsim(ft3,u2,ts);
%Grafica
figure
plot(ts,X3(:,2),'-b',ts,ylin+Eq3(2),'--r','LineWidth',2);
title('Comparación modelo Lineal vs No Lineal (caso 3)');
xlabel('tiempo (min)');
ylabel('C_b mol/litro');
legend('No Lineal','Lineal','Location','SouthEast')

Función del Modelo:

  function x = vvmodel(~,x,par)
%  Modelo dinamico del Reactor de Van de Vusse
%  x    = Estados
%  x(1) = concentración de A
%  x(2) = concentración de B
%  t    = tiempo
%
%  fov  = Tasa de disolución
%  ca   = concentración de a
%  cb   = concentración de b
%  caf  = Concentración de A en la entrada del reactor
%
%  Esquema de la Reacción
%
%  A ---> B ---> C
%    (k1)   (k2)
% 2A ---> D
%    (k3)
%
k1  = par.k1; % (min^-1)
k2  = par.k2; % (min^-1)
k3  = par.k3; % (mol/(liter min)
caf = par.cafs; %  mol/liter
fov = par.fsov; % (min^-1)
%Definición de las Variables de Estado
ca   = x(1);
cb   = x(2);
%
%  Sistema de Ecuaciones
%
dcadt = fov*(caf - ca) -k1*ca - k3*ca*ca;
dcbdt = -fov*cb + k1*ca - k2*cb;
%
%  Vector de Estados
%
x = [dcadt;dcbdt];

Bibliografia

  • Bequette, B. (1998). Process Dynamics: Modeling, Analysis, and Simulation. New Jersey: Prentice Hall

Eso es todo por la entrada del dia de hoy, espero les haya gustado y hayan aprendido algo nuevo. Si te ha servido el contenido de esta entrada, de los videos y los códigos de implementación y deseas apoyar mi trabajo invitandome a un café super barato, puedes hacerlo en el siguiente link:

👉 Invitar a Sergio a un Café ☕️

Que esten muy bien, nos vemos en la siguiente entrada.

Entradas relacionadas

Deja una respuesta

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *

Este sitio usa Akismet para reducir el spam. Aprende cómo se procesan los datos de tus comentarios.

Comentarios (2)

oi Sérgio td bem? fiz a disciplina do Maurício com vc esse semestre na Ufrj mas aluna da puc. Eu estava conferindo seu código com o meu pois tive alguns erros nos meus trabalhos e queria entender melhor algo: no seu código, o k1, k2 e k3 são constantes, certo? isso é, não variam com a temperatura como no exercício da disciplina, certo? ou esses k1, k2 e k3 do seu código se referem aos k1, k2 e k3 iniciais?

Responder

Olá Sarah, tudo certo!!. É correto teu analise. Neste exemplo aqui, são consideradas as constantes cineticas K1, K2, e K3 como constantes que são aqualas que aparecem no código e não variam com a Temperatura. O modelo do exemplo desta entrada, não está conseiderando-se o balaço de energia da jaqueta.
Na outra entrada, está o reator de VV para o caso MIMO, que é igual ao visto na disciplina, ali sim consideramos a jaqueta e as constantes variam com relação à temperatura. Podes dar click aqui, caso tu não tenhas visto o outro exemplo. Qualquer dúvida é só me procurar. Grande abraço.

Responder