Flujo de Flujo de Aparece dentro Desaparece dentro
Variación de la masa que masa que del Sistema por del Sistema por
masa dentro del entra al sale del crecimiento muerte
Sistema Sistema dm dm
sistema,
dm
x 0 , F
dt x1 , F dt G dt D
de la cual se obtiene la siguiente ecuación:
dm dm
x0 F x1 F (1)
dt G dt D
m x1V (2)
Vdx1 dm dm
x0 F x1 F V V (3)
dt dt G dt D
donde:
dm
G x1 (4)
dt G
dm
D x1 (5)
dt D
Vdx1
x0 F x1 F V G x1 V D x1 (6)
dt
dx1 F F
x0 x1 G x1 D x1 (7)
dt V V
Ahora definiendo D= F/V como velocidad de
dilución cuya unidad es [L/h], se obtiene al
sustituir en la ecuación anterior y agrupando
términos, la ecuación que modelará el
comportamiento de la concentración de
(m.o.) en el transcurso del tiempo:
dx1
x0 D G D D x1 (8)
dt
dx1
x0 D G D x1 (9)
dt
0 x0 D G D x1 (10)
D
x1 x0 (11)
G D
D
X 1 (s) X 0 (s) (12)
s D G
D
x1 xo (13)
D G
dx1
dt x0 D G D x1 Tanque # 1
dx2 x D D x Tanque # 2
dt 1 G 2
:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: (14)
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
dxn
dt xn 1 D G D xn Tanque # n
Aplicando la transformada de Laplace,
obtenemos el sistema de ecuaciones que modela
el comportamiento de la concentración de (m.o.)
en cada tanque en función de transferencia:
D
1 0 s D 0
X s x X s Tanque # 1
G
D
X 2 s x0 X1 s Tanque # 2
s G D
:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::: (15)
::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
D
X n s x0 s D X n 1 s Tanque # n
G
Se Puede simular el sistema de ecuaciones
(15) en Simulink mediante el diagrama de
bloque mostrado a continuación:
Fig 4. Diagrama de bloque para la simulación de
una batería de “n” bio reactores.
Clase N°
Desarrollo de un modelo
para una batería por 3 bio
reactores de CSTR de flujo
continuo
Para una batería conformada por tres bio
reactores de fermentación continua suponga
μG = 0,5 h-1 y D = 5 h-1, si x0 = 10 g/L, determine:
a) La respuesta dinámica (variación de la
concentración de los m.o.).
b) La respuesta dinámica para D = 6 h-1.
c) envié las variables al ambiente de Matlab y
presente sus graficas.
Para dar respuesta al inciso a) según la teoría, se
monta el diagrama de bloques con D = 5 h-1 y
μG = 0,5 h-1 mostrado en figura 7, donde el
bloque escalón simula la concentración de
(m.o.) en el efluente al tanque, en los bloques, la
función de transferencia, se han sustituidos los
valores de D y μG que poseen el mismo valor en
cada tanque y con los bloques visualizadores
denominados x0, x1, x2 y x3 se grafica la
concentración de (m.o.) de entrada y de salida
de cada tanque.
Fig. 7 Esquema de simulación de la batería de
tres bioreactores.
La figura 8 a, b, c muestra la respuesta dinámica
de la concentración de (m.o.) de salida de
cada tanque dada por los visualizadores x1, x2 y
x3, de estas gráficas se observa, como que el
valor de la respuesta a estado estacionario de
la salida de cada tanque es 1,1 el valor a
estado estacionario de la entrada a este y el
tiempo de establecimiento del sistema es
relativamente corto, de aproximadamente 2 h.
Para la respuesta del inciso b), en cada función
de transferencia cambiamos el valor de D a 0,6
h-1 en figura 7 y las respuesta correspondientes a
x1, x2 y x3, se grafican en figura 8 d, e, f, de esta
se observa que el valor de la concentración de
(m.o.) a estado estacionario es 6 veces el valor
de la respuesta a estado estacionario del
tanque anterior y que el tiempo de
establecimiento ahora es mucho mayor para el
ultimo tanque, alrededor de las 100 horas.
Luego, nuevamente se observa que hacer D
pequeño se traduce en un aumento de la
concentración de ( m.o.) de varias veces la
entrada, pero la respuesta dinámica del
sistema en este caso es mucho más lenta.
Por tanto si es conveniente en algún proceso
aumentar la concentración de (m.o.) debe
existir un compromiso con el tiempo en el
que se logra obtener el régimen estacionario.
Fig 8. Gráficos de los valores instantáneos de la
concentración de m.o
Fig 8. Gráficos de los valores instantáneos de la
concentración de m.o
Fig 8. Gráficos de los valores instantáneos de la
concentración de m.o
Fig 8. Gráficos de los valores instantáneos de la
concentración de m.o
Fig 8. Gráficos de los valores instantáneos de la
concentración de m.o
Fig 8. Gráficos de los valores instantáneos de la
concentración de m.o
Entre las ventajas de la simulación tenemos
que se nos permite visualizar:
♦ La respuesta del sistema ante un estímulo.
♦ La respuesta en cada bio reactor de un
batería de bio reactores en serie, donde la
salida de un tanque constituye el estímulo
del siguiente, lo que hace que la solución
analítica sea más compleja a medida que
el número de tanques aumenta.
La variación del comportamiento de la
respuesta (concentración de m.o.) a la
salida de un bio reactor como en cada
tanque de una batería de bio reactores en
función de la velocidad de dilución y de la
velocidad de crecimiento específica si
desea comparar diferentes microorganismos
El montaje del esquema de simulación en
Simulink es fácil, rápido y sistemático.
El aumento del número de tanques en una
batería no aumenta excesivamente la
complejidad del esquema de simulación,
se traduce en añadir más bloques función
de transferencia y visualizadores y esto no
crea una demora mucho mayor en el
tiempo de simulación.
.
.
.
.
.
.
.
