Balance de Masa y Energia Modelo Matematico
Balance de Masa y Energia Modelo Matematico
Balance de Masa y Energia Modelo Matematico
Modelo Semifsico de Base Fenomenolgica del Proceso Continuo de Fermentacin Alcohlica Ortega
Vol. 27(1), 21-32 (2016)
doi: 10.4067/S0718-07642016000100004
Recibido Jul. 7, 2015; Aceptado Jul. 27, 2015; Versin final Oct. 6, 2015, Publicado Feb. 2016
Resumen
Este trabajo presenta la construccin de un modelo semifsico de base fenomenolgica del proceso
continuo de fermentacin alcohlica con Saccharomyces cerevisiae. La estructura del modelo es obtenida
por los principios de conservacin de masa y energa aplicados al proceso. Los parmetros del modelo
fueron obtenidos por aplicacin de las ecuaciones de fenmenos de transporte y cintica de los procesos
microbianos. El modelo tiene en cuenta las inhibiciones por concentracin de biomasa, concentracin de
sustrato y concentracin de producto. Adicionalmente tiene en cuenta el efecto de la temperatura sobre
los parmetros asociados al metabolismo de la levadura. La simulacin y la validacin muestran que el
modelo obtenido puede predecir de forma satisfactoria el comportamiento dinmico del proceso
estudiado.
INTRODUCCIN
Los bioprocesos son el fundamento de muchas industrias de alimentos, qumicas y farmacuticas. Estos
hacen uso de clulas o partes de ellas, para manufacturar nuevos productos, destruir desechos peligrosos,
entre otras labores. El uso de microorganismos para transformar materiales biolgicos tiene su origen en la
antigedad y desde entonces los bioprocesos han sido desarrollados para enorme variedad de productos
comerciales, desde materiales relativamente econmicos tales como alcohol industrial, levaduras y
solventes orgnicos, hasta sustancias muy costosas como los antibiticos, enzimas y vacunas (Doran,
2013; Dutta, 2008). Debido a la enorme importancia de los bioprocesos en el mundo contemporneo, los
investigadores se plantean la necesidad de mejorarlos, para lo cual requieren de modelos matemticos que
describan con gran nivel de exactitud lo que est pasando en el proceso. Una vez definidos los modelos,
que pueden ser utilizados para describir el proceso de produccin bajo diferentes condiciones de proceso,
tales como temperatura, pH, aireacin y mezcla. Los modelos matemticos pueden facilitar la reduccin de
experimentacin innecesaria, la reduccin de los costes del proceso y aumentar la calidad del producto, y se
han utilizado en la comprensin y la mejora de los sistemas biolgicos, as como en la optimizacin y control
automtico de los bioprocesos (Dodi et al., 2012; Koutinas et al., 2012).
Con el aumento del inters en la aplicacin industrial de la fermentacin alcohlica, se han examinado
varios modelos cinticos para el crecimiento microbiano, la formacin de producto y consumo de sustrato,
encaminadas a mejorar la productividad en trminos de etanol, ya que este metabolito es muy codiciado en
la industria debido a su amplio uso comercial. Dentro de estas propuestas de modelos sobresalen los
Modelos Semifsicos de Base Fenomenolgica, los cuales han logrado describir los fenmenos que suceden
en la fermentacin (Dodi et al., 2012; Ribas et al., 2011; Gmez et al., 2008; Duarte, 1995; Martnez et al.,
1993). Lastimosamente, estas propuestas dejan a un lado algunas dinmicas que pueden ser
fundamentales en el diseo de sistemas de control y los modelos que representan la cintica microbiana no
incluyen efecto inhibitorio por concentracin de producto, por concentracin de sustrato y por efecto de la
temperatura. Por lo tanto, en la presente propuesta se incluyen las dinmicas de temperatura del fluido
utilizado para controlar la temperatura del fluido en el tanque y se adicionan los efectos inhibitorios que no
tienen en cuenta otras propuestas de modelado de la fermentacin alcohlica.
Este trabajo propone un Modelo Semifsico de Base Fenomenolgica (MSBF) del proceso continuo de
fermentacin alcohlica con Saccharomyces cerevisiae. Inicialmente se presenta el MSBF del proceso
continuo de fermentacin alcohlica obtenido aplicando los diez pasos de la metodologa de lvarez et al.
(2009). Luego se muestran los resultados de simulacin del modelo y se discute su validacin cualitativa. Se
termina con la seccin de conclusiones y las referencias.
Un MSBF se dice que es de base fenomenolgica porque toma su estructura de los balances de materia,
energa y cantidad de movimiento del proceso y es semifsico porque adiciona a la estructura ecuaciones
empricas para varios de sus parmetros llamadas ecuaciones constitutivas (lvarez et al., 2009). El
planteamiento de MSBF en los bioprocesos debe su importancia a que an existen fenmenos que no estn
totalmente explicados, por ejemplo el estrs de una clula con respecto a una variacin de su medio
ambiente, ya sea por pH, temperatura o gradiente de velocidad, recurrindose al empirismo para obtener
ecuaciones matemticas que permitan predecir estos cambios. A continuacin, se plantea un Modelo
Semifsico de Base Fenomenolgica del proceso continuo de fermentacin alcohlica aplicando los diez
pasos de la metodologa de lvarez et al. (2009).
El modelo contestar la pregunta: Cmo cambia la produccin de etanol a la salida del fermentador, ante
cambio en las variables externas como el caudal, temperatura y concentracin de glucosa en la corriente de
entrada del fermentador?
Objetivo del modelo. Predecir dinmicas de volumen, concentracin de glucosa, concentracin de biomasa,
concentracin de etanol y temperaturas en el fermentador y en la chaqueta. Con esta informacin, el modelo
podr usarse en el control, en el escalado y optimizacin del proceso de fermentacin.
Hiptesis del modelo. El proceso de fermentacin alcohlica se realiza por la metabolizacin de la glucosa
por la levadura Saccharomyces cerevisiae para producir etanol como metabolito principal. La velocidad de
crecimiento de la levadura se representa por un modelo cintico con efecto inhibitorio de concentracin de
sustrato, concentracin celular y de concentracin de producto. La velocidad de formacin de etanol se
asocia al crecimiento celular y puede ser expresada por la ecuacin de Luedeking-Piret (Clarke, 2013;
Ccopa et al., 2006). Durante el metabolismo de la levadura se genera calor, el cual es retirado por un fluido
de enfriamiento por el mecanismo de la conveccin de calor.
Supuestos del modelo: El jarabe posee slo glucosa como fuente de carbono; se garantiza mezcla perfecta
dentro del tanque; el CO2 liberado desplaza el O2 residual al inicio de la fermentacin; las clulas y el caldo
pueden considerarse una sola fase lquida; el etanol es despreciable en la corriente gaseosa de salida; los
parmetros asociados con propiedades de los fluidos (densidad, viscosidad absoluta y calor especfico) son
constantes; la masa del gas-vapor retenida es muy pequea comparada con el lquido; el equilibrio trmico
entre las fases es rpido; el pH es constante durante el proceso con un valor de 4 (pHptimo); debido a que el
tamao de una levadura es menor de 10 micrmetros, no se considera muerte por estrs hidrodinmico
(Trujillo y Valdez, 2006).
Se definen tantos Sistemas de Proceso (SdeP) sobre el proceso que se modelar como lo exija el nivel de
detalle y representar la relacin de todos los SdeP en un diagrama de bloques.Como el modelo ser de
carcter macroscpico, se pueden tomar dos sistemas de proceso: 1. Fluido contenido en la chaqueta y 2.
Lquido dentro del reactor. Como se observa en la Fig. 2, la interaccin entre los dos SdeP es energtica y
se da a travs de la transferencia de calor entre el lquido en el fermentador y el agua de enfriamiento
contenida en la chaqueta.
Paso 4. Aplicar el principio de Conservacin sobre cada uno de los sistemas de proceso.
Sistema de Proceso I:
dVw
w =w V -w V wout (1)
dt win
Como la densidad es considerada constante y no hay acumulacin de masa en la chaqueta, la ecuacin (1)
se reduce a una ecuacin algebraica:
dTw 1
= [ c V (T -T )+Q intercambio ] (3)
dt w cpw Vw w pw win win w
dV
=V in -V out (4)
dt
dX 1 dV
= (X V -XV out +Vrx -X ) (5)
dt V in in dt
dS 1 dV
= (Sin V in -SV out -Vrs -S ) (6)
dt V dt
dP 1 dV
= (-PV out +Vrp -P ) (7)
dt V dt
dT 1 dV
= [cp V in (Tin -Tref )-cp V out (T-Tref )-cp (T-Tref ) +Q gen -Q intercambio +Q agitador ] (8)
dt cp V dt
Paso 5. Seleccionar de las EDB aquellas con informacin valiosa para cumplir con el objetivo del modelo.
Todas las EDB desde (3) a la (8) son importantes y son las que brindan toda la informacin sobre la
dinmica de inters del proceso.
Paso 6. Definir para las EDB esenciales, los parmetros, variables y constantes conocidas en cada SdeP.
(i) Parmetros. Los parmetros que requieren ecuaciones constitutivas para su clculo estn ligados a la
cintica del proceso (rx, rs, rp, Q gen ), al flujo de calor intercambiado con la chaqueta (Q intercambio ) y a la
potencia suministrada por el agitador que se disipa en energa calrica (Q agitador ); (ii) Estados. S, X, P, T, Tw,
V; (iii) Variables de entrada (acciones de control o perturbaciones).Sin, Xin, Tin, Twin, V win, V in; y (iv)
Constantes. La densidad y el calor especfico de los fluidos en el biorreactor y en la chaqueta, que se
asumen que varan poco con los cambios de concentracin y temperaturas de operacin. El volumen de la
chaqueta.
Se refiere a encontrar las ecuaciones que permitan calcular el mayor nmero de parmetros en cada
Sistema de Proceso. Para la velocidad de crecimiento celular: Se utiliza una cintica que contiene efectos
inhibitorios de concentracin de sustrato, concentracin de producto y concentracin de biomasa (Clarke;
2013; Akpa, 2012; Alt y Markov; 2012; Dodi et al., 2012; Tian et al., 2011; Ccopa et al., 2006; Schgerl y
Bellgardt, 2000):
S -k S X m P n
rx =max e i (1- ) (1- ) X-Kd X (9)
ks +S Xmax Pmax
Para la velocidad de produccin de etanol: Se considera que la produccin de etanol se debe al crecimiento
junto con el mantenimiento del microorganismo, este fenmeno se describe por la ecuacin de Luedeking-
Piret (Liu, 2013; Ccopa et al., 2006; Schgerl y Bellgardt, 2000):
rx
rs = +mx X (11)
Yx
De acuerdo a la descripcin de las ecuaciones (9)-(11), hay 12 parmetros que tienen que ser estimados,
dentro de los cuales hay cinco parmetros dependientes de la temperatura, max, Xmax, Pmax, Yx y Ypx , y esta
dependencia puede ser descrita por la ecuacin (12):
B D
Parmetro(T)= Ae T +Ce T (12)
Parmetro A B C D
Para la potencia de agitacin: Para el clculo de la potencia de agitacin se utiliza la ecuacin (13), la cual
est definida para valores de Re (D2i N/) mayores a 10000. Esta potencia se considera que se transforma
Para la transferencia de calor: La velocidad de transferencia de calor entre el fluido del biorreactor y el fluido
de la chaqueta es:
Para el clculo del Coeficiente Global de Transferencia de Calor, se utiliza la expresin de circuito
equivalente de resistencias entre la chaqueta y el fluido en el interior del biorreactor:
1
U= 1 e 1 (15)
+ +
hi kmetal he
Para calcular los valores de los coeficientes convectivos de la ecuacin (15) se utilizan las ecuaciones
descritas en (16) y (18). Estas ecuaciones dependen de relaciones entre parmetros de diseo y de la
geometra del biorreactor, las cuales estn debidamente reportadas en Mohan et al. (1992) y Katoh y
Yoshida, (2009), por tanto no sern mencionadas en este trabajo, al lector interesado se invita a
consultarlas.
Coeficiente convectivo de calor externo (chaqueta) (Mohan et al., 1992): La chaqueta se puede considerar
como un nulo que tiene un dimetro equivalente (De ).
De vw
Para valores de Re ( ) mayores a 10000, el coeficiente convectivo se puede estimar por la ecuacin
w
(16):
0.68 0.14
he De De vw cpw w 1/3 w
=0.81 ( ) ( ) ( ) (16)
kw w kw wwall
Para valores de Re menores a 2100, el coeficiente convectivo se puede estimar por la ecuacin (17):
1/3 0.14
he De De vw cpw w 1/3 De
=1.86 ( ) ( ) ( )( w ) (17)
kw w kw L wwall
Donde:
D22 -D21
De = (18)
D1
Coeficiente convectivo de calor interno (superficie lateral interna del tanque), (Katoh and Yoshida, 2009):
Di 2 N
Para 40<Re= < 300000, el coeficiente convectivo se puede estimar por la ecuacin (19):
2
1 0.14
hi DT Di 2 N 3
cp 3 (19)
=0.54 ( ) ( ) ( )
k k wall
Para el flujo de salida (Smith y Corripio, 1997): El flujo de salida se puede estimar por la expresin que
representa el flujo a travs de una vlvula:
GL = #Ecuaciones - # Incgnitas. El MSBF final est constituido por las ecuaciones diferenciales (3)-(8), y
las ecuaciones constitutivas (9)-(11), (13), (14), (20) y (21). El modelo contiene 13 incgnitas: los estados X,
S, P, T, Tw y V; los flujos de calor representados en Q agitador , Q gen , Q intercambio ; las velocidades , y y el
).Con lo cual se puede decir que el modelo es coherente y con GL=0.
flujo de salida del biorreactor (
La resolucin de las ecuaciones diferenciales y algebraicas simultneas se hizo con Matlab 2014a, usando
el mtodo de resolucin de Runge-Kutta con paso de 0.01h. Las diferentes condiciones para la simulacin
se determinaron buscando que ayudaran en la validacin del modelo, tal y como se discute en la siguiente
seccin.
Se realiza la simulacin del modelo definiendo algunas caractersticas de diseo: tipo de tanque, volumen y
dimetro del tanque, cantidad y tipo de deflectores, tipo, tamao y velocidad del agitador. Se definen
adems el tipo de enfriamiento y otras caractersticas.
Para realizar la simulacin del modelo, se tomaron las siguientes caractersticas de diseo: Tanque
3
cilndrico de lote cuadrado con volumen de lquido en estado estacionario de 1m , dimetro del tanque de
1.08 m, con 4 deflectores de 0.09 m de ancho, un agitador tipo turbina de 0.3613 m de dimetro girando a
3
60 rpm. La chaqueta del tanque tiene un volumen de 0.3329 m donde fluye agua de enfriamiento a 25 C y
5/2
un dimetro equivalente de 0.3436 m. La vlvula de salida tiene un de 0.1921 m /s. El resto de los
parmetros utilizados en las simulaciones se presentan en la Tabla 2 (Ribas et al., 2011; McCabe et al.,
2007; Ccopa et al., 2006).
Parmetro ks ki Kd mp mx m n HR Np
Valor 4.1 0.002 0.0018 0.1 0.2 1 1.5 557330 6
Con el modelo totalmente identificado, se realizaron varias simulaciones en Matlab. La Fig. 3 muestra el
resultado para una fermentacin que se realiza en 72 horas. Se inicia la simulacin con las siguientes
3
condiciones iniciales: Tasa de dilucin (D) de 0.2 /h, un valor de biomasa X= 5 kg/m , un valor de sustrato
3 3
S=150 kg/m . La concentracin del sustrato a la entrada es Sin= 150 kg/m con una temperatura de entrada
T=30C. El comportamiento mostrado por el modelo es tpico de los procesos continuos de fermentacin
reales, en el cual se tiene un sustrato que se consume para obtener un producto con crecimiento de
biomasa y existen dinmicas rpidas de temperaturas de los fluidos dentro del tanque y en la chaqueta
como consecuencia del calor producido por la fermentacin y retirado por el fluido de servicio para
enfriamiento (Ccopa et al., 2006; Schgerl y Bellgardt, 2000).
3
Si se realiza una perturbacin tipo escaln en la concentracin de sustrato a un valor de 300 Kg/m a las 10
h de iniciado el proceso, se obtiene el resultado mostrado en la Fig. 4. Se observa que la concentracin de
biomasa disminuye al igual que la concentracin de producto, lo que implica que una cantidad apreciable
del sustrato no se consume, este fenmeno se presenta por el efecto de inhibicin por sustrato, el cual es
descrito en la ecuacin (9).
CONCLUSIONES
El modelo contempla acercamientos al proceso real como son las inhibiciones por biomasa, sustrato y
producto, y efectos de la temperatura sobre los parmetros asociados al metabolismo de la levadura.
Se logr el objetivo del modelo, como lo demostraron la simulacin y la validacin cualitativa: predecir las
dinmicas de volumen, concentracin de sustrato (glucosa), concentracin de biomasa (Saccharomyces
cerevisiae), concentracin de producto (etanol) y temperatura en el fermentador y en la chaqueta.
NOTACIN
2 3
A: rea de transferencia de calor [m ] rs : Velocidad de consumo de sustrato [Kg/(m s)]
3
cp : Calor especfico del caldo de rx : Velocidad de formacin de biomasa [Kg/(m s)]
fermentacin [J/(Kg K)]
cpw : Calor especfico del fluido en la chaqueta S: Concentracin de sustrato en el biorreactor
3
[J/(Kg K)] [Kg/m ]
C'v: Coeficiente de la vlvula en la salida Sin: Concentracin de sustrato a la entrada al
5/2 3
[m /s] biorreactor [Kg/m ]
De : Dimetro equivalente de la chaqueta [m] T: Temperatura en el caldo de fermentacin [K]
Di : Dimetro del agitador [m] Tin : Temperatura en la entrada al biorreactor [K]
DT : Dimetro del biorreactor [m] Tout : Temperatura a la salida del biorreactor [K]
D1 : Dimetro interno de la chaqueta [m] Tref : Temperatura de referencia [K]
2 : Dimetro externo de la chaqueta [m] Twin : Temperatura en la entrada a la chaqueta [K]
e: Espesor metlico del biorreactor [m] Twout : Temperatura a la salida de la chaqueta [K]
REFERENCIAS
Akpa, J., Modeling of a bioreactor for the fermentation of palmwine by Saccaharomyce cerevisiae (yeast)
and lactobacillus (bacteria). American Journal of Scientific and Industrial Research: 3(4), 231-240 (2012)
Alt, R. y Markov, S., Theoretical and computational studies of some bioreactor models. Computers and
Mathematics with Applications: 64, 350-360 (2012)
lvarez, H. y otros tres autores, Metodologa para la Obtencin de Modelos Semifsicos de Base
Fenomenolgica Aplicada a una Sulfitadora de Jugo de Caa de Azcar. Revista Iberoamericana de
Automtica e Informtica Industrial: 6(3), 10-20 (2009)
Ccopa, E. y otros cinco autores, Evaluation of optimization techniques for parameter estimation: Application
to etanol fermentation considering the effect of temperature. Process Biochemistry: 41, 1682-1687 (2006)
Clarke, K., Bioprocess engineering. An introductory engineering and life science approach. Wood head
Publishing, Cambridge, UK (2013)
Dodi, J. y otros cinco autores, Kinetic modelling of batch ethanol production from sugar beet raw juice.
Applied Energy: 99,192-197 (2012)
Doran, P., Bioprocess Engineering Principles. Second edition. Academic Press. Oxford, Reino Unido (2013)
Duarte, A., Introduccin a la ingeniera bioqumica. Universidad Nacional, Facultad de Ingenieras. Bogot,
Colombia (1995)
Dutta, R.,Fundamentals of Biochemical Engineering. Springer-Ane Books India. Nueva Delhi, India (2008)
Gmez, C., Caldern, Y. y lvarez, H., Construccin de modelos semifsicos de base fenomenolgica. Caso
proceso de fermentacin. Revista Facultad de Ciencias Agropecuarias: 6(2), 28-39 (2008)
Katoh, S. y Yoshida, F., Biochemical Engineering. A Textbook for Engineers, Chemists and Biologists.
WILEY-VCH. Weinheim, Alemania (2009)
Koutinas, M. y otros tres autores, Bioprocess systems engineering: transferring traditional process
engineering principles to industrial biotechnology. Computational and structural Biotechnology Journal: 3(4),
1-9 (2012)
Liu, S., Bioprocess engineering. Kinetics, biosystems, sustainability, and reactor design. Elsevier.
Amsterdam, Holanda (2013)
Martnez, L. y otros tres autores, Fermentacin etanlica continua anaerbica con Saccharomyces
cerevisiae resistentes al alcohol. Anales de qumica: 89 (5-6), 579-583 (1993)
McCabe, W., Smith, J. y Harriot, P., Operaciones Unitarias en Ingeniera Qumica. McGraw-Hill
Interamericana. Mxico D.F., Mxico (2007)
Mohan, P., Emery, A. y Al-Hassan, T., Heat transfer to Newtonian Fluids in mechanically agitated vessels.
Experimental Thermal and Fluid Science: 5(6), 861-883 (1992)
Ribas, M. y otros cuatro autores, Metodologa para la modelacin matemtica de procesos. Caso de estudio,
fermentacin alcohlica. ICIDCA. Sobre los derivados de la caa de azcar: 45(1), 37-47 (2011)
Schgerl, K. y Bellgardt, K., Bioreaction Engineering. Modeling and Control. Springer. Berln, Alemania.
(2000)
Tian, Y.y otros tres autores, Theoretical approach to modelling and analysis of the bioprocess with product
inhibition and impulse effect. BioSystems: 104, 77-86 (2011)
Trujillo, M y Valdez, N., El estrs hidrodinmico: Muerte y dao celular en cultivos agitados. Revista
Latinoamericana de Microbiologa: 48 (3-4), 269-280 (2006)