Turbulence Modelling of Mixing Layers under Anisotropic Strain

Questo studio indaga l'impatto dei tassi di deformazione anisotropi sugli strati di miscelamento turbolenti utilizzando il modello K-L RANS, dimostrando che una chiusura della deformazione trasversale migliora l'accuratezza predittiva rispetto all'approccio isotropo predefinito e suggerendo le corrispondenti modifiche per i modelli K-$\epsilon$ e K-$\omega$.

L'essenziale

Immagina di avere una ciotola con due liquidi di colori diversi, come olio e acqua, appoggiata su un tavolo. Se scuoti il tavolo, il confine tra i due diventa disordinato e inizia a mescolarsi. Questo è simile a ciò che accade in uno strato di miscelazione turbolenta in fisica: due fluidi di densità diversa vengono spinti insieme, creando un caos vorticoso.

Questo articolo riguarda la comprensione di ciò che accade quando non si scuote solo il tavolo, ma si allunga o si comprime l'intera stanza in cui avviene la miscelazione.

Ecco una panoramica della storia dell'articolo, utilizzando semplici analogie:

1. L'ambientazione: Allungare la stanza

In molti scenari reali, come l'esplosione di una stella (supernova) o la compressione di una bomba a fusione nucleare, lo spazio in cui i fluidi si mescolano non è semplicemente fermo. Lo spazio stesso si espande o si contrae.

• L'analogia: Immagina che lo strato di miscelazione sia un pezzo di pasta da impastare. Di solito, gli scienziati studiano come la pasta si mescola quando la si spinge semplicemente. Ma in questo articolo, gli autori chiedono: "Cosa succede se, mentre impasti, qualcuno sta anche tirando il tavolo su cui si trova la pasta, allungandolo in lunghezza o comprimendolo in larghezza?"
• Il problema: Lo "stiramento" (deformazione) non è lo stesso in tutte le direzioni. Se tiri un elastico, diventa più lungo in una direzione ma più sottile nelle altre. Questo è chiamato deformazione anisotropa. La maggior parte dei modelli informatici utilizzati per prevedere queste miscele assume che lo stiramento sia uguale in ogni direzione (come gonfiare un palloncino perfetto), il che non corrisponde alla realtà.

2. Lo strumento: Il modello "K-L"

Per prevedere come si mescolano i fluidi, gli autori utilizzano un programma informatico chiamato modello di turbolenza K-L.

• L'analogia: Pensa a questo modello come a un libro di ricette per prevedere il caos. Ha due ingredienti principali che tiene traccia:
  1. Quanta energia c'è nei vortici (Energia Cinetica Turbolenta).
  2. Quanto sono grandi i vortici (Scala di Lunghezza Turbolenta).
• Il modello cerca di indovinare quanto diventeranno grandi i vortici mentre i fluidi si mescolano. La parte difficile è una regola nella ricetta chiamata termine di "compressione globale". Questa regola dice al modello come cambia la dimensione dei vortici quando l'intera stanza viene compressa o stirata.

3. L'esperimento: Testare tre regole diverse

Gli autori hanno eseguito simulazioni al computer per vedere quale "regola" per la compressione globale funzionava meglio quando la stanza veniva stirata in direzioni specifiche. Hanno testato tre versioni della ricetta:

1. La regola "Media": Assume che lo stiramento sia uguale in tutte le direzioni (l'impostazione predefinita).
2. La regola "Longitudinale": Assume che la dimensione dei vortici cambi in base solo a quanto la stanza si stirata lungo la direzione di miscelazione.
3. La regola "Trasversale": Assume che la dimensione dei vortici cambi in base a quanto la stanza si stirata attraverso la direzione di miscelazione (perpendicolare al flusso).

4. I risultati: Vince la regola "Trasversale"

Gli autori hanno confrontato le loro previsioni al computer con simulazioni ad alta risoluzione estremamente dettagliate (che fungono da riferimento "perfetto").

• La scoperta: La regola "Media" predefinita era accettabile, ma non ottima. La regola "Longitudinale" ha effettivamente reso le previsioni peggiori.
• Il vincitore: La regola "Trasversale" (utilizzando la deformazione trasversale) è stata la più accurata.
• Perché? Gli autori spiegano che quando si stirra uno strato di miscelazione, i grandi "vortici" (eddies) si comportano in modo diverso a seconda della direzione. Si scopre che la dimensione di questi vortici è più sensibile a come lo spazio cambia trasversalmente (di lato) rispetto a come cambia longitudinalmente. Utilizzando lo stiramento trasversale per regolare la dimensione dei vortici nella ricetta, il modello ha previsto la larghezza e l'energia della miscelazione con molta più accuratezza.

5. Il quadro generale: Una nuova ricetta "a tre parti"

L'articolo ha anche esaminato come semplificare queste equazioni complesse in un modello di "Galleggiamento-Resistenza" (un modo più semplice per pensare alla miscelazione).

• Hanno realizzato che la "larghezza della miscela" e la "dimensione dei vortici" reagiscono effettivamente a forze diverse. La larghezza si stirra con la trazione longitudinale, ma la dimensione dei vortici reagisce alla compressione trasversale.
• La conclusione: Per ottenere la previsione migliore, è necessario un modello che tratti queste due cose separatamente. Invece di una regola per tutto, serve un modello a tre parti che evolva la larghezza e la dimensione dei vortici in modo indipendente.

Riepilogo

In breve, questo articolo riguarda la correzione di un modello informatico utilizzato per prevedere come i fluidi si mescolano quando lo spazio che li circonda viene distorto. Gli autori hanno scoperto che il modo standard di calcolare come i "vortici" si restringono o crescono era errato per queste condizioni specifiche. Cambiando la regola per guardare come lo spazio si stirata trasversalmente invece di semplicemente mediare, hanno reso il modello molto più accurato. Questo aiuta gli scienziati a comprendere meglio eventi complessi come le esplosioni stellari o gli esperimenti sull'energia da fusione, dove i fluidi vengono costantemente compressi e stirati in modi disuguali.

Sommario tecnico

1. Enunciato del Problema

Gli strati di miscelazione turbolenta, come quelli generati dalle instabilità di Rayleigh-Taylor (RTI) e Richtmyer-Meshkov (RMI), sono fondamentali per la fisica delle alte densità di energia, inclusa la fusione a confinamento inerziale (ICF) e fenomeni astrofisici come le supernove. In questi scenari, il flusso spesso sperimenta tassi di deformazione anisotropi a causa del moto radiale in geometrie convergenti (implosioni) o divergenti (esplosioni).

I modelli di turbolenza attuali basati sulle equazioni di Reynolds-Averaged Navier-Stokes (RANS), in particolare il ampiamente utilizzato modello K-L (un modello a due equazioni che risolve per l'energia cinetica turbolenta $K$ e la scala di lunghezza turbolenta $L$), si basano tipicamente su ipotesi di chiusura derivate per deformazione isotropa o flussi incomprimibili. La chiusura standard per il termine di compressione globale nell'equazione di trasporto della scala di lunghezza assume che la scala di lunghezza evolva in base al tasso di deformazione isotropo medio (divergenza della velocità divisa per 3). Tuttavia, in scenari anisotropi realistici (ad esempio, implosioni sferiche), i tassi di deformazione nelle direzioni radiale (assiale) e circonferenziale (trasversale) differiscono significativamente. Il documento affronta la mancanza di una definizione teorica e di una modellazione accurata su come la compressione globale influisca sulla scala di lunghezza turbolenta in queste condizioni anisotrope.

2. Metodologia

Gli autori hanno condotto un'indagine sistematica utilizzando simulazioni RANS monodimensionali per isolare gli effetti della direzione della deformazione sull'evoluzione dello strato di miscelazione.

• Quadro Numerico: Le simulazioni hanno utilizzato il modello di turbolenza K-L implementato nel codice FLAMENCO. Il modello risolve le equazioni di trasporto per l'energia cinetica turbolenta mediata di Favre ($\tilde{K}$) e la scala di lunghezza turbolenta ($L$).
• Casi di Test: Lo studio si è basato sul caso di instabilità di Richtmyer-Meshkov "quarter-scale $\theta$-group" (configurazione da pesante a leggero, numero di Atwood $A_t=0.5$, numero di Mach dell'urto $M_a \approx 1.84$).
• Applicazione della Deformazione: Per isolare gli effetti dell'anisotropia, sono stati applicati tassi di deformazione normale uniformi al dominio di simulazione in due direzioni distinte:
  • Deformazione Assiale ($S_A$): Deformazione applicata nella direzione della miscelazione (normale all'interfaccia).
  • Deformazione Trasversale ($S_T$): Deformazione applicata nelle direzioni parallele all'interfaccia (piano omogeneo).
  • Sono stati testati sia profili di velocità costante che profili di tasso di deformazione costante.
• Variazioni di Chiusura: Sono state confrontate tre diverse chiusure per il termine di compressione globale ($\bar{\rho} L S_\phi$) nell'equazione della scala di lunghezza:
  1. Chiusura Isotropa (Predefinita): Scala $L$ con il tasso di deformazione medio ($S_I = \nabla \cdot \mathbf{u} / 3$).
  2. Chiusura Assiale: Scala $L$ con il tasso di deformazione assiale ($S_A$).
  3. Chiusura Trasversale: Scala $L$ con il tasso di deformazione trasversale ($S_T$).
• Validazione: I risultati RANS sono stati confrontati con dati ad alta fedeltà della Simulazione a Grande Vortice Implicita (ILES) provenienti da studi precedenti (Pascoe et al.), che hanno funto da verità di riferimento.

3. Contributi Chiave

• Identificazione delle Carenze di Chiusura: Il documento dimostra che l'ipotesi di chiusura isotropa standard non riesce a prevedere accuratamente l'evoluzione degli strati di miscelazione turbolenta sotto deformazione anisotropa, in particolare per quanto riguarda l'energia cinetica turbolenta (TKE) e la larghezza integrale.
• Proposta di Chiusura Trasversale: Gli autori propongono e validano una nuova strategia di chiusura in cui l'evoluzione della scala di lunghezza turbolenta sotto compressione globale è guidata dal tasso di deformazione trasversale ($S_T$) piuttosto che dalla deformazione isotropa media.
• Derivazione di un Modello di Galleggiamento-Resistenza a Tre Equazioni: Attraverso un'analisi di autosimilarità, gli autori hanno derivato un modello semplificato di galleggiamento-resistenza. Hanno scoperto che una singola scala di lunghezza non può catturare simultaneamente la scalatura della larghezza integrale (che segue la deformazione assiale) e la scala di lunghezza turbolenta (che segue la deformazione trasversale). Di conseguenza, hanno sviluppato un modello a tre equazioni che evolve separatamente la larghezza integrale e la scala di lunghezza turbolenta.
• Generalizzazione ad Altri Modelli: Lo studio estende i risultati ad altri popolari modelli a due equazioni (K-$\epsilon$ e K-$\omega$), derivando modifiche equivalenti per i loro termini di compressione globale basate sulla chiusura della deformazione trasversale.

4. Risultati

• Larghezza Integrale ($W$):
  • Sotto deformazione assiale, la larghezza dello strato di miscelazione viene direttamente allungata o compressa dal gradiente di velocità del flusso medio. Tutte le chiusure hanno catturato questa tendenza, ma la chiusura trasversale ha fornito la previsione più accurata del tasso di crescita turbolenta, riducendo l'Errore Percentuale Assoluto Medio (MAPE) di un fattore quattro rispetto alla chiusura isotropa.
  • Sotto deformazione trasversale, la chiusura isotropa ha mostrato quasi nessuna variazione nella larghezza (cancellando erroneamente gli effetti di produzione di taglio e compressione), mentre la chiusura trasversale ha previsto correttamente il lieve aumento del tasso di crescita osservato nei dati ILES.
• Energia Cinetica Turbolenta (TKE):
  • La chiusura trasversale ha superato significativamente gli altri metodi nella previsione della TKE.
  • La chiusura assiale ha ottenuto i risultati peggiori, prevedendo spesso il segno sbagliato del cambiamento per la TKE sotto espansione.
  • La chiusura isotropa ha sovrastimato l'influenza della deformazione sulla TKE.
• Analisi di Autosimilarità:
  • L'analisi ha rivelato che sotto deformazione anisotropa, il rapporto tra il picco della scala di lunghezza turbolenta e la larghezza integrale ($\beta_W$) non è costante, violando l'ipotesi standard di autosimilarità utilizzata in molte calibrazioni K-L.
  • Ciò ha reso necessario lo sviluppo del modello di galleggiamento-resistenza a tre equazioni, in cui la larghezza integrale scala con la deformazione assiale ($S_A$) e la scala di lunghezza turbolenta scala con la deformazione trasversale ($S_T$).
• Equivalenza dei Modelli: Le modifiche derivate per il modello K-L sono state mappate con successo sui modelli K-$\epsilon$ e K-$\omega$, suggerendo che anche il termine di compressione globale in questi modelli dovrebbe essere aggiustato per tenere conto della deformazione trasversale piuttosto che solo della divergenza media.

5. Significato

Questo lavoro fornisce un avanzamento critico nella modellazione della miscelazione turbolenta comprimibile nelle applicazioni ad alta densità di energia.

• Migliorata Capacità Predittiva: Correggendo la chiusura della compressione globale per tenere conto dell'anisotropia, il modello K-L (e i suoi equivalenti) può ora prevedere con maggiore precisione la crescita degli strati di miscelazione e la TKE in scenari di implosione/esplosione, che sono centrali per la Fusione a Confinamento Inerziale (ICF) e le simulazioni astrofisiche (ad esempio, supernove).
• Comprensione Fisica: Lo studio chiarisce i ruoli fisici distinti delle deformazioni assiali e trasversali: la deformazione assiale guida principalmente lo stiramento geometrico dello strato di miscelazione, mentre la deformazione trasversale governa l'evoluzione dei vortici turbolenti (scala di lunghezza) e la dissipazione.
• Applicazione Pratica: Le modifiche proposte sono relativamente semplici da implementare nei codici RANS esistenti, ma offrono miglioramenti sostanziali in termini di accuratezza senza i costi computazionali della LES o della DNS. Il modello di galleggiamento-resistenza a tre equazioni derivato offre uno strumento computazionalmente efficiente per applicazioni ingegneristiche dove le simulazioni RANS complete sono troppo costose.

In conclusione, il documento sostiene che per ambienti con deformazione anisotropa, la scala di lunghezza turbolenta dovrebbe essere modellata come evolutiva con il tasso di deformazione trasversale, una scoperta che sfida l'ipotesi di lunga data della scalatura isotropa nelle chiusure di turbolenza standard.