Two pathways to diapycnal mixing in strongly stratified flows with no initial vertical shear

Questo studio combina la teoria lineare e le simulazioni numeriche dirette per rivelare che nei flussi fortemente stratificati privi di sforzo di taglio verticale iniziale, le instabilità di taglio orizzontale guidano inevitabilmente il rimescolamento diapicnale attraverso due percorsi distinti — o tramite l'emergere diretto di uno sforzo di taglio verticale o attraverso un'evoluzione non lineare in vortici colonnari — entrambi i quali innescano infine piccole instabilità di Kelvin-Helmholtz ma producono diverse efficienze di rimescolamento a causa della loro eccitazione di diverse scale verticali.

L'essenziale

Immaginate l'oceano o l'atmosfera come una gigantesca torta a strati. Gli strati sono fatti di fluidi con densità diverse (come diversi gusti di torta) e non amano mescolarsi facilmente. Di solito, gli scienziati studiano cosa succede quando si spingono questi strati lateralmente e verso l'alto o verso il basso contemporaneamente (taglio verticale). Ma questo articolo pone una domanda diversa: cosa succede se si spingono gli strati solo lateralmente, senza un movimento iniziale su e giù, in un ambiente fortemente stratificato?

I ricercatori hanno scoperto che, anche se si parte con un flusso perfettamente piatto e orizzontale, la natura possiede due "ricette" o percorsi distinti per creare il caos (turbolenza) e mescolare gli strati. La scelta della ricetta da parte della natura dipende interamente da come si "inocula" l'esperimento — essenzialmente, dal piccolo colpetto che si dà al fluido all'inizio.

Ecco la suddivisione dei due percorsi utilizzando analogie semplici:

L'Incipit: Il Fiume Calmo

Immaginate un fiume ampio e calmo che scorre orizzontalmente. L'acqua è stratificata come una pila di pancake (forte stratificazione). All'inizio, il flusso è fluido e bidimensionale (si muove solo a destra e a sinistra, non su e giù).

Percorso 1: L'Effetto "Stanza Affollata" (La Via Diretta)

Come inizia: Date al fiume un piccolo, casuale colpetto ovunque nello stesso momento (come lanciare un pugno di coriandoli in aria).
Cosa succede:

1. L'Ondulazione: A causa degli strati, il fluido non si limita a incresparsi a destra e sinistra; inizia immediatamente a incresparsi su e giù in molte dimensioni diverse contemporaneamente. Pensate a una folla di persone in una stanza che cercano di muoversi tutte insieme, creando un movimento caotico e multidirezionale.
2. Il Taglio: Queste increspature creano forti correnti verticali (taglio) molto rapidamente.
3. La Rottura: Queste correnti verticali diventano così forti da spezzarsi, creando piccoli vortici violenti (come minuscoli mulinelli). Questa è l'instabilità di "Kelvin-Helmholtz", che assomiglia alle onde che si infrangono quando il vento soffia sull'acqua.
   Il Risultato: Il mescolamento avviene in modo efficiente. Poiché l'energia è distribuita su molte dimensioni diverse di increspature, l' "attrito" (dissipazione viscosa) è minore, rendendo il processo di mescolamento relativamente efficiente.

Percorso 2: L'Effetto "Danza Sincronizzata" (La Via Indiretta)

Come inizia: Date al fiume una spinta molto specifica e organizzata (come un direttore d'orchestra che agita la bacchetta per far muovere tutti secondo uno schema preciso).
Cosa succede:

1. Il Vortice: Invece di increspature caotiche, il fluido si organizza in lunghe colonne verticali di acqua rotante (come enormi tornado che stanno in piedi nel fiume). Per molto tempo, il flusso rimane perfettamente bidimensionale, composto solo da queste grandi colonne rotanti.
2. L'Oscillazione: Alla fine, queste enormi colonne diventano instabili. Iniziano a oscillare in un modo molto specifico e ad alta frequenza. I ricercatori chiamano questo fenomeno "instabilità iperbolica". Immaginate una trottola che inizia a oscillare violentemente proprio prima di cadere.
3. La Rottura: Questa oscillazione violenta crea strati molto sottili e netti di taglio verticale. Questi strati sottili poi si spezzano in piccoli vortici violenti, proprio come nel Percorso 1.
   Il Risultato: Il mescolamento avviene, ma è meno efficiente. Perché? Perché questo percorso crea strati estremamente sottili e netti. Richiede molta energia (attrito) per creare e rompere questi strati così piccoli e sottili. È come cercare di tagliare un grosso blocco di formaggio con un coltello smussato (Percorso 1) rispetto a una lametta (Percorso 2); la lametta crea un taglio molto più netto e dispendioso in termini di energia.

La Grande Conclusione

L'articolo dimostra che il taglio verticale (movimento su e giù) non deve essere presente all'inizio. È un prodotto inevitabile del taglio orizzontale in fluidi fortemente stratificati, a patto che il fluido sia abbastanza denso (alto numero di Reynolds).

• Se si parte con rumore casuale: Si ottiene il Percorso 1 (Mescolamento diretto ed efficiente).
• Se si parte con un modello specifico: Si ottiene il Percorso 2 (Percorso indiretto, meno efficiente).

I ricercatori hanno utilizzato potenti simulazioni al computer per dimostrare che questi due percorsi sono reali, distinti e che la "ricetta" scelta all'inizio determina quanta energia viene sprecata sotto forma di calore e quanta viene invece utilizzata per mescolare effettivamente gli strati.

In breve: Anche in un fluido stratificato perfettamente calmo, una spinta orizzontale creerà inevitabilmente un caos verticale. Ma a seconda di come si dà la spinta, quel caos sarà diverso e mescolerà gli strati con diversi livelli di efficienza.

Sommario tecnico

Problematica
La stabilità dei flussi orizzontali su larga scala e fortemente stratificati rispetto a perturbazioni tridimensionali (3D) è un meccanismo critico per il trasferimento verso il basso dell'energia e il trasporto diapicnale nei sistemi geofisici e astrofisici. Mentre i flussi stratificati con taglio verticale sono ampiamente studiati, i loro controparti con taglio orizzontale (senza un iniziale taglio verticale) hanno ricevuto meno attenzione. Una questione centrale rimane: come fanno i flussi di taglio orizzontale inizialmente bidimensionali (2D) e fortemente stratificati ($Fr \ll 1$) a transire verso la turbolenza e generare uno shear verticale? La letteratura precedente suggerisce due percorsi distinti: uno che coinvolge la crescita diretta di modi primari 3D che portano a instabilità di Kelvin-Helmholtz (KH), e un altro che coinvolge l'evoluzione non lineare in vortici colonnari seguiti da instabilità secondarie. Tuttavia, non è stato chiaro se questi percorsi siano mutuamente esclusivi a causa di differenze nelle condizioni iniziali (ad esempio, tangente iperbolica rispetto a profili sinusoidali) o nella formulazione del problema (Problema di Valore Iniziale vs. forzatura corporea), o se rappresentino rotte alternative disponibili all'interno dello stesso sistema fisico. Inoltre, l'impatto del percorso selezionato sull'efficienza del mixing diapicnale non è stato sistematicamente quantificato.

Metodologia
Gli autori combinano la teoria della stabilità lineare con Simulazioni Numeriche Dirette (DNS) per investigare queste transizioni in un dominio triplamente periodico. Lo studio si concentra su un flusso di Kolmogorov piano, orizzontale e verticalmente invariante ($\bar{u}(y) = \sin(y)\mathbf{e}_x$) in un fluido linearmente stratificato con frequenza di galleggiamento costante $N$. Le simulazioni sono condotte a numeri di Reynolds ($Re = 10.000$) e di Péclet ($Pe = 10.000$) elevati, con numeri di Froude bassi ($Fr = 0,15$e$Fr = 0,1$), che risultano in numeri di Reynolds di galleggiamento ($Re_b = Fr^2 Re$) pari a 225 e 100, rispettivamente.

Per isolare l'effetto delle condizioni iniziali sulla selezione del percorso, vengono eseguiti due set di DNS:

1. DNS1 (Percorso P1): Inizializzato con rumore bianco a scala di griglia, permettendo a tutti i modi instabili della instabilità di taglio orizzontale primaria di crescere.
2. DNS2 (Percorso P2): Inizializzato con rumore a scala di griglia più un "seme" specifico per il modo eigen più veloce del 2D ($k_z=0$) della instabilità primaria.

Inoltre, gli autori eseguono un'analisi di stabilità lineare di un modello analitico "congelato" del flusso colonnare che emerge in DNS2. Questo modello approssima il flusso 2D evoluto non linearmente come un array di vortici di Taylor-Green sovrapposto a flussi sinusoidali, consentendo il calcolo degli eigenmodi e dei tassi di crescita delle perturbazioni 3D secondarie.

Contributi Chiave e Risultati

1. Conferma di Due Percorsi in un Singolo Setup:
   Lo studio dimostola che entrambi i percorsi verso la turbolenza possono verificarsi all'interno dello stesso setup di modello, distinti esclusivamente dallo spettro delle perturbazioni iniziali:

• Percorso P1 (Crescita 3D Diretta): Quando inizializzato con rumore bianco, gli eigenmodi con $k_z \neq 0$ della instabilità di taglio orizzontale primaria crescono direttamente. Questi modi generano uno shear verticale su una vasta gamma di scale, che diventa successivamente instabile a instabilità KH secondarie su piccola scala quando $Re_b$ è sufficientemente grande. Questo percorso si allinea alle osservazioni nei flussi di Kolmogorov forzati da un corpo esterno.
• Percorso P2 (Evoluzione Colonnare): Quando viene seminato il modo $k_z=0$, esso domina il flusso, guidando lo sfondo verso un flusso colonnare (vorticoso) 2D a lungo termine e tempo-dipendente. Questo stato 2D è successivamente instabile a modi 3D secondari. Gli autori identificano questi modi secondari come instabilità iperboliche derivanti dai punti di stasi iperbolici nell'array di vortici. Queste instabilità presentano alti numeri d'onda verticali ($k_z \sim Fr^{-1}$) e guidano l'emergere dello shear verticale, che alla fine innesca instabilità KH terziarie.

2. Caratterizzazione dell'Instabilità Secondaria nel Percorso P2:
   Attraverso l'analisi di stabilità lineare del flusso colonnare congelato, gli autori mostrano che l'instabilità secondaria è un'instabilità iperbolica ad alto $k_z$. Il tasso di crescita di questa instabilità scala con il tasso di deformazione locale dei punti iperbolici ($\Delta_h$) ed è approssimativamente $\lambda \approx 0,46 \Delta_h$, ampiamente indipendente dalla forza della stratificazione per ampiezze elevate. L'instabilità si manifesta come "increspature" nel campo di velocità verticale centrate attorno ai punti iperbolici, distinte dalla tradizionale instabilità "zigzag" osservata in coppie di vortici isolate.

3. Impatto sull'Efficienza di Mixing:
   Lo studio quantifica il mixing diapicnale e l'efficienza di mixing ($\eta$) per entrambi i percorsi.

• Dissipazione del Galleggiamento ($\chi$): Il tasso massimo di dissipazione della varianza del galleggiamento è paragonabile in entrambi i percorsi, poiché è controllato da moti a piccola scala ($k_z \gtrsim Fr^{-1}$) in entrambi i casi.
• Dissipazione Viscosa ($\epsilon$): Il Percorso P2 presenta una dissipazione viscosa significativamente più alta rispetto al Percorso P1. Ciò accade perché l'instabilità iperbolica in P2 eccita uno spettro stretto ad alto $k_z$, portando a una dissipazione intensa in sottili strati di shear verticale. Al contrario, P1 eccita uno spettro ampio di scale verticali, distribuendo la dissipazione in modo più esteso.
• Efficienza di Mixing: Di conseguenza, l'efficienza di picco del mixing è significativamente più alta in P1 ($\eta_{max} \approx 0,25$) rispetto a P2 ($\eta_{max} \approx 0,165$). La differenza deriva dal fatto che P2 dissipa più energia cinetica rispetto alla quantità di galleggiamento miscelato.

Significatività e Rivendicazioni
L'articolo sostiene che l'emergere dello shear verticale che guida le instabilità KH a piccola scala è un sottoprodotto inevitabile delle instabilità di taglio orizzontale in flussi fortemente stratificati a un $Re_b$ sufficientemente grande, a condizione che il flusso possieda un punto di inflessione. La selezione tra i due percorsi è determinata dalle condizioni iniziali: specificamente, il rapporto tra l'ampiezza iniziale del modo $k_z=0$ rispetto ai modi $k_z \neq 0$.

Gli autori suggeriscono che la prevalenza del Percorso P2 nella letteratura precedente (specificamente negli studi sugli strati di taglio a tangente iperbolica) possa essere un artefatto delle pratiche di inizializzazione in cui un seme $k_z=0$ viene spesso aggiunto per accelerare i calcoli, sopprimendo così la crescita diretta dei modi 3D (Percorso P1). Al contrario, la mancanza di P1 in alcuni contesti può essere dovuta a limiti di stratificazione più deboli dove l'intervallo di modi $k_z \neq 0$ instabili è limitato.

In definitiva, il lavoro stabilisce che le instabilità di taglio orizzontale sono attori chiave nella generazione della turbolenza stratificata e del mixing diapicnale, anche in assenza di un preesistente shear verticale. Evidenzia come la specifica rotta verso la turbolenza (crescita 3D diretta vs. transizione 2D-a-3D) alteri fondamentalmente lo spettro di dissipazione dell'energia e, di conseguenza, l'efficienza del mixing, un fattore che deve essere considerato quando si parametrizzano la turbolenza nei modelli climatici e geofisici. Gli autori notano che, sebbene i loro risultati siano robusti per flussi idealizzati, l'interazione con la rotazione e i campi d'onda ambientali in contesti geofisici realistici rimane una questione aperta.