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

Resumo Técnico: Solvers de Alta Ordem Rápidos para a Equação de Lippmann–Schwinger em Meios Heterogêneos de Suavidade Segmentada

Enunciado do Problema
O artigo aborda a solução numérica da equação de Lippmann–Schwinger bidimensional decorrente da dispersão de ondas harmônicas temporais por meios heterogêneos penetráveis e de suavidade segmentada (piecewise-smooth). O desafio principal reside em alcançar simultaneamente dois objetivos que são frequentemente mutuamente exclusivos na literatura existente:

1. Precisão de alta ordem: Resolver com precisão as soluções quando o índice de refração (contraste) é apenas de suavidade segmentada, exibindo descontinuidades de salto através de interfaces de materiais.
2. Convergência iterativa robusta à frequência: Resolver os sistemas lineares densos resultantes com um número de iterações de subespaço de Krylov (GMRES) que permanece limitado conforme o número de onda aumenta.

Os métodos de alta ordem existentes tipicamente dependem de grades cartesianas uniformes, que perdem precisão quando as interfaces não são conformes à grade. Por outro lado, métodos conformes à interface (como esquemas de Nyström de malha ajustada) oferecem precisão de alta ordem, mas frequentemente sofrem com contagens de iteração que crescem rapidamente em altas frequências por carecerem de precondicionadores compatíveis e robustos à frequência.

Metodologia
Os autores propõem um solver que acopla uma discretização de alta ordem conforme à interface com um precondicionador de grade estruturada via operadores de transferência.

1. Discretização Não Estruturada (O Solver):

   • O domínio $\Omega$ (suporte do contraste) é discretizado usando uma malha triangular ajustada que é conforme às interfaces de materiais.
   • Um método de Nyström de alta ordem é empregado utilizando nós de quadratura de Vioreanu–Rokhlin (VR).
   • O potencial de volume é avaliado usando o Método de Interpolação de Densidade de Volume (VDIM). O VDIM regulariza o núcleo fracamente singular construindo interpolantes polinomiais da densidade em cada triângulo, permitindo integração de alta precisão mesmo através de descontinuidades.
   • O sistema linear resultante é denso, mas permite a aplicação matrix-free do operador de potencial de volume usando o Método de Multipolo Rápido (FMM) ou compressão H-matrix, resultando em complexidade quase linear por iteração.

2. Estratégia de Precondicionamento:

   • Para acelerar a convergência, o sistema não estruturado é precondicionado usando um "precondicionador de esparsificação" (especificamente o método de Ying [46]) projetado para grades cartesianas uniformes.
   • Operadores de Transferência: Dois operadores fazem a ponte entre os nós de quadratura não estruturados e a grade cartesiana uniforme:
     • De Cartesiano para Quadratura ($T^Q_C$): Uma prolongação bilinear que mapeia valores da grade para os nós de quadratura.
     • De Quadratura para Cartesiano ($T^C_Q$): O artigo propõe e analisa duas construções para este operador de restrição:
       • Abordagem I: Uma projeção $L^2$ local sobre funções linearmente segmentadas.
       • Abordagem II: Uma transposta renormalizada do operador de prolongação (padrão em multigrid geométrico).
   • O precondicionador $P$ é construído como $P = T^Q_C S^{(C)} T^C_Q + (I - T^Q_C T^C_Q)$, onde $S^{(C)}$ é o precondicionador de esparsificação na grade cartesiana. O termo de correção $(I - T^Q_C T^C_Q)$ é essencial para compensar a deficiência de posto do termo de imersão.

3. Análise Teórica:

   • Os autores estabelecem condições suficientes para a invertibilidade do precondicionador $P$.
   • Para a Abordagem I, a invertibilidade é provada sob a condição de que a razão de malha $r_h = h/h_C$ (razão entre o tamanho da malha não estruturada e o passo cartesiano) seja suficientemente pequena.
   • Para a Abordagem II, a invertibilidade é provada para qualquer tamanho de malha $h$, desde que o precondicionador de esparsificação $S^{(C)}$ satisfaça uma condição espectral (evitando o eixo real não positivo em um produto interno específico), que está ligada ao balanço de energia do problema de dispersão contínuo.

Principais Contribuições

1. Novo Solver: O primeiro solver de Nyström rápido e de alta ordem para a equação de Lippmann–Schwinger 2D em meios de suavidade segmentada que mantém alta precisão nas interfaces enquanto utiliza solvers baseados em convolução rápida.
2. Precondicionamento Híbrido: Uma nova estratégia que acopla discretizações de malha ajustada com precondicionadores de grade estruturada via operadores de transferência, permitindo o uso de solvers robustos à frequência (como o precondicionador de esparsificação) em malhas não estruturadas.
3. Provas de Invertibilidade: Garantias teóricas rigorosas para a invertibilidade do precondicionador resultante sob condições explícitas de refinamento de malha e propriedades do solver cartesiano subjacente.
4. Validação Numérica: Experimentos extensivos demonstrando que o sistema precondicionado requer significativamente menos iterações de GMRES do que os equivalentes não precondicionados. As contagens de iteração mostram-se quase independentes tanto do tamanho da malha quanto do número de onda, mesmo para dispersão de alta frequência e índices de refração de alto contraste com suavidade segmentada.

Resultos e Significância
O artigo afirma resolver a tensão de longa data entre a flexibilidade geométrica de alta ordem e a convergência iterativa robusta à frequência. Os resultados numéricos ilustram que o método proposto:

• Alcança taxas de convergência de alta ordem para contrastes de suavidade segmentada com descontinuidades de salto.
• Exibe comportamento "robusto à frequência", onde o número de iterações de GMRES permanece limitado conforme o número de onda aumenta, uma propriedade tipicamente perdida em métodos de alta ordem não estruturados.
• Mantém estabilidade em relação ao refinamento da malha.

Os autores posicionam este trabalho como um passo significativo para a resolução de problemas de dispersão de alta frequência em estruturas complexas criadas pelo homem (ex: nanofotônica e óptica integrada), onde as interfaces de materiais são nítidas e os índices de refração variam continuamente dentro dos subdomínios. O método aproveita os melhores aspectos da precisão de malha ajustada e da eficiência de solvers de grade estruturada sem sacrificar nenhum dos dois.