Convergence of finite element right-hand-side computation from finite difference data

Sintesi Tecnica: Convergenza del Calcolo del Termine Sorgente Destro degli Elementi Finiti da Dati di Differenze Finite

Definizione del Problema
Nelle simulazioni aeroacustiche (CAA) ibride e nei problemi di campo accoppiato, è richiesto un trasferimento accurato dei campi tra griglie strutturate di differenze finite (FD) (spesso utilizzate per la simulazione del flusso) e mesh del metodo degli elementi finiti (FEM) non strutturate (utilizzate per la propagazione acustica). La sfida principale risiede nell'assemblaggio del vettore del termine destro (RHS) del FEM, definito come $b_i = \int_{\Omega} N_i(x)f(x) d\Omega$, dove $f(x)$ è un campo sorgente definito sulla griglia FD.

Gli approcci standard spesso si affidano a una semplice interpolazione, il che può introdurre errori significativi, particolarmente quando i termini sorgente presentano transizioni nette o forti oscillazioni spaziali (ad esempio, termini sorgente di tipo Lighthill). Il documento affronta la necessità di metodi di integrazione che garantiscano accuratezza, conservazione dell'energia (globale e locale) e convergenza all'aumentare della risoluzione della mesh, senza essere limitati dalle assunzioni di regolarità (smoothness) del campo sorgente.

Metodologia
Il documento propone e analizza due distinte strategie di integrazione per il trasferimento dei dati FD al RHS del FEM:

1. Integrazione Gaussiana ad Alto Ordine con Interpolazione del Campo:

   • Meccanismo: Questo metodo valuta l'integrale utilizzando la quadratura di Gauss all'interno di ogni elemento FEM. I valori del campo FD vengono prima interpolati nei punti di quadratura di Gauss utilizzando B-Spline ad alto ordine (ordini $p \in \{1, \dots, 5\}$) o polinomi Lagrangiani (ordini $p \in \{1, 3\}$).
   • Analisi dell'Errore: L'errore totale viene scomposto in errore di interpolazione ($\varepsilon_{interp}$), errore di quadratura ($\varepsilon_{quad}$) ed errore di discretizzazione FEM ($\varepsilon_{fem}$). Lo studio dimostra che, affinché il metodo converga, l'ordine di interpolazione e l'ordine di quadratura devono essere aumentati in modo coerente. Se l'ordine di quadratura viene aumentato senza un corrispondente aumento dell'ordine di interpolazione, l'errore di interpolazione diventa il fattore limitante, impedendo la convergenza alla precisione di macchina.

2. Integrazione Cut-Cell (Supermesh) con Clipping Geometrico:

   • Meccanismo: Questo approccio evita gli errori di interpolazione a livello di griglia-mesh trattando il problema come un'intersezione geometrica di due mesh non conformi. Le celle FD e gli elementi FEM vengono sovrapposti per creare "poligoni tagliati" (supermesh).
   • Implementazione: Il dominio viene decomposto in poligoni di intersezione, che vengono poi tessellati in triangoli. L'integrazione viene eseguita tramite quadratura di Gauss su questi sotto-triangoli nello spazio fisico.
   • Complessità: Una sfida computazionale significativa è la mappatura inversa necessaria per valutare le funzioni di forma FEM nei punti di Gauss fisici situati all'interno delle regioni tagliate. Per mappature geometriche non lineari, ciò richiede la risoluzione di un sistema non lineare (ad esempio, tramite iterazione Newton-Raphson) per ogni punto di quadratura.
   • Ottimizzazione: L'algoritmo utilizza l'indicizzazione AABB (Axis-Aligned Bounding Box) e la ricerca binaria per identificare efficientemente le celle FD intersecanti, seguita dal clipping dei poligoni (ad esempio, Sutherland-Hodgman) per definire i domini di integrazione.

Risultati Chiave

• Comportamento di Convergenza:
  • Il metodo basato sull'interpolazione converge all'aumentare dell'ordine delle B-Spline e dell'ordine di quadratura, a condizione che la risoluzione della mesh FD sia coerente. Tuttavia, presenta un limite sistematico di errore determinato dall'ordine di ricostruzione dell'interpolazione. Fatica a raggiungere la precisione di macchina per termini sorgente altamente oscillanti, a meno di non utilizzare ordini estremamente elevati.
  • Il metodo supermesh elimina gli errori di interpolazione e di approssimazione all'interfaccia griglia-mesh. Esso converge alla precisione di macchina in tutte le risoluzioni di mesh, indipendentemente dal fatto che il campo sorgente sia regolare o fortemente oscillante.
• Sensibilità del Termine Sorgente: Studi comparativi utilizzando termini sorgente di tipo Lighthill (varianti spazialmente, con transizioni nette) e termini dell'equazione d'onda convettiva perturbata (più regolari) mostrano che l'approccio supermesh è universalmente robusto. I metodi basati sull'interpolazione mostrano errori più elevati per i termini di tipo Lighthill più complessi.
• Prestazioni: Sebbene il metodo supermesh comporti una maggiore complessità geometrica e il costo della mappatura inversa, esso è altamente ottimizzato. Per problemi che coinvolgono molti elementi, può superare i metodi di interpolazione non ottimizzati. Lo scaling debole dell'integrazione è lineare per entrambi i metodi.

Significato e Rivendicazioni
Il documento sostiene che l'approccio di integrazione supermesh fornisce un'alternativa superiore per le simulazioni CAA ibride dove la conservazione dell'energia del termine sorgente è critica. Eliminando il passaggio di interpolazione tra la griglia FD e i punti di Gauss del FEM, il metodo assicura che l'intera energia dal campo di flusso venga trasferita alla mesh acustica senza riduzione o smoothing artificiale.

Gli autori sottolineano che, sebbene l'interpolazione ad alto ordine possa ridurre gli errori, non può garantire una convergenza sufficiente in pratica se le quantità di campo non sono "sufficientemente regolari". Pertoco, l'approccio supermesh viene presentato come una soluzione robusta che è "universale" sia per quantità di campo regolari che fortemente oscillanti. Il lavoro contribuisce con un'analisi comparativa rigorosa di questi schemi di integrazione, evidenziando come la scelta del metodo dipenda fortemente dalla natura del termine sorgente e dal requisito di una stretta conservazione.

Conclusione
Lo studio conclude che, per un trasferimento accurato del campo nelle aeroacustiche ibride, il metodo di integrazione supermesh è l'approccio preferibile quando la conservazione e l'alta accuratezza sono fondamentali, poiché rimuove gli errori sistematici inerenti agli schemi basati sull'interpolazione. Il metodo basato sull'interpolazione rimane vitale per termini sorgente più regolari, ma richiede un attento bilanciamento tra gli ordini di interpolazione e di quadratura per minimizzare l'errore.