On weak solutions for the stationary Cahn-Hillard-Navier-Stokes equations with singular potential

Technische Samenvatting: Zwakke Oplossingen voor de Stationaire Cahn-Hilliard-Navier-Stokes-vergelijkingen met Singulier Potentieel

Probleemstelling
Dit artikel onderzoekt het bestaan van zwakke oplossingen voor het stationaire Navier–Stokes–Cahn-Hilliard-systeem dat een driedimensionale, samendrukbare, tweefasige vloeistofmengsel met een diffuus interface beschrijft. Het systeem wordt geregeerd door de behoudswetten voor massa, impuls en concentratie, gekoppeld aan een vergelijking voor het chemisch potentieel. Een onderscheidend kenmerk van het model is de opname van een singuliere vrije-energiedichtheid van het Flory–Huggins-type (logaritmisch). Deze keuze zorgt ervoor dat de massafraction $c$ binnen het fysische bereik $(-1, 1)$ blijft en de aanwezigheid van vacuümtoestanden ($\rho = 0$) toelaat.

Het specifieke randwaardeprobleem omvat een begrensd domein $\Omega \subset \mathbb{R}^3$ met Dirichlet-randvoorwaarden voor de snelheid en Neumann-randvoorwaarden voor de concentratie en het chemisch potentieel. De voornaamste wiskundige uitdagingen vloeien voort uit twee bronnen:

1. De Stationaire Instelling: In tegenstelling tot evolutionaire problemen mist het stationaire systeem een a priori energie-ongelijkheid ($E(t) \leq E(0)$) gedreven door tijdsontwikkeling, wat doorgaans de nodige controle biedt over het singuliere potentieel en het concentratieveld.
2. Singuliere Potentieel en Vacuüm: De logaritmische singulariteit in de vrije energie verhindert standaard Taylor-expansie-regularisaties, en de potentiële degeneratie van de dichtheid nabij vacuümtoestanden bemoeilijkt de compactheidsargumenten die nodig zijn voor het nemen van de limiet.

Methodologie
De auteurs stellen het bestaan vast via een tweeledig limietproces toegepast op een zorgvuldig geconstrueerd benaderingssysteem.

1. Gespecialiseerde Regularisatie van het Potentieel:
   Om het gebrek aan $L^2$-grenzen voor $c$ en het falen van standaard Taylor-expansies aan te pakken, introduceren de auteurs een gespecialiseerde regularisatie $f^\delta_2(c)$ voor de singuliere logaritmische term. In tegenstelling tot eerdere benaderingen wordt deze benadering zo geconstrueerd dat het verschil $f^\delta_2(c) - \frac{\theta_c}{2}c^2$ slechts lineair groeit naarmate $|c| \to \infty$. Deze specifieke constructie is ontworpen om de kwadratische groei veroorzaakt door de anti-diffusie-term te annuleren, waardoor de nodige compactheid wordt hersteld en uniforme schattingen mogelijk worden gemaakt.

2. Tweeledig Benaderingsschema:

   • Niveau 1 ($\epsilon$-limiet): De auteurs introduceren eerst een kunstmatige viscositeit en regularisatieparameters ($\epsilon$) om een geregulariseerd systeem te construeren. Met behulp van een vast-puntargument (Schaefer's Vast-Puntstelling) bewijzen ze het bestaan van sterke oplossingen voor dit benaderingssysteem. Uniforme schattingen die onafhankelijk zijn van $\epsilon$ worden afgeleid, waardoor de limiet $\epsilon \to 0$ kan worden genomen, wat resulteert in een zwakke oplossing voor een tussenliggend systeem met de geregulariseerde vrije energie.
   • Niveau 2 ($\delta$-limiet): De tweede stap betreft het verwijderen van de regularisatie van het potentieel ($\delta \to 0$). Een grote moeilijkheid hier is het beheersen van de $L^{3/2}$-integreerbaarheid van de dichtheid. De auteurs maken gebruik van een specifieke keuze voor kunstmatige druk en een subtiel interpolatieargument. Dit argument maakt gebruik van de structuur van het geregulariseerde potentieel, met name de observatie dat $\frac{\partial f^\delta_2}{\partial c} - \theta_c c = O(\ln \frac{1}{\delta})$ voor $|c| > 1+\delta$. Dit stelt hen in staat de schattingen voor de dichtheidsnorm te sluiten en uniforme grenzen vast te stellen die onafhankelijk zijn van $\delta$.

3. Sterke Convergentie en Fysische Grenzen:
   Het nemen van de limiet $\delta \to 0$ vereist het bewijzen van de sterke convergentie van de dichtheid om de niet-monotoone drukterm te behandelen. Door technieken uit de theorie van stationaire samendrukbare Navier–Stokes-vergelijkingen aan te passen, stellen de auteurs de sterke convergentie van de dichtheid vast. Bovendien bewijzen ze dat de limietconcentratie $c$ strikt voldoet aan de fysische grenzen $-1 < c < 1$ bijna overal op de drager van de dichtheid ($\{\rho > 0\}$), waardoor de singuliere term goed gedefinieerd is.

Belangrijkste Resultaten
De hoofdstelling (Stelling 1.1) stelt het bestaan vast van een zwakke oplossing $(\rho, u, \mu, c)$ voor het stationaire probleem onder de volgende voorwaarden:

• Domein: Een begrensd, enkelvoudig samenhangend domein met een $C^2$-rand.
• Krachten: Externe krachten $g_1, g_2$ zijn begrensd.
• Adiabatische Exponent ($\gamma$):
  • $\gamma > 3/2$ als de externe kracht irrotatief is ($\nabla \times g_1 = 0$).
  • $\gamma > 5/3$ als de externe kracht rotatief is ($\nabla \times g_1 \neq 0$).
• Massabeperkingen: De totale massa $m_1 > 0$ en relatieve massa $m_2$ voldoen aan $m_2 \in (-m_1, m_1)$.

De oplossing voldoet aan de geregenormaliseerde continuïteitsvergelijking, de zwakke formulering van de impuls- en concentratievergelijkingen, en de energie-ongelijkheid. Cruciaal is dat de concentratie voldoet aan $-1 < c(x) < 1$ bijna overal waar $\rho(x) > 0$.

Betekenis en Claims
De auteurs claimen dat dit werk het eerste existentie-resultaat biedt voor zwakke oplossingen van het stationaire samendrukbare Navier–Stokes–Cahn-Hilliard-systeem dat gelijktijdig omvat:

1. Een singuliere vrije-energiepotentieel (Flory–Huggins-type), wat fysisch consistent is voor binaire mengsels maar wiskundig uitdagend is vanwege het gebrek aan globale grenzen voor de concentratie.
2. De aanwezigheid van vacuümtoestanden ($\rho = 0$), wat degeneratie in de vergelijkingen introduceert.

Het artikel merkt op dat hoewel het bereik van toelaatbare adiabatische exponenten $\gamma$ mogelijk uitbreidbaar is op basis van resultaten voor enkelvoudige vloeistof-stationaire samendrukbare Navier–Stokes-vergelijkingen, het huidige bereik wordt beperkt door de moeilijkheid om gewogen schattingen voor de concentratie nabij de rand te verkrijgen vanwege de singulariteit van de vrije energie. Het werk vormt een belangrijke stap in het overbruggen van de kloof tussen evolutionaire modellen (waar energie-ongelijkheden controle bieden) en stationaire modellen met singuliere potentialen.