Miscelazione in fluidi viscoelastici utilizzando la turbolenza elastica

Informazioni di base

• ID articolo: 2409.06391
• Titolo: Mixing in viscoelastic fluids using elastic turbulence
• Autori: Reinier van Buel, Holger Stark (Technische Universität Berlin)
• Classificazione: physics.flu-dyn
• Data di pubblicazione: Settembre 2024
• Link articolo: https://arxiv.org/abs/2409.06391

Riassunto

Questo articolo esamina l'effetto della turbolenza elastica sulla miscelazione del campo di concentrazione scalare in fluidi viscoelastici all'interno di una geometria Taylor-Couette bidimensionale, utilizzando soluzioni numeriche del modello Oldroyd-B. Lo stato del flusso è determinato mediante un parametro d'ordine del flusso secondario, che indica la regione di turbolenza elastica. Quando si parte da uno stato turbolento e successivamente si riduce il numero di Weissenberg, al di sotto di Wic appare un flusso debolmente caotico. L'advection sia nello stato turbolento che in quello debolmente caotico induce miscelazione, illustrata attraverso l'evoluzione temporale della deviazione standard della concentrazione di soluto rispetto alla distribuzione uniforme. In particolare, la miscelazione è forte nello stato di turbolenza elastica, quantificata mediante il tasso di miscelazione, il tempo di miscelazione e l'efficienza di miscelazione. Tutte e tre le quantità seguono leggi di scala. Significativamente, il parametro d'ordine è fortemente correlato al tasso di miscelazione, e quindi rappresenta anche un buon indicatore della miscelazione all'interno del fluido.

Contesto di ricerca e motivazione

Contesto del problema

1. Sfide di miscelazione a basso numero di Reynolds: Nei fluidi newtoniani, a basso numero di Reynolds il flusso è laminare e il processo di miscelazione è lento, controllato principalmente dalla diffusione. Questo rappresenta un importante limite nei dispositivi microfluidici.
2. Scoperta della turbolenza elastica: L'aggiunta di polimeri ad alto peso molecolare al fluido lo rende viscoelastico, consentendo l'esibizione di fenomeni di turbolenza elastica anche quando l'inerzia è trascurabile.
3. Importanza della miscelazione a microscala: Nei dispositivi lab-on-a-chip, una miscelazione fluida efficiente è cruciale, in particolare a scala micrometrica.

Motivazione della ricerca

• Esplorare il potenziale della turbolenza elastica nel migliorare l'efficienza di miscelazione
• Comprendere la natura fisica dei meccanismi di miscelazione nei fluidi viscoelastici
• Fornire guida teorica per applicazioni microfluidiche

Contributi principali

1. Identificazione della natura della transizione: Identificazione della transizione da flusso laminare a turbolenza elastica come transizione subcritica, con scoperta dell'esistenza di uno stato di flusso debolmente caotico
2. Stabilimento di metodi di caratterizzazione della miscelazione: Introduzione di tre indicatori quantitativi—tasso di miscelazione, tempo di miscelazione ed efficienza di miscelazione—che seguono tutti leggi di scala
3. Scoperta della forte correlazione tra parametro d'ordine e miscelazione: Dimostrazione della forte correlazione tra il parametro d'ordine del flusso secondario e il tasso di miscelazione, utilizzabile come buon indicatore del grado di miscelazione del fluido
4. Verifica del meccanismo di miscelazione di Batchelor: Osservazione di decadimento esponenziale nello stato di turbolenza elastica, corrispondente alla regione di miscelazione di Batchelor, in accordo con i risultati sperimentali

Dettagli metodologici

Configurazione geometrica e equazioni di controllo

Lo studio impiega una geometria Taylor-Couette bidimensionale contenente due cilindri concentrici:

• Raggio cilindro esterno: $r_o = 10\mu m$
• Raggio cilindro interno: $r_i = 2.5\mu m$
• Il cilindro esterno ruota in senso antiorario con velocità angolare $\Omega = 2\pi s^{-1}$

Il moto del fluido è descritto dall'equazione di Navier-Stokes generalizzata: $\rho\left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u}\right) = -\nabla p + \eta_s \nabla^2 \mathbf{u} + \nabla \cdot \boldsymbol{\tau}$

Modello costitutivo Oldroyd-B

Il tensore di stress polimerico $\boldsymbol{\tau}$ è modellato secondo il modello Oldroyd-B: $\boldsymbol{\tau} + \lambda \stackrel{\nabla}{\boldsymbol{\tau}} = \eta_p[\nabla \otimes \mathbf{u} + (\nabla \otimes \mathbf{u})^T]$

dove $\stackrel{\nabla}{\boldsymbol{\tau}}$ è la derivata convettiva superiore.

Modellazione della miscelazione

L'evoluzione del campo di concentrazione scalare $c(\mathbf{r}, t)$ è governata dall'equazione di advection-diffusione: $\frac{\partial c}{\partial t} + (\mathbf{u} \cdot \nabla)c = D\nabla^2 c$

In questo studio, $D = 0$, trascurando la diffusione molecolare.

Parametri adimensionali chiave

• Numero di Reynolds: $Re = \rho\Omega r_o^2/\eta_s = 0.628 \times 10^{-3}$
• Numero di Weissenberg: $Wi = \lambda\Omega$
• Rapporto di viscosità: $\beta = \eta_p/\eta_s = 1.5$

Configurazione sperimentale

Metodo numerico

• Utilizzo del risolutore open-source OpenFOAM® a volume finito
• Impiego del risolutore specializzato RheoTool per il flusso viscoelastico
• Griglia: mesh radiale con 100 elementi radiali e 120 elementi angolari
• Passo temporale: $\delta t = 10^{-5}s$

Definizione del parametro d'ordine

Definizione dell'intensità del flusso secondario come parametro d'ordine: $\Phi = \sigma = \sqrt{\langle[\mathbf{u}(\mathbf{r},t) - \mathbf{u}_0(\mathbf{r})]^2\rangle}/u_0^{\max}$

Indicatori di caratterizzazione della miscelazione

Deviazione standard normalizzata: $\Delta c = \sqrt{\frac{\langle[c(\mathbf{r},t) - c_\infty]^2\rangle}{\langle[c(\mathbf{r},0) - c_\infty]^2\rangle}}$

dove $c_\infty$ è la concentrazione nello stato di miscelazione completa.

Risultati sperimentali

Transizione dello stato di flusso

1. Transizione subcritica: Scoperta che la transizione da flusso laminare a turbolenza elastica presenta carattere subcritico, con isteresi presente
2. Numero di Weissenberg critico:
   • Ramo ascendente: $Wi_c = 9.95$
   • Ramo discendente: $Wi_c^I = 9.26$
3. Stato debolmente caotico: Osservazione di uno stato di flusso debolmente caotico al di sotto di $Wi_c^I$

Analisi delle caratteristiche di miscelazione

Stato di flusso debolmente caotico ($5.03 < Wi \leq 8.17$)

• Decadimento della deviazione standard secondo legge di potenza con esponente prossimo a $-1/6$
• Lo spettro di potenza mostra legge di scala con esponente $-5$ o superiore

Stato di turbolenza elastica ($Wi \geq 9.42$)

• Decadimento iniziale secondo legge di potenza (esponente circa $1/5$), seguito da decadimento esponenziale
• Forma di decadimento esponenziale: $\Delta c \sim \exp(-\alpha_{mix}t)$
• Lo spettro di potenza presenta caratteristiche dello spettro di Batchelor (scala $k^{-1}$)

Relazioni di scala

1. Tempo di miscelazione: $t_{mix} \sim (Wi - Wi_c^I)^{-0.52}$
2. Concentrazione regione II: $c_{II}/c_\infty \sim (Wi - Wi_c^I)^{0.28}$
3. Relazione di scala: $c_{II} \sim t_{mix}^{-1/2}$

Correlazione tra parametro d'ordine e tasso di miscelazione

I risultati sperimentali mostrano una forte correlazione tra il parametro d'ordine $\Phi_d$ e il tasso di miscelazione $\alpha_{mix}$, indicando che l'aumento delle fluttuazioni di velocità migliora direttamente l'efficienza di miscelazione del fluido.

Lavori correlati

Ricerca sulla turbolenza elastica

• Prima scoperta della turbolenza elastica in soluzioni polimeriche da parte di Groisman & Steinberg
• Analisi di stabilità lineare e non lineare in geometria Taylor-Couette
• Osservazione di modelli di flusso elicoidale e a bande nelle simulazioni numeriche tridimensionali

Ricerca sull'aumento della miscelazione

• Miscelazione di flusso estensionale in laminatoi a quattro rulli
• Instabilità elastica in geometria a croce
• Aumento del trasporto di calore convettivo in canali a serpentina

Conclusioni e discussione

Conclusioni principali

1. Natura della transizione: Conferma della natura subcritica della transizione di turbolenza elastica nel flusso Taylor-Couette
2. Meccanismo di miscelazione: Verifica che la turbolenza elastica possa realizzare miscelazione efficiente a numero di Reynolds estremamente basso
3. Regione di Batchelor: Osservazione di caratteristiche di miscelazione di Batchelor nello stato di turbolenza elastica
4. Indicatore predittivo: Il parametro d'ordine può servire come indicatore affidabile dell'efficienza di miscelazione

Limitazioni

1. Limitazione bidimensionale: Lo studio è limitato a geometria bidimensionale; gli effetti tridimensionali potrebbero produrre risultati diversi
2. Dipendenza dalla griglia: Il comportamento di isteresi scompare con griglie più fini, indicando l'influenza della precisione numerica
3. Sensibilità alle condizioni iniziali: Le caratteristiche di miscelazione dipendono dalla scelta della distribuzione di concentrazione iniziale

Direzioni future

1. Effetto delle onde elastiche: Studio dell'influenza delle onde elastiche su efficienza, tasso e tempo di miscelazione
2. Estensione tridimensionale: Verifica delle leggi di scala scoperte in geometrie tridimensionali
3. Applicazioni pratiche: Esplorazione del potenziale di applicazione pratica nei dispositivi microfluidici

Valutazione approfondita

Punti di forza

1. Ricerca sistematica: Analisi completa delle caratteristiche di miscelazione durante la transizione da flusso laminare a turbolenza elastica
2. Caratterizzazione quantitativa: Introduzione di molteplici indicatori quantitativi per caratterizzare il grado e l'efficienza di miscelazione
3. Verifica teorica: I risultati numerici sono coerenti con la teoria di miscelazione di Batchelor e le osservazioni sperimentali
4. Valore pratico: Fornisce guida teorica importante per applicazioni microfluidiche

Carenze

1. Problemi di precisione numerica: La dipendenza dalla griglia del comportamento di isteresi suggerisce limitazioni potenziali del metodo numerico
2. Comprensione del meccanismo: La spiegazione del meccanismo fisico del decadimento secondo legge di potenza nello stato debolmente caotico non è sufficientemente approfondita
3. Intervallo di parametri: L'intervallo di parametri studiato è relativamente limitato, mancando un'esplorazione più ampia dello spazio dei parametri

Impatto

1. Contributo teorico: Fornisce importanti intuizioni per la comprensione dei meccanismi di miscelazione nei fluidi viscoelastici
2. Prospettive di applicazione: Ha significato guida importante per la progettazione di dispositivi microfluidici
3. Valore metodologico: I metodi di caratterizzazione stabiliti possono essere generalizzati ad altre geometrie di flusso

Scenari applicabili

• Ottimizzazione della miscelazione nei dispositivi microfluidici
• Controllo della miscelazione nei processi di lavorazione dei polimeri
• Miscelazione a microscala in applicazioni biomediche

Bibliografia

L'articolo cita 57 riferimenti correlati, coprendo molteplici aspetti tra cui turbolenza elastica, stabilità del flusso viscoelastico e teoria della miscelazione, fornendo una solida base teorica per la ricerca.