Weak-strong uniqueness for bi-fluid compressible system with algebraic closure

A continuación presento un resumen técnico detallado del artículo "Weak–strong uniqueness for bi-fluid compressible system with algebraic closure" (Unicidad débil-fuerte para un sistema de bi-fluidos compresibles con cierre algebraico), escrito por Yang Li, Mária Lukáčová-Medvid'ová, Milan Pokorný y Ewelina Zatorska.

1. Planteamiento del Problema

El artículo aborda un sistema de ecuaciones de fluidos compresibles viscosos de dos fases (bi-fluidos) en un dominio acotado tridimensional $\Omega$. El modelo se caracteriza por:

• Campo de velocidad común: Ambas fases comparten la misma velocidad $u$.
• Cierre algebraico: Las fases están acopladas mediante una relación constitutiva implícita, específicamente la igualdad de presiones $p_+(\rho_+) = p_-(\rho_-)$, donde las presiones siguen leyes barotrópicas de potencia ($p \sim \rho^\gamma$).
• Variables: Las incógnitas son las fracciones volumétricas $\alpha$ (y $1-\alpha$), las densidades parciales $\rho_+$ y $\rho_-$, y la velocidad $u$.

El desafío principal:
A diferencia de los modelos de cierre diferencial (donde la fracción volumétrica $\alpha$ obedece una ecuación de transporte pura), el cierre algebraico implica que $\alpha$ satisface una ecuación de transporte con un término no lineal adicional dependiente de la divergencia de la velocidad ($\text{div } u$). Esta estructura no lineal rompe las estimaciones estándar de entropía relativa utilizadas en modelos más simples, generando términos "problemáticos" que no se cancelan fácilmente.

El objetivo es probar el principio de unicidad débil-fuerte: demostrar que cualquier solución débil de energía finita (existente globalmente, según Novotný y Pokorný) coincide con la solución fuerte local (existente según Piasecki y Zatorska) siempre que ambas partan de los mismos datos iniciales, durante el tiempo de vida de la solución fuerte.

2. Metodología

Los autores emplean el método de entropía relativa, una técnica robusta para estudiar la estabilidad y unicidad en sistemas de ecuaciones en derivadas parciales hiperbólicas y parabolicas.

A. Funcional de Entropía Relativa

Definen un funcional de energía relativa $\mathcal{E}$ que mide la distancia entre la solución débil $(\alpha, \rho_+, \rho_-, u)$ y la solución fuerte de referencia $(\beta, \tilde{\rho}_+, \tilde{\rho}_-, v)$:
$$\mathcal{E} = \int_\Omega \left[ \frac{1}{2}(R+Q)|u-v|^2 + \alpha \mathcal{H}_+(\rho_+|\tilde{\rho}_+) + (1-\alpha)\mathcal{H}_-(\rho_-|\tilde{\rho}_-) + \frac{1}{2}(\alpha-\beta)^2 \right] dx$$
donde $R=\alpha\rho_+$ y $Q=(1-\alpha)\rho_-$ son las masas parciales, y $\mathcal{H}$ representa la entropía termodinámica (potencial de Helmholtz).

B. Derivación de la Desigualdad

1. Ecuaciones para $\alpha$ y densidades: Derivan rigurosamente las ecuaciones de evolución para $\alpha$, $\rho_+$ y $\rho_-$ a partir de las ecuaciones de continuidad y la relación algebraica de presiones. Se obtiene una ecuación de transporte para $\alpha$ con un coeficiente no lineal $\omega_\alpha$ que depende de $\alpha$ y $\gamma$.
2. Descomposición del residuo: Al aplicar la desigualdad de entropía relativa, surgen términos residuales ($J_1, \dots, J_8$) que deben ser acotados.
3. Tratamiento de los términos problemáticos:
   • La mayor dificultad radica en los términos que involucran la diferencia de fracciones volumétricas $(\alpha - \beta)$ y la divergencia de la velocidad.
   • Los autores introducen un mecanismo adicional para controlar el defecto en la fracción volumétrica. Aprovechan la estructura implícita específica de la ley de presiones para derivar una estimación de estabilidad para las fracciones volumétricas que es compatible con la desigualdad de energía relativa.
   • Utilizan una descomposición del dominio en conjuntos "esenciales" (donde las densidades están acotadas) y "residuales" (donde pueden ser muy grandes o pequeñas) para manejar la no linealidad de las presiones.

C. Argumento de Gronwall

Combinan la estimación de la entropía relativa con una estimación auxiliar para la diferencia de las fracciones volumétricas $(\alpha - \beta)^2$ (Lema 5.2). Al sumar ambas desigualdades y aplicar el lema de Gronwall, demuestran que si la entropía inicial es cero (mismos datos), la entropía permanece cero en todo el intervalo de tiempo.

3. Contribuciones Clave

1. Resolución de la obstrucción del cierre algebraico: A diferencia de trabajos anteriores (como Jin y Novotný, 2020) que trataban el cierre diferencial, este trabajo resuelve la complejidad añadida por el cierre algebraico ($p_+ = p_-$). El término no lineal en la ecuación de transporte de $\alpha$ genera nuevos términos en la desigualdad de entropía que requieren un tratamiento específico.
2. Eliminación de restricciones de acotamiento: Una contribución significativa es que no se requiere que la fracción volumétrica $\alpha$ esté acotada estrictamente lejos de 0 y 1 (es decir, no se asume $0 < c \leq \alpha \leq 1-c < 1$). El análisis es válido incluso cuando $\alpha$ puede acercarse a los extremos, lo cual es físicamente más realista y matemáticamente más desafiante.
3. Unicidad sin convexidad estricta: El método evita depender de suposiciones de convexidad estricta del potencial de Helmholtz libre respecto a las variables conservadas, un requisito que en trabajos previos limitaba la aplicabilidad del resultado.
4. Generalización del método: Demuestran que el método de entropía relativa puede adaptarse exitosamente a sistemas de bi-fluidos con relaciones de cierre algebraico, abriendo la puerta a estudios de límites singulares (como el límite de bajo número de Mach o límites incompresibles) para este tipo de modelos.

4. Resultados Principales

• Teorema 2.1 (Unicidad Débil-Fuerte): Se establece que la solución débil global de energía finita (construida en [35]) coincide con la solución fuerte local (construida en [36]) en el intervalo de tiempo donde existe la solución fuerte, siempre que compartan los mismos datos iniciales.
• Estabilidad: Se demuestra que la solución fuerte es estable frente a perturbaciones en el sentido de la entropía relativa. Si los datos iniciales son cercanos, las soluciones permanecen cercanas.
• Estimaciones de Residuos: Se logra cerrar todas las estimaciones de los términos residuales ($J_1$ a $J_8$) utilizando la estructura algebraica del sistema y las propiedades de regularidad de la solución fuerte.

5. Significado e Impacto

Este trabajo es fundamental en el análisis matemático de fluidos multifásicos por varias razones:

• Rigor Matemático: Cierra una brecha importante en la teoría de existencia y unicidad para modelos de bi-fluidos compresibles con cierre algebraico, un modelo ampliamente utilizado en ingeniería (flujo de petróleo-gas, reactores nucleares) pero difícil de analizar matemáticamente debido a la no linealidad implícita.
• Validación de Modelos: Al probar la unicidad débil-fuerte, se valida que las soluciones numéricas o asintóticas que convergen a una solución fuerte son únicas, lo que es crucial para la confiabilidad de las simulaciones computacionales.
• Herramienta para Límites Singulares: La desigualdad de entropía relativa derivada proporciona un mecanismo de estabilidad robusto que permite justificar rigurosamente límites formales (como la incompresibilidad) para sistemas de dos fluidos, algo que anteriormente era un problema abierto para este tipo de cierres.
• Avance en Teoría de EDPs: El tratamiento de la ecuación de transporte con término fuente no lineal acoplado a la divergencia de la velocidad ofrece nuevas técnicas que pueden ser aplicadas a otros sistemas de fluidos multifásicos o reactivos.

En resumen, el artículo logra un hito técnico al superar las dificultades analíticas impuestas por el cierre algebraico en sistemas de dos fluidos, estableciendo un principio de unicidad robusto sin restricciones artificiales en las fracciones volumétricas.