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

Resumo Técnico: Uma Estratégia de Decomposição de Domínio para a Imposição Natural de Condições de Contorno Mistas em Sistemas Port-Hamiltonianos

Definição do Problema
O artigo aborda o desafio de impor condições de contorno mistas (condições de Dirichlet e Neumann simultâneas) em sistemas port-hamiltonianos (pH) hiperbólicos utilizando métodos de elementos finitos. Em formulações padrão de elementos finitos, as condições de contorno são tipicamente impostas ou fortemente (por meio da restrição dos espaços de funções discretas) ou fracamente (via termos de contorno naturais na formulação variacional). A imposição forte frequentemente leva a equações diferenciais algébrico-diferenciais (DAEs), que são computacionalmente mais difíceis de resolver do que equações diferenciais ordinárias (ODEs). Por outro lado, a imposição fraca geralmente requer a introdução de multiplicadores de Lagrange para impor as condições de contorno essenciais, aumentando a dimensão e a complexidade do sistema. Os autores observam que, embora o arcabouço port-hamiltoniano trate naturalmente os fluxos de contorno via estrutura Stokes-Dirac, os esquemas de discretização existentes muitas vezes têm dificuldade em incorporar condições mistas explicitamente sem o uso de multiplicadores de Lagrange, sem comprometer a estrutura geométrica subjacente.

Metodologia
A estratégia proposta emprega uma abordagem de decomposição de domínio combinada com Cálculo Exterior de Elementos Finitos (FEEC) e uma formulação primal-dual. A metodologia central envolve as seguintes etapas:

1. Decomposição de Domínio: O domínio espacial $\Omega$ é dividido em dois subdomínios não sobrepostos, $\Omega_1$ e $\Omega_2$, separados por uma interface interna arbitrária $\Gamma_{int}$.
2. Representação de Campo Dual: Duas formulações distintas de elementos finitos mistos são aplicadas simultaneamente:
   • Sistema 1 (em $\Omega_1$): Utiliza uma formulação onde a condição de contorno em $\Gamma_1$ é imposta naturalmente (fracamente), e a condição na interface é tratada como uma entrada. Esta formulação utiliza espaços de elementos finitos derivados de um subcomplexo de Hilbert onde o espaço da variável dual é um subconjunto da imagem do espaço primal sob o operador adjunto.
   • Sistema 2 (em $\Omega_2$): Utiliza a formulação dual onde a condição de contorno em $\Gamma_2$ é imposta naturalmente, e a condição na interface é tratada como uma entrada. Os espaços de elementos finitos são selecionados a partir do subcomplexo de Hilbert adjunto.
3. Interconexão de Feedback: Os dois subdomínios são acoplados na interface $\Gamma_{int}$ via uma interconexão de girador (um tipo específico de feedback em sistemas port-hamiltonianos). Esta interconexão impõe a continuidade dos espaços de elementos finitos de maneira fraca. Especificamente, a entrada de um subdomínio é definida como igual à saída (variável conjugada) do outro, e vice-versa.
4. Discretização: O método utiliza subcomplexos de elementos finitos (por exemplo, complexos de de Rham) para garantir estabilidade e preservação de estrutura. Por exemplo, na equação da onda, $\Omega_1$ pode utilizar elementos de Galerkin Descontínuo (DG) e Raviart-Thomas (RT), enquanto $\Omega_2$ utiliza elementos de Galerkin Contínuo (CG) e Nédélec.
5. Integração Temporal: O artigo avalia dois esquemas de integração temporal:
   • Ponto Médio Implícito: Preserva a energia global e a estrutura de Poisson do sistema, mas requer a resolução de um sistema monolítico.
   • Störmer-Verlet: Permite o desacoplamento dos dois subdomínios (resolvendo-os separadamente), mas não preserva estritamente a energia total, embora mantenha o balanço de potência dentro de cada subdomínio.

Principais Contribuições
O artigo formaliza e demonstra as seguintes contribuições:

• BCs Mistas Livres de Lagrange: A estratégia impõe com sucesso condições de contorno mistas explicitamente sem introduzir multiplicadores de Lagrange, convertendo o problema em um sistema de ODEs em vez de DAEs.
• Decomposição de Domínio Primal-Dual: Diferente dos métodos clássicos de Schwarz alternados Dirichlet-Neumann, que dependem de acoplamento iterativo, este método utiliza uma formulação primal-dual simultânea para acoplar os subdomínios via uma interconexão de feedback direta.
• Preservação de Estrutura: O uso de subcomplexos de Hilbert garante que a discretização preserve a estrutura geométrica subjacente (estrutura Stokes-Dirac) e as propriedades físicas (ex: balanço de potência) do sistema port-hamiltoniano.
• Robustez ao Bloqueio de Cisalhamento (Shear Locking): A formulação mista mostra-se capaz de evitar o fenômeno de shear locking em problemas de mecânica estrutural, um problema comum em discretizações padrão de elementos finitos para placas e vigas.
• Extensibilidade: A metodologia demonstra ser aplicável a sistemas semilineares contendo não linearidades algébricas (onde a não linearidade não interfere nos operadores diferenciais) e é teoricamente extensível a problemas elípticos estáticos e fenômenos eletromagnéticos.

Resultados Numéricos
Os autores validam o método através de quatro exemplos numéricos:

1. Viga Intrínseca Geometricamente Exata (1D): Um modelo semilinear demonstrando a capacidade do método de lidar com não linearidades algébricas e evitar o shear locking em um cenário de enrolamento quase estático.
2. Equação da Onda 2D: Um estudo de convergência confirma as taxas de convergência teórica ($h^k$). A análise espectral mostra que o método proposto produz frequências próprias com precisão comparável ou superior às discretizações clássicas de elementos Lagrange lineares. O balanço de potência é preservado até a precisão de máquina dentro de cada subdomínio.
3. Elastodinâmica Linear: A análise modal em um quadrado unitário demonstra que a proposta de discretização mista alcança precisão comparável aos elementos Lagrange quadráticos clássicos.
4. Placa de Mindlin: Um problema de referência utilizado para testar o shear locking. Os resultados mostram que a proposta de discretização mista supera significativamente as discretizações primais clássicas em termos de precisão para frequências próprias, eliminando efetivamente o shear locking mesmo com malhas grosseiras.

Significância e Alegações
O artigo afirma que esta estratégia oferece uma alternativa robusta para modelar sistemas multifísicos complexos onde a descrição modular e a preservação de estrutura são críticas. Ao evitar multiplicadores de Lagrange, o método simplifica a estrutura algébrica das equações resultantes, o que é benéfico não apenas para a simulação, mas também para a otimização numérica, onde restrições algébricas podem ser problemáticas. Os autores posicionam o trabalho como uma "prova de conceito" demonstrando que formulações mistas primal-duais podem ser usadas simultaneamente para lidar com condições de contorno mistas de forma natural. Eles sugerem que, embora a implementação atual exija a definição de uma interface, a abordagem poderia ser estendida para métodos de Galerkin descontínuos completos, onde a interconexão ocorreria entre cada elemento finito, potencialmente reduzindo o custo computacional de resolver grandes sistemas monolíticos. O trabalho é projetado primariamente para sistemas port-hamiltonianos hiperbólicos, mas nota explicitamente seu potencial de extensão para problemas estáticos e semilineares.