La diffusione transitoria solido-liquido è governata dalle leggi di Fick, dalla conservazione della massa e dalla differenza di concentrazione variabile tra le fasi solida e liquida. All’interno del solido, il soluto si muove per diffusione verso la superficie; l’agitazione trasporta il soluto attraverso il liquido e mantiene il liquido di massa approssimativamente uniforme. Poiché il volume del liquido è finito, la sua concentrazione aumenta durante l’estrazione, riducendo la forza motrice fino al raggiungimento dell’equilibrio o di un recupero prestabilito.

Concetto fondamentale: Un modello di estrazione su scala pilota deve accoppiare la diffusione transitoria all’interno del solido con un bilancio dinamico della massa nella fase liquida. La diffusione interna, la dimensione delle particelle, la diffusività effettiva, l’equilibrio interfacciale, il volume del liquido e l’agitazione determinano insieme la velocità di estrazione e il tempo di batch richiesto.

Il modello fisico governante

La prima legge di Fick determina il flusso locale

La prima legge di Fick mette in relazione il flusso di soluto con il gradiente di concentrazione:

[ \mathbf{N}A=-D{\mathrm{eff}}\nabla C_A ]

Qui, (D_{\mathrm{eff}}) è la diffusività effettiva del soluto (A) attraverso la matrice solida, mentre (C_A) è la sua concentrazione all’interno del solido.

Il segno negativo indica che il soluto si sposta dalle regioni a concentrazione elevata verso quelle a concentrazione bassa. In una particella porosa, (D_{\mathrm{eff}}) riflette gli effetti della struttura dei pori, della porosità, della tortuosità, della saturazione, della temperatura e della viscosità del liquido.

La seconda legge di Fick descrive l’esaurimento transitorio

Per una diffusività costante e in assenza di reazioni all’interno del solido, il campo di concentrazione transitorio soddisfa:

[ \frac{\partial C_A}{\partial t}

D_{\mathrm{eff}}\nabla^2 C_A ]

Per una particella sferica, una geometria comune nei modelli di lisciviazione, si ottiene:

[ \frac{\partial C_A}{\partial t}

D_{\mathrm{eff}} \left[ \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2\frac{\partial C_A}{\partial r} \right) \right] ]

La geometria è importante perché lastre piane, cilindri e sfere presentano lunghezze del percorso di diffusione e profili di concentrazione diversi.

La condizione iniziale rappresenta il solido caricato

Il modello inizia con una distribuzione iniziale del soluto, spesso approssimata come uniforme:

[ C_A(r,0)=C_{A,0} ]

I solidi più complessi possono presentare un profilo radiale di concentrazione causato da precedenti operazioni di essiccamento, impregnazione, reazione o precipitazione. Tale condizione iniziale può influenzare significativamente la velocità di estrazione nelle fasi iniziali.

La simmetria vincola il centro della particella

Per una particella sferica, nessun flusso netto attraversa il centro:

[ \left. \frac{\partial C_A}{\partial r} \right|_{r=0}=0 ]

Questa condizione impedisce una singolarità artificiale al centro e riflette la simmetria fisica della particella.

Come entra nel modello il volume finito del liquido

Il liquido di massa viene spesso trattato come perfettamente miscelato

Con un’agitazione vigorosa, la concentrazione del liquido di massa può essere approssimata come uniforme nello spazio:

[ C_L=C_L(t) ]

Questa ipotesi elimina la necessità di risolvere un campo completo di concentrazione nella fase liquida. Non elimina il trasferimento di massa; significa che i gradienti di concentrazione sul lato liquido sono rappresentati attraverso una resistenza interfacciale oppure si presume che siano sufficientemente piccoli.

Il bilancio del liquido accoppia entrambe le fasi

Il soluto rilasciato dal solido si accumula nel volume finito del liquido:

[ V_L\frac{dC_L}{dt}

\dot{m}_{A,\mathrm{solid\ to\ liquid}} ]

Per (N_p) particelle identiche con area superficiale (A_p), ciò può essere scritto come:

[ V_L\frac{dC_L}{dt}

N_p A_p N_A ]

dove (N_A) è il flusso di soluto verso l’esterno sulla superficie della particella.

Il bilancio totale del soluto fornisce un importante controllo di coerenza:

[ M_{A,0}

M_A^{\mathrm{solid}}(t) + V_L C_L(t) ]

assumendo che non vi siano perdite di soluto, degradazione o rimozione tramite campionamento.

Il volume finito crea una forza motrice decrescente

Una condizione interfacciale comunemente utilizzata è:

[ -D_{\mathrm{eff}} \left. \frac{\partial C_A}{\partial r} \right|_{r=R}

k_L\left(C_{A,i}-C_L\right) ]

dove (C_{A,i}) è la concentrazione interfacciale equivalente nel liquido e (k_L) è il coefficiente di trasferimento di massa sul lato liquido.

Se l’equilibrio locale è rappresentato da una relazione di ripartizione,

[ C_{A,i}=K_{\mathrm{eq}}C_{A,s} ]

allora la concentrazione superficiale nel solido e la concentrazione nel liquido sono collegate attraverso (K_{\mathrm{eq}}). All’aumentare di (C_L), la differenza tra la concentrazione del liquido all’equilibrio e (C_L) diminuisce, per cui l’estrazione rallenta.

I principali principi adimensionali

Il numero di Fourier misura il progresso della diffusione

Il tempo di diffusione adimensionale è rappresentato dal numero di Fourier:

[ \mathrm{Fo}=\frac{D_{\mathrm{eff}}t}{R^2} ]

Qui, (R) è il raggio caratteristico della particella. Valori piccoli indicano che la diffusione ha penetrato solo una parte del solido, mentre valori più elevati indicano un esaurimento più esteso.

Il tempo caratteristico di diffusione scala come:

[ t_D\sim \frac{R^2}{D_{\mathrm{eff}}} ]

Questa dipendenza quadratica rende la dimensione delle particelle particolarmente importante nella progettazione di impianti pilota. Dimezzare il raggio delle particelle può ridurre il tempo caratteristico di diffusione interna di circa un fattore quattro, purché le altre limitazioni rimangano invariate.

Il numero di Biot confronta la resistenza esterna e quella interna

Un numero di Biot di trasferimento di massa confronta la resistenza sul lato liquido con la diffusione interna nel solido:

[ \mathrm{Bi}m=\frac{k_L R}{D{\mathrm{eff}}} ]

Un valore elevato di (\mathrm{Bi}_m) indica generalmente che la resistenza esterna è relativamente piccola rispetto alla resistenza alla diffusione interna. Un valore ridotto di (\mathrm{Bi}_m) indica che lo strato limite liquido può controllare o influenzare significativamente la velocità di estrazione.

La definizione esatta può variare in funzione della geometria e delle convenzioni di ripartizione, ma la domanda di fondo è la stessa: dove si verifica la maggior parte della caduta di concentrazione?

Le soluzioni di similitudine spiegano il comportamento iniziale

Per la diffusione seminfinità all’interno o all’esterno di un solido, la variabile di similitudine è:

[ \eta=\frac{x}{2\sqrt{D_{\mathrm{eff}}t}} ]

L’equazione della diffusione può quindi produrre profili di concentrazione basati sulla funzione errore, come:

[ \Theta(\eta)=\operatorname{erfc}(\eta) ]

Queste soluzioni mostrano che la profondità di penetrazione cresce con la radice quadrata del tempo:

[ \delta(t)\propto \sqrt{D_{\mathrm{eff}}t} ]

Sono utili per comprendere l’estrazione iniziale da solidi di grandi dimensioni o da strati superficiali sottili, sebbene siano necessarie soluzioni per particelle finite quando il fronte di diffusione raggiunge il centro della particella.

Il flusso iniziale spesso decresce come (t^{-1/2})

Per un solido sem infinito inizialmente uniforme esposto a una concentrazione superficiale costante, il flusso superficiale ha la forma:

[ N_A\propto \left(\frac{D_{\mathrm{eff}}}{\pi t}\right)^{1/2} \Delta C ]

La dipendenza (t^{-1/2}) descrive la caratteristica diminuzione del flusso di diffusione transitorio. In un estrattore a volume finito, la diminuzione è solitamente più accentuata o comunque modificata, poiché la concentrazione del liquido cambia e la differenza di concentrazione disponibile diminuisce.

Applicazione del modello agli impianti pilota

Il tempo di estrazione è una variabile progettuale accoppiata

Un impianto pilota deve garantire un tempo di contatto sufficiente affinché la frazione target di soluto abbandoni il solido. Tale tempo dipende da (D_{\mathrm{eff}}), dal raggio delle particelle, dalla temperatura, dalle proprietà del solvente, dal rapporto liquido-solido, dal comportamento all’equilibrio e dall’agitazione.

Una semplice stima basata soltanto su (R^2/D_{\mathrm{eff}}) è utile per il dimensionamento iniziale, ma non può prevedere il recupero finale quando l’accumulo nel liquido e le limitazioni all’equilibrio sono importanti.

L’agitazione controlla principalmente il trasporto esterno

L’agitazione riduce lo strato limite liquido che circonda le particelle e aumenta il coefficiente di trasferimento di massa sul lato liquido (k_L). Può inoltre mantenere le particelle in sospensione, prevenire la sedimentazione, migliorare l’uniformità della temperatura e ridurre i gradienti di concentrazione nel recipiente.

Quando la resistenza esterna diventa ridotta, un’ulteriore agitazione produce miglioramenti decrescenti, poiché la diffusione interna nel solido è diventata la resistenza dominante.

La dimensione delle particelle influenza fortemente la produttività

La riduzione della dimensione delle particelle accorcia il percorso caratteristico di diffusione e aumenta l’area superficiale esterna disponibile per unità di massa. Questi cambiamenti generalmente accelerano l’estrazione, ma possono anche aumentare il consumo energetico, generare particelle fini, complicare la filtrazione e causare problemi di movimentazione o sicurezza.

È quindi necessario effettuare prove pilota per individuare la dimensione minima pratica delle particelle, invece di presumere che la macinazione massima sia ottimale.

Il volume del liquido controlla la forza motrice

Un volume maggiore di solvente può mantenere più bassa la concentrazione del liquido di massa per una determinata quantità di soluto estratto. Ciò mantiene una forza motrice maggiore e può aumentare il recupero finale quando la ripartizione all’equilibrio limiterebbe altrimenti l’estrazione.

Il compromesso consiste in un maggiore utilizzo di solvente, in un carico di concentrazione a valle più elevato e potenzialmente in un volume maggiore delle apparecchiature. Il rapporto liquido-solido appropriato è quindi una variabile di ottimizzazione, non una costante universale.

I dati pilota collegano la teoria allo scale-up

Una campagna pilota efficace misura la concentrazione del liquido di massa in funzione del tempo, registrando al contempo la dimensione delle particelle, la temperatura, la velocità di agitazione, il volume del solvente e il carico di solido. Questi dati possono essere adattati a un modello di diffusione transitoria per stimare (D_{\mathrm{eff}}), (k_L), i parametri di equilibrio o loro combinazioni.

La sola velocità di agitazione non è un criterio trasferibile per lo scale-up. Lo scale-up dovrebbe preservare il comportamento di trasporto dominante, utilizzando grandezze pertinenti come il tempo di miscelazione, la potenza immessa, la qualità della sospensione, il numero di Reynolds o coefficienti di trasferimento di massa determinati sperimentalmente.

Comprendere i compromessi

L’ipotesi di perfetta miscelazione ha dei limiti

Un liquido può sembrare vigorosamente miscelato pur contenendo gradienti locali di concentrazione vicino alle particelle, alle pareti del recipiente o in regioni non adeguatamente sospese. Se il volume del liquido è grande o la sospensione è densa, un (C_L(t)) uniforme nello spazio può essere un’approssimazione inaccurata.

Quando l’uniformità del liquido di massa è incerta, si possono utilizzare misure del tempo di miscelazione, prove con traccianti oppure modelli a compartimenti e di fluidodinamica computazionale.

La diffusione interna potrebbe non essere l’unico meccanismo

La velocità apparente di estrazione può includere il rapido lavaggio superficiale, la diffusione nei pori, il rigonfiamento, la frammentazione delle particelle, la dissoluzione, la reazione chimica, l’adsorbimento o il desorbimento. Trattare ogni deviazione dal comportamento fickiano ideale come una variazione della diffusività può produrre parametri adattati fuorvianti.

Un modello di reazione-diffusione o un modello a particella in restringimento può essere necessario quando la struttura del solido cambia durante l’estrazione.

La diffusività effettiva è solitamente una proprietà apparente

(D_{\mathrm{eff}}) non è necessariamente la diffusività molecolare misurata in un liquido puro. Spesso incorpora la geometria dei pori, la tortuosità, la costrizione, la saturazione, la ripartizione e le interazioni con la matrice solida.

Il suo valore può cambiare con la temperatura, il contenuto di umidità, la concentrazione del soluto o il grado di esaurimento. Un coefficiente adattato alla scala di laboratorio dovrebbe quindi essere convalidato in condizioni pilota.

L’equilibrio può limitare il recupero finale

Aumentare l’agitazione non può superare una relazione di equilibrio sfavorevole. Se il liquido si avvicina alla concentrazione di equilibrio, la velocità di estrazione diventa bassa anche quando la miscelazione è eccellente.

Il funzionamento in controcorrente, il rinnovo del solvente, l’estrazione a stadi, la regolazione del pH o una reazione chimica selettiva possono migliorare la forza motrice complessiva quando questi approcci sono compatibili con il prodotto e il processo.

I sistemi multicomponente richiedono maggiore attenzione

Per l’estrazione diluita di un singolo soluto, un modello fickiano o pseudobinario è spesso adeguato. Nei sistemi concentrati o multicomponente, le interazioni tra specie, l’adsorbimento competitivo, la termodinamica non ideale e la diffusione accoppiata possono richiedere descrizioni del trasporto multicomponente secondo Stefan-Maxwell o altri approcci.

Il modello dovrebbe essere adeguato all’intervallo di composizione e al livello di accuratezza richiesto per la decisione relativa all’impianto pilota.

Come applicare questo al proprio progetto

Un modello di estrazione transitoria dovrebbe essere sviluppato a partire da proprietà misurabili e convalidato mediante dati di concentrazione in funzione del tempo. Le seguenti priorità aiutano ad allineare il modello all’obiettivo ingegneristico:

Se l’obiettivo principale è il tempo di estrazione: Stimare (R^2/D_{\mathrm{eff}}), misurare le distribuzioni delle dimensioni delle particelle e adattare il modello completo di diffusione transitoria nel solido per determinare il tempo di batch necessario al recupero richiesto.
Se l’obiettivo principale è valutare le prestazioni dell’agitazione: Variare la velocità dell’agitatore o le condizioni di fluidizzazione misurando al contempo la velocità di estrazione, la qualità della sospensione e il tempo di miscelazione, per distinguere la resistenza sul lato liquido dalla diffusione interna.
Se l’obiettivo principale è il consumo di solvente: Utilizzare il bilancio del liquido a volume finito e la relazione di equilibrio per selezionare il rapporto liquido-solido o il numero di stadi di estrazione che garantisca il recupero con un consumo di solvente accettabile.
Se l’obiettivo principale è lo scale-up dall’impianto pilota alla produzione: Determinare se il processo è controllato dalla diffusione interna, dal trasferimento di massa esterno, dalla miscelazione o dall’equilibrio, quindi preservare il comportamento adimensionale dominante e convalidarlo su scala maggiore.
Se l’obiettivo principale è l’affidabilità del modello: Verificare la conservazione totale della massa, misurare la temperatura e la viscosità, caratterizzare la porosità e la dimensione delle particelle e verificare che (D_{\mathrm{eff}}) e (k_L) rimangano validi nell’intervallo operativo.

Il compito ingegneristico centrale consiste nell’identificare la resistenza dominante e progettare l’impianto pilota in modo che le sue misure distinguano gli effetti della diffusione, della miscelazione e dell’equilibrio.

Tabella riepilogativa:

Principio	Ruolo nella progettazione dell’estrazione
Prima legge di Fick	Calcola il flusso di soluto a partire dai gradienti di concentrazione.
Seconda legge di Fick	Prevede l’esaurimento transitorio del soluto all’interno delle particelle.
Bilancio del liquido finito	Segue l’aumento della concentrazione di massa e la diminuzione della forza motrice.
Numero di Fourier	Collega il progresso della diffusione al tempo e alla dimensione delle particelle.
Numero di Biot per il trasferimento di massa	Confronta la resistenza esterna del liquido con la diffusione interna.
Raggio delle particelle	Influenza fortemente il tempo di estrazione, approssimativamente proporzionale a R^2/D_eff.
Agitazione	Aumenta il trasferimento di massa sul lato liquido, la sospensione e l’uniformità della temperatura.
Relazione di equilibrio	Determina il recupero massimo raggiungibile.
Dati pilota	Convalidano D_eff, k_L, i parametri di equilibrio e le ipotesi di scale-up.

Migliora il tuo impianto pilota di estrazione

Devi determinare il tempo di estrazione, il volume di solvente, i requisiti di agitazione o le condizioni di scale-up? Contatta LABPARK per discutere di un impianto pilota per operazioni unitarie a scopo didattico o professionale, personalizzato per il tuo progetto.

LABPARK aiuta università, istituti di ricerca e imprese a valutare l’estrazione solido-liquido utilizzando apparecchiature pratiche su scala pilota per applicazioni di ingegneria chimica, bioprocessi e biotecnologie, nonché ambientali e di trattamento delle acque. I nostri sistemi supportano studi controllati su dimensione delle particelle, temperatura, rapporto liquido-solido, agitazione, miscelazione, trasferimento di massa ed equilibrio, affinché il tuo team possa identificare la resistenza dominante e generare dati affidabili di concentrazione in funzione del tempo.

Collabora con LABPARK per collegare la teoria della diffusione transitoria alle prestazioni misurabili del processo, migliorare l’efficienza nell’uso di solvente ed energia, convalidare i parametri di scale-up e selezionare condizioni operative che raggiungano il recupero target.