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

Résumé Technique : Synthèse de la Diffusion Temporelle Efficace via la Soustraction de Singularités dans le Domaine Fréquentiel

Énoncé du Problème
L'article traite du défi que représente la simulation précise de la diffusion d'ondes temporelles par des obstacles présentant des géométries de « piégeage » (trapping). Bien que les méthodes « Hybrides Fréquence-Temps » (FTH) basées sur la transformée de Fourier offrent des solutions précises et sans dispersion pour les problèmes de diffusion généraux en résolvant les équations de Helmholtz à des fréquences discrètes et en effectuant une transformée de Fourier inverse, elles souffrent d'une dégradation sévère des performances en présence d'obstacles piégeants. Dans ces cas, la solution de Helmholtz dans le domaine fréquentiel présente des résonances complexes quasi-réelles (pôles). Ces pôles se manifestent sous la forme de pics acérés le long de l'axe des fréquences réelles, rendant l'intégrande de la transformée de Fourier inverse presque singulière. Par conséquent, la quadrature numérique standard nécessite un échantillonnage fréquentiel extrêmement dense pour atteindre la convergence, ce qui rend les simulations de longue durée temporellement infaisables. De plus, la décroissance lente des champs diffusés dans les régions de piégeage pose des défis aux méthodes traditionnelles du domaine temporel, telles que les méthodes de différences finies ou d'éléments finis.

Méthodologie
Les auteurs proposent une méthode FTH avec Soustraction de Singularités (FTH-SS). La stratégie centrale consiste à régulariser l'intégrande de la transformée de Fourier inverse en soustrayant les contributions des résonances complexes dans le domaine fréquentiel, puis en les réincorporant via une expansion asymptotique efficace pour les temps longs. La méthodologie se compose de trois composantes principales :

1. Algorithme d'Excitation par Incidence (IE) :
   Pour éviter le coût élevé de l'évaluation des résolvants sur des contours de fréquences complexes, l'article introduit un nouvel algorithme adaptatif qui identifie les résonances complexes et leurs résidus en utilisant uniquement des données de fréquences réelles.

   • La méthode utilise l'algorithme d'approximation rationnelle AAA (Adaptive Antoulas-Anderson) appliquée aux densités d'équations intégrales excitées par un champ incident spécifique.
   • En utilisant une approche de « sketching aléatoire » pour les fonctions à valeurs vectorielles, l'algorithme construit des approximants rationnels des solutions de densité.
   • Une stratégie récursive adaptative partitionne l'intervalle de fréquence pour s'assurer que tous les pôles pertinents (ceux dont les résidus sont suffisamment importants pour causer des pics dans la solution) sont capturés.
   • Crucialement, les approximants rationnels générés lors de l'identification des pôles sont réutilisés pour calculer les résidus et les valeurs de champ sans nécessiter d'inversions supplémentaires et coûteuses des intégrales de bord.

2. Soustraction de Singularités :
   Une fois les résonances complexes pertinentes ($\rho_n$) et leurs résidus ($d_{p,n}$) identifiés, la solution dans le domaine fréquentiel $U_p(r, \omega)$ est décomposée :
   $$U_p(r, \omega) = U^s_{p,h}(r, \omega) + \sum_{n} \frac{d_{p,n}(r)}{\omega - \rho_n}$$
   où $U^s_{p,h}$ est le terme soustrait de la singularité (régularisé). La transformée de Fourier inverse est divisée en deux parties :

   • $I_1(r, t)$ : L'intégrale du terme régularisé. Puisque les pics acérés sont supprimés, ce terme peut être évalué efficacement en utilisant des règles de quadrature FTH haute fréquence standards, même pour des temps longs.
   • $I_2(r, t)$ : L'intégrale des termes singuliers soustraits.

3. Évaluation Asymptotique des Contributions Singulières :
   Pour la partie singulière $I_2$, l'article utilise une expansion asymptotique à grand temps dérivée par déformation de contour.

   • Aux temps pré-asymptotiques, $I_2$ est évalué numériquement à l'aide d'une quadrature spécialisée qui gère les quasi-singularités.
   • Pour les temps longs ($t > T_{inc}$), l'intégrale est approximée par une « expansion de singularité » :
     $$E^I_h(r, t) = -i \sum_{n} d_{p,n}(r) A(\sigma_n) e^{-i\sigma_n t}$$
   • L'analyse théorique (Théorème 1) fournit des bornes d'erreur montrant que la différence entre l'intégrale exacte et l'expansion asymptotique décroît exponentiellement avec le temps, à condition que l'intervalle de fréquence soit suffisamment large.

Contributions Clés

• Nouvel Algorithme IE : Une méthode de calcul efficace pour déterminer les résonances complexes excitées par incidence et leurs résidus en utilisant uniquement des inversions d'intégrales de bord à fréquences réelles, évitant ainsi l'intégration de contours de fréquences complexes ou la recherche de racines dans le plan complexe.
• Régularisation par Soustraction : Une technique qui transforme une transformée de Fourier inverse quasi-singulière en une intégrale lisse et une somme de termes exponentiels connus, permettant une simulation temporelle longue et efficace pour des structures hautement piégeantes.
• Réutilisation des Approximants Rationnels : La démonstration que les approximants rationnels utilisés pour la recherche de pôles peuvent être directement utilisés pour calculer les résidus et les valeurs de champ, réduisant considérablement le coût computationnel de l'algorithme global.
• Validation de l'Expansion de Singularité : Bien que la validité des expansions de singularité pour les obstacles piégeants ne soit pas rigoureusement prouvée dans la littérature, l'article fournit des preuves numériques solides montrant que l'expansion décrit avec précision le comportement à temps long des champs diffusés pour des géométries de piégeage arbitraires.

Résultats
Les expériences numériques présentées dans l'article démontrent :

• Précision et Efficacité : La méthode FTH-SS atteint une haute précision pour la diffusion par des obstacles hautement piégeants (ex: petits cercles à ouverture étroite, cavités en forme de fusée, structures de type "whispering-gallery") là où les méthodes FTH standards échouent ou nécessitent des ressources de calcul excessives.
• Régularisation : La soustraction des résonances lisse efficacement l'intégrande dans le domaine fréquentiel, permettant une convergence rapide de la quadrature numérique.
• Capture Adaptative des Pôles : L'algorithme IE identifie avec succès toutes les résonances pertinentes pour un champ incident spécifique, filtrant celles ayant des résidus négligeables.
• Comportement à Temps Long : La méthode capture avec précision la décroissance lente des champs diffusés dans les régions de piégeage. L'expansion asymptotique montre qu'elle converge rapidement vers la solution réelle à mesure que le temps augmente, même pour des géométries complexes.
• Comparaison : Dans les cas non piégeants, FTH-SS coïncide avec le FTH standard. Dans les cas piégeants, FTH-SS surpasse significativement le FTH standard, atteignant une précision de l'ordre de la machine avec des ordres de grandeur de moins d'évaluations de fréquences.

Signification
L'article affirme que l'approche FTH-SS proposée permet des simulations de diffusion temporelle précises et efficaces pour des durées arbitrairement longues, même pour des structures hautement piégeantes. En découplant le besoin d'un échantillonnage fréquentiel dense de la présence de résonances complexes, la méthode surmonte un goulot d'étranglement majeur dans les simulations de diffusion d'ondes. De plus, les résultats numériques fournissent des preuves convaincantes soutenant la validité pratique des expansions de singularité pour les obstacles piégeants, un régime où la justification théorique était auparavant évasive. Ce travail établit un cadre robuste pour traiter les pôles « quasi-réels » qui tourmentent typiquement les approches de domaine fréquentiel pour les problèmes temporels.