An LES model with finite-rate phase change and subgrid spray based on a thermodynamically consistent four-equation multiphase model

Este trabajo propone un modelo de simulación LES que integra un modelo multifásico de cuatro ecuaciones con cambios de fase de tasa finita y un modelo de spray de subrejilla basado en la ecuación $\Sigma$, logrando predicciones precisas y termodinámicamente consistentes validadas con datos experimentales de combustión.

Lo esencial

El Gran Problema: ¿Cómo simular una lluvia de combustible en un motor?

Imagina que quieres diseñar el motor de un cohete o de un coche de carreras súper avanzado. Para que el motor sea eficiente, necesitas saber exactamente cómo el combustible (que sale como un chorro líquido a toda presión) se convierte en gas para quemarse.

El problema es que esto es un caos absoluto. El combustible se rompe en gotas diminutas, esas gotas se evaporan, algunas explotan por el cambio de presión (como cuando abres una lata de refresco muy rápido) y todo ocurre en un espacio tan pequeño que las computadoras actuales "sufren" para calcularlo todo con precisión. Si intentas calcular cada molécula, la computadora tardaría años. Si lo haces de forma muy simple, el motor que diseñes podría explotar en la vida real porque tus cálculos fallaron.

La Solución: "El truco de la aproximación inteligente"

Los científicos de este estudio (de Stanford y la Universidad de KTH) han creado un nuevo modelo matemático. Para entenderlo, vamos a usar tres analogías:

1. El Modelo de las Cuatro Ecuaciones (La "Regla de la Paz")

Imagina que tienes una fiesta con dos grupos de personas: los "Líquidos" y los "Gases". En una simulación normal, tendrías que vigilar la velocidad, la temperatura y la presión de cada persona por separado. ¡Es imposible!

Lo que estos científicos hacen es aplicar una regla de "paz inmediata": asumen que, aunque sean grupos distintos, la presión y la temperatura se equilibran casi al instante. Es como decir: "No me importa qué tan rápido corra cada persona, asumo que todo el salón tiene la misma temperatura y la misma presión". Esto ahorra una cantidad enorme de energía a la computadora, permitiéndole trabajar mucho más rápido sin perder la esencia de lo que pasa.

2. El Modelo de Spray Sub-rejilla (El "Efecto Niebla")

Las computadoras dividen el mundo en pequeños cubitos (como si fueran piezas de LEGO). Si una gota de combustible es más pequeña que un cubito de LEGO, la computadora no puede "verla". Es como intentar ver un grano de arena desde un avión; para ti, es solo una mancha borrosa.

En lugar de intentar dibujar cada granito de arena (que sería imposible), los científicos inventaron un modelo de "Niebla Inteligente" ($\Sigma$). En lugar de decir "aquí hay una gota", la computadora dice: "No veo la gota, pero sé que hay una cantidad de 'superficie de contacto' aquí, como si fuera una niebla espesa". Esto permite que la computadora sepa que hay combustible ahí y cuánto se puede evaporar, aunque no pueda dibujar la gota individualmente.

3. Cambio de Fase de Tasa Finita (El "Semáforo de la Evaporación")

Antes, los modelos decían que el combustible pasaba de líquido a gas de forma instantánea (como si fuera un interruptor de luz: o está encendido o está apagado). Pero en la realidad, la naturaleza es más lenta; es más como un semáforo.

Los autores crearon un modelo que respeta las leyes de la termodinámica (la "policía de la energía"). El combustible no puede evaporarse más rápido de lo que la física permite. Han puesto un "límite de velocidad" para que la simulación no haga cosas locas o imposibles, como hacer aparecer gas de la nada.

¿Para qué sirve esto en la vida real?

Para validar su invento, probaron su modelo con un experimento real llamado ECN Spray A (un chorro de combustible n-dodecano). Los resultados fueron casi idénticos a lo que ocurre en el laboratorio.

En resumen: Han creado un "mapa inteligente" que permite a las computadoras predecir cómo se comporta el combustible en un motor de forma rápida y muy precisa, sin necesidad de ser tan detallistas que la computadora se bloquee, pero sin ser tan descuidados que el resultado sea mentira.

Esto ayudará a crear motores más limpios, cohetes más potentes y sistemas de extinción de incendios más efectivos.

Resumen técnico

Resumen Técnico: Modelo LES con Cambio de Fase de Tasa Finita y Spray de Subrejilla

Título original: An LES model with finite-rate phase change and subgrid spray based on a thermodynamically consistent four-equation multiphase model

1. El Problema (Contexto y Motivación)

La simulación de sistemas multifásicos y multicomponentes a escala de ingeniería (como motores de combustión o propulsión de cohetes) es extremadamente compleja debido a la necesidad de capturar la mezcla intrafase, interfaces con altas relaciones de densidad y procesos de transferencia de masa.

Los modelos actuales enfrentan un dilema de eficiencia y precisión:

• Modelos de no-equilibrio total: Son físicamente precisos pero computacionalmente prohibitivos para escalas de ingeniería.
• Modelos de equilibrio homogéneo (HEM): Son eficientes pero sobreestiman la tasa de transferencia de masa al asumir que el equilibrio es instantáneo, lo cual es incorrecto para estructuras de interfaz no resueltas por la malla (subrejilla).
• Modelos de tasa finita existentes: A menudo carecen de una base termodinámica sólida o no integran correctamente la dependencia de la transferencia de masa con el área superficial de la interfaz en entornos de simulación de grandes remolinos (LES).

2. Metodología

Los autores proponen un marco de trabajo basado en un modelo multifásico de cuatro ecuaciones que asume equilibrio mecánico (presión) y térmico (temperatura) entre fases, lo que reduce drásticamente el costo computacional. Para compensar la restricción del equilibrio, introducen dos innovaciones principales:

• Cambio de fase de tasa finita termodinámicamente acotado: En lugar de usar tasas arbitrarias, proponen un modelo basado en la cinética de Hertz-Knudsen. Para evitar resultados no físicos (como la transferencia de masa más allá del equilibrio), el modelo está "acotado" por la energía libre de Gibbs. La tasa de transferencia de masa se define mediante un tiempo de relajación ($\tau$) que garantiza que el sistema siempre tienda hacia el estado de equilibrio químico definido por el modelo HEM.
• Modelo de Spray de Subrejilla ($\Sigma$): Para capturar el área interfacial que la malla no puede resolver, utilizan una ecuación de transporte para la densidad de área superficial ($\Sigma'$). Introducen una formulación "phase-confined" (confinada a la fase) para evitar la "fuga" (leakage) de la subrejilla a través de las interfaces, asegurando que el área de la subrejilla solo exista dentro del volumen de la fase gaseosa.
• Esquema Numérico: Utilizan una técnica híbrida de alta resolución que combina esquemas que preservan la energía cinética y la entropía con esquemas de Godunov para capturar gradientes agudos (choques e interfaces).

3. Contribuciones Clave

1. Consistencia Termodinámica: Desarrollo de un modelo de cambio de fase de tasa finita que utiliza el potencial químico como fuerza impulsora, asegurando que la evaporación, condensación y flashing (evaporación súbita por caída de presión) sean físicamente consistentes.
2. Modelo $\Sigma$ para LES: Una nueva formulación de la ecuación de transporte de área superficial que es robusta frente a la difusión numérica y que permite modelar procesos de ruptura turbulenta (dense spray) y ruptura secundaria (dilute spray), así como la coalescencia.
3. Integración de Regímenes: Un sistema de indicadores (basado en la fracción de volumen y la escala de Hinze) que permite al modelo alternar automáticamente entre regímenes de spray denso y diluido, y entre procesos de flashing y evaporación convencional.

4. Resultados y Validación

El modelo se validó utilizando el caso de estudio ECN Spray A (inyección de n-dodecano en nitrógeno) en dos condiciones:

• Caso No-Evaporante: Se validó la dinámica del spray y el modelo de área superficial ($\Sigma$). Los resultados de la densidad de masa proyectada (PMD) y el diámetro medio de Sauter (SMD) mostraron una excelente concordancia con los experimentos, demostrando que el modelo captura correctamente la expansión del chorro y la ruptura de gotas.
• Caso Evaporante: Se validó el acoplamiento completo (dinámica de spray + cambio de fase). Las predicciones de la longitud de penetración del líquido y la extensión de la nube de gas evaporado coincidieron estrechamente con las mediciones experimentales, validando la capacidad del modelo para manejar la transferencia de masa compleja.

5. Significancia

Este trabajo representa un avance importante en la simulación de fluidos multifásicos porque demuestra que es posible utilizar un modelo de cuatro ecuaciones (más rápido y económico) sin sacrificar la precisión predictiva, siempre y cuando se incorporen modelos de subrejilla que respeten las leyes de la termodinámica. Esto abre la puerta a simulaciones de alta fidelidad de motores de combustión interna y sistemas de propulsión que antes eran computacionalmente inalcanzables.