A domain decomposition strategy for natural imposition of mixed boundary conditions in port-Hamiltonian systems

Resumen Técnico: Una Estrategia de Descomposición de Dominios para la Imposición Natural de Condiciones de Contorno Mixtas en Sistemas Port-Hamiltonianos

Planteamiento del Problema
El artículo aborda el desafío de imponer condiciones de contorno mixtas (condiciones de Dirichlet y Neumann simultáneas) en sistemas port-hamiltonianos (pH) hiperbólicos mediante métodos de elementos finitos. En las formulaciones estándar de elementos finitos, las condiciones de contorno se imponen típicamente de forma fuerte (mediante la restricción de los espacios de funciones discretas) o débil (a través de términos de contorno naturales en la formulación variacional). La imposición fuerte suele derivar en ecuaciones diferencial-algebraicas (DAE), que son computacionalmente más difíciles de resolver que las ecuaciones diferenciales ordinarias (ODE). Por el contrario, la imposición débil suele requerir la introducción de multiplicadores de Lagrange para imponer las condiciones de contorno esenciales, lo que aumenta la dimensión y la complejidad del sistema. Los autores señalan que, si bien el marco port-hamiltoniano gestiona naturalmente los flujos de contorno mediante la estructura Stokes-Dirac, los esquemas de discretización existentes a menudo tienen dificultades para incorporar condiciones mixtas explícitamente sin recurrir a multiplicadores de Lagrange, preservando al mismo tiempo la estructura geométrica subyacente.

Metodología
La estrategia propuesta emplea un enfoque de descomposición de dominios combinado con Cálculo Exterior de Elementos Finitos (FEEC) y una formulación primal-dual. La metodología central consiste en los siguientes pasos:

1. Descomposición de Dominios: El dominio espacial $\Omega$ se divide en dos subdominios no superpuestos, $\Omega_1$ y $\Omega_2$, separados por una interfaz interna arbitraria $\Gamma_{int}$.
2. Representación de Campo Dual: Se aplican simultáneamente dos formulaciones distintas de elementos finitos mixtos:
   • Sistema 1 (en $\Omega_1$): Utiliza una formulación donde la condición de contorno en $\Gamma_1$ se impone naturalmente (débilmente), y la condición en la interfaz se trata como una entrada. Esta formulación utiliza espacios de elementos finitos derivados de un subcomplejo de Hilbert donde el espacio de la variable dual es un subconjunto de la imagen del espacio primal bajo el operador adjunto.
   • Sistema 2 (en $\Omega_2$): Utiliza la formulación dual donde la condición de contorno en $\Gamma_2$ se impone naturalmente, y la condición en la interfaz se trata como una entrada. Los espacios de elementos finitos se seleccionan del subcomplejo de Hilbert adjunto.
3. Interconexión de Retroalimentación: Los dos subdominios se acoplan en la interfaz $\Gamma_{int}$ mediante una interconexión de girador (un tipo específico de retroalimentación en sistemas port-hamiltonianos). Esta interconexión fuerza la continuidad de los espacios de elementos finitos de manera débil. Específicamente, la entrada de un subdominio se establece igual a la salida (variable conjugada) del otro, y viceversa.
4. Discretización: El método utiliza subcomplejos de elementos finitos (por ejemplo, complejos de de Rham) para asegurar la estabilidad y la preservación de la estructura. Por ejemplo, en la ecuación de onda, $\Omega_1$ podría utilizar elementos de Galerkin Descontinuos (DG) y Raviart-Thomas (RT), mientras que $\Omega_2$ utiliza elementos de Galerkin Continuos (CG) y Nédélec.
5. Integración Temporal: El artículo evalúa dos esquemas de integración temporal:
   • Punto Medio Implícito: Preserva la energía global y la estructura de Poisson del sistema, pero requiere resolver un sistema monolítico.
   • Störmer-Verlet: Permite el desacoplamiento de los dos subdominios (resolviéndolos por separado), pero no preserva estrictamente la energía total, aunque mantiene el balance de potencia dentro de cada subdominio.

Contribuciones Clave
El artículo formaliza y demuestra las siguientes contribuciones:

• Condiciones de Contorno Mixtas Libres de Lagrange: La estrategia impone con éxito condiciones de contorno mixtas de forma explícita sin introducir multiplicadores de Lagrange, convirtiendo el problema en un sistema de ODE en lugar de DAE.
• Descomposición de Dominios Primal-Dual: A diferencia de los métodos clásicos de Schwarz alternante Dirichlet-Neumann que dependen de un acoplamiento iterativo, este método utiliza una formulación primal-dual simultánea para acoplar los subdominios mediante una interconexión de retroalimentación directa.
• Preservación de la Estructura: El uso de subcomplejos de Hilbert asegura que la discretización preserve la estructura geométrica subyacente (estructura Stokes-Dirac) y las propiedades físicas (por ejemplo, el balance de potencia) del sistema port-hamiltoniano.
• Robustez ante el Bloqueo por Cizalladura (Shear Locking): Se demuestra que la formulación mixta evita los fenómenos de bloqueo por cizalladura en problemas de mecánica estructural, un problema común en las discretizaciones estándar de elementos finitos de placas y vigas.
• Extensibilidad: Se demuestra que la metodología es aplicable a sistemas semilineales que contienen no linealidades algebraicas (donde la no linealidad no interfiere con los operadores diferenciales) y es teóricamente extensible a problemas elípticos estáticos y fenómenos electromagnéticos.

Resultados Numéricos
Los autores validan el método mediante cuatro ejemplos numéricos:

1. Viga Intrínseca Geométricamente Exacta (1D): Un modelo semilineal que demuestra la capacidad del método para manejar no linealidades algebraicas y evitar el bloqueo por cizalladura en un escenario de enrollamiento cuasi-estático.
2. Ecuación de Onda 2D: Un estudio de convergencia confirma las tasas de convergencia teórica ($h^k$). El análisis espectral muestra que el método propuesto produce frecuencias propias con una precisión comparable o superior a las discretizaciones clásicas de elementos Lagrange lineales. El balance de potencia se preserva hasta la precisión de la máquina dentro de cada subdominio.
3. Elastodinámica Lineal: El análisis modal en un cuadrado unitario demuestra que la presente discretización mixta logra una precisión comparable a la de los elementos Lagrange cuadráticos clásicos.
4. Placa de Mindlin: Un problema de referencia utilizado para probar el bloqueo por cizalladura. Los resultados muestran que la presente discretización mixta supera significativamente a las discretizaciones primales clásicas en términos de precisión para las frecuencias propias, eliminando eficazmente el bloqueo por cizalladura incluso con mallas gruesas.

Significancia y Reivindicaciones
El artículo sostiene que esta estrategia ofrece una alternativa robusta para modelar sistemas multifísicos complejos donde la descripción modular y la preservación de la estructura son críticas. Al evitar los multiplicadores de Lagrange, el método simplifica la estructura algebraica de las ecuaciones resultantes, lo cual es beneficioso no solo para la simulación, sino también para la optimización numérica donde las restricciones algebraicas pueden ser problemáticas. Los autores posicionan este trabajo como una "prueba de concepto" que demuestra que las formulaciones mixtas primal-dual pueden utilizarse simultáneamente para manejar condiciones de contorno mixtas de forma natural. Sugieren que, si bien la implementación actual requiere definir una interfaz, el enfoque podría extenderse a métodos de Galerkin Descontinuos totales donde la interconexión ocurra entre cada elemento finito, reduciendo potencialmente el costo computacional de resolver grandes sistemas monolíticos. El trabajo está diseñado primordialmente para sistemas port-hamiltonianos hiperbólicos, pero señala explícitamente su potencial extensión a problemas estáticos y semilineales.