Predicción del comportamiento hidrodinámico en el escalado de un reactor de tanque agitado para procesos aerobios, mediante CFD
Ricardo Gelves1 Aldo Benavides2 Juan Carlos Quintero3
1Grupo de Biotecnología. Sede de Investigación Universitaria-SIU. Universidad de Antioquia. Carrera 53 No 61-30. Medellín, Colombia. E-mail: egrgz832@udea.edu.co
2Grupo de Bioinstrumentación e Ingeniería Clínica. Universidad de Antioquia. Calle 67 N. 53-108. Medellín, Colombia. E-mail: aldogbmoran@gmail.com
3Grupo de Bioprocesos. Universidad de Antioquia. Calle 67 N. 53-108. Medellín, Colombia. E-mail: jcquinte@udea.edu.co
RESUMEN
En el presente trabajo se analiza mediante simulaciones en CFD (Dinámica de Fluidos Computacional) el comportamiento de la hidrodinámica líquido-gaseosa en el escalado de un biorreactor, aplicando criterios empíricos comúnmente utilizados en procesos biotecnológicos (Número de Reynolds Re y Potencia por unidad de volumen P/V) en reactores de tanque agitado de 0,050 y 0,500 m3. Los efectos de la turbulencia, flujo rotativo, ruptura y coalescencia de burbujas son simulados mediante los modelos k-?, MRF (Multiple Reference Frame) y PBM (Population Balance Model). El criterio de escalado P/V presenta mejores similitudes hidrodinámicas que el criterio de escalado Re respecto del reactor de 0,050 m3; sin embargo, el alto consumo de potencia y el posible daño microbiano generado por la turbulencia en cultivos de células animales pueden ser desventajas en prototipos a gran escala.
Palabras clave: Biorreactor, dinámica de fluidos computacional (CFD), marco de referencia múltiple (MRF), método de balance poblacional (PBM), escalado.
ABSTRACT
In this paper the effects of some empirical scale up criteria (Re y P/V) on hydrodynamics are analyzed by using CFD simulations (Computational Fluid Dynamics) in 0.050 and 0.500 m3 stirred tank bioreactors. The effects of turbulence, rotating flow, bubbles breakage and coalescence were simulated by using the k-?, MRF (Multiple Reference Frame) and PBM (Population Balance Model), respectively. The P/V scale up criteria showed better similar hydrodynamics than Re scale up criteria with respect to the 0.050 m3reactor. However, the high power consumption and cellular damage due to turbulence in animal cell cultures could be disadvantages in large scale prototypes of bioreactors.
Keywords: Bioreactor, computational fluid dynamics (CFD), multiple reference frame (MRF), population balance model (PBM), scale up.
INTRODUCCIÓN
El éxito de muchos procesos biotecnológicos a escala industrial depende de una agitación y aireación efectiva de los fluidos presentes. La hidrodinámica desarrollada en un biorreactor promueve la transferencia de masa y calor entre las fases y la distribución uniforme de los nutrientes durante un proceso fermentativo. Durante el cambio de escala surgen problemas asociados con la hidrodinámica de biorreactores que consecuentemente reducen, en algunos casos, la productividad de un bioproceso, respecto de escalas de laboratorio [1-3]. El alejamiento de las condiciones de mezcla completa con el incremento de la escala, genera gradientes que conducen a un distanciamiento de las condiciones óptimas encontradas en condiciones de laboratorio. Gradientes de gases, esfuerzos cortantes, entre otros, en fluidos multifásicos que comúnmente son utilizados en procesos biotecnológicos, son causados por una deficiente mezcla generada por la adopción de métodos empíricos como estrategias de escalado. Todos los criterios empíricos de escalado como Re (Número de Reynolds constante) y P/V (Potencia por unidad de volumen constante), entre otros, se basan en aproximaciones ideales para mantener un parámetro de operación constante en ambas escalas, dejando de lado la hidrodinámica y los fenómenos de transporte que gobiernan el proceso de escalado. Conocer el comportamiento hidrodinámico durante la operación de biorreactores a diferentes escalas de operación permite identificar el grado de alejamiento de la mezcla perfecta que viene asociado cuando se emplean los diferentes criterios empíricos en el escalado. La dinámica de fluidos computacional (CFD del inglés Computational Fluid Dynamics) es una muy útil herramienta para analizar en detalle el fenómeno de la turbulencia en biorreactores de tanque agitado, los cuales son comúnmente usados en procesos de fermentaciones microbianas. Para la simulación en CFD del escalado de un proceso aerobio es importante considerar fenómenos como la ruptura y coalescencia de burbujas a partir de modelos de balance poblacional, ya que en el biorreactor estos fenómenos desencadenan la formación de zonas poco oxigenadas y gradientes de gases que pueden afectar la productividad durante el escalado debido a que: i) El comportamiento de las burbujas tiene una influencia significativa en la transición del régimen de flujo porque el fenómeno más importante que causa esta transición desde un régimen homogéneo a un régimen heterogéneo es la ocurrencia de grandes burbujas, ii) El tamaño de las burbujas tiene una amplia distribución en un régimen heterogéneo que comúnmente se presenta en un proceso de cambio de escala [4], iii) Las velocidades radiales, axiales y tangenciales de las burbujas y la disipación de la energía cinética turbulenta influyen significativamente en los fenómenos de ruptura y coalescencia que como consecuencia definen los perfiles de gradientes de aire en un biorreactor. En los últimos años ha habido un rápido avance en el modelado a través de CFD (Dinámica de Fluidos Computacional) de sistemas de agitación de reactores, incluyendo desde modelos de rotación de fluidos en una fase [5-9] hasta la utilización de modelos de balance poblacional en sistemas de agitación líquido-gaseosa [10-27]. Sin embargo, no existen en la actualidad investigaciones que analicen mediante predicciones en CFD (Dinámica de Fluidos Computacional), los efectos que podrían generarse en fenómenos como la ruptura y coalescencia de burbujas, como consecuencia de comportamientos no ideales en la hidrodinámica generados por la aplicación de criterios empíricos comúnmente utilizados durante la etapa de escalado de biorreactores de tanque agitado. Por tal motivo en el presente trabajo se pretende identificar comportamientos no ideales generados durante estudios de escalado entre 0,050 y 0,500 m3 empleando dos criterios empíricos de escalado utilizados comúnmente en bioprocesos. Para esto se realizaron simulaciones con el software de CFD (Dinámica de Fluidos Computacional) ANSYS-FLUENT v13 que permite obtener información sobre parámetros hidrodinámicos como perfiles de dispersión de gases, turbulencia, distribución de tamaños de burbujas, etc. Esta estrategia de predicción de problemas que se pueden generar durante la etapa de escalado de un bioproceso puede evitar la necesidad de construir prototipos costosos e ineficientes, lo que traería como ventaja un ahorro en tiempo de experimentación y presupuesto para una empresa del sector biotecnológico.
MATERIALES Y MÉTODOS
Características físicas del reactor de 0,050 m3
Para la etapa inicial (volumen de referencia) se utilizó un biorreactor de tanque agitado New Brunswick BioFlo 5000 a escala piloto (0,080 m3 de volumen total) con 0,050 m3 de agua como volumen de trabajo (temperatura de 293 K y presión de 101325 Pa). El tanque está caracterizado por las siguientes dimensiones: diámetro DT: 0,380 m, altura del líquido HL: 0,380 m, equipado con cuatro deflectores espaciados en un ángulo de 90° con una anchura Wb: 0,1 Dt instalados a 0,010 m de la pared del tanque. El equipo cuenta con un impulsor tipo turbina Rushton (disco con seis paletas rectas y perpendiculares) con un diámetro Di: 0,1265 (Di: 0,3 DT), ancho de la paleta Wi: 0,030 m (Wi: 0,25 Di), altura de la paleta WL: 0,040 m (WL: 0,33 Di:). El aire es suministrado a través de un difusor de 0,125 m de diámetro localizado debajo de la turbina a una altura de 0,09 m. Para la etapa de escalado virtual simulada en CFD se dimensionó un biorreactor de 0,500 m3 (DT: 0,75 m y HL: 1,1 m)con una similitud geométrica respecto de biorreactor de 0,050 m3. Las propiedades utilizadas para la fase primaria (agua) fueron: ?L: 998,2 kgm-3; ?L: 0,001 kgm-1 s-1. ?: 0,073 Nm1. Las Propiedades de la fase secundaria (aire) fueron: ?G: 1,225 kgm-3, ?G = 1,789 10-5 kgm-1 s-1[28].
Criterios empíricos de escalado
Las condiciones de operación para el proceso de escalado a 0,500 m3 se calcularon con base en los criterios Re y P/V teniendo en cuenta la velocidad de agitación Ni del biorreactor de 0,050 m3 (400 rpm) y manteniendo constante el flujo de aire a 1 vvm en ambas escalas.
Criterio número de Reynolds constante Re. El escalado con base en el número de Reynolds asume que la magnitud relativa de los esfuerzos inerciales respecto de los esfuerzos viscosos se mantiene constante en ambas escalas. Si los esfuerzos inerciales dominan respecto de los esfuerzos viscosos se dice que el flujo no puede ser laminar. A la pérdida de laminaridad se le denomina turbulencia y esta es muy importante para la transferencia de masa y dispersión de gases. En esta investigación se utilizó el criterio Número de Reynolds constante, ya que por su carácter teórico y empírico define una supuesta semejanza en el patrón de flujo hidrodinámico, lo que significa que adoptando este criterio se trataría de mantener el mismo grado de turbulencia, en ambas escalas [29]. Con el incremento de escala, la velocidad de agitación requerida para mantener Re constante se reduce significativamente y el grado de heterogeneidad puede aumentar, por lo que se generarían comportamientos no ideales que podrían ser identificados a través de Dinámica de Fluidos Computacional.
Criterio potencia por unidad de volumen constante P/V. El criterio de escalado P/V se basa en mantener constante el consumo de potencia por unidad de volumen de un proceso líquido-gaseoso debido a la agitación [29]. P/V requiere de un mayor grado de agitación respecto de la velocidad requerida para mantener el criterio Re constante. Es por tal motivo que su desventaja radica en el alto consumo de potencia, ya que la velocidad en la punta del agitador
se incrementa con la escala y, por lo tanto, el diseño de un motor de agitación a gran escala resulta poco práctico en los bioprocesos. Como consecuencia del incremento en la velocidad en la punta del agitador, el esfuerzo cortante puede ocasionar daños en la morfología de células microbianas. Sin embargo, una agitación y aireación efectiva de los fluidos presentes en un biorreactor promueven la dispersión de gases y tamaño de burbujas para la transferencia de masa. Lo anterior soporta las grandes diferencias significativas en la hidrodinámica líquido-gaseosa que podrían presentarse utilizando los criterios P/V y Re. Por esta razón en esta investigación se han incluido estos criterios de escalado, ya que se busca predecir a través de CFD tendencias hidrodinámicas opuestas generadas según el criterio de escalado empleado (P/V y Re).
Modelo de CF
Ecuaciones de Navier-Stokes
La fase gaseosa y la fase líquida son tratadas como series continuas de interpenetración y las ecuaciones de conservación de masa y momentum se resuelven para cada fase. Debido a que el volumen de una fase no puede ser ocupado por las otras fases se utiliza el concepto de fracción de volumen que son funciones continuas en el espacio y el tiempo. Las ecuaciones de conservación para cada fase se derivan para obtener un sistema de ecuaciones, que tienen estructura similar para todas las fases [17-22].
El modelo euleriano es el más complejo de los modelos multifásicos en ANSYS FLUENT v13. Soluciona un sistema de n-ecuaciones de momentum y continuidad para cada fase presente. El acoplamiento se alcanza a través de la presión y de los coeficientes de intercambio interfacial. La ecuación de conservación de masa para cada fase se muestra a continuación:
 |
(1) |
Donde 
representan la densidad, la fracción de volumen y la velocidad de la fase i (G o L). Se asume que la fase líquida L y la fase gaseosa G comparten el espacio en proporción a su volumen, de manera que el total de las fracciones suman un valor igual una unidad.
 |
(2) |
La ecuación de momentum para la fase i se describe a continuación:
 |
(3) |
? es la presión compartida por ambas fases, 
es la gravedad y 
representa el intercambio de momentum interfaz. El término 
representa las fuerzas de Coriolis y centrífuga expresada en el modelo para flujos rotativos MRF y se representa como:
 |
(4) |
es la velocidad angular, 
es la posición vectorial. El tensor de esfuerzos de Reynolds 
se relaciona con los gradientes medios de velocidad mediante la hipótesis de Boussinesq [17]:
 |
(5) |
?lam, i es la viscosidad molecular de la fase i, 
es el tensor de deformación, ?t, i es la viscosidad turbulenta de la fase i.ki es la energía cinética turbulenta.
Intercambio de momentum interfacial
La fuerza más importante de intercambio entre las fases es la fuerza de arrastre RG que actúa sobre las burbujas que funciona como una fuerza de cerradura para ambas fases en las ecuaciones de momentum. Esta fuerza depende de la fricción, de la presión, de la cohesión y de otros efectos hidrodinámicos [30].
 |
(6) |
K es el coeficiente de intercambio de las fases líquida y gaseosa y se determina con la ecuación:
 |
(7) |
d es el diámetro de burbuja y Cd es el coeficiente de arrastre que se define como una función de número de Reynolds Rep:
 |
(8) |
Para el cálculo del coeficiente de arrastre se utiliza la correlación estándar [31]:
 |
(9) |
Ecuaciones para el modelo de turbulencia
El modelo de turbulencia k-e disperso es el modelo multifásico de turbulencia estándar. Representa la extensión del modelo monofásico k-? y es aplicable cuando las concentraciones de la fase secundaria son diluidas. Las ecuaciones de k y ? que describen este modelo son como sigue:
En estas ecuaciones, Gk, l representa la generación de energía cinética turbulenta kL de la fase líquida debido a los gradientes de la velocidad media, ?L es la disipación de la energía cinética turbulenta. 
k, L y?, L representan la influencia de la fase dispersa en la fase continua y son modeladas con las ecuaciones de Elgobashi y Rizk [32].
La viscosidad turbulenta ?t, L es calculada de:
 |
(12) |
Los valores de las constantes utilizados en este experimento fueron C1?: 1,44 C2?: 1,92 C?: 0,09 ?k: 1,00 y ??: 1,30. ?k y ?? representan el número de Prandtl turbulento para k y ?, respectivamente.
Método de balance poblacional
En este trabajo se utilizó el método discreto [33-35] para solucionar las ecuaciones de balance poblacional por su robustez numérica, ya que la población de burbujas es discretizada en un número finito de intervalos de diámetros. Las ecuaciones de balance poblacional para las diferentes clases de burbujas en el reactor se pueden escribir como [27, 36-37]:
 |
(13) |
Donde ni es el número de clases de burbujas i, BiC y BiB son las velocidades de nacimiento debidas a la coalescencia y rotura, respectivamente, Dic y DiB son las tasas de mortalidad.
Los términos de rotura y coalescencia son:
?(v, v') es la tasa de coalescencia entre las burbujas de tamaño v y v'; g(v) es la tasa de ruptura de burbujas de tamaño v; g(v') es la frecuencia de ruptura de burbujas v' y ?(v|v') es la probabilidad de la función densidad de las burbujas rotas desde el volumen v' a una burbuja de volumen v; p es el número de burbujas formadas por cada rompimiento; n es el número de clases de burbujas y ?v es el volumen de una burbuja v. La ruptura de burbujas es analizada en términos de la interacción de las burbujas con los remolinos turbulentos. Estos remolinos turbulentos incrementan la energía de superficie de las burbujas hasta causar la deformación. La ruptura ocurre si el incremento en la energía de superficie alcanza un valor crítico. La tasa de ruptura se define como [6]:
Donde d es el diámetro de partícula, ? es el tamaño adimensional de remolinos, f es la frecuencia de rompimiento, ? es la tensión superficial y ? es la fracción de partículas. La coalescencia de burbujas se modela considerando la colisión de burbujas debida a la turbulencia, flotabilidad y corte laminar. La tasa de coalescencia se define como el producto entre la frecuencia de colisión ?ag(vi, vj) y la probabilidad de coalescencia Pag(vi, vj) y se define como [6]:
 |
(21) |
La frecuencia de colisión se define como:
 |
(22) |
Donde 
es la velocidad característica de colisión de dos partículas con diámetro di y dj y densidad de burbujas ni y ni
 |
(23) |
Donde c1 es una constante de primer orden, xij = di/dJ ?1 y ?2 son las densidades de las fases primaria (agua) y secundaria (aire).
El número de clases de burbujas i, ni se relaciona con la fracción de volumen de la fase gaseosa, de la siguiente manera:
 |
(24) |
La sumatoria de las fracciones en volumen de cada grupo de burbujas es igual a la fracción en volumen de la fase dispersa:
 |
(25) |
La fracción en volumen de cada grupo de tamaños se expresa en términos de la fracción total de la fase dispersa:
 |
(26) |
Por lo tanto, incluyendo la ecuación (26) en la ecuación (13), resulta:
 |
(27) |
Para acoplar las ecuaciones de balance poblacional con la hidrodinámica del reactor se utiliza el diámetro de sauter d32 como valor de entrada en la ecuación:
 |
(28) |
Especificaciones de los reactores y técnica numérica
El dominio de solución para esta investigación se presenta en la Figura 1. Se utilizó el software AnsysFluent 13 [38] para generar la geometría y mallado de los reactores. La malla híbrida en 3D es transformada a celdas poliédricas y está constituida por 450,095 celdas para el reactor de 0,05 m3 y 700,000 celdas para el reactor de 0,500 m3. La calidad del mallado se analiza mediante el criterio skewness [38]. En esta investigación todas las celdas presentan un valor skewness inferior a 0,70, lo que indica que el mallado es aceptable [38]. Se utilizó la técnica de volúmenes finitos [30] para convertir las ecuaciones de Navier-Stokes en ecuaciones algebraicas que pueden ser resueltas numéricamente. Las paredes del tanque, las superficies del agitador y los bafles en los biorreactores se tratan como condiciones de no deslizamiento y funciones de pared estándar.
Figura 1. Mallado en 3D del biorreactor de 0,05 m3 con 450,095 celdas.
En la superficie líquida de los biorreactores se adicionó una pequeña zona de aire (Pressure outlet) para evitar que el líquido se escape del tanque [11], por lo tanto, solo se permite que se escape el gas. El flujo de gas en el difusor se define como una condición de entrada (velocity inlet) con una fracción de volumen de gas igual a una unidad.
Para solucionar las ecuaciones diferenciales parciales se utilizó el algoritmo PC SIMPLE que acopla la presión y velocidad. También se utilizó una discretización implícita de primer orden para la integración del término transitorio y un esquema Upwind de segundo orden para los términos espaciales. Se asume que la solución converge cuando los residuales escalados permanecen con valores menores a 10-5.
Técnicas experimentales y correlaciones utilizadas
Para verificar la validez de las simulaciones se determinó experimentalmente el valor de kL a escala de 0,05 m3 y se comparó con valores obtenidos por CFD. Adicionalmente se utilizaron correlaciones de la bibliografía para verificar las simulaciones a escala de 0,5 m3. Para la determinación experimental de kLa, el biorreactor de 0,05 m3 es gasificado con nitrógeno (N2) hasta alcanzar una concentración mínima de oxígeno disuelto. En esta etapa se adiciona aire al biorreactor a través del difusor y se toman mediciones temporales de oxígeno disuelto CL hasta que se alcance la saturación. La concentración de oxígeno disuelto se mide utilizando un sensor polarográfico de METTLER TOLEDO. La velocidad de cambio de oxígeno en la fase líquida es determinada usando la siguiente ecuación:
 |
(29) |
Donde C* es la concentración de saturación de oxígeno disuelto. Integrando la ecuación (29) es posible obtener la siguiente ecuación:
 |
(30) |
Graficando el término izquierdo de la ecuación (30) en función del tiempo se obtiene el valor experimental de kL?.
Para calcular la velocidad en la punta del impulsor NDa se utilizó la siguiente correlación:
 |
(31) |
Para calcular la potencia consumida de una turbina Rushton con aireación, Nagata [43] determinó la siguiente correlación a partir de datos experimentales:
 |
(32) |
Donde Pg es la potencia consumida con aireación, P potencia consumida sin aireación, v velocidad del gas y Q flujo volumétrico del gas. Las demás variables ya han sido definidas anteriormente.
Se calcularon los valores del torque M del impeler para determinar la potencia consumida mediante CFD, para compararlos con los valores obtenidos por las correlaciones antes mencionadas y verificar la validez de las simulaciones utilizando la siguiente expresión matemática:
 |
(33) |
La teoría de la penetración de Hibies fue usada para calcular los valores de kL? a partir de CFD:
 |
(34) |
 |
(35) |
kL es el coeficiente de transferencia de la fase líquida, DCL es la difusividad del oxígeno, a es el área interfacial de las burbujas, d32 es el diámetro de sauter y ?L es la viscosidad molecular de la fase líquida.
RESULTADOS Y DISCUSIONES
El objetivo principal de esta investigación fue predecir el comportamiento hidrodinámico en el escalado de un proceso aerobio utilizando CFD (Dinámica de Fluidos Computacional). Se hizo un énfasis especial en la hidrodinámica y la transferencia de masa líquido-gaseosa. Las condiciones de operación y criterios de escalado simulados en este trabajo se basan en valores comúnmente usados en cultivos de microorganismos aerobios. En la Figura 2 se observan los contornos de distribución de la fase dispersa (aire) en los biorreactores de 0,050 y 0,500 m3 respectivamente, calculados mediante CFD aplicando los criterios de escalado Re y P/V. Usando el criterio de escalado Re en el biorreactor de 0,500 m3 (Figura 2 b) se observa poca distribución del gas, especialmente en el fondo del biorreactor, en comparación con el biorreactor de 0,050 m3 (Figura 2 a) debido a las bajas fuerzas centrífugas generadas por la turbina Rushton que no son suficientes para superar el fenómeno de flotabilidad del gas (reflejado en los perfiles de velocidad del aire, Figura 4 b) ocasionando de esta forma ambientes heterogéneos (poca dispersión de aire en el fondo del biorreactor) y gradientes que pueden afectar el crecimiento de un microorganismo aerobio. Teniendo en cuenta la transferencia de masa simulada mediante CFD (Figura 6 b) se observa que en estas zonas del biorreactor, la transferencia de oxígeno es muy pobre. Este fenómeno de flotabilidad es más pronunciado con el criterio de escalado Re debido a la baja velocidad de agitación requerida y al incremento en la velocidad superficial del aire para tratar de mantener en ambas escalas el mismo flujo de aire por unidad de volumen del biorreactor (restricción para los dos criterios de escalado).
Figura 2. Fracción en volumen de la fase dispersa (aire) para los biorreactores de (a) 0,050 m3 y 0,500 m3 aplicando los criterios de escalado (b) Re y (c) P/V.
Figura 3.1. Diámetro medio de sauter d32 (m) de la fase dispersa (aire) del biorreactor 0,050 m3. Corte transversal a 0,12 m de altura.
Figura 3.2. Diámetro medio de sauter d32 (m) de la fase dispersa (aire) para los biorreactores de (a) 0,050 m3 y 0,500 m3 aplicando los criterios de escalado (b) Re y (c) P/V.
Figura 4. Perfiles de velocidad (m/s) de la fase dispersa (aire) para los biorreactores de (a) 0,050 m3 y 0,500 m3 aplicando los criterios de escalado (b) Re y (c) P/V.
Figura 5. Perfiles del diámetro de Kolmogórov (m) para los biorreactores de (a) 0,050 m3 y 0,500 m3 aplicando los criterios de escalado (b) Re y (c) P/V.
Figura 6. Perfiles del coeficiente de transferencia de masa kLa(s-1)para los biorreactores de (a) 0,050 m3 y 0,500 m3 aplicando los criterios de escalado (b) Re y (c) P/V.
Las grandes diferencias en cuanto a los perfiles (34) de aire utilizando el criterio P/V (Figura 2 c) se deben a que la velocidad requerida para mantener este criterio constante en ambas escalas es mayor respecto del criterio Re; por lo tanto, las fuerzas centrífugas se incrementan, y el aire es capturado detrás de las aspas del agitador, fenómeno que disminuye el escape de las burbujas del biorreactor, lo que genera una mejor dispersión de las burbujas de aire y un ambiente hidrodinámico gobernado por fenómenos de ruptura de burbujas que trae como consecuencia una mejor transferencia de oxígeno (valores altos de kL?).
La dispersión del gas se produce principalmente en las zonas cercanas a la turbina. Este es arrastrado a las zonas de menor presión estática existentes detrás de las aspas del agitador. El gas que sale del difusor junto con una gran parte de gas recirculado es arrastrado a estas cavidades.
En una fermentación aeróbica, el aire fluye de forma continua a través del reactor de tanque agitado. Es por tal razón que la distribución de tamaño de burbujas resulta de la interacción de mecanismos físicos como la ruptura y coalescencia. El parámetro más importante para analizar estos mecanismos hidrodinámicos es el diámetro medio de sauter. En las Figuras 3.1 y 3.2 se presentan los contornos del diámetro medio de sauter en los biorreactores de 0,050 (a) y 0,500 m3, aplicando los criterios de escalado (b) Re y (c) P/V. Detrás de las aspas se acumulan las burbujas para luego ser dispersadas (Figura 3.1), lo que confirma el fenómeno de coalescencia en estas aéreas. Este comportamiento también ha sido reportado por Kerdouss [11].
Como se explicó anteriormente, el fenómeno de flotabilidad encontrado en el biorreactor de 0,500 m3, aplicando el criterio de escalado Re constante (Figura 3.2 b), influye negativamente en la dispersión de burbujas por lo que el aire permanece menor tiempo en el biorreactor y por esta razón la interacción entre las burbujas debida a los fenómenos de rompimiento, es débil, lo que genera zonas muertas en el fondo del biorreactor, donde las burbujas de aire no alcanzan a ser dispersadas. Por lo tanto estas zonas tienden a ser poco oxigenadas según los bajos valores de kL? simulados (Figura 6 b). Este fenómeno también ha sido reportado por Wang [4] en columnas de burbujeo a escala de laboratorio.
Con el incremento de escala, el fenómeno de coalescencia puede ser importante en un proceso líquido gaseoso. Tal es el caso que se presenta con el criterio de escalado P/V (Figura 3.2 c) donde el comportamiento no ideal desencadenado por la pérdida de homogeneidad debida al incremento de escala observado en el biorreactor de 0,500 m3 es causado por los fenómenos de ruptura y coalescencia, ya que en las zonas comprendidas entre el fondo del biorreactor, la turbina de agitación y los bafles deflectores, donde la disipación de la energía cinética turbulenta es más alta y los tamaños de remolinos son más bajos (Figura 5 c), el fenómeno de ruptura reduce significativamente el tamaño de las burbujas, causando de esta manera la deformación de las mismas hasta causar el rompimiento, como se observa gradualmente desde el fondo del biorreactor hasta la superficie a excepción del centro de rotación; por lo tanto, existe una mejor transferencia de oxígeno que en las zonas cercanas a la superficie del fluido líquido donde el fenómeno de coalescencia generado por las altas velocidades tangenciales, la baja disipación de la energía cinética turbulenta y el incremento en el tamaño de remolinos ocasionan comportamientos no ideales en el centro de agitación que deben evitarse durante un proceso de escalado mediante la instalación de una segunda turbina en esta zona del biorreactor.
En la Figura 4 se presentan los perfiles de velocidad del gas calculados por CFD en los biorreactores de (a) 0,050 y 0,500 m3, aplicando los criterios de escalado (b) Re y (c) P/V.
Se observan perfiles de velocidad desde 0,1 a 1 m/s para el caso del biorreactor a escala de 0,050 m3 (Figura 4 a) mientras que los perfiles encontrados para el caso del escalado con los criterios empíricos Re y P/V (Figuras 4 b y 4 c) los perfiles de velocidad del aire se reducen observándose intervalos desde 0 a 0,55 y desde 0,2 a 0,8 m/s.
El criterio P/V presenta mayores similitudes con el biorreactor de 0,050 m3 en cuanto a los perfiles de velocidad de aire en el fondo del biorreactor y en las zonas cercanas a la turbina, en su comparación con los perfiles simulados con el criterio Re. Estos altos valores de velocidad encontrados en el fondo del biorreactor conducen a una mejor recirculación de aire y como consecuencia la transferencia de oxígeno kL? se incrementa en estas zonas (Figura 6 c). Sin embargo la pérdida de homogeneidad con el incremento de escala genera comportamientos no ideales en zonas alejadas de la turbina de agitación que se presentan en los dos criterios de escalado (Re y P/V).
Lo anterior, también es reflejado en la simulación del tamaño de remolinos mediante el diámetro de Kolmogórov (Figura 5) que muestra en todos los casos los valores más bajos en las zonas ubicadas entre la turbina y el difusor de aire y también en las paredes de los bafles deflectores. La ruptura de burbujas, como se explicó anteriormente, es analizada en términos de la interacción de las burbujas con los remolinos turbulentos. Estos remolinos turbulentos incrementan la energía de superficie de las burbujas hasta causar la deformación. Es por eso que en las simulaciones se observa que la ruptura ocurre si el incremento en la energía de superficie, generada según el tamaño de remolinos, alcanza un valor crítico (Figuras 3.2 y 5). Por lo tanto, estas zonas son candidatas a una transferencia de masa aceptable como se observa en la Figura 6 para el crecimiento de microorganismos aerobios. Sin embargo, la teoría del daño celular debida a la turbulencia sugiere que si una célula es menor al diámetro de remolinos de Kolmogórov, la entidad biológica no sufrirá daño celular [45]. El tamaño promedio de algunas células animales se encuentra en el orden de 15-18 ?m y teniendo en cuenta los valores mínimos de tamaño de remolinos simulados mediante CFD (4 ?m) aplicando el criterio de escala P/V, se puede afirmar que se espera que ocurra daño microbiano en su estructura celular; por lo tanto, es necesario utilizar otro sistema de agitación cuando el objetivo es el cultivo de células animales.
El kL? es el parámetro que se usa comúnmente para describir la eficiencia de un sistema de aireación. Para calcular los valores de kL? a través de simulaciones en CFD se utilizó la teoría de la penetración de Higbies, en la cual el valor de kL? se obtiene como el producto del coeficiente de transferencia de masa de la fase líquida kL y el área interfaz ? [37]. El área interfaz se expresa como una función de la fracción local de volumen de gas y el diámetro de sauter local. Se observan similitudes entre los valores calculados por CFD a escala de 0,050 m3 y los datos obtenidos experimentalmente (Tabla 1 y Figura 6). Lo que significa que los términos de ruptura y coalescencia del modelo de balance poblacional predicen con buen acercamiento los fenómenos de transferencia de masa involucrados en las condiciones estudiadas a escala de 0,050 m3. Partiendo de esta información se simularon los valores de kL? para el biorreactor de 0,500 m3 con base en los criterios de escalado Re y P/V (Figura 6). Los resultados demuestran el efecto proporcional que ocurre en los valores de kL? con el incremento de la potencia consumida en cada criterio de escalado. Acontecimiento que ya ha sido verificado por investigaciones experimentales [39-41].
Tabla 1. Comparación de valores de CFD (ND?: velocidad en la punta del agitador; Pg: Potencia con aireación y kL?: Coeficiente de transferencia de masa) con datos obtenidos por correlaciones y técnicas experimentales (Exp).
El aire en contacto con las aspas de la turbina tiende a disminuir el coeficiente de retardo o rozamiento asociado con el giro de la turbina y la consecuente disminución del consumo de potencia [42].
En la Tabla 1 se presentan los valores de potencia consumida con aireación Pg y velocidad en la punta del impulsor NDa calculados por CFD y comparados por las correlaciones mencionadas en el apartado de materiales y métodos para los biorreactores de 0,050 y 0,500 m3. Se observan similitudes entre los valores calculados por CFD (Tabla 1) y los datos obtenidos por las correlaciones.
Lo anterior argumenta que el modelo de balance poblacional acoplado a CFD predice con buen acercamiento los fenómenos de turbulencia en las condiciones hidrodinámicas estudiadas. En este trabajo se encontró que la transferencia de masa tiene una muy clara estructura definida por la hidrodinámica y los fenómenos de coalescencia y rompimiento de burbujas. El conocimiento acerca de la transferencia de masa en cultivos microbianos a gran escala es importante para determinar las zonas muertas poco oxigenadas que pueden afectar las actividades metabólicas de un cultivo celular, que se pueden presentar como consecuencia de un proceso de escalado.
Esta estrategia presentada en esta investigación de predicción de problemas que se pueden generar durante la etapa de escalado de un cultivo de microorganismos aerobios mediante CFD, puede ser una herramienta poderosa para comprender y solucionar problemas de transferencia líquido-gaseosa que comúnmente ocurren en cultivos celulares. De esta manera, la información obtenida a partir de simulaciones en CFD de las heterogeneidades en cultivos de microorganismos aerobios podría ser usada para optimizar el escalado de un bioproceso.
CONCLUSIONES
El modelo de CFD aplicado en este trabajo simula los efectos hidrodinámicos generados por los criterios de escalado evaluados (Re y P/V) y predice tendencias hidrodinámicas que afectan un proceso de escalado. La velocidad en la punta del agitador, la potencia consumida con aireación y el coeficiente de transferencia de oxígeno calculados mediante CFD han sido verificados satisfactoriamente en correspondencia con datos experimentales
obtenidos en esta investigación y correlaciones de la bibliografía. El criterio de escalado P/V utilizando una turbina Rushton presenta mejores similitudes hidrodinámicas que el criterio de escalado Re con referencia al biorreactor de 0,050 m3; sin embargo, el alto consumo de potencia y el posible daño microbiano generado por la turbulencia en cultivos de células animales son desventajas en prototipos a gran escala.
REFERENCIAS
[1] F. Bylund, E. Collet, S.O. Enfors and G. Larsson. "Substrate gradient formation in the large-scale bioreactor lowers cell yield and increases by-product formation". Bioprocess and Biosystems Engineering. Vol. 18, Issue 3, pp. 171-180. 1998. DOI: 10.1007/s004490050427.
[2] E. Dunlop. "Micromixing in fermentors: Metabolic changes in Saccharomyces cerevisiae and their relationship to fluid turbulence". Biotechnology and Bioengineering. Vol. 36 N° 8, pp. 854-864. October 1990. DOI: 10.1002/bit.260360816.
[3] C.J. Hewitt, G.N. Caron, B. Axelsson, C.M. Mcfarlane and A.W. Nienow. "Studies related to the scale-up of high-cell-density E. coli fed-batch fermentations using multiparameter flow cytometry: Effect of a changing microenvironment with respect to glucose and dissolved oxygen concentration". Biotechnology and Bioengineering. Vol. 70 N° 4, pp. 381-390. 2000. ISSN: 1097-0290.
[4] T. Wang, J. Wang and Y. Jin. "CFD-PBM Coupled Model for Gas-Liquid Flows". AIChE Journal. Vol. 52 N° 1, pp. 125-140. 2006. DOI: 10.1002/aic.10611.
[5] M. Jenne and M. Reuss. "A critical assessment on the use of k-e turbulence models for simulation of the turbulent liquid flow induced by Rushton turbine in baffled stirred-tank reactors". Chemical Engineering Science. Vol. 54 N° 17, pp. 3921-3941. September 1999. DOI: 10.1016/S0009-2509(99)00093-7.
[6] J.Y. Luo, R.I. Issa and A.D. Gosman. "Prediction of impeller induced flows in mixing vessels using multiple frames of reference". In IchemE Symposium Series. Institution of Chemical Engineers. Vol. 136, pp. 549-556. 1994.
[7] G. Micale, A. Brucato and F. Grisafi. "Prediction of flow fields in a dual-impeller stirred tank". AICHE Journal. Vol. 45 N° 3, pp. 445-464. 1999.
[8] K. Rutherford, K.C. Lee, S.M.S. Mahmoudi and M. Yianneskis. "Hydrodynamic characteristics of dual rushton impeller stirred vessels. AICHE Journal". Vol. 42 N° 2, pp. 332-346. 1996. DOI: 10.1002/aic.690420204.
[9] A.D. Tabor, G. Gosman and R.I. Issa. "Numerical simulation of the flow in a mixing vessel stirred by a Rushton Turbine". IChemE Symposium Series. Vol. 140, pp. 25-34. 1996. DOI:10.1016/j.cep.2005.06.006.
[10] A. Bakker and H.E.A. Van den Akker. "A computational model for the gas-liquid flow in stirred reactors". Transactions of IChemE. Vol. 72 N° A4, pp. 594-606. July, 1994. ISSN: 0888-5885.
[11] F. Kerdouss, A. Bannari, P. Proulx, R. Bannari, M. Skrga and Y. Labrecque."Two-phase mass transfer coefficient prediction in stirred vessel with a CFD model". Computers and Chemical Engineering. Vol. 32 N° 22, pp. 1943-1955. August 2008. DOI:10.1016/j.compchemeng. 2007.10.01.
[12] F. Kerdouss, L. Kiss, P. Proulx, J.F. Bilodeau and C. Dupuis. "Mixing characteristics of an axial flow rotor: Experimental and numerical study". International Journal of Chemical Reactor Engineering. Vol. 3 N° A35. 2005. ISSN: 1542-6580.
[13] F. Kerdouss, A. Bannari, P. Proulx, R. Bannari, M. Skrga and Y. Labrecque. "Two-phase mass transfer coefficient prediction in stirred vessel with a CFD model. Computers and Chemical Engineering". Vol. 3 N° 8, pp. 1-13. 2007. ISSN: 0098-1354.
[14] G.L. Lane, M.P. Schwarz and G.M. Evans. "Numerical modeling of gas-liquid flow in stirred tank. Chemical Engineering Science". Vol. 60 N° 8-9, pp. 2203-2214. 2005. DOI: 10.1016/j.ces.2004.11.046.
[15] C.H. Venneker and H. Derksen Van Den Akker. "Population balance modeling of aerated stirred vessels based on CFD". AICHE Journal. Vol. 48 N° 4, pp. 673-684. 2002. DOI: 10.1002/aic.690480404.
[16] F. Scargiali, A. D'Orazio, F. Grisafi, A. Brucato."Modelling and Simulation of Gas-Liquid Hydrodynamics in Mechanically Stirred Tanks". Chemical Engineering Research and Design, Vol. 85 N° 5, pp. 637-646. 2007. DOI: 10.1205/cherd0624.
[17] G.R. Kasat, A.B. Pandit, y V.V. Ranade. "CFD Simulation of Gas-Liquid Flows in a Reactor Stirred by Dual Rushton Turbines". International Journal of Chemical Reactor Engineering. Vol. 6, pp. 1-28. 2008. ISSN: 1542-6580.
[18] Lu. Wei-Ming. "Multiple Impeller Gas-Liquid Contactors". Informe técnico. National Taiwan University. 2004. ISBN: 957-564-947-7.
[19] V.V. Ranade and J.B. Dommeti. "Computational Snapshot of flow generated by axial impellers in baffled stirred vessels". I Chem. Vol. 74 N° 4, pp. 476-484. 1990. ISSN 0263-8762.
[20] V.V. Ranade. "Computational Flow Modeling for Chemical Reactor Engineering". Process Systems Engineering Series. Academic Press. 2002. ISBN: 0-12-576960-1.
[21] M. Jahoda, L. Tomásková, M. Mostêk. "CFD prediction of liquid homogenization in a gas-liquid stirred tank". Chemical Engineering Research and Design. Vol. 87 N° 4, pp. 460-467. April 2009. DOI: 10.1016/j.cherd.2008.12.006.
[22] E.M. Marshall and A. Bakker, "Computational fluid mixing". Informe técnico. Fluent-Inc. 2002. ISBN: 0-9719532-0-1.
[23] M. Laakkonen, P. Moilanen, V. Alopaeus and J. Aittamaa. "Modelling Local Gas-Liquid Mass Transfer in Agitated Vessels". Chemical Engineering Research and Design. Vol. 85 N° 5, pp. 665-675. 2007. DOI: 10.1205/cherd0617
[24] R. Panneerselvam and S. Savithri. "Investigations on hydrodynamics and mass transfer in gas-liquid stirred reactor using computational fluid dynamics". Chemical Engineering Science. Vol. 66 N° 14-15, pp. 3108-3124. July 2011. DOI: 10.1016/j.ces.2011.03.00.
[25] M. Martinov and S.D. Vlaev. "Increasing Gas-Liquid Mass Transfer in Stirred Power Law Fluids by Using a New Saving Energy Impeller". Chem. Biochem. Eng. Vol. 16 N° 1, pp. 1-6. 2002. ISSN: 0352-9568.
[26] W. Tiefeng and W. Jinfu. "Numerical simulations of gas-liquid mass transfer in bubble columns with a CFD-PBM coupled model". Chemical Engineering Science. Vol. 62 N° 24, pp. 7107-7118. December, 2007.
[27] P. Chen, M.P. Dudukovic and J. Sanyal. "Numerical simulation of bubble columns flows: effect of different breakup and coalescence closures". Chemical Engineering Science. Vol. 60 N° 4, pp. 1085-1101. 2005. ISSN 0009-2509.
[28] CRC Handbook. Handbook of Chemistry and Physics (83rd ed.). Boca Raton, FL, USA: CRC Press LLC. 2002.
[29] A. Junker. "Scale-Up Methodologies for Escherichia coli and Yeast Fermentation Processes". Journal of Bioscience and Bioengineering. Vol. 97 N° 6, pp. 347-364. 2004. ISSN: 1347-4421.
[30] Ansys fluent v13. User's Manual. Fluent Inc. Centrera Resource Park, 10 Cavendish Court. Lebanon, USA. 2011.
[31] M. Ishii and N. Zuber. "Drag coefficient and relative velocity in bubbly, droplet or particulate flows". AIChE Journal. Vol. 25 N° 5, pp. 843-855. June 1979. DOI: 10.1002/aic.690250513.
[32] S.E. Elgobashiand M.A. Rizk. "A two-equation turbulence model for dispersed dilute confined two-phase flows". International Journal of Multiphase Flow. Vol. 15 N° 1, pp. 119-133. February, 1989. DOI: 10.1016/0301-9322(89)90089-X
[33] M.J. Hounslow, R.L. Ryall and V.R. Marschall. "A Discretized Population Balance for Nucleation, Growth and Aggregation". AIChE Journal. Vol. 34 N° 11, pp. 1821-1832. November, 1988. DOI: 10.1002/aic.690341108.
[34] J.D. Litster, D.J. Smit and M.J. Hounslow. "Adjustable Discretization Population Balance for Growth and Aggregation". AIChE Journal. Vol. 41 N° 3, pp. 591-603. 1995.DOI: 10.1002/aic.690410317.
[35] Ramkrishna, Population Balances: Theory and Applications to Particulate Systems in Engineering. Academic Press, San Diego, CA. 2000.ISBN: 0125769709. DOI: 10.1016/j.ces.2007.08.03.
[36] L. Hagesaether, H.A. Jakobsen, K. Hjarbo and H.F. Svendsen. "A coalescence and breakup module for implementation in CFD codes. In European Symposium on Computer Aided Process Engineering". Vol. 1, pp. 367-372. 2000.DOI:10.1016/S1570-7946(00)80063-2.
[37] J. Sanyal, D.L. Marchisio, R.O. Fox and K. Dhanasekharan. "On the comparison between population balance models for CFD simulation of bubble columns. Industrial Engineering and Chemical Research". Vol. 44, pp. 5063-5072. 2005. DOI: 10.1021/ie049555j.
[38] Gambit 2. User's Manual to GAMBIT 2. Centrera Resource Park, 10 Cavendish Court. Fluent Inc. Lebanon, USA. 2004.
[39] F. Flórez, F. Pérez and M. de la Torre. "Scale-up of Bacillus thuringiensis fermentation based on oxygen transfer. Journal of Fermentation and Bioengineering". Vol. 83 N° 6, pp. 561-564. 1997. DOI: 10.1016/S0922-338X(97)81137-3.
[40] S. J. Arjunwadkar, K. Sarvanan, P. Kulkarni, and A.B. Pandit. "Gas-Liquid Mass Transfer in Dual Impeller Bioreactor". Biochem. Eng. Vol. 1 N° 2, pp. 99-106. March 1998. ISSN: 1369-703X.
[41] V.B. Shukla, U Parasu, P.R. Kulkarni and A.B. Pandit. "Scale-up Biotransformation Process in Stirred Tank Reactor using Dual Impeller Bioreactor". Biochem. Eng. Vol. 8, pp. 19-29. 2001.DOI: 10.1016/S1369-703X(00)00130-3.
[42] P.M. Doran. "Bioprocess Engineering Principles". Elsevier Science and Technology Books. University of New South Wales Sydney. Australia. 1995. ISBN: 978-0-12-220856-0.
[43] S. Nagata. "Mixing: Principles and Application". New York, NY. John Wiley & Sons, pp. 59-62. 1975.
[44] K.R. Westerterp, L.L. van Dierendonck and J.A. Kraa. "Interfacial areas in agitated gas-liquid contactors". Chem. Eng. Sci. Vol. 18, pp. 157-176. 1963. DOI:10.1016/0009-2509(63)85002-2.
[45] V.V. Ranade. "Computational Flow Modeling for Chemical Reactor Engineering". Process Systems Engineering Series. Academic Press. 2002.
Recibido 6 de diciembre de 2011, aceptado 22 de abril de 2013
