Upscaling the Navier-Stokes-Cahn-Hilliard model for incompressible multiphase flow in inhomogeneous porous media

Este trabajo presenta un modelo macroscópico riguroso para el flujo bifásico en medios porosos inhomogéneos, derivado mediante promediado volumétrico de las ecuaciones de Navier-Stokes y Cahn-Hilliard, que incorpora formalmente el comportamiento de mojado en el potencial químico promediado y valida su capacidad mediante simulaciones numéricas.

Lo esencial

¡Claro que sí! Imagina que este artículo científico es como un manual de instrucciones para predecir cómo se mezclan dos líquidos (como aceite y agua) cuando intentan pasar a través de un laberinto de roca (como un suelo poroso o una esponja).

Aquí tienes la explicación, traducida a un lenguaje sencillo y con analogías divertidas:

1. El Problema: El Laberinto Invisible

Imagina que tienes una esponja gigante llena de millones de agujeros diminutos (poros). Si viertes agua y aceite sobre ella, ¿cómo se mueven?

• A nivel microscópico (el tamaño de los poros): Es un caos. Las gotas de agua y aceite chocan, se estiran, se rompen y luchan por el espacio. Los científicos pueden simular esto con superordenadores, pero es como intentar seguir a cada persona en una multitud de un millón de personas. ¡Es demasiado lento y costoso para proyectos grandes como limpiar un acuífero o extraer petróleo!
• A nivel macroscópico (el tamaño de la esponja entera): Necesitamos una regla simple que nos diga: "Si inyecto agua aquí, ¿cuánto petróleo saldrá allá?". El problema es que las reglas actuales (llamadas Ley de Darcy) son como mapas antiguos: funcionan bien en terreno llano, pero fallan cuando el terreno es complejo o cuando los líquidos tienen comportamientos extraños (como si el agua "se pegara" a las paredes de la roca).

2. La Solución: El "Traductor" de Escalas

Los autores de este paper han creado un traductor matemático. Han tomado las leyes físicas complejas que gobiernan el movimiento de los fluidos en cada agujero diminuto y las han "promediado" para crear una nueva ley que funcione para toda la esponja.

Piensa en esto como si quisieras describir el tráfico en una ciudad:

• En lugar de seguir a cada coche individualmente (nivel microscópico), quieres saber el flujo promedio de coches por hora (nivel macroscópico).
• Pero, a diferencia de los coches, los líquidos en la roca tienen una "personalidad": se pegan a las paredes (mojabilidad) y se repelen entre sí.

3. La Magia: El "Imán" Invisible (La Mojabilidad)

Lo más genial de este estudio es cómo manejan la mojabilidad.

• La analogía: Imagina que las paredes de los poros de la roca son como imanes. Si la roca es "hidrófila" (le gusta el agua), actúa como un imán que atrae al agua y empuja al aceite. Si es "hidrófoba", hace lo contrario.
• El problema anterior: Las fórmulas viejas ignoraban este "imán" o lo trataban de forma muy tosca.
• La innovación de este paper: Han creado una nueva ecuación que incluye un "potencial químico promedio". Imagina que este potencial es como un termómetro de deseo. Mide cuánto "quieren" los líquidos estar en contacto con la roca o entre ellos. Al incluir esto en la ecuación principal, el modelo puede predecir con mucha más precisión cómo se comportará el fluido, incluso si la roca tiene formas extrañas o si el líquido se pega a las paredes.

4. El Método: El "Tamiz" Matemático

Para hacer esto, usaron una técnica llamada promedio de volumen.

• Imagina que tienes un cubo de hielo con trozos de fruta dentro. Si quieres saber la "densidad promedio" de la fruta en todo el cubo, no cuentas cada trozo uno por uno. Tomas el cubo, lo agitas y calculas un promedio.
• Los autores hicieron esto con las ecuaciones de movimiento. Pero al promediar, surgieron "huecos" en la matemática (términos desconocidos). Para llenar esos huecos, resolvieron pequeños problemas matemáticos dentro de un solo "agujero" representativo y luego aplicaron esa solución a todo el sistema. Es como resolver un rompecabezas en una pieza pequeña y luego usar esa pieza para armar todo el cuadro.

5. Las Pruebas: ¿Funciona de verdad?

Para asegurarse de que su nuevo "traductor" no estaba loco, lo pusieron a prueba en tres escenarios:

1. Flujo simple: Como agua corriendo por un tubo lleno de arena. El modelo coincidió perfectamente con la teoría conocida.
2. La transición: Como pasar de un río libre a un río que entra en una zona de roca. El modelo manejó bien el cambio de velocidad.
3. El caos (Dedos viscosos): Imagina inyectar agua en aceite dentro de la roca. A veces, el agua avanza como "dedos" que se ramifican y rompen. El modelo logró predecir cómo se forman estos dedos y, lo más importante, cómo cambia la forma de los dedos si la roca es más o menos "pegajosa" (mojabilidad).

En Resumen

Este trabajo es como haber creado un GPS de alta precisión para fluidos en el subsuelo.
Antes, los ingenieros tenían que adivinar o usar reglas aproximadas para saber cómo se movería el petróleo o el agua en el suelo. Ahora, gracias a este modelo, tienen una herramienta matemática rigurosa que entiende la "personalidad" de los líquidos (cómo se pegan a la roca) y puede predecir su movimiento desde el nivel de un grano de arena hasta el de una montaña entera.

¿Por qué importa?
Esto es vital para:

• Limpiar el agua subterránea contaminada.
• Extraer petróleo de forma más eficiente.
• Almacenar CO2 bajo tierra para combatir el cambio climático.

Básicamente, han enseñado a las matemáticas a "ver" los detalles pequeños para poder predecir el comportamiento grande con mucha más inteligencia.

Resumen técnico

A continuación presento un resumen técnico detallado del artículo en español, estructurado según los puntos solicitados:

Título: Escalamiento del modelo Navier-Stokes-Cahn-Hilliard para flujo multifásico incompresible en medios porosos inhomogéneos

1. Planteamiento del Problema

El flujo de fluidos multifásicos a través de medios porosos es fundamental en aplicaciones como la recuperación mejorada de petróleo, la remediación de acuíferos y el secuestro de CO₂. Sin embargo, modelar estos procesos es un desafío debido a la complejidad geométrica de los poros y a la amplia gama de escalas temporales y espaciales involucradas.

• Limitaciones actuales: Los enfoques tradicionales a escala macroscópica (escala Darcy) se basan en leyes empíricas (como la ley de Darcy modificada) que requieren funciones de permeabilidad relativa y presión capillary ajustadas experimentalmente. Estas formulaciones a menudo fallan en capturar la física subyacente, especialmente fenómenos de histéresis, efectos de mojabilidad y dinámicas interfaciales complejas, ya que no derivan rigurosamente de las ecuaciones de gobierno a escala de poro.
• Necesidad: Existe una necesidad crítica de un marco teórico que integre rigurosamente la dinámica a escala de poro (donde ocurren los fenómenos de interfaz y mojabilidad) con las ecuaciones de transporte a escala macroscópica, evitando la dependencia excesiva de parámetros empíricos.

2. Metodología

Los autores emplean el método de promediado volumétrico para escalar las ecuaciones de gobierno desde la escala de poro hasta la escala de Darcy (o Volumen Elemental Representativo, REV).

• Ecuaciones a escala de poro: Se utiliza el sistema acoplado Navier-Stokes-Cahn-Hilliard (NSCH).
  • La evolución de la interfaz entre los dos fluidos inmiscibles se describe mediante la ecuación de Cahn-Hilliard, basada en un funcional de energía libre que incluye términos de energía volumétrica, energía interfacial y energía de pared (mojabilidad).
  • El flujo de momento se rige por las ecuaciones de Navier-Stokes, incluyendo fuerzas interfaciales derivadas del potencial químico.
• Proceso de Promediado:
  • Se aplican teoremas de promediado espacial y temporal sobre un REV.
  • Se utiliza la descomposición de Gray, separando las cantidades físicas en una componente promediada (lenta) y una desviación espacial (rápida).
  • Se derivan ecuaciones para las variables promediadas, lo que genera términos "no cerrados" (desviaciones espaciales) que deben modelarse.
• Cierre del Modelo:
  • Se resuelven problemas de cierre locales dentro del REV para expresar las desviaciones en función de las variables promediadas.
  • Innovación clave en el cierre: Se incorpora formalmente la condición de contorno de mojabilidad en el potencial químico promediado. Esto permite que el efecto de la mojabilidad (ángulo de contacto) aparezca naturalmente en la ecuación macroscópica sin necesidad de funciones empíricas de presión capilar.
  • La fuerza de arrastre (drag) se modela utilizando una formulación que incluye términos de Darcy (baja velocidad) y Forchheimer (alta velocidad), junto con viscosidad efectiva.

3. Contribuciones Clave

1. Modelo de Fluido Único Riguroso: A diferencia de los enfoques de "dos fluidos" que requieren ecuaciones de momento separadas para cada fase, este trabajo deriva un modelo de fluido único a escala macroscópica. Esto simplifica la estructura matemática y evita la necesidad de condiciones de contorno interfaciales empíricas explícitas.
2. Derivación Formal de la Mojabilidad: El aporte más significativo es la inclusión rigurosa de la mojabilidad dentro del potencial químico promediado. La condición de contorno de la pared a escala de poro se transfiere directamente al término de fuerza capilar macroscópico, vinculando intrínsecamente el ángulo de contacto con la dinámica del flujo.
3. Marco Unificado de Escalas: El modelo actúa como un puente teórico que conecta la física a escala de poro con la escala de Darcy. Se demuestra que, bajo ciertas condiciones (porosidad uniforme, baja velocidad), el modelo se reduce a las ecuaciones de Darcy-Brinkman clásicas, pero con términos de tensión superficial y mojabilidad derivados teóricamente en lugar de empíricamente.
4. Implementación Numérica: Se desarrolló un método de Lattice Boltzmann (LBM) adaptado para resolver las nuevas ecuaciones promediadas, validando la viabilidad computacional del modelo.

4. Resultados

Los autores validaron el modelo mediante simulaciones numéricas de varios casos de prueba:

• Flujo de Poiseuille en medio poroso homogéneo: Los resultados numéricos mostraron una excelente concordancia con las soluciones analíticas para perfiles de velocidad, confirmando la precisión de las ecuaciones hidrodinámicas promediadas.
• Flujo en canal fluido-poroso: Se simuló el flujo en una interfaz entre un fluido libre y un medio poroso. El modelo capturó correctamente la discontinuidad en el esfuerzo cortante y la continuidad de la velocidad, validando la capacidad del modelo para manejar medios inhomogéneos y transiciones de porosidad.
• Problema de Buckley-Leverett: Se simuló el desplazamiento de un fluido inmiscible por otro. El modelo reprodujo con precisión la formación de la onda de choque (frente de saturación) y su velocidad, comparándose favorablemente con la solución analítica clásica, a pesar de pequeñas oscilaciones numéricas típicas del método LBM.
• Dedos Viscosos (Viscous Fingering): Se estudió la inestabilidad interfacial en medios porosos. Los resultados demostraron que el modelo puede capturar la formación de dedos viscosos y, crucialmente, cómo la mojabilidad afecta la dinámica:
  • En medios hidrofílicos (ángulo de contacto bajo), los dedos se dispersan y se inhibe su crecimiento.
  • En medios hidrofóbicos, los dedos son más pronunciados.
  • Esto confirma que el modelo integra correctamente los efectos de mojabilidad en la dinámica macroscópica.

5. Significado e Impacto

Este trabajo representa un avance teórico significativo en la mecánica de fluidos en medios porosos:

• Superación de Empirismo: Proporciona una alternativa fundamentada físicamente a las correlaciones empíricas tradicionales para la permeabilidad relativa y la presión capilar, derivándolas de primeros principios termodinámicos y mecánicos.
• Predicción de Fenómenos Complejos: Al incluir la mojabilidad de manera formal en el potencial químico, el modelo tiene el potencial de predecir comportamientos de flujo en condiciones donde las leyes empíricas fallan (e.g., cambios de mojabilidad, histéresis compleja).
• Herramienta para Ingeniería: Ofrece un marco robusto para simular procesos de recuperación de hidrocarburos, almacenamiento de CO₂ y flujo en reactores catalíticos, permitiendo una mejor optimización de procesos industriales al capturar la física interfacial real sin necesidad de resolver costosas simulaciones a escala de poro en dominios grandes.
• Base para Futuras Investigaciones: Establece una base para el desarrollo de modelos más precisos que puedan incorporar efectos de temperatura, tensión superficial variable y geometrías de poro más complejas mediante la refinación de los coeficientes de cierre.

En resumen, el artículo presenta un marco unificado y riguroso que cierra la brecha entre la dinámica microscópica de interfases y el transporte macroscópico, ofreciendo una herramienta teórica y numérica superior para el estudio de flujos multifásicos en medios porosos.