Goal-Oriented Time Adaptivity for Linear Port-Hamiltonian Differential-Algebraic Equations of Index~1

Technische Samenvatting: Doelgerichte Tijdsadaptiviteit voor Lineaire Port-Hamiltoniaanse DAE's van Index 1

Probleemstelling
Port-Hamiltoniaanse (pH) systemen bieden een gestructureerd raamwerk voor het modelleren van fysieke systemen door energieopslag, dissipatie en interconnectie te coderen. Een fundamentele uitdaging bij de numerieke simulatie van deze systemen, met name over lange tijdsintervallen, is het behoud van de energiebalansvergelijking op het discrete niveau. Hoewel structuurbehoudende integratoren (bijvoorbeeld discrete gradiëntmethoden) een exacte energiebalans afdwingen, beperken ze vaak de keuze van tijdsintegratieschema's.

Dit werk behandelt lineaire port-Hamiltoniaanse differentiaal-algebraïsche vergelijkingen (DAE's) van index 1. De specifieke uitdaging ligt in het beheersen van de schending van de energiebalans zonder gebonden te zijn aan een specifieke integrator. De auteurs beogen een tijdsadaptieve strategie te ontwikkelen die de schending van de energiebalans behandelt als een "hoeveelheid van belang" (QoI), waarbij het tijdsrooster alleen wordt verfijnd waar nodig om deze schending te reduceren tot een voorgeschreven tolerantie.

Methodologie
De voorgestelde aanpak maakt gebruik van de Dual Weighted Residual (DWR) methode voor doelgerichte a posteriori foutenschatting. De methodologie verloopt via de volgende stappen:

1. Modelreductie en Formulering:
   De index-1 pH-DAE wordt getransformeerd naar een semi-expliciete vorm met behulp van een kern-splitsingspaar. Dit maakt eliminatie van algebraïsche variabelen mogelijk via een Schur-complementreductie, resulterend in een gereduceerde evolutievergelijking uitsluitend voor de differentiaalvariabelen. Deze reductie is cruciaal voor het formuleren van het variatieprobleem en het daaropvolgende geadjungeerde systeem.

2. Discretisatie:
   Het gereduceerde primaire probleem wordt gediscetiseerd met een discontinu Galerkin-methode van orde nul in de tijd (dG(0)), wat overeenkomt met de impliciete Euler-methode. Het discrete geadjungeerde probleem wordt afgeleid door de blokkenstructuur van het primaire systeem te transponeren, wat resulteert in een recursie achterwaarts in de tijd.

3. Doelfunctionaal en Geadjungeerd Probleem:
   De hoeveelheid van belang wordt gedefinieerd als het residu van de globale energiebalans, wat de som is van de kwadratische lokale schendingen van de energiebalans over de tijdsintervallen. Omdat deze functionaal quartisch is in de toestandsvariabele, omvat het geadjungeerde probleem springcondities op de tijdsknooppunten. Het continue geadjungeerde probleem is een achterwaarts in de tijd lopende ODE met stuksgewijs gladde forcing en specifieke springcondities afgeleid van de afgeleide van de doelfunctionaal.

4. Foutrepresentatie en Indicatoren:
   Met behulp van het DWR-raamwerk wordt een foutidentiteit afgeleid die de fout in de doelfunctionaal relateert aan het product van het primaire residu en de geadjungeerde fout. De auteurs leiden berekenbare lokale foutindicatoren af door de geadjungeerde fout te benaderen met een stuksgewijs lineair interpolant van de discrete geadjungeerde oplossing. Deze indicatoren combineren het lokale primaire residu (afwijking tussen de rechterkant en de discrete operator) met de geadjungeerde gevoeligheid (toestandsprongen gewogen met geadjungeerde verschillen).

5. Parallelle Geadjungeerde Benadering:
   Het sequentieel oplossen van het discrete geadjungeerde systeem is rekenkundig duur vanwege de achterwaartse koppeling. De auteurs stellen een paralleliseerbare Block-Jacobi-benadering voor. Door gebruik te maken van de dissipatieve structuur van het pH-systeem bewijzen ze dat de versterkingsmatrix van de Jacobi-iteratie een strikte contractie is in de energienorm (op voorwaarde dat het symmetrische deel van het Schur-complement coercief is). Dit maakt het ontkoppelen van intervallen mogelijk, waardoor parallelle berekening van de geadjungeerde benadering wordt mogelijk gemaakt.

Belangrijkste Bijdragen
Het artikel levert vier primaire bijdragen:

• Op Maat Gemaakte Adaptieve Systeem: Definitie van een doelgerichte tijdsadaptieve schema specifiek voor lineaire index-1 pH-DAE's, waarbij gebruik wordt gemaakt van Schur-complementeliminatie om volledig binnen de differentiaalvariabelen te werken.
• Foutrepresentatie: Afleiding van een op DWR gebaseerde a posteriori foutrepresentatie voor de energiebalansfout, inclusief de formulering van berekenbare lokale foutindicatoren die rekening houden met de quartische aard van de doelfunctionaal.
• Parallelle Geadjungeerde Oplosser: Voorstellen van een Block-Jacobi-benadering voor het achterwaartse geadjungeerde probleem. De auteurs bewijzen dat de dissipativiteit van de pH-structuur garandeert dat deze benadering een strikte contractie is, wat het gebruik rechtvaardigt voor efficiënte, geparalleliseerde berekening.
• Numerieke Validatie: Demonstratie van de doeltreffendheid van de methode door simulaties van elektrische schakelmodelen, waarbij aanzienlijke winst in efficiëntie wordt aangetoond ten opzichte van uniforme tijdstappen.

Resultaten
Numerieke experimenten werden uitgevoerd op twee benchmarks: een academisch 3-dimensionaal systeem en een transmissielijnmodel (RCL-ladder).

• Efficiëntie: De DWR-adaptieve strategie bereikte dezelfde nauwkeurigheid in de schending van de energiebalans als uniforme verfijning, maar met aanzienlijk minder tijdstappen. In de transmissielijn-benchmark reduceerde de methode het aantal vereiste intervallen met tot wel 89% voor specifieke foutdoelen.
• Convergentie: Zowel uniforme als adaptieve schema's vertoonden de verwachte asymptotische convergentiesnelheid van $O(N^{-3})$ voor de energiebalansfout, consistent met de discretisatie van de toestand van de eerste orde en de kwadratische doelfunctionaal. De adaptieve methode bereikte dit met een aanzienlijk gereduceerde foutconstante.
• Stabilisatie van Geadjungeerde: Het onderzoek naar de parallelle Block-Jacobi-iteratie onthulde dat de set intervallen die zijn gemarkeerd voor verfijning (de Dörfler-set) goed stabiliseert voordat de geadjungeerde oplossing volledig convergeert. Dit maakt vroegtijdige beëindiging van de Jacobi-iteraties mogelijk, wat verdere rekenkundige versnellingen biedt, met name in de vroege en tussenliggende fasen van de adaptieve lus.
• Effectiviteit: De foutenschatting bleek begrensd, met een effectiviteitsindex die doorgaans varieert tussen 2 en 5. De schatter overschat consistent de grootte van de fout, waardoor een conservatief stopcriterium wordt gewaarborgd.

Betekenis en Beweringen
Het artikel beweert dat de primaire betekenis ligt in het bieden van een rigoureus, doelgericht raamwerk voor het beheersen van schendingen van de energiebalans in pH-DAE's zonder de keuze van tijdsintegrator te beperken. Door de energiebalans te behandelen als een berekenbare hoeveelheid van belang, richt de methode rekenkundige middelen op de tijdsintervallen waar ze het meest effectief zijn.

De auteurs benadrukken dat hoewel de DWR-methode bekendstaat voor optimale besturing en andere toepassingen, de toepassing ervan op de specifieke structuur van pH-DAE's – met name het hanteren van index-1 beperkingen en de quartische energiefunctionaal – de specifieke reducties en geadjungeerde formuleringen vereist die hier worden gepresenteerd. Het werk benadrukt ook het praktische nut van de parallelle Block-Jacobi-benadering, die gebruik maakt van de fysieke dissipativiteit van het systeem om de geadjungeerde berekening te ontkoppelen, waardoor het adaptieve algoritme rekenkundig haalbaar wordt voor grotere systemen.

Het artikel concludeert met de opmerking dat hoewel het huidige werk zich richt op lineaire systemen en dG(0)-discretisatie, het theoretische raamwerk zich uitstrekt tot methoden van hogere orde en niet-lineaire systemen, hoewel de stabiliteitsanalyse voor de parallelle oplosser verdere generalisatie zou vereisen.