Stabilization of Parabolic Time-Varying PDEs using Certified Reduced-Order Receding Horizon Control

Technische Samenvatting: Stabilisatie van Paraboolse Tijdsvariabele PDE's met Gecertificeerde Gereduceerde Orde Receding Horizon Control

Probleemstelling
Het artikel behandelt de stabilisatie van lineaire, tijdsvariabele paraboolse partiële differentiaalvergelijkingen (PDE's) met behulp van eindig-dimensionale receding horizon control (RHC), ook bekend als Model Predictive Control (MPC). Een primaire uitdaging bij het toepassen van RHC op PDE's zijn de rekenkosten die gepaard gaan met het herhaaldelijk oplossen van hoog-dimensionale open-loop optimale controleproblemen over overlappende tijdsintervallen. Hoewel gereduceerde orde-modellen (ROM's) een weg bieden om deze berekeningen te versnellen, kan het rechtstreeks toepassen van uit ROM's afgeleide besturingen op het volledig-orde model (FOM) de gesloten-lus stabiliteit en suboptimaliteitsgaranties van het oorspronkelijke systeem in gevaar brengen. Bovendien behandelt het artikel niet-gladde kostenfunctionalen, specifiek die welke de gekwadrateerde $\ell_1$-norm omvatten om sparsiteit in de besturingsinput af te dwingen, wat complexiteit toevoegt aan de optimalisatie en foutanalyse.

Methodologie
De auteurs stellen een tweestapskader voor dat een rigoureuze stabiliteitsanalyse voor het volledig-orde systeem combineert met een gecertificeerde modelreductiestrategie.

1. Analyse van het Volledig-Orde Model (FOM):
   De auteurs stellen eerst de exponentiële stabiliteit en suboptimaliteit vast van een continue-tijd RHC-scheme voor tijdsvariabele paraboolse systemen in Hilbertruimten. In tegenstelling tot eerdere werken die vertrouwen op terminale kosten of beperkingen om stabiliteit te waarborgen, maakt dit kader gebruik van een specifiek concateneringsschema en een observabiliteitsongelijkheid op het eindtijdstip. Deze aanpak maakt het mogelijk om zwakkere stapskosten (specifiek een $L^2(\Omega)$-trackingterm) en niet-gladde besturingskosten (inclusief sparsiteitsbevorderende regularisatie) te gebruiken zonder vereisten voor observabiliteitsschattingen op het begintijdstip.

2. Gereduceerde Orde Modellering (ROM) en Foutschatting:
   Om de rekenkosten aan te pakken, maken de auteurs gebruik van projectiegebaseerde Galerkin-modelordereeductie, specifiek met behulp van Proper Orthogonal Decomposition (POD). Er wordt een rigoureuze a posteriori foutanalyse ontwikkeld voor de optimale controleproblemen met eindige horizon. Dit omvat het afleiden van foutschattingen voor de toestand, de toegevoegde toestand, de optimale besturing en de waardenfunctie. Cruciaal is dat deze schattingen rekening houden met niet-gladde convexe besturingskosten en onnauwkeurige beginwaarden.

3. Gecertificeerd ROM-RHC Algoritme:
   De kernbijdrage is een adaptief ROM-RHC-algoritme (Algoritme 2). In plaats van gereduceerde besturingen blind toe te passen, gebruikt het algoritme de afgeleide a posteriori foutschattingen om online een "Ontspannen Dynamisch Programmeren" (RDP) voorwaarde te verifiëren.

   • Op elk tijdstap lost het algoritme het gereduceerde probleem op en berekent het prestatiegrenzen ($\bar{\alpha}^r_k$) op basis van de ROM-waardenfunctie en foutschattingen.
   • Als de prestatiegrens een door de gebruiker gespecificeerde drempel $\tilde{\alpha}$ voldoet, wordt de gereduceerde besturing toegepast op het volledige systeem.
   • Als aan de voorwaarde niet wordt voldaan, wordt de gereduceerde basis verfijnd (bijvoorbeeld door nieuwe POD-snapshots van de volledig-orde oplossing toe te voegen) en wordt de stap herhaald.
     Dit adaptieve proces zorgt ervoor dat de gereduceerde besturing het volledig-orde systeem stabiliseert met gegarandeerde prestatiegrenzen ten opzichte van het oneindige-horizon FOM-probleem.

Belangrijkste Bijdragen
Het artikel levert vier primaire bijdragen:

1. Stabiliteitsanalyse voor FOM: Het stelt exponentiële stabiliteit en suboptimaliteit vast voor continue-tijd RHC van tijdsvariabele lineaire paraboolse systemen met behulp van een $L^2(\Omega)$-trackingterm en niet-gladde kosten. Deze analyse berust op een observabiliteitsongelijkheid op het eindtijdstip in plaats van op het begintijdstip, en biedt een nieuwe expliciete bovengrens voor de suboptimaliteitsparameter.
2. Rigoureuze Foutschattingen: De auteurs leiden a posteriori foutschattingen af voor de ROM-benadering van de waardenfunctie en de kostenfunctionaal. Deze schattingen zijn van toepassing op algemene niet-gladde convexe besturingskosten en onnauwkeurige beginwaarden, en ze vestigen interpolatie-eigenschappen en asymptotische convergentie van Galerkin-ROM's onder versoepelde regulariteitsaannames.
3. Gecertificeerd ROM-RHC Schema: Er wordt een nieuw algoritme voorgesteld dat deze foutschattingen combineert met een variant van het principe van Ontspannen Dynamisch Programmeren. Dit schema garandeert exponentiële stabilisatie van het FOM en suboptimaliteit ten opzichte van het oneindige-horizon FOM-probleem. De methode zorgt ervoor dat de ROM zo kan worden geconstrueerd dat de RDP-ongelijkheid wordt voldaan met een eindig-dimensionale basis.
4. Numerieke Validatie: De aanpak wordt gevalideerd op systemen met exponentieel instabiele vrije dynamiek en niet-gladde $\ell_1$-regularisatie. De resultaten tonen aan dat ROM-RHC aanzienlijke rekenversnellingen bereikt terwijl de stabiliteits- en suboptimaliteitsgaranties van de volledig-orde methode behouden blijven.

Resultaten
Numerieke experimenten op een 2D paraboolse PDE met 13 actuators en een niet-gladde kostenfunctionaal bevestigen de theoretische bevindingen:

• Stabiliteit: Zowel FOM- als ROM-RHC-schema's stabiliseren het systeem succesvol met exponentiële vervaltempo's die consistent zijn met de theoretische grenzen.
• Nauwkeurigheid: De gereduceerde-orde oplossingen vertonen extreem lage relatieve fouten ten opzichte van de volledig-orde oplossingen (bijvoorbeeld relatieve fouten in de kostenfunctie in de orde van $10^{-7}$ tot $10^{-8}$).
• Efficiëntie: De ROM-RHC-schema's bereiken een versnelling van ongeveer 10–13 keer in rekentijd ten opzichte van FOM-RHC. Dit wordt bereikt door het aantal volledig-orde gradiëntevaluaties met een factor van ongeveer 20 te verminderen, aangezien volledig-orde oplossingen alleen worden geactiveerd wanneer het model verfijning vereist.
• Foutschatting: De a posteriori foutschattingen blijken betrouwbaar, hoewel ze de werkelijke fout de neiging hebben te overschatten, met name voor grote voorspellingshorizons en kleine regularisatieparameters. De schattingen vangen echter correct de exponentiële afname van de fout naarmate de dimensie van de gereduceerde basis toeneemt.

Betekenis en Beweringen
Het artikel beweert een robuust kader te bieden voor de efficiënte en stabiele besturing van tijdsvariabele paraboolse PDE's. Door rigoureuze a posteriori foutanalyse te integreren met adaptieve modelreductie, overwint de methode de rekenonuitvoerbaarheid van standaard RHC voor hoog-dimensionale PDE's zonder theoretische garanties op te offeren. De auteurs benadrukken dat hun aanpak geen terminale kosten of beperkingen vereist om stabiliteit te waarborgen, waardoor het toepasbaar is op een bredere klasse van problemen, inclusief die met niet-gladde doelstellingen. Het certificeringsmechanisme zorgt ervoor dat de gereduceerde-orde controller niet slechts een benadering is, maar een geverifieerde stabilisator voor het volledig-orde systeem, met kwantificeerbare prestatiegrenzen. Het werk wordt gepresenteerd als een stap richting het behandelen van complexere systemen, met toekomstig werk voorgesteld voor niet-lineaire paraboolse systemen, parametrische onzekerheden en adaptieve horizonselectie.