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

Resumen Técnico: Adaptatividad Temporal Orientada a Objetivos para EDAE Lineales de Port-Hamiltoniano de Índice 1

Enunciado del Problema
Los sistemas de Port-Hamiltoniano (pH) proporcionan un marco estructurado para modelar sistemas físicos codificando el almacenamiento de energía, la disipación y la interconexión. Un desafío fundamental en la simulación numérica de estos sistemas, particularmente sobre horizontes temporales largos, es preservar la ecuación de balance de energía a nivel discreto. Aunque los integradores que preservan la estructura (por ejemplo, métodos de gradiente discreto) imponen un balance de energía exacto, a menudo restringen la elección de los esquemas de integración temporal.

Este trabajo aborda ecuaciones diferenciales-algebraicas (EDAE) lineales de Port-Hamiltoniano de índice 1. El desafío específico radica en controlar la violación del balance de energía sin estar restringido a un integrador específico. Los autores buscan desarrollar una estrategia adaptativa temporal que trate la violación del balance de energía como una "cantidad de interés" (QoI), refinando la malla temporal solo donde sea necesario para reducir dicha violación a una tolerancia prescrita.

Metodología
El enfoque propuesto utiliza el método del Residuo Ponderado Dual (DWR) para la estimación de error a posteriori orientada a objetivos. La metodología procede a través de los siguientes pasos:

1. Reducción de Modelo y Formulación:
   La EDAE de pH de índice 1 se transforma en una forma semi-explicita utilizando un par de descomposición del núcleo. Esto permite la eliminación de variables algebraicas mediante una reducción de complemento de Schur, resultando en una ecuación de evolución reducida para las variables diferenciales únicamente. Esta reducción es crucial para formular el problema variacional y el sistema adjunto subsiguiente.

2. Discretización:
   El problema primal reducido se discretiza utilizando un método de Galerkin discontinuo de orden cero en el tiempo (dG(0)), que corresponde al método de Euler implícito. El problema adjunto discreto se deriva transponiendo la estructura de bloques del sistema primal, resultando en una recurrencia hacia atrás en el tiempo.

3. Funcional Objetivo y Problema Adjunto:
   La cantidad de interés se define como el residuo global del balance de energía, que es la suma de las violaciones cuadráticas locales del balance de energía sobre los intervalos de tiempo. Debido a que este funcional es cuártico en la variable de estado, el problema adjunto implica condiciones de salto en los nodos temporales. El problema adjunto continuo es una EDO hacia atrás en el tiempo con forzamiento suave a trozos y condiciones de salto específicas derivadas de la derivada del funcional objetivo.

4. Representación del Error e Indicadores:
   Utilizando el marco DWR, se deriva una identidad de error que relaciona el error en el funcional objetivo con el producto del residuo primal y el error adjunto. Los autores derivan indicadores de error locales computables aproximando el error adjunto con un interpolante lineal a trozos de la solución adjunta discreta. Estos indicadores combinan el residuo primal local (desajuste entre el lado derecho y el operador discreto) con la sensibilidad adjunta (saltos de estado ponderados por las diferencias adjuntas).

5. Aproximación Adjunta Paralela:
   Resolver el sistema adjunto discreto secuencialmente es computacionalmente costoso debido al acoplamiento hacia atrás. Los autores proponen una aproximación Block-Jacobi paralelizable. Aprovechando la estructura disipativa del sistema pH, demuestran que la matriz de amplificación de la iteración de Jacobi es una contracción estricta en la norma de energía (siempre que la parte simétrica del complemento de Schur sea coerciva). Esto permite desacoplar los intervalos, habilitando el cálculo paralelo de la aproximación adjunta.

Contribuciones Clave
El artículo presenta cuatro contribuciones principales:

• Esquema Adaptativo Específico: Definición de un esquema adaptativo temporal orientado a objetivos específicamente para EDAE lineales de pH de índice 1, utilizando la eliminación de complemento de Schur para trabajar enteramente dentro de las variables diferenciales.
• Representación del Error: Derivación de una representación de error a posteriori basada en DWR para el error de balance de energía, incluyendo la formulación de indicadores de error locales computables que tienen en cuenta la naturaleza cuártica del funcional objetivo.
• Solver Adjunto Paralelo: Propuesta de una aproximación Block-Jacobi para el problema adjunto hacia atrás. Los autores demuestran que la disipatividad de la estructura pH asegura que esta aproximación sea una contracción estricta, justificando su uso para un cálculo eficiente y paralelizado.
• Validación Numérica: Demostración de la eficacia del método mediante simulaciones de modelos de circuitos eléctricos, mostrando ganancias significativas de eficiencia sobre la escalada temporal uniforme.

Resultados
Se realizaron experimentos numéricos en dos puntos de referencia: un sistema académico de 3 dimensiones y un modelo de línea de transmisión (escalera RCL).

• Eficiencia: La estrategia adaptativa DWR logró la misma precisión en la violación del balance de energía que el refinamiento uniforme, pero con significativamente menos pasos de tiempo. En el punto de referencia de la línea de transmisión, el método redujo el número de intervalos requeridos hasta en un 89% para objetivos de error específicos.
• Convergencia: Tanto los esquemas uniformes como adaptativos exhibieron la tasa de convergencia asintótica esperada de $O(N^{-3})$ para el error de balance de energía, consistente con la discretización de estado de primer orden y el funcional objetivo cuadrático. El método adaptativo logró esto con una constante de error sustancialmente reducida.
• Estabilización Adjunta: El estudio sobre la iteración Block-Jacobi paralela reveló que el conjunto de intervalos marcados para refinamiento (el conjunto de Dörfler) se estabiliza bien antes de que la solución adjunta converja completamente. Esto permite la terminación temprana de las iteraciones de Jacobi, proporcionando mayores aceleraciones computacionales, particularmente en las etapas inicial e intermedia del bucle adaptativo.
• Efectividad: Se encontró que el estimador de error está acotado, con un índice de efectividad que típicamente oscila entre 2 y 5. El estimador consistentemente sobreestima la magnitud del error, asegurando un criterio de parada conservador.

Significado y Afirmaciones
El artículo afirma que su significado principal radica en proporcionar un marco riguroso y orientado a objetivos para controlar las violaciones del balance de energía en EDAE de pH sin restringir la elección del integrador temporal. Al tratar el balance de energía como una cantidad de interés computable, el método dirige los recursos computacionales a los intervalos de tiempo donde son más efectivos.

Los autores enfatizan que, aunque el método DWR es conocido para el control óptimo y otras aplicaciones, su aplicación a la estructura específica de las EDAE de pH —particularmente el manejo de las restricciones de índice 1 y el funcional de energía cuártico— requiere las reducciones específicas y las formulaciones adjuntas presentadas aquí. El trabajo también destaca la utilidad práctica de la aproximación Block-Jacobi paralela, que explota la disipatividad física del sistema para desacoplar el cálculo adjunto, haciendo así viable computacionalmente el algoritmo adaptativo para sistemas más grandes.

El artículo concluye señalando que, aunque el trabajo actual se centra en sistemas lineales y discretización dG(0), el marco teórico se extiende a métodos de orden superior y sistemas no lineales, aunque el análisis de estabilidad para el solver paralelo requeriría una generalización adicional.