Exact expressions for the unresolved stress in a finite-volume based large-eddy simulation

Sintesi Tecnica: Espressioni Esatte per lo Stress Non Risolto in una Simulazione Large-Eddy Basata su Volumi Finiti

Definizione del Problema
La simulazione Large-Eddy (LES) mira a modellare i flussi turbolenti risolvendo le scale grandi e modellando gli effetti delle sottoscale non risolte. Tradizionalmente, il problema della chiusura nella LES si concentra sull'errore di commutatore derivante dalla non commutazione tra il filtraggio spaziale e i termini di advezione non lineari, che risulta nel classico tensore dello stress residuo (RST), $\tau_{ij} = \overline{u_i u_j} - \bar{u}_i \bar{u}_j$. Tuttavia, nelle simulazioni pratiche, la LES viene risolta utilizzando schemi di discretizzazione numerica, comunemente il Metodo dei Volumi Finiti (FVM).

Un problema critico è che la larghezza del filtro $\Delta$ nella LES è spesso dello stesso ordine di grandezza della spaziatura della griglia $h$ (ovvero, $\Delta \approx h$). In tali regimi, gli errori di discretizzazione — derivanti dai flussi numerici e dagli operatori di divergenza discreti — sono dello stesso ordine degli stress sub-filtro fisici. Gli approcci classici tipicamente definiscono il modello di chiusura basandosi su equazioni continue e poi discretizzano, oppure trattano gli errori di discretizzazione come problemi di convergenza separati. Questa separazione porta a un disallineamento tra i dati target (il residuo vero) e il modello, specialmente quando la discretizzazione stessa introduce effetti significativi non locali o non simmetrici che i modelli di chiusura locali e simmetrici standard (come lo Smagorinsky) non possono catturare.

Metodologia
Gli autori propongono un quadro matematico unificato denominato LES-FVM, che integra direttamente il processo di filtraggio LES con la discretizzazione FVM. La metodologia centrale prevede:

1. Formalismo Unificato: Invece di applicare il filtraggio LES alle equazioni continue e poi discretizzare, gli autori applicano il filtro FVM (una media volumetrica sui volumi di controllo) alle equazioni filtrate LES. Ciò crea un sistema "doppio-filtrato" dove la larghezza del filtro $\Delta$ e la spaziatura della griglia $h$ sono trattate come parametri indipendenti.
2. Derivazione dei Residui Esatti: Utilizzando una proprietà di "scambio del filtro discreto" (dove l'operatore di differenza finita è equivalente a una derivata esatta filtrata), gli autori derivano espressioni esatte per il flusso residuo nelle equazioni discrete.
3. Decomposizione dei Residui: Il flusso residuo totale $\tau_{\Delta,h}$ è scomposto in tre componenti distinte:
   • Termine Classico: L'errore di commutatore standard associato al doppio filtro ($\Delta$ e $h$).
   • Termine di Flusso: Errori derivanti dalla differenza tra il flusso esatto e il flusso numerico ($r$ rispetto a $r_h$).
   • Termine di Divergenza: Errori derivanti dalla differenza tra la divergenza esatta e la divergenza discreta.
4. Estensione al Flusso Incomprimibile 3D: Il framework viene esteso alle equazioni di Navier-Stokes incomprimibili 3D. Ciò richiede la gestione del termine di pressione tramite operatori di proiezione per garantire che il campo di velocità rimanga privo di divergenza nell'ambito discreto. Gli autori derivano un RST proiettato che tiene conto dell'interazione tra la proiezione della pressione e la discretizzazione.

Contributi Chiave

• RST Esatto Informato dalla Discretizzazione: Il documento deriva un'espressione esatta per il tensore dello stress non risolto in una LES basata su FVM. A differenza del classico RST, questa nuova espressione include esplicitamente i contributi dei flussi numerici e degli operatori di divergenza discreti.
• Non-Simmetria e Non-Località: Per i flussi incomprimibili, gli autori dimostrano che il proposto RST LES-FVM è intrinsecamente non simmetrico e non locale (a causa della proiezione della pressione). Ciò contrasta nettamente con i classici RST della LES, che sono tipicamente assunti come simmetrici e locali.
• Errore A-Posteriori Nullo: In esperimenti "assistiti da DNS" (dove il residuo esatto è calcolato da dati di Simulazione Diretta dei Numeri e iniettato nell'equazione LES), la proposta espressione di RST produce un errore a-posteriori nullo. Al contrario, le espressioni RST classiche producono errori che si accumulano nel tempo.
• Miglioramento del Modellazione della Chiusura: Quando il nuovo RST viene utilizzato come dato target per adattare il coefficiente di un modello di Smagorinsky, il coefficiente risultante è più elevato rispetto a quando viene adattato al classico RST. Ciò porta a prestazioni migliorate in termini di spettri di energia e errore relativo, particolarmente per bassi rapporti filtro-griglia ($\Delta/h$).

Risultati
Gli autori validano la loro teoria utilizzando due casi di test: l'equazione di Burgers viscosa 1D e la turbolenza isotropa decadente 3D.

• Equazione di Burgers 1D:

  • Analisi A-Priori: Lo studio mostra che per comuni rapporti di filtro ($\Delta \le 4h$), gli errori di discretizzazione (termini di flusso e divergenza) costituiscono una parte significativa (fino al 70-80%) del residuo totale. Trascurare questi termini porta a una rappresentazione scarsa del vero residuo.
  • Analisi A-Posteriori: L'utilizzo del residuo esatto informato dalla discretizzazione come termine di chiusura produce una ricostruzione "perfetta" della soluzione DNS filtrata.
  • Modello di Smagorinsky: L'adattamento del coefficiente di Smagorinsky al nuovo RST produce un valore di coefficiente più alto rispetto al classico RST. Le simulazioni che utilizzano questo coefficiente informato mostrano errori relativi inferiori e spettri di energia che corrispondono meglio alla DNS filtrata, specialmente ai bassi rapporti $\Delta/h$.

• Equazioni di Navier-Stokes 3D:

  • La derivazione conferma che l'RST 3D è non simmetrico e non locale.
  • Esperimenti assistiti da DNS su turbolenza decadente 3D confermano che il framework proposto cattura correttamente lo stress residuo, validando l'estensione della teoria 1D al flusso completo incomprimibile.

Significato e Rivendicazioni
Il documento sostiene che la tradizionale separazione tra "modellare la fisica" (chiusura LES) e "modellare la numerica" (discretizzazione) è fallace quando la larghezza del filtro e la spaziatura della griglia sono comparabili. Unificando questi processi nel framework LES-FVM, gli autori forniscono una definizione matematicamente rigorosa dello stress non risolto che è coerente con lo schema numerico.

Il significato primario risiede nella dimostrazione che gli errori indotti dalla discretizzazione non sono meri rumori numerici, ma sono parti integrali dello stress residuo che devono essere considerate nel modello di chiusura. Gli autori sostengono che l'uso del classico RST come target per modelli di chiusura basati su dati o sintonizzati introduce un'incoerenza fondamentale, portando a coefficienti di modello subottimali ed errori accumulati. Il loro RST "informato dalla discretizzazione" fornisce il corretto target di dati, consentendo lo sviluppo di modelli di chiusura più accurati e stabili, particolarmente per la LES implicita o per scenari con griglie grossolane dove $\Delta \approx h$.

Gli autori rimangono modesti riguardo all'applicazione immediata di queste scoperte a flussi industriali complessi, osservando che, sebbene il nuovo RST migliori le prestazioni del semplice modello di Smagorinsky, il pieno potenziale di questo framework sarà probabilmente realizzato con modelli di chiusura più espressivi, non locali e non simmetrici (ad esempio, approcci basati sul machine learning) che possano sfruttare appieno la struttura del residuo derivato.