Extrapolation of solvability of the parabolic $L^p$ Neumann problem on bounded Lipschitz cylinders

1. Probleemstelling en Context

Het artikel richt zich op de parabool $L^p$-Neumann-probleem voor partiële differentiaalvergelijkingen (PDE's) van het type:
$$Lu = -\partial_t u + \text{div}(A\nabla u) = 0$$
op een domein $\Omega = O \times \mathbb{R}$, waarbij $O$ een begrensde Lipschitz-domein is in de ruimtelijke variabelen. De tijdvariabele loopt van $-\infty$ tot $+\infty$. De coëfficiëntenmatrix $A$ wordt verondersteld te zijn reëel, begrensd, meetbaar en voldoet aan een uniforme ellipticiteitsvoorwaarde, maar mag afhankelijk zijn van de tijd.

Het centrale vraagstuk:
Een eerdere studie (Dindoš, Li, Pipher) had bewezen dat voor ongebonden grafische domeinen ($\Omega = \{(x', x_n) : x_n > \phi(x')\} \times \mathbb{R}$) de oplosbaarheid van het $L^p$-Neumann-probleem kan worden "geëxtrapoleerd". Dit betekent: als het probleem oplosbaar is voor een zekere $p \in (1, \infty)$ (en het bijbehorende $L^{p'}$-Dirichlet-probleem voor de geadjungeerde operator $L^*$), dan is het ook oplosbaar voor alle $q$ met $1 < q < p$.

Deze paper vult een belangrijke lacune aan door te onderzoeken of deze extrapolatie ook geldt voor begrensde Lipschitz-cilinders ($O \times \mathbb{R}$). Het blijkt dat de bewijstechnieken voor ongebonden domeinen hier niet direct toepasbaar zijn, omdat de geometrie en de randvoorwaarden fundamenteel anders zijn. De auteurs bewijzen dat de extrapolatie wel degelijk geldt voor dit geval, zonder extra aannames over de kleinheid van de Lipschitz-norm of de coëfficiënten (zoals Carleson-maatvoorwaarden) die in eerdere werken vaak nodig waren.

2. Methodologie en Bewijsstrategie

De kern van het bewijs is het aantonen van de oplosbaarheid voor de randwaarde $p=1$ in de zin van de atomische Hardy-ruimte $\mathcal{H}^1_{\text{ato}}(\partial\Omega)$. Als dit lukt, volgt de oplosbaarheid voor alle $1 < q < p$ door middel van reële interpolatie.

De bewijsstrategie verloopt als volgt:

1. Reductie tot Atomen:
   Het doel is om te laten zien dat voor elke $L^\infty$-atom $g$ (een functie met draagverzameling in een kleine parabolische kubus, gemiddelde nul en begrensd door de inverse van het volume) de oplossing $u$ voldoet aan:
   $$\|\tilde{N}(\nabla u)\|_{L^1(\partial\Omega)} \leq C$$
   waarbij $\tilde{N}$ de gemodificeerde niet-tangentiële maximale functie is.

2. Scheiding van Gevallen:
   De auteurs onderscheiden twee gevallen gebaseerd op de grootte $r$ van de atoom-draagverzameling ten opzichte van de diameter van het domein $O$:

   • Grote Atomen ($r \gtrsim 1$): Hier wordt gebruikgemaakt van de bekende $L^p$-oplosbaarheid en exponentiële afname van de oplossing in de tijd (via een variant van de Poincaré-ongelijkheid en Caccioppoli-ongelijkheden). Omdat het Neumann-gegeven nul is voor grote $t$, vervaagt de oplossing exponentieel.
   • Kleine Atomen ($r < 1$): Dit is het technisch meest uitdagende deel. De auteurs introduceren een hulpdomein.

3. De "Unbounded Domain" Truc (Hulpdomein):
   Voor kleine atomen wordt het lokale probleem op het begrensde domein $\Omega$ gekoppeld aan een probleem op een ongebonden grafisch domein $\tilde{\Omega} = \{(x', x_n) : x_n > \tilde{\phi}(x')\} \times \mathbb{R}$.

   • De coëfficiënten $A$ en het domein worden lokaal uitgebreid tot dit ongebonden domein.
   • Er wordt een oplossing $w$ geconstrueerd op $\tilde{\Omega}$ met dezelfde randvoorwaarde $g$.
   • Omdat $u$ (op $\Omega$) en $w$ (op $\tilde{\Omega}$) dezelfde PDE en randvoorwaarde delen op het ondersteunende gebied, is het verschil $u-w$ een oplossing met nul Neumann-data op dat gebied.
   • Door de bekende resultaten voor ongebonden domeinen (uit [7]) toe te passen op $w$, en de Caccioppoli-ongelijkheid op $u-w$, kunnen de auteurs de $L^1$-norm van de gradient van $u$ controleren.

4. Reflectie en Schattingen:
   Een cruciale stap is het gebruik van een even reflectie van de oplossing over de rand van het ongebonden domein (buiten de steun van de atoom). Dit maakt het mogelijk om de oplossing te behandelen als een oplossing van een parabolische vergelijking op een groter domein, waarbij de maximumprincipe en interne regulariteitseigenschappen kunnen worden gebruikt om afnemende gedrag (decay) van de oplossing en zijn gradient op grote afstand te bewijzen.

5. Absorptie van Termen:
   Door de schattingen te combineren en een geschikte parameter $K$ (gerelateerd aan de tijdsduur van de atoom) te kiezen, kunnen de auteurs een "sluitende" ongelijkheid aantonen waarbij de rechterkant kan worden geabsorbeerd door de linkerkant, wat leidt tot de gewenste uniforme $L^1$-schatting.

3. Belangrijkste Resultaten

• Hoofdstelling (Theorem 1.3): Laat $\Omega = O \times \mathbb{R}$ met $O$ een begrensde Lipschitz-domein. Als het $L^p$-Neumann-probleem voor $L$ en het $L^{p'}$-Dirichlet-probleem voor de geadjungeerde operator $L^*$ oplosbaar zijn voor een zekere $p \in (1, \infty)$, dan is het $L^q$-Neumann-probleem oplosbaar voor alle $1 < q < p$.
• Hardy-ruimte Oplosbaarheid: Het probleem is ook oplosbaar voor Neumann-data in de atomische Hardy-ruimte $\mathcal{H}^1_{\text{ato}}(\partial\Omega)$ (het geval $p=1$).
• Geen Extra Aannames: Het resultaat vereist geen kleineheid van de Lipschitz-norm van het domein en geen Carleson-maatvoorwaarden voor de coëfficiënten $A$. Dit is een significant vooruitgang ten opzichte van eerdere resultaten die dergelijke beperkende aannames nodig hadden voor begrensde domeinen.

4. Technische Bijdragen

• Overbrugging van Geometrie: De paper lost het probleem op dat de eerdere bewijzen voor ongebonden grafische domeinen niet direct generaliseerbaar waren naar begrensde cilinders. De introductie van het lokale ongebonden hulpdomein $\tilde{\Omega}$ is een innovatieve methode om lokale eigenschappen te koppelen aan globale theorie.
• Extrapolatie zonder Interpolatie naar $L^\infty$: Voor parabolische Dirichlet-problemen wordt vaak geïnterpoleerd tussen $L^p$ en $L^\infty$ (via het maximumprincipe). Voor Neumann-problemen is dit niet direct mogelijk. De auteurs gebruiken in plaats daarvan interpolatie tussen $L^p$ en de Hardy-ruimte $\mathcal{H}^1$, wat een "schoner" resultaat oplevert zonder extra aannames.
• Uitgebreide Regulariteitstechnieken: Het bewijs maakt gebruik van een combinatie van Caccioppoli-ongelijkheden, niet-tangentiële maximale functies, Sobolev-embeddings en de theorie van fundamentele oplossingen (Aronson-schattingen) op een zeer verfijnde manier.

5. Significatie

Deze paper sluit een belangrijke kennislacune in de theorie van parabolische randwaardeproblemen.

• Het complementeert de theorie voor elliptische en parabolische Dirichlet- en Regulariteitsproblemen.
• Het bevestigt dat de structuur van de oplossing (oplosbaarheid voor een interval van $p$) fundamenteel wordt bepaald door de ellipticiteit en de Lipschitz-geometrie, en niet door de grootte van het domein of kleine coëfficiënt-variaties.
• Het resultaat is van toepassing op een breed scala aan fysische modellen die worden beschreven door parabolische vergelijkingen op complexe, begrensde gebieden, en biedt een robuust theoretisch fundament voor numerieke methoden en verdere analyse op deze domeinen.

Kortom, Dindoš en Liu leveren een bewijs dat de oplosbaarheid van het parabolische Neumann-probleem op begrensde Lipschitz-cilinders even robuust is als op ongebonden grafische domeinen, mits de basisvoorwaarden voor ellipticiteit en de geadjungeerde Dirichlet-probleem zijn vervuld.