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 Articolo: 2511.06847
- Titolo: Global well-posedness of strong solutions to a bulk-surface Navier-Stokes-Cahn-Hilliard model with non-degenerate mobilities in two dimensions
- Autore: Jonas Stange (Universität Regensburg, Germania)
- Classificazione: math.AP (Analisi di EDP), math-ph (Fisica Matematica), math.MP (Fisica Matematica)
- Data di Sottomissione: 10 novembre 2025
- Link Articolo: https://arxiv.org/abs/2511.06847v1
Il presente articolo studia un modello di interfaccia diffusa termodinamicamente coerente per descrivere miscele di fluidi viscosi bulk-surface. Il modello è costituito dalle equazioni di Navier-Stokes-Cahn-Hilliard nel dominio bulk accoppiate con il sistema di Navier-Stokes-Cahn-Hilliard di superficie al contorno. L'autore dimostra la buona posizione globale di soluzioni forti nel caso bidimensionale, coprendo il caso fisicamente significativo di funzioni di mobilità non degeneri. Inoltre, vengono provati l'unicità della soluzione forte e la dipendenza continua dai dati iniziali. L'approccio si basa su una nuova teoria di buona posizione e regolarità per l'equazione di Cahn-Hilliard convettiva bulk-surface con mobilità non degeneri, nonché sulla teoria di regolarità per l'equazione di Stokes bulk-surface con coefficienti non costanti.
- Modellazione di Fenomeni Fisici: La ricerca riguarda modelli di interfaccia diffusa per miscele di fluidi bifasici dove i fluidi hanno densità diverse e gli effetti dinamici viscosi sono presenti sia nel dominio bulk che al contorno. Tali modelli trovano importanti applicazioni nella scienza dei materiali, nella dinamica delle membrane biologiche e in altri campi.
- Sfide Matematiche: Il modello accoppia:
- Equazioni di Navier-Stokes-Cahn-Hilliard nel dominio bulk (descrivono il moto del fluido e la separazione di fase)
- Sistema di Navier-Stokes-Cahn-Hilliard di superficie al contorno (descrive la dinamica di superficie)
- Condizioni al contorno dinamiche (permettono lo scambio di massa tra bulk e superficie)
- Limitazioni della Teoria Esistente:
- Per modelli di interfaccia diffusa con condizioni al contorno dinamiche, la letteratura esistente si concentra principalmente sull'esistenza di soluzioni deboli
- Lo studio dell'esistenza di soluzioni forti è meno frequente, specialmente per il caso di funzioni di mobilità non degeneri
- La teoria globale di soluzioni forti nel caso tridimensionale rimane incompleta
- Ragionevolezza Fisica: Le funzioni di mobilità non degeneri (cioè mobilità con limiti inferiori e superiori positivi in tutto il dominio) sono più coerenti con la realtà fisica
- Rigore Matematico: Le soluzioni forti forniscono stime di regolarità più raffinate, cruciali per comprendere il comportamento a lungo termine delle soluzioni
- Completamento della Teoria: Colma il vuoto nella teoria delle soluzioni forti per sistemi di Navier-Stokes-Cahn-Hilliard con condizioni al contorno dinamiche
I principali contributi dell'articolo includono:
- Esistenza di Soluzioni Forti Globali: Nel caso bidimensionale, viene provata l'esistenza di soluzioni forti globali per il sistema bulk-surface di Navier-Stokes-Cahn-Hilliard con mobilità non degeneri (Teorema 3.1)
- Unicità e Dipendenza Continua: Viene provata l'unicità della soluzione forte e la dipendenza continua dai dati iniziali (Teorema 3.4), con la condizione di restrizione L∈(0,∞]
- Nuova Teoria di Regolarità: Vengono stabilite due teorie ausiliarie fondamentali:
- Teoria di buona posizione e regolarità per l'equazione di Cahn-Hilliard convettiva bulk-surface con mobilità non degeneri (Appendice A)
- Teoria di regolarità per l'equazione di Stokes bulk-surface con coefficienti non costanti (Appendice B)
- Disuguaglianza Energetica: Viene provato che la soluzione soddisfa la legge di dissipazione energetica:
dtdEtot(v,w,ϕ,ψ)+termini di dissipazione=0
- Proprietà di Separazione: Viene provato che le funzioni di fase rimangono rigorosamente nell'intervallo fisico (−1,1): ∣ϕ∣<1 q.o. in Q, ∣ψ∣<1 q.o. su Σ
Si consideri un dominio limitato bidimensionale Ω⊂R2 e il suo contorno Γ=∂Ω. Si ponga Q=Ω×(0,∞), Σ=Γ×(0,∞).
Sistema di Equazioni di Controllo (equazione 1.1) include:
- Equazione di Quantità di Moto nel Dominio Bulk:
∂t(ρ(ϕ)v)+div(v⊗(ρ(ϕ)v+J))=divT,divv=0 in Q
- Equazione di Quantità di Moto di Superficie:
∂t(σ(ψ)w)+divΓ(w⊗(σ(ψ)w+K))=divΓTΓ+Z,divΓw=0 on Σ
- Equazione di Cahn-Hilliard nel Dominio Bulk:
∂tϕ+div(ϕv)=div(mΩ(ϕ)∇μ),μ=−Δϕ+F′(ϕ) in Q
- Equazione di Cahn-Hilliard di Superficie:
∂tψ+divΓ(ψw)=divΓ(mΓ(ψ)∇Γθ)−βmΩ(ϕ)∂nμ on Σθ=−ΔΓψ+G′(ψ)+α∂nϕ on Σ
- Condizioni al Contorno:
w=v∣Γ,v⋅n=0 on ΣK∂nϕ=αψ−ϕ,LmΩ(ϕ)∂nμ=βθ−μ on Σ
Dove:
- v,w sono i campi di velocità nel bulk e sulla superficie
- ϕ,ψ sono le variabili di fase nel bulk e sulla superficie
- μ,θ sono i potenziali chimici
- ρ(ϕ),σ(ψ) sono funzioni di densità dipendenti dalla fase
- mΩ,mΓ sono funzioni di mobilità (soddisfano 0<m∗≤mΩ,mΓ≤m∗)
La Strategia di Prova utilizza il metodo di approssimazione semi-Galerkin:
- Si introduce il sistema di autofunzioni {(v~j,w~j)}j∈N dell'operatore di Stokes bulk-surface A
- Si definisce il sottospazio a dimensione finita Vk=span{(v~1,w~1),…,(v~k,w~k)}
- Si costruisce la soluzione approssimata (vk,wk,ϕk,ψk,μk,θk) soddisfacente la forma debole Galerkin delle equazioni
Per un campo di velocità fissato (v∗,w∗), si risolvono:
- L'equazione di Cahn-Hilliard convettiva bulk-surface (4.8), ottenendo (ϕk,ψk,μk,θk)
- L'approssimazione Galerkin dell'equazione di Navier-Stokes (4.12), ottenendo (vk,wk)
Mediante il teorema del punto fisso di Schauder si prova che la mappa Λ:(v∗,w∗)↦(vk,wk) ammette un punto fisso.
Si stabiliscono stime uniformi cruciali (indipendenti da k):
- Stima Energetica (equazioni 4.35-4.37):
∥(vk,wk)∥L∞(0,∞;L2)+∥(vk,wk)∥L2(0,∞;H1)≤C∥(ϕk,ψk)∥L∞(0,∞;H1)≤C
- Stime di Ordine Superiore (equazioni 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
- Stime di Ordine Superiore del Campo di Velocità (equazione 4.61):
supt≥0∫tt+1∥(vk,wk)∥H22+∥(∂tvk,∂twk)∥L22+∥(pk,qk)∥H12ds≤C
Utilizzando argomenti di compattezza (lemma di Aubin-Lions) e tecniche di monotonia, si estrae una sottosuccessione convergente quando k→∞, ottenendo la soluzione forte.
- Trattamento della Mobilità Non Degenere: Diversamente dal caso di mobilità costante, la mobilità non degenere introduce termini non lineari aggiuntivi. La tecnica chiave è utilizzare la teoria di regolarità ellittica (equazioni 2.13-2.15) per controllare questi termini.
- Stime Energetiche di Ordine Superiore: Mediante differenziazione del funzionale energetico (equazione 4.46), combinata con la teoria di regolarità di Stokes bulk-surface (Teorema B.2), si stabilisce la stima H2 del campo di velocità:
dtd(energia)+dissipazione≤ϖ∥(vk,wk)∥H22+termini di ordine inferiore
Scegliendo opportunamente ϖ, si assorbe la norma H2.
- Preservazione della Proprietà di Separazione: Utilizzando la singolarità del potenziale logaritmico e una variante del principio del massimo, si prova che le funzioni di fase rimangono rigorosamente in (−1,1).
- Prova di Unicità: Per la differenza di due soluzioni (v,w,ϕ,ψ)=(v1−v2,w1−w2,ϕ1−ϕ2,ψ1−ψ2), si stabilisce una disuguaglianza differenziale (equazioni 5.2-5.3):
dtd(∥(v,w)∥L22+∥(ϕ,ψ)∥K,α2)+dissipazione≤P(t)×(norma della differenza)
dove P(t)∈L1(0,T), applicando il lemma di Gronwall si ottiene l'unicità.
Nota: Il presente articolo è un lavoro di matematica teorica pura e non contiene esperimenti numerici. I risultati principali sono prove matematiche di esistenza e unicità.
L'articolo verifica la teoria attraverso i seguenti strumenti matematici:
- Configurazione degli Spazi Funzionali:
- Spazio dei campi di velocità: H0,div1={(v,w)∈H01:divv=0,divΓw=0}
- Spazio delle fasi: VL,β1 (dipendente dai parametri di accoppiamento L,β)
- Spazio dei potenziali chimici: HL,β1
- Requisiti di Regolarità:
- Dati iniziali: (v0,w0)∈Hdiv1, (ϕ0,ψ0)∈H1 soddisfacenti ∥ϕ0∥L∞≤1, ∥ψ0∥L∞≤1
- Condizioni di compatibilità: esistono (μ0,θ0)∈HL,β1 soddisfacenti specifiche relazioni integrali
- Restrizioni sui Parametri:
- Dimensione spaziale: d=2
- Parametri di accoppiamento: K∈(0,∞), L∈[0,∞] (l'unicità richiede L∈(0,∞])
- Mobilità: 0<m∗≤mΩ,mΓ≤m∗<∞
(A1-A7) includono:
- Regolarità del dominio: Ω⊂R2 con contorno C3
- Funzioni potenziale: F(s)=F0(s)−2cFs2, dove F0 è fortemente convessa e singolare per s→±1
- Condizioni di crescita: F0′′(s)≥Θ>0 o tassi di crescita specifici ((A7.1)-(A7.2))
Teorema 3.1 (Esistenza di Soluzioni Forti Globali): Sotto le ipotesi (A1)-(A7), per K∈(0,∞), L∈[0,∞], esiste una soluzione forte globale (v,w,ϕ,ψ,μ,θ) soddisfacente:
- Regolarità:
- (v,w)∈BC([0,∞);H0,div1)∩Luloc2([0,∞);H2)∩Huloc1(0,∞;Ldiv2)
- (ϕ,ψ)∈L∞(0,∞;H3)
- (F′(ϕ),G′(ψ)),(F′′(ϕ),G′′(ψ))∈L∞(0,∞;Lp) per ogni 2≤p<∞
- Proprietà di Separazione: ∣ϕ∣<1 q.o. in Q, ∣ψ∣<1 q.o. su Σ
- Disuguaglianza Energetica:
Etot(v(t),w(t),ϕ(t),ψ(t))+∫0t(dissipazione)ds=Etot(v0,w0,ϕ0,ψ0)
Teorema 3.4 (Unicità): Sotto le ipotesi del Teorema 3.1, se L∈(0,∞], allora la soluzione forte è unica.
L'articolo stabilisce le seguenti stime quantitative (equazioni 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
- Limitazione del Caso L=0: Quando L=0, la condizione al contorno (1.1h)2 non consente di riformulare la derivata normale ∂nμ, causando il fallimento della prova di unicità. L'autore sottolinea che ciò richiede ipotesi strutturali aggiuntive (Osservazione 3.5).
- Sfide nel Caso Tridimensionale: L'articolo affronta solo il caso bidimensionale, poiché nel caso tridimensionale:
- Si perde la propagazione della proprietà di separazione
- Sono necessarie tecniche di approssimazione più raffinate per controllare i termini non lineari
- Si può sperare solo in soluzioni forti locali piuttosto che globali
- Modelli Senza Condizioni al Contorno Dinamiche:
- Modello AGG (Abels-Garcke-Grün): 2, 4, 16, 17, 21, 22
- Teoria matura di soluzioni forti già disponibile
- Modelli con Condizioni al Contorno Dinamiche:
- Principalmente focalizzati sull'esistenza di soluzioni deboli: 9, 14, 15, 19, 30
- Quasi-soluzioni forti: sistema Allen-Cahn-Navier-Stokes-Voigt 13
- Regolarità dell'Equazione di Cahn-Hilliard:
- Caso di mobilità costante: 20, 28, 29
- Mobilità non degenere: 11, 39 (lavori precedenti dell'autore)
Rispetto ai lavori esistenti, il presente articolo:
- Per la prima volta affronta le soluzioni forti del sistema di Navier-Stokes-Cahn-Hilliard bulk-surface con condizioni al contorno dinamiche
- Per la prima volta copre il caso di mobilità non degeneri
- Stabilisce una nuova teoria di regolarità per l'equazione di Stokes bulk-surface (Appendice B)
- Sviluppa tecniche di stime di ordine superiore per l'equazione di Cahn-Hilliard convettiva (Appendice A)
- Nel caso bidimensionale, viene provata l'esistenza e l'unicità di soluzioni forti globali per il sistema bulk-surface di Navier-Stokes-Cahn-Hilliard con mobilità non degeneri (quando L∈(0,∞])
- La soluzione soddisfa la proprietà di separazione rigorosa e la legge di dissipazione energetica
- Vengono stabilite due importanti strutture teoriche ausiliarie, applicabili indipendentemente ad altri problemi
- Restrizione sulla Dimensione: I risultati si applicano solo al caso bidimensionale. Il caso tridimensionale richiede:
- Stime non lineari più raffinate
- Possibilmente solo soluzioni forti locali
- Nel caso di mobilità costante, l'unicità della soluzione debole e l'esistenza della soluzione forte sono note
- Restrizioni sui Parametri:
- K=0 non è fattibile (caso di mobilità non degenere)
- L'unicità per L=0 rimane irrisolta, a meno di ipotesi strutturali aggiuntive
- Condizioni di Compatibilità: Quando L=0, la funzione di densità deve soddisfare la relazione speciale β(σ2−σ1)=ρ2−ρ1 (equazione 3.2)
- Generalizzazione Tridimensionale: Sviluppare tecniche applicabili al caso tridimensionale, possibilmente richiedendo:
- Teoria di esistenza locale
- Globalità per piccoli dati iniziali o strutture speciali
- Ipotesi di mobilità costante
- Caso L=0: Ricercare nuovi metodi per affrontare questo caso o determinarne l'infeasibilità
- Comportamento a Lungo Termine: Studiare il comportamento asintotico delle soluzioni, esistenza di attrattori
- Verifica Numerica: Sviluppare schemi numerici efficienti per verificare le previsioni teoriche
- Rigore Teorico:
- Prove complete e dettagliate, trattamento tecnico rigoroso
- Spiegazioni chiare delle varie difficoltà tecniche
- Enunciati di lemmi e teoremi precisi
- Innovatività:
- Primo stabilimento della teoria di soluzioni forti per questa classe di sistemi
- Il trattamento della mobilità non degenere presenta sfide significative
- Le teorie ausiliarie nell'appendice hanno valore indipendente
- Struttura Chiara:
- I risultati principali sono chiaramente enunciati nella Sezione 3
- La strategia di prova è ben articolata
- I dettagli tecnici sono collocati nell'appendice, non ostacolando la lettura della linea principale
- Ragionevolezza Fisica:
- Il modello è termodinamicamente coerente
- La mobilità non degenere è più realistica
- La conservazione della massa e la legge di dissipazione energetica sono garantite
- Restrizione sulla Dimensione: I risultati bidimensionali hanno applicabilità limitata; il caso tridimensionale è più importante ma rimane irrisolto
- Restrizioni sui Parametri: Le limitazioni K∈(0,∞) e L∈(0,∞] escludono alcuni scenari fisici
- Assenza di Verifica Numerica: Come lavoro puramente teorico, non fornisce esempi numerici a supporto dei risultati teorici
- Complessità Tecnica: La prova è altamente tecnica, richiedendo vaste conoscenze preliminari, potenzialmente limitando il pubblico di lettori
- Dipendenza delle Costanti: Molte stime contengono costanti C dipendenti da molteplici parametri, senza espressioni esplicite
- Contributi Teorici:
- Colma un importante vuoto nella teoria delle soluzioni forti per sistemi di Navier-Stokes-Cahn-Hilliard bulk-surface con condizioni al contorno dinamiche
- Fornisce strumenti tecnici e quadri teorici per ricerche successive
- Valore Metodologico:
- La combinazione del metodo semi-Galerkin con argomenti del punto fisso è applicabile ad altri sistemi accoppiati
- Le tecniche di stime di ordine superiore hanno carattere generale
- Potenziale Applicativo:
- Fornisce supporto teorico per metodi numerici
- Applicabile alla modellazione in scienza dei materiali, dinamica di membrane biologiche e altri campi
- Scienza dei Materiali: Separazione di fase e dinamica di interfaccia in materiali fluidi bifasici
- Biofisica: Modello a mosaico fluido di membrane cellulari, dinamica di doppi strati lipidici
- Dinamica della Linea di Contatto: Movimento e bagnabilità di gocce su superfici solide
- Ricerca Teorica:
- Teoria di regolarità per equazioni differenziali alle derivate parziali non lineari
- Analisi di buona posizione per sistemi accoppiati
- Teoria matematica delle condizioni al contorno dinamiche
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.
Valutazione Complessiva: Questo è un articolo di alta qualità nel campo dell'analisi matematica delle equazioni differenziali alle derivate parziali, che apporta contributi significativi. Nonostante le limitazioni sulla dimensione e sui parametri, l'innovazione tecnica e la profondità teorica meritano riconoscimento. Per i matematici e gli scienziati applicati che lavorano in campi correlati, questo articolo fornisce preziosi strumenti teorici e paradigmi di ricerca.