A Nonlocal $p$-Laplacian Interface Model with Sharp Interface

Technische Samenvatting: Een niet-lokaal p-Laplace interface model met een scherpe interface

Probleemstelling
Het artikel behandelt de uitdaging van het formuleren van interfaceproblemen binnen een niet-lokaal kader. Klassieke interfaceproblemen betreffen partiële differentiaalvergelijkingen (PDE's) die gedefinieerd zijn op een domein $\Omega$ dat wordt verdeeld door een interface $\Gamma$, waarbij de oplossing of de afgeleiden ervan discontinuïteiten vertonen. Specifiek beschouwen de auteurs de uitbreiding van het Poisson interfaceprobleem en het $p$-Laplace interfaceprobleem met membraancondities naar niet-lokale settings. Een primaire moeilijkheid in niet-lokale formuleringen is dat klassieke lokale interfacecondities (zoals continuïteit van de oplossing en voorgeschreven sprongen in de normale flux) niet direct kunnen worden opgelegd. Bestaande benaderingen benaderen de scherpe interface $\Gamma$ vaak met een volumetrische laag van dikte $\delta$, wat het modelleren van onafhankelijke diffusie in subdomeinen en het afhandelen van randgegevens bemoeilijkt.

Methodologie
De auteurs stellen een nieuw energiegebaseerd niet-lokaal benaderingsmodel voor dat een scherpe interface behoudt. In plaats van de geometrie van het domein te wijzigen, gebruikt het model specifieke energiefunctionalen waarbij de interfacecondities impliciet worden gecodeerd.

1. Energiefunctie-formulering:

   • Voor het lineaire geval ($p=2$) wordt de lokale energiefunctie $E(u)$ benaderd door een niet-lokale functie $E_\delta(u)$. De niet-lokale diffusieterm wordt geconstrueerd met behulp van een kernel $\Lambda(x,y)$ die interacties beperkt tot binnen de subdomeinen $\Omega_1$ en $\Omega_2$, wat onafhankelijke niet-lokale diffusie in elke regio waarborgt.
   • De Dirichlet interfaceconditie (sprong in flux) wordt afgehandeld via een strafterm die betrokken is bij niet-lokale gemiddelden ($K_{\delta,i} * u$) langs de interface $\Gamma$.
   • De Neumann randvoorwaarde wordt op natuurlijke wijze geïncorporeerd via de variationele structuur.
   • Het raamwerk wordt uitgebreid naar het niet-lineaire $p$-Laplace geval met membraancondities, waarbij de sprong in de oplossing over de interface wordt bestraft. De resulterende niet-lokale functie $F_\delta(u)$ heeft dezelfde structurele vorm als $E_\delta(u)$.

2. Theoretische Analyse:

   • Bestaan en Uniciteit: De auteurs stellen het bestaan en de uniciteit van minimizers voor de niet-lokale functies $E_\delta$ en $F_\delta$ vast met behulp van de directe methode van de calculus van variaties. Dit steunt op een nieuw afgeleide niet-lokale Poincaré-ongelijkheid die rekening houdt met de interface-sprongstrafterm.
   • $\Gamma$-convergentie: De kern van de theoretische bijdrage is het bewijs van $\Gamma$-convergentie. De auteurs tonen aan dat wanneer de niet-lokale horizon $\delta$ verdwijnt ($\delta \to 0$), de niet-lokale functies $\Gamma$-convergeren naar hun lokale tegenhangers ($E$ en $F$).
   • Compactheid: Gecombineerd met compactheidsresultaten voor niet-lokale $p$-Laplace operatoren, impliceert de $\Gamma$-convergentie dat de minimizers van de niet-lokale problemen convergeren naar de minimizers van de lokale interfaceproblemen in de $L^p$-norm.

3. Numerieke Implementatie:

   • De auteurs maken gebruik van een efficiënte eindige elementenmethode gebaseerd op de tensorproductstructuur van Gaussische kernels en multi-kubische polynoombasissen.
   • Deze aanpak ontbindt hoog-dimensionale integralen in producten van één-dimensionale integralen, wat de computationele kosten aanzienlijk vermindert.
   • Speciale kwadratuurregels (Gauss-Legendre) zijn afgestemd om de randlagen te behandelen die voortvloeien uit de kubische basisfuncties nabij cel-interfaces.

Belangrijkste Resultaten

• Behoud van de Scherpe Interface: In tegen tegenovergaan van eerdere volume-constraint benaderingen, behoudt het voorgestelde model een scherpe interface, wat onafhankelijke niet-lokale diffusie in disjuncte subdomeinen mogelijk maakt zonder de interfacegegevens in de bulk uit te breiden.
• Convergentiebewijs: Rigoureuze bewijzen bevestigen dat de minimizers van de niet-lokale functies convergeren naar de lokale minimizers wanneer $\delta \to 0$. De convergentiesnelheid wordt waargenomen als eerste-orde.
• Numerieke Validatie: Numerieke experimenten in één en twee dimensies verifiëren de theoretische convergentie. De $L^2$-fout tussen de niet-lokale oplossing en de grondwaarheid lokale oplossing neemt lineair af met $\delta$, en de resultaten zijn robuust over verschillende mesh-resoluties.

Betekenis en Claims
Het artikel claimt dat het voorgestelde scherpe-interface niet-lokale model een natuurlijkere en flexibelere methode biedt voor interfaceproblemen vergeleken met bestaande volumetrische benaderingen. Door het probleem via energiefuncties te formuleren, worden de moeilijk te behandelen normale flux interfacecondities impliciet voldaan. De auteurs benadrukken dat deze aanpak:

1. De uitbreiding naar complexere interfaceproblemen faciliteert, zoals aangetoond door de succesvolle toepassing op de niet-lineaire $p$-Laplace met membraancondities.
2. Het voorkomen van de noodzaak om interfacegegevens (bijv. $g$) van $L^2(\Gamma)$ naar $H^1(\Omega)$ uit te breiden.
3. Een wiskundig rigoureuze link biedt tussen niet-lokale en lokale interfacemodellen via $\Gamma$-convergentie, wat garandeert dat het niet-lokale model een geldige benadering is van het klassieke fysieke probleem naarmate de horizon verdwijnt.

Het werk stelt geen nieuwe fysieke toepassingen voor, maar biedt een robuust wiskundig en numeriek instrument voor het modelleren van bestaande interfacefenomenen in velden zoals vloeistofdynamica en materiaalkunde waar niet-lokale effecten relevant zijn.