Fast high-order solvers for the Lippmann--Schwinger equation in piecewise-smooth heterogeneous media

Resumen Técnico: Solucionadores de Alto Orden Rápidos para la Ecuación de Lippmann–Schwinger en Medios Heterogéneos de Suavidad Segmentada

Planteamiento del Problema
El artículo aborda la solución numérica de la ecuación de Lippmann–Schwinger bidimensional que surge de la dispersión de ondas armónicas en el tiempo por medios heterogéneos de suavidad segmentada (piecewise-smooth) y penetrables. El desafío principal radica en alcanzar simultáneamente dos objetivos que suelen ser mutuamente excluyentes en la literatura existente:

1. Precisión de alto orden: Resolver con precisión las soluciones cuando el índice de refracción (contraste) es solo de suavidad segmentada, exhibiendo discontinuidades de salto a través de las interfaces de los materiales.
2. Convergencia iterativa robusta frente a la frecuencia: Resolver los sistemas lineales densos resultantes con un número de iteraciones de subespacio de Krylov (GMRES) que permanezca acotado a medida que aumenta el número de onda.

Los métodos de alto orden existentes suelen depender de mallas cartesianas uniformes, que pierden precisión cuando las interfaces no coinciden con la malla. Por el contrario, los métodos que se ajustan a la interfaz (como los esquemas de Nyström de malla ajustada) ofrecen una precisión de alto orden, pero a menudo sufren de conteos de iteración que crecen rápidamente a altas frecuencias debido a que carecen de precondicionadores compatibles y robustos frente a la frecuencia.

Metodología
Los autores proponen un solucionador que acopla una discretización de alto orden que se ajusta a la interfaz con un precondicionador de malla estructurada mediante operadores de transferencia.

1. Discretización No Estructurada (El Solucionador):

   • El dominio $\Omega$ (soporte del contraste) se discretiza utilizando una malla triangular ajustada que se conforma a las interfaces de los materiales.
   • Se emplea un método de Nyström de alto orden utilizando nodos de cuadratura de Vioreanu–Rokhlin (VR).
   • El potencial de volumen se evalúa utilizando el Método de Interpolación de Densidad de Volumen (VDIM). El VDIM regulariza el núcleo débilmente singular construyendo interpolantes polinómicos de la densidad en cada triángulo, lo que permite una integración de alta precisión incluso a través de las discontinuidades.
   • El sistema lineal resultante es denso, pero permite la aplicación libre de matrices (matrix-free) del operador de potencial de volumen utilizando el Método de Multipolos Rápidos (FMM) o compresión H-matrix, lo que produce una complejidad cuasi-lineal por iteración.

2. Estrategia de Precondicionamiento:

   • Para acelerar la convergencia, el sistema no estructurado se precondiciona utilizando un "precondicionador de esparcimiento" (específicamente el método de Ying [46]) diseñado para mallas cartesianas uniformes.
   • Operadores de Transferencia: Dos operadores conectan los nodos de cuadratura no estructurados y la malla cartesiana uniforme:
     • De Cartesiana a Cuadratura ($T^Q_C$): Una prolongación bilineal que mapea los valores de la malla a los nodos de cuadratura.
     • De Cuadratura a Cartesiana ($T^C_Q$): El artículo propone y analiza dos construcciones para este operador de restricción:
       • Enfoque I: Una proyección $L^2$ local sobre funciones lineales por partes.
       • Enfoque II: Una transpuesta renormalizada del operador de prolongación (estándar en multigrid geométrico).
   • El precondicionador $P$ se construye como $P = T^Q_C S^{(C)} T^C_Q + (I - T^Q_C T^C_Q)$, donde $S^{(C)}$ es el precondicionador de esparcimiento en la malla cartesiana. El término de corrección $(I - T^Q_C T^C_Q)$ es esencial para compensar la deficiencia de rango del término de incrustación.

3. Análisis Teórico:

   • El artículo establece condiciones suficientes para la invertibilidad del precondicionador $P$.
   • Para el Enfoque I, la invertibilidad se demuestra bajo la condición de que la razón de malla $r_h = h/h_C$ (razón entre el tamaño de la malla no estructurada y el paso cartesiano) sea suficientemente pequeña.
   • Para el Enfoque II, la invertibilidad se demuestra para cualquier tamaño de malla $h$, siempre que el precondicionador de esparcimiento $S^{(C)}$ satisfaga una condición espectral (evitando el eje real no positivo en un producto interno específico), la cual está vinculada al balance de energía del problema de dispersión continuo.

Contribuciones Clave

1. Solucionador Novedoso: El primer solucionador rápido de alto orden de Nyström para la ecuación de Lippmann–Schwinger 2D en medios de suavidad segmentada que mantiene una alta precisión en las interfaces mientras utiliza solucionadores rápidos basados en convolución.
2. Precondicionamiento Híbrido: Una nueva estrategia que acopla discretizaciones de malla ajustada con precondicionadores de malla estructurada mediante operadores de transferencia, permitiendo el uso de solucionadores robustos frente a la frecuencia (como el precondicionador de esparcimiento) en mallas no estructuradas.
3. Pruebas de Invertibilidad: Garantías teóricas rigurosas para la invertibilidad del precondicionador resultante bajo condiciones explícitas de refinamiento de malla y propiedades del solucionador cartesiano subyacente.
4. Validación Numérica: Experimentos extensos que demuestran que el sistema precondicionado requiere significativamente menos iteraciones de GMRES que sus contrapartes sin precondicionar. Los conteos de iteraciones muestran ser casi independientes tanto del tamaño de la malla como del número de onda, incluso para la dispersión de alta frecuencia y los índices de refracción de alto contraste con suavidad segmentada.

Resultados y Significado
El artículo afirma resolver la tensión de larga data entre la flexibilidad geométrica de alto orden y la convergencia iterativa robusta frente a la frecuencia. Los resultados numéricos ilustran que el método propuesto:

• Logra tasas de convergencia de alto orden para contrastes de suavidad segmentada con discontinuidades de salto.
• Exhibe un comportamiento "robusto frente a la frecuencia", donde el número de iteraciones de GMRES permanece acotado a medida que aumenta el número de onda, una propiedad que típicamente se pierde en los métodos no estructurados de alto orden.
• Mantiene la estabilidad con respecto al refinamiento de la malla.

Los autores posicionan este trabajo como un paso significativo hacia la resolución de problemas de dispersión de alta frecuencia en estructuras complejas creadas por el hombre (por ejemplo, nanofotónica y óptica integrada) donde las interfaces de materiales son nítidas y los índices de refracción varían continuamente dentro de los subdominios. El método aprovecha los mejores aspectos tanto de la precisión de malla ajustada como de la eficiencia de los solucionadores de malla estructurada sin sacrificar ninguno de los dos.