Laminar-to-Turbulent Transition of Yield-Stress Fluids in Pipe and Channel Flows

Questo studio presenta le prime simulazioni numeriche dirette (DNS) che risolvono l'intero processo di transizione da flusso laminare a turbolento per fluidi con soglia di snervamento (modello Herschel-Bulkley) in condotti e canali, rivelando come lo stress di Reynolds locale debba superare la tensione di snervamento per innescare la turbolenza e definendo i confini del regime di transizione in accordo con dati sperimentali su fluidi Carbopol.

L'essenziale

🌊 Il Mistero del "Fango" che diventa Turbina: Una Storia di Fluidi Strani

Immagina di avere due tipi di liquidi molto diversi:

1. L'acqua: Scorre sempre, si muove appena la spingi. È un fluido "normale".
2. Il Dentifricio (o la maionese): Se lo lasci fermo nel tubo, rimane lì, immobile come un solido. Devi premere forte (superare una certa "soglia di forza") per farlo uscire. Una volta uscito, però, diventa più fluido e scorre veloce. Questo è un fluido con sforzo di snervamento (in gergo tecnico: yield-stress).

Gli scienziati di questo studio (dalla Georgia Tech, negli USA) si sono chiesti: "Cosa succede quando questi fluidi 'testardi' come il dentifricio vengono spinti così forte da diventare turbolenti, proprio come l'acqua in una tempesta?"

Hanno usato supercomputer per simulare questo fenomeno, perché è molto difficile vederlo chiaramente in un esperimento reale.

🧱 Il Concetto Chiave: Il "Tappo" di Dentifricio

Quando spingi un fluido come il dentifricio in un tubo, succede una cosa strana:

• Al centro: Il liquido è così "pigro" che non si muove affatto. Si forma un tappo solido (chiamato plug) che scorre tutto insieme, come un blocco di gelato che scivola nel cono.
• Ai bordi: Solo vicino alle pareti del tubo il liquido viene "schiacciato" abbastanza da muoversi e scivolare.

Finché spingi piano, il tappo centrale rimane intatto e il flusso è calmo (laminare). Ma cosa succede se aumenti la pressione?

🚀 La Transizione: Da "Tappo" a "Tempesta"

Lo studio ha scoperto che c'è una battaglia silenziosa tra due forze:

1. La forza dell'inerzia: La spinta che vuoi dare al fluido per farlo andare veloce.
2. La forza del "tappo": La resistenza interna del fluido che vuole restare solido.

Ecco le tre fasi che hanno osservato, usando un'analogia con un traffico automobilistico:

1. Fase Laminare (Il Tappo Intatto - $Re_G < 1735$):
   Immagina un'autostrada dove tutte le auto sono bloccate in un unico grande blocco solido. Non c'è movimento relativo tra le auto. Il fluido scorre, ma al centro è un blocco unico. È tutto calmo, nessun caos.

2. Fase di Transizione (Il Blocco si Rompe - $1735 < Re_G < 2920$):
   Qui succede la magia. Spingi sempre di più. La forza diventa così forte che inizia a "rompere" il blocco centrale.

   • Il tappo di dentifricio inizia a sgretolarsi.
   • Appaiono dei vortici (piccoli tornado) vicino alle pareti, come se le auto iniziarono a fare manovre strane e a cambiare corsia.
   • È un momento di caos: a volte il tappo si riforma, a volte si rompe. È la zona di transizione dove il fluido decide se rimanere calmo o diventare turbolento.

3. Fase Turbolenta (Il Caos Totale - $Re_G > 2920$):
   Hai spinto abbastanza forte da distruggere completamente il tappo centrale. Ora il fluido è tutto mescolato, con vortici ovunque, proprio come l'acqua in una cascata o l'aria in un uragano. Il "dentifricio" si comporta finalmente come un fluido normale e turbolento.

🔍 Cosa hanno scoperto gli scienziati?

Usando simulazioni al computer (chiamate DNS, che sono come "filmati" ultra-precisi di ogni singola molecola), hanno capito tre cose fondamentali:

• La regola d'oro: La turbolenza può nascere solo se la spinta è abbastanza forte da superare la resistenza del "tappo" centrale. Se la forza non è sufficiente, il fluido rimane calmo anche se spingi forte, perché il centro resiste.
• Conferma sperimentale: Hanno confrontato i loro calcoli con esperimenti reali fatti con un gel chiamato Carbopol (usato nei dentifrici e nei gel). I risultati combaciano perfettamente: il computer ha previsto esattamente quando il fluido sarebbe diventato turbolento.
• Asimmetria: Hanno notato che, durante la fase di transizione, il fluido non si comporta sempre in modo simmetrico. A volte il caos inizia più da un lato che dall'altro, come se il traffico si bloccasse più a destra che a sinistra.

🏁 Perché è importante?

Questa ricerca non è solo teoria. Capire come questi fluidi si comportano è cruciale per:

• Industria: Trasportare fanghi, vernici, cosmetici o cibo (come il ketchup) attraverso tubi senza che si blocchino o si rovinino.
• Medicina: Capire come scorre il sangue (che ha proprietà simili) nelle arterie.
• Geologia: Prevedere come si muovono le frane o i flussi di fango durante le alluvioni.

In sintesi, questo studio ci ha dato la "mappa" per capire esattamente quando un fluido "testardo" smette di comportarsi come un solido e inizia a diventare una tempesta turbolenta, aiutandoci a progettare tubi e pompe più efficienti per il futuro.

Sommario tecnico

Ecco una sintesi tecnica dettagliata del documento presentato, redatta in italiano.

Titolo: Transizione da Laminare a Turbolento per Fluidi con Soglia di Snervamento in Flussi in Tubo e Canale

1. Problema e Contesto

Il lavoro affronta la complessa dinamica della transizione dal regime laminare a quello turbolento per i fluidi non newtoniani che seguono il modello di Herschel-Bulkley (HB). Questi fluidi, caratterizzati da una soglia di snervamento ($\tau_y$) e da un comportamento di assottigliamento al taglio (shear-thinning), sono onnipresenti in applicazioni industriali e geofisiche (es. trasporto di fanghi, prodotti cosmetici come il Carbopol, processi di stampa).

La sfida principale risiede nel fatto che la presenza di una soglia di snervamento crea regioni "non snervate" (core a plug) all'interno del flusso, dove la velocità è costante e il gradiente di velocità è nullo. Questo fenomeno altera drasticamente i meccanismi di trasporto della quantità di moto e l'insorgenza della turbolenza rispetto ai fluidi newtoniani. Nonostante studi precedenti su correlazioni empiriche e simulazioni parziali, mancava una visione unificata basata su Simulazioni Numeriche Dirette (DNS) che risolvesse l'intera sequenza di transizione (laminare, transizione, turbolento) sia per geometrie cilindriche (tubi) che piane (canali rettangolari), includendo la formazione e la rottura del "plug".

2. Metodologia

Gli autori hanno condotto simulazioni DNS ad alta fedeltà utilizzando un approccio numerico avanzato:

• Modello Reologico: È stato implementato il modello di Herschel-Bulkley regolarizzato a doppia viscosità (bi-viscosity regularization) per gestire la singolarità numerica nelle regioni a basso gradiente di velocità (dove $\dot{\gamma} \to 0$). La viscosità apparente è limitata da un valore massimo ($\nu_0$) sotto una soglia di taglio critica ($\dot{\gamma}_0$).
• Geometrie e Griglie:
  • Flusso in Tubo: Dominio circolare con raggio $R=25.4$ mm e lunghezza $L=4\pi D$. Griglia strutturata a topologia O-grid con circa $4 \times 10^6$ volumi di controllo, con clustering radiale verso la parete per risolvere i gradienti.
  • Flusso in Canale: Dominio rettangolare con lunghezza streamwise $L_x = 0.2$ m. Griglia strutturata $128 \times 129 \times 128$(circa$4.2$ milioni di celle).
• Parametri Adimensionali: È stata utilizzata una definizione di Numero di Reynolds Generalizzato ($Re_G$) basata sulla viscosità alla parete, permettendo di coprire un ampio intervallo da $Re_G = 378$ a $5300$. Sono stati inoltre calcolati i numeri di Reynolds di attrito ($Re_\tau$) e i numeri di Reynolds di Metzner-Reed ($Re_{MR}$) per confronto con la letteratura esistente.
• Validazione: I risultati sono stati confrontati con dati sperimentali (LDV) su soluzioni di Carbopol 0.1% riportati da Guzel et al. [20], focalizzandosi su profili di velocità, intensità di turbolenza e strutture vorticali.

3. Contributi Chiave

• Prima DNS Completa: Questo studio rappresenta la prima simulazione DNS che risolve l'intera sequenza di transizione (laminare $\to$ transizione $\to$ turbolento) per fluidi HB sia in tubi che in canali, fornendo una visione coerente tra le due geometrie.
• Meccanismo di Transizione: È stato identificato e quantificato il meccanismo fisico fondamentale: la transizione avviene solo quando gli sforzi di Reynolds locali superano la soglia di snervamento del fluido. Questo permette la deformazione e la successiva rottura del nucleo a "plug".
• Mappatura dei Regimi: Definizione precisa dei confini tra i regimi di flusso basati su $Re_G$, validati sia per tubi che per canali.
• Analisi Strutturale: Visualizzazione dettagliata delle strutture coerenti (vortici) tramite il criterio Q, mostrando come l'attività vorticale si evolva da strutture organizzate e sparse (laminare) a una ricca attività a piccola scala (turbolento).

4. Risultati Principali

A. Flusso in Tubo:

• Regimi Identificati:
  • $Re_G < 1735$: Regime laminare con un forte nucleo a plug e turbolenza trascurabile.
  • $1735 < Re_G < 2920$: Regime di transizione caratterizzato da un aumento netto dell'intensità di turbolenza e da un comportamento intermittente.
  • $Re_G > 2920$: Regime pienamente turbolento con il collasso del plug e dominanza della turbolenza vicino alla parete.
• Intensità di Turbolenza: L'intensità di turbolenza ($I$) mostra un picco nella regione di transizione, specialmente a $r/R \approx \pm 0.75$ (fuori dal centro), mentre il centrolinea rimane a bassa intensità finché il plug non si rompe.
• Asimmetria: È stata osservata una leggera asimmetria nei profili di turbolenza durante la transizione, coerente con le osservazioni sperimentali su fluidi con soglia di snervamento.

B. Flusso in Canale:

• Le simulazioni confermano la stessa fisica di transizione: formazione di un plug centrale, crescita delle strutture turbolente vicino alla parete (visibili tramite il criterio Q) e progressivo rafforzamento delle fluttuazioni di velocità ($u'_{rms}$) e dell'energia cinetica turbolenta all'aumentare di $Re_G$.
• I profili di velocità media passano da una forma fortemente "appiattita" (plug-dominated) a una forma più piena e tipica della turbolenza di parete.

C. Confronto Sperimentale:

• I risultati DNS mostrano un accordo eccellente con i dati sperimentali di Guzel et al. per quanto riguarda l'andamento dell'intensità di turbolenza in funzione di $Re_G$ e la forma dei profili di velocità, validando l'implementazione numerica del modello HB.

5. Significato e Implicazioni

Questo lavoro fornisce un quadro unificato per comprendere come la viscosità dipendente dal taglio e la soglia di snervamento influenzino la stabilità del flusso e la transizione alla turbolenza.

• Modellazione Industriale: I risultati offrono basi fisiche solide per migliorare i modelli di progettazione di pompe, sistemi di trasporto di fanghi e processi industriali che coinvolgono fluidi complessi.
• Fisica Fondamentale: Dimostra che la turbolenza in fluidi HB non è semplicemente una versione modificata della turbolenza newtoniana, ma è governata da un bilancio locale tra trasporto di quantità di moto turbolento e resistenza allo snervamento.
• Validazione Futura: Il lavoro prepara il terreno per ulteriori validazioni sperimentali (previste con PIV in un canale acrilico di 3.6 m) per quantificare le soglie di transizione dipendenti dalla geometria.

In sintesi, lo studio colma una lacuna significativa nella letteratura sulla fluidodinamica computazionale, offrendo una descrizione dettagliata e fisicamente coerente della transizione laminare-turbolenta per una classe di fluidi critici per l'ingegneria moderna.