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

Titel: Oplossing van een tweecomponentendruppel op een macrofase door oppervlaktespanningseffect

Auteur: Alexey Kabalnov (Ink Splat Company)

1. Probleemstelling

Het artikel onderzoekt het fysisch-chemische probleem van de oplossing (dissolutie) van een bolvormige vloeibare druppel bestaande uit twee volledig mengbare componenten in een omringend medium. De drijvende kracht voor deze oplossing is de excess Laplace-druk (capillaire druk) binnen de druppel, veroorzaakt door oppervlaktespanning.

De kernvraag is hoe de aanwezigheid van een moeilijk oplosbare component (additief) de oplossingssnelheid beïnvloedt. Het is bekend dat dergelijke additieven Ostwald-rijping (groei van grote deeltten ten koste van kleine) kunnen vertragen of stoppen, maar de onderliggende kinetiek voor een mengsel met een eindige oplosbaarheid van het additief was niet volledig gekwantificeerd. Het doel is om een model te ontwikkelen dat de oplossingssnelheid voorspelt over het volledige bereik van samenstellingen (van zuivere component 1 tot zuivere component 2), rekening houdend met niet-ideale oplossingen en verschillen in molaire volumes.

2. Methodologie

De auteur hanteert een combinatie van analytische modellering en numerieke simulaties:

• Diffusie- en Chemische Potentiaaltheorie: Het model begint met de diffusievergelijkingen voor twee componenten (J1 en J2) rond een bolvormige deeltje. De concentratieverschillen aan het oppervlak worden gekoppeld aan de chemische potentiaalverschillen ($\Delta\mu$) tussen de druppel en de omringende macrofase.
• Laplace- en Raoult-effecten: De chemische potentiaal wordt ontbonden in twee bijdragen:
  1. Het Laplace-effect (positief, drijft oplossing door capillaire druk).
  2. Het Raoult-effect (negatief voor de goed oplosbare component, veroorzaakt door de toename van de concentratie van de slecht oplosbare component in de druppel).
• Lock-in Benadering: Een cruciale aanname is de "lock-in" (vergrendeling) toestand. Hierin wordt aangenomen dat het Raoult-effect de Laplace-druk voor de goed oplosbare component bijna volledig compenseert ($\Delta\mu_1 \approx 0$). Hierdoor wordt de diffusie van de moeilijk oplosbare component (component 2) de snelheidsbepalende stap.
• Numerieke Validatie: De analytische afleidingen worden getoetst door numerieke simulaties van de gekoppelde differentiaalvergelijkingen voor verschillende scenario's (variatie in oppervlaktespanning, molaire volumes en initiële samenstelling).

3. Belangrijkste Bijdragen en Concepten

A. Drie Stadia van Oplossing

Het paper identificeert drie distincte fasen in het levensverloop van de druppel:

1. Pre-lock-in: Een korte initiële fase waarin de concentratie van de moeilijk oplosbare component snel toeneemt. De Laplace-druk domineert nog.
2. Lock-in (Steady State): De belangrijkste fase. De toename van de concentratie van component 2 zorgt ervoor dat het Raoult-effect de Laplace-druk voor component 1 neutraliseert. De oplossingssnelheid wordt nu volledig bepaald door de trage "bloeding" (diffusie) van component 2. De samenstelling van de druppel blijft hierbij nagenoeg constant.
3. Late lock-in: Aan het einde van de levensduur van het deeltje (kleine straal) wordt de Laplace-druk zo groot dat het Raoult-effect niet meer volledig kan compenseren. De samenstelling verandert sneller en de kinetiek wijkt af van de lineaire trend.

B. Lock-in Aantallen (Dimensionless Numbers)

De auteur introduceert twee dimensieloze grootheden om de geldigheid van het model te kwantificeren:

• Eerste Lock-in Getal ($L_1$): Bepaalt of het systeem de lock-in toestand bereikt en hoe lang deze duurt. $L_1 \gg 1$ betekent dat de samenstelling tijdens de meeste levensduur constant blijft en de klassieke kubieke wet geldt.
• Tweede Lock-in Getal ($L_2$): Gerelateerd aan de verhouding van de diffusie en oplosbaarheid, bepaalt het de mate waarin de "bloeding" van component 2 de stabiliteit beïnvloedt.

C. Verbeterde Empirische Vergelijking

De paper presenteert een nieuwe, uitgebreide vergelijking voor de oplossingssnelheid ($w$) die de hele samenstellingsrange dekt (van $x_2=0$ tot $x_2=1$):
$$w = \frac{4\pi D_2 C_2^* V_m}{\left( \frac{\phi_1}{C_1^* V_{m1}} + \frac{\phi_2}{C_2^* V_{m2}} \right)^{-1}}$$
(Let op: De exacte vorm in vergelijking 31 combineert de bijdragen van beide componenten en correcties voor molaire volumes en activiteitscoëfficiënten).

Deze vergelijking is een uitbreiding van eerdere modellen (zoals die van Webster en Cates) en is geldig voor:

• Niet-ideale oplossingen.
• Verschillende molaire volumes van de componenten.
• Het volledige bereik van volumefracties.

4. Resultaten

• Lineariteit van Volume-tijd: Numerieke simulaties tonen aan dat voor $L_1 > 10$ de relatie tussen het deeltjesvolume en de tijd bijna perfect lineair is (wat overeenkomt met de klassieke kubieke wet voor de straal: $r^3 \sim t$).
• Invloed van Oppervlaktespanning: Bij hoge oppervlaktespanning (en dus lage $L_1$) wordt de "pre-lock-in" fase langer en treden er afwijkingen op in de late fase, maar de algemene trend blijft voorspelbaar.
• Validatie van de Empirische Formule: De nieuwe formule (Eq. 31) voorspelt de oplossingssnelheid met een nauwkeurigheid van ongeveer 20% zelfs bij zeer asymmetrische molaire volumes (factor 8 verschil). Het oude model (Eq. 36) faalde hierbij met afwijkingen tot een factor 5.
• Bleeding Effect: Zelfs als de "bloeding" van de slecht oplosbare component significant is (lage $L_2$), blijft de kinetiek kubisch, maar verandert de effectieve oplossingssnelheid. De stabilisatie treedt niet op door een verandering in de kinetische wet, maar door een verandering in de snelheidsconstante.

5. Significantie en Toekomstperspectief

• Fundamenteel Begrip: Het paper levert een volledig theoretisch raamwerk voor het begrijpen van hoe additieven de oplossing van druppels vertragen, wat essentieel is voor de stabiliteit van emulsies, inkten en farmaceutische formuleringen.
• Ostwald-rijping: De resultaten vormen de basis voor een vervolgpapier over Ostwald-rijping. Het model suggereert dat voor $L_1 \gg 1$ de rijpingssnelheid volgt uit de afgeleide vergelijkingen, terwijl voor $L_1 \approx 1$ een bimodale deeltjesgrootteverdeling kan ontstaan.
• Praktische Toepassing: De afgeleide vergelijkingen maken het mogelijk om formuleringen te ontwerpen die optimaal bestand zijn tegen coalescentie en rijping, door de juiste balans te kiezen tussen oplosbaarheid, oppervlaktespanning en molaire volumes van de componenten.

Conclusie:
Kabalnov heeft een robuust model ontwikkeld dat de complexe interactie tussen Laplace-druk en Raoult-effecten kwantificeert. Door de introductie van de "lock-in" toestand en de bijbehorende dimensieloze getallen, kan de oplossingssnelheid van tweecomponentendruppels nauwkeurig worden voorspeld over een breed scala aan condities, wat een belangrijke stap voorwaarts is in de colloïdechemie en formuleringstechnologie.