Investigation of Shock-Capturing with Bound-Preserving Limiters for the Nonlinearly Stable Flux Reconstruction Method

Sintesi Tecnica: Investigazione del Capturing di Shock con Limiter Preservanti i Limiti per il Metodo di Ricostruzione del Flusso Nonlineamente Stabile

Definizione del Problema
Gli schemi numerici ad alto ordine per i flussi comprimibili, come il metodo Discontinuous Galerkin (DG) e la Ricostruzione del Flusso (FR), offrono un'elevata accuratezza ma spesso faticano in termini di robustezza in presenza di shock forti. Sebbene lo schema Entropy-Stable Flux Reconstruction (ESFR) e la sua estensione, Nonlinearly Stable Flux Reconstruction (NSFR), forniscano una stabilità non lineare e una stabilità entropica provabili, queste proprietà da sole non garantiscono la preservazione dei limiti fisici, come densità e pressione positive. Gli schemi standard ad alto ordine spesso producono oscillazioni spurie o valori non fisici in prossimità degli shock, portando alla divergenza del solutore. Le strategie esistenti per mitigare questi problemi, come la viscosità artificiale o il limiting subcell, possono alterare l'ordine di accuratezza, introdurre limitazioni sul passo temporale o richiedere implementazioni complesse. Inoltre, le specifiche condizioni di Courant-Friedrichs-Lewy (CFL) richieste per mantenere la positività per i flussi a due punti utilizzati nei metodi FR non erano state pienamente stabilite.

Metodologia
Il documento investiga l'applicazione di un limiter modificato che preserva i limiti al metodo NSFR per la risoluzione delle equazioni di Eulere comprimibili. La metodologia si concentra su tre aree primarie:

1. Limiter Modificato Preservante i Limiti: Gli autori estendono il limiter Zhang-Shu preservante i limiti. Il limiter originale e i suoi successivi miglioramenti apportati da Wang e Shu garantivano la positività in specifici punti di quadratura, ma mancavano di garanzie di robustezza ai nodi della soluzione, permettendo potenzialmente valori di densità o pressione negativi in tali nodi. La modifica proposta incorpora i valori minimi di densità e pressione ai nodi della soluzione (oltre ai punti di quadratura standard) quando si determinano i fattori di scala ($\theta_1$ e $\theta_2$) utilizzati per limitare la soluzione. Ciò garantisce la positività in tutto l'elemento, inclusi i nodi della soluzione.
2. Configurazione dello Schema NSFR: Lo studio utilizza lo schema NSFR con una forma di differenziazione del flusso a due punti skew-symmetric. L'implementazione impiega il flusso a due punti Chandrashekar-Ranocha (CHRA), che è conservativo dell'entropia, preservante l'energia cinetica e preservante l'equilibrio di pressione. La discretizzazione temporale utilizza il metodo Runge-Kutta a Stabilità Forte (SSPRK3) del terzo ordine.
3. Investigazione Parametrica: Il documento analizza sistematicamente l'impatto di vari parametri dello schema sulle capacità di shock-capturing:
   • Flussi a due punti: Confronto tra i flussi CHRA, Chandrashekar (CH), Ismail-Roe (IR) e Kennedy-Gruber (KG).
   • Nodi di Quadratura: Confronto tra i nodi Gauss-Legendre (GL) e Gauss-Lobatto-Legendre (GLL) per l'evaluazione del flusso.
   • Parametro di Ricostruzione del Flusso (FR) ($c$): Investigazione di diversi valori del parametro di correzione $c$, che variano dal valore standard DG ($c_{DG}=0$) al differenziale spettrale ($c_{SD}$), all' $g_2$ di Huynh ($c_{HU}$), e al valore massimo stabile ($c_+$), nonché valori superiori a $c_+$.

Contributi Chiave

• Robustezza del Limiter Migliorata: Il contributo primario è la modifica del limiter Zhang-Shu per tenere conto dei valori dei nodi della soluzione, il che si è dimostrato necessario per la simulazione con successo di problemi 2D dominati da shock forti senza fallimenti.
• Verifica della Condizione CFL: Il documento stabilisce numericamente che i flussi a due punti utilizzati in NSFR aderiscono a specifiche condizioni CFL per la preservazione della positività. Nello specifico, viene riscontrato che una condizione CFL basata sulla distanza minima tra i nodi di quadratura GLL (Eq. 29) è sufficiente per la preservazione della positività attraverso vari gradi polinomiali sia in 1D che in 2D.
• Capacità di Shock-Capturing di NSFR: Lo studio dimostra che lo schema NSFR, quando accoppiato con il limiter modificato, può risolvere shock forti con soluzioni essenzialmente prive di oscillazioni, particolarmente quando si utilizzano parametri di ricostruzione del flusso più forti (ad esempio, $c_+$ o superiori).
• Analisi Comparativa: Il lavoro fornisce un confronto esaustivo tra NSFR e i metodi standard Strong DG, evidenziando la superiore robustezza e i vantaggi nel passo temporale dell'approccio NSFR.

Risultati
La metodologia proposta è stata validata attraverso una serie di casi di test 1D e 2D:

• Accuratezza: Il limiter modificato preserva l'elevata accuratezza dello schema, come verificato dagli studi di convergenza sul caso Low Density.
• Robustezza e Vantaggi del Passo Temporale: Nei problemi Sod Shock Tube e Shu-Osher, lo schema NSFR con il limiter modificato ha dimostrato una robustezza significativamente maggiore rispetto al metodo Strong DG. Ad esempio, nel Sod Shock Tube, NSFR poteva operare senza il limiter a un CFL di 0,2, mentre lo Strong DG falliva a CFL anche di 0,01 senza un limiter. Nello Strong Vortex Shock Wave interaction, NSFR ha permesso un CFL di 0,58 con il limiter, rispetto a 0,27 senza di esso, mantenendo l'accuratezza della soluzione.
• Mitigazione delle Oscillazioni: È stato riscontrato che l'aumento del parametro di correzione FR ($c$) mitiga efficacemente le oscillazioni spurie e gli overshoots nella scia degli shock. Ad esempio, nei casi Leblanc Shock Tube e Double Mach Reflection, lo schema $c_+$ ha prodotto soluzioni più pulite con rumore ridotto rispetto allo schema $c_{DG}$.
• Casi Estremi: Lo schema ha simulato con successo benchmark estremi, inclusi i casi di Astrophysical Jet Mach 80 e Mach 2000. Nel caso Mach 2000, un parametro di ricostruzione del flusso dieci volte superiore a $c_+$ ha permesso alla simulazione di girare con un CFL di 0,0005, preservando la positività dove i metodi DG standard fallivano.
• Sensibilità di Flusso e Nodo: Il flusso Kennedy-Gruber (KG) è risultato insufficientemente dissipativo per shock forti, portando a crescita di entropia e oscillazioni. Mentre i nodi GL offrivano una quadratura più forte, i nodi GLL hanno fornito un vantaggio significativo nei passi temporali ammissibili grazie alla ridotta sensibilità nelle proiezioni di entropia ai bordi.

Significatività
Il documento sostiene che lo schema NSFR modificato, combinato con il limiter potenziato che preserva i limiti, offre un framework robusto per simulare flussi comprimibili non stazionari con shock forti. La significatività risiede nella capacità di ottenere soluzioni essenzialmente prive di oscillazioni senza ricorrere a limiter TVD o viscosità artificiale, che spesso degradano l'accuratezza. Lo studio evidenzia come lo schema NSFR non solo mantenga la stabilità entropica, ma fornisca anche vantaggi pratici in termini di passi temporali ammissibili e fedeltà di shock-capturing. Gli autori suggeriscono che la capacità di scalare il parametro FR in base alle discontinuità locali rappresenta una direzione promettente per la ricerca futura, mirando a massimizzare l'efficienza del passo temporale preservando l'accuratezza d'alto ordine nelle regioni fluide (smooth). Il lavoro stabilisce una base per l'applicazione di questi schemi a una gamma più ampia di flussi, inclusi problemi di shock viscosi e interazioni shock-turbolenza.