Leveraging unstructured grids for direct numerical simulations of wall turbulence

Questo articolo introduce la {\eta}-griglia, un framework di generazione di griglie non strutturate per simulazioni numeriche dirette della turbolenza di parete che scala le dimensioni della griglia con la scala locale di Kolmogorov, ottenendo un'accuratezza paragonabile a quella delle griglie cartesiane convenzionali riducendo al contempo in modo significativo il costo computazionale, specialmente ad alti numeri di Reynolds e su geometrie complesse a riblet.

L'essenziale

Immagina di provare a simulare come l'aria scorre sopra un'auto o come l'acqua si muove lungo lo scafo di una nave. Per farlo con precisione su un computer, gli scienziati utilizzano una tecnica chiamata Simulazione Numerica Diretta (DNS). Pensa alla DNS come alla creazione di un gigantesco microscopio digitale tridimensionale che scompone il fluido (aria o acqua) in milioni di minuscoli cubi invisibili (una griglia). Il computer calcola quindi come ogni singolo cubo si muove e interagisce con i suoi vicini.

Il problema è che il flusso del fluido vicino a una superficie (come il fianco di una nave) è incredibilmente caotico e dettagliato. Per ottenere un quadro chiaro, è necessario un numero massiccio di questi minuscoli cubi proprio accanto alla superficie. Tuttavia, man mano che ci si allontana dalla superficie, il caos si attenua e non sono più necessari cubi così piccoli.

Il Vecchio Metodo: Il "Muro di Mattoni Rigidi"

Tradizionalmente, gli scienziati utilizzavano una griglia cartesiana. Immagina di costruire un muro con mattoni identici e rigidi.

• Il Problema: Per vedere i minuscoli dettagli vicino alla superficie, devi rendere i mattoni alla base molto piccoli. Ma poiché questi mattoni sono rigidi e collegati in linea retta, sei costretto a utilizzare quegli stessi mattoni minuscoli fino alla cima del tuo muro, anche dove i dettagli non sono importanti.
• Il Risultato: Finisci con un muro composto da miliardi di mattoni minuscoli, la maggior parte dei quali è inutile. Questo rende la simulazione al computer incredibilmente lenta e costosa, come cercare di contare ogni singolo granello di sabbia su una spiaggia solo per misurare la marea.

La Nuova Soluzione: La "Rete Intelligente ed Elastica"

Questo articolo introduce un nuovo metodo chiamato griglia $\eta$ (griglia eta). Invece di mattoni rigidi, immagina una rete da pesca intelligente ed elastica.

• Come funziona: Gli autori hanno progettato un sistema in cui la dimensione dei fori della rete cambia automaticamente in base alla quantità di dettaglio necessaria.
  • Vicino alla superficie (Il "Livello Interno"): La rete ha fori molto piccoli e stretti per catturare i minuscoli vortici caotici del fluido.
  • Più lontano (Il "Livello Esterno"): Man mano che il fluido si calma, la rete si allunga automaticamente, rendendo i fori molto più grandi.
• L'Ingrediente Segreto: La dimensione di questi fori si basa su qualcosa chiamato scala di Kolmogorov (indicata come $\eta$). Pensa a questa come al "vortice più piccolo possibile" che può esistere nel fluido a una data altezza. La nuova griglia dice semplicemente: "Rendi la dimensione del foro appena abbastanza grande da catturare il vortice più piccolo a questa specifica altezza, e non di più."

Perché è una Grande Notizia

Gli autori hanno testato questa "rete intelligente" su due diversi tipi di codici informatici (uno simile a un metodo agli elementi spettrali, l'altro simile a un metodo ai volumi finiti) e l'hanno confrontata con il vecchio metodo dei "mattoni rigidi".

1. Precisione: I risultati erano quasi identici. La "rete intelligente" ha catturato la fisica esattamente quanto i "mattoni rigidi", con una differenza inferiore all'1% nelle misurazioni chiave come attrito e velocità.
2. Risparmi Massicci: È qui che avviene la magia.
   • Per superfici lisce (come una parete piana), la nuova griglia ha ridotto il numero di "mattoni" richiesti (punti della griglia) di circa il 90% ad alte velocità.
   • Per superfici ruvide (come una parete con minuscole scanalature chiamate "riblet" progettate per ridurre la resistenza), i risparmi sono stati ancora più drammatici: fino al 97% in meno di punti della griglia.

L'Analogia del "Muro Scanalato"

Per comprendere la parte relativa ai riblet, immagina un muro coperto da minuscole scanalature parallele (come la texture di una palla da golf o la pelle di uno squalo).

• Il Vecchio Metodo: Per simulare questo, il metodo dei mattoni rigidi doveva mantenere i mattoni minuscoli ovunque perché le scanalature costringevano la griglia a essere fine fino in alto. Era come cercare di contare ogni singolo filo di un maglione, anche le parti lontane dal tessuto.
• Il Nuovo Metodo: La "rete intelligente" sa che una volta che ci si trova a pochi centimetri sopra quelle minuscole scanalature, il flusso diventa di nuovo liscio. Allarga immediatamente i fori appena sopra le scanalature, ignorando i dettagli minuscoli che non contano più.

Il Conclusione

Gli autori hanno creato un framework che agisce come una lente di zoom intelligente per le simulazioni di fluidi. Concentra la sua potenza di calcolo esattamente dove serve (vicino alla parete) e si rilassa dove non serve.

• Per Pareti Lisce: Scala lo sforzo molto più lentamente man mano che la simulazione diventa più grande.
• Per Pareti Ruvide: Scala ancora meglio, rendendo fattibili su computer che precedentemente non potevano gestirle le simulazioni di superfici complesse che riducono la resistenza.

In sintesi, hanno trovato un modo per ottenere lo stesso lavoro di alta qualità con una frazione della potenza di calcolo, trasformando un compito che potrebbe richiedere a un supercomputer un mese in uno che potrebbe essere completato in pochi giorni.

Sommario tecnico

Sintesi Tecnica: Sfruttamento di Griglie Non Strutturate per Simulazioni Numeriche Dirette della Turbolenza Parete

Enunciato del Problema
La Simulazione Numerica Diretta (DNS) della turbolenza parete è computazionalmente costosa, principalmente a causa della necessità di risolvere la scala di lunghezza di Kolmogorov ($\eta$) in tutto il dominio. I solver convenzionali basati su griglie cartesiane, ampiamente utilizzati per superfici lisce e non lisce (ad esempio tramite Metodi di Frontiera Immersa), impongono tipicamente una risoluzione della griglia streamwise ($\Delta x^+$) e spanwise ($\Delta z^+$) fissa o lentamente variabile in tutto il dominio. Questo approccio porta a una sovrarisoluzione significativa nelle regioni esterne, dove la scala di Kolmogorov cresce, creando un collo di bottiglia computazionale. Ciò è particolarmente acuto per le simulazioni che coinvolgono superfici complesse come i riblet, dove la griglia deve essere estremamente fine vicino alla superficie per risolvere i micro-solchi, eppure questa risoluzione fine viene mantenuta inutilmente per tutta lo spessore dello strato limite nei framework cartesiani.

Metodologia
Gli autori propongono un nuovo framework di generazione di griglie non strutturate denominato $\eta$-grid. Il principio fondamentale consiste nel impostare le dimensioni della griglia nella direzione normale alla parete ($\Delta y^+$) e spanwise ($\Delta z^+$) proporzionali alla scala locale di Kolmogorov ($\eta^+$), consentendo un alleggerimento della griglia nelle regioni del flusso esterno dove le scale turbolente sono più grandi.

1. Formulazione della Griglia:

   • Strato Interno: Uno strato sottile ($\sim 50$ unità viscose) adiacente alla parete utilizza dimensioni della griglia scalate viscosamente simili alla DNS convenzionale ($0.3 \lesssim \Delta y^+ \lesssim 4$, $\Delta z^+ \simeq 5$). Per superfici non lisce (ad esempio riblet), la griglia risolve la lunghezza d'onda superficiale più piccola ($\ell^+$) con $\ell^+/30 \lesssim \Delta y^+, \Delta z^+ \lesssim 4$.
   • Regione Esterna: Al di sopra dello strato interno, le dimensioni della griglia aumentano proporzionalmente alla scala locale di Kolmogorov ($\Delta y^+ \simeq \Delta z^+ \simeq 2\eta^+$).
   • Leggi di Scalatura: Gli autori affinano adattamenti semi-empirici per $\eta^+$ sia nella regione logaritmica che nella regione di scia per il flusso turbolento in canale e per strati limite turbolenti (TBL) a gradiente di pressione zero (ZPG). Propongono adattamenti a legge di potenza a tratti per $\eta^+$ che transizionano a $y^+ \approx \delta^+/2$.
   • Adattamento ai Riblet: Per le geometrie a riblet, il framework introduce un "sottostrato a riblet" ($\delta_r^+$) riempito di celle quadrate per risolvere i flussi secondari e i rotori di Kelvin-Helmholtz, transizionando in modo fluido verso la scalatura basata su $\eta$ al di sopra.

2. Implementazione:

   • Il framework utilizza il pacchetto open-source Gmsh per generare griglie non strutturate multi-blocco.
   • Le griglie sono testate utilizzando due solver distinti: SOD2D (Metodo agli Elementi Spettrali) e OpenFOAM (Metodo ai Volumi Finiti).
   • Le simulazioni coprono il flusso turbolento in canale e TBL ZPG su pareti lisce e varie geometrie a riblet (solchi triangolari) fino a un numero di Reynolds di attrito $\delta_0^+ = 1000$.

Contributi Chiave

• Sviluppo del Framework: Una strategia robusta di generazione di griglie non strutturate multi-blocco che adatta dinamicamente la spaziatura della griglia alla scala locale di Kolmogorov, applicabile sia a superfici lisce che a superfici ruvide complesse.
• Leggi di Scalatura Affinate: Adattamenti semi-empirici aggiornati per la scala di Kolmogorov nella regione di scia dei TBL e dei flussi in canale, migliorando l'accuratezza del dimensionamento della griglia nello strato esterno.
• Indipendenza dal Solver: Dimostrazione che il framework $\eta$-grid produce risultati coerenti e ad alta fedeltà attraverso diversi schemi di discretizzazione numerica (SEM e FVM).
• Formulazione Specifica per Riblet: Una strategia di griglia su misura per i riblet che risolve la fisica del sottostrato (flussi secondari, rotori KH) consentendo al contempo un rapido alleggerimento nel flusso esterno, superando i limiti degli approcci Cartesiani/IBM.

Risultati

• Accuratezza: La $\eta$-grid produce risultati con una differenza inferiore all'1% rispetto alle griglie cartesiane convenzionali e ai dati di riferimento DNS/sperimentali. Le metriche includono i coefficienti di attrito superficiale ($C_f$), i profili di velocità media ($U^+$), gli stress turbolenti e gli spettrogrammi.
• Efficienza della Griglia:
  • Pareti Lisce: Il numero di punti della griglia ($N_\eta$) scala come $\propto \delta_0^{+2.5}$, rispetto a $\propto \delta_0^{+3.0}$ per le griglie cartesiane con stiramento tangente iperbolico. A $\delta_0^+ = 6000$, la $\eta$-grid richiede solo $\sim 10\%$ dei punti di griglia di una griglia cartesiana standard ($N_\eta/N_{Tanh} \simeq 0.1$).
  • Riblet: Il guadagno di efficienza è ancora più pronunciato. Per i riblet triangolari a riduzione di resistenza, $N_\eta$ scala come $\propto \delta_0^{+2.0}$, mentre le griglie cartesiane scalano come $\propto \delta_0^{+3.0}$. A $\delta_0^+ = 6000$, la $\eta$-grid richiede solo $\sim 3\%$ dei punti di griglia di una griglia cartesiana ($N_\eta/N_{Tanh} \simeq 0.03$).
• Fisica dei Riblet: Il framework cattura con successo meccanismi di flusso complessi sopra i riblet, inclusi l'origine virtuale, i flussi secondari e i rotori di Kelvin-Helmholtz, con convergenza della griglia ottenuta utilizzando parametri di griglia moderati ($C_y=2.0, C_z=2.5, y_{in}^+=50$).
• Validazione TBL: La DNS di TBL ZPG su riblet ha mostrato un buon accordo con i dati sperimentali (Baron & Quadrio; Choi & Orchard) riguardo alle percentuali di riduzione della resistenza e ai profili di velocità, nonostante le discrepanze intrinseche tra sperimentale e DNS nelle regioni di scia.

Significato
Il documento afferma che il framework $\eta$-grid offre un "enorme risparmio di griglia", rendendo computazionalmente fattibile la DNS ad alto numero di Reynolds della turbolenza parete su superfici complesse. Allineando la risoluzione della griglia alle scale fisiche locali (scala di Kolmogorov) invece di mantenere risoluzioni fisse, il metodo riduce drasticamente il costo computazionale (punti di griglia e vincoli sul passo temporale) senza sacrificare l'accuratezza. Gli autori sottolineano che questo approccio rende economicamente sostenibile condurre DNS di TBL su superfici complesse a numeri di Reynolds che corrispondono alle condizioni sperimentali, un compito precedentemente ostacolato dal costo proibitivo dei metodi basati su griglie cartesiane. Il lavoro convalida l'efficacia delle griglie non strutturate per simulazioni di turbolenza ad alta fedeltà, superando i vincoli delle mesh cartesiane strutturate.