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

Technische Zusammenfassung: Existenz, Eindeutigkeit und Regularität schwacher Lösungen für ein thermodynamisch konsistentes Zweiphasenströmungsmodell

Problemstellung
Diese Arbeit befasst sich mit der mathematischen Analyse eines kürzlich vorgeschlagenen thermodynamisch konsistenten Modells für inkompressible und ungemischte Zweiphasenströmungen in porösen Medien, das ursprünglich von Kou et al. [25, 35] eingeführt wurde. Während bereits energiestabile numerische Methoden für dieses Modell entwickelt wurden, fehlte bisher eine fundierte mathematische Grundlage hinsichtlich der Existenz, Eindeutigkeit und Regularität seiner schwachen Lösungen. Das Modell charakterisiert Kapillareffekte durch eine logarithmische freie Energiefunktion, wodurch die Einhaltung des zweiten Hauptsatzes der Thermodynamik gewährleistet wird und eine elegante mathematische Symmetrie entsteht. Die primäre Herausforderung liegt in der hochgradig nichtlinearen Natur des Systems, insbesondere in der Kopplung zwischen Sättigung, chemischem Potenzial und Druck, ohne dabei auf künstliche Variablen zurückzugreifen, denen die physikalische Bedeutung fehlt.

Methodik
Die Autoren verwenden einen mehrstufigen analytischen Ansatz, um die Wohlgestelltheit des Modells zu etablieren:

1. Diskretisierung und Approximation: Die Studie beginnt mit der Konstruktion einer voll impliziten zeitsemi-diskreten Approximation und einer voll diskreten Approximation (unter Verwendung von Galerkin-Methoden) für die schwache Formulierung des Modells. Die primären Variablen sind das chemische Potenzial ($\mu_w$) und der reale Phasendruck ($p$), wobei die Verwendung künstlicher globaler oder komplementärer Drücke vermieden wird, die in der klassischen Analyse oft erforderlich sind.
2. Existenz mittels Fixpunktsätzen: Um die Existenz von Lösungen für das voll diskrete Schema zu beweisen, nutzen die Autoren den Satz über die „Nullstellen eines Vektorfeldes“ (eine Variante des Brouwerschen Fixpunktsatzes). Dies beinhaltet die Definition einer stetigen Abbildung und den Nachweis, dass diese spezifische Kohärenzbedingungen auf dem Rand einer ausreichend großen Kugel im endlichen dimensionalen Raum erfüllt.
3. Grenzwertbildung (Passage to the Limit): Die Existenz von Lösungen für die semi-diskreten und kontinuierlichen Modelle wird durch Techniken der schwachen Konvergenz abgeleitet. Entscheidend für diesen Schritt ist die Herleitung einheitlicher Energiestabilitätsabschätzungen und der Beschränktheit der Sättigung. Der Aubin-Lions-Simon-Lemma wird angewendet, um stark konvergente Teilfolgen zu extrahieren, wodurch der Grenzwert der diskreten Schemata die kontinuierliche schwache Formulierung erfüllt.
4. Analyse der Eindeutigkeit: Die Eindeutigkeit wird unter der Annahme etabliert, dass das chemische Potenzial im Raum Lipschitz-stetig ist (mit einer zeitabhängigen Lipschitz-Konstante in $L^2$). Der Beweis nutzt die Differenz zweier potenzieller Lösungen und wendet die Grönwall-Ungleichung an, um zu zeigen, dass die Differenz verschwinden muss.
5. Regularitätstheorie: Um eine höhere Regularität zu etablieren, führen die Autoren künstliche und komplementäre Drücke ein (folgend dem Rahmenwerk von Chen [19]), um das System in ein gekoppeltes elliptisch-parabolisches System zu transformieren. Durch die Nutzung der Regularitätstheorie für elliptische partielle Differentialgleichungen (PDEs) und unter der Annahme vollständiger Neumann-Randbedingungen leiten sie eine $H^2$-Regularität für den Druck und das chemische Potenzial her.

Wesentliche Beiträge und Ergebnisse

• Existenz schwacher Lösungen: Die Arbeit beweist die Existenz einer schwachen Lösung $(S_w, \mu_w, p)$ für das neue thermodynamisch konsistente Modell. Ein signifikantes Ergebnis ist, dass dieser Beweis keine Annahmen über die Beschränktheit höherer Ablekungsterme erfordert – eine häufige Einschränkung in klassischen Modellen, die oft den Einsatz nicht-physikalischer komplementärer Drücke notwendig macht.
• Eindeutigkeit: Die Autoren etablieren die Eindeutigkeit der schwachen Lösung. Im Gegensatz zu früheren Arbeiten, die die Lipschitz-Stetigkeit eines künstlichen globalen Drucks erforderten (was eine Hölder-Stetigkeit der Sättigung impliziert), setzt diese Arbeit die Lipschitz-Stetigkeit des chemischen Potenzials $\mu_w$ voraus.
• Regularität: Unter spezifischen Annahmen bezüglich der Domänengeometrie (konvex mit $C^{1,\beta}$-Rand), der Anfangsdaten und der Materialparameter stellt die Arbeit fest, dass die schwache Lösung eine höhere Regularität besitzt. Speziell gehört die Sättigung $S_w$ zu $L^\infty(0, T; H^1(\Omega)) \cap H^1(0, T; L^2(\Omega))$, während das chemische Potenzial $\mu_w$ und der Druck $p$ zu $L^2(0, T; H^2(\Omega)) \cap L^\infty(0, T; H^1(\Omega))$ gehören.
• Sättigungsgrenzen: Die Analyse demonstriert rigoros, dass die Sättigung innerhalb physikalisch sinnvoller Grenzen bleibt ($S_\epsilon \le S_w \le 1 - S_\epsilon$), indem die Monotonie des chemischen Potenzials sowie spezifische Bedingungen auf die Randdaten und Quellterme genutzt werden.

Bedeutung
Die Arbeit beansprucht, die erste Arbeit zu sein, die die Wohlgestelltheit (Existenz, Eindeutigkeit und Regularität) des spezifischen thermodynamisch konsistenten Zweiphasenströmungsmodells aus [25, 35] rigoros behandelt. Durch die Vermeidung künstlicher Variablen und die Verwendung des physikalischen chemischen Potenzials als primäre Variable bietet die Analyse eine direktere mathematische Rechtfertigung für die Struktur des Modells. Die Ergebnisse bestätigen, dass das Modell mathematisch wohlgestellt ist, was die zuvor entwickelten energiestabilen numerischen Methoden validiert und eine theoretische Basis für zukünftige Erweiterungen unter Berücksichtigung von Gesteinskompressibilität oder poro-viskoelastischen Medien schafft.