Efficient time-domain scattering synthesis via frequency-domain singularity subtraction

Resumo Técnico: Síntese de Espalhamento no Domínio do Tempo Eficiente via Subtração de Singularidades no Domínio da Frequência

Definição do Problema
O artigo aborda o desafio de simular com precisão o espalhamento de ondas no domínio do tempo por obstáculos que exibem geometrias de "aprisionamento" (trapping). Embora os métodos "Híbridos de Frequência-Tempo" (FTH), baseados na transformada de Fourier, ofereçam soluções precisas e livres de dispersão para problemas gerais de espalhamento ao resolver equações de Helmholtz em frequências discretas e realizar uma transformada inversa de Fourier, eles sofrem uma degradação severa de desempenho na presença de obstáculos de aprisionamento. Nesses casos, a solução de Helmholtz no domínio da frequência exibe ressonâncias complexas quase reais (polos). Esses polos manifestam-se como características agudas, do tipo pico, ao longo do eixo de frequência real, fazendo com que o integrando da transformada inversa de Fourier se torne quase singular. Consequentemente, a quadratura numérica padrão exige um amostragem de frequência excessivamente densa para alcançar a convergência, tornando as simulações de longo tempo computacionalmente inviáveis. Além disso, o decaimento temporal lento dos campos espalhados em regiões de aprisionamento impõe desafios para métodos tradicionais no domínio do tempo, como as técnicas de diferenças finitas ou elementos finitos.

Metodologia
Os autores propõem um método FTH com Subtração de Singularidade (FTH-SS). A estratégia central envolve regularizar o integrando da transformada inversa de Fourier subtraindo as contribuições das ressonâncias complexas no domínio da frequência e, em seguida, reincorporando-as via uma expansão assintótica eficiente para tempos longos. A metodologia consiste em três componentes primários:

1. Algoritmo de Excitação-Incidência (IE):
   Para evitar o alto custo de avaliar resolventes sobre contornos de frequência complexos, o artigo introduz um novo algoritmo adaptativo que identifica ressonâncias complexas e seus resíduos utilizando apenas dados de frequência real.

   • O método utiliza o algoritmo de aproximação racional AAA (Adaptive Antoulas-Anderson) aplicado às densidades de equações integrais excitadas por um campo incidente específico.
   • Utilizando uma abordagem de "esboço aleatório" (random sketching) para funções de valores vetoriais, o algoritmo constrói aproximantes racionais para as soluções de densidade.
   • Uma estratégia recursiva adaptativa particiona o intervalo de frequência para garantir que todos os polos relevantes (aqueles com resíduos significativos o suficiente para causar picos na solução) sejam capturados.
   • Crucialmente, os aproximantes racionais gerados durante a identificação de polos são reutilizados para calcular os resíduos e os valores de campo sem exigir inversões adicionais e custosas de integrais de contorno.

2. Subtração de Singularidade:
   Uma vez identificados os ressonâncias complexas relevantes ($\rho_n$) e seus resíduos ($d_{p,n}$), a solução no domínio da frequência $U_p(r, \omega)$ é decomposta:
   $$U_p(r, \omega) = U^s_{p,h}(r, \omega) + \sum_{n} \frac{d_{p,n}(r)}{\omega - \rho_n}$$
   onde $U^s_{p,h}$ é o termo subtraído da singularidade (regularizado). A transformada inversa de Fourier é dividida em duas partes:

   • $I_1(r, t)$: A integral do termo regularizado. Como os picos agudos foram removidos, este termo pode ser avaliado eficientemente usando regras de quadratura FTH de alta frequência padrão, mesmo para tempos longos.
   • $I_2(r, t)$: A integral dos termos singulares subtraídos.

3. Avaliação Assintótica de Contribuições Singulares:
   Para a parte singular $I_2$, o artigo utiliza uma expansão assintótica de grande tempo derivada via deformação de contorno.

   • Em tempos pré-assintóticos, $I_2$ é avaliado numericamente usando uma quadratura especializada que lida com as quase-singularidades.
   • Para tempos longos ($t > T_{inc}$), a integral é aproximada por uma "expansão de singularidade":
     $$E^I_h(r, t) = -i \sum_{n} d_{p,n}(r) A(\sigma_n) e^{-i\sigma_n t}$$
   • A análise teórica (Teorema 1) fornece limites de erro mostrando que a diferença entre a integral exata e a expansão assintótica decai exponencialmente com o tempo, desde que o intervalo de frequência seja suficientemente grande.

Principais Contribuições

• Novo Algoritmo IE: Um método computacionalmente eficiente para determinar ressonâncias complexas excitadas pela incidência e seus resíduos usando apenas inversões de integral de contorno de frequência real, evitando a necessidade de integração de contorno de frequência complexa ou busca de raízes no plano complexo.
• Regularização via Subtração: Uma técnica que transforma uma transformada de Fourier inversa quase singular em uma integral suave e uma soma de termos exponenciais conhecidos, permitindo a simulação eficiente de longo tempo para estruturas altamente de aprisionamento.
• Reuso de Aproximantes Racionais: A demonstração de que os aproximantes racionais usados para a busca de polos podem ser diretamente utilizados para calcular resíduos e valores de campo, reduzindo significativamente o custo computacional de todo o algoritmo.
• Validação da Expansão de Singularidade: Embora a validade das expansões de singularidade para obstáculos de aprisionamento não tenha sido rigorosamente provada na literatura, o artigo fornece evidências numéricas fortes de que a expansão descreve com precisão o comportamento de tempo tardio de campos espalhados para geometrias de aprisionamento arbitrárias.

Resultados
Experimentos numéricos apresentados no artigo demonstram:

• Precisão e Eficiência: O método FTH-SS alcança alta precisão para o espalhamento por obstáculos de alto aprisionamento (ex: círculos de pequena abertura, cavidades em forma de foguete, estruturas de galerias sussurrantes/whispering-gallery) onde os métodos FTH padrão falham ou exigem recursos computacionais excessivos.
• Regularização: A subtração de ressonâncias suaviza efetivamente o integrando no domínio da frequência, permitindo uma convergência rápida da quadratura numérica.
• Captura Adaptativa de Polos: O algoritmo IE identifica com sucesso todas as ressonâncias relevantes para um campo incidente específico, filtrando aquelas com resíduos negligenciáveis.
• Comportamento de Longo Tempo: O método captura com precisão o decaimento lento de campos espalhados em regiões de aprisionamento. A expansão assintótica mostra-se capaz de convergir rapidamente para a solução real conforme o tempo aumenta, mesmo para geometrias complexas.
• Comparação: Em casos não de aprisionamento, o FTH-SS coincide com o FTH padrão. Em casos de aprisionamento, o FTH-SS supera significativamente o FTH padrão, alcançando precisão de ordem de máquina com ordens de magnitude a menos de avaliações de frequência.

Significância
O artigo afirma que a abordagem FTH-SS proposta permite simulações de espalhamento no domínio do tempo precisas e eficientes para tempos arbitrariamente longos, mesmo para estruturas de alto aprisionamento. Ao desacoplar a necessidade de amostragem densa de frequência da presença de ressonâncias complexas, o método supera um grande gargalo nas simulações de espalhamento de ondas. Além disso, os resultados numéricos fornecem evidências convincentes que apoiam a validade prática das expansões de singularidade para obstáculos de aprisionamento, um regime onde a justificativa teórica tem sido anteriormente elusiva. O trabalho estabelece um arcabouço robusto para lidar com os polos "quase reais" que tipicamente assolam as abordagens de domínio da frequência para problemas de domínio do tempo.