Diseño de

Biorreactores
Alfredo Palomino Infante
Profesor Principal, FQIQ,UNMSM
[email protected]
 Introducción
En este curso se presentan los fundamentos del
diseño de biorreactores, basado en la experiencia del
autor en la enseñanza del curso de Ingeniería de
Diseño, de la maestría de Biotecnología en la Facultad
de Farmacia y Bioquímica de la Universidad Nacional
Mayor de San Marcos.
Esta presentación se inicia con una aproximación
jerárquica al modelamiento de biorreactores,
continuando con los conceptos de diseño y el análisis
de biorreactores ideales.
Para ayudar la comprensión de los modelos derivados,
se han incluido problemas resueltos, así como una
lista de abreviaturas usadas en el trabajo.
 Abreviaturas Usadas

C Concentración de células de levaduras suspendidas en

el líquido(g/L)
ks Constante de saturación(g/L)
k’s Constante de saturación(g/L)
P Concentración de producto(g/L)
v Velocidad superficial(cm/min)
S Concentración de sustrato(g/L)
S0 Concentración de entrada sustrato(g/L)
X Concentración de células o biomasa (g/L)
Yx/s Factor de rendimiento de producto(g células/g sustrato)
YP/s Factor de rendimiento de producto(g producto/g sustrato)
Z Longitud del reactor tubular(cm)
k1 Constante empírica de inhibición de producto(L/g)
rx velocidad de formación de células de levaduras(g/L h)
rp velocidad de formación de productos(g/L h)
rs velocidad de consumo de sustrato(g/L h)
t tiempo (horas)
ΚLa coeficiente global de transferencia de masa (1/hora)
ε Fracción vacía de la columna empacada
μ velocidad específica de crecimiento celular(1/h)
μm velocidad específica máxima de crecimiento celular(1/h)
ν velocidad específica de formación de producto(1/h)
νm velocidad específica máxima de formación de producto(1/h)
β factor de rendimiento= 1/Yxs
τ tiempo espacial( h ó min)
D Coeficiente de dilución, (1/h).
 ESTRUCTURA JERARQUICA DE SISTEMAS
BIOLOGICOS
Primer nivel: células Segundo nivel: agregados
de células

Mecanismo de
crecimiento-factores Efectos de interacción
condicionantes intercelular

Tercer nivel
Efectos de
transporte de
masa, energía y
cantidad de
movimiento
ESTRUCTURA JERARQUICA DE SISTEMAS BIOLOGICOS

Cuarto nivel: Micro

hidrodinámica, tipo de
flujo local
aglomerados

Quinto nivel: Macro

hidrodinámica, tipo de
flujo global

Zonas
muertas
 Biorreacciones

Fermentación de bacterias, levaduras ,

hongos, etc, para la producción de
células, productos intracelulares y
extracelulares.
Reacciones enzimáticas

Sistemas basados en células de plantas y

animales.
 DE QUE TRATA EL DISEÑO DE
BIORREACTORES

Elegir la mejor configuración/proceso/programa

de operación (batch, fed-batch, continuo) para
un objetivo particular.

Estimar y optimizar el requerimiento de tamaño/

tiempo para un rendimiento máximo o
productividad.
Consideraciones en el diseño

• ¿Cuál es la naturaleza del producto?

• Naturaleza de las células o enzima.

• Células o enzimas inmovilizadas o
libres
• Fases presentes (ej. solventes)
• Esterilidad requerida.
• Aireación.
• Control.
 ALTERNATIVAS DE DISEÑO

Elegir condiciones de
operación

1. Modo de operación: Contínuo o batch

2 Mezclado/agitación, remoción de calor
3. Sistemas de aireación
4 Retención y mejoramiento de biomasa
5. Escala de operación, reglas de ingeniería
Pasos del diseño

1.Estequiometría
2. Cinética
3. Medición y Balances
4. Aspectos de ingeniería
 ANÁLISIS DE BIORREACTORES

Considerando que los modelos cinéticos simples para

procesos enzimáticos y microbianos, pueden ser
usados para diseñar los tipos básicos de biorreactores.
Tomando en cuenta las tres configuraciones básicas
de biorreactores y modos operacionales (batch, mezcla
completa de flujo (CSTBR) y flujo tapón (PFBR). En
estos modelos una suposición clave es que la célula
microbiana o enzima está disuelta y son físicamente
indistinguibles del líquido; es decir se asume un medio
pseudo homogéneo
 FERMENTADOR BATCH
CON CINETICA DE MONOD

Estamos interesados
solamente en las células (X)
y sustrato (S) . En
consecuencia deseamos
calcular como varían las
concentraciones de X y S
con el tiempo.
Considerando un
fermentador de volumen
constante y haciendo un
balance de células tenemos:
• Un balance de
sustrato suponiendo
rendimiento, Yx/s
constante da:

• Estas ecuaciones se
resuelven con las
siguientes condiciones:
 Solución bajo restricciones
Si la relación (ks/
S)<<1, se
demuestra que la
ecuación(1) queda
como sigue:
Las ecuaciones (4) y
(3) se resuelven y
dan como
resultado:
 Resolviendo
• La solución
general de la
ecuación se
obtiene
primero
resolviendo la
ecuación (2):

• Reemplazando
(6) en (1) y
aplicando (3)
se demuestra
queDonde:
la solución
es
 Ejemplo de Consumo de Sustrato

Los datos adimensionalizados

X/Xo S/So
de consumo de sustrato que se μm t,
presentan a continuación, ( si es
muestran que la concentración logarítmico)

de sustrato permanece 2.7 1 1

prácticamente constante 148 1 5
22,026 0.98 10
durante los primeros estadios de 59,874 0.94 11
crecimiento y luego decae 162,755 0.837 12
442,413 0.558 13
rápidamente. Consideramos por 729,416 0.271 13.5
ejemplo un sistema con X∞ /
Xo=10 6:
 Estimación de la velocidad de
reacción
En la primera fase de crecimiento, X aumenta
rápidamente, pero S no varía apreciablemente, de modo
que
rx aumenta con X hasta que disminuya el sustrato,

máximo
 FERMENTADOR DE FLUJO DE TAPON
Este biorreactor se
Fig. Esquema del biorreactor
caracteriza por mantener PFBR
flujo de tapón y no
admite retromezcla ni v 0, X 0 S 0 v0, X. S
gradiente radial a lo
largo del tubo. Se
considera flujo El balance de células da
estacionario de modo
que las composiciones a
lo largo del biorreactor
siempre son las mismas,
es decir, no varían con el
tiempo.
La relación vo/A es la velocidad superficial media (u) del
fluido, la cual se considera independiente de la posición
axial (Z). Si hacemos dZ=udt la ecuación (9) se
convierte en una ecuación equivalente a la ecuación (1),
toda vez que (Z/u) = (L/u) = t = tiempo de retención
global. Además, un balance de sustrato entre la entrada
y salida Z proporciona

De este modo, la evolución de la fermentación de un flujo

de tapón, en términos de coordenadas espaciales
transformadas, viene a ser idéntico al tiempo de
evolución en un reactor batch.
FERMENTADOR DE FLUJO TIPO TANQUE
AGITADO, CSTBR (Quimiostato)

Es un biorreactor de flujo Asumiendo estado estacionario

continuo de mezcla y volumen constante. El
completa (CSTBR), donde balance Global de sustrato da
el fluido que sale del
biorreactor tiene la misma
composición que el fluido
que se halla dentro, dada
la asunción de mezcla
Un balance de células
perfecta.
da
 Formas Alternas
En la expresión anterior, rx es la velocidad de
producción celular por unidad de volumen. La ecuación
(12) puede escribirse de las siguientes formas:

donde

Sabemos por definición que

• Desde que el proceso de crecimiento celular es
auto catalítico, es perfectamente posible
sustentar el crecimiento celular sin un
suministro continuo de células, esto es con una
alimentación estéril ó X0=0, se tiene

Esto significa que la velocidad específica de

crecimiento
debe ser igual a la velocidad de dilución para
estado estable,
independiente de la cinética precisa
Así, asumiendo cinética de Monod

donde S es la concentración del sustrato dentro del

fermentador continuo, igual a la concentración que deja el
reactor. Resolviendo para S la ecuación (15) obtenemos:
 Fig. Forma de la
solución Reemplazando la
ecuación (16) en
S, X (11) y haciendo
X X0=0 obtenemos:
S0

X0=Yx/sS0

S
μm
D*

D
Implica que hay una conversión completa se sustrato a
células
Hay un límite superior para D (este es el límite
inferior del tiempo de retención) desde X O
cuando

Para la mayoría de los

sistemas
..

Este es el punto de lavado o ( WASHOUT) donde la

velocidad de remoción de células excede la velocidad
máxima de producción en el fermentador. Debería
notarse que en la mayoría de las situaciones y las
curvas mostradas arriba tienen pendientes muy
pronunciadas en la región donde D=μm.

El resultado de la ecuación (16) es remarcable, puesto

que la concentración de salida de sustrato no depende de
la concentración de entrada So, pero sí de los parámetros
cinéticos y de D. Se dice que el proceso es limitado por
nutrientes. La concentración celular depende de S0 y
cuando
La productividad máxima se presenta en dP/dD=0.
A partir de la ecuación (18) se demuestra que
el valor de D óptimo viene dado por al ecuación:
Combinando la ecuación (19) con la
ecuación (17)
se demuestra también:

En la práctica, ks>>So. Por tanto:

 COMPARACIÓN DE UN FERMENTADOR BATCH
CON UN FERMENTADOR CONTINUO

A partir de las ecuaciones (1) y

(2),
considerando rendimiento, Y,
constante tenemos:

Para una operación batch integramos la ecuación (21),

donde tb es el tiempo de retención para una extensión de
reacción dada:
Para una operación en CSTBR el balance de
materia da:
(i)

(ii)

Usando (ii) junto con (3a), se tiene

1/r B 1/ B
batc r
h C B’
C
B’

CSTBR
D A
D A
S S0 S S
Fig. b 0

área ABCD ecuación 22; área AB’CD ecuación 23

Fig. Comparación de las ecuaciones 22 y 23

 Interpretación de las
ecuaciones (22) y (23)
• Es conveniente considerar aquí la Ec(22)
como la integral bajo la curva (1/r) Vs. S, es
decir el área ABCD en la figura
correspondiente. La Ec. (23) es el área del
rectángulo AB CD, es decir 1/r versus S-So.
Estos dos ejemplos muestran claramente
que AB CD es menor que ABCD. Es decir el
proceso continuo requiere un tiempo de
retención ( o volumen para una producción
dada; el otro representa un caso donde lo
inverso es cierto
 Cinética de Monod
A partir de la ecuación
(1):
(24
)

reemplazando en esta ecuación el

valor de X
deducido de la ecuación (6) o (12)
obtenemos:
(25
)
El cual tiene la forma ilustrada.
Para una operación CSTBR el balance de
materia da:

dond
e
El máximo se obtiene en dr/dw =0.
Para X0 = 0, se demuestra que
el valor óptimo es:
Si 1/k >> 1 entonces
obtenemos:

Si k << 1,
obtenemos:
 EJEMPLO DE APLICACIONES

Ejemplo 1. Un cultivo de E. Coli se deja crecer sobre

lactosa en un reactor de flujo con mezcla (V=1 litro).
Usando varias velocidades de flujo y una alimentación
constante de S0= 160 mg/L se obtienen los siguientes
resultados:

Flujo(L/h) Sustrato S(mg/L) Conc. células

0.2 4 15.6
0.4 10 15.0
0.8 40 12.0
1.0 100 6.0

Encuentre la expresión de velocidad que representa este

crecimiento para X0=0.
Solución: Para X0=0 hemos deducido las ecuaciones
14,16 y 17. La ecuación 17 representa el crecimiento celular,
para el cual determinaremos los parámetros Yxs, D, ks, μm .
La ecuación 12 la llevamos a la forma de un línea recta:
X= Yxs(S0-S)

X
S0-S X
156 15.6
150 15.0
Pendiente=Yxs= 0.1
120 12.0
60 6
S0-
S
La ecuación (16) es transformada a la siguiente forma:

A partir de la información dada elaboramos la

siguiente tabla:
1/D 1/S
0.06 0.25
0.07 0.10
0.08 Será posible obtener ks
0.025 y mm a partir de esta
información ?
0.17 0.01
 Ejemplo 2.

La E. Coli crece sobre manitol con la siguiente

velocidad

donde Yx/s= 0.1 g. de células/ g. de manitol.

Determine la concentración de salida de las

células del reactor, cuando el flujo al reactor
CSTBR (V=5 m3) es de 1 m3/hora con una
concentración So= 6 g./ m3.
 Solución

En condiciones estacionarias y si X0=0 obtenemos

de la ecuación 15:

Por analogía con la expresión dada en este

problema tenemos ks=2 y μm= 1.2. Luego, de la
ecuación (17):
X=0.57
 Ejemplo 3.
En su tesis publicada en 1948, Monod propuso su
ecuación sustentado con datos experimentales de cuatro
corridas batch para el crecimiento de un cultivo
bacteriano de lactosa. Ajuste los siguientes datos a la
ecuación de Monod.
Medic Δt,horas S promedio X X promedio

1 0.54 137 15.5-23.0 19.5

2 0.36 114 23.0-30.0 26.5
3 0.33 90 30.0- 38.8 34.4
4 0.35 43 38.8-48.5 43.65
5 0.37 29 48.5-58.3 53.4
6 0.38 9 58.3-61.3 59.8
7 0.37 2 61.3-62.5 61.9
 Solución:

La ecuación 1 puede escribirse en la forma lineal siguiente

t ΔX/Δt dX/dt (del

De donde podemos método)
1 19.5 17.0
determinar los 2 26.5 23.5
parámetros del 3 34.4 31.0
modelo, para ello 4 43.65 39.0
determinamos previa- 5 53.4 49.0
mente las derivadas 6 59.8 57.0
dX/dt 7 61.9 61.0
utilizando el método
de
Para la gráfica de la ecuación B construimos la siguiente tabla:

1/X X/(dX/dt) 1/S

Se obtiene los parámetros:
0.051 0.867
0.007 0.038 0.893 μm= 1.156 ks= 1.965
0.008 0.029 0.899 La ecuación de Monod es :
0.011 0.023 0.897
0.023 0.019 0.931
0.034 0.017 0.969
0.111 0.016 0.976 T S μ
0.500 1 137 1.14
2 114 1.13
3 90
1.13 4 43
1.10 5 29
1.08 6 9
0.95 7 2
Se indicó que:

Llevando esta ecuación a la forma lineal se

tiene:
t, horas CA dCA/dt dCA/dt
0 1.0 1
0.84 -0.16 -0.16 2 Del método
0.68 -0.16 -0.155 3 de las áreas
0.53 -0.15 -0.15 4 equivalentes
0.38 -0.15 -0.145 5
0.27 -0.11 -0.11 6
0.16 -0.11 -0.077 7
0.09 -0.07 -0.063 8
0.04 -0.05 -0.028 9
0.018 -0.022 -0.017 10
0.006 -0.012 -0.007 11
0.0023 -0.0037 -0.002
La siguiente tabla muestra los resultados para la
gráfica de la ecuación (A)

1/(dCA/dt) 1/CA
-6.25 1 Del gráfico se obtiene
-6.50 1.2
-6.70 1.5 pendiente= -1
-6.90 1.9 intercepto= -5
-9.10 2.6
-13.0 3.7 De donde se estiman las
-15.9 11.1 constantes
-35.70 25 k3= 20 con CE= 0.01
-58.80 55.6
km= 20 * 0.01 = 0.2
-142.90
166.7 -500
434.8
 Ejempl
o
La bacteria E coli crece constantemente en glucosa según :

con Yx/s= 0.1. Que velocidad de alimentación de Glucosa

(S0= 6 kg/m3) a un reactor CSTBR (V= 1 m3) podría originar
una velocidad máxima de consumo de glucosa y la máxima
producción de células de E coli.
Solución: En este caso aplicamos la ecuación 18 deducida
anteriormente. Con ks=0.4 y μm= 40
Reemplazando datos obtenemos:
D= 40*0.75=30; luego: vo = D*V = 30*1 = 30 m3/h
 Referencias
• Doran Pauline M , Principios de Ingeniería de los Bioprocesos. Edit. Acribia
S.A. España.
• Bailey, J. y Ollis D. Biochemical Engineering Fundamentals, 2nd, edit. Mc
Graw Hill Book Co, 984, USA, (1986).
• Palomino, Alfredo, Ingeniería de Reacciones Bioquímicas, Apuntes, FQIQ,
UNMSM, 1999.
• Fogler, Scott. Elements of Chemical Reaction Engineering, Edit. Prentice
Hall, USA, 1985.
• Schuler, M. and Weigand, W.,Editors, Biochemical Engineering V, Annals
of the New York Academy of Sciences, Vol 506, USA, 1987.
• Biotechnology Progress, Edited by American Institute of Chemical
Engineers, AIChE, USA, 1999. Bimontly periodical.
• Kafarov, V., Vinarov A, y Gordeev L. ICE, Part. 1 pp. 615, (1987)
• Kafarov V., Vinarov A, y Gordeev L. Int. Chem Eng, Part 2, pp. 14 Jan
(1988).
• Ciurlizza G. y Delgado J., Chem Engng., pp 115, 1997.
• Gencer M. And Mutharasan. Biotech and Bioeng . pp. 2243, 1983.
• Plowman R. Chem Eng. pp. 74, Mech (1997).
• Davis, Mark . Numerical Methods and Modeling for Chemical Engineers,
John Wiley and sons USA, (1984).
• Jenson V, and Je reys. Mathematical Methods in Chemical Engineering,
• Hanna O. y Sanball O. Computational Methods in Chemical Engineering,
Prentice Hall Inc. USA, 1995.
• Rieder W. y Busby H. Introductory Engineering Modeling, USA, 1986. John
wiley & Sons.
• Math Cad User´s Guide, Math Cad 7, Professional Mathsoft, Inc. USA.
http:/www Mathsoft con USA, 1997.
• Wylie Ray, Matemáticas Superiores para Ingeniería, 4ta. Edición. Mc
Graw Hill Book Co, México, 1982.
• Kreyszig E. Advanced Engineering Mathematics, Sixth. Edit, John Wiley
& Sons, USA, 1988.
• Marquardt, W. Int. Chem Eng., vol. 34, Nº 1, pp. 28, (1994).
• Reverchon E. AIChE Journal, Jun 1996, vol 42, Nº 6 pp. 1765 (1996)
• Shamel R. And chow J. Chem Eng. Prog. Nº 85 Nº 12, pp. 33 (1989).
• Colt W. Chem Eng. Prog. Vol 90, Nº 12, pp. 23, (1994).
• Buchh volz, K. Int Chem Eng., Vol 32 Nº 1 pp 1, (1992)
• Kunigita, E. , Nishitani H. y Wan Y. Int. Chem Engng., vol 33, Nº 2, pp. 289,
(1993).
• Izumoto E, Fukuda H. y Nakanishi, E. Int. Chem Eng., vol 32, Nº 3, pp. 570,
(1992).
• Schmidt Ch. y Schuler H. Int. Chem. Eng. Vol 33 Nº 4, pp 559, (1993).
• Dylag M y Talaga J., Int. Chem Eng., Vol 34 Nº 4, pp. 539, (1994).
• Dylag M y Talaga J., Int. Chem. Eng. Vol, 34, pp. 552 (1994).