Discrete FEM-BEM coupling with the Generalized Optimized Schwarz Method

Resumen Técnico: Acoplamiento Discreto FEM-BEM con el Método de Schwarz Optimizado Generalizado

Planteamiento del Problema
El artículo aborda la solución numérica de problemas de propagación de ondas acústicas modelados por la ecuación de Helmholtz tanto en dominios acotados como no acotados. El objetivo es desarrollar un resolvedor versátil y eficiente capaz de manejar materiales heterogéneos y geometrías complejas, incluyendo aquellas con puntos de cruce (puntos compartidos por más de dos subdominios). Los autores se centran en acoplar el Método de Elementos Finitos (FEM), adecuado para regiones acotadas heterogéneas, con el Método de Elementos de Contorno (BEM), óptimo para dominios homogéneos por partes y no acotados. Un desafío específico abordado es la potencial maldad de la formulación clásica de FEM-BEM en "resonancias espurias" (números de onda donde la formulación variacional admite un núcleo no trivial a pesar de que el problema físico sea bien planteado).

Metodología
Los autores proponen un enfoque de Descomposición de Dominio (DD) no superpuesto basado en el Método de Schwarz Optimizado Generalizado (GOSM). La metodología implica:

1. Subestructuración: El dominio se particiona en subdominios no superpuestos. Las incógnitas del resolvedor iterativo están asociadas a las interfaces (el esqueleto) en lugar del volumen.
2. Marco Discreto: El análisis extiende la teoría continua presentada en [20] a un entorno totalmente discreto utilizando aproximaciones de Galerkin (FEM y BEM).
3. Formulación de Operadores: El método reformula el problema global en una ecuación subestructurada de la forma $(Id + \Pi S)q = \text{rhs}$, donde:
   • $S$ es un operador de dispersión que resuelve problemas locales con condiciones de contorno de impedancia optimizadas.
   • $\Pi$ es un operador de intercambio no local (que reemplaza al operador de intercambio estándar utilizado en los Métodos de Schwarz Optimizados clásicos) diseñado para manejar puntos de cruce.
   • $q$ representa las trazas de impedancia en el esqueleto.
4. Acoplamientos FEM-BEM: El estudio analiza específicamente tres acoplamientos FEM-BEM clásicos: Costabel, Johnson-Nédélec y Bielak-MacCamy. Investiga la "suposición de signo" (relacionada con la disipación física) para los operadores de frontera $A_\Gamma$ asociados con estos acoplamientos.
5. Resolvedor Iterativo: Cuando se cumple la suposición de signo (específicamente para el acoplamiento Costabel), el operador $\Pi S$ es una contracción. Esto permite el uso de un esquema iterativo de Richardson con convergencia geométrica garantizada. Para los acoplamientos donde la suposición de signo no se cumple (Johnson-Nédélec, Bielak-MacCamy), se propone un resolvedor GMRes.
6. Operadores de Transmisión: Los experimentos numéricos utilizan varios operadores de transmisión, incluyendo operadores locales de Després, operadores hipersingulares de Yukawa no locales y operadores de complemento de Schur derivados de un mapa de Dirichlet-a-Neumann positivo.

Contribuciones Clave

• Bien Planteamiento Discreto: El artículo establece que las formulaciones subestructuradas de los acoplamientos FEM-BEM clásicos son bien planteadas para números de onda diferentes de las resonancias espurias clásicas.
• Relación de Dimensión del Núcleo: Se establece una relación explícita entre las dimensiones de los núcleos de la formulación variacional inicial, los problemas locales y la formulación subestructurada (Teorema 7.3).
• Robustez ante Resonancias Espurias: Un resultado teórico significativo es la prueba de que la formulación subestructurada del acoplamiento FEM-BEM de Costabel permanece bien planteada incluso en resonancias espurias. Esto se logra porque la dimensión del núcleo del operador subestructurado coincide con la del problema local, lo que permite un núcleo trivial en el sistema subestructurado bajo condiciones específicas.
• Teoría de Convergencia: Los autores demuestran que para el acoplamiento de Costabel (que satisface la suposición de signo), el operador subestructurado toma la forma de "identidad + contracción". Esto garantiza la convergencia geométrica del resolvedor de Richardson.
• Manejo de Puntos de Cruce: Se demuestra que el método es robusto ante puntos de cruce sin requerir incógnitas adicionales por punto de cruce o el acoplamiento de problemas locales en estos puntos, una dificultad común en otros métodos de DD.
• Robustez h: Los estimadores teóricos y los experimentos numéricos demuestran que se puede lograr una convergencia geométrica que es robusta con respecto al tamaño de la malla ($h$), siempre que se utilicen operadores de transmisión no locales (como los tipos Yukawa o complemento de Schur).

Resultados
Se realizaron experimentos numéricos en 2D para problemas homogéneos y heterogéneos, así como para problemas que involucran obstáculos impenetrables y cavidades abiertas (introduciendo puntos de cruce).

• Conteos de Iteración: Se analizó el número de iteraciones requeridas para la convergencia frente al tamaño de la malla ($h$) y el número de onda ($\kappa$).
  • Las configuraciones que utilizan operadores de transmisión locales (Després) mostraron un crecimiento cuasi-cuadrático en las iteraciones a medida que la malla se refinaba ($O(h^{1.69} \dots O(h^{1.78})$).
  • Las configuraciones que utilizan operadores de transmisión no locales (Yukawa, complemento de Schur) exhibieron conteos de iteración cuasi-constantes con respecto al refinamiento de la malla, confirmando la robustez h.
  • Con respecto al número de onda, las configuraciones no locales generalmente mostraron un crecimiento cuasi-constante o sublineal, mientras que las configuraciones locales mostraron un crecimiento sublineal.
• Comparación de Acoplamientos: Aunque las garantías teóricas para la convergencia geométrica se derivaron principalmente para el acoplamiento de Costabel, los resultados numéricos indicaron que los acoplamientos Johnson-Nédélec y Bielak-MacCamy produjeron conteos de iteración similares en los escenarios probados.
• Robustez de Puntos de Cruce: En el caso de prueba de la cavidad rectangular que involucra puntos de cruce, el método mantuvo la robustez h, con conteos de iteración estables a medida que la malla se refinaba.
• Selección de Operador de Transmisión: La configuración "Y-S" (operador Yukawa para BEM, complemento de Schur para FEM) generalmente produjo los conteos de iteración más bajos. El estudio destaca que emparejar el tipo de operador de transmisión con la naturaleza del subdominio (disperso vs. denso) es beneficioso para la eficiencia de la implementación.

Significancia y Reivindicaciones
El artículo afirma proporcionar el primer algoritmo iterativo h-robusto para el acoplamiento FEM-BEM con límites de convergencia garantizados. Al reformular los acoplamientos FEM-BEM clásicos en una forma subestructurada, los autores demuestran que estos métodos clásicos pueden resolverse eficientemente utilizando técnicas de descomposición de dominio incluso en presencia de puntos de cruce y resonancias espurias. El trabajo cierra la brecha entre la teoría de descomposición de dominios continua y las implementaciones discretas de FEM-BEM, ofreciendo una estrategia sistemática para acoplar estos métodos que es tanto teóricamente sólida como numéricamente eficiente para aplicaciones de alto rendimiento. Los autores enfatizan que su enfoque evita la complejidad de acoplar problemas locales en los puntos de cruce, lo cual es una ventaja distintiva sobre otras formulaciones FEM-BEM de multidominio.