Global weak solutions to a compressible Navier--Stokes/Cahn--Hilliard system with singular entropy of mixing

Résumé technique : Solutions faibles globales pour un système Navier–Stokes compressible / Cahn–Hilliard avec entropie de mélange singulière

Énoncé du problème
Ce travail traite de l'analyse mathématique d'un système Navier–Stokes/Cahn–Hilliard compressible modélisant l'évolution d'un mélange binaire de fluides visqueux subissant une séparation de phases. Le système est défini sur un domaine borné $\Omega \subset \mathbb{R}^3$ sur un intervalle de temps $(0, T)$. Les équations de conservation décrivent la conservation de la masse, de la quantité de mouvement et l'évolution de la variable de phase (différence de concentration de masse) couplée à un potentiel chimique.

La nouveauté et le défi principal de cette étude résident dans le choix de la densité d'énergie libre. Contrairement à la littérature existante sur les mélanges compressibles, qui utilise généralement des potentiels réguliers à croissance polynomiale, ce travail intègre l'entropie de mélange de Flory-Huggins (logarithmique), physiquement pertinente. Ce potentiel est singulier aux limites de l'intervalle physique $[-1, 1]$, imposant la contrainte que la variable de phase $c$ reste dans $(-1, 1)$. Le système est soumis à des conditions aux limites de non-glissement pour la vitesse et des conditions de Neumann homogènes pour la phase et le potentiel chimique.

Méthodologie
Les auteurs emploient une méthode de compacité basée sur un schéma de régularisation-approximation pour prouver l'existence de solutions faibles globales dans le temps pour des données initiales arbitrairement grandes. La méthodologie centrale procède comme suit :

1. Schéma d'approximation : Le potentiel singulier de Flory-Huggins $F(c)$ est approximé par une séquence de potentiels réguliers $F_\varepsilon(c)$ définis sur toute la droite réelle, présentant au plus une croissance quadratique à l'infini. Cela permet l'application des résultats d'existence existants pour les systèmes régularisés (spécifiquement ceux établis par Abels et Feireisl pour les mélanges compressibles avec des potentiels réguliers).
2. Estimations uniformes : Un composant critique de la preuve est la dérivation d'estimations uniformes pour les solutions approchées $(\varrho_\varepsilon, u_\varepsilon, c_\varepsilon)$ qui sont indépendantes du paramètre de régularisation $\varepsilon$.
   • Estimations d'énergie : Des inégalités d'énergie standards fournissent des bornes sur la densité, la vitesse et le gradient de la variable de phase.
   • Nouvelles estimations pour le Cahn-Hilliard dépendant de la densité : Les auteurs dérivent de nouvelles estimations pour le potentiel chimique $\mu_\varepsilon$ et le terme $\varrho_\varepsilon F'_\varepsilon(c_\varepsilon)$ sous des hypothèses minimales sur la densité (non négative, $\gamma$-intégrable avec $\gamma > 3/2$). En exploitant le contrôle de $\nabla_x \mu_\varepsilon$ à partir du bilan d'énergie et en utilisant une fonction test spécifique impliquant la différence de masse totale, ils établissent des bornes uniformes en $L^2$ sur $\int_\Omega \varrho_\varepsilon |F'_\varepsilon(c_\varepsilon)| dx$.
   • Estimations de pression : Pour gérer la compressibilité, les auteurs adaptent la théorie de Lions-Feireisl. Ils dérivent des estimations d'intégrabilité améliorées pour la pression $p(\varrho_\varepsilon)$ en testant l'équation de quantité de mouvement avec une fonction test spécifique basée sur l'opérateur de Bogovskii. Cela produit des bornes $L^q$ sur la pression, où $q$ dépend de l'exposant adiabatique $\gamma$.
3. Compacité et passage à la limite : En utilisant les estimations uniformes, les auteurs extraient des sous-suites convergentes. La convergence forte de la densité $\varrho_\varepsilon$ est établie via l'équation de continuité renormalisée et l'argument du flux visqueux effectif. La convergence forte de la variable de phase $c_\varepsilon$ est obtenue en adaptant les arguments de compacité d'études antérieures sur les cas incompressibles, en tirant parti de la convergence forte de la densité et des bornes uniformes sur le potentiel chimique.

Contributions clés et résultats
Le résultat principal, énoncé dans le Théorème 1.1, établit l'existence de solutions faibles globales dans le temps pour $\gamma > 3/2$ dans des domaines bornés en trois dimensions. La solution satisfait les propriétés suivantes :

• Régularité : La solution $(\varrho, u, c, \mu)$ possède des classes de régularité spécifiques, incluant $\varrho \in C_{weak}([0, T]; L^\gamma(\Omega))$ et $c \in L^\infty(0, T; W^{1,2}(\Omega)) \cap L^2(0, T; W^{2, \frac{2\gamma}{\gamma+1}}(\Omega))$.
• Contraintes physiques : Un accomplissement significatif est la preuve que la variable de phase $c$ respecte les bornes physiques presque partout. Plus précisément :
  • $-1 \leq c \leq 1$ presque partout dans $(0, T) \times \Omega$.
  • $-1 < c < 1$ presque partout sur l'ensemble où la densité est positive ($\{\varrho > 0\}$).
• Données initiales générales : Le résultat est valable pour des données initiales possédant une énergie totale finie, sans nécessiter que la densité soit strictement positive ou éloignée de zéro, ni exiger que la concentration initiale soit strictement comprise dans $(-1, 1)$ partout (seule la masse relative $M_c/M \in (-1, 1)$ est requise).

Signification et affirmations
L'article affirme fournir le premier résultat d'existence de solutions faibles globales pour le système Navier–Stokes/Cahn–Hilliard compressible avec le potentiel de Flory-Huggins physiquement pertinent. Les auteurs soulignent que les travaux précédents sur les mélanges compressibles étaient limités à des potentiels réguliers à croissance polynomiale, qui ne parviennent pas à imposer la contrainte physique $c \in [-1, 1]$.

La signification de ce travail réside dans :

1. Combler le fossé : Il traite le fossé analytique entre les modèles de champ de phase incompressibles et compressibles en gérant l'entropie singulière dans le régime compressible.
2. Hypothèses minimales : La dérivation des estimations pour l'équation de Cahn-Hilliard dépendant de la densité ne repose pas sur la positivité stricte ou la bornitude de la densité, une condition souvent requise dans les études précédentes (par exemple, dans le cas incompressible par le second auteur et Temam).
3. Cohérence physique : La preuve démontre rigoureusement que la solution limite respecte les bornes physiques de la variable de phase, assurant la cohérence du modèle avec la thermodynamique sous-jacente de la séparation de phases.

Les auteurs notent que leur approche repose sur une méthode de compacité et que la preuve de la convergence forte des densités suit des techniques établies (théorie de Lions-Feireisl) adaptées à ce système couplé spécifique. Ils mentionnent également que les résultats peuvent être étendus à des domaines en deux dimensions avec une régularité améliorée et que des potentiels singuliers plus généraux satisfaisant des conditions de croissance similaires peuvent être traités.