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

Resumen Técnico: Síntesis de Dispersión en el Dominio del Tiempo Eficiente mediante la Sustracción de Singularidades en el Dominio de la Frecuencia

Planteamiento del Problema
El artículo aborda el desafío de simular con precisión la dispersión de ondas en el dominio del tiempo por obstáculos que presentan geometrías de "atrapamiento" (trapping). Aunque los métodos de "Híbrido Frecuencia-Tiempo" (FTH) basados en la transformada de Fourier ofrecen soluciones precisas y sin dispersión para problemas de dispersión generales al resolver las ecuaciones de Helmholtz en frecuencias discretas y realizar una transformada de Fourier inversa, sufren una degradación severa del rendimiento ante obstáculos con propiedades de atrapamiento. En tales casos, la solución de Helmholtz en el dominio de la frecuencia exhibe resonancias (polos) complejas casi reales. Estos polos se manifiestan como características de tipo pico muy agudas a lo largo del eje de frecuencia real, lo que provoca que el integrando de la transformada de Fourier inversa sea casi singular. En consecuencia, la cuadratura numérica estándar requiere un muestreo de frecuencia prohibitivamente denso para lograr la convergencia, lo que hace que las simulaciones de tiempos largos sean computacionalmente inviables. Además, el lento decaimiento temporal de los campos dispersados en regiones de atrapamiento plantea desafíos para los métodos tradicionales en el dominio del tiempo, como los de diferencias finitas o elementos finitos.

Metodología
Los autores proponen un método FTH con Sustracción de Singularidades (FTH-SS). La estrategia central consiste en regularizar el integrando de la transformada de Fourier inversa mediante la sustracción de las contribuciones de las resonancias complejas en el dominio de la frecuencia y su posterior reincorporación mediante una expansión asintótica eficiente para tiempos largos. La metodología consta de tres componentes principales:

1. Algoritmo de Incidencia-Excitación (IE):
   Para evitar el alto costo de evaluar los resolventes sobre contornos de frecuencia complejos, el artículo introduce un nuevo algoritmo adaptativo que identifica las resonancias complejas y sus residuos utilizando únicamente datos de frecuencia real.

   • El método emplea el algoritmo de aproximación racional AAA (Adaptive Antoulas-Anderson) aplicado a las densidades de la ecuación integral excitadas por un campo incidente específico.
   • Utilizando un enfoque de "sketching aleatorio" para funciones de valor vectorial, el algoritmo construye aproximantes racionales a las soluciones de densidad.
   • Una estrategia recursiva adaptativa particiona el intervalo de frecuencia para asegurar que todos los polos relevantes (aquellos con residuos lo suficientemente significativos como para causar picos en la solución) sean capturados.
   • Crucialmente, los aproximantes racionales generados durante la identificación de polos se reutilizan para calcular los residuos y los valores de campo sin requerir inversiones adicionales y costosas de la integral de frontera.

2. Sustracción de Singularidades:
   Una vez identificados los resonantes complejos relevantes ($\rho_n$) y sus residuos ($d_{p,n}$), la solución en el dominio de la frecuencia $U_p(r, \omega)$ se descompone en:
   $$U_p(r, \omega) = U^s_{p,h}(r, \omega) + \sum_{n} \frac{d_{p,n}(r)}{\omega - \rho_n}$$
   donde $U^s_{p,h}$ es el término sustracto de la singularidad (regularizado). La transformada de Fourier inversa se divide en dos partes:

   • $I_1(r, t)$: La integral del término regularizado. Dado que los picos agudos han sido eliminados, este término puede evaluarse eficientemente utilizando reglas de cuadratura FTH de alta frecuencia, incluso para tiempos largos.
   • $I_2(r, t)$: La integral de los términos singulares sustraídos.

3. Evaluación Asintótica de las Contribuciones Singulares:
   Para la parte singular $I_2$, el artículo utiliza una expansión asintótica de tiempo largo derivada mediante la deformación de contorno.

   • En tiempos pre-asintóticos, $I_2$ se evalúa numéricamente utilizando una cuadratura especializada que maneja las casi-singularidades.
   • Para tiempos largos ($t > T_{inc}$), la integral se aproxima mediante una "expansión de singularidad":
     $$E^I_h(r, t) = -i \sum_{n} d_{p,n}(r) A(\sigma_n) e^{-i\sigma_n t}$$
   • El análisis teórico (Teorema 1) proporciona límites de error que muestran que la diferencia entre la integral exacta y la expansión asintótica decae exponencialmente con el tiempo, siempre que el intervalo de frecuencia sea suficientemente grande.

Contribuciones Clave

• Nuevo Algoritmo IE: Un método computacionalmente eficiente para determinar las resonancias excitadas por incidencia y sus residuos utilizando únicamente inversiones de integral de frontera en frecuencia real, evitando la necesidad de integración de contorno en frecuencia compleja o búsqueda de raíces en el plano complejo.
• Regularización vía Sustracción: Una técnica que transforma una transformada de Fourier inversa casi singular en una integral suave y una suma de términos exponenciales conocidos, permitiendo simulaciones de tiempo largo eficientes para estructuras altamente atrapantes.
• Reutilización de Aproximantes Racionales: La demostración de que los aproximantes racionales utilizados para la búsqueda de polos pueden utilizarse directamente para calcular residuos y valores de campo, reduciendo significativamente el costo computacional del algoritmo global.
• Validación de la Expansión de Singularidad: Aunque la validez de las expansiones de singularidad para obstáculos de atrapamiento no ha sido probada rigurosamente en la literatura, el artículo proporciona una fuerte evidencia numérica de que la expansión describe con precisión el comportamiento de tiempo largo de los campos dispersados para geometrías de atrapamiento arbitrarias.

Resultados
Los experimentos numéricos presentados en el artículo demuestran:

• Precisión y Eficiencia: El método FTH-SS logra una alta precisión para la dispersión por obstáculos altamente atrapantes (por ejemplo, círculos de pequeña apertura, cavidades con forma de cohete, estructuras de galería susurrante) donde los métodos FTH estándar fallan o requieren recursos computacionales excesivos.
• Regularización: La sustracción de resonancias suaviza eficazmente el integrando en el dominio de la frecuencia, permitiendo una convergencia rápida de la cuadratura numérica.
• Captura Adaptativa de Polos: El algoritmo IE identifica con éxito todas las resonancias relevantes para un campo incidente específico, filtrando aquellas con residuos insignificantes.
• Comportamiento de Tiempo Largo: El método captura con precisión el lento decaimiento de los campos dispersados en regiones de atrapamiento. Se muestra que la expansión asintótica converge rápidamente a la solución real a medida que el tiempo aumenta, incluso para geometrías complejas.
• Comparación: En casos no atrapantes, el FTH-SS coincide con el FTH estándar. En casos de atrapamiento, el FTH-SS supera significativamente al FTH estándar, logrando una precisión de orden de máquina con órdenes de magnitud menos evaluaciones de frecuencia.

Significancia
El artículo afirma que el enfoque FTH-SS propuesto permite simulaciones de dispersión en el dominio del tiempo precisas y eficientes para tiempos arbitrariamente largos, incluso para estructuras altamente atrapantes. Al desacoplar la necesidad de un muestreo de frecuencia denso de la presencia de resonancias complejas, el método supera un cuello de botella importante en las simulaciones de dispersión de ondas. Además, los resultados numéricos proporcionan evidencia convincente que respalda la validez práctica de las expansiones de singularidad para obstáculos de atrapamiento, un régimen donde la justificación teórica ha sido previamente esquiva. El trabajo establece un marco robusto para manejar los polos "casi reales" que típicamente afectan a los enfoques de dominio de la frecuencia para problemas del dominio del tiempo.