FACULTAD DE INGENIERÍA

Departamento de Ingeniería Química y Ambiental

SISTEMA DE CONTROL DE UN PROCESO DE FERMENTACIÓN AERÓBICA PARA LA

OBTENCIÓN DE ETANOL

En el presente documento se busca el diseño de un sistema de control mediante lazos

de tipo realimentado, cuyo objetivo es mantener las especificaciones de una corriente
etanol-agua, producto de una reacción de fermentación.

En primer lugar, la fermentación es un proceso que puede ser definido a partir de la

obtención de un producto a través del empleo de un agente biológico. El proceso de
fermentación más ampliamente conocido es la fermentación alcohólica, de acuerdo a la cual se
tiene una corriente de mosto o sustrato que es rica en azúcares. Dicha corriente se pone en
contacto con un cultivo generalmente bacteriano para el cual un subproducto metabólico de la
degradación del azúcar sea el alcohol deseado.

Previo a la fermentación pueden tenerse operaciones de prensado o molienda que

permitan obtener el mosto. Lo anterior puede efectuarse a partir de cultivos de maíz, papa,
yuca, sorgo, remolacha y caña de azúcar. Éstos se seleccionan debido principalmente a su alta
concentración de carbohidratos que los entes biológicos son capaces de degradar. El producto
de fermentación es generalmente una mezcla entre biomasa, agua y etanol. La biomasa
corresponde a los cuerpos celulares cuya densidad y composición los hace susceptibles a ser
retirados mediante filtraciones. Una gran proporción de la mezcla resultante tras remover los
cuerpos celulares está conformada por agua, por lo que se tiene entonces una solución
etanol-agua que generalmente es sometida a rectificación en una columna. Dependiendo de
las, condiciones de proceso, pueden alcanzarse concentraciones de etanol mayores o
menores de acuerdo al destino del alcohol. El proceso de producción de etanol utilizado
actualmente en Colombia aún no es el más eficiente, la etapa del proceso que generalmente
limita la eficiencia es la fermentación, puesto que es la etapa con mayores pérdidas energéticas
y mayor producción. La baja eficiencia se debe a la técnica utilizada en el cultivo, lo que implica
la selección de la cepa bacteriana y control de las bacterias para la fermentación, las cinéticas
de reacción para procesos de fermentación suelen ser de carácter inestable y oscilante, lo que
dificulta el control del proceso. Un sistema de control apropiado debe encargarse de mantener
en los niveles adecuados, la temperatura, la adición de sustratos y el nivel de acidez en el
tanque de fermentación, garantizando un ambiente ideal para el crecimiento de las bacterias
encargadas de fermentar el jugo, aumentando la concentración y disminuyendo el tiempo de
producción de etanol. Lo anterior plantea la necesidad de mejorar las estrategias de control del
proceso de fermentación para hacerlo más eficiente. [1]

1. Objetivos específicos
1.1​. Diseñar un controlador PID para regular la temperatura en el interior del fermentador en un
valor igual a 32 °C que opere en un lazo realimentado.
1.2​. Diseñar un controlador PID para regular la razón de reflujo en el rectificador que opere en
un lazo realimentado.

La temperatura óptima de fermentación de la ​Saccharomyces Cerevisiae ​en condiciones

aeróbicas es 32 °C y el pH óptimo es de 4.5. En general, si se conocen los parámetros de
ajuste de la cinética de la reacción de fermentación, existen correlaciones que permiten derivar
el potencial de hidrógeno a través del tiempo.

2. Variables asociadas
Una vez se tiene el volumen de control para el cual se desarrollarán los lazos de control, deben
definirse las variables involucradas y su carácter dentro de la simulación. Las variables
aparecen registradas a continuación. Como se definió en los objetivos del proyecto, las
variables controladas corresponderán a la temperatura del medio y el nivel del reactor.

Variable Descripción Manipulables Perturbaciones

h Altura de la masa de reacción

T​r Temperatura del reactor

Reactor c​P Concentración de etanol

c​S Concentración de sustrato

c​X Concentración de biomasa

c​O2 Concentración de oxígeno

Corriente

F​i Flujo de alimentación

1 c​S, in Concentración del sustrato

T​in Temperatura de alimentación

 F​ag Flujo de refrigerante
2
T​in,ag Temperatura de ingreso del
refrigerante

3 T​ag Temperatura del refrigerante

en la chaqueta

4 F​e Flujo de salida

A continuación se presenta el análisis de grados de libertad del sistema, que permite

determinar el número de variables controladas en el establecimiento de los lazos de control.

N D = N V − N E = 13 − 7 = 6
NM = ND − NP = 6 − 3 = 3

Donde:
N​V.​ Número de variables
N​E.​ Número de ecuaciones
N​P. ​ Número de perturbaciones
N​D.
N​M. ​ Número de variables manipuladas

3. MODELO DINÁMICO

3.1. Parámetros

A1​ ​ = 9,5 x 108​ KO2​ ​ = 8,86 mg l-1 ​ V = 10 l

A2​ ​ = 2,55 x 1033 ​ KP​ ​ = 0,139 g l-1​ V​j​ = 0,5 l
AT​ ​ = 0,2 m2​ KP1 ​ ​ = 0,070 g l-1​ Y​O2​ = 0,970 mg / mg
-1​
​
Cheat,ag ​ = 4,18 J g​ K​
-1
KS​ ​ = 1,030 g l​ -1
ΔHr​ ​ = 518 kJ / mol O2​
-1​
Cheat,r​ ​ = 4,18 J g​ K​ -1
KS1 ​ ​ = 1,680 g l-1 ​ μ​O2​ = 0,5 h-1
​
Ea1​ ​ = 55000 J mol-1 ​ KT​ ​ = 3,6 x 10​5​ J h​-1​ m​-2​ K​-1 μ​P​ = 1,790 h​-1
Ea2 ​ ​ = 220000 J mol-1 ​ R​SP​ = 0,435 ρ​ag​ = 1000 g l-1​
(k​l​a)0​ ​ = 38 h-1
​ R​SX​ = 0,607 ρ​r​ = 1080 g l-1
​

3.2. Entradas del sistema

F i = F e = 0, 5 l h−1
T in = F e = 25 °C
cS,in = 60 g/l
T in,ag = 15 °C

3.3. Control de temperatura

La dependencia de la concentración de oxígeno de equilibrio con la temperatura en agua
destilada está dada por la siguiente ecuación empírica ( ).

c*O = 14, 6 − 0, 3943T r + 0, 007714T r2 − 0, 0000646T r3

2

(k l a) = (k l a)0 (1, 024)T r −20

La velocidad de consumo de oxígeno es:

cO2
rO2 = μO2 Y 1 cx K
O2 O + cO
2 2

La expresión de la velocidad de crecimiento máximo involucra la velocidad de incremento con

la temperatura y el efecto de la desnaturalización por calor.

μx = A1 exp [ E a1
R (T r +273) ] − A exp [
2
E a2
R (T r +273) ]
Los balances de materia para la biomasa, el producto, el sustrato y el oxígeno disuelto son
expresados por las ecuaciones

dcX cS Fe
dt = μx cx K e−K p cp − V cx
S + cS

dcp cS Fe
dt = μp cx K e−K p1 cp − V cp
S1 + cS

dcS 1 cS cS Fi Fe
dt =− μ cx
RSX x K S + cS e−K p cp − 1
RSP μ cx K e−K p1 cp − V cS,in − V cS
p S1 + cS

La concentración de oxígeno disuelto en el medio de reacción depende de la cantidad de

oxígeno que entra por transferencia de masa y la velocidad de consumo de oxígeno de las
levaduras para la reacción de fermentación

d cO2
dt = (k l a)(cO2* − cO 2 ) − rO 2

Los balances de energía para el reactor y la chaqueta están dados por las ecuaciones ( ) y ( ),
respectivamente.

dT r Fi Fe rO2 ΔH r K T AT (T r −T ag )
dt = V (T in + 273) − V (T r + 273) + 32 ρr C heat,r + V ρr C heat,r
dT ag F ag K T AT (T r −T ag )
dt = V j (T in,ag − T ag ) + V j ρag C heat,ag

3.4. Control de nivel

dm ︿ ︿
dt = mi − me

d (ρe Ah)
dt = ρi F i − ρe F e

Se supone densidad constante a lo largo del tiempo para el flujo de entrada (glucosa) y el flujo
de salida para obtener la expresión de la variación de altura controlada por el flujo de entrada.

dh ρi F i −ρe F e
dt = ρe h

3.5. Análisis dinámico del proceso en lazo abierto

Luego de haber introducido las ecuaciones constitutivas del sistema, definir las variables de
control y definir las variables de proceso, procedemos a realizar el análisis de los resultados
que obtenemos con las variables de proceso en estado estable sin perturbación, de esta
manera encontramos los valores que introducimos como condiciones iniciales en el bloque
integrator ​del sistema en ​Simulink, ​para la temperatura:
La temperatura constante alcanzada es de 26,04 °C, como se mencionó antes cuando los
valores de entrada del flujo de sustrato y de refrigerante son constantes.

La altura constante alcanzada es de 0.2646 cm, como se mencionó antes cuando los valores
de entrada del flujo de sustrato y de refrigerante son constantes.

4. P&ID
El diagrama P&ID del proceso se presenta en el anexo A.
4.1. Válvulas
Las válvulas empleadas para el control de la temperatura y el nivel en el interior del reactor, TV
100 y LV 101 respectivamente, responden a una señal neumática. Son válvulas de globo de
conexión estándar sin divisor de flujo con característica lineal. Ver anexo C.

Diámetro nominal de la tubería (NPS) = 1 in

Ganancia de la válvula (K​VS​) = 0,1

Figura

La válvula contenida en la corriente 4 es manual y se acciona en caso de necesitar restringir el

flujo al tanque de almacenamiento (por llenado). Las tres válvulas involucradas en el proceso
son normalmente cerradas.

Figura .

4.2. Sensor transmisor

5. Análisis dinámico del proceso

Realizamos el bloque del controlador,insertando un PID que recibe el error que será la
diferencia entre el valor de referencia y la respuesta de la variable controlada resultante del
sistema de ecuaciones obtenido de las ecuaciones constitutivas del biorreactor. La señal
transmitida por el PID ingresa a un actuador, este bloque es el que indica el porcentaje de
apertura que se realizará en la válvula de control para regular los valores normales alterados
por la perturbación; la salida del actuador se multiplica por la ganancia de la válvula que genera
los valores que pretender regular las variables a controlar.

Generamos la perturbación en el flujo de entrada del refrigerante, y posteriormente en el flujo

de entrada del sustrato como se muestra a continuación.

Iniciamos el controlador para la altura con una perturbación generada por el flujo de alimento
de sustrato, observado la respuesta del controlador y verificando que se ajusta a la referencia
esperada, cabe aclarar que en el bloque del controlador escogimos un PID que actúa con la
acción PI pues la acción derivativa es cero.

1. PERTURBACIÓN DEL FLUJO DE ALIMENTACIÓN DE SUSTRATO

Respuesta del controlador en el nivel del reactor cuando se reduce el flujo de alimento
de 0,5 g/h a 0,25 g/h
Temperatura del reactor con la perturbación generada.

Concentración de sustrato en el tiempo a 50 h

Concentración de etanol en el tiempo a 50 h

2. PERTURBACIÓN DEL FLUJO DE ALIMENTACIÓN DE REFRIGERANTE

Respuesta del controlador a la temperatura del reactor aumenta el flujo de refrigerante
de 0,5 g/h a 1 g/h

Nivel del tanque cuando se genera la perturbación

6. Estructura alternativa de control
Otra forma de controlar el nivel del reactor es mediante la manipulación de la válvula a la salida
del mismo, por lo tanto se cambia el lazo, tomando como elemento final de control la válvula ya
mencionada. Por su parte, la válvula de acción manual se sitúa a la entrada del reactor. Ver
anexo 5.

7. Conclusiones:
● El proceso puede ser controlado mediante lazos realimentados empleando
controladores tipo PID. Para las condiciones seleccionadas se lograron estados
estacionarios alrededor de las referencias definidas.
● Existen alternativas con respecto al tipo de controlador empleado (como el controlador
difuso presentado por Hurtado et. al [1]) Sin embargo, su complejidad y dificultad en el
modelamiento debido a la inestabilidad del sistem, confieren al control PID ser una
solución sencilla y económica para un amplio número de plantas.
Nomenclatura

A​1​, A2​ Factor pre-exponencial en la ecuación K​P1 Constante de la inhibición de

de Arrhenius fermentación por el etanol (g/l)

​
Cheat,ag Capacidad calorífica del agente de K​S Constante en el término de sustrato para
enfriamiento (J g​-1​ K​-1​) el crecimiento (g/l)

​
Cheat,r Capacidad calorífica de la masa de K​S1 Constante en el término de sustrato para
reacción (J g​-1​ K​-1​) la producción de etanol (g/l)

​
cO2 Concentración de oxígeno en la fase K​T Coeficiente de transferencia de calor
líquida (mg/l) (J h​-1​ m​-2​ K​-1​)

​ ​*
cO2 Concentración de equilibrio del oxígeno ​
rO2 Velocidad de consumo de oxígeno
en la fase líquida (mg/l) (mg l​-1​ h​-1​)

cp​ Concentración del producto, etanol (g/l) R Constante universal de los gases
(8,314 J mol​-1​ K​-1​)
cS​ Concentración del sustrato, glucosa (g/l)
​
RSP Relación de etanol producido por
​
cS,in Concentración de glucosa en el flujo de glucosa consumida para la fermentación
alimento (g/l)
​
RSX Relación de células producidas por
cX​ Concentración de biomasa, levadura glucosa consumida para el crecimiento
(g/l)
t Tiempo
E​a1​, Ea2
​ Energía de activación aparente para la
reacción de crecimiento y la de ​
Tag Temperatura del agente refrigerante en
desnaturalización, respectivamente la chaqueta (°C)

​
Fag Flujo de agente de enfriamiento (l h​-1​) Tin​ Temperatura del flujo de sustrato que
ingresa al reactor (°C)
Fe​ Flujo de salida del reactor ( l h​-1​)
​
Tin,ag Temperatura del agente refrigerante que
Fi​ Flujo de sustrato entrante al reactor ingresa en la chaqueta (°C)
(l h​-1​)
Tr​ Temperatura en el reactor
(kl​ a)
​ Producto del coeficiente de
transferencia de masa para el oxígeno y V Volumen de la masa de reacción (l)
el área específica de la fase gaseosa
(h​-1​) V​j Volumen de la chaqueta (l)

(kl​ a)​
​ 0 Producto del coeficiente de Y​O2 Factor de rendimiento de la biomasa en
transferencia de masa a 20 °C para el el oxígeno (mg/mg)
oxígeno y el área específica de la fase
gaseosa (h​-1​)

K​O2 Constante de consumo de oxígeno (g/l)

K​P Constante de la inhibición del

 crecimiento por el etanol (g/l)

Símbolos griegos

μ​O2 Velocidad específica de consumo de ΔH​r Calor de reacción de la fermentación

oxígeno máxima (h​-1​) (kJ/ mol O​2​ consumido)

μ​P Velocidad específica de fermentación ρ​ag Densidad del agente de enfriamiento (g/l)
máxima (h​-1​)
ρ​r Densidad de la masa de reacción (g/l)
μ​X Velocidad específica de crecimiento
máxima (h​-1​)

REFERENCIAS
[1]