TECNOLÓGICO NACIONAL DE MÉXICO

INSTITUTO TECNOLÓGICO DE LA PAZ

Ingeniería de Biorreactores

Serie 4

Isaac Mackliz Villegas

Dania Dayali Montaño Manríquez

Benito Ramírez Barrera

Ing. Jorge Humberto Silva Ojeda

Solución

tiempo (min) DO (mg/L) dCL/dt dCl/dt + qo2X

-1 3.3

Air off 0 3.3

1 2.4

2 1.3

3 0.3

4 0.1

5 0

Air on 6 0
 7 0.3 0.3 0.9971

8 1 0.7 1.3971

9 1.6 0.9 1.5971

10 2 1.1 1.7971

11 2.4 1.3 1.9971

12 2.7 1.4 2.0971

13 2.9 1.5 2.1971

14 3 1.5 2.1971

15 3.1 1.6 2.2971

16 3.2 1.6 2.2971

Gráfico del proceso

3.5

2.5

2
DO

1.5

0.5

0
-2 -1 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
t

De la etapa del proceso donde se corta el suministro de aire se puede calcular la

pendiente, para obtener el valor así obtener el valor de la demanda de oxigeno:

𝑑𝐶𝐿
= −𝑞𝑂2𝑋
𝑑𝑡
 Aire off
3.5 y = -0.6971x + 2.9762
R² = 0.9132
3
2.5
2
1.5
DO

1
0.5
0
0 1 2 3 4 5 6
-0.5
-1
t

El valor de la pendiente y por lo tanto, de la demanda de oxigeno es qO2X = 0.6971

m/L.min, este dato se utiliza para graficar la densidad de oxigeno VS la transferencia total
de oxígeno y así obtener el valor de KLa:

Aire on
3.5

2.5
y = 2.2598x - 2.0395
2 R² = 0.9956
OD

1.5

0.5

0
0 0.5 1 1.5 2 2.5
dCl/dt + qO2X

1 𝑑𝐶𝐿
𝐶𝐿 = − ( ⁄𝑑𝑡 + 𝑞𝑂2𝑋) + 𝐶 ∗
𝐾𝐿𝑎
1 −1
− 𝐾𝐿𝑎 = −9.219 ∴ 𝐾𝐿𝑎 = −9.219
→ 𝐾𝐿𝑎 = 0.1085 𝑚𝑖𝑛−1
Solución

Datos Formulas

V (L) 800 La ecuación global de transferencia en el steady state:

𝐾𝐿 𝑎 (𝐶 ∗ − 𝐶𝐿 ) = 𝑞𝑂2 𝑥

x (g/L) 20 𝑞𝑂2 𝑚𝑎𝑥 𝐶𝐿

𝑞𝑂2 = 𝑚𝑔
0.2 𝐿 + 𝐶𝐿

KLa (h-1) 120 Igualando las ec. anteriores:

% O2 entrada 80 𝐾𝐿 𝑎 (𝐶 ∗ − 𝐶𝐿 ) 𝑞𝑂2 𝑚𝑎𝑥 𝐶𝐿

= 𝑚𝑔
𝑥 0.2 + 𝐶𝐿
𝐿

C* (mg/L) 28 Igualando a 0:

𝑞𝑂2𝑚𝑎𝑥 (mg. O2/g*h) 240 −𝐾𝐿 𝑎𝐶𝐿2 − 0.2𝐾𝐿 𝑎 + 𝐾𝐿 𝑎𝐶 ∗ 𝐶𝐿 − 𝑞𝑂2𝑚𝑎𝑥 𝐶𝐿 𝑥 + 0.2𝐾𝐿 𝑎𝐶 ∗ = 0

Resultado Sustitución
𝑚𝑔 −(120 ∗ 𝐶𝐿2 ) − (0.2 ∗ 120 ∗ 𝐶𝐿 ) + (120 ∗ 28 ∗ 𝐶𝐿 ) − (240 ∗ 20 ∗ 𝐶𝐿 )
𝐶𝐿 = 0.443
𝑙 + (0.2 ∗ 120 ∗ 28) = 0

−120𝐶𝐿2 − 1464𝐶𝐿 + 672 = 0

Por medio de una resolución ecuación cuadrática se encuentra

𝐶𝐿
 1. Durante la fermentación de Pseudomonas ovalis B 1486 a distintas velocidades de
agitación se obtuvieron los datos mostrados abajo. El volumen del fermentador fue
de 3 litros. El flujo de aire en el recipiente se mantuvo a 1 vvm. La presión del aire
fue de 3 atm y la temperatura de 29 °C.

Tiempo de fermentación (h) 300 rpm 500 rpm

CL F0 CL F0

0 5.9 0.21 5.9 0.21

4 5.6 0.209

5 5.3 0.209

6 5.2 0.208

7 4.7 0.208 4.9 0.207

8 4.1 0. 207 4.4 0.206

9 3.4 0.206 4.0 0.205

10 3.4 0.206 4.0 0.205

11 3.5 0.207 4.2 0.205

Determine Kla para cada velocidad de agitación.

Solución

Datos: Formulas

V = 3L ρa (Q)
𝐶𝐿 = − (𝑓𝑖 − 𝑓0) + 𝐶 ∗
Q = 1 vvm 𝐾𝐿𝑎 (𝑉)

ρa (3 atm, 29°C) = 3.5047 kg/m3 𝑄 = 𝑣 (𝑉)

60 𝑚𝑖𝑛
𝑄 = (1 𝑣𝑣𝑚)(3 𝐿) × ( ) → 𝑄 = 180 𝐿/ℎ
1ℎ
kg 𝑔 1000 𝑚𝑔 1 𝑚3 𝑚𝑔
ρa = 3.5047 = 3.5047 𝐸03 × × = 3504.7
m3 𝑚3 1𝑔 1000 𝐿 𝐿

ρa = 3504.7 𝑝𝑝𝑚
 300 rpm 500 rpm
Tiempo de fermentación
(h) Cl (ppm) f0 fi-f0 Cl (ppm) f0 fi-f0

0 5.9 0.21 0 5.9 0.21 0

4 5.6 0.209 0.001

5 5.3 0.209 0.001

6 5.2 0.208 0.002

7 4.7 0.208 0.002 4.9 0.207 0.003

8 4.1 0.207 0.003 4.4 0.206 0.004

9 3.4 0.206 0.004 4 0.205 0.005

10 3.4 0.206 0.004 4 0.205 0.005

11 3.5 0.207 0.003 4.2 0.205 0.005

300 rpm
7
6
5
4
Cl

3 y = -625x + 5.925
R² = 0.999
2
1
0
0 0.0005 0.001 0.0015 0.002 0.0025 0.003 0.0035 0.004 0.0045
fi-f0
 500 rpm
7
6
5
4
Cl

3 y = -386.88x + 5.9625
R² = 0.9953
2
1
0
0 0.001 0.002 0.003 0.004 0.005 0.006
fi-f0

300 rpm 500 rpm

𝑦 = −625𝑥 + 5.925 𝑦 = −386.88𝑥 + 5.9625

𝑟 = 0.999 𝑟 = 0.995

ρa (Q) (3504.7 ppm)(180 𝐿/ℎ)

𝐾𝐿𝑎 = − 𝐾𝐿𝑎 = −
−𝑚 (𝑉) −(386.83 𝑝𝑝𝑚)(3 𝐿)

(3504.7 ppm)(180 𝐿/ℎ) 𝐾𝐿𝑎 = 543.5326 ℎ−1

𝐾𝐿𝑎 = −
−(625 𝑝𝑝𝑚)(3 𝐿)

𝐾𝐿𝑎 = 336.4512 ℎ−1

 2. Calcule los requerimientos específicos de potencia (P/V) para una configuración
standard, con bafles y equipada con turbinas Rushton conteniendo agua a 250°C.
El diámetro del tanque es de 0.5 m, la velocidad de rotación del agitador es 300 rpm.

Solución

Datos de tablas Formulas Sustitución NRe

(H2O a 25°C)

ρ (kg/m3) 996.95 𝑛𝐷𝑖2 𝜌 (5)(0.5)2 (996.95) 1380.0526

𝑁𝑅𝑒 = 𝑁𝑅𝑒 =
𝜇 0.903

μ 0.903 𝑃𝑔 Valor es obtenido de la curva Np

𝑁𝑝 =
(N*s/m2) 𝑛3 𝐷𝑖5 𝜌 de potencia (al ser turbinas
Rushton, se habla de una
turbina de pala plana)
Datos 5.182
operacionales del
fermentador

Dt (m) 0.5 𝑁𝑝 𝑛3 𝐷𝑖5 𝜌 (5.182)(5)3 (0.5)5 (996.95) P(

𝑚∗𝑘𝑔
)
𝑃= 𝑃= 𝑠
𝑔 9.81
n (s-1) 5 2057.13

g (m/s2) 9.81 P (hp)

27.035

kW

20.16

Nota: Datos de ρ y μ obtenidos interpolando el valor de 25°C para agua en tablas de

Potter, MC., y Scott, E.P. (2006). Termodinámica. México: Thomson Editores, S.A. de
C.V.
 3. Se tiene un fermentador de 30 m3, equipado con 1 set de turbinas de paletas planas
y 4 deflectores. La geometría del fermentador es:

HL = DT HL: Altura del líquido en el tanque

Di = 1/3* DT DT: Diámetro del tanque

Di: Diámetro del agitador

Las características del caldo de cultivo son una densidad de 1200 kg/m3 y una
viscosidad de 0.02 kg/m sec. La velocidad de aireación es de 0.16 vvm. Si la
potencia entregada al sistema aireado (Pg) es de 80.84 kgf m/s y se supone que el
fermentador opera bajo régimen turbulento.

Determine la velocidad máxima de agitación a la cual se puede encontrar sometido

el fermentador.

Solución

𝑘𝑔𝑚 𝑚
𝑃 = 𝑃𝑔(𝑔𝑐) = 80.84 (9.81 ) = 792.77 𝑊
𝑠 𝑠2

3 30 (4)
𝐷𝑡 = √ = 3.3678 𝑚
𝜋

3.3678 𝑚
𝐷𝑖 = = 1.1226 𝑚
3

Si se supone que el número de potencia tiene un valor de 6, al tratarse de un flujo turbulento, la

velocidad de rotación se calcularía

𝑘𝑔𝑚
3 𝑃 3 (792.7655 𝑠 ) 60 𝑠𝑒𝑔
𝑁𝑝 = √ = √ × = 23.712 𝑟𝑝𝑚
𝜌(𝑁𝑝)(𝐷𝑖)5 𝑘𝑔 5 1 𝑚𝑖𝑛
(1200 𝑚3)(6)(1.1226 𝑚)
 4. Las dimensiones de un fermentador equipado con un agitador de turbina y 4
deflectores, son las siguientes:

Diámetro del agitador = 1,5 m

Diámetro del fermentador = 4,5 m

Anchura de cada tabique deflector = 0,45 m

Altura del líquido = 5 m.

La densidad del caldo de cultivo es de 1200 kg/m3 y su viscosidad de 0,02 kg/ms.

La velocidad de giro del agitador es de 60 r.p.m. mientras que la velocidad de
aireación es de 0,4 vvm.

Calcular:

a) Potencia requerida por el agitador en ausencia de aireación.

b) Potencia requerida cuando el sistema se encuentra aireado.

Solución

Datos Formulas

gc = 9.81 m/s2 ρ = 1200 kg/m3 (𝑁𝑖)(𝐷𝑖 2 )(ρ)

𝑅𝑒 =
μ
Dt = 4.5 m μ = 0.02 kg/ms
(𝑃)(𝑔𝑐)
Di = 1.5 m Ni = 60 rpm 𝑁𝑝 =
(𝑁𝑖 3 (𝐷𝑖 5 )(ρ)
Wb = 0.45 m vel = 0.4 vvm 𝑃𝑟𝑒𝑎𝑙 = 𝑃(𝑓𝑐)
Hl = 5 m ∗
(𝐷𝑙⁄𝐷𝑖 ) (𝐻𝑙⁄𝐷𝑖 )∗
𝑓𝑐 = √
(𝐷𝑙⁄𝐷𝑖 )(𝐻𝑙⁄𝐷𝑖 )

𝑃𝑔
= 𝑓(𝑁𝑎) ; 𝐹 𝑣𝑒𝑙(𝑉)
𝑃
𝐹
𝑁𝑎 =
(𝑁𝑖)(𝐷𝑖 3 )
a)

𝑘𝑔 1
(1200 𝑚3)(𝑠𝑒𝑔)2 (1.5 𝑚)
𝑅𝑒 = = 90,000 → 𝑇𝑢𝑟𝑏𝑢𝑙𝑒𝑛𝑡𝑜 ∴ 𝑁𝑝 = 6
𝑘𝑔
(0.02 𝑚𝑠𝑒𝑔)

𝑘𝑔 1 3
(6)(1200 )( ) (1.5 𝑚)5
𝑚3 𝑠𝑒𝑔
𝑃= 𝑚 = 5573.3945 𝑘𝑔. 𝑚/𝑠𝑒𝑔
9.81 2
𝑠𝑒𝑔

4.5 5
( )( )
𝑓𝑐 = √ 1.5 1.5 = 1.0541 ∴ 𝑃𝑟𝑒𝑎𝑙 = (5573.3945 𝑘𝑔𝑚) (1.0541) = 5874.9151 𝑘𝑔𝑚/𝑠𝑒𝑔
(3)(3) 𝑠𝑒𝑔

b)

𝑑 2 3𝑚 2
𝑉 = 𝜋( ) ℎ = 𝜋 ( ) (5 𝑚) = 35.3429 𝑚3
2 2
1 𝑚𝑖𝑛
𝐹 = 𝑣(𝑉) = (0.4 𝑣𝑣𝑚)(35.3429 𝑚3) × = 0.2356 𝑚3/𝑠𝑒𝑔
60 𝑠𝑒𝑔

0.2356 𝑚3/𝑠𝑒𝑔 𝑃𝑔
𝑁𝑎 = = 0.0698 ∴ = 0.68
1 𝑃0
( )(1.5 𝑚)
𝑠𝑒𝑔

𝑘𝑔𝑚 𝑘𝑔𝑚
𝑃𝑔 = 𝑃0(0.68) = 5874.9151 (0.68) = 3994.9423
𝑠𝑒𝑔 𝑠𝑒𝑔
 5. Un tanque cilíndrico de 1.22 metros de diámetro se llena con agua para operar a un
nivel igual al diámetro del tanque. El tanque se encuentra equipado con cuatro
dispersores regularmente espaciados cuya amplitud es un décimo del diámetro del
tanque. El tanque es agitado con seis turbinas de hoja plana con diámetro de 0.30
m. la velocidad de rotación del impulsor es de 2.8 rps. El aire entra a través de una
tubería situada abajo del impulsor y tiene un flujo de 0.00416 m3 /s a 1.08 atm y 25
C.

a) calcule los requerimientos de potencia para el sistema no aireado

b) para el sistema aireado

La viscosidad y densidad del agua del agua a 25 C son:

μ = 8.904x 10-4 Kg/m s 𝜌 = 997.08 Kg/m3

Solución