On weak solutions for the stationary Cahn-Hillard-Navier-Stokes equations with singular potential

Resumen Técnico: Soluciones Débiles para las Ecuaciones de Cahn-Hilliard-Navier-Stokes Estacionarias con Potencial Singular

Planteamiento del Problema
Este artículo investiga la existencia de soluciones débiles para el sistema estacionario de Navier–Stokes–Cahn-Hilliard que describe una mezcla fluida bifásica compresible tridimensional con una interfaz difusa. El sistema está gobernado por la conservación de la masa, el momento y la concentración, acoplados con una ecuación de potencial químico. Una característica distintiva del modelo es la inclusión de una densidad de energía libre singular de tipo Flory–Huggins (logarítmica). Esta elección asegura que la fracción másica $c$ permanezca dentro del rango físico $(-1, 1)$ y permite la presencia de estados de vacío ($\rho = 0$).

El problema de valores en la frontera específico involucra un dominio acotado $\Omega \subset \mathbb{R}^3$ con condiciones de Dirichlet para la velocidad y condiciones de Neumann para la concentración y el potencial químico. Los principales desafíos matemáticos surgen de dos fuentes:

1. El Marco Estacionario: A diferencia de los problemas evolutivos, el sistema estacionario carece de una desigualdad de energía a priori ($E(t) \leq E(0)$) impulsada por la evolución temporal, la cual típicamente proporciona el control necesario sobre el potencial singular y el campo de concentración.
2. Potencial Singular y Vacío: La singularidad logarítmica en la energía libre impide las regularizaciones estándar mediante expansión de Taylor, y la posible degeneración de la densidad cerca de los estados de vacío complica los argumentos de compacidad requeridos para pasar al límite.

Metodología
Los autores establecen la existencia mediante un proceso de límite de dos niveles aplicado a un sistema aproximado cuidadosamente construido.

1. Regularización Especializada del Potencial:
   Para abordar la falta de cotas $L^2$ en $c$ y el fracaso de las expansiones de Taylor estándar, los autores introducen una regularización especializada $f^\delta_2(c)$ para el término logarítmico singular. A diferencia de enfoques anteriores, esta aproximación se construye de modo que la diferencia $f^\delta_2(c) - \frac{\theta_c}{2}c^2$ crezca solo linealmente cuando $|c| \to \infty$. Esta construcción específica está diseñada para cancelar el crecimiento cuadrático inducido por el término anti-difusivo, restaurando así la compacidad necesaria y permitiendo estimaciones uniformes.

2. Esquema de Aproximación de Dos Niveles:

   • Nivel 1 (Límite $\epsilon$): Los autores introducen primero una viscosidad artificial y parámetros de regularización ($\epsilon$) para construir un sistema regularizado. Utilizando un argumento de punto fijo (Teorema del Punto Fijo de Schaefer), prueban la existencia de soluciones fuertes para este sistema aproximado. Se derivan estimaciones uniformes independientes de $\epsilon$, lo que permite tomar el límite $\epsilon \to 0$, obteniendo una solución débil para un sistema intermedio con la energía libre regularizada.
   • Nivel 2 (Límite $\delta$): El segundo paso implica eliminar la regularización del potencial ($\delta \to 0$). Una dificultad mayor aquí es controlar la integrabilidad $L^{3/2}$ de la densidad. Los autores emplean una elección específica de presión artificial y un argumento de interpolación delicado. Este argumento explota la estructura del potencial regularizado, específicamente la observación de que $\frac{\partial f^\delta_2}{\partial c} - \theta_c c = O(\ln \frac{1}{\delta})$ para $|c| > 1+\delta$. Esto les permite cerrar las estimaciones para la norma de la densidad y establecer cotas uniformes independientes de $\delta$.

3. Convergencia Fuerte y Cotas Físicas:
   Pasar al límite $\delta \to 0$ requiere probar la convergencia fuerte de la densidad para manejar el término de presión no monótono. Adaptando técnicas de la teoría de las ecuaciones de Navier–Stokes compresibles estacionarias, los autores establecen la convergencia fuerte de la densidad. Además, prueban que la concentración límite $c$ satisface las estrictas cotas físicas $-1 < c < 1$ casi en todas partes en el soporte de la densidad ($\{\rho > 0\}$), asegurando que el término singular esté bien definido.

Resultados Clave
El teorema principal (Teorema 1.1) establece la existencia de una solución débil $(\rho, u, \mu, c)$ para el problema estacionario bajo las siguientes condiciones:

• Dominio: Un dominio acotado y simplemente conexo con frontera $C^2$.
• Fuerzas: Las fuerzas externas $g_1, g_2$ están acotadas.
• Exponente Adiabático ($\gamma$):
  • $\gamma > 3/2$ si la fuerza externa es irrotacional ($\nabla \times g_1 = 0$).
  • $\gamma > 5/3$ si la fuerza externa es rotacional ($\nabla \times g_1 \neq 0$).
• Restricciones de Masa: La masa total $m_1 > 0$ y la masa relativa $m_2$ satisfacen $m_2 \in (-m_1, m_1)$.

La solución satisface la ecuación de continuidad renormalizada, la formulación débil de las ecuaciones de momento y concentración, y la desigualdad de energía. Crucialmente, la concentración satisface $-1 < c(x) < 1$ casi en todas partes donde $\rho(x) > 0$.

Significado y Afirmaciones
Los autores afirman que este trabajo proporciona el primer resultado de existencia para soluciones débiles del sistema estacionario de Navier–Stokes–Cahn-Hilliard compresible que incorpora simultáneamente:

1. Un potencial de energía libre singular (tipo Flory–Huggins), que es físicamente consistente para mezclas binarias pero matemáticamente desafiante debido a la falta de cotas globales en la concentración.
2. La presencia de estados de vacío ($\rho = 0$), lo que introduce degeneración en las ecuaciones.

El artículo señala que, aunque el rango de exponentes adiabáticos admisibles $\gamma$ podría ser extensible basándose en resultados para las ecuaciones de Navier–Stokes compresibles estacionarias de un solo fluido, el rango actual está limitado por la dificultad de obtener estimaciones ponderadas para la concentración cerca de la frontera debido a la singularidad de la energía libre. El trabajo representa un paso significativo para cerrar la brecha entre los modelos evolutivos (donde las desigualdades de energía proporcionan control) y los modelos estacionarios con potenciales singulares.