Dissolution of a two-component drop onto macrophase due to surface tension effect

Titolo

Dissoluzione di una goccia bifase su una macrofase dovuta all'effetto della tensione superficiale

1. Il Problema

Il lavoro affronta il problema fisico della dissoluzione di una goccia liquida sferica composta da due componenti completamente miscibili (Componente 1 e Componente 2) immersa in un mezzo continuo (macrophase).

• Contesto: Si studia come gli additivi scarsamente solubili possano rallentare o inibire l'ingrossamento di Ostwald (Ostwald ripening) nei sistemi dispersi.
• Ipotesi di partenza: Il Componente 2 è significativamente meno solubile nel mezzo rispetto al Componente 1 ($C_2^* \ll C_1^*$).
• Meccanismo: La dissoluzione è guidata dalla pressione di Laplace eccessiva all'interno della goccia (dovuta alla tensione superficiale), che aumenta il potenziale chimico dei componenti, spingendoli a diffondere verso la macrofase.
• Obiettivo: Analizzare la cinetica di dissoluzione considerando l'intero intervallo di composizioni (frazioni molari da 0 a 1), i volumi molari dei componenti e i coefficienti di attività (soluzioni non ideali), superando i limiti dei modelli precedenti che assumevano additivi completamente insolubili o soluzioni ideali.

2. Metodologia

L'autore sviluppa un modello teorico basato su equazioni di diffusione e termodinamica, integrato da simulazioni numeriche per la validazione.

• Equazioni Fondamentali:
  • Si utilizzano le equazioni di flusso di massa stazionario per una sfera (legge di Fick modificata per la geometria sferica).
  • I potenziali chimici ($\Delta\mu$) sono decomposti in due contributi: l'effetto Laplace (pressione capillare) e l'effetto Raoult (variazione di composizione).
  • Si applica l'equazione di Gibbs-Duhem per relazionare i potenziali chimici dei due componenti senza assumere idealità della soluzione.
• Approssimazione "Lock-in" (Blocco):
  • Si identifica uno stato stazionario in cui l'effetto Raoult (dovuto all'accumulo del componente meno solubile) bilancia quasi perfettamente l'effetto Laplace per il componente più solubile.
  • In questo stato, la dissoluzione è controllata dalla lenta "sanguinamento" (bleeding) del componente scarsamente solubile.
• Simulazioni Numeriche:
  • Vengono risolte direttamente le equazioni differenziali accoppiate per i potenziali chimici e i flussi di massa, senza assumere a priori che $\Delta\mu_1 = 0$.
  • Vengono testati diversi valori di tensione superficiale, frazioni volumetriche iniziali, volumi molari e solubilità per validare le approssimazioni analitiche.

3. Contributi Chiave

Il paper introduce concetti teorici nuovi e generalizza modelli esistenti:

1. Identificazione di Tre Stadi Cinetici:

   • Pre-lock-in: Fase iniziale rapida in cui la concentrazione del componente scarsamente solubile aumenta velocemente per bilanciare la pressione di Laplace.
   • Lock-in (Stato di blocco): Fase stazionaria in cui la cinetica di dissoluzione è controllata dal componente meno solubile. La composizione della goccia rimane quasi costante e la cinetica segue una legge cubica classica.
   • Late lock-in: Fase finale in cui, quando il raggio diventa molto piccolo, le variazioni relative di concentrazione diventano significative, portando a una deviazione dalla linearità (termine quadratico nell'equazione di volume).

2. Numeri di Lock-in ($L_1$ e $L_2$):

   • $L_1$ (Primo Numero di Lock-in): Un gruppo adimensionale che determina la validità dell'approssimazione di composizione costante. Se $L_1 \gg 1$, la deviazione dalla legge cubica è trascurabile.
   • $L_2$ (Secondo Numero di Lock-in): Un rapporto che caratterizza il confine tra i regimi di dissoluzione dominati dal componente 1 o dal componente 2.

3. Nuova Equazione Empirica Generalizzata:

   • Viene proposta una nuova equazione (Eq. 31 nel testo) per il tasso di dissoluzione che unifica i regimi di composizione pura e mista.
   • Questa equazione estende il lavoro precedente [2] includendo esplicitamente le differenze nei volumi molari e la non idealità della soluzione (coefficienti di attività), aspetti che i modelli precedenti non gestivano correttamente.

4. Risultati

Le simulazioni numeriche e l'analisi teorica confermano quanto segue:

• Validità della Legge Cubica: Se la concentrazione iniziale del componente scarsamente solubile è superiore a una certa soglia (corrispondente a $L_1 > 10$), la cinetica di dissoluzione segue quasi perfettamente la legge cubica classica ($V \propto t$), nonostante le variazioni di composizione.
• Effetto della Tensione Superficiale:
  • Con alta tensione superficiale (basso $L_1$), si osservano fasi di pre-lock-in e late-lock-in più lunghe e deviazioni dalla linearità nel grafico Volume-Tempo.
  • Con bassa tensione superficiale (alto $L_1$), il sistema entra rapidamente nello stato di lock-in e la linearità è quasi perfetta.
• Robustezza dell'Equazione Proposta: L'equazione empirica (Eq. 31) mostra un accordo eccellente (entro il 20%) con le simulazioni numeriche anche per volumi molari molto asimmetrici (fattore di 8 di differenza), dove i modelli precedenti fallivano con errori fino a un fattore 5.
• Soglia di Blocco: Il "blocco" completo si verifica solo se la frazione molare iniziale del componente meno solubile è sufficientemente alta (circa > 0.1 nel caso simulato). Al di sotto di questa soglia, il componente "sangina" troppo velocemente per stabilire un equilibrio di lock-in completo.

5. Significato e Implicazioni

Questo studio è fondamentale per la scienza dei colloidi e delle emulsioni:

• Stabilizzazione delle Dispersioni: Fornisce una base teorica rigorosa per spiegare come gli additivi scarsamente solubili inibiscano l'ingrossamento di Ostwald, permettendo di progettare formulazioni più stabili (es. inchiostri, cosmetici, farmaci).
• Generalizzazione del Modello: Estende la teoria di Ostwald ripening a sistemi reali, non ideali e con componenti di diverse dimensioni molecolari, colmando un vuoto nella letteratura scientifica.
• Prospettive Future: Il modello getta le basi per il prossimo articolo della serie, che applicherà questi risultati alla teoria completa dell'ingrossamento di Ostwald. In particolare, suggerisce che per $L_1 \approx 1$ potrebbe emergere una distribuzione bimodale delle dimensioni delle particelle, un fenomeno non previsto dalle teorie classiche (Lifshits-Slezov-Wagner).

In sintesi, il paper offre un quadro completo e matematicamente solido per prevedere la vita utile e la stabilità cinetica di gocce bifase, introducendo parametri adimensionali critici per il controllo della formulazione.