Trefftz methods with evanescent plane waves

1. Il Problema: Instabilità Numerica nei Metodi Trefftz Classici

Il lavoro affronta le limitazioni fondamentali dei metodi Trefftz classici applicati all'equazione di Helmholtz ($\Delta u + \kappa^2 u = 0$).

• Contesto: I metodi Trefftz approssimano le soluzioni delle PDE utilizzando spazi discreti composti da funzioni che soddisfano esattamente l'equazione differenziale all'interno di ogni sottodomino. Per l'equazione di Helmholtz, le basi tradizionali utilizzano onde piane propagative (PPW) della forma $e^{i\kappa d \cdot x}$ con $d \in \mathbb{R}^n$.
• La Criticità: Sebbene le PPW offrano tassi di convergenza rapidi, le loro implementazioni soffrono di forti instabilità numeriche.
• Analisi Teorica: Il paper cita un risultato teorico (Teorema 1, basato su [11] e [5]) che dimostra come, per certe soluzioni di Helmholtz (in particolare le "modi evanescenti" o onde ad alto indice), un'approssimazione accurata tramite PPW richieda necessariamente coefficienti di ampiezza enorme.
  • Se l'errore target è $\delta$, i coefficienti $\mu$ necessari possono avere norma $\|\mu\| \ge M$, dove $M$ può superare l'inverso della precisione di macchina.
  • Questo porta a una cancellazione numerica inevitabile nell'aritmetica dei computer, rendendo impossibile rappresentare stabilmente l'approssimazione, indipendentemente dalle tecniche di algebra lineare (precondizionamento, decomposizioni) utilizzate.

2. Metodologia: Onde Piane Evanescenti (EPW)

Per superare l'instabilità delle PPW, gli autori propongono l'uso di Onde Piane Evanescenti (EPW) come funzioni di base.

• Definizione: Le EPW generalizzano le PPW permettendo vettori di propagazione complessi $d \in \mathbb{C}^n$ tali che $d \cdot d = 1$.
  • La forma è data da: $EW(x) = e^{i\kappa \zeta p \cdot x} e^{-\kappa \eta e \cdot x}$.
  • Qui $p$ è la direzione di propagazione, $e$ è la direzione di decadimento (ortogonale a $p$), $\eta$ è la forza di decadimento e $\zeta = \sqrt{1+\eta^2}$ è il numero d'onda apparente.
• Vantaggio Teorico: A differenza delle PPW, le EPW permettono una rappresentazione di Herglotz continua di qualsiasi soluzione di Helmholtz in $H^1(B)$ con coefficienti limitati. Questo garantisce l'esistenza di approssimazioni accurate con coefficienti piccoli, risolvendo il problema di instabilità alla radice.

3. Strategia di Selezione delle Basi (Ricetta di Campionamento)

Un contributo pratico significativo è la proposta di una strategia semplice e concreta per selezionare i parametri delle EPW in uno spazio discreto di dimensione $P$ (DOF - gradi di libertà).

• Parametri: Per ogni elemento della mesh, si campionano:
  • Direzione di propagazione $\theta$ (uniforme su $[0, 2\pi)$).
  • Direzione di decadimento $\phi$ (uno dei due vettori ortogonali).
  • Forza di decadimento $\eta$ (o $\zeta$).
• La "Ricetta" (in 2D):
  • Si definisce un livello di troncamento di Fourier $L \approx P/4$.
  • Si campiona un parametro $\xi \in [0, 1]$ e si calcola $\zeta_p = \max\left(1, \frac{2L}{\kappa \cdot \text{diam}(K)} \xi_p\right)$.
  • Se $2L \le \kappa \cdot \text{diam}(K)$, tutte le basi sono PPW ($\zeta=1$).
  • Per $P$ più grandi, una frazione crescente di basi diventa EPW ($\zeta > 1$), permettendo di catturare le componenti evanescenti della soluzione.
• Estensione 3D: Viene proposta una ricetta analoga con $L = \sqrt{P/8}$.

4. Formulazione Variazionale e Regularizzazione

Le basi EPW sono integrate nella Formulazione Variazionale Ultraweak (UWVF), un caso particolare dei metodi Discontinuous Galerkin (DG).

• Sovracampionamento (Oversampling): Per garantire la stabilità anche quando le funzioni di base sono quasi linearmente dipendenti, si utilizza un approccio di sovracampionamento (es. 10% in più di funzioni di base rispetto ai gradi di libertà finali).
• Regularizzazione: Si utilizza la SVD tronca della matrice diagonale $D$ (che rappresenta i termini locali). Si inverte solo la parte della matrice con valori singolari superiori a una soglia $\epsilon$ (es. $10^{-14}$), eliminando le componenti che causano instabilità numerica. Questo sistema regolarizzato viene risolto efficientemente in parallelo per ogni elemento.

5. Risultati Numerici

Gli esperimenti confrontano le PPW e le EPW su tre scenari:

1. Soluzione Fondamentale (Punto Sorgente):
   • Con un numero moderato di DOF, PPW ed EPW convergono in modo simile.
   • All'aumentare dei DOF, l'errore delle PPW si stabilizza (stallo) a causa dell'instabilità numerica.
   • Le EPW continuano a convergere, raggiungendo errori fino a $10^{-6}$ volte inferiori rispetto alle PPW.
2. Influenza del Numero d'Onda ($\kappa$):
   • Mantenendo il numero di onde per cella costante ($P \propto \kappa$), l'errore delle EPW migliora all'aumentare di $\kappa$, mentre le PPW mostrano degradazione.
3. Problema di Scattering (Ostacolo Non Convesso):
   • In presenza di singolarità agli angoli e di una cavità, le PPW falliscono nel raggiungere alta precisione.
   • Le EPW forniscono una precisione significativamente superiore, gestendo meglio le complessità geometriche e le componenti evanescenti generate dagli angoli.

6. Significato e Conclusioni

Il lavoro dimostra che l'instabilità dei metodi Trefftz per l'equazione di Helmholtz non è un difetto intrinseco del metodo, ma una conseguenza della scelta della base (PPW).

• Innovazione: L'introduzione delle onde piane evanescenti (EPW) con una strategia di campionamento sistematica risolve il problema della "cancellazione numerica" permettendo approssimazioni stabili con coefficienti piccoli.
• Impatto: Questo approccio rende i metodi Trefftz (in particolare UWVF) molto più robusti per problemi ad alto numero d'onda e geometrie complesse, superando i limiti delle implementazioni tradizionali basate solo su onde propagative.
• Disponibilità: I dettagli completi e ulteriori esempi sono previsti nella pubblicazione di riferimento [6] (citata nel testo).