Mixing of miscible liquids: Dimensionless scaling for intermediate-to-large density differences in a stirred tank

Questo studio utilizza simulazioni numeriche di un serbatoio agitato con liquidi miscibili per dimostrare che, sebbene il tempo di miscelazione correli positivamente con il numero di Richardson, una legge di scala esponenziale derivata basata sui numeri di Potenza, Froude e Richardson riesce a far convergere tutti i dati su un'unica curva maestra per differenze di densità da intermedie a grandi.

L'essenziale

Immagina di dover mescolare un enorme serbatoio contenente due liquidi diversi: uno sciroppo pesante e denso sul fondo e un succo più leggero e fluido in superficie. Fai cadere all'interno una gigantesca paletta rotante (una girante) per mescolarli insieme.

Nel mondo reale, questo è un compito comune nelle fabbriche che producono di tutto, dai farmaci al trattamento delle acque reflue. Ma ecco il punto cruciale: poiché i liquidi hanno pesi diversi (densità), quello pesante vuole rimanere sul fondo mentre quello leggero vuole galleggiare in superficie. Questo crea una "lotta" tra la paletta rotante che cerca di mescolarli e la gravità che cerca di mantenerli separati.

Questo articolo è come una storia investigativa in cui gli scienziati hanno utilizzato potenti simulazioni al computer per capire esattamente quanto tempo ci vuole per mescolare perfettamente questi due liquidi e come prevedere quel tempo senza dover costruire un serbatoio fisico ed eseguire costosi test ogni volta.

La Preparazione: Una Cucina Sperimentale Digitale

I ricercatori hanno costruito una versione virtuale di un serbatoio di miscelazione industriale standard.

• Il Serbatoio: È un grande cilindro con pareti e quattro alette verticali (paratie) per impedire al liquido di girare semplicemente in tondo come un fiume pigro.
• La Paletta: Una lama rotante al centro.
• I Liquidi: Hanno simulato un miscuglio 50/50 di un liquido pesante e uno più leggero. Non hanno utilizzato sostanze chimiche reali; li hanno semplicemente trattati come fluidi "pesanti" e "leggeri" con la stessa viscosità (spessore).
• Il Metodo: Invece di utilizzare le equazioni matematiche standard, hanno usato un trucco intelligente chiamato Metodo di Boltzmann su Reticolo. Immagina questo come simulare il liquido non come un blocco continuo, ma come miliardi di minuscole biglie da biliardo invisibili che rimbalzano e collidono. Questo ha permesso loro di vedere esattamente come si comportava la turbolenza (il vortice caotico).

La Grande Domanda: Quanto Velocemente Possiamo Mescolare?

L'obiettivo principale era trovare una "formula magica" per prevedere il tempo di miscelazione.

• Le Variabili: Hanno modificato due cose principali:
  1. La velocità di rotazione della paletta (Numero di Reynolds): Una rotazione più veloce significa solitamente più turbolenza e una miscelazione più rapida.
  2. La differenza di peso (Numero di Richardson): Se i liquidi hanno quasi lo stesso peso, si mescolano facilmente. Se uno è molto più pesante, la gravità contrasta la miscelazione, creando strati difficili da rompere.

La Scoperta: La Battaglia "Gravità contro Rotazione"

I ricercatori hanno trovato alcuni schemi interessanti:

1. Quando la Gravità Non Conta (Stesso Peso):
   Se i due liquidi hanno esattamente lo stesso peso, il tempo di miscelazione è sorprendentemente costante. Non importa quanto velocemente giri la paletta (entro un certo intervallo), il "tempo di miscelazione adimensionale" (un modo elegante per dire "quanti giri della paletta sono necessari") rimane costante a circa 20 giri. È come una regola della natura: una volta che l'acqua è sufficientemente agitata, girarla più velocemente non la fa mescolare più velocemente in termini di giri della paletta.

2. Quando la Gravità Contrattacca (Pesi Diversi):
   Non appena i liquidi hanno pesi diversi, quello pesante vuole rimanere sul fondo. Più grande è la differenza di peso, più difficile è mescolare.

   • La Tendenza: Più grande è la differenza di peso, più tempo ci vuole per mescolare.
   • La Sorprendente Svolta: Se mantieni costante la "differenza di peso" e giri semplicemente la paletta più velocemente, il tempo di miscelazione non diminuisce sempre. A volte, girare più velocemente fa effettivamente impiegarci più tempo per raggiungere un punto specifico di miscelazione.
   • Perché? Immagina che il liquido pesante sia come una coperta spessa. Se giri la paletta troppo velocemente, crei molta energia, ma il liquido pesante forma un "coperchio" stabile che il liquido più leggero non riesce a penetrare. L'energia viene sprecata nel vorticare lo strato superiore mentre lo strato inferiore rimane sigillato. È come cercare di mescolare una pentola di zuppa in cui le verdure pesanti si sono depositate in un blocco solido sul fondo; girare il cucchiaio più velocemente fa solo schizzare il brodo sopra senza rompere il blocco di verdure.

La Soluzione: Una Nuova "Curva Maestra"

Il più grande successo del team è stato creare una singola formula semplice che combina tutti questi fattori. Hanno scoperto che se si osserva il tempo di miscelazione attraverso la lente di tre numeri specifici (Potenza, Froude e Richardson), tutti i loro dati disordinati si allineano su una singola curva esponenziale e liscia.

Pensala così: prima, gli ingegneri dovevano indovinare o eseguire centinaia di test per vedere come si sarebbe mescolato un nuovo liquido. Ora, hanno una "ricetta". Se si dice loro la differenza di peso e la velocità di rotazione prevista, questa formula prevede il tempo di miscelazione con alta precisione.

Il Punto Cruciale

L'articolo conclude che per questi specifici serbatoi industriali:

• La turbolenza è fondamentale: Una volta che il liquido è completamente agitato, il comportamento di miscelazione è prevedibile.
• La gravità è il capo: Se i liquidi hanno densità diverse, la gravità crea una "stratificazione" (a strati) che resiste alla miscelazione.
• Più veloce non è sempre meglio: In sistemi con grandi differenze di densità, aumentare semplicemente la velocità del motore non garantisce una miscelazione più rapida; a volte crea solo una separazione più stabile.

Gli autori forniscono questa nuova formula per aiutare gli ingegneri a progettare processi di miscelazione migliori senza dover costruire costosi prototipi in anticipo. Hanno intenzione di testare questa formula su diverse forme di serbatoi e tipi di pale in futuro, ma per ora funziona perfettamente per i serbatoi standard che hanno simulato.

Sommario tecnico

Sintesi Tecnica: Scalatura Adimensionale per Differenze di Densità da Intermedie a Elevate in una Vasca Agitata

Enunciato del Problema
La miscelazione di liquidi miscibili è un'operazione unitaria critica nei processi industriali, che spazia dalla produzione farmaceutica al trattamento delle acque reflue. Sebbene la miscelazione turbolenta in sistemi monofase e omogenei sia ben caratterizzata, la presenza di differenze di densità tra i fluidi introduce una stratificazione guidata dalla spinta di galleggiamento. Questo fenomeno crea regioni stratificate verticalmente in cui le forze gravitazionali contrastano il trasporto convettivo, potenzialmente sopprimendo la miscelazione verticale e aumentando il tempo necessario per raggiungere l'omogeneità. Sebbene esistano correlazioni empiriche per la miscelazione di traccianti monofase e per casi stratificati specifici, un quadro generalmente valido che tenga conto coerentemente degli effetti gravitazionali, inerziali, viscosi e di galleggiamento attraverso diversi regimi di flusso e geometrie rimane sfuggente. Questo studio affronta la sfida di prevedere i tempi di miscelazione in vasche agitate con deflettori contenenti due liquidi miscibili con differenze di densità da intermedie a elevate.

Metodologia
Gli autori hanno impiegato la Dinamica dei Fluidi Computazionale (CFD) utilizzando il Metodo di Boltzmann su Reticolo (LBM) per simulare la miscelazione turbolenta in una vasca agitata con deflettori. Il sistema consisteva in un rapporto volumetrico 50/50 di due liquidi miscibili con viscosità dinamiche identiche ma densità variabili ($\rho_l$ e $\rho_h$), dove il fluido più leggero era inizialmente stratificato sopra il fluido più pesante.

• Approccio Numerico: Le simulazioni hanno risolto le equazioni di Navier-Stokes per fluidi incomprimibili utilizzando uno schema reticolare D3Q19. Per gestire alti numeri di Reynolds ($Re > 10^5$), è stato implementato un modello di Simulazione delle Grandi Scie (S-L LES) Smagorinsky-Lilly per modellare la turbolenza sottogriglia.
• Trasporto Scalare: È stato adottato un modello di trasporto scalare basato sul flusso di massa per tracciare la concentrazione del liquido più leggero. Questo modello ha accoppiato il campo scalare all'equazione della quantità di moto tramite una forza di corpo gravitazionale dipendente dalla concentrazione, consentendo la simulazione degli effetti di galleggiamento derivanti dalle variazioni locali di densità.
• Spazio di Simulazione: È stato condotto uno studio sistematico dei parametri variando la velocità di rotazione dell'agitatore ($N$) e la densità della fase più leggera ($\rho_l$). Ciò ha portato a un'ampia gamma di numeri adimensionali: numeri di Reynolds ($Re$) da 50.000 a 125.000 e numeri di Richardson ($Ri$) da 0 a 14.
• Validazione: L'accuratezza numerica è stata verificata attraverso un'analisi dell'Indice di Convergenza della Griglia (GCI), che ha stimato l'incertezza numerica a circa il 2,39% sulla griglia di produzione. Inoltre, i risultati della simulazione sono stati validati rispetto a dati sperimentali di traccianti dalla letteratura, mostrando un buon accordo nella risposta delle sonde e nelle definizioni del tempo di miscelazione.

Risultati Chiave
Lo studio ha quantificato il tempo di miscelazione ($t_m$), definito come il tempo necessario per raggiungere il 90% di omogeneità, e lo ha espresso in forma adimensionale ($t^*_m = N \cdot t_m$).

1. Baseline Omogenea ($Ri = 0$): In assenza di differenze di densità, il tempo di miscelazione adimensionale collassa su un valore costante di $t^*_m \approx 20$ su tutto l'intervallo di numeri di Reynolds investigato. Ciò conferma che il flusso è completamente turbolento e che la miscelazione inerziale domina, in accordo con la letteratura consolidata per vasche con deflettori.
2. Effetto della Stratificazione ($Ri > 0$): All'aumentare del numero di Richardson, indicando differenze di densità maggiori, il tempo di miscelazione adimensionale aumenta sistematicamente. I dati seguono una chiara tendenza esponenziale su tutti i numeri di Reynolds.
3. Legge di Scalatura Proposta: Gli autori hanno derivato una curva maestra che collassa tutti i dati di simulazione su una singola relazione di scalatura esponenziale:
   $$t^*_m = 20 \cdot e^{0.58 \cdot N_p \cdot Fr \cdot Ri}$$
   Dove $N_p$ è il numero di potenza, $Fr$ è il numero di Froude e $Ri$ è il numero di Richardson. Il prefattore di 20 rappresenta la baseline inerziale, mentre il termine esponenziale cattura l'effetto stabilizzante del galleggiamento.
4. Comportamento Non Monotono: Analizzando i tempi di miscelazione dimensionali ($t_m$) a $Ri$ costante, lo studio ha osservato una dipendenza non monotona da $Re$ per numeri di Richardson più elevati ($Ri \geq 8$). Sebbene l'aumento di $Re$ (e quindi della velocità dell'agitatore) riduca inizialmente il tempo di miscelazione potenziando la turbolenza, ulteriori aumenti portano a una stratificazione di densità più forte (a causa dell'accoppiamento dei parametri nella configurazione della simulazione), che sopprime il trasporto verticale. Questa interazione comporta un tempo di miscelazione minimo a numeri di Reynolds intermedi.
5. Struttura del Flusso: L'analisi del campo di velocità ha rivelato che ad alti $Ri$, la fase più leggera può diventare "sigillata" rispetto all'apporto di potenza dell'agitatore, con strutture verticali che appaiono principalmente nel liquido più denso. L'analisi della Densità Spettrale di Potenza (PSD) ha indicato che, sebbene il galleggiamento influenzi l'organizzazione spaziale del flusso, la dinamica guidata dall'agitatore domina ancora la cascata di energia.

Significato e Affermazioni
Il documento afferma di fornire una correlazione compatta e predittiva per i tempi di miscelazione turbolenta in sistemi inizialmente stratificati per densità all'interno del regime completamente turbolento. La legge di scalatura derivata, basata sull'influenza combinata dei numeri di Potenza, Froude e Richardson, offre una semplice espressione analitica per ingegneri ed esperti di applicazioni.

Gli autori notano che, per quanto loro sia noto, non è stata precedentemente riportata alcuna correlazione comparabile per sistemi stratificati nel regime completamente turbolento. Il significato di questo lavoro risiede nella sua capacità di snellire la progettazione dei processi di miscelazione riducendo la dipendenza da costosi test nel mondo reale. Il modello proposto consente la stima dei tempi di miscelazione in sistemi in cui gli effetti di galleggiamento sono significativi, catturando la complessa competizione tra turbolenza potenziata e stratificazione rafforzata. Gli autori suggeriscono modestamente che il lavoro futuro coinvolgerà il test di questa correlazione su diverse geometrie di vasca e tipi di agitatori, inclusi agitatori angolari, per valutare la robustezza della scalatura proposta.