The zero capillarity limit for the Euler-Korteweg system with no-flux boundary conditions

Titel: Der Grenzwert verschwindender Kapillarität für das Euler-Korteweg-System mit No-Flux-Randbedingungen
Autoren: Paolo Antonelli und Yuri Cacchio

1. Problemstellung

Der Artikel untersucht das asymptotische Verhalten des Euler-Korteweg-Systems (EK) im Grenzwert verschwindender Kapillarität ($\varepsilon \to 0$). Das System beschreibt die Zeitentwicklung einer kompressiblen, reibungsfreien Strömung unter dem Einfluss kapillarer Kräfte (diffuse Grenzflächen) in einem Gebiet $\Omega \subset \mathbb{R}^d$ mit glattem Rand.

Die Gleichungen lauten:
$$\begin{cases} \partial_t \rho + \text{div}(\rho u) = 0 \\ \partial_t (\rho u) + \text{div}(\rho u \otimes u) + \nabla p(\rho) = \varepsilon^2 \rho \nabla \left( \text{div}(k(\rho)\nabla\rho) - \frac{1}{2}k'(\rho)|\nabla\rho|^2 \right) \end{cases}$$
wobei $\rho$ die Dichte, $u$ die Geschwindigkeit und $k(\rho)$ der Kapillarkoeffizient ist.

Randbedingungen:
Das System wird mit No-Flux-Randbedingungen (kein Durchfluss) betrachtet:
$$u \cdot n|_{\partial\Omega} = 0, \quad \partial_n \rho|_{\partial\Omega} = 0.$$

Ziel:
Die Konvergenz der schwachen Lösungen des EK-Systems (mit $\varepsilon > 0$) gegen die starken Lösungen des klassischen kompressiblen Euler-Systems (mit $\varepsilon = 0$) zu beweisen. Das limitierende Euler-System erfüllt jedoch nur die Bedingung $u_E \cdot n|_{\partial\Omega} = 0$.

Herausforderung:
Es besteht eine Inkompatibilität der Randbedingungen: Das EK-System erfordert $\partial_n \rho = 0$, während das limitierende Euler-System dies nicht garantiert. Dies führt zur Bildung einer Grenzschicht (boundary layer) nahe dem Rand, was die Konvergenzanalyse, insbesondere für den Dichtegradienten, erheblich erschwert.

2. Methodik

Die Autoren verwenden eine Relative-Energie-Methode (auch Relativ-Entropie-Methode), ein etabliertes Werkzeug zur Untersuchung der Stabilität und von Singulargrenzwerten bei PDEs.

Schlüsseltechniken:

1. Erweiterung des Relativ-Energie-Funktionals:
   Anstatt nur die Standard-Energie zu betrachten, führen die Autoren ein augmentiertes Relativ-Energie-Funktional $E_h$ ein. Dies geschieht durch die Interpretation der kapillaren Energie als kinetische Energie eines Hilfsgeschwindigkeitsfeldes $v(\rho)$.
   $$v(\rho) = \sqrt{\frac{k(\rho)}{\rho}} \nabla \rho = \nabla \theta(\rho)$$
   Damit wird die kapillare Energie $\frac{\varepsilon^2}{2} k(\rho)|\nabla\rho|^2$ als $\frac{\varepsilon^2}{2} \rho |v|^2$ interpretiert. Dies erlaubt die Behandlung des Dichtegradienten als eine Art "Geschwindigkeit" in einem erweiterten System.

2. Randlagen-Korrektur (Boundary Layer Correction):
   Da die Lösung des limitierenden Euler-Systems ($\rho_E, u_E$) im Allgemeinen die Neumann-Randbedingung $\partial_n \rho_E = 0$ (und somit $v(\rho_E) \cdot n = 0$) nicht erfüllt, kann sie nicht direkt als Testfunktion im Relativ-Energie-Funktional verwendet werden.
   Die Autoren konstruieren eine Randlagen-Korrektur $v_{bl}$, die in einer dünnen Schicht der Dicke $\delta$ nahe dem Rand definiert ist. Die korrigierte Testfunktion $v_{bl}^{corr} = v(\rho_E) - v_{bl}$ erfüllt dann exakt die erforderliche homogene Neumann-Bedingung.

3. Schwache Lösungen und Vakuum:
   Die Analyse erlaubt das Auftreten von Vakuum ($\rho = 0$). Daher werden die Variablen in Dichte $\rho$ und Impuls $J = \rho u$ (bzw. Pseudo-Geschwindigkeit $\Lambda = \sqrt{\rho}u$) formuliert, um Singularitäten zu vermeiden. Es werden finite-energy schwache Lösungen betrachtet.

4. Abschätzung der Restterme:
   Durch Einsetzen der korrigierten Testfunktionen in die relative Energie-Ungleichung entstehen Restterme, die aus der Diskrepanz zwischen der exakten Lösung und der korrigierten Testfunktion resultieren. Die Autoren zeigen, dass diese Terme (insbesondere die durch die Grenzschicht verursachten) kontrolliert werden können, ohne zusätzliche Annahmen über das Fehlen anomaler Energiekonzentrationen treffen zu müssen.

3. Hauptbeiträge und Ergebnisse

• Konvergenz der hydrodynamischen Variablen (Satz 2.6):
  Es wird gezeigt, dass die schwachen Lösungen des EK-Systems gegen die starken Lösungen des Euler-Systems in den natürlichen Topologien konvergieren:
  $$\|\rho^\varepsilon - \rho^E\|_{L^\infty(L^1 \cap L^\gamma)} + \|\Lambda^\varepsilon - \sqrt{\rho^\varepsilon} u^E\|_{L^\infty(L^2)} \to 0.$$
  Dieser Teil der Analyse ist unempfindlich gegenüber der Grenzschichtbildung, da er nur die Konvergenz der Dichte und des Impulses betrachtet.

• Konvergenz des Dichtegradienten (Satz 2.12):
  Dies ist das zentrale neue Ergebnis. Durch das augmentierte Funktional $E_h$ und die Randlagen-Korrektur wird die Konvergenz des Dichtegradienten bewiesen:
  $$\varepsilon \|\nabla \beta(\rho^\varepsilon) - \sqrt{\rho^\varepsilon} v^E\|_{L^\infty(L^2)} \to 0.$$
  Hier ist $\beta(\rho)$ eine Funktion, die mit der Kapillarkraft verknüpft ist ($\nabla \beta \sim \sqrt{k(\rho)} \nabla \rho$).

• Vergleich mit dem viskosen Grenzwert:
  Im Gegensatz zum Grenzwert verschwindender Viskosität (Navier-Stokes zu Euler), bei dem starke Annahmen über das Fehlen anomaler Dissipation an der Grenze notwendig sind, ist die Grenzschicht im kapillaren Fall schwächer. Die relative Energie-Methode absorbiert die Grenzschichtfehler automatisch, sodass keine zusätzlichen Bedingungen an die Kapillar-Energie am Rand erforderlich sind.

• Allgemeinheit:
  Die Ergebnisse gelten für eine allgemeine Klasse von Kapillarkoeffizienten $k(\rho) = c_\alpha \rho^\alpha$ und für Gebiete mit nicht-trivialen Rändern (beschränkt oder unbeschränkt), nicht nur für periodische Gebiete oder den gesamten Raum.

4. Bedeutung und Fazit

Dieser Artikel liefert eine rigorose mathematische Begründung für den Übergang von Modellen mit diffusen Grenzflächen (Euler-Korteweg) zu Modellen mit scharfen Grenzflächen (klassisches Euler-System) in Anwesenheit von Wänden.

• Methodischer Fortschritt: Die Einführung der Randlagen-Korrektur im Kontext der relativen Energie für kapillare Flüssigkeiten ist ein innovativer Ansatz, der die Lücke zwischen der physikalischen Erwartung einer Grenzschicht und der mathematischen Konvergenzanalyse schließt.
• Physikalische Relevanz: Die Ergebnisse sind relevant für die Modellierung von Quantenflüssigkeiten (Bose-Einstein-Kondensate, Superfluid Helium), wo das EK-System als Quanten-Hydrodynamik (QHD) auftritt, sowie für Flüssigkeits-Dampf-Übergänge.
• Robustheit: Die Methode ist weniger restriktiv als frühere Arbeiten, da sie weniger Regularität der Anfangsdaten erfordert und keine speziellen Annahmen über die Energiekonzentration am Rand macht.

Zusammenfassend demonstriert das Papier, dass die relative Energie-Methode, erweitert um eine geschickte Behandlung der Randbedingungen durch Korrekturterme, ein leistungsfähiges Werkzeug ist, um singuläre Grenzwerte in kompressiblen Strömungsmodellen mit nicht-trivialen Randbedingungen zu analysieren.