Discrete FEM-BEM coupling with the Generalized Optimized Schwarz Method

Resumo Técnico: Acoplamento Discreto FEM-BEM com o Método de Schwarz Otimizado Generalizado

Enunciado do Problema
O artigo aborda a solução numérica de problemas de propagação de ondas acústicas modelados pela equação de Helmholtz tanto em domínios limitados quanto ilimitados. O objetivo é desenvolver um solver versátil e eficiente capaz de lidar com materiais heterogêneos e geometrias complexas, incluindo aquelas com pontos de cruzamento (cross-points — pontos compartilhados por mais de dois subdomínios). Os autores focam no acoplamento do Método dos Elementos Finitos (FEM), adequado para regiões limitadas heterogêneas, com o Método dos Elementos de Contorno (BEM), ideal para domínios de homogeneidade parcial e ilimitados. Um desafio específico abordado é a potencial mal-definição das formulações clássicas FEM-BEM em "ressonâncias espúrias" (números de onda onde a formulação variacional admite um núcleo não trivial, apesar de o problema físico ser bem definido).

Metodologia
Os autores propõem uma abordagem de Decomposição de Domínio (DD) não sobreposta baseada no Método de Schwarz Otimizado Generalizado (GOSM). A metodologia envolve:

1. Subestruturação: O domínio é particionado em subdomínios não sobrepostos. As incógnitas do solver iterativo estão associadas às interfaces (o esqueleto) em vez do volume.
2. Estrutura Discreta: A análise estende a teoria contínua apresentada em [20] para um cenário totalmente discreto usando aproximações de Galerkin (FEM e BEM).
3. Formulação de Operador: O método reformula o problema global em uma equação subestruturada da forma $(Id + \Pi S)q = \text{rhs}$, onde:
   • $S$ é um operador de espalhamento que resolve problemas locais com condições de contorno de impedância otimizadas.
   • $\Pi$ é um operador de troca não local (substituindo o operador de troca padrão usado nos Métodos de Schwarz Otimizados clássicos) projetado para lidar com pontos de cruzamento.
   • $q$ representa os traços de impedância no esqueleto.
4. Acoplamentos FEM-BEM: O estudo analisa especificamente três acoplamentos FEM-BEM clássicos: Costabel, Johnson-Nédélec e Bielak-MacCamy. Investiga-se a "suposição de sinal" (relacionada à dissipação física) para os operadores de contorno $A_\Gamma$ associados a esses acoplamentos.
5. Solver Iterativo: Quando a suposição de sinal é satisfeita (especificamente para o acoplamento Costabel), o operador $\Pi S$ é uma contração. Isso permite o uso de um esquema iterativo de Richardson com convergência geométrica garantida. Para acoplamentos onde a suposição de sinal não é satisfeita (Johnson-Nédélec, Bielak-MacCamy), um solver GMRes é proposto.
6. Operadores de Transmissão: Os experimentos numéricos utilizam vários operadores de transmissão, incluindo operadores locais de Després, operadores hipersingulares de Yukawa não locais e operadores de complemento de Schur derivados de um mapa de Dirichlet-para-Neumann positivo.

Principais Contribuições

• Bem-definição Discreta: O artigo estabelece que as formulações subestruturadas dos acoplamentos FEM-BEM clássicos são bem definidas para números de onda diferentes de ressonâncias espúrias clássicas.
• Relação de Dimensão do Núcleo: Uma relação explícita é estabelecida entre as dimensões dos núcleos da formulação variacional inicial, dos problemas locais e da formulação subestruturada (Teorema 7.3).
• Robustez em Ressonâncias Espúrias: Um resultado teórico significativo é a prova de que a formulação subestruturada do acoplamento FEM-BEM de Costabel permanece bem definida mesmo em ressonâncias espúrias. Isso é alcançado porque a dimensão do núcleo do operador subestruturado coincide com a do problema local, permitindo um núcleo trivial no sistema subestruturado sob condições específicas.
• Teoria de Convergência: Os autores provam que, para o acoplamento Costabel (que satisfaz a suposição de sinal), o operador subestruturado assume a forma "identidade + contração". Isso garante a convergência geométrica do solver de Richardson.
• Tratamento de Pontos de Cruzamento: O método mostra-se robusto a pontos de cruzamento sem exigir incógnitas adicionais por ponto de cruzamento ou o acoplamento de problemas locais nesses pontos, uma dificuldade comum em outros métodos de DD.
• Robustez em $h$: Estimativas teóricas e experimentos numéricos demonstram que a convergência geométrica pode ser alcançada de forma robusta em relação ao tamanho da malha ($h$-robusta), desde que sejam utilizados operadores de transmissão não locais (como os tipos Yukawa ou complemento de Schur).

Resultos
Experimentos numéricos foram conduzidos em 2D para problemas homogêneos e heterogêneos, bem como para problemas envolvendo obstáculos impenetráveis e cavidades abertas (introduzindo pontos de cruzamento).

• Contagens de Iteração: O número de iterações necessárias para a convergência foi analisado em relação ao tamanho da malha ($h$) e ao número de onda ($\kappa$).
  • Configurações utilizando operadores de transmissão locais (Després) mostraram um crescimento quase quadrático nas iterações conforme a malha era refinada ($O(h^{1.69} \dots O(h^{1.78})$).
  • Configurações utilizando operadores de transmissão não locais (Yukawa, complemento de Schur) exibiram contagens de iteração quase constantes com relação ao refinamento da malha, confirmando a $h$-robustez.
  • Em relação ao número de onda, as configurações não locais geralmente mostraram crescimento quase constante ou sublinear, enquanto as configurações locais mostraram crescimento sublinear.
• Comparação de Acoplamento: Embora as garantias teóricas para a convergência geométrica tenham sido derivadas primariamente para o acoplamento Costabel, os resultados numéricos indicaram que os acoplamentos Johnson-Nédélec e Bielak-MacCamy produziram contagens de iteração semelhantes nos cenários testados.
• Robustez de Pontos de Cruzamento: No caso de teste da cavidade retangular envolvendo pontos de cruzamento, o método manteve a $h$-robustez, com as contagens de iteração permanecendo estáveis conforme a malha era refinada.
• Seleção de Operador de Transmissão: A configuração "Y-S" (operador de Yukawa para BEM, complemento de Schur para FEM) geralmente apresentou o menor número de iterações. O estudo destaca que combinar o tipo de operador de transmissão com a natureza do subdomínio (esparso vs. denso) é benéfico para a eficiência de implementação.

Significância e Alegações
O artigo afirma fornecer o primeiro algoritmo iterativo $h$-robusto para acoplamento FEM-BEM com limites de convergência garantidos. Ao reformular acoplamentos FEM-BEM clássicos em uma forma subestruturada, os autores demonstram que esses métodos clássicos podem ser resolvidos eficientemente usando técnicas de decomposição de domínio, mesmo na presença de pontos de cruzamento e ressonâncias espúrias. O trabalho une a lacuna entre a teoria de decomposição de domínio contínua e as implementações discretas de FEM-BEM, oferecendo uma estratégia sistemática para acoplar esses métodos que é tanto teoricamente sólida quanto numericamente eficiente para aplicações de computação de alto desempenho. Os autores enfatizam que sua abordagem evita a complexidade de acoplar problemas locais em pontos de cruzamento, o que é uma vantagem distinta sobre outras formulações de múltiplos domínios de FEM-BEM.