We examine a thermodynamically consistent diffuse interface model for bulk-surface viscous fluid mixtures. This model consists of a Navier--Stokes--Cahn--Hilliard model in the bulk coupled to a surface Navier--Stokes--Cahn--Hilliard system on the boundary. In this paper, we address the global well-posedness of strong solutions in the two-dimensional setting, also covering the physically meaningful case of non-degenerate mobility functions. Lastly, we prove the uniqueness of the corresponding strong solutions and their continuous dependence on the initial data. Our approach hinges upon new well-posedness and regularity theory for a convective bulk-surface Cahn--Hilliard equation with non-degenerate mobilities, as well as a bulk-surface Stokes equation with non-constant coefficients.
- ID de l'article: 2511.06847
- Titre: Global well-posedness of strong solutions to a bulk-surface Navier-Stokes-Cahn-Hilliard model with non-degenerate mobilities in two dimensions
- Auteur: Jonas Stange (Universität Regensburg, Allemagne)
- Classification: math.AP (Analyse des EDP), math-ph (Physique mathématique), math.MP (Physique mathématique)
- Date de soumission: 10 novembre 2025
- Lien de l'article: https://arxiv.org/abs/2511.06847v1
Cet article étudie un modèle d'interface diffuse thermodynamiquement cohérent pour décrire les mélanges de fluides visqueux en volume-surface. Le modèle est constitué des équations de Navier-Stokes-Cahn-Hilliard dans le volume couplées à un système de Navier-Stokes-Cahn-Hilliard de surface sur la frontière. L'auteur démontre la bonne posture globale des solutions fortes en dimension deux, couvrant le cas physiquement pertinent des fonctions de mobilité non dégénérées. De plus, l'unicité des solutions fortes et leur dépendance continue par rapport aux données initiales sont établies. L'approche repose sur une nouvelle théorie de bonne posture et de régularité pour les équations de Cahn-Hilliard convectives en volume-surface avec mobilités non dégénérées, ainsi que sur les équations de Stokes en volume-surface avec coefficients non constants.
- Modélisation des phénomènes physiques: L'étude porte sur les modèles d'interface diffuse pour les mélanges de fluides biphasiques où les fluides ont des densités différentes et où les effets de viscosité dynamique sont présents à la fois dans le domaine volumique et sur la frontière. Ces modèles trouvent des applications importantes en science des matériaux et en dynamique des membranes biologiques.
- Défis mathématiques: Le modèle couple:
- Les équations de Navier-Stokes-Cahn-Hilliard dans le domaine volumique (décrivant le mouvement du fluide et la séparation de phase)
- Le système de Navier-Stokes-Cahn-Hilliard de surface sur la frontière (décrivant la dynamique de surface)
- Les conditions de frontière dynamiques (permettant l'échange de masse entre le volume et la surface)
- Limitations de la théorie existante:
- Pour les modèles d'interface diffuse avec conditions de frontière dynamiques, la littérature existante se concentre principalement sur l'existence de solutions faibles
- L'étude de l'existence de solutions fortes est moins développée, particulièrement pour le cas des fonctions de mobilité non dégénérées
- La théorie des solutions fortes globales en dimension trois reste incomplète
- Pertinence physique: Les fonctions de mobilité non dégénérées (c'est-à-dire avec des bornes positives inférieures et supérieures sur tout le domaine) correspondent mieux à la réalité physique
- Rigueur mathématique: Les solutions fortes fournissent des estimations de régularité plus fines, essentielles pour comprendre le comportement à long terme des solutions
- Complétude théorique: Combler les lacunes dans la théorie des solutions fortes pour les systèmes de Navier-Stokes-Cahn-Hilliard avec conditions de frontière dynamiques
Les principales contributions de cet article sont:
- Existence de solutions fortes globales: En dimension deux, l'existence de solutions fortes globales pour le système de Navier-Stokes-Cahn-Hilliard en volume-surface avec mobilités non dégénérées est démontrée (Théorème 3.1)
- Unicité et dépendance continue: L'unicité des solutions fortes et leur dépendance continue par rapport aux données initiales sont établies (Théorème 3.4), sous la condition de restriction L∈(0,∞]
- Nouvelle théorie de régularité: Deux théories auxiliaires clés sont établies:
- Théorie de bonne posture et de régularité pour les équations de Cahn-Hilliard convectives en volume-surface avec mobilités non dégénérées (Appendice A)
- Théorie de régularité pour les équations de Stokes en volume-surface avec coefficients non constants (Appendice B)
- Inégalités énergétiques: Les solutions satisfont la loi de dissipation énergétique:
dtdEtot(v,w,ϕ,ψ)+termes de dissipation=0
- Propriété de séparation: Les fonctions de phase restent strictement dans l'intervalle physique (−1,1): ∣ϕ∣<1 p.p. dans Q, ∣ψ∣<1 p.p. sur Σ
Considérons un domaine borné bidimensionnel Ω⊂R2 et sa frontière Γ=∂Ω. Posons Q=Ω×(0,∞), Σ=Γ×(0,∞).
Système d'équations gouvernantes (équation 1.1) comprend:
- Équation de quantité de mouvement volumique:
∂t(ρ(ϕ)v)+div(v⊗(ρ(ϕ)v+J))=divT,divv=0 dans Q
- Équation de quantité de mouvement de surface:
∂t(σ(ψ)w)+divΓ(w⊗(σ(ψ)w+K))=divΓTΓ+Z,divΓw=0 sur Σ
- Équation de Cahn-Hilliard volumique:
∂tϕ+div(ϕv)=div(mΩ(ϕ)∇μ),μ=−Δϕ+F′(ϕ) dans Q
- Équation de Cahn-Hilliard de surface:
∂tψ+divΓ(ψw)=divΓ(mΓ(ψ)∇Γθ)−βmΩ(ϕ)∂nμ sur Σθ=−ΔΓψ+G′(ψ)+α∂nϕ sur Σ
- Conditions de frontière:
w=v∣Γ,v⋅n=0 sur ΣK∂nϕ=αψ−ϕ,LmΩ(ϕ)∂nμ=βθ−μ sur Σ
Où:
- v,w sont les champs de vitesse volumique et de surface
- ϕ,ψ sont les variables de phase volumique et de surface
- μ,θ sont les potentiels chimiques
- ρ(ϕ),σ(ψ) sont les fonctions de densité dépendant de la phase
- mΩ,mΓ sont les fonctions de mobilité (satisfaisant 0<m∗≤mΩ,mΓ≤m∗)
La stratégie de preuve utilise une méthode semi-Galerkin:
- Introduire le système de fonctions propres {(v~j,w~j)}j∈N de l'opérateur de Stokes en volume-surface A
- Définir le sous-espace de dimension finie Vk=span{(v~1,w~1),…,(v~k,w~k)}
- Construire la solution approchée (vk,wk,ϕk,ψk,μk,θk) satisfaisant la forme faible Galerkin des équations
Pour un champ de vitesse fixé (v∗,w∗), résoudre:
- L'équation de Cahn-Hilliard convective en volume-surface (4.8), obtenant (ϕk,ψk,μk,θk)
- L'approximation Galerkin de l'équation de Navier-Stokes (4.12), obtenant (vk,wk)
L'existence d'un point fixe de l'application Λ:(v∗,w∗)↦(vk,wk) est démontrée via le théorème du point fixe de Schauder.
Des estimations uniformes clés (indépendantes de k) sont établies:
- Estimations énergétiques (équations 4.35-4.37):
∥(vk,wk)∥L∞(0,∞;L2)+∥(vk,wk)∥L2(0,∞;H1)≤C∥(ϕk,ψk)∥L∞(0,∞;H1)≤C
- Estimations d'ordre supérieur (équations 4.43-4.45):
∥(ϕk,ψk)∥L∞(0,∞;W2,p)+∥(F′(ϕk),G′(ψk))∥L∞(0,∞;Lp)≤K1,p∥(μk,θk)∥L∞(0,∞;H1)+∥(μk,θk)∥Luloc2(0,∞;H3)≤K2
- Estimations d'ordre supérieur du champ de vitesse (équation 4.61):
supt≥0∫tt+1∥(vk,wk)∥H22+∥(∂tvk,∂twk)∥L22+∥(pk,qk)∥H12ds≤C
Utilisant des arguments de compacité (lemme d'Aubin-Lions) et des techniques de monotonie, on extrait une sous-suite convergente quand k→∞, obtenant une solution forte.
- Traitement des mobilités non dégénérées: Contrairement au cas des mobilités constantes, les mobilités non dégénérées introduisent des termes non linéaires supplémentaires. La technique clé est d'utiliser la théorie de régularité elliptique (équations 2.13-2.15) pour contrôler ces termes.
- Estimations énergétiques d'ordre supérieur: En différentiant la fonctionnelle énergétique (équation 4.46) et en combinant avec la théorie de régularité de Stokes en volume-surface (Théorème B.2), on établit l'estimation H2 du champ de vitesse:
dtd(eˊnergie)+dissipation≤ϖ∥(vk,wk)∥H22+termes d’ordre infeˊrieur
En choisissant ϖ approprié, on absorbe la norme H2.
- Préservation de la propriété de séparation: Utilisant la singularité du potentiel logarithmique et une variante du principe du maximum, on démontre que les fonctions de phase restent strictement dans (−1,1).
- Preuve d'unicité: Pour la différence de deux solutions (v,w,ϕ,ψ)=(v1−v2,w1−w2,ϕ1−ϕ2,ψ1−ψ2), on établit une inégalité différentielle (équations 5.2-5.3):
dtd(∥(v,w)∥L22+∥(ϕ,ψ)∥K,α2)+dissipation≤P(t)×(norme de la diffeˊrence)
où P(t)∈L1(0,T), et on applique le lemme de Gronwall pour obtenir l'unicité.
Remarque: Cet article est un travail mathématique purement théorique et ne contient pas d'expériences numériques. Les résultats principaux sont des preuves mathématiques d'existence et d'unicité.
L'article vérifie la théorie via les outils mathématiques suivants:
- Configuration des espaces fonctionnels:
- Espace des champs de vitesse: H0,div1={(v,w)∈H01:divv=0,divΓw=0}
- Espace des phases: VL,β1 (dépendant des paramètres de couplage L,β)
- Espace des potentiels chimiques: HL,β1
- Exigences de régularité:
- Données initiales: (v0,w0)∈Hdiv1, (ϕ0,ψ0)∈H1 satisfaisant ∥ϕ0∥L∞≤1, ∥ψ0∥L∞≤1
- Conditions de compatibilité: existence de (μ0,θ0)∈HL,β1 satisfaisant certaines relations intégrales
- Restrictions de paramètres:
- Dimension spatiale: d=2
- Paramètres de couplage: K∈(0,∞), L∈[0,∞] (l'unicité nécessite L∈(0,∞])
- Mobilités: 0<m∗≤mΩ,mΓ≤m∗<∞
(A1-A7) incluent:
- Régularité du domaine: Ω⊂R2 avec frontière C3
- Fonctions potentielles: F(s)=F0(s)−2cFs2, où F0 est fortement convexe et singulière quand s→±1
- Conditions de croissance: F0′′(s)≥Θ>0 ou taux de croissance spécifique (conditions (A7.1)-(A7.2))
Théorème 3.1 (Existence de solutions fortes globales): Sous les hypothèses (A1)-(A7), pour K∈(0,∞), L∈[0,∞], il existe une solution forte globale (v,w,ϕ,ψ,μ,θ) satisfaisant:
- Régularité:
- (v,w)∈BC([0,∞);H0,div1)∩Luloc2([0,∞);H2)∩Huloc1(0,∞;Ldiv2)
- (ϕ,ψ)∈L∞(0,∞;H3)
- (F′(ϕ),G′(ψ)),(F′′(ϕ),G′′(ψ))∈L∞(0,∞;Lp) pour tout 2≤p<∞
- Propriété de séparation: ∣ϕ∣<1 p.p. dans Q, ∣ψ∣<1 p.p. sur Σ
- Inégalité énergétique:
Etot(v(t),w(t),ϕ(t),ψ(t))+∫0t(dissipation)ds=Etot(v0,w0,ϕ0,ψ0)
Théorème 3.4 (Unicité): Sous les hypothèses du Théorème 3.1, si L∈(0,∞], alors la solution forte est unique.
L'article établit les estimations quantitatives suivantes (équations A.10-A.12):
∥(ϕ,ψ)∥L∞(0,∞;H1)2+∥(∂tϕ,∂tψ)∥L2(0,∞;(HL,β1)′)2≤C(1+∥(v,w)∥L2(0,∞;L2)2)
∥(ϕ,ψ)∥Luloc2([0,∞);W2,p)2+∥(F′(ϕ),G′(ψ))∥Luloc2([0,∞);Lp)2≤Cp(1+∥(v,w)∥L2(0,∞;L2)2)
∥(ϕ,ψ)∥Luloc4([0,∞);H2)4≤C(1+∥(v,w)∥L2(0,∞;L2)2)2
- Limitation du cas L=0: Quand L=0, la condition de frontière (1.1h)2 ne permet pas de reformuler la dérivée normale ∂nμ, ce qui rend la preuve d'unicité inefficace. L'auteur note que cela nécessite des hypothèses structurelles supplémentaires (Remarque 3.5).
- Défis en dimension trois: L'article traite uniquement le cas bidimensionnel car en dimension trois:
- La propriété de séparation n'est pas préservée
- Des techniques d'approximation plus fines sont nécessaires pour contrôler les termes non linéaires
- Seules les solutions fortes locales plutôt que globales peuvent être attendues
- Modèles sans conditions de frontière dynamiques:
- Modèle AGG (Abels-Garcke-Grün): 2, 4, 16, 17, 21, 22
- Théorie mature des solutions fortes existante
- Modèles avec conditions de frontière dynamiques:
- Accent principal sur l'existence de solutions faibles: 9, 14, 15, 19, 30
- Solutions quasi-fortes: système Allen-Cahn-Navier-Stokes-Voigt 13
- Régularité de l'équation de Cahn-Hilliard:
- Cas des mobilités constantes: 20, 28, 29
- Mobilités non dégénérées: 11, 39 (travaux antérieurs de l'auteur)
Par rapport aux travaux existants, cet article:
- Pour la première fois traite les solutions fortes du système de Navier-Stokes-Cahn-Hilliard avec conditions de frontière dynamiques
- Pour la première fois couvre le cas des mobilités non dégénérées
- Établit une nouvelle théorie de régularité pour les équations de Stokes en volume-surface (Appendice B)
- Développe des techniques d'estimations d'ordre supérieur pour l'équation de Cahn-Hilliard convective (Appendice A)
- En dimension deux, l'existence et l'unicité (quand L∈(0,∞]) des solutions fortes globales pour le système de Navier-Stokes-Cahn-Hilliard en volume-surface avec mobilités non dégénérées sont démontrées
- Les solutions satisfont une propriété de séparation stricte et une loi de dissipation énergétique
- Deux cadres théoriques auxiliaires importants sont établis, applicables indépendamment à d'autres problèmes
- Restriction dimensionnelle: Les résultats s'appliquent uniquement en dimension deux. Le cas tridimensionnel nécessite:
- Des estimations non linéaires plus fines
- Possiblement seulement des solutions fortes locales
- Pour le cas des mobilités constantes, l'unicité des solutions faibles et l'existence de solutions fortes sont connues
- Restrictions de paramètres:
- K=0 n'est pas viable (cas des mobilités non dégénérées)
- L'unicité pour L=0 n'est pas résolue, sauf avec des hypothèses structurelles supplémentaires
- Conditions de compatibilité: Quand L=0, les fonctions de densité doivent satisfaire la relation spéciale β(σ2−σ1)=ρ2−ρ1 (équation 3.2)
- Généralisation tridimensionnelle: Développer des techniques applicables en dimension trois, possiblement nécessitant:
- Théorie d'existence locale
- Globalité pour petites données initiales ou structures spéciales
- Hypothèse de mobilités constantes
- Cas L=0: Chercher de nouvelles méthodes pour traiter ce cas ou établir son infaisabilité
- Comportement à long terme: Étudier le comportement asymptotique des solutions, existence d'attracteurs
- Vérification numérique: Développer des schémas numériques efficaces pour vérifier les prédictions théoriques
- Rigueur théorique:
- Preuve complète et détaillée, traitement technique conforme aux normes
- Explication claire des diverses difficultés techniques
- Énoncés des lemmes et théorèmes précis
- Caractère innovant:
- Première établissement de la théorie des solutions fortes pour cette classe de systèmes
- Le traitement des mobilités non dégénérées présente des défis significatifs
- Les théories auxiliaires dans les appendices ont une valeur indépendante
- Structure claire:
- Les résultats principaux sont énoncés clairement à la Section 3
- La stratégie de preuve est bien organisée hiérarchiquement
- Les détails techniques en appendice ne perturbent pas la lecture du fil principal
- Pertinence physique:
- Modèle thermodynamiquement cohérent
- Les mobilités non dégénérées correspondent mieux à la réalité
- Conservation de la masse et loi de dissipation énergétique garanties
- Restriction dimensionnelle: Les résultats en dimension deux ont une applicabilité limitée, le cas tridimensionnel étant plus important mais non résolu
- Restrictions de paramètres: Les restrictions K∈(0,∞) et L∈(0,∞] excluent certains cas physiques
- Absence de vérification numérique: En tant que travail purement théorique, aucun exemple numérique n'est fourni pour soutenir les résultats théoriques
- Complexité technique: La preuve est hautement technique, nécessitant des connaissances préalables substantielles, ce qui peut limiter le lectorat
- Dépendance des constantes: De nombreuses estimations contiennent des constantes C dépendant de plusieurs paramètres, sans expressions explicites
- Contributions théoriques:
- Comble une lacune importante dans la théorie des solutions fortes pour les systèmes de Navier-Stokes-Cahn-Hilliard avec conditions de frontière dynamiques
- Fournit des outils techniques et un cadre théorique pour les recherches ultérieures
- Valeur méthodologique:
- La combinaison de la méthode semi-Galerkin et de l'argument du point fixe peut s'appliquer à d'autres systèmes couplés
- Les techniques d'estimations d'ordre supérieur ont un caractère général
- Potentiel applicatif:
- Fournit un soutien théorique aux méthodes numériques
- Applicable à la modélisation en science des matériaux, dynamique des membranes biologiques, etc.
- Science des matériaux: Séparation de phase et dynamique d'interface des matériaux fluides biphasiques
- Biophysique: Modèle de mosaïque fluide des membranes cellulaires, dynamique des bicouches lipidiques
- Dynamique de la ligne de contact: Mouvement et mouillage des gouttelettes sur surfaces solides
- Recherche théorique:
- Théorie de régularité des équations aux dérivées partielles non linéaires
- Analyse de bonne posture des systèmes couplés
- Théorie mathématique des conditions de frontière dynamiques
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.
Évaluation globale: Ceci est un article mathématique théorique de haute qualité qui apporte des contributions importantes au domaine de l'analyse des équations aux dérivées partielles. Malgré les restrictions dimensionnelles et paramétriques, son innovation technique et sa profondeur théorique méritent reconnaissance. Pour les mathématiciens et scientifiques appliqués travaillant dans les domaines connexes, cet article fournit des outils théoriques précieux et un paradigme de recherche.