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

Título: Dissolução de uma gota de dois componentes em uma macrofase devido ao efeito da tensão superficial

Autor: Alexey Kabalnov (Ink Splat Company)

1. Problema Investigado

O artigo aborda o problema fundamental da dissolução de uma gota líquida esférica composta por dois componentes miscíveis em um meio contínuo (macrofase). O foco principal é entender como a presença de um componente pouco solúvel (aditivo) afeta a cinética de dissolução induzida pela pressão de Laplace excessiva (efeito de tensão superficial).

Este estudo é uma extensão e revisão de modelos anteriores sobre a estabilização contra o Ostwald ripening (amadurecimento de Ostwald). Enquanto trabalhos anteriores focaram em aditivos completamente insolúveis, este artigo investiga o caso onde o aditivo possui solubilidade finita, cobrindo todo o espectro de composições (fração molar de 0 a 1), volumes molares diferentes e coeficientes de atividade não ideais.

2. Metodologia

A abordagem combina análise teórica baseada em termodinâmica e difusão com simulações numéricas para validação.

• Formulação Termodinâmica:
  • O sistema é modelado considerando o equilíbrio de potenciais químicos entre a gota e a macrofase.
  • A força motriz para a dissolução é a diferença de potencial químico ($\Delta\mu$) causada pela pressão de Laplace ($2\sigma/r$) e pela variação na composição (efeito de Raoult).
  • São estabelecidas equações de fluxo de massa estacionário ($J_1, J_2$) baseadas nos coeficientes de difusão ($D$) e nas diferenças de concentração na superfície da partícula e no infinito.
• Aproximação "Lock-in" (Bloqueio):
  • O modelo assume que, após uma fase inicial rápida, o sistema entra em um estado de "lock-in". Neste estado, o efeito de Raoult (devido ao aumento da concentração do componente menos solúvel) quase que completamente compensa o efeito da pressão de Laplace para o componente mais solúvel.
  • Consequentemente, o potencial químico do componente mais solúvel torna-se quase zero, e a taxa global de dissolução passa a ser controlada pela lenta "sangria" (difusão) do componente menos solúvel.
• Simulações Numéricas:
  • Foram realizadas simulações resolvendo diretamente as equações diferenciais acopladas para validar as aproximações analíticas.
  • Os parâmetros variaram incluem tensão interfacial ($\sigma$), frações molares iniciais, volumes molares ($V_m$) e solubilidades.

3. Contribuições Principais e Resultados

A. Identificação de Três Estágios de Dissolução

O estudo caracteriza a cinética de dissolução em três fases distintas:

1. Pré-lock-in (Pré-bloqueio): Uma fase rápida onde a concentração do componente pouco solúvel aumenta rapidamente para compensar a pressão de Laplace.
2. Lock-in (Bloqueio): Um estado estacionário onde a dissolução é controlada pela difusão do componente menos solúvel. A cinética segue uma lei cúbica clássica (volume vs. tempo), e a composição da gota permanece quase constante.
3. Late lock-in (Bloqueio tardio): No final da vida da partícula, quando o raio se torna muito pequeno, a pressão de Laplace domina novamente sobre o efeito de Raoult, levando a um aumento rápido na concentração do componente pouco solúvel e desvios da linearidade no volume.

B. Números Adimensionais (Lock-in Numbers)

O autor introduz dois números adimensionais críticos para prever o comportamento do sistema:

• $L_1$ (Primeiro Número de Lock-in): Relacionado à magnitude da mudança relativa na concentração do componente menos solúvel ao entrar no estado de lock-in.
  • Se $L_1 \gg 1$, a aproximação de composição constante é válida e a cinética segue a lei cúbica clássica.
  • Se $L_1 \approx 1$ ou menor, ocorrem desvios significativos e a fase de "late lock-in" afeta a vida útil total da partícula.
• $L_2$ (Segundo Número de Lock-in): Relacionado à razão entre as taxas de dissolução dos componentes puros e a eficiência do bloqueio.

C. Nova Equação Empírica Unificada

O principal resultado quantitativo é a proposta de uma equação melhorada para a taxa de dissolução ($w$) que cobre todo o intervalo de composições (de $x_{02}=0$ a $1$) e considera volumes molares diferentes e não idealidades:

$$w = \frac{4\pi D_2 C_2^* V_{m2}}{r} \left( \frac{\phi_1}{\phi_1^* \gamma_1} + \frac{\phi_2}{\phi_2^* \gamma_2} \right)^{-1}$$

(Nota: A equação exata na forma apresentada no texto é a Eq. 31, que combina as taxas de dissolução individuais ponderadas pelas frações de volume e coeficientes de atividade).

• Esta equação generaliza a fórmula anterior do autor (Ref [2]), que falhava em casos de volumes molares assimétricos.
• Simulações mostram que esta nova equação prevê taxas de dissolução com desvios de apenas ~20% mesmo para volumes molares com fatores de assimetria de 8x, enquanto a fórmula antiga poderia ter erros de até 5x.

D. Validação Numérica

• As simulações confirmaram que, para valores altos de $L_1$ (ex: >10), a relação entre o volume da partícula e o tempo é altamente linear (comportamento cúbico), validando a aproximação de estado estacionário.
• Para valores baixos de $L_1$ (ex: <2), observa-se uma curvatura significativa no final da dissolução, indicando que a composição da gota muda substancialmente antes da dissolução completa.

4. Significância e Implicações

• Estabilização de Dispersões: O trabalho fornece uma base teórica robusta para entender como aditivos pouco solúveis (mas não totalmente insolúveis) inibem o Ostwald ripening em emulsões e dispersões coloidais.
• Aplicação em Formulações: A equação desenvolvida permite prever a estabilidade de formulações complexas (como tintas, cosméticos ou fármacos) sem a necessidade de simulações computacionais pesadas, cobrindo desde aditivos traço até fases majoritárias.
• Conexão com o Amadurecimento de Ostwald: O modelo serve como pré-requisito para a próxima etapa da pesquisa do autor, que aplicará esses princípios à teoria de amadurecimento de Ostwald em sistemas com distribuição de tamanhos de partículas. O autor sugere que, para $L_1 \gg 1$, a distribuição de tamanhos segue a teoria clássica de Lifshits-Slezov-Wagner, enquanto para $L_1 \approx 1$, pode surgir uma distribuição bimodal.

Em resumo, este artigo avança significativamente a teoria de estabilidade de gotas bifásicas, fornecendo uma ferramenta matemática unificada que integra efeitos termodinâmicos (Raoult), capilares (Laplace) e difusivos, superando limitações de modelos anteriores que assumiam idealidade ou insolubilidade total.