Mixing of miscible liquids: Dimensionless scaling for intermediate-to-large density differences in a stirred tank

Este estudio utiliza simulaciones numéricas de un tanque agitado con líquidos miscibles para demostrar que, si bien el tiempo de mezcla se correlaciona positivamente con el número de Richardson, una ley de escalado exponencial derivada basada en los números de potencia, Froude y Richardson colapsa con éxito todos los datos en una única curva maestra para diferencias de densidad intermedias a grandes.

Lo esencial

Imagina que estás intentando mezclar un tanque gigante de dos líquidos diferentes: un jarabe pesado y espeso en el fondo y un jugo más ligero y delgado en la parte superior. Dejas caer una paleta giratoria gigante (un impulsor) para mezclarlos.

En el mundo real, esta es una tarea común en fábricas que producen todo, desde medicamentos hasta tratamiento de aguas residuales. Pero aquí está el truco: como los líquidos tienen diferentes pesos (densidades), el pesado quiere quedarse en el fondo y el ligero quiere flotar en la parte superior. Esto crea una "lucha" entre la paleta giratoria que intenta mezclarlos y la gravedad que intenta mantenerlos separados.

Este artículo es como una historia de detectives donde los científicos utilizaron potentes simulaciones por computadora para determinar exactamente cuánto tiempo tarda en mezclarse perfectamente estos dos líquidos, y cómo predecir ese tiempo sin tener que construir un tanque físico y realizar pruebas costosas cada vez.

El Escenario: Una Cocina de Pruebas Digital

Los investigadores construyeron una versión virtual de un tanque de mezcla industrial estándar.

• El Tanque: Es un gran cilindro con paredes y cuatro aletas verticales (deflectores) para evitar que el líquido simplemente gire en círculo como un río perezoso.
• La Paleta: Una hoja giratoria en el medio.
• Los Líquidos: Simularon una mezcla 50/50 de un líquido pesado y uno más ligero. No utilizaron químicos reales; simplemente los trataron como fluidos "pesados" y "ligeros" con el mismo espesor (viscosidad).
• El Método: En lugar de usar ecuaciones matemáticas estándar, utilizaron un truco inteligente llamado Método de Boltzmann en Red. Piensa en esto como simular el líquido no como una masa continua, sino como miles de millones de bolas de billar diminutas e invisibles rebotando y chocando entre sí. Esto les permitió ver exactamente cómo se comportaba la turbulencia (el remolino caótico).

La Gran Pregunta: ¿Qué Tan Rápido Podemos Mezclar?

El objetivo principal era encontrar una "fórmula mágica" para predecir el tiempo de mezcla.

• Las Variables: Cambiaron dos cosas principales:
  1. Qué tan rápido gira la paleta (número de Reynolds): Girar más rápido generalmente significa más turbulencia y mezcla más rápida.
  2. Qué tan diferentes son los pesos (número de Richardson): Si los líquidos tienen casi el mismo peso, se mezclan fácilmente. Si uno es mucho más pesado, la gravedad lucha contra la mezcla, creando capas que son difíciles de romper.

El Descubrimiento: La Batalla "Gravedad vs. Giro"

Los investigadores encontraron algunos patrones interesantes:

1. Cuando la Gravedad No Importa (Mismo Peso):
   Si los dos líquidos tienen exactamente el mismo peso, el tiempo de mezcla es sorprendentemente consistente. No importa qué tan rápido gires la paleta (dentro de cierto rango), el "tiempo de mezcla adimensional" (una forma elegante de decir "cuántas vueltas de la paleta se necesitan") se mantiene constante en aproximadamente 20 vueltas. Es como una regla de la naturaleza: una vez que el agua está revolviendo lo suficiente, girarla más rápido no la hace mezclarse más rápido en términos de vueltas de la paleta.

2. Cuando la Gravedad Lucha (Pesos Diferentes):
   Tan pronto como los líquidos tienen pesos diferentes, el pesado quiere quedarse en el fondo. Cuanto mayor sea la diferencia de peso, más difícil es mezclar.

   • La Tendencia: Cuanto más diferentes son los pesos, más tiempo tarda en mezclarse.
   • El Giro Sorprendente: Si mantienes constante la "diferencia de peso" y simplemente giras la paleta más rápido, el tiempo de mezcla no siempre disminuye. A veces, girar más rápido en realidad hace que tome más tiempo alcanzar un punto específico de mezcla.
   • ¿Por qué? Imagina que el líquido pesado es como una manta gruesa. Si giras la paleta demasiado rápido, creas mucha energía, pero el líquido pesado forma una "tapa" estable que el líquido ligero no puede penetrar. La energía se desperdicia en revolver la capa superior mientras la capa inferior permanece sellada. Es como intentar revolver una olla de sopa donde los vegetales pesados se han asentado en un bloque sólido en el fondo; girar la cuchara más rápido solo salpica el caldo de arriba sin romper el bloque de vegetales.

La Solución: Una Nueva "Curva Maestra"

El mayor logro del equipo fue crear una sola fórmula simple que combine todos estos factores. Descubrieron que si observas el tiempo de mezcla a través de la lente de tres números específicos (Potencia, Froude y Richardson), todos sus datos desordenados colapsan en una sola curva suave y exponencial.

Piénsalo así: Antes, los ingenieros tenían que adivinar o realizar cientos de pruebas para ver cómo se mezclaría un nuevo líquido. Ahora, tienen una "receta". Si les dices la diferencia de peso y qué tan rápido planeas girar, esta fórmula predice el tiempo de mezcla con alta precisión.

La Conclusión

El artículo concluye que para estos tanques industriales específicos:

• La turbulencia es clave: Una vez que el líquido está revolviendo completamente, el comportamiento de mezcla es predecible.
• La gravedad es el jefe: Si los líquidos tienen diferentes densidades, la gravedad crea una "estratificación" (capas) que resiste la mezcla.
• Más rápido no siempre es mejor: En sistemas con grandes diferencias de densidad, simplemente aumentar la velocidad del motor no garantiza una mezcla más rápida; a veces solo crea una separación más estable.

Los autores proporcionan esta nueva fórmula para ayudar a los ingenieros a diseñar mejores procesos de mezcla sin necesidad de construir prototipos costosos primero. Planean probar esta fórmula en diferentes formas de tanques y tipos de paletas en el futuro, pero por ahora, funciona perfectamente para los tanques estándar que simularon.

Resumen técnico

Resumen Técnico: Escalamiento Adimensional para Diferencias de Densidad Intermedias a Grandes en un Tanque Agitado

Enunciado del Problema
La mezcla de líquidos miscibles es una operación unitaria crítica en procesos industriales, que abarca desde la producción farmacéutica hasta el tratamiento de aguas residuales. Si bien la mezcla turbulenta en sistemas monofásicos y homogéneos está bien caracterizada, la presencia de diferencias de densidad entre los fluidos introduce estratificación impulsada por flotabilidad. Este fenómeno crea regiones estratificadas verticalmente donde las fuerzas gravitatorias contrarrestan el transporte convectivo, suprimiendo potencialmente la mezcla vertical y aumentando el tiempo requerido para alcanzar la homogeneidad. Aunque existen correlaciones empíricas para la mezcla de trazadores en fase única y casos estratificados específicos, sigue siendo esquivo un marco generalmente válido que cuente de manera consistente con los efectos gravitatorios, inerciales, viscosos y de flotabilidad a través de diferentes regímenes de flujo y geometrías. Este estudio aborda el desafío de predecir los tiempos de mezcla en tanques agitados con baffles que contienen dos líquidos miscibles con diferencias de densidad intermedias a grandes.

Metodología
Los autores emplearon Dinámica de Fluidos Computacional (CFD) utilizando el Método de Boltzmann en Red (LBM) para simular la mezcla turbulenta en un tanque agitado con baffles. El sistema consistió en una relación de volumen 50/50 de dos líquidos miscibles con viscosidades dinámicas idénticas pero densidades variables ($\rho_l$ y $\rho_h$), donde el fluido más ligero estaba inicialmente estratificado sobre el fluido más pesado.

• Enfoque Numérico: Las simulaciones resolvieron las ecuaciones de Navier-Stokes incompresibles utilizando un esquema de red D3Q19. Para manejar altos números de Reynolds ($Re > 10^5$), se implementó un modelo de Simulación de Grandes Remolinos (S-L LES) de Smagorinsky-Lilly para modelar la turbulencia subred.
• Transporte Escalar: Se adoptó un modelo de transporte escalar basado en flujo másico para rastrear la concentración del líquido más ligero. Este modelo acopló el campo escalar a la ecuación de momento mediante una fuerza de cuerpo gravitatoria dependiente de la concentración, permitiendo la simulación de efectos de flotabilidad derivados de variaciones locales de densidad.
• Espacio de Simulación: Se realizó un estudio sistemático de parámetros variando la velocidad de rotación del impulsor ($N$) y la densidad de la fase ligera ($\rho_l$). Esto resultó en un amplio rango de números adimensionales: números de Reynolds ($Re$) de 50.000 a 125.000 y números de Richardson ($Ri$) de 0 a 14.
• Validación: La precisión numérica se verificó mediante un análisis del Índice de Convergencia de Malla (GCI), que estimó una incertidumbre numérica de aproximadamente 2,39 % en la malla de producción. Además, los resultados de la simulación se validaron contra datos experimentales de trazadores de la literatura, mostrando un buen acuerdo en la respuesta de la sonda y en las definiciones del tiempo de mezcla.

Resultados Clave
El estudio cuantificó el tiempo de mezcla ($t_m$), definido como el tiempo necesario para alcanzar el 90 % de homogeneidad, y lo expresó en forma adimensional ($t^*_m = N \cdot t_m$).

1. Línea Base Homogénea ($Ri = 0$): En ausencia de diferencias de densidad, el tiempo de mezcla adimensional colapsó a un valor constante de $t^*_m \approx 20$ en todo el rango de números de Reynolds investigado. Esto confirma que el flujo es completamente turbulento y que la mezcla inercial domina, consistente con la literatura establecida para tanques con baffles.
2. Efecto de la Estratificación ($Ri > 0$): A medida que aumentó el número de Richardson, indicando mayores diferencias de densidad, el tiempo de mezcla adimensional aumentó sistemáticamente. Los datos siguieron una tendencia exponencial clara en todos los números de Reynolds.
3. Ley de Escalamiento Propuesta: Los autores derivaron una curva maestra que colapsa todos los datos de simulación en una única relación de escalamiento exponencial:
   $$t^*_m = 20 \cdot e^{0.58 \cdot N_p \cdot Fr \cdot Ri}$$
   Donde $N_p$ es el número de potencia, $Fr$ es el número de Froude y $Ri$ es el número de Richardson. El prefactor de 20 representa la línea base inercial, mientras que el término exponencial captura el efecto estabilizador de la flotabilidad.
4. Comportamiento No Monótono: Al analizar los tiempos de mezcla dimensionales ($t_m$) a $Ri$ constante, el estudio observó una dependencia no monótona con respecto a $Re$ para números de Richardson más altos ($Ri \geq 8$). Si bien aumentar $Re$ (y por tanto la velocidad del impulsor) inicialmente reduce el tiempo de mezcla al mejorar la turbulencia, aumentos adicionales conducen a una estratificación de densidad más fuerte (debido al acoplamiento de parámetros en la configuración de la simulación), lo que suprime el transporte vertical. Esta interacción resulta en un tiempo de mezcla mínimo en números de Reynolds intermedios.
5. Estructura del Flujo: El análisis del campo de velocidades reveló que a altos $Ri$, la fase ligera puede quedar "sellada" de la entrada de potencia del impulsor, apareciendo estructuras verticales principalmente en el líquido más denso. El análisis de la Densidad Espectral de Potencia (PSD) indicó que, si bien la flotabilidad afecta la organización espacial del flujo, la dinámica impulsada por el agitador sigue dominando la cascada de energía.

Significado y Afirmaciones
El artículo afirma proporcionar una correlación predictiva compacta para los tiempos de mezcla turbulenta en sistemas inicialmente estratificados por densidad dentro del régimen completamente turbulento. La ley de escalamiento derivada, basada en la influencia combinada de los números de Potencia, Froude y Richardson, ofrece una expresión analítica simple para ingenieros y expertos en aplicaciones.

Los autores señalan que, hasta donde ellos conocen, no se ha reportado previamente una correlación comparable para sistemas estratificados en el régimen completamente turbulento. El significado de este trabajo radica en su capacidad para agilizar el diseño de procesos de mezcla al reducir la dependencia de costosas pruebas en el mundo real. El modelo propuesto permite estimar los tiempos de mezcla en sistemas donde los efectos de flotabilidad son significativos, capturando la compleja competencia entre la turbulencia mejorada y la estratificación reforzada. Los autores sugieren modestamente que el trabajo futuro implicará probar esta correlación en diferentes geometrías de tanques y tipos de impulsores, incluidos agitadores angulados, para evaluar la robustez del escalamiento propuesto.