Optimal error estimates for a discontinuous Galerkin method on curved boundaries with polygonal meshes

Technische Samenvatting: Optimale Foutschattingen voor een Discontinuous Galerkin-methode op Gekromde Grenzen met Polygonaal Mesh

Probleemstelling
Het artikel behandelt de numerieke oplossing van randwaardeproblemen, specifs twee-dimensionale lineaire advectie-diffusie-reactievergelijkingen, op domeinen met willekeurige gladde gekromde grenzen ($\Omega$). Een primaire uitdaging in deze context is het verlies van convergentieorde dat optreedt wanneer het fysieke domein wordt benaderd door een polygonaal mesh ($\Omega_h$). Standaard Discontinuous Galerkin (DG) methoden, wanneer deze worden toegepast op dergelijke polygonaal benaderingen waarbij randvoorwaarden direct op de computationele grens ($\partial\Omega_h$) worden opgelegd, lijden over het algemeen aan suboptimale convergentie. De geometrische fout die wordt geïntroduceerd door de mismatch tussen de gekromde fysieke grens en de stuksgewijs lineaire computationele grens kan de oplossing vervuilen, waardoor de nauwkeurigheid beperkt wordt tot maximaal tweede orde, ongeacht de gebruikte polynoomgraad.

Methodologie
De auteurs stellen de DG–ROD (Discontinuous Galerkin – Reconstruction for Off-site Data) methode voor en analyseren deze. Deze aanpak breidt een strategie uit die eerder is ontwikkeld voor Finite Element (FE) en Finite Volume (FV) methoden naar het DG-raamwerk.

• Mesh en Ruimte: Het fysieke domein wordt gediscretiseerd met behulp van standaard polygonaal (triangulair) mesh. De methode maakt gebruik van gebroken Sobolev-ruimten en discontinue stuksgewijze polynoomruimten ($S_{h,N}$).
• Behandeling van Randvoorwaarden: In plaats van de voorgeschreven randvoorwaarden direct op de polygonaal grens $\partial\Omega_h$ op te leggen, maakt de DG–ROD methode gebruik van een polynoomreconstructie. Voor elk randelement wordt een geconstrueerde minimalisatieprobleem opgelost om een lokale polynoom te bepalen die de afstand tot de numerieke oplossing minimaliseert, terwijl deze voldoet aan de exacte randvoorwaarden op specifieke punten op de ware fysieke grens $\partial\Omega$.
• Variatierepresentatie: De methode combineert het Symmetric Interior Penalty (SIP) schema voor diffusie met een upwind-schema voor advectie-reactietermen. De gereconstrueerde randgegevens ($g_{ROD}$) worden direct in de primaire variatierepresentatie opgenomen, wat effectief de geometrische mismatch binnen de zwakke vorm corrigeert.
• Implementatie: De auteurs hanteren een globale formulering waarbij de reconstructie en de DG-solver in één enkele stap worden opgelost, in plaats van een iteratieve alternerende aanpak.

Belangrijkste Bijdragen en Theoretische Analyse
Het artikel biedt een rigoureuze theoretische analyse van de DG–ROD methode en stelt de existentie, uniciteit en optimale foutschattingen vast voor zowel convexe als niet-convexe domeinen.

1. Welgesteldheid: De auteurs bewijzen de begrensdheid van de bilineaire vorm en stellen een inf-sup conditie vast (via de Banach–Nečas–Babuška stelling), wat de existentie en uniciteit van de discrete oplossing garandeert voor een voldoende kleine meshgrootte $h$ en een voldoende grote strafparameter $\eta$.
2. Convexe Domeinen: Voor convexe domeinen bewijzen de auteurs dat de Galerkin-orthogonaliteit standhoudt. Onder geschikte regulariteitsveronderstellingen op de exacte oplossing leiden zij af:
   • Optimale convergentie van orde $O(h^N)$ in de DG-norm (energienorm).
   • Optimale convergentie van orde $O(h^{N+1})$ in de $L^2$-norm met behulp van dualiteitsargumenten.
3. Niet-convexe Domeinen: In niet-convexe omgevingen verbreekt de geometrische benadering de standaard Galerkin-orthogonaliteit. De auteurs introduceren een residu-term die voortvloeit uit de extensie van de oplossing naar een groter glad domein $\tilde{\Omega}$. Zij leiden foutschattingen af die een residu-term bevatten die afhankelijk is van de geometrische mismatch.
   • Theoretische analyse suggereert dat voor niet-convexe domeinen met $N > 2$, optimale convergentie in de $L^2$-norm mogelijk niet wordt bereikt, tenzij het computationele domein het fysieke domein met hogere nauwkeurigheid benadert (orde $h^q, q>2$) of de oplossing een regelmatige extensie toelaat die aan specifieke voorwaarden voldoet.
   • De analyse benadrukt dat het verlies van orthogonaliteit de belangrijkste technische uitdaging vormt in het niet-convexe geval.

Numerieke Resultaten
De theoretische bevindingen worden gevalideerd door numerieke experimenten op drie benchmark-geometrieën: een cirkel (convex), een annulus (niet-convex) en een rozenvormig domein (niet-convex).

• Vergelijking: De resultaten vergelijken de klassieke DG-methode (het opleggen van condities op $\partial\Omega_h$) met de DG–ROD methode.
• Prestaties: De klassieke DG-methode vertoont consequent een verslechtering van de nauwkeurigheid, waarbij de convergentiesnelheden beperkt blijven tot ongeveer tweede orde, ongeacht de gebruikte polynoomgraad $N$.
• Effectiviteit van DG–ROD: De DG–ROD methode slaagt er succesvol in om optimale convergentiesnelheden ($O(h^{N+1})$ in $L^2$) te herstellen voor $N=2, 3, 4$ over alle geteste geometrieën en diverse advectie-reactiecoëfficiënten. Opmerkelijk genoeg, zelfs in het niet-convexe rozenvormige domein, bereikt de methode optimale snelheden, wat de auteurs toeschrijven aan het gebruik van gladde 'manufactured solutions' die regelmatige extensies mogelijk maken, waardoor de theoretische worst-case beperkingen voor algemene niet-convexe gevallen effectief worden omzeild.

Betekenis en Claims
Het artikel claimt de eerste rigoureuze theoretische analyse van de DG–ROD methode te bieden, waarbij het theoretische raamwerk van Finite Elements wordt uitgebreid naar Discontinuous Galerkin-methoden.

• Herstel van Nauwkeurigheid: De primaire betekenis is de demonstratie dat optimale convergentiesnelheden kunnen worden hersteld op eenvoudige polygonaal meshes zonder de noodzaak van complexe gekromde elementen of isoparametrische mapping, mits de randgegevens worden gereconstrueerd via de ROD-techniek.
• Robuustheid: De numerieke resultaten bevestigen de robuustheid van de methode bij het afhandelen van zowel convexe als niet-convexe domeinen met gekromde grenzen.
• Bescheidenheid over Beperkingen: De auteurs erkennen dat hoewel de theoretische foutschattingen voor niet-convexe domeinen potentiële sub-optimaliteit suggereren voor $N > 2$ zonder specifieke regulariteitsvoorwaarden, de numerieke benchmarks met gladde oplossingen laten zien dat de methode in de praktijk optimaal presteert. Zij merken op dat uitbreidingen naar niet-lineaire vergelijkingen, tijdsafhankelijke problemen en andere randvoorwaarden (Neumann, Robin) openstaande uitdagingen blijven voor toekomstig onderzoek. Het werk suggereert ook dat de aanpak aanpasbaar is aan drie-dimensionale problemen, verwijzend naar gerelateerd werk in klassieke Lagrange finite elements.