Invariant domain preserving limiting of time explicit and time implicit discretizations for systems of conservation laws

Sintesi Tecnica: Limitazione che Preserva i Domini Invarianti per Discretizzazioni Temporali Esplicite e Implicite

Enunciato del Problema
Il lavoro affronta la sfida di progettare approssimazioni numeriche ad alto ordine (HO) per sistemi iperbolici non lineari di leggi di conservazione che preservino i domini invarianti (ID). Sebbene i metodi HO catturino con precisione le soluzioni lisce, spesso non soddisfano i vincoli fisici a livello discreto, come la positività della densità e dell'energia interna o il rispetto delle disuguaglianze di entropia. Al contrario, gli schemi a basso ordine (LO) (ad esempio, Godunov, Lax-Friedrichs) sono noti per preservare tali domini invarianti, ma soffrono di eccessiva diffusione numerica. L'obiettivo è costruire una tecnica di limitazione che vincoli la soluzione HO a rimanere all'interno dell'insieme convesso degli stati ammissibili $D$ (ad esempio, densità positiva, energia interna positiva, limiti di entropia) senza sacrificare la conservazione o l'accuratezza ad alto ordine nelle regioni lisce.

Metodologia
Gli autori propongono un quadro generale di limitazione applicabile a qualsiasi discretizzazione spaziale conservativa (Volume Finito, Galerkin Discontinuo) e a un'ampia gamma di schemi di integrazione temporale, inclusi Runge-Kutta (RK) espliciti e impliciti, metodi multistep e integrazione temporale Galerkin Discontinuo (DG).

1. Concetto Fondamentale: Il metodo limita la soluzione HO attorno a una soluzione LO che è garantita essere preservante i domini invarianti. Generalizza il quadro del Trasporto Corretto di Flusso (FCT) limitando i flussi antidiffusivi (la differenza tra i residui HO e LO), ma definisce i coefficienti di limitazione in modo da esprimere la soluzione finale come una combinazione convessa di quantità preservanti i domini invarianti, simile al quadro di limitazione convessa.
2. Coefficienti di Limitazione: A differenza dell'FCT standard che impone limiti sulle componenti della soluzione, questo approccio definisce coefficienti $l_{\kappa,\kappa'}^{i,j}$ tali che la soluzione limitata $u_h^{n+1}$ sia una combinazione convessa della soluzione LO e dell'aggiornamento HO. I coefficienti sono determinati trovando i più grandi fattori di scala $k$ tali che gli stati intermedi rimangano all'interno del dominio invariante $D$.
3. Procedura Iterativa: Il limitatore è incorporato in un algoritmo iterativo. Partendo dalla soluzione LO, il metodo affina iterativamente la soluzione limitata riapplicando il limitatore alla differenza tra la soluzione HO e la soluzione limitata corrente. Questo processo continua fino al raggiungimento di una tolleranza o di un numero massimo di iterazioni, mirando a recuperare la soluzione HO il più possibile mantenendo le proprietà di preservazione dei domini invarianti (IDP).
4. Accelerazione: Viene introdotto un fattore di accelerazione euristico $\beta \geq 1$ per scalare i coefficienti di limitazione, accelerando la convergenza del processo iterativo verso la soluzione HO.
5. Adattamento all'Integrazione Temporale: Il quadro si adatta a vari schemi temporali definendo uno schema LO a singolo passo (Euler in avanti per l'esplicito, Euler all'indietro per l'implicito) per calcolare i residui necessari e i flussi antidiffusivi, evitando la necessità di ricalcolare le soluzioni LO a ogni stadio dei metodi multistadio o multistep.
6. Discretizzazione Spaziale: Vengono fornite derivazioni specifiche per i flussi antidiffusivi per gli schemi Volume Finito (FV), gli schemi DG generali e il Metodo agli Elementi Spettrali Galerkin Discontinuo (DGSEM). Per il DGSEM, il metodo combina la limitazione basata sulla media cellulare con un limitatore di scala lineare per garantire che anche i valori nodali soddisfino i vincoli del dominio invariante.

Contributi Chiave

• Generalizzazione ai Sistemi e agli Schemi Impliciti: Il metodo estende i concetti di limitazione convessa ai sistemi di leggi di conservazione e, crucialmente, si applica alle discretizzazioni temporali implicite, un dominio in cui i precedenti approcci di limitazione convessa erano limitati agli schemi espliciti.
• Limitazione Basata sui Flussi senza Stati "Bar": L'approccio definisce coefficienti di limitazione per formare combinazioni convesse di soluzioni LO e HO senza fare affidamento su stati "bar" (stati intermedi utilizzati in alcune formulazioni di limitazione convessa). Ciò consente l'uso di qualsiasi flusso numerico coerente e conservativo nello schema HO.
• Raffinamento Iterativo: L'introduzione di una procedura iterativa permette alla soluzione limitata di convergere verso la soluzione HO, recuperando l'accuratezza nelle regioni lisce mentre preserva rigorosamente i domini invarianti.
• Euristica di Accelerazione: Viene proposto un semplice fattore di scala $\beta$ per ridurre significativamente il numero di iterazioni richieste per la convergenza.
• Ampia Applicabilità: La derivazione è presentata per discretizzazioni conservative generali, con implementazioni specifiche dettagliate per FV, DG e DGSEM, accoppiate con integratori temporali RK espliciti/impliciti, multistep e DG.

Risultati Numerici
Il lavoro valida il metodo attraverso esperimenti numerici su problemi 1D e 2D che coinvolgono equazioni iperboliche scalari e le equazioni di Eulero per gas comprimibili:

• Problemi di Riemann (Onde d'Urto Forti): Per le equazioni di Eulero comprimibili, lo schema HO senza limitazione fallisce a causa di pressione/densità negativa. Il limitatore proposto risolve con successo le onde d'urto forti in modo netto mantenendo la positività e i limiti di entropia sia per l'avanzamento temporale esplicito (ERK3) che implicito (DIRK33), anche a numeri CFL elevati (fino a 10).
• Instabilità di Kelvin-Helmholtz: Nelle simulazioni 2D con grandi passi temporali, lo schema LO smorza eccessivamente le scale piccole. Lo schema limitato recupera l'alta risoluzione dello schema HO, preservando le strutture a scala fine.
• Diffrazione d'Urto e Flusso su Gradino: Il limitatore dimostra robustezza in geometrie 2D complesse, mantenendo l'ammissibilità fisica dove lo schema HO illimitato fallisce. Raggiunge una risoluzione paragonabile a quella degli schemi espliciti con piccoli passi temporali, ma con la stabilità dei metodi impliciti.
• Equazioni Scalari: I test su trasporto lineare, l'equazione di Burgers e il problema KPP (flusso non convesso) mostrano che il limitatore rimuove efficacemente le oscillazioni spurie e impone il principio del massimo e le disuguaglianze di entropia.
• Convergenza: Si dimostra che il limitatore iterativo con il fattore di accelerazione $\beta=2$ converge alla soluzione HO significativamente più velocemente dell'approccio standard ($\beta=1$).

Significato e Affermazioni
Gli autori affermano che questo lavoro fornisce un limitatore abbastanza generale in grado di gestire qualsiasi legge di conservazione iperbolica su mesh generali con un'ampia gamma di schemi di integrazione temporale. Il significato principale risiede nel colmare il divario tra la robustezza degli schemi impliciti a basso ordine e l'accuratezza dei metodi ad alto ordine. Garantendo che lo schema numerico preservi i domini invarianti a livello discreto, il metodo consente l'uso di grandi passi temporali nelle simulazioni implicite senza il rischio di produrre soluzioni non fisiche (ad esempio, densità negativa). Il lavoro sottolinea che il limitatore mantiene la conservazione e può essere applicato iterativamente per affinare la soluzione, offrendo uno strumento pratico per simulazioni robuste e ad alto ordine di problemi complessi di fluidodinamica. Il lavoro conclude notando che le future estensioni considereranno l'avanzamento temporale di Euler all'indietro linearizzato per problemi di stato stazionario.