A HDG method with transmission variables for time-harmonic wave propagation problems with constant coefficients

Sintesi Tecnica: Un Metodo HDG con Variabili di Trasmissione per la Propagazione di Onde in Regime di Tempo-Armonico

Definizione del Problema
I solutori iterativi a elementi finiti per problemi di propagazione ondosa in regime tempo-armonico sono notoriamente lenti nella convergenza a causa delle proprietà intrinseche dei sistemi lineari sottostanti. Sebbene i metodi Discontinuous Galerkin (DG) siano ben adatti per soluzioni ad alta fedeltità in geometrie complesse e regimi ad alta frequenza, i sistemi lineari sparsi e non strutturati risultanti pongono sfide significative per i solutori diretti su architetture parallele. Sebbene i solutori iterativi offrano efficienza di memoria, i loro tassi di convergenza sono spesso insufficienti. I metodi Hybridizable Discontinuous Galerkin (HDG) esistenti migliorano il condizionamento del sistema eliminando le incognite fisiche all'interno degli elementi, lasciando un sistema più piccolo definito sullo scheletro della mesh. Tuttavia, le formulazioni HDG standard, che utilizzano tracce numeriche come variabili ibride, non sempre garantiscono la convergenza più rapida per gli schemi iterativi nei contesti di propagazione ondosa.

Metodologia
Gli autori propongono una variante del metodo HDG, denominata CHDG (Characteristic-based HDG), che sostituisce la traccia numerica standard con variabili di trasmissione derivate dall'analisi delle caratteristiche dei sistemi iperbolici governanti.

1. Quadro Matematico: Il metodo si rivolge a sistemi iperbolici lineari, simmetrici, con coefficienti costanti nel dominio della frequenza. Le equazioni governanti vengono analizzate tramite decomposizione delle caratteristiche, suddividendo la soluzione in variabili di trasmissione uscenti ($G^+$) e entranti ($G^-$) sulla base degli autovalori della matrice del flusso.
2. Formulazione DG: Viene costruita una formulazione DG utilizzando flussi upwind espressi esplicitamente in termini di queste variabili di trasmissione. Il flusso numerico su un'interfaccia utilizza la variabile uscente dell'elemento corrente e la variabile entrante dell'elemento vicino.
3. Ibridazione: L'innovazione principale risiede nella strategia di ibridazione. Invece di eliminare le variabili fisiche tramite una traccia numerica, il metodo introduce le variabili di trasmissione entranti ($G^-$) come incognite globali sullo scheletro della mesh.
   • Le variabili fisiche all'interno di ogni elemento vengono eliminate risolvendo problemi locali dove $G^-$ è trattata come dato al contorno.
   • Ciò produce un sistema ibridato globale della forma $(I - \Pi S)G^- = b$, dove $S$ è un operatore di scattering (che mappa le variabili entranti in uscenti) e $\Pi$ è un operatore di scambio (che gestisce l'accoppiamento delle interfacce e le condizioni al contorno).
4. Proprietà Teoriche: Gli autori dimostrano che, per condizioni al contorno passive, l'operatore $\Pi S$ è una contrazione stretta. Di conseguenza, il sistema ibridato è ben posto ed è suscettibile di iterazione a punti fissi.

Applicazione all'Aeroacustica
Il framework è applicato specificamente alle equazioni di Euler linearizzate (LEE) in regime tempo-armonico per modellare la propagazione ondosa in un flusso medio subsonico uniforme. Questa applicazione introduce complicazioni specifiche:

• Il sistema supporta sia onde sonore che onde di vorticità.
• Il numero di variabili di trasmissione entranti e uscenti su una faccia dell'elemento non è costante; dipende dall'orientamento della faccia rispetto alla velocità del flusso medio ($u_0 \cdot n$).
• Il framework generale proposto gestisce naturalmente questa variabilità nel numero di condizioni al contorno e di variabili di trasmissione attraverso il dominio.

Risultati Numerici
Il metodo viene validato attraverso benchmark numerici bidimensionali che coinvolgono onde sonore e di vorticità in flussi medi uniformi, confrontando la convergenza di tre schemi iterativi (Fixed Point, CGNR e GMRES) applicati sia ai sistemi DG standard che ai sistemi CHDG:

• Propagazione di Onde Piane: Test con diverse direzioni di propagazione rispetto al flusso medio e diverse condizioni al contorno (impedenza vs. pressione).
• Sorgente Puntiforme: Un problema di sorgente singolare in un dominio circolare con diversi numeri di Mach.
• Onde di Vorticità: Propagazione in una guida d'onda rigida con modi di ordine superiore.

Risultati Chiave:

• L'iterazione a punti fissi (Fixed Point) applicata al sistema CHDG converge costantemente, allineandosi con la dimostrazione teorica della contrazione.
• Gli schemi iterativi (CGNR e GMRES) applicati al sistema CHDG convergono significativamente più velocemente (meno iterazioni) rispetto a quando applicati ai sistemi DG standard.
• Nel contesto CHDG, il metodo a punti fissi performa quasi quanto il GMRES e supera costantemente il CGNR.
• La presenza di onde di vorticità e l'orientamento variabile del flusso non ostacolano le proprietà di convergenza dell'approccio CHDG.

Significatività e Rivendicazioni
Il lavoro sostiene di estendere la metodologia CHDG, precedentemente stabilita per i problemi di Helmholtz, a un framework generale che copre una vasta classe di problemi di propagazione ondosa definiti da sistemi iperbolici simmetrici con coefficienti costanti.

• Gestione Unificata: Il framework unifica il trattamento di diversi tipi di onde (acustiche e di vorticità) all'interno di una singola struttura matematica.
• Robustezza: Gestisce efficacemente la complessità del numero variabile di variabili di trasmissione causato dall'orientamento del flusso medio, una caratteristica non gestita naturalmente dalle precedenti implementazioni specifiche.
• Efficienza: La riformulazione accelera significativamente la convergenza dei solutori iterativi rispetto agli approcci DG e HDG standard, rendendolo un candidato promettente per simulazioni di propagazione ondosa ad alta frequenza su larga scala.

Gli autori osservano che, sebbene i risultati attuali si basino su script MATLAB, il lavoro preliminare sulle onde elettromagnetiche suggerisce benefici simili, e il lavoro futuro mira a estendere questo framework alle onde elastiche e ai mezzi non omogenei.