BALANCES DE MATERIA
PARA SISTEMAS DE CULTIVO
I. ECUACION DE BALANCE DE MATERIA EN EL BIORREACTOR
BIORREACTOR O FERMENTADOR
Es el equipo o recipiente donde se realiza el proceso fermentativo o bioprocesos en general.
Provee todos los servicios que son necesarios para el cultivo: mezclado, termostatizacin,
suministro de oxgeno, entradas para adicin nutrientes, control del pH, etc.
ESQUEMA DE UN BIORREACTOR
Se indica los caudales (F) y concentraciones del componente i (Ci) a la entrada (1) y salida (2).
La flecha que rodea al eje del agitador significa que el cultivo est perfectamente mezclado.
F1
Ci1
F2
Ci2
Figura 1
En la Figura 1 para un componente i cualquiera del cultivo (Ci): i = sustrato,
biomasa o producto, se puede plantear el siguiente balance de materia en el
biorreactor
Velocidad
de
acumulacin
(1)
d (V Ci)
dt
V
F1
F2
Ci1
Ci2
rf i
Velocidad
de
ingreso
Velocidad
de
salida
F1 Ci1
F2 Ci2
Velocidad
de
formacin
V rf i
Velocidad
de
consumo
V rci
Volumen de cultivo,
Caudal de alimentacin, (Lh-1)
Caudal de salida, (Lh-1)
Concentracin del componente "i" en la alimentacin
Concentracin del componente "i" en el caudal de salida
Si el cultivo esta bien mezclado, se puede asumir idntica a la que
hay dentro del biorreactor.
Velocidad de formacin del componente "i"
Velocidad de consumo del componente "i"
rci
EJERCICIO 1:
Aplicando la Ecuacin (1):
Deduzca la Ecuacin de Balance de Materia para las Clulas (Ci = X)
(1)
II.
d (V Ci)
dt
F1 Ci1
F2 Ci2
V rf i
V rci
d (V X)
dt
F1 X1
F2 X2
V rx
V rd
d (X)
dt
F1 X1 / V
F2 X2 / V +
V rx / V
V rd / V
d (X)
dt
F1 X1 / V
F2 X2 / V +
rx
rd
d (X)
dt
F1 X1 / V
F 2 X2 / V +
SISTEMAS DE CULTIVO Y ECUACIONES DE BALANCE DE
MATERIA
Los sistemas de cultivo o mtodos de cultivo, hace referencia al modo de
operar de un biorreactor, esto es en forma continua o discontinua.
Suponiendo que la densidad del cultivo (d = m / V) y de su alimentacin son
iguales resulta:
o
El volumen del cultivo (V) variar en el tiempo segn sea F1 y
(2)
dV
dt
F1 - F2
F1 Caudal de entrada
F2 Caudal de salida
Dependiendo de como sea F1 y F2 surgen tres sistemas de cultivo:
F2
1. Cultivo Continuo
F1
Ambos caudales son iguales
Por la ec. (2): V es constante
Por lo tanto: la ec. (1) se reduce a ec. (3)
F2
F1 = F2
(1)
d (V Ci)
dt
F1 Ci1
F2 Ci2
(3)
V d (Ci)
Dt
F ( Ci1
Ci2 )
V rf i
V ( rf i
V rci
rci )
2. Cultivo Batch Alimentado
El caudal de salida F2 , es nulo
F1
Por lo que V aumentar en el tiempo en
funcin del caudal de entrada:
Ec. (4)
F2 = 0
dV
dt
F1
En la Ec. (1) se anula el trmino:
(1)
d (V Ci)
dt
F1 Ci1
(5)
d (V Ci)
dt
F Ci1
F2 Ci2
F2 Ci2
V rf i
V ( rf i
V rci
rci )
3. Cultivo Batch
Ambos caudales son nulos F1 = F2 = 0
V es constante
En la Ec. (1) se anulan los trminos:
F1 Ci1
(1)
(5)
F2 Ci2
V rf i
V rci
V ( rf i
rci )
( rf i
rci )
d (V Ci)
dt
V d (Ci)
dt
d (Ci)
dt
F1 Ci1
F2 Ci2
Su duracin:
Tambin es limitada en el tiempo y depende de las condiciones iniciales del
cultivo
Una vez inoculado el medio:
La Biomasa aumenta a expendas de los nutrientes y cuando el substrato que
limita el crecimiento se agota, finaliza el batch.
II.
BALANCES DE MATERIA DE LOS COMPONENTES
( X, P, S ) PARA SISTEMAS DE CULTIVO
Analizaremos ahora con ms detalle cada uno de los sistemas vistos
aplicando en particular los balances de materia a la biomasa (X), al producto
(P), y al sustrato limitante de crecimiento (S).
Para la biomasa (X), y el producto (P) slo es necesario considerar la cintica
de formacin
rx
rp
respectivamente, con lo cual rci = 0 en ambos casos.
Para el sustrato (S) se tiene el caso inverso y slo deber considerarse la
velocidad de consumo, rs .
Si por las caractersticas del proceso, la lisis celular o la descomposicin del
producto son importantes, deber incluirse en el balance un trmino adicional
que contemple este aspecto.
II.1. Cultivo continuo
II.1.1. Operacin:
Para poner en marcha un cultivo continuo, se realiza previamente un
cultivo batch
En un momento dado se comienza a alimentar con medio fresco a
caudal F y por rebalse se mantiene el volumen constante.
3.1.2. Componentes del caudal de salida:
Clulas
Medio de cultivo parcialmente agotado y
Eventualmente algn producto.
3.1.3. Consideraciones:
Si alimentamos con medio fresco significa que:
X1 = 0
P1 = 0
Por lo que, slo se debe considerar la concentracin de sustrato
limitante del crecimiento (S1).
3.1.4. Balance de Materia para X, S, P
En base a las consideraciones anteriores y de acuerdo a Ec. (3) los
balances de materia para X, S y P ser la: Ec. (7) , Ec. (8) Ec. (9)
CUANDO: Ci = X
(7)
VdX
dt
DEDUCCION MATEMATICA DE EC. (7)
F X
V rX
(3)
CUANDO: Ci = S
(8)
VdS
dt
V d (Ci)
dt
F (Ci1 -Ci2) + V (rf i - rci )
Cuando: Ci = X
F (S1 S2) - V rS
VdX
dt
F ( 0 - X ) + V (rX - 0 )
CUANDO: Ci = P
(9)
VdP
dt
F P
V rP
(7)
VdX
dt
F X
V rX
3.1.5. En Estado Estacionario
Las concentraciones dentro del biorreactor permanecern constantes
en el tiempo, lo que significa igualar a cero las ecuaciones (7), (8) y (9).
3.1.6. Relacin entre la Velocidad de Dilucin ( D ) y Crecimiento Celular ( )
De la primera Ec. (7) y teniendo en cuenta que rx =
(10)
VELOCIDAD DE DILUCIN (D)
(10)
F
V
D =
X, resulta Ec.
DEDUCCION MATEMATICA EC. (10)
(7)
u
(10)
VdX
dt
F X
V rX
F X
V rX
F X
F
V
D
uX
Si se reemplaza D Ec. (10) en la ec. (19) resulta Ec. (11)
CONCENTRACION DE SUBSTRATO EN ESTADO ESTACIONARIO ()
(11)
= KS D / (u m - D)
DEDUCCI
(10)
(19)
Ejercicio de Aplicacin:
KS = 0.05 g L-1
D
= 0, 0.25, 0.5, 0.75, 1.0
h -1
K
S
D
u m = 1.0 h-1
1. Calcule los valores de
2. Traze la curva de evolucin de
Rpta: = 0.017 g/L; cuando D = 0.25
(11)
317 `Rlaci "ents a Vlcidad e Diui n>(D )
e Sbstra ( s
eha Ecc ,8)!$sr e c. (1
Consm
VEOIA E CON MO E
USRTO )
!
4
(12
) $ !
! xr " =
$ !
D !(S
DEDUCCION MATEMATICA Ec. (11)
(8)
VdS
dt
0
F (S1 S) - V rC
F (S1 S) - V rS
!!
ezcicio de AplicAcin:
E
D
= 0, 0.25, 0.5, 0.75, 1.0
h -1
(12)
V bS
F (S1 S)
rS
F (S1 S)
rS
D (S1 )
/ V
S1 = 5 g L-1
= valores del ejercicio anterior
1. Calcule rS para los valores de
Rpta:
Por la Ec. (26) :
VELOCIDAD DE CONSUMO DE
SUBSTRATO
(26)
rs = u X / Yx/s
DEDUCCION MATEMATICA Ec. (26)
(26)
Y x/s
rx / rs
rs =
rx / Y x/s
rs =
uX / Yx/s
Por lo tanto, comparando ec. (26) y dc. (12) y reemplazando u = D se tiene ec.(1#)
CONCENTRACION DE BIOMASA EN
DEDUCCION MATEMATICA Ec. (13)
ESTADO ESTACIONARIO (x)
(12)
rS
D (S1 )
(13)
Yx/s ( S1 - )
Ejebcicio de Aplicacin:
(26)
rS = u X / Yx/s
u X / Yx/s = D (S1 )
(26) y (12)
Reemphazando u = D
Yp-s = 0.5
D X = Yx/s D (S1 - )
S1 = 5 g L-1
= valores del ejercicio anterior
(13)
X = Yx/s (S1 - )
1. Calcule X
2. Gr!fi1ue la evolucin de la biomasa
Rpta:
X = 2.49 g/L;
cuando = 0. 17 g/L
Si en particular S es la fuente de carbono y energa, r S vendr dada por la Ec. (29)
del Cap. 5 y X ser:
CONCENTRACION DE BIOMASA EN
ESTADO ESTACIONARIO (x)
(13 B)
DEDUCCION MATEMATICA Ec. (13 B)
(11)
(13)
X = Yx/s ( S1 Ks D / u m D)
(13B)
= KS D / (u m - D)
=
Yx/s (S1 - )
Introduciendo (11) en (13)
resulta (13 B)
X = Yx/s ( S1 Ks D / u m D)
CONCENTRACION DE BIOMASA EN
ESTADO ESTACIONARIO (x)
(14)
D Yx/s ( S1 - )
(D + ms Yx/s )
DEDUCCION MATEMATICA Ec. (14)
En ec. (29) reemplazar rx = X
(29)
rs = rx / Yx/s + ms X
=
X / Yx/s
Factorizando en X
+ ms X
Ejercicio de Aplicacin:
D
= 0, 0.25, 0.5, 0.75, 1.0
rs = X [ / Yx/s + ms ]
h-1
Introducir la ec. (26):
rs = X /Yx/s
Yx/s = 0.5
X /Yx/s
S1 = 5 g L-1
= valores del ejercicio anterior
X [ / Yx/s + ms ]
Introducir la ec. (13):
ms = 0.06 g / g h
X =Yx/s (S1 - )
Yx/s (S1 - ) / Yx/s = X [ /Yx/s + ms ]
-------------------------------------------------1. Calcule X
2. Grafique la evolucin
Despejando X
X [ /Yx/s + ms ] .Yx/s = Yx/s (S1 - )
Rpta:
------------
Yx/s ( S1 - )
[ /Yx/s + ms ] .Yx/s
D
Yx/s ( S1 - )
D + ms . Yx/s
Cuando ms = 0 , la ecuacin (14) se reduce a la ecuacin (13).
Figura 16.
Cultivo continuo.
Muestra como vara la concentracin () en estado estacionario
en funcin de la velocidad de dilucin (D)
La curva superior continua corresponde a Ec. (13) y la inferior Ec. (14)
(13)
ms = 0
Yx/s ( S1 - )
(14)
ms = 0.08 g-1g-1h-1
D Yx/s ( S1 - )
(D + ms Yx/s)
S1 = 5 gL-1
Ks = 0.05 gL-1
Dc = 0.99 h-1
m = 1 h-1
Concentracin de biomasa (X) y de substrato limitante
(S1) en estado estacionario a distintos valores de D.
Parmetros:
m = 1 h-1
Ks = 0.05 gL-1 ;
Yx/s = 0.5
S1 = 5 gL-1
En el caso de la Ec.(14) el efecto del mantenimiento celular se hace notable a
bajas velocidades de dilucin.
En ambos, casos puede apreciarse que existe un valor D por encima del cual es
X = 0, con lo cual por la ec. (13) o (14) es = S1
(14)
D Yx/s ( S1 - )
(D + ms Yx/s)
(13)
Yx/s ( S1 - )
Yx/s ( S1 - )
D Yx/s ( S1 - )
(D + ms Yx/s)
= S1
D Yx/s ( S1 - )
= S1
Si se reemplaza este valor en la Ec. (11) se obtiene Ec. (15)
VELOCIDAD DE DILUCION CRITICA
Dc
DEDUCCION MATEMATICA Ec. (15)
(11)
= KS D / (u m D)
S1
(15)
Dc
S1 = KS D / (u m D)
m
Ks + S1
KS D = S1 . (u m D)
(15)
KS D =
u m S1
S1 D
KS D + S1 D =
u m S1
D (KS + S1 )
u m S1
S1
D
Ks + S1
Si como ocurre normalmente S1 >> Ks se tiene Dc m lo cual es un criterio
muy til en el momento de seleccionar un valor D apropiado, ya que debe
cumplirse que D < m
Caso contrario, ocurrir el lavado del cultivo debido a que la velocidad de
salida de las clulas del biorreactor ser mayor que la del crecimiento.
Debe tomarse en cuenta que la Fig. 16 se representa en caso que esta dado por
la ecuacin de Monod, pero si en el medio de cultivo hay sustancias inhibidoras
del crecimiento ( o son formadas por microorganismos ) o el substrato limitante
es el inhibidor, debern emplearse las expresiones de correspondientes a las
dadas en el capitulo 5.
3.1.8. Relacin entre la Velocidad Dilucin y Formacin de Producto
De la tercera Ec. (9) en estado estacionario se reduce a:
CONCENTRACION DE PRODUCTO EN
ESTADO ESTACIONARIO (P)
DEDUCCION MATEMATICA EC. (10)
(9)
(10)
O bien:
(17)
P =
rP / D
rP = qp . X
_
P = qp . X / D
VdP
dt
V rP
0 = -F P
F P = V rP
P = V rP / F
V rP
P =
F P
rP / D
Donde P : concentracin de producto en estado estacionario
Dependiendo de cmo sea la cintica de formacin de producto ser la forma de
la curva P vs D.
El cultivo continuo, alternativamente puede emplearse, y de hecho se emplea,
para elucidar que tipo de relacin existe entre q p y , ya que como se mencion
en el capitulo 6 es uno de los factores a tener en cuenta para planear la estrategia
de produccin, como se ver mas adelante cuando se trate la produccin de
penicilina.
DETERMINACION DE LOS PARAMETROS DE CRECIMIENTO
El cultivo continuo es sumamente til para determinar parmetros de crecimiento
tales como Ks, u m o Yx/s.
As reordenando la ec. (11) se tiene la ec.(15)
(11)
(15)
= KS D / (u m D)
KS D =
. (u m D)
KS D =
um
KS D +
D (KS +
D =
=
um
um
m
Ks +
La Ecuacin de Monod se puede linearizar tomando los recprocos para obtener
la expresin
m
D
m
+
Ks
m
D
Ks
1
Ks
1
+
Es una lnea recta
Graficando 1 / D en funcin de 1 / , los puntos debern ajustarse a una lnea
recta cuya interseccin con el eje 1/D
dar el valor de 1/ m y la pendiente
Ks / m , si es que el cultivo puede ser representado por una cintica como la de
Monod.
Y
a=Ks/m
1/D
b= 1/m
1/
Por otra parte reordenando la ec. (14) resulta ec. (19)
(14)
D Yx/s ( S1 - )
(D + ms Yx/s)
D Yx/s ( S1 - )
D ( S1 - )
(19)
X (D + ms Yx/s)
X (D + ms Yx/s) / Yx/s
D ( S1 - ) /
D / Yx/s
ms Yx/s / Yx/s
D ( S1 - ) /
D / Yx/s
ms
D ( S1 - ) /
1 / Yx/s . D
De donde la grfica de:
1 / Yx/s y ms
D ( S1 - ) / X
ms
Es una lnea recta
en funcin de D permitir estimar
Obviamente en esta caso el substrato considerado es la fuente de carbono y
energa
A = 1 / Yx/s
D ( S1 - ) / X
B = ms
3.2. Cultivo Batch Alimentado
3.2.1. Operacin:
Se inicia la alimentacin del cultivo cuando el substrato limitante se ha
agotado. No es un requisito indispensable, pero simplifica el tratamiento
matemtico.
3.2.2. Objetivo:
Controlar la
alimentacin
velocidad
de
crecimiento
mediante
la
velocidad
de
3.2.3. Consideraciones:
Supondremos que alimentamos con medio fresco, es decir:
X1 = 0
P1 = 0
Por lo que, slo deberemos considerar la concentracin de sustrato
limitante del crecimiento, S1 .
3.2.4. Balance de Materia para X , S , P
Segn Ec.(5)
(5)
d (V Ci)
dt
F Ci1
V ( rf i
rci )
Se tiene:
CUANDO: Ci = X
(20)
d XV
dt
DEDUCCION MATEMATICA DE EC.(20)
V rX
(5)
CUANDO: Ci = S
(21)
d SV
dt
d (VCi)
dt
F (Ci1 )
Cuando: Ci = X
F S1 V r S
dVX
dt
V (rf i - rci )
donde X = 0
F ( 0 ) + V (rX - 0 )
CUANDO: Ci = P
(22)
d PV
dt
V rP
(20)
d XV
dt
V rX
Si en la ec. (21) se reemplaza rS por rx / Yx/s que resulta de Yx/s = rx / rs y
reemplazando V rx por su equivalencia de la ec. (20).
DEDUCCION MATEMATICA DE EC.(23)
(23)
d SV
dt
F S1 1 / Yx/s . dXV/dt
(21)
d SV
dt
F S1 V rX / Yx/s
d SV
dt
F S1 V rS
F S1 1 / Yx/s . V rX
Pero como Ec. (20):
d XV
dt
(23)
d SV
dt
V rX
F S1 1 / Yx/s . dXV/dt
Si deseamos que la velocidad de crecimiento d(XV)/dt este controlada por la
alimentacin d(SV)/dt , esta deber ser tal que en todo momento S 0 y por lo
tanto d(SV) / dt 0
Esto equivale decir que el sustrato es consumido totalmente ni bien ingresa al
birreactor
Luego ec. (23) es:
(23)
1 / Yx/s
(24)
F S1 1 / Yx/s
dXV/dt
dXV/dt
dXV/dt
F S1
Yx/s . F S1
E integrando:
XV
dXV
XoVo
t
=
Yx/s . F S1
to
dt
(25)
XV
XV - Xo Vo
Yx/s . F S1 t
XV
Xo Vo + Yx/s . F S1 t
Xo Vo + Yx/s . F S1 t
Donde:
Xo = Concentracin de biomasa en el momento de iniciar la alimentacin.
Vo = Volumen del cultivo en el momento de iniciar la alimentacin
La variacin de V con el tiempo se obtiene integrando la ec. (4): dv/dt = F
V
dv
Vo
V - Vo
(26)
V - Vo
dt
to
F t
F t
UNIVERSIDAD NACIONAL DE TRUJILLO
ESCUELA POST GRADO
MAESTRIA: BIOTECNOLOGIA Y BIOINGENIERA
DOCENTE : MSc CARLOS E. VILLANUEVA AGUILAR
Universidad Nacional Pedro Ruiz Gallo Lambayeque 074-283610
villabal@[Link]
EJERCICIO N 1
Para un cultivo continuo se tiene los sgtes. parmetros:
m = 1 h-1 ; Yx/s = 0.5 ; Ks = 0.05 gL-1 ; S1 = 5 gL-1
a) Calcule la velocidad de dilucin crtica Dc
b) Si como ocurre normalmente S1 >> Ks . Calcule Dc.
c) Que criterio deber cumplirse en el momento de seleccionar un valor D apropiado En caso contrario que ocurrira?
EJERCICIO N 2
Aplicando la ecuacin general de balance de materia determine las ecuaciones de balance de un sistema de cultivo
batch o discontinuo para los componentes biomasa (X) , producto (P) y substrato (S).
EJERCICIO N 3
Tomando como base el ejemplo 4.6. Produccin de goma xantano, realice los clculos de balance similares a los
solicitados para cuando la cantidad de gas de salida recogido durante todo el proceso es de 1350 kg.
EJERCICIO N 4
Si considera el esquema de datos para el biorreactor indicado deduzca:
1.
2.
3.
4.
La ecuacin general de balance de materiales, tomando como patrn el La ecuacin general de balance
La ecuacin de balance cuando el componente es biomasa (X)
La ecuacin de balance cuando el componente es substrato (S)
La ecuacin de balance cuando el componente es producto (P)
EJERCICIO N 5
Deducida la ecuacin de balance de materia para biomasa (Ci = X) y suponiendo las siguientes condiciones para el
biorreactor:
Estado estacionario
Alimentacin estril
Velocidad de muerte despreciable
Establezca la ecuacin resultante para dichas condiciones