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

Sintesi Tecnica: Sintesi di Scattering nel Dominio del Tempo Efficiente tramite Sottrazione di Singolarità nel Dominio della Frequenza

Problematica
Il documento affronta la sfida di simulare accuratamente lo scattering ondulatorio nel dominio del tempo da parte di ostacoli che presentano geometrie di tipo "trapping" (intrappolanti). Sebbene i metodi "Frequency-Time Hybrid" (FTH), basati sulla trasformata di Fourier, offrano soluzioni accurate e prive di dispersione per problemi di scattering generali risolvendo le equazioni di Helmholtz a frequenze discrete e performing una trasformata di Fourier inversa, essi soffrono di un grave degrado delle prestazioni in presenza di ostacoli intrappolanti. In tali casi, la soluzione di Helmholtz nel dominio della frequenza presenta risonanze (poli) quasi reali e complesse. Questi poli si manifestano come picchi acuti lungo l'asse della frequenza reale, rendendo l'integrando della trasformata di Fourier quasi singolare. Di conseguenza, la quadratura numerica standard richiede un campionamento di frequenza proibitivamente denso per raggiungere la convergenza, rendendo computazionalmente impraticabili le simulazioni a lungo termine. Inoltre, il lento decadimento temporale dei campi diffusi nelle regioni di intrappolamento pone sfide per le tecniche tradizionali nel dominio del tempo, come i metodi a differenze finite o a elementi finiti.

Metodologia
Gli autori propongono un metodo FTH con Sottrazione di Singolarità (FTH-SS). La strategia centrale consiste nel regolarizzare l'integrando della trasformata di Fourier inversa sottraendo i contributi delle risonanze complesse nel dominio della frequenza e reincorporandoli tramite un'efficiente espansione asintotica per tempi lunghi. La metodologia si compone di tre componenti primarie:

1. Algoritmo di Incidence-Excitation (IE):
   Per evitare l'alto costo della valutazione dei risolventi su contorni di frequenza complessi, il documento introduce un nuovo algoritmo adattivo che identifica le risonanze complesse e i loro residui utilizzando solo dati a frequenza reale.

   • Il metodo impiega l'algoritmo di approssimazione razionale AAA (Adaptive Antoulas-Anderson) applicato alle densità dell'equazione integrale eccitate da un campo incidente specifico.
   • Utilizzando un approccio di "sketching casuale" per funzioni vettoriali, l'algoritmo costruisce approssimanti razionali per le soluzioni di densità.
   • Una strategia ricorsiva adattiva partiziona l'intervallo di frequenza per garantire che tutti i poli rilevanti (quelli con residui sufficientemente significativi da causare picchi nella soluzione) siano catturati.
   • Fondamentalmente, gli approssimanti razionali generati durante l'identificazione dei poli vengono riutilizzati per calcolare i residui e i valori del campo senza richiedere ulteriori e costose inversioni di integrali di bordo.

2. Sottrazione della Singolarità:
   Una volta identificati i relativi poli complessi ($\rho_n$) e i loro residui ($d_{p,n}$), la soluzione nel dominio della frequenza $U_p(r, \omega)$ viene decomposta come:
   $$U_p(r, \omega) = U^s_{p,h}(r, \omega) + \sum_{n} \frac{d_{p,n}(r)}{\omega - \rho_n}$$
   dove $U^s_{p,h}$ è il termine sottratto della singolarità (regolarizzato). La trasformata di Fourier inversa è divisa in due parti:

   • $I_1(r, t)$: L'integrale del termine regolarizzato. Poiché i picchi acuti sono stati rimossi, questo termine può essere valutato efficientemente utilizzando regole di quadratura FTH ad alta frequenza standard, anche per tempi lunghi.
   • $I_2(r, t)$: L'integrale dei termini singolari sottratti.

3. Valutazione Asintotica dei Contributi Singolari:
   Per la parte singolare $I_2$, il documento utilizza un'espansione asintotica per tempi lunghi derivata dalla deformazione del contorno.

   • Per tempi pre-asintotici, $I_2$ viene valutato numericamente utilizzando una quadratura specializzata che gestisce le quasi-singolarità.
   • Per tempi lunghi ($t > T_{inc}$), l'integrale è approssimato da un' "espansione di singolarità":
     $$E^I_h(r, t) = -i \sum_{n} d_{p,n}(r) A(\sigma_n) e^{-i\sigma_n t}$$
   • L'analisi teorica (Teorema 1) fornisce limiti di errore dimostrando che la differenza tra l'integrale esatto e l'espansione asintotica decade esponenzialmente con il tempo, a condizione che l'intervallo di frequenza sia sufficientemente ampio.

Contributi Chiave

• Nuovo Algoritmo IE: Un metodo computazionalmente efficiente per determinare le risonanze complesse eccitate dall'incidenza e i relativi residui utilizzando solo inversioni di integrali di bordo a frequenza reale, evitando la necessità di integrazione su contorni di frequenza complessi o la ricerca di radici nel piano complesso.
• Regolarizzazione tramite Sottrazione: Una tecnica che trasforma una trasformata di Fourier inversa quasi singolare in un integrale liscio e in una somma di termini esponenziali noti, consentendo simulazioni efficienti a lungo termine per strutture altamente intrappolanti.
• Riutilizzo degli Approssimanti Razionali: La dimostrazione che gli approssimanti razionali utilizzati per la ricerca dei poli possono essere direttamente utilizzati per calcolare i residui e i valori del campo, riducendo significativamente il costo computazionale dell'intero algoritmo.
• Validazione dell'Espansione di Singolarità: Sebbene la validità delle espansioni di singolarità per ostacoli intrappolanti non sia stata rigorosamente provata in letteratura, il documento fornisce forti prove numeriche che l'espansione descrive accuratamente il comportamento a lungo termine dei campi diffusi per arbitrarie geometrie di intrappolamento.

Risultati
Gli esperimenti numerici presentati nel documento dimostrano:

• Accuratezza ed Efficienza: Il metodo FTH-SS raggiunge un'elevata accuratezza per lo scattering da parte di ostacoli altamente intrappolanti (ad esempio, cerchi a piccola apertura, cavità a forma di razzo, strutture a galleria di sussurro) dove i metodi FTH standard falliscono o richiedono risorse computazionali eccessive.
• Regolarizzazione: La sottrazione delle risonanze regolarizza efficacemente l'integrando nel dominio della frequenza, permettendo una rapida convergenza della quadratura numerica.
• Cattura Adattiva dei Poli: L'algoritmo IE identifica con successo tutte le risonanze rilevanti per un determinato campo incidente, filtrando quelle con residui trascurabili.
• Comportamento a Lungo Termine: Il metodo cattura accuratamente il lento decadimento dei campi diffusi nelle regioni di intrappolamento. Si osserva che l'espansione asintotica converge rapidamente alla soluzione reale all'aumentare del tempo, anche per geometrie complesse.
• Confronto: Nei casi non intrappolanti, il FTH-SS coincide con il FTH standard. Nei casi di intrappolamento, il FTH-SS supera significativamente il FTH standard, raggiungendo un'accuratezza vicina alla precisione di macchina con ordini di grandezza in meno di valutazioni di frequenza.

Significatività
Il documento sostiene che l'approccio FTH-SS proposto consente simulazioni di scattering nel dominio del tempo accurate ed efficienti per tempi arbitrariamente lunghi, anche per strutture altamente intrappolanti. Decoupling (disaccoppiando) la necessità di un campionamento denso della frequenza dalla presenza di risonanze complesse, il metodo supera un importante collo di bottiglia nelle simulazioni di scattering ondulatorio. Inoltre, i risultati numerici forniscono prove convincenti a supporto della validità pratica delle espansioni di singolarità per gli ostacoli intrappolanti, un regime in cui la giustificazione teorica è stata precedentemente elusiva. Il lavoro stabilisce un framework robusto per gestire i poli "quasi reali" che tipicamente affliggono gli approcci nel dominio della frequenza per problemi nel dominio del tempo.