The Existence, uniqueness, and regularity of weak solutions for a thermodynamically consistent two-phase flow model in porous media

Résumé Technique : Existence, Unicité et Régularité des Solutions Faibles pour un Modèle de Flux Bifasique Thermodynamiquement Cohérent

Énoncé du Problème
Cet article traite de l'analyse mathématique d'un modèle récemment proposé pour le flux bifasique incompressible et immiscible en milieu poreux, qui est thermodynamiquement cohérent, introduit initialement par Kou et al. [25, 35]. Bien que des méthodes numériques stables sur le plan de l'énergie aient été développées pour ce modèle, une base mathématique rigoureuse concernant l'existence, l'unicité et la régularité de ses solutions faibles faisait défaut jusqu'à présent. Le modèle caractérise les effets de capillarité par une fonction d'énergie libre logarithmique, assurant le respect du second principe de la thermodynamique et possédant une élégante symétrie mathématique. Le défi principal réside dans la nature hautement non linéaire du système, particulièrement le couplage entre la saturation, le potentiel chimique et la pression, sans recourir à des variables artificielles dépourvues de signification physique.

Méthodologie
Les auteurs emploient une approche analytique en plusieurs étapes pour établir la bien-poséité du modèle :

1. Discrétisation et Approximation : L'étude commence par la construction d'une approximation semi-discrète temporelle entièrement implicite et d'une approximation entièrement discrète (utilisant des méthodes de Galerkin) pour la formulation faible du modèle. Les variables primaires choisies sont le potentiel chimique ($\mu_w$) et la pression réelle de la phase ($p$), évitant ainsi l'utilisation de pressions globales ou complémentaires artificielles souvent requises dans les analyses classiques.
2. Existence via les Théorèmes de Point Fixe : Pour prouver l'existence de solutions pour le schéma entièrement discret, les auteurs utilisent le théorème des « zéros d'un champ de vecteurs » (une variante du théorème du point fixe de Brouwer). Cela implique de définir une application continue et de démontrer qu'elle satisfait des conditions de coercivité spécifiques sur le bord d'une boule suffisamment grande dans l'espace de dimension finie.
3. Passage à la Limite : L'existence de solutions pour les modèles semi-discrets et continus est dérivée par des techniques de convergence faible. La clé de cette étape est la dérivation d'estimations de stabilité énergétique uniformes et le caractère borné de la saturation. Le lemme d'Aubin-Lions-Simon est appliqué pour extraire des sous-suites de convergence forte, permettant à la limite des schémas discrets de satisfaire la formulation faible continue.
4. Analyse de l'Unicité : L'unicité est établie sous l'hypothèse que le potentiel chimique est lipschitzien dans l'espace (avec une constante de Lipschitz dépendante du temps dans $L^2$). La preuve utilise la différence de deux solutions potentielles, appliquant l'inégalité de Grönwall pour montrer que la différence doit s'annuler.
5. Théorie de la Régularité : Pour établir une régularité supérieure, les auteurs introduisent des pressions artificielles et complémentaires (suivant le cadre de Chen [19]) afin de transformer le système en un système couplé elliptique-parabolique. En exploitant la théorie de la régularité pour les équations aux dérivées partielles (EDP) elliptiques et en supposée des conditions de Neumann complètes, ils dérivent une régularité $H^2$ pour la pression et le potentiel chimique.

Contributions Clés et Résultats

• Existence de Solutions Faibles : L'article prouve l'existence d'une solution faible $(S_w, \mu_w, p)$ pour le nouveau modèle thermodynamiquement cohérent. Un résultat significatif est que cette preuve ne nécessite pas d'hypothèses sur le caractère borné des termes de dérivées d'ordre supérieur, une contrainte courante dans les modèles classiques qui nécessite souvent l'usage de pressions complémentaires non physiques.
• Unicité : Les auteurs établissent l'unicité de la solution faible. Contrairement aux travaux précédents qui nécessitaient la continuité de Lipschitz d'une pression globale artificielle (ce qui implique la continuité de Hölder de la saturation), ce travail suppose la continuité de Lipschitz du potentiel chimique $\mu_w$.
• Régularité : Sous des hypothèses spécifiques concernant la géométrie du domaine (convexe avec un bord $C^{1,\beta}$), les données initiales et les paramètres matériels, l'article établit que la solution faible possède une régularité supérieure. Plus précisément, la saturation $S_w$ appartient à $L^\infty(0, T; H^1(\Omega)) \cap H^1(0, T; L^2(\Omega))$, tandis que le potentiel chimique $\mu_w$ et la pression $p$ appartiennent à $L^2(0, T; H^2(\Omega)) \cap L^\infty(0, T; H^1(\Omega))$.
• Bornes de Saturation : L'analyse démontre rigoureusement que la saturation reste dans des limites physiquement significatives ($S_\epsilon \le S_w \le 1 - S_\epsilon$) en utilisant la monotonicité du potentiel chimique et des conditions spécifiques sur les données aux limites et les termes sources.

Signification
L'article affirme être le premier travail à traiter rigoureusement la bien-poséité (existence, unicité et régularité) de ce modèle spécifique de flux bifasique thermodynamiquement cohérent proposé dans [25, 35]. En évitant les variables artificielles et en s'appuyant sur le potentiel chimique physique comme variable primaire, l'analyse fournit une justification mathématique plus directe de la structure du modèle. Les résultats confirment que le modèle est mathématiquement bien posé, validant ainsi les méthodes numériques stables sur le plan de l'énergie précédemment développées pour celui-ci, et fournissant une base théorique pour de futures extensions impliquant la compressibilité de la roche ou les milieux poro-viscoélastiques.