Time integration as filtering: a space-time discretization-aware LES formulation

Sintesi Tecnica: L'integrazione temporale come filtraggio nella discretizzazione spazio-temporale consapevole della LES

Definizione del problema
Nella Large Eddy Simulation Implicita (ILES), il filtro spaziale è spesso identificato con la griglia computazionale stessa. Le tradizionali formulazioni LES "consapevoli della discretizzazione" affrontano il problema della chiusura della turbolenza derivando uno stress residuo esatto basandosi sulla proprietà di "scambio del filtro" (filter-swap): la derivata centrale discreta grossolana viene riconosciuta come una derivata spaziale esatta filtrata con un top-hat. Tuttavia, le formulazioni esistenti trattano l'integrazione temporale come un passaggio separato e di post-processing. Di conseguenza, l'errore di discretizzazione temporale rimane al di fuori della garanzia di stress residuo esatto. Man mano che i passi temporali aumentano (numeri CFL più elevati), le chiusure standard basate solo sullo spazio degradano perché non riescono a tenere conto dell'errore specifico introdotto dall'integratore temporale, specialmente quando si utilizzano schemi espliciti come il Forward Euler.

Metodologia
Il documento estende il concetto di scambio del filtro spaziale allo spazio-tempo, osservando che la differenza temporale del Forward-Euler è essa stessa una derivata temporale esatta filtrata (top-hat, unidirezionale).

1. Scambio del filtro spazio-temporale: Applicando un filtro spazio-temporale $(\cdot)_{h,\tau} = (\cdot)_\tau \circ (\cdot)_h$ alla legge di conservazione scalare, gli autori derivano una legge discreta esatta in cui sia le derivate spaziali che quelle temporali sono sostituite dalle loro controparti filtrate.
2. Decomposizione del flusso: Il flusso spazio-temporale esatto è decomposto nel flusso numerico più un flusso residuo, $r_{h,\tau}$. Questo residuo è suddiviso in quattro termini di errore distinti:
   • $r^{LES}_h$: Il classico commutatore LES (filtraggio della griglia rispetto alla non linearità).
   • $r^{num}_h$: L'errore del flusso numerico.
   • $r^{div}_h$: Una correzione per la divergenza discreta.
   • $r^{time}_\tau$: Un nuovo termine che rappresenta l'errore di quadratura temporale del flusso (la differenza tra il flusso mediato nel passo temporale e il flusso istantaneo).
3. Modellazione: Gli autori propongono un modello di viscosità di eddy consapevole della discretizzazione spazio-temporale. Espandono il termine di errore di quadratura temporale, $r^{time}_\tau$, dimostrando che si comporta come un flusso diffusivo proporzionale a $\tau (f'(u))^2 \partial_x u$. Questo viene identificato come una diffusione di tipo Lax–Wendroff, che è costruita per essere non negativa e scala linearmente con il passo temporale $\tau$. La chiusura finale aumenta il modello di Smagorinsky standard con questo termine temporale.

Risultati Chiave
Il framework è stato testato su una turbolenza di Burgers viscosa utilizzando uno schema a volumi finiti Forward-Euler con numeri CFL variabili (da 0,06 a 0,89).

• Analisi a priori: Nella decomposizione del flusso residuo, il termine temporale ($r^{time}_\tau$) cresce significativamente con il numero CFL, diventando il singolo termine di errore più grande ai valori di CFL pratici (ad esempio, 0,59 del residuo al CFL 0,89 per i flussi centrali). Al contrario, i termini spaziali mantengono proporzioni fisse indipendentemente dal passo temporale.
• Analisi a posteriori:
  • Flusso Centrale: Le chiusure standard (senza modello, Smagorinsky classico e Smagorinsky basato solo sullo spazio) degradano rapidamente all'aumentare del numero CFL, con errori che salgono da 0,1 a ~0,85. La chiusura consapevole dello spazio-tempo proposta mantiene un errore relativo basso e costante (0,09) attraverso l'intero intervallo di numeri CFL.
  • Flusso Upwind: Il flusso numerico upwind è intrinsecamente sovradisciplinante. In questo caso, l'errore temporale (che è non dissipativo) compensa parzialmente la sovradisciplina del flusso upwind. Sebbene il modello spazio-temporale rimanga accurato, è leggermente sovradisciplinante per passi temporali ampi, poiché modella gli effetti di cancellazione in modo indipendente invece di sfruttare la fortuita cancellazione.

Significatività e Rivendicazioni
Il documento sostiene che riconoscere la differenza temporale del Forward-Euler come un filtro top-hat permette una decomposizione esatta del flusso discreto nello spazio-tempo. Il contributo primario è l'identificazione dell'errore di quadratura temporale del flusso come un termine distinto e dominante nel residuo a passi temporali grossolani.

• Gli autori dimostrano che una chiusura aumentata con il termine principale di questo errore (diffusione di tipo Lax–Wendroff) permette alle simulazioni Forward-Euler di rimanere accurate anche a passi temporali elevati dove le chiusure basate solo sullo spazio falliscono.
• Il lavoro stabilisce che l'integratore temporale contribuisce esattamente con un nuovo termine al residuo e nessun termine di commutazione temporale separato.
• Il documento nota con modestia che gran parte dell'errore di quadratura temporale è non dissipativo e quindi invisibile ai modelli standard di viscosità di eddy, suggerendo che le chiusure apprese tramite machine learning, addestrate sul target esatto o su integratori impliciti/multi-stadio, siano naturali estensioni future.