Globale Wohlgestelltheit starker Lösungen für ein Bulk-Surface Navier-Stokes-Cahn-Hilliard-Modell mit nicht-degenerierten Mobilitäten in zwei Dimensionen

Grundinformationen

• Paper-ID: 2511.06847
• Titel: Global well-posedness of strong solutions to a bulk-surface Navier-Stokes-Cahn-Hilliard model with non-degenerate mobilities in two dimensions
• Autor: Jonas Stange (Universität Regensburg, Deutschland)
• Klassifizierung: math.AP (Analysis von PDEs), math-ph (Mathematische Physik), math.MP (Mathematische Physik)
• Einreichungsdatum: 10. November 2025
• Paper-Link: https://arxiv.org/abs/2511.06847v1

Zusammenfassung

Die vorliegende Arbeit untersucht ein thermodynamisch konsistentes Diffuse-Interface-Modell zur Beschreibung von Bulk-Surface-Viskosflüssigkeitsmischungen. Das Modell besteht aus der Navier-Stokes-Cahn-Hilliard-Gleichung im Bulkbereich, gekoppelt mit einem Surface-Navier-Stokes-Cahn-Hilliard-System auf dem Rand. Der Autor beweist die globale Wohlgestelltheit starker Lösungen im zweidimensionalen Fall für physikalisch bedeutsame nicht-degenerierte Mobilitätsfunktionen. Darüber hinaus werden Eindeutigkeit starker Lösungen sowie stetige Abhängigkeit von Anfangsdaten nachgewiesen. Der Ansatz basiert auf neuer Wohlgestelltheits- und Regularitätstheorie für konvektive Bulk-Surface-Cahn-Hilliard-Gleichungen mit nicht-degenerierten Mobilitäten sowie für Bulk-Surface-Stokes-Gleichungen mit nicht-konstanten Koeffizienten.

Forschungshintergrund und Motivation

Problemhintergrund

1. Modellierung physikalischer Phänomene: Die Untersuchung betrifft Diffuse-Interface-Modelle für zweiphasige Flüssigkeitsmischungen mit unterschiedlichen Dichten, wobei viskose Dynamikeffekte sowohl im Bulkbereich als auch auf dem Rand auftreten. Solche Modelle finden wichtige Anwendungen in der Materialwissenschaft und Biodynamik von Membranen.
2. Mathematische Herausforderungen: Das Modell koppelt:
   • Navier-Stokes-Cahn-Hilliard-Gleichungen im Bulkbereich (beschreiben Flüssigkeitsbewegung und Phasenseparation)
   • Surface-Navier-Stokes-Cahn-Hilliard-Systeme auf dem Rand (beschreiben Oberflächendynamik)
   • Dynamische Randbedingungen (ermöglichen Massenaustausch zwischen Bulk und Surface)
3. Limitierungen bestehender Theorien:
   • Für Diffuse-Interface-Modelle mit dynamischen Randbedingungen konzentriert sich die Literatur hauptsächlich auf Existenz schwacher Lösungen
   • Existenzresultate für starke Lösungen sind selten, besonders für nicht-degenerierte Mobilitätsfunktionen
   • Die Theorie globaler starker Lösungen im dreidimensionalen Fall ist noch unvollständig

Forschungsmotivation

1. Physikalische Plausibilität: Nicht-degenerierte Mobilitätsfunktionen (d.h. Mobilitäten mit positiven oberen und unteren Schranken im gesamten Bereich) entsprechen der physikalischen Realität besser
2. Mathematische Strenge: Starke Lösungen liefern feinere Regularitätsabschätzungen, die für das Verständnis des Langzeitverhaltens von Lösungen entscheidend sind
3. Theoretische Vervollständigung: Schließung von Lücken in der Theorie starker Lösungen für Navier-Stokes-Cahn-Hilliard-Systeme mit dynamischen Randbedingungen

Kernbeiträge

Die Hauptbeiträge dieser Arbeit umfassen:

1. Existenz globaler starker Lösungen: Im zweidimensionalen Fall wird die Existenz globaler starker Lösungen für das Bulk-Surface-Navier-Stokes-Cahn-Hilliard-System mit nicht-degenerierten Mobilitäten nachgewiesen (Theorem 3.1)
2. Eindeutigkeit und stetige Abhängigkeit: Eindeutigkeit starker Lösungen sowie stetige Abhängigkeit von Anfangsdaten werden bewiesen (Theorem 3.4), erfordern jedoch die Bedingung $L \in (0,\infty]$
3. Neue Regularitätstheorie: Zwei wesentliche Hilfstheorien werden etabliert:
   • Wohlgestelltheits- und Regularitätstheorie für konvektive Bulk-Surface-Cahn-Hilliard-Gleichungen mit nicht-degenerierten Mobilitäten (Anhang A)
   • Regularitätstheorie für Bulk-Surface-Stokes-Gleichungen mit nicht-konstanten Koeffizienten (Anhang B)
4. Energieungleichungen: Es wird nachgewiesen, dass Lösungen die Energiedissipationsgesetze erfüllen: $\frac{d}{dt}E_{tot}(v,w,\phi,\psi) + \text{Dissipationsterme} = 0$
5. Separationseigenschaft: Phasenfeldvariablen bleiben strikt im physikalischen Intervall $(-1,1)$ enthalten: $|\phi| < 1$ f.ü. in $Q$, $|\psi| < 1$ f.ü. auf $\Sigma$

Methodische Details

Aufgabendefinition

Betrachtet wird ein zweidimensionales beschränktes Gebiet $\Omega \subset \mathbb{R}^2$ mit Rand $\Gamma = \partial\Omega$. Setze $Q = \Omega \times (0,\infty)$, $\Sigma = \Gamma \times (0,\infty)$.

Kontrollgleichungssystem (Gleichung 1.1) umfasst:

1. Bulk-Impulsgleichung: $\partial_t(\rho(\phi)v) + \text{div}(v \otimes (\rho(\phi)v + J)) = \text{div} T, \quad \text{div} v = 0 \text{ in } Q$
2. Surface-Impulsgleichung: $\partial_t(\sigma(\psi)w) + \text{div}_\Gamma(w \otimes (\sigma(\psi)w + K)) = \text{div}_\Gamma T_\Gamma + Z, \quad \text{div}_\Gamma w = 0 \text{ on } \Sigma$
3. Bulk-Cahn-Hilliard-Gleichung: $\partial_t\phi + \text{div}(\phi v) = \text{div}(m_\Omega(\phi)\nabla\mu), \quad \mu = -\Delta\phi + F'(\phi) \text{ in } Q$
4. Surface-Cahn-Hilliard-Gleichung: $\partial_t\psi + \text{div}_\Gamma(\psi w) = \text{div}_\Gamma(m_\Gamma(\psi)\nabla_\Gamma\theta) - \beta m_\Omega(\phi)\partial_n\mu \text{ on } \Sigma$$\theta = -\Delta_\Gamma\psi + G'(\psi) + \alpha\partial_n\phi \text{ on } \Sigma$
5. Randbedingungen: $w = v|_\Gamma, \quad v \cdot n = 0 \text{ on } \Sigma$$K\partial_n\phi = \alpha\psi - \phi, \quad Lm_\Omega(\phi)\partial_n\mu = \beta\theta - \mu \text{ on } \Sigma$

Dabei:

• $v, w$ sind Bulk- und Surface-Geschwindigkeitsfelder
• $\phi, \psi$ sind Bulk- und Surface-Phasenfeldvariablen
• $\mu, \theta$ sind chemische Potentiale
• $\rho(\phi), \sigma(\psi)$ sind phasenfeldabhängige Dichtenfunktionen
• $m_\Omega, m_\Gamma$ sind Mobilitätsfunktionen (erfüllen $0 < m_* \leq m_\Omega, m_\Gamma \leq m^*$)

Modellarchitektur

Beweistrategie verwendet die Semi-Galerkin-Approximationsmethode:

Schritt 1: Konstruktion des Approximationsproblems

1. Einführung des Eigenfunktionssystems $\{(\tilde{v}_j, \tilde{w}_j)\}_{j\in\mathbb{N}}$ des Bulk-Surface-Stokes-Operators $A$
2. Definition endlichdimensionaler Unterräume $V_k = \text{span}\{(\tilde{v}_1, \tilde{w}_1), \ldots, (\tilde{v}_k, \tilde{w}_k)\}$
3. Konstruktion von Näherungslösungen $(v_k, w_k, \phi_k, \psi_k, \mu_k, \theta_k)$ erfüllend die Galerkin-Form der schwachen Gleichungen

Schritt 2: Fixpunktargument

Für festes Geschwindigkeitsfeld $(v_*, w_*)$ wird gelöst:

• Die konvektive Bulk-Surface-Cahn-Hilliard-Gleichung (4.8), ergebend $(\phi_k, \psi_k, \mu_k, \theta_k)$
• Die Galerkin-Approximation der Navier-Stokes-Gleichung (4.12), ergebend $(v_k, w_k)$

Durch den Schauder-Fixpunktsatz wird die Existenz eines Fixpunkts der Abbildung $\Lambda: (v_*, w_*) \mapsto (v_k, w_k)$ nachgewiesen.

Schritt 3: A-priori-Abschätzungen

Etablierung von Schlüsselabschätzungen (unabhängig von $k$):

1. Energieabschätzungen (Gleichungen 4.35-4.37): $\|(v_k, w_k)\|_{L^\infty(0,\infty; L^2)} + \|(v_k, w_k)\|_{L^2(0,\infty; H^1)} \leq C$$\|(\phi_k, \psi_k)\|_{L^\infty(0,\infty; H^1)} \leq C$
2. Höherordnungsabschätzungen (Gleichungen 4.43-4.45): $\|(\phi_k, \psi_k)\|_{L^\infty(0,\infty; W^{2,p})} + \|(F'(\phi_k), G'(\psi_k))\|_{L^\infty(0,\infty; L^p)} \leq K_{1,p}$$\|(\mu_k, \theta_k)\|_{L^\infty(0,\infty; H^1)} + \|(\mu_k, \theta_k)\|_{L^2_{uloc}(0,\infty; H^3)} \leq K_2$
3. Höherordnungsabschätzungen für Geschwindigkeitsfelder (Gleichung 4.61): $\sup_{t\geq 0} \int_t^{t+1} \|( v_k, w_k)\|_{H^2}^2 + \|(\partial_t v_k, \partial_t w_k)\|_{L^2}^2 + \|(p_k, q_k)\|_{H^1}^2 \, ds \leq C$

Schritt 4: Grenzübergang

Unter Verwendung von Kompaktitätsargumenten (Aubin-Lions-Lemma) und Monotoniemethoden wird beim Grenzübergang $k \to \infty$ eine konvergente Teilfolge extrahiert, die eine starke Lösung ergibt.

Technische Innovationen

1. Behandlung nicht-degenerierter Mobilitäten: Im Gegensatz zum Fall konstanter Mobilitäten führen nicht-degenerierte Mobilitäten zusätzliche nichtlineare Terme ein. Die Schlüsseltechnik besteht in der Verwendung elliptischer Regularitätstheorie (Gleichungen 2.13-2.15) zur Kontrolle dieser Terme.
2. Höherordnungs-Energieabschätzungen: Durch Differentiation des Energiefunktionals (Gleichung 4.46), kombiniert mit der Bulk-Surface-Stokes-Regularitätstheorie (Theorem B.2), werden $H^2$-Abschätzungen für Geschwindigkeitsfelder etabliert: $\frac{d}{dt}\text{(Energie)} + \text{Dissipation} \leq \varpi\|(v_k, w_k)\|_{H^2}^2 + \text{niederordnige Terme}$ Durch Wahl geeigneter $\varpi$ wird die $H^2$-Norm absorbiert.
3. Erhaltung der Separationseigenschaft: Unter Verwendung der Singularität logarithmischer Potentiale und einer Variante des Maximumprinzips wird nachgewiesen, dass Phasenfeldvariablen strikt in $(-1, 1)$ bleiben.
4. Eindeutigkeitsbeweis: Für die Differenz $(v, w, \phi, \psi) = (v_1 - v_2, w_1 - w_2, \phi_1 - \phi_2, \psi_1 - \psi_2)$ zweier Lösungen wird eine Differentialungleichung etabliert (Gleichungen 5.2-5.3): $\frac{d}{dt}\left(\|(v, w)\|_{L^2}^2 + \|(\phi, \psi)\|_{K,\alpha}^2\right) + \text{Dissipation} \leq P(t) \times \text{(Norm der Differenz)}$ wobei $P(t) \in L^1(0,T)$. Anwendung des Gronwall-Lemmas ergibt Eindeutigkeit.

Experimentelle Einrichtung

Hinweis: Diese Arbeit ist eine rein theoretische mathematische Arbeit ohne numerische Experimente. Die Hauptergebnisse sind mathematische Beweise für Existenz und Eindeutigkeit.

Theoretisches Verifikationsgerüst

Der Artikel verifiziert die Theorie durch folgende mathematische Werkzeuge:

1. Funktionsraumeinstellung:
   • Geschwindigkeitsfeldraum: $H^1_{0,\text{div}} = \{(v,w) \in H^1_0: \text{div} v = 0, \text{div}_\Gamma w = 0\}$
   • Phasenfeldraum: $V^1_{L,\beta}$ (abhängig von Kopplungsparametern $L, \beta$)
   • Chemisches-Potential-Raum: $H^1_{L,\beta}$
2. Regularitätsanforderungen:
   • Anfangsdaten: $(v_0, w_0) \in H^1_{\text{div}}$, $(\phi_0, \psi_0) \in H^1$ erfüllend $\|\phi_0\|_{L^\infty} \leq 1$, $\|\psi_0\|_{L^\infty} \leq 1$
   • Kompatibilitätsbedingungen: Existenz von $(\mu_0, \theta_0) \in H^1_{L,\beta}$ erfüllend spezifische Integralrelationen
3. Parameterbeschränkungen:
   • Raumdimension: $d = 2$
   • Kopplungsparameter: $K \in (0, \infty)$, $L \in [0, \infty]$ (Eindeutigkeit erfordert $L \in (0, \infty]$)
   • Mobilitäten: $0 < m_* \leq m_\Omega, m_\Gamma \leq m^* < \infty$

Schlüsselannahmen

(A1-A7) umfassen:

• Gebietsregularität: $\Omega \subset \mathbb{R}^2$ mit $C^3$-Rand
• Potentialfunktionen: $F(s) = F_0(s) - \frac{c_F}{2}s^2$, wobei $F_0$ stark konvex und singulär bei $s \to \pm 1$
• Wachstumsbedingungen: $F_0''(s) \geq \Theta > 0$ oder spezifische Wachstumsratenbeschränkungen (A7.1)-(A7.2)

Experimentelle Ergebnisse

Haupttheoretische Resultate

Theorem 3.1 (Existenz globaler starker Lösungen): Unter den Annahmen (A1)-(A7), für $K \in (0, \infty)$, $L \in [0, \infty]$, existiert eine globale starke Lösung $(v, w, \phi, \psi, \mu, \theta)$ erfüllend:

1. Regularität:
   • $(v, w) \in BC([0,\infty); H^1_{0,\text{div}}) \cap L^2_{uloc}([0,\infty); H^2) \cap H^1_{uloc}(0,\infty; L^2_{\text{div}})$
   • $(\phi, \psi) \in L^\infty(0,\infty; H^3)$
   • $(F'(\phi), G'(\psi)), (F''(\phi), G''(\psi)) \in L^\infty(0,\infty; L^p)$ für alle $2 \leq p < \infty$
2. Separationseigenschaft: $|\phi| < 1$ f.ü. in $Q$, $|\psi| < 1$ f.ü. auf $\Sigma$
3. Energieungleichung: $E_{tot}(v(t), w(t), \phi(t), \psi(t)) + \int_0^t \text{(Dissipation)} \, ds = E_{tot}(v_0, w_0, \phi_0, \psi_0)$

Theorem 3.4 (Eindeutigkeit): Unter den Annahmen von Theorem 3.1, falls $L \in (0, \infty]$, ist die starke Lösung eindeutig.

Schlüsselabschätzungen

Die Arbeit etabliert folgende quantitative Abschätzungen (Gleichungen A.10-A.12):

$\|(\phi, \psi)\|_{L^\infty(0,\infty; H^1)}^2 + \|(\partial_t\phi, \partial_t\psi)\|_{L^2(0,\infty; (H^1_{L,\beta})')}^2 \leq C(1 + \|(v,w)\|_{L^2(0,\infty; L^2)}^2)$

$\|(\phi, \psi)\|_{L^2_{uloc}([0,\infty); W^{2,p})}^2 + \|(F'(\phi), G'(\psi))\|_{L^2_{uloc}([0,\infty); L^p)}^2 \leq C_p(1 + \|(v,w)\|_{L^2(0,\infty; L^2)}^2)$

$\|(\phi, \psi)\|_{L^4_{uloc}([0,\infty); H^2)}^4 \leq C(1 + \|(v,w)\|_{L^2(0,\infty; L^2)}^2)^2$

Überwindung technischer Schwierigkeiten

1. Einschränkung für $L = 0$-Fall: Wenn $L = 0$, erlaubt die Randbedingung (1.1h)$_2$ keine Umformulierung der Normalableitung $\partial_n\mu$, was zum Scheitern des Eindeutigkeitsbeweises führt. Der Autor weist darauf hin, dass dies zusätzliche Strukturannahmen erfordert (Bemerkung 3.5).
2. Herausforderungen im dreidimensionalen Fall: Die Arbeit behandelt nur den zweidimensionalen Fall, da im dreidimensionalen Fall:
   • Die Separationseigenschaft verloren geht
   • Feinere Approximationstechniken zur Kontrolle nichtlinearer Terme erforderlich sind
   • Nur lokale starke Lösungen statt globaler starker Lösungen zu erwarten sind

Verwandte Arbeiten

Hauptforschungsrichtungen

1. Modelle ohne dynamische Randbedingungen:
   • Abels-Garcke-Grün (AGG)-Modell: 2, 4, 16, 17, 21, 22
   • Bereits etablierte Theorie starker Lösungen
2. Modelle mit dynamischen Randbedingungen:
   • Hauptsächlich Fokus auf Existenz schwacher Lösungen: 9, 14, 15, 19, 30
   • Quasi-starke Lösungen: Allen-Cahn-Navier-Stokes-Voigt-Systeme 13
3. Regularität von Cahn-Hilliard-Gleichungen:
   • Konstante-Mobilitäts-Fall: 20, 28, 29
   • Nicht-degenerierte Mobilitäten: 11, 39 (Vorarbeiten des Autors)

Innovativität dieser Arbeit

Im Vergleich zu bestehenden Arbeiten:

• Erstmals Behandlung starker Lösungen für Navier-Stokes-Cahn-Hilliard-Systeme mit dynamischen Randbedingungen
• Erstmals Abdeckung des Falls nicht-degenerierter Mobilitäten
• Etablierung neuer Bulk-Surface-Stokes-Gleichungs-Regularitätstheorie (Anhang B)
• Entwicklung höherordniger Abschätzungstechniken für konvektive Cahn-Hilliard-Gleichungen (Anhang A)

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

1. Im zweidimensionalen Fall wird die Existenz und Eindeutigkeit globaler starker Lösungen für das Bulk-Surface-Navier-Stokes-Cahn-Hilliard-System mit nicht-degenerierten Mobilitäten nachgewiesen (wenn $L \in (0, \infty]$)
2. Lösungen erfüllen strikte Separationseigenschaften und Energiedissipationsgesetze
3. Zwei wichtige Hilfstheorien werden etabliert, die unabhängig auf andere Probleme anwendbar sind

Limitierungen

1. Dimensionsbeschränkung: Ergebnisse gelten nur für den zweidimensionalen Fall. Der dreidimensionale Fall erfordert:
   • Feinere nichtlineare Abschätzungen
   • Möglicherweise nur lokale starke Lösungen
   • Für konstante Mobilitäten sind Eindeutigkeit schwacher Lösungen und Existenz starker Lösungen bekannt
2. Parameterbeschränkungen:
   • $K = 0$ ist nicht möglich (nicht-degenerierter Mobilitäts-Fall)
   • Eindeutigkeit für $L = 0$ ist ungelöst, außer unter zusätzlichen Strukturannahmen
3. Kompatibilitätsbedingungen: Wenn $L = 0$, müssen Dichtenfunktionen die spezielle Relation $\beta(\sigma_2 - \sigma_1) = \rho_2 - \rho_1$ erfüllen (Gleichung 3.2)

Zukünftige Richtungen

1. Dreidimensionale Verallgemeinerung: Entwicklung von Techniken für den dreidimensionalen Fall, möglicherweise erfordernde:
   • Lokale Existenztheorie
   • Globale Existenz für kleine Anfangswerte oder unter speziellen Strukturannahmen
   • Konstante-Mobilitäts-Annahme
2. $L = 0$-Fall: Suche nach neuen Methoden zur Behandlung dieses Falls oder Bestimmung seiner Unmöglichkeit
3. Langzeitverhalten: Untersuchung asymptotischen Verhaltens von Lösungen, Existenz von Attraktoren
4. Numerische Verifikation: Entwicklung effizienter numerischer Schemata zur Verifikation theoretischer Vorhersagen

Tiefgehende Bewertung

Stärken

1. Theoretische Strenge:
   • Vollständige und detaillierte Beweise mit standardisierter technischer Behandlung
   • Klare Erklärung verschiedener technischer Schwierigkeiten
   • Präzise Formulierung von Lemmata und Theoremen
2. Innovativität:
   • Erstmalige Etablierung der Theorie starker Lösungen für diese Systemklasse
   • Behandlung nicht-degenerierter Mobilitäten stellt eine Herausforderung dar
   • Hilfstheorien in den Anhängen haben unabhängigen Wert
3. Klare Struktur:
   • Hauptergebnisse sind in Abschnitt 3 klar dargelegt
   • Beweisstrategien sind hierarchisch strukturiert
   • Technische Details in Anhängen beeinträchtigen nicht die Lesbarkeit des Haupttextes
4. Physikalische Plausibilität:
   • Modell ist thermodynamisch konsistent
   • Nicht-degenerierte Mobilitäten entsprechen der Realität besser
   • Energiedissipation und Massenerhaltung sind garantiert

Schwächen

1. Dimensionsbeschränkung: Zweidimensionale Ergebnisse haben begrenzte praktische Anwendbarkeit; der dreidimensionale Fall ist wichtiger, aber ungelöst
2. Parameterbeschränkungen: Die Einschränkungen $K \in (0, \infty)$ und $L \in (0, \infty]$ schließen bestimmte physikalische Szenarien aus
3. Fehlende numerische Verifikation: Als rein theoretische Arbeit fehlen numerische Beispiele zur Unterstützung der theoretischen Ergebnisse
4. Technische Komplexität: Der Beweis ist hochgradig technisiert und erfordert umfangreiches Vorwissen, was die Leserschaft einschränken kann
5. Konstantenabhängigkeit: Viele Abschätzungen enthalten Konstanten $C$, die von mehreren Parametern abhängen, ohne explizite Ausdrücke zu geben

Auswirkungen

1. Theoretische Beiträge:
   • Schließung wichtiger Lücken in der Theorie starker Lösungen für Navier-Stokes-Cahn-Hilliard-Systeme mit dynamischen Randbedingungen
   • Bereitstellung von Techniken und theoretischen Rahmen für nachfolgende Forschung
2. Methodologischer Wert:
   • Die Kombination von Semi-Galerkin-Methode und Fixpunktargument ist auf andere gekoppelte Systeme anwendbar
   • Höherordnungs-Abschätzungstechniken haben allgemeine Gültigkeit
3. Anwendungspotential:
   • Theoretische Unterstützung für numerische Methoden
   • Anwendbarkeit auf Modellierung in Materialwissenschaft und Biophysik von Membranen

Anwendungsszenarien

1. Materialwissenschaft: Phasenseparation und Grenzflächendynamik in zweiphasigen Flüssigkeitsmaterialien
2. Biophysik: Flüssig-Mosaik-Modell von Zellmembranen, Dynamik von Lipiddoppelschichten
3. Kontaktliniendynamik: Bewegung und Benetzungsphänomene von Flüssigkeitstropfen auf festen Oberflächen
4. Theoretische Forschung:
   • Regularitätstheorie nichtlinearer partieller Differentialgleichungen
   • Wohlgestelltheitsanalyse gekoppelter Systeme
   • Mathematische Theorie dynamischer Randbedingungen

Ausgewählte Referenzen

2 H. Abels, H. Garcke, and A. Giorgini. Global regularity and asymptotic stabilization for the incompressible Navier-Stokes-Cahn-Hilliard model with unmatched densities. Math. Ann., 2024.

11 M. Conti, P. Galimberti, S. Gatti, and A. Giorgini. New results for the Cahn-Hilliard equation with non-degenerate mobility. Calc. Var. PDE, 2025.

19 A. Giorgini and P. Knopf. Two-phase flows with bulk-surface interaction: thermodynamically consistent Navier-Stokes-Cahn-Hilliard models. J. Math. Fluid Mech., 2023.

30 P. Knopf and J. Stange. A thermodynamically consistent model for bulk-surface viscous fluid mixtures. Preprint, 2025.

39 J. Stange. Well-posedness and long-time behavior of a bulk-surface Cahn-Hilliard model with non-degenerate mobility. Preprint, 2025.

Gesamtbewertung: Dies ist eine hochwertige theoretische mathematische Arbeit, die bedeutende Beiträge zur Analyse partieller Differentialgleichungen leistet. Trotz Dimensions- und Parameterbeschränkungen verdienen die technischen Innovationen und theoretische Tiefe Anerkennung. Für Mathematiker und angewandte Wissenschaftler in verwandten Bereichen bietet diese Arbeit wertvolle theoretische Werkzeuge und Forschungsparadigmen.