Tres intérpretes y la danza del vórtice ✛

‹›

Un tubo vertical con expansión repentina es un campo de pruebas clásico para la turbulencia. El flujo de aire, girado como un remolino, irrumpe en una cámara ancha y genera vórtices, corrientes de retorno y zonas donde partículas de vidrio de 45 μm de diámetro, la mitad del grosor de un cabello, se mueven en trayectorias complejas. Martin Sommerfeld y su colega Hui Qiu midieron las velocidades en ese flujo ya en la década de 1990, y desde entonces su experimento se ha convertido en un referente para validar modelos de flujos bifásicos. Tres versiones populares del modelo k-ε (una forma de estimar las fluctuaciones turbulentas a través de la energía y la disipación) pretenden describir con precisión esta danza de vórtices. Cada una es como un intérprete que escucha la misma música pero la transmite a su manera.

Los intérpretes difieren en los detalles. El modelo k-ε estándar, el más antiguo y simple, resultó ser el más fiel: sus perfiles de velocidad coincidieron con el experimento en las cuatro secciones (x = 3, 52, 112 y 155 mm). El modelo realizable, creado posteriormente para corregir las deficiencias del estándar, falló en la velocidad radial: lejos de los chorros predijo un flujo hacia adentro excesivamente fuerte, y cerca de la pared hacia afuera, mientras que en el experimento el flujo estaba dirigido hacia adentro. El modelo RNG, derivado de la mecánica estadística, generó zonas de recirculación adicionales al comienzo de la burbuja central, vórtices falsos que no aparecieron en el experimento. En cuanto al tiempo de cómputo para un intervalo de 0,5 s, el modelo estándar resultó ser también el más rápido: 17,13 horas frente a 19,44 del RNG y 24,81 del realizable (un 45% más).

Las partículas también se pueden describir de dos maneras: seguir cada una por separado (enfoque lagrangiano, que rastrea cada partícula a lo largo de su trayectoria) o agruparlas en paquetes, llamados parcels, donde cien partículas idénticas se mueven como una sola. La lógica sugiere que agrupar ahorrará tiempo, ya que habrá menos ecuaciones de movimiento. Pero en esta prueba el ahorro fue solo del 4,2% (31,15 horas frente a 29,84), casi nada. La razón es simple: la mayor carga computacional recae en resolver las ecuaciones de la fase gaseosa, y hay pocas partículas en el flujo. Los perfiles promedio de velocidad del gas no cambiaron en absoluto al agrupar hasta 100 partículas por paquete; la diferencia es imperceptible.

Este trabajo ofrece a los ingenieros una guía clara: para flujos bifásicos giratorios con carga moderada de partículas, el modelo k-ε estándar sigue siendo un compromiso razonable entre precisión y costo computacional. El modelo realizable, a pesar de sus mejoras teóricas, puede dar una velocidad radial cualitativamente incorrecta, y el RNG, estructuras de vórtices falsas. El enfoque de agrupación de partículas solo se justifica con alta concentración de la fase dispersa, cuando el ahorro en el número de ecuaciones de movimiento de partículas se vuelve notable. Hoy en día, este conocimiento se aplica al diseño de cámaras de combustión de carbón pulverizado, separadores ciclónicos y sistemas de gasificación con captura de CO2, donde es importante predecir la distribución de partículas y las zonas de recirculación.

El siguiente paso natural es repetir la comparación con una carga de partículas más alta (superior a 0,034 kg/kg) y un giro más fuerte (número de giro mayor que 0,6), característico de los quemadores reales. Quizás sea necesario adaptar las funciones de pared (una descripción aproximada del flujo cerca de la pared) o pasar a modelos que resuelvan la subcapa viscosa para describir mejor las corrientes secundarias. Pero ya está claro: en el mundo de los flujos bifásicos turbulentos, un intérprete confiable vale más que uno de moda.