Evolving finite elements for advection diffusion with an evolving interface

Technische Samenvatting: Evoluerende Eindige Elementen voor Advektie-Diffusie met een Evoluerend Interface

Probleemstelling
Dit artikel behandelt de numerieke benadering van parabolische vergelijkingen gedefinieerd op een stationair domein $\Omega$ dat een bewegend interface $\Gamma(t)$ bevat. Het interface verdeelt het domein in twee subdomeinen, $\Omega_1(t)$ en $\Omega_2(t)$. De besturende vergelijkingen modelleren het transport en de diffusie van scalaire velden $u_1$ en $u_2$ binnen deze subdomeinen, onderworpen aan:

• Advektie-diffusie-reactiedynamiek binnen elk subdomein met tijdsafhankelijke coëfficiënten ($A_i, B_i, C_i$).
• Een Dirichlet-randvoorwaarde op de buitenrand $\partial\Omega$.
• Continuïteit van de oplossing over het interface ($u_1 = u_2$).
• Een voorgeschreven fluxsprong over het interface ($A_1 \partial_{\nu} u_1 - A_2 \partial_{\nu} u_2 = g$).

De belangrijkste geïdentificeerde uitdagingen zijn de evolutie van de subdomeinen in de tijd en de aanwezigheid van de fluxsprong over het bewegende interface. Hoewel toepassingen bestaan in de stromingsleer, materiaalkunde en celbiologie, richt het artikel zich op de wiskundige en numerieke formulering van het probleem.

Methodologie
De auteurs stellen een Arbitrary Lagrangian-Eulerian (ALE)-benadering voor die gebruikmaakt van Evoluerende Eindige Elementruimtes (ESFEM), gebaseerd op het abstracte raamwerk zoals vastgesteld in [21].

1. Functioneel-analytische Setting: Het continue probleem wordt geformuleerd binnen de context van evoluerende Hilbertruimtes. De auteurs definiëren compatibele paren ruimtes $(X(t), \phi_t)$ waarbij $\phi_t$ een stroomafbeelding is die functies transporteert van een referentiedomein naar het evoluerende domein. Dit maakt de definitie van zwakke materiële afgeleiden mogelijk en leidt tot een goedgestelde zwakke formulering (Stelling 2.16).
2. Meshconstructie:
   • Initieel Mesh: Het domein wordt initieel getrianguleerd met stuksgewijs lineaire simplexen. Cruciaal is dat de hoekpunten van interfacefacetten precies op het initiële interface $\Gamma(0)$ liggen (Aanname M5).
   • Isoparametrische Afbeelding: Om hogere-orde nauwkeurigheid te bereiken, wordt het lineaire mesh afgebeeld op een isoparametrisch mesh met behulp van een diffeomorfisme $\Psi_h$. Deze afbeelding duwt de interfacefacetten naar het exacte gebogen interface $\Gamma(0)$, terwijl de continuïteit van het mesh behouden blijft.
   • Tijds-evolutie: Het mesh evolueert door de Lagrangiaanse knopen te verplaatsen volgens een bekende, gladde globale snelheidsveld $w$. De discrete stroom $\Phi^h$ wordt geconstrueerd als een polynoominterpolatie van de continue stroom, zodat het mesh te allen tijde een conformerende toelaatbare onderverdeling blijft.
3. Lift-afbeelding: Een lift-operator $\Lambda_h$ wordt geïntroduceerd om functies af te beelden van het discrete computationele domein (dat het gebogen interface benadert) terug naar het exacte continue domein. Dit is essentieel voor het vergelijken van de discrete oplossing met de exacte oplossing in de foutanalyse.
4. Discreet Schema: De methode maakt gebruik van isoparametrische eindige elementen van orde $k$. Het discrete variatieprobleem wordt geformuleerd met behulp van de discrete materiële afgeleide en de geliftte coëfficiënten. Het schema wordt opgelost met een tijdstapmethode (specifiek Backward Difference Formulas, BDF, in de numerieke experimenten).

Belangrijkste Bijdragen

• Goedgesteldheid: Het artikel biedt een rigoureuze functioneel-analytische setting die het bestaan en de uniciteit van de oplossing van het continue probleem bewijst (Stelling 2.16), inclusief aanvullende regulariteitsresultaten onder specifieke aannames over de data.
• Isoparametrisch ESFEM-raamwerk: De auteurs breiden het abstracte raamwerk voor evoluerende eindige elementen uit om tweezijdige interfaceproblemen met fluxsprongen te behandelen. Ze beschrijven in detail de constructie van isoparametrische elementen die het gebogen interface nauwkeurig benaderen, waardoor convergentie van hogere orde mogelijk wordt.
• Optimale Foutgrenzen: Het belangrijkste theoretische resultaat (Stelling 5.1) stelt optimale foutgrenzen vast. Specifiek wordt de fout in de $L^2$-norm begrensd door $O(h^{k+1})$, en de fout in de $H^1$-semi-norm door $O(h^k)$, waarbij $h$ de meshgrootte is en $k$ de polynoomgraad. Deze grenzen gelden tot op een constante en komen overeen met de interpolatiefout van een bekende gladde oplossing.
• Implementatie en Verificatie: Het artikel bevat een robuuste implementatie met behulp van het Firedrake-pakket en GMSH. Numerieke experimenten in 2D en 3D verifiëren de theoretische convergentiesnelheden voor verschillende polynoomordes.

Resultaten
De numerieke experimenten tonen de effectiviteit van de methode aan op een modelprobleem met een bekende exacte oplossing die continu is maar niet differentieerbaar over het interface.

• Convergentie: Voor ruimtelijke discretisatie van orde $k$ bevestigen de numerieke resultaten convergentie van orde $k+1$ in de $L^2$-norm en orde $k$ in de $H^1$-semi-norm.
• Robuustheid: De resultaten gelden voor zowel 2D- als 3D-domeinen, zelfs wanneer de meshes onafhankelijk worden gegenereerd (in plaats van via uniforme verfijning), wat de stabiliteit van de methode aangeeft.
• Beschikbaarheid van Code: De auteurs verstrekken de simulatiecode om reproduceerbaarheid en verder gebruik door de gemeenschap mogelijk te maken.

Betekenis en Reikwijdte
Het artikel claimt een volledig theoretisch en praktisch raamwerk te bieden voor het oplossen van advektie-diffusieproblemen met bewegende interfaces met behulp van isoparametrische eindige elementen van hoge orde. De betekenis ligt in het overbruggen van de kloof tussen de theorie van abstracte evoluerende ruimtes en de praktische implementatie voor interfaceproblemen met fluxsprongen.

De auteurs houden een bescheiden reikwijdte aan met betrekking tot het snelheidsveld: ze veronderstellen dat een globale, gladde snelheidsveld a priori gegeven is. Ze stellen expliciet dat het werk niet ingaat op hoe zo'n snelheidsveld berekend moet worden (bijvoorbeeld via herparametrisatie om meshverstrengeling te voorkomen), en merken op dat dergelijke methoden in de literatuur bestaan [8, 18, 19]. Bovendien, hoewel de methode theoretisch van toepassing zou kunnen zijn op gekoppelde PDE-systemen waarbij domeinevolutie deel uitmaakt van de oplossing, bestrijkt de huidige analyse dergelijke gekoppelde problemen niet. De focus blijft strikt gericht op de benadering van het parabolische transmissieprobleem gegeven een bekende domeinevolutie.