Vol. 4 (1) pp.

15-23, enero-junio 2014

RESEARCH ARTICLE

Anlisis de estabilidad de un sistema de fermentacin

acetona-butanol-etanol (ABE) a partir de glucosa em-
pleando Clostridium acetobutylicum ATCC 824
Laura Johana Castellanos Surez,*1 Luis Gernimo Matallana Prez*2
y Luis Javier Lpez Giraldo3
1
Universidad Industrial de Santander, Facultad de Ingenieras Fisicoqumicas, Escuela de Ingeniera Qumica, Carrera 27 calle 9 Ciudad Universitaria, Bucaraman-
ga, Colombia. *Autor para correspondencia: [email protected]

2
Universidad de Antioquia, Facultad de Ingeniera, Calle 67 N 53-108, Medelln, Colombia *Autor para correspondencia: [email protected]

3
Universidad Industrial de Santander, Facultad de Ingenieras Fisicoqumicas, Escuela de Ingeniera Qumica, Carrera 27 calle 9 Ciudad Universitaria, Bucara-
manga, Colombia.

Stability analysis of acetone-butanol-ethanol (ABE) fermentation by Clostridium aceto-

butylicum ATCC 824 with glucose feeding

Abstract
In the present study we analyzed the dynamics and stability in production of butanol by fermentation on a con-
tinuous bioreactor. The dynamic response of the system is studied as an evaluation criterion selecting the con-
centration of butanol obtained and using a kinetic model statistically validated, which describe the continuous
cultures oscillatory behavior taking into account the influence of pH, inhibition by products and influencing para-
meters in the system. The bifurcation parameters chosen were the dilution rate (D) and feed substrate concentra-
tion (Sf). We used a calculation routine developed in Matlab and MATCONT to generate diagrams for dynamic
behavior, steady state values and bifurcation diagrams. The results show that it is convenient to operate at sma-
ller dilution rates (D <0.1 [h ^ (-1)]) and feed substrate concentrations (Sf) in the range of 40 to 90 [g / L] for higher
concentrations of butanol. The practical implications of the phenomena discovered are analyzed in order to apply
an appropriate methodology for this type of study and the implementation of it in bioprocesses.

Keywords: Bifurcation, Butanol, Dynamic simulation, Nonlinear dynamics, Stability Analysis.

Editora: Pataquiva-Mateus, A. Y.
Resumen
Citation: Castellanos, L. J., Matallana, L. G. y Lpez, L. J. (2014). Anlisis de Se analiza la dinmica y estabilidad en la produccin
estabilidad de un sistema de fermentacin Acetona-Butanol-Etanol (ABE) a
partir de glucosa empleando Clostridium Acetobutylicum ATCC 824. Revista
de butanol por fermentacin en un biorreactor conti-
Mutis 4(1); pag. 15-23 nuo. La respuesta dinmica del sistema se estudi eli-
giendo como criterio de evaluacin la concentracin
Received: Jun 30, 2014; accepted: Junio 8, 2014; Published on line: Ju-
nio 30, 2014
de butanol, en un modelo cintico validado estadsti-
camente, que describe el comportamiento oscilatorio
Copyright: 2014 Castellanos et al. This is an open-access article, which de cultivos continuos teniendo en cuenta la influen-
permits unrestricted use, distributions and reproduction in any medium, cia del pH, inhibicin por productos y los parmetros
provided the original author and source are credited
influyentes en el sistema. Se tomaron como parme-
Competing Interests: The authors have no conflict of interest.
MUTIS, Journal of the Faculty of Sciences and Engineering, Jorge Tadeo Lozano University, is licensed under the Creative Commons 4.0: Attribution - Noncommercial -
No Derivative Works
 Fermentacin ABE empleando Clostridium A.
16

tros de bifurcacin la dilucin (D) y la concentracin anterior debido a que la relacin estequiomtrica ai-
de sustrato de alimentacin (Sf). Se emple una ruti- re:combustible es menor y en consecuencia es posi-
na de clculo en Matlab y la herramienta Matcont ble inyectar mayor cantidad de combustible (Patakova
para generar diagramas de comportamiento dinmi- et al., 2011). Adems, estas mezclas tienen mejores
co, valores de estado estacionario y diagramas de bi- propiedades fisicoqumicas que las de etanol en cuan-
furcacin. Los resultados muestran que es convenien- to a la estabilidad de fase en presencia de agua, me-
te operar a velocidades de dilucin pequeas (D<0,1 jores propiedades a baja temperatura, menor suscep-
) y concentraciones de sustrato de alimentacin tibilidad a la oxidacin durante el almacenamiento y
(Sf) en un rango de 40 a 90 [g/L] para obtener las con- una menor presin de vapor.
centraciones ms altas de butanol. Las implicaciones
prcticas de los fenmenos descubiertos se analizan Se infiere entonces que el butanol es una buena op-
con el objetivo de aplicar una metodologa adecuada cin para diversificar la oferta energtica de combus-
para este tipo de estudios y para la puesta en marcha tibles lquidos y tambin que es el producto ms sig-
de este tipo de bioprocesos. nificativo de la fermentacin Acetona-Butanol-Etanol
(ABE). Diversos estudios han demostrado que su pro-
Palabras clave: anlisis estabilidad, bifurcacin, buta- duccin por va fermentativa se puede desarrollar em-
nol, dinmica no lineal, simulacin. pleando diferentes tipos de bacterias, pero es de in-
ters para el presente trabajo, citar el gnero Clostri-
dium especie acetobutylicum debido a que es el ms
Introduccin empleado industrialmente: presenta mayor potencial
La produccin de alcoholes empleando fermentacin para la produccin eficiente del butanol (Jang et al.,
de sustratos azucarados tiene un amplio mercado, 2012; Kamiski et al., 2011; Patakova et al., 2013; Tas-
que abarca desde la industria farmacutica hasta su hiro et al., 2010).
empleo como aditivos oxigenantes de combustibles
Por otra parte, el incremento en la demanda de al-
para reducir su nivel contaminante. La fermentacin
coholes obtenidos por va fermentativa para ser em-
ABE (Acetona-Butanol-Etanol) es un ejemplo tpico y
pleados como biocombustibles, ha generado en las
fue uno de los primeros procesos fermentativos de-
industrias la necesidad de aumentar la productividad
sarrollados industrialmente. No obstante, el aumen-
de estos procesos. Adems, debido a los altos precios
to de los costos de los sustratos y la posibilidad de
de las materias primas, se hace necesario el diseo de
obtener materias primas y productos ms econmi-
procesos eficientes con el fin de minimizar los costos
cos desde la industria petroqumica, fueron las causas
de operacin y maximizar los rendimientos de los pro-
para que este proceso pasara a un segundo plano. En
ductos de inters. Por lo tanto, se hace imprescindi-
los ltimos aos, el uso de biocombustibles ha sido
ble plantear estudios y anlisis para el entendimiento
ampliamente desarrollado y su produccin ha crecido
detallado y optimizacin de los procesos fermentati-
exponencialmente (Donaldson et al., 2007; Kirschner,
vos para la obtencin de alcoholes. Por ejemplo, Paz
2006; Lee et al., 2008). Este renovado inters se basa
(Paz, 2011a; Paz, 2011b) propone que en el diseo de
casi exclusivamente en la preocupacin por el agota-
procesos en general, se debe incluir un anlisis rigu-
miento de los combustibles fsiles y la incertidumbre
roso de la estabilidad de los posibles estados estacio-
de los suministros de petrleo.
narios del sistema. De esta manera, ser viable prede-
As las cosas, es necesario la diversificacin de las cir y con mayor certeza el comportamiento dinmico
fuentes de energa al tiempo que se mejoren las ca- del mismo, y as establecer condiciones de operacin
ractersticas de los biocombustibles existentes o se adecuadas para el desempeo eficiente del proceso.
propongan nuevos con mejores caractersticas fisico-
qumicas. En ese sentido, el butanol presenta ventajas
con respecto al etanol y mejora considerablemente Mtodo de un anlisis de estabilidad para
las propiedades de los combustibles fsiles utilizados. un bioproceso
Por ejemplo, las mezclas butanol/gasolina hasta una Contrario a lo esperado, el diseo de los bioprocesos
relacin del 30% (v:v) no requieren modificaciones en a escala industrial an es desarrollado empricamente
el motor y adems se mejora la potencia de este. Lo a pesar de los estudios presentados por la comunidad

Revista Digital de la Facultad de Ciencias Naturales e Ingeniera de la UJTL

Castellanos et al., 2014
17

cientfica, ignorando los lmites de estabilidad caractersti-

cos de los sistemas. Por esta razn, el diseo de estos pro- no depende explcitamente de la variable independien-
cesos ha sido clasificado como un arte emprico: preocu- te . As, se dice que el sistema descrito anteriormente
pado principalmente por saber hacer sin tener en cuen- en la Ecuacin 1, es autnomo o invariante en el tiempo
ta la importancia de conocimientos cientficos que expli- (Ajbar et al., 2012; Strogatz, 1994; Seydel, 2010). En re-
quen el por qu de estos (Ajbar et al., 2012; Stockar et al., lacin con el nmero de ecuaciones , corresponden
2003). Lo anterior, ha llevado a elevar los costos globales a la cantidad de variables de estado y tambin denotan
de los bioprocesos, y en la mayora de los casos, ha sido un la dimensin del sistema.
factor diferenciador y de gran desventaja frente a los pro- La seleccin del modelo se realiz con base en el anlisis
cesos tradicionales. Es entonces evidente que la naturale- de las referencias bibliogrficas. As las cosas, el modelo
za dinmica del proceso fermentativo debe ser estudiada de Yang X. (1990) y Tsao G. (Yang et al., 1994) fue selec-
con el propsito de entender el efecto de las variables de cionado porque: (i) describe el comportamiento de inhi-
diseo sobre los fenmenos de estabilidad que afectan la bicin por producto, (ii) plantea una expresin cintica
productividad. para el crecimiento celular, (iii) considera el efecto del pH
Por consiguiente, el anlisis de un bioproceso, desde un sobre la cintica de crecimiento celular, (iv) se tiene cono-
punto de vista de estabilidad necesitar mnimo: cimiento de cada una de las constantes cinticas, (v) est
validado experimental y estadsticamente y, (vi) conside-
1. Modelo cintico que describa el proceso fermentativo. ra la recuperacin in situ de solventes (modelo que consi-
dera la adicin de resinas).
2. Parmetros a evaluar en el anlisis de estabilidad.

3. Conceptos necesarios para realizar un anlisis de esta- Parmetros a evaluar en el anlisis

bilidad del sistema.
de estabilidad
La velocidad de dilucin (D), y la concentracin de sus-
Modelo cintico trato alimentado (Sf) son los parmetros elegidos para
Durante las ltimas dos dcadas se han logrado signi- el anlisis de bifurcacin debido a que estos son los ms
ficativos avances en el rea de investigacin relaciona- fciles de manipular durante la operacin de un sistema
da con el modelamiento matemtico de la fermentacin de fermentacin, y por lo tanto es realista pensar que
ABE (Mayank et al., 2012). Los modelos de rutas meta- los resultados puedan ser empleados en la puesta en
blicas pueden clasificarse en 2 categoras: los mode- marcha de este (Garhyana et al., 2003; Mahecha-Botero
los estequiomtricos y modelos de simulacin cinticos et al., 2005). En la Figura 1 se presenta un diagrama es-
(Tashiro et al., 2010); Los ltimos describen el compor- quemtico de los biorreactores junto con sus variables
tamiento dinmico de los metabolitos y permiten desa- de entrada y salida.
rrollar estrategias para la operacin de biorreactores.

En el anlisis de estabilidad de biorreactores continuos,

los sistemas son modelados por medio de conjuntos
de ecuaciones diferenciales ordinarias de primer orden
(EDO). Las ecuaciones diferenciales usadas para modelar
los sistemas, expresadas de forma vectorial, son:

(1)

Donde t representa el tiempo, es el

vector de variables de estado, es el
vector de las funciones no lineales,
es el conjunto de parmetros del modelo, y es el vec- Figura 1. Diagrama esquemtico del fermentador continuo
tor de las condiciones iniciales. Tambin se asume que

Vol. 4 N. 1 pp. 15-23, enero-junio 2014

 Fermentacin ABE empleando Clostridium A.
18

Anlisis de estabilidad
(3)
La simulacin de biorreactores continuos a travs de
modelos cinticos genera el reto de determinar cules (4)
de los mltiples estados estables posibles son los ms
adecuados para la operacin del proceso (Fradkov et (5)
al., 1999). Como se ha mencionado en prrafos an-
teriores, este anlisis debe ser desarrollado basado (6)
en criterios ingenieriles y no en criterios empricos. A
continuacin se presentan brevemente los conceptos (7)
fundamentales y el orden metodolgico asociados a
(8)
un anlisis de estabilidad.
(9)
Estados estables o estacionarios
Los estados estacionarios de un sistema dinmico, re-
presentan las soluciones invariantes en el tiempo del
sistema dinmico representado por (1) y deben cum-
plir que: Donde

(2) ,

Siendo (2) un sistema de ecuaciones algebraico, gene-

ralmente de carcter no lineal. Estos puntos de estado
estacionario, t les llama puntos fijos, puntos de equi-
librio, o puntos crticos (Ajbar et al., 2012; Strogatz,
1994). Las soluciones de estado estacionario sa-
tisfacen todo el conjunto de ecuaciones algebraicas,
donde los parmetros de bifurcacin tambin tienen
valores fijos determinados.

En los procesos fermentativos, los modelos se carac-

terizan por ser altamente no lineales y en consecuen-
cia conducirn a soluciones (estados estacionarios)
que pueden ser o no estables. Esta situacin genera la
necesidad de analizar cules de las soluciones condu-
cen a regiones de estabilidad (Paz et al., 2011a; Seydel
2010; Strogatz, 1994). Resolviendo todas las ecuacio-
nes de carcter no lineal (2), se obtienen las concen-
traciones que satisfacen la condicin, siendo estas los
valores de estado estacionario de los modelos.

Para obtener los valores de estado estacionario se

emplea el software MATLAB. Se desarroll una ru-
tina de clculo que permite la solucin simultnea de
las ecuaciones diferenciales para el modelo cintico
elegido (3-10), solucionando los balances en estado
no estacionario y al mismo tiempo, permite obtener
las grficas del comportamiento dinmico.

Revista Digital de la Facultad de Ciencias Naturales e Ingeniera de la UJTL

Castellanos et al., 2014
19

Estabilidad local de estados estacionarios cin de . El objetivo de la teora de la bifurcacin es

Para la operacin de un proceso se hace necesario te- caracterizar cualitativamente los cambios en el compor-
ner una medida cuantitativa de su estabilidad. Este tipo tamiento dinmico de los sistemas no lineales mientras
de informacin puede ser obtenida mediante la linea- los parmetros principales o elegidos cambian.
lizacin alrededor de un punto fijo, correspondiente a
los estados o valores estacionarios del proceso (Seydel Los diagramas de bifurcacin son de gran utilidad. Las
2010; Strogatz, 1994). De esa manera, la estabilidad lo- soluciones continuas se denominan ramas. Las ramas
cal del estado estacionario (Xs) puede determinarse es- son curvas unidimensionales, cuando varan los valo-
tudiando las propiedades de la Ecuacin (1), cuando esta res de , el nmero de soluciones y en consecuencia,
se linealiza alrededor de Xs. Para este propsito se desa- de ramas, pueden cambiar. Nuevas ramas pueden surgir,
rrolla la matriz jacobiana. La estabilidad del equilibrio es pueden terminar o pueden cruzarse con otras en ciertos
caracterizada por los eigenvalores de esta puntos llamados puntos de bifurcacin de estado esta-
matriz. Para el caso en estudio, como la matriz jacobia- ble (Ajbar et al., 2012).
na est constituida por n eigenvectores linealmente in-
dependientes , la estabilidad local de la Los diagramas de bifurcacin se realizaron con ayuda
Ec. (1) puede ser concluida fcilmente. De acuerdo con del software Matcont, el cual permite hacer anlisis de
el resultado atribuido a Liapunov, las partes reales de los continuidad y determinar las bifurcaciones posibles del
eigenvalores de la matriz jacobiana , evaluada en la so- sistema de ecuaciones diferenciales ordinarias plantea-
lucin estacionaria ( ), determinan su estabilidad: si la do (Dhooge et al., 2006).
parte real de los eigenvalores son menores que cero, el
estado estacionario evaluado es considerado estable, de
lo contrario el estado estacionario es inestable.
Puntos de bifurcacin de estado estable
Los puntos comunes encontrados en diagramas de bi-
furcacin se presentan en la Tabla 1.
Diagrama de bifurcacin
Despus de contar con valores escalares del vector es
posible representar grficamente las soluciones de la
ecuacin (1) en un diagrama que muestra a en fun-

Tabla 1. Puntos de bifurcacin comunes en sistemas continuos

Punto de bifurcacin Principales caractersticas

Punto de inflexin En este punto 2 ramas de una solucin de estado estacionario se

Bifurcacin Fold unen formando tangentes limitantes.
Punto lmite estacionario Dos ramas nacen o se anulan entre s.
Nodo de silla Localmente no hay soluciones a un lado del punto y al otro lado hay
(SLP) dos soluciones.

Punto de bifurcacin estacionario En este punto dos y solo dos curvas con diferentes tangentes se
(SPB) o (BP) cruzan entre s.

Punto cspide estacionario Es el punto de contacto entre dos curvas de estado estacionario con
(SCP) la misma tangente.
Asociado con la aparicin de dos eigenvalores complejos conjuga-
Punto Hopf
dos , con > 0.
(HP)
Solo puede aparecer en sistemas de dimensin 2 o ms.
Adaptado de: Ajbar et al. (2012) y Seydel (2010).

Vol. 4 N. 1 pp. 15-23, enero-junio 2014

 Fermentacin ABE empleando Clostridium A.
20

Resultados y discusin estacionario (E.E) para cada modelo cintico. Adems,

se incluye el estado estacionario de lavado, donde se
Determinacin y seleccin de los puntos estacionarios
produce arrastre de la biomasa al flujo de salida del bio-
La solucin del modelo cintico descrito en la seccin rreactor y no se obtienen productos, los cuales sern
2.3.1, bajo condiciones de estado estacionario se mues- necesarios para acotar el anlisis de estabilidad.
tra en la tabla 2. Del mismo modo, en la tabla 3 se mues- Para este caso el transcurso de la fermentacin se fa-
tran las tres condiciones iniciales con las que se obtie- vorece a dilucin y concentracin de sustrato de 0,03
nen las mayores concentraciones de butanol en estado y 40 [g/L], respectivamente. Con respecto a la

Tabla 2. Condiciones iniciales y concentracin de butanol en estado estacionario

Condicin X S Caa Cba Cb Ca Ce Sf Cb (E.E.)

H D[h1]
Inicial [g/L] [g/L] [g/L] [g/L] [g/L] [g/L] [g/L] [g/L] [g/L]
1 0 40 0 0 0 0 0 7,9433E-05 0,03 40 0
2 0,10 40 0 0 0 0 0 7,9433E-05 0,10 40 8,58
3 0,10 20 0 0 0 0 0 7,9433E-05 0,03 20 8,95
4 0,01 40 0 0 0 0 0 7,9433E-05 0,03 40 10,63

Tabla 3. Valores de estado estacionario obtenidos para cada una de las

variables evaluadas

Variable E.E Lavado estables obtenibles cuando se vara la velocidad de di-

lucin y se mantiene constante la concentracin de
X [g/L] 5,3080 0 sustrato de alimentacin en 40 g/L. En la Figura 2 se
S [g/L] 0,0882 40 presentan los diagramas de bifurcacin.
Caa[g/L] 2,1996 0
Cba[g/L] 2,4522 0 En general, se puede observar que los diagramas de
Ca [g/L] 4,6642 0 bifurcacin exhiben un comportamiento similar y se
Cb [g/L] 10,6312 0 obtienen bifurcaciones comunes para sistemas conti-
nuos (Seydel, 2010). En efecto las figuras tienen dos
Ce [g/L] 0,3254 0
tipos de puntos: puntos Hopf (HP) y puntos de bifur-
H 1,0001E-07 1,0001E-07 cacin (SBP). Las grficas fueron divididas en 5 regio-
pH 6,9999 6,9999 nes para realizar el anlisis, donde las lneas continuas
simbolizan las ramas estables del sistema, y las pun-
teadas o discontinuas las inestables. De igual modo,
condicin de lavado, la velocidad de dilucin crtica de se analizar con base en 4 diagramas de bifurcacin
0,54 h-1 corresponde al rango (0,4-0,7 h-1) de los valo- representativos de todo el sistema: (i) biomasa, (ii)
res obtenidos para este tipo de sistemas (Bai, 2007; sustrato, (iii) butanol para el caso de los solventes y,
Godin et al., 1990; Mayank et al., 2012; Mahecha-Bo- (iv) cido butrico para el caso de los cidos.
tero et al., 2005; Zhang et al., 2001).
La primera regin est limitada por velocidades de di-
lucin . En esta zona existe multi-
Anlisis de estabilidad y diagramas plicidad de soluciones para la velocidad de dilucin,
constituidas por dos ramas, una estable y la otra ines-
de bifurcacin table. En la estable, las concentraciones mximas de
Se eligieron los diagramas y resultados ms represen- butanol son mayores a 10 g/L. Para esta misma re-
tativos, correspondientes a aquellos cuyo parmetro gin, la biomasa alcanza concentraciones mximas
de bifurcacin es la tasa de dilucin (D). En la tabla 2 del orden de 9 g/L.
(seccin 3.1) se resumen los valores necesarios para
iniciar este anlisis. En este caso en particular el inte-
rs es de determinar cuntos y cules son los estados

Revista Digital de la Facultad de Ciencias Naturales e Ingeniera de la UJTL

Castellanos et al., 2014
21

Por ejemplo para una concentracin de butanol igual a

10,02 g/L, la concentracin de biomasa es de 9,412 g/L.
Este comportamiento corresponde a los resultados re-
portados en estudios de biorreactores continuos (Godin
et al, 1990). As, en la mayora de los casos y con la cepa
Clostridium cuando se est en un biorreactor de clulas
libres se opera a bajas diluciones ( ) alcanzando
las mejores productividades de solventes. Anlogamen-
te, el sustrato alcanza conversiones mayores al 90%.

En contraste, en la segunda regin

(a) se puede observar una disminucin del 20% en la con-
centracin de butanol para un . Por otra
parte, se aumenta la concentracin de cidos en el me-
dio, en un porcentaje similar al del producto Cba=5,509
g/L, lo que afecta adems el crecimiento de la bioma-
sa, donde para estas mismas condiciones corresponde
a X=6,926 g/L. Adems la conversin del sustrato dis-
minuye a un 67% aproximadamente. En la rama inesta-
ble, todas las concentraciones son nulas, a excepcin de
la concentracin de sustrato que corresponde al mismo
valor de alimentacin. Esta rama describe el comporta-
(b) miento del sistema cuando se encuentra en condicin
de lavado.

En la tercera regin, para velocidades de dilucin com-

prendidas entre ; el comportamiento
de la rama estable e inestable siguen la tendencia en-
contrada en la anterior regin analizada. No obstante,
la rama inestable experimenta la aparicin de un pun-
to Hopf. Este punto confirma su inestabilidad y permite
pensar que para este valor de dilucin, aunque sea ines-
table, podran existir soluciones peridicas o ciclos lmi-
tes que deben considerarse en un estudio de inestabili-
(c) dad ms profundo. Para este punto, todos los valores de
la parte real de los eigenvalores tienen el mismo valor
que la tasa de dilucin ( ) y todos
son positivos menos el primero.

Finalmente, en la ltima regin ;

la rama estable disminuye hasta llegar a cero, a excepcin
de la del sustrato. Con el aumento de la velocidad de di-
lucin, se llega al valor de dilucin crtica ( )
para el sistema y al mismo tiempo a un punto de bifurca-
cin. En este punto, convergen las dos ramas que se ana-
lizaron a lo largo de las 4 regiones, donde una desaparece
(d) (la rama estable) y la otra pasa a ser estable, describiendo
de nuevo el lavado pero en esta ocasin, como una rama
Figura 2. Diagramas de bifurcacin para biorreactor empleando como
parmetro la velocidad de dilucin (a) concentracin de biomasa, (b) con- estable del sistema.
centracin de butanol, (c) concentracin de sustrato y (d) concentracin de
cido butrico

Vol. 4 N. 1 pp. 15-23, enero-junio 2014

 Fermentacin ABE empleando Clostridium A.
22

Conclusiones Godin, C., & Engasser, J. M. (1990). Two-stage contin-

uous fermentation of Clostridium acetobutylicum:
Con referencia a la velocidad de dilucin, el anlisis
effects of pH and dilution rate.Applied Microbiol-
de estabilidad demostr que la zona de estados es-
ogy and Biotechnology,33(3), 269-273.
tables posibles que favorecen el aumento en la con-
centracin de butanol est delimitado para aquellos Jang, Y.-S., Malaviya, A., Cho C., Lee J., Lee S. (2012).
valores menores a 0,078 h-1. Adems, el mismo anli- Butanol production from renewable biomass by
sis permite concluir que el valor ms adecuado de ve- clostridia. Bioresource Technology 123, 653663.
locidad de dilucin es de 0,03 sin sacrificar la eficacia
del proceso. Kamiski, W., Tomczak, E., Grak, A. (2011). Biobuta-
nol-Production and purification methods. Ecologi-
El anlisis de estabilidad permite concluir que su uti- cal Chemistry and Engineering. 18(1), 31-37.
lizacin simplifica y reduce el nmero de experimen-
tos a desarrollar, con miras a identificar las zonas es- Kirschner M. (2006). n-Butanol. Chemical Market Re-
tables del proceso de fermentacin ABE que permiten porter January 30February 5, ABI/INFORM Glob-
la produccin de altas concentraciones de solventes. al, p. 42.

Finalmente, la zona de estados estables encontrada Lee, S. Y., Park, J. H., Jang, S. H., Nielsen, L. K., Kim, J.,
concuerda con los valores experimentales reportados Jung, K. S. (2008). Fermentative Butanol Produc-
por varios autores. tion by Clostridia. Biotechnology and Bioengineer-
ing. 101(2), 209-28.

Referencias Mahecha-Botero, A., Garhyan, P., Elnashaie, S.S.E.H.

(2005). Bifurcation, Stabilization, and Ethanol Pro-
Ajbar, A., & Alhumaizi, K. (2011).Dynamics of the che-
ductivity Enhancement for a Membrane Fermen-
mostat: a bifurcation theory approach. CRC Press.
tor. Mathematical and Computer Modelling. 41,
Bai, F. W. (2007). Process oscillations in continuous 391-406.
ethanol fermentation with Saccharomyces cerevi-
Patakova, P., Linhova, M., Rychtera, M., Paulova, L., &
siae. Tesis doctoral. Universidad de Waterloo. Wa-
Melzoch, K. (2013). Novel and neglected issues of
terloo. 218 p.
acetonebutanolethanol (ABE) fermentation by
Dhooge, A., Govaerts,W., Y. A., K., Mestrom,W., y Riet, clostridia: Clostridium metabolic diversity, tools
A. M, a. S. B. Matcont and C. L. Matcont. (2006): for process mapping and continuous fermentation
Continuation toolboxes in Matlab. University Gent, systems.Biotechnology Advances,31(1), 58-67.
Blgica.
Paz I., Cardona C. (2011a). Influencia de las regiones
Donaldson, G. K., Huang, L. L., Maggio-Hall, L. A., de estabilidad del bioreactor en la produccin efi-
Nagarajan, V., Nakamura, C. E., Suh, W. (2007). Fer- ciente de alcohol carburante. Revista Colombiana
mentative production of four carbon alcohols. In- de Biotecnologa 8(1), 23-32.
ternational Patent WO2007/041269.
Paz, I., Cardona C. (2011b) Importance of stability
Fradkov, A. L., Miroshnik, I. V., & Nikiforov, V. O. study of continuous systems for ethanol produc-
(1999).Nonlinear and adaptive control of complex tion. Journal of Biotechnology 151, 43-55.
systems(Vol. 491). Springer.
Seydel, R. (2010). Practical Bifurcation and Stability
Garhyan, P., Elnashaie, S. S. E. H., Al-Haddad, S. M., Analysis.Third Edition. Springer, New York.
Ibrahim, G., & Elshishini, S. S. (2003). Exploration
Stockar, U., Valentinotti, S., Marison, I., Cannizzaro,
and exploitation of bifurcation/chaotic behavior of
C., Herwig, C. (2003). Know-how and know-why in
a continuous fermentor for the production of eth-
biochemical engineering. Biotechnology Advances
anol. Chemical Engineering Science, 58(8), 1479-
21, 417430.
1496.

Revista Digital de la Facultad de Ciencias Naturales e Ingeniera de la UJTL

Castellanos et al., 2014
23

Strogatz, S. H. (1996). Nonlinear dynamics and chaos, Yang, X., Tsao G. (1994). Mathematical modeling of in-
1994.Reading: Perseus Books. hibition kinetics in acetone-butanol fermentation by
Clostridium acetobutylicum. Biotechnology Progress
Tashiro, Y., & Sonomoto, K. (2010). Advances in butanol 10(5), 532-538.
production by clostridia.Current Research, Technolo-
gy and Education Topics in Applied Microbiology and Zhang, Y., & Henson, M. A. (2001). Bifurcation analysis of
Microbial Biotechnology (Microbiology Book Series, continuous biochemical reactor models.Biotechnol-
number# 2). FORMATEX, Spain, 1383-1394. ogy Progress,17(4), 647-660.

Yang, X. (1990). A novel extractive fermentation for the

production of Acetone-butanol by Clostridium aceto-
butylicum. PurdueUniversity.

Vol. 4 N. 1 pp. 15-23, enero-junio 2014