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

Technische Samenvatting: Exacte expressies voor de onopgeloste spanning in een op eindige volumes gebaseerde Large-Eddy Simulatie

Probleemstelling
Large-eddy simulatie (LES) beoogt turbulente stromingen te modelleren door grote schalen op te lossen en de effecten van onopgeloste sub-grid schalen te modelleren. Traditioneel richt de sluitingsproblematiek bij LES zich op de commutatorfout die voortvloeit uit de niet-commutativiteit van ruimtelijke filtering en niet-lineaire advectietermen, wat resulteert in de klassieke residuele spanningstensor (RST), $\tau_{ij} = \overline{u_i u_j} - \bar{u}_i \bar{u}_j$. Echter, in praktische simulaties wordt LES opgelost met numerieke discretisatieschema's, meestal de Eindige Volume Methode (FVM).

Een cruciaal probleem is dat de filterbreedte $\Delta$ in LES vaak van dezelfde orde grootte is als de roosterafstand $h$ (d.w.z. $\Delta \approx h$). In dergelijke regimes zijn discretisatiefouten — voortkomend uit numerieke fluxen en discrete divergentie-operatoren — van dezelfde orde grootte als de fysieke sub-filter spanningen. Klassieke benaderingen definiëren het sluitingsmodel doorgaans op basis van continue vergelijkingen en discretiseren deze vervolgens, of behandelen discretisatiefouten als aparte convergentieproblemen. Deze scheiding leidt tot een mismatch tussen de doelgegevens (de ware residu) en het model, met name wanneer de discretisatie zelf significante niet-lokale of niet-symmetrische effecten introduceert die standaard symmetrische, lokale sluitingsmodellen (zoals Smagorinsky) niet kunnen opvangen.

Methodologie
De auteurs stellen een verenigd wiskundig kader voor, getiteld LES-FVM, dat het LES-filteringsproces direct integreert met de FVM-discretisatie. De kern van de methodologie omvat:

1. Verenigde Formalisme: In plaats van LES-filtering toe te passen op continue vergelijkingen en deze vervolgens te discretiseren, passen de auteurs de FVM-filter (een volumegemiddelde over controlevolumes) toe op de LES-gefilterde vergelijkingen. Dit creëert een "dubbel-gefilterd" systeem waarbij de filterbreedte $\Delta$ en de roosterafstand $h$ als onafhankelijke parameters worden behandeld.
2. Afleiding van Exacte Residuen: Door gebruik te maken van een "discrete filter-swap" eigenschap (waarbij de eindige differentie-operator gelijk is aan een gefilterde exacte afgeleide), leiden de auteurs exacte expressies af voor de residuele flux in de gediscretiseerde vergelijkingen.
3. Decompositie van Residuen: De totale residuele flux $\tau_{\Delta,h}$ wordt gedecomposeerd in drie afzonderlijke componenten:
   • Klassieke Term: De standaard commutatorfout geassocieerd met de dubbele filter ($\Delta$ en $h$).
   • Flux Term: Fouten voortkomend uit het verschil tussen de exacte flux en de numerieke flux ($r$ versus $r_h$).
   • Divergentie Term: Fouten voortkomend uit het verschil tussen de exacte divergentie en de discrete divergentie.
4. Extensie naar 3D Incompressibele Stroming: Het kader wordt uitgebreid naar de 3D incompressibele Navier-Stokes vergelijkingen. Hiervoor moet de drukterm worden afgehandeld via projectie-operatoren om te waarborgen dat het snelheidveld in de discrete setting divergentievrij blijft. De auteurs leiden een geprojecteerde RST af die rekening houdt met de interactie tussen de drukprojectie en de discretisatie.

Belangrijkste Bijdragen

• Exacte Discretisatie-geïnformeerde RST: Het artikel leidt een exacte expressie af voor de onopgeloste spanningstensor in een FVM-gebaseerde LES. In tegen tegenstelling tot de klassieke RST, bevat deze nieuwe expressie expliciet bijdragen van numerieke fluxen en discrete divergentie-operatoren.
• Niet-Symmetrie en Niet-Lokaliteit: Voor incompressibele stromingen demonstreren de auteurs dat de voorgestelde LES-FVM RST inherent niet-symmetrisch en niet-lokaal is (door de drukprojectie). Dit staat scherp in contrast met klassieke LES RST's, die doorgaans als symmetrisch en lokaal worden beschouwd.
• Nul A-Posteriori Fout: In "DNS-aided" experimenten (waarbij het exacte residu wordt berekend uit Direct Numerical Simulation data en in de LES-vergelijking wordt geïnjecteerd), levert de voorgestelde RST-expressie een a-posteriori fout van nul op. In tegenstelling hiertoe resulteren klassieke RST-expressies in fouten die over de tijd accumuleren.
• Verbeterde Sluitingsmodellering: Wanneer de nieuwe RST wordt gebruikt als doelgegevens om de coëfficiënt van een Smagorinsky-model aan te passen, is de resulterende coëfficiënt hoger dan wanneer deze wordt aangepast aan de klassieke RST. Dit leidt tot verbeterde prestaties wat betreft energiespectra en relatieve fout, met name voor kleine filter-naar-rooster ratio's ($\Delta/h$).

Resultaten
De auteurs valideren hun theorie met twee testgevallen: de 1D viskeuze Burgers-vergelijking en 3D vervallende isotrope turbulentie.

• 1D Burgers Vergelijking:

  • A-Priori Analyse: De studie laat zien dat voor veelvoorkomende filterratio's ($\Delta \le 4h$), discretisatiefouten (flux- en divergentietermen) een aanzienlijk deel (tot 70-80%) van het totale residu vormen. Het negeren van deze termen leidt tot een slechte representatie van het ware residu.
  • A-Posteriori Analyse: Het gebruik van het exacte discretisatie-geïnformeerde residu als een sluitingsterm resulteert in een "perfecte" reconstructie van de gefilterde DNS-oplossing.
  • Smagorinsky Model: Het aanpassen van de Smagorinsky-coëfficiënt aan de nieuwe RST levert een hogere coëfficiëntwaarde op vergeleken met het aanpassen aan de klassieke RST. Simulaties met deze geïnformeerde coëfficiënt vertonen lagere relatieve fouten en energiespectra die beter overeenkomen met de gefilterde DNS, met name bij lage $\Delta/h$ ratio's.

• 3D Navier-Stokes Vergelijkingen:

  • De afleiding bevestigt dat de 3D RST niet-symmetrisch en niet-lokaal is.
  • DNS-aided experimenten op 3D vervallende turbulentie bevestigen dat het voorgestelde kader de residuele spanning correct vastlegt, wat de extensie van de 1D-theorie naar volledige incompressibele stroming valideert.

Betekenis en Claims
Het artikel beweert dat de traditionele scheiding tussen "het modelleren van de fysica" (LES-sluiting) en "het modelleren van de numeriek" (discretisatie) gebrekkig is wanneer de filterbreedte en de roosterafstand vergelijkbaar zijn. Door deze processen te verenigen in het LES-FVM kader, bieden de auteurs een wiskundig rigoureuze definitie van de onopgeloste spanning die consistent is met het numerieke schema.

De primaire betekenis ligt in de demonstratie dat discretisatie-geïnduceerde fouten niet louter numerieke ruis zijn, maar integrale onderdelen van de residuele spanning die in de sluitingsmodellen moeten worden meegenomen. De auteurs argumenteren dat het gebruik van de klassieke RST als doel voor data-gedreven of getunede sluitingsmodellen een fundamentele inconsistentie introduceert, wat leidt tot suboptimale sluitingscoëfficiënten en accumulerende fouten. Hun voorgestelde "discretisatie-geïnformeerde" RST biedt de correcte doelgegevens, wat de ontwikkeling van nauwkeurigere en stabielere sluitingsmodellen mogelijk maakt, met name voor implicit LES of scenario's met grove roosters waar $\Delta \approx h$.

De auteurs blijven bescheiden over de directe toepassing van deze bevindingen op complexe industriële stromingen, waarbij zij opmerken dat hoewel de nieuwe RST de prestaties van het eenvoudige Smagorinsky-model verbetert, het volledige potentieel van dit kader waarschijnlijk pas gerealiseerd zal worden met meer expressieve, niet-lokale en niet-symmetrische sluitingsmodellen (bijv. machine learning-gebaseerde benaderingen) die de structuur van de afgeleide residu volledig kunnen exploiteren.