A Temperature-Coupled Cahn-Hilliard-Stokes-Heat Model for Thermally Driven Phase Separation

Resumen Técnico: Modelo Acoplado de Cahn–Hilliard–Stokes–Calor

1. Planteamiento del Problema

El trabajo aborda la modelización matemática y numérica de la separación de fases en mezclas viscosas e incompresibles impulsada por gradientes de temperatura. El contexto motivador es la condensación de humedad dentro de los alojamientos de baterías en vehículos eléctricos, un fenómeno que causa corrosión y degradación del aislamiento.

El desafío principal radica en capturar la interacción compleja entre:

• La separación de fases (transición de vapor a líquido).
• El flujo de fluido (movimiento del fluido debido a diferencias de densidad y viscosidad).
• El transporte de calor (que controla la termodinámica de la transición).

A diferencia de modelos anteriores que asumen isotermia o acoplamientos ad-hoc, este estudio propone un sistema acoplado donde la temperatura modifica directamente la energía libre del sistema, alterando la estabilidad de las fases.

2. Metodología y Formulación Matemática

A. El Modelo Físico-Matemático
Los autores desarrollan un modelo de interfaz difusa que combina tres subsistemas acoplados:

1. Ecuación de Cahn–Hilliard (Convección-Difusión):

   • Describe la evolución del parámetro de orden $c$ (fracción volumétrica o concentración).
   • Incluye un término convectivo $(u \cdot \nabla)c$ debido al flujo del fluido.
   • Potencial Termodinámico: Se introduce un potencial de Landau dependiente de la temperatura:
     $$f(c, \Theta) = 2c^4 - (1 - \beta(\Theta))c^2$$
     Donde $\beta(\Theta) = \Theta / \Theta_S$ es la temperatura adimensionalizada respecto a la temperatura de saturación $\Theta_S$.
     • Si $\Theta > \Theta_S$ ($\beta > 1$): El potencial es de un solo pozo (fase homogénea estable).
     • Si $\Theta < \Theta_S$ ($\beta < 1$): El potencial es de doble pozo (inestabilidad espinodal, separación de fases).

2. Sistema de Stokes (Flujo):

   • Describe el flujo de fluido en régimen de bajo número de Reynolds (flujo lento/creep).
   • Incluye densidad $\rho(c)$ y viscosidad $\eta(c)$ que varían linealmente con el parámetro de orden.
   • Incorpora un término de fuerza de flotabilidad (gravedad) $G \rho(c) \hat{g}$.
   • Condiciones de frontera mixtas: no deslizamiento en algunas paredes y tensión libre en otras.

3. Ecuación de Calor:

   • Ecuación de advección-difusión para la temperatura $\Theta$.
   • La temperatura alimenta el potencial de Cahn–Hilliard, cerrando el ciclo de retroalimentación.

B. Escalamiento Adimensional
El sistema se adimensionaliza utilizando números clave:

• Número de Péclet de masa ($Pe$): Relación entre transporte convectivo y difusivo de masa.
• Número de Péclet térmico ($Pe_\Theta$): Relación entre transporte convectivo y difusivo de calor.
• Número de Cahn ($Ch$): Relación entre el grosor de la interfaz y el tamaño del dominio.
• Parámetro gravitacional ($G$): Relación entre fuerzas gravitatorias y viscosas.

3. Contribuciones Clave

El artículo ofrece contribuciones tanto analíticas como numéricas:

A. Teoría de Existencia (Análisis Matemático)

• Existencia Local: Se prueba la existencia de soluciones débiles locales en el tiempo para el sistema acoplado regularizado en dimensiones 2 y 3.
• Estrategia de Punto Fijo: Se utiliza el Teorema del Punto Fijo de Schauder. El sistema se descompone en tres operadores:
  1. Operador de Calor ($H$): Mapea velocidad a temperatura.
  2. Operador de Cahn–Hilliard ($C$): Mapea velocidad y temperatura al parámetro de orden.
  3. Operador de Stokes ($S$): Mapea concentración a velocidad.
• Regularización: Para garantizar la existencia, se introducen regularizaciones analíticas (un término de derivada temporal en Stokes y una suavización espacial de la temperatura en el potencial químico) que no se utilizan en la simulación numérica final.
• Marco Funcional: El análisis se realiza en espacios de Sobolev ($H^2 \to L^2 \to H^{-2}$) para manejar la falta de cotas $L^\infty$ globales para el parámetro de orden cuando se usa un potencial polinómico.

B. Esquema Numérico

• Discretización: Se propone un esquema de elementos finitos totalmente discreto.
• Estrategia de Descomposición (Splitting):
  1. Paso de Calor: Discretización implícita (Euler hacia atrás).
  2. Paso de Stokes: Resolución implícita con coeficientes variables (densidad/viscosidad).
  3. Paso de Cahn–Hilliard: Discretización convex-concava (convex splitting) estabilizada. Esto asegura la estabilidad energética incondicional para el subsistema de difusión isotérmica.
• Implementación: Se utiliza el software FreeFEM++.

4. Resultados de las Simulaciones Numéricas

Los autores realizaron cuatro experimentos numéricos en 2D para validar el modelo:

1. Validación con Flujo de Cavidad (Benchmark):

   • Comparación con un caso de cavidad impulsada por una tapa (lid-driven cavity) con interfaz difusa.
   • Resultado: El esquema reproduce correctamente los perfiles de velocidad y la deformación de la interfaz, validando la acoplamiento hidrodinámico.

2. Separación de Fases Isotérmica:

   • Se estudia la evolución sin flujo ni gradiente térmico.
   • Resultado: Se confirma la transición controlada por temperatura:
     • $\Theta > \Theta_S$: Mezcla homogénea (un solo pozo).
     • $\Theta < \Theta_S$: Separación de fases rápida seguida de coarsening (crecimiento de dominios).

3. Separación de Fases por Gravedad:

   • Se activa la flotabilidad con densidades y viscosidades diferentes.
   • Resultado: Las fases densas migran hacia abajo y se acumulan en la pared inferior, demostrando la capacidad del modelo para capturar la sedimentación y la formación de capas bajo condiciones de frontera mixtas.

4. Fase Cambiante Inducida Térmicamente (Condensación/Evaporación):

   • Se aplica un gradiente térmico en los límites.
   • Caso Calentado: Una pared caliente suprime la separación de fases localmente, disolviendo las estructuras de fase.
   • Caso Enfriado: Una pared fría induce la nucleación de fases densas (condensación) cerca de la pared, que luego crecen.
   • Resultado: El modelo captura cualitativamente la nucleación de gotas, el crecimiento y la interacción con las paredes, reproduciendo fenómenos relevantes para la gestión térmica de baterías.

5. Significado y Conclusión

Este trabajo es significativo por varias razones:

• Marco Unificado: Proporciona un marco matemático riguroso que acopla termodinámica, flujo de fluidos y transporte de calor en un solo sistema de ecuaciones, superando las limitaciones de los modelos isotérmicos.
• Aplicabilidad Práctica: Es directamente aplicable a problemas de ingeniería como la gestión de humedad en baterías de vehículos eléctricos, donde la condensación es un riesgo crítico.
• Robustez Numérica: El esquema propuesto demuestra ser una herramienta robusta para simular geometrías confinadas y condiciones de frontera complejas, mostrando cómo parámetros adimensionales (Péclet, Cahn, Gravedad) influyen en la formación de patrones.
• Fundamento Teórico: La prueba de existencia de soluciones débiles para este sistema acoplado con condiciones de frontera mixtas llena un vacío en la literatura matemática sobre modelos de interfaz difusa no isotérmicos.

En resumen, el artículo establece una base sólida tanto teórica como computacional para estudiar la separación de fases impulsada térmicamente en sistemas complejos, ofreciendo una herramienta predictiva para el diseño y la seguridad de sistemas electroquímicos.