Emulsion Stability

Quick Start

Those who make emulsions recognise that systematically changing a surfactant blend changes the stability, for reasons that involve either handwaving or incomprehension.

In a tour-de-force, Kiran and Acosta came up with a model that describes emulsion stability in terms of HLD-NAC, starting with the observation that emulsions are least stable when HLD=0.

This expert-level app includes not just the basic theory but the effects of creaming and Ostwald ripening. Take your time; it's heavy going, but rewarding.

Emulsion Stability

Tail length L Å
Rigidity Er kT
T °C
Surf φ
MWt
A Ų
Oil fraction φ
DLVO Å
c mole %
rpm
D cm
ρOil
ρWater
μ cP
Sugg. Er
r nm
Min.Surf.
Sep. Time
rd
γOW
Ionic
Ethoxylate
Extended
NAC Viscosity Correction

The stability of emulsions depends on the initial drop size distribution. The initial drop size depends on the interfacial tension (γOW) and the mixing conditions. The interfacial tension can be calculated using the HLD-NAC model. To calculate the drop size there are different equations that have been proposed in the literature, some that are more appropriate for laminar flow and others for turbulent flow. The equation used in this application is derived from the Hinze-Kolmogorov theory of drop breakup in turbulent flow. The equation to calculate the drop radius (rd) from the Sauter diameter (d32) can set in terms of the density of the continuous phase (ρ), the diameter of the impeller of the mixer (D), and the speed of agitation (rpm):

rd=0.063.D[ρ(rpm/60)²D³/γOW]-0.6

Although this equation has its limitations and others might be chosen, it is adequate for our purposes. Note that in addition to the obvious importance of rpm and impeller diameter, γOW which, of course, is controlled via HLD, is very important.

Drop Size

rpm
D cm
γOW mN/m
ρ g/cc
Rd µm

This equation is valid for systems where the ratio between the viscosity of the dispersed drop and the viscosity of the continuous phase is approximately 1. Systems that have high or low viscosity ratios are more difficult to emulsify and the resulting drop sizes are substantially larger than the values predicted here.

The equations included in this app to calculate emulsion stability assumes that this separation is controlled by coalescence (default case). The coalescence time is calculated based on the collision-coalescence-separation (CCS) separation model of Kiran (Sumit Kiran PhD thesis: Application of the HLD and NAC models to the formation and stability of emulsions, 2013) developed in the Acosta group at U. Toronto.

The CCS model assumes that the drops are small enough to coalesce following the classical model of Davies and Rideal, which considers the probability of coalescence resulting from the Brownian collision of drops. This is the model used in the Coalescence app. The CCS further assumes that after coalescing, the resulting drop quickly separates from the emulsion. Kiran assumes that the activation energy of the coalescence process has two terms, one associated with the energy cost of removing the surfactant from the interface (proportional to the interfacial rigidity, Er), and another term associated with the energy required to form a “neck” connecting the drops. The energy cost associated with the formation of the neck depends on the interfacial tension and the critical separation between the colliding drops. The model calculates this critical separation considering the size of the drops, the Debye length (DLVO slider) and the solvation distance (also introduced through the DLVO slider). In summary, the inputs are:

  • The tail length of the surfactant, L
  • The Rigidity of the surfactant, in units of kT which can be assumed as 1 for a typical ionic and 4-5 for a typical non-ionic and 5-10 for extended surfactants, depending on the tail length. A rigidity estimator is provided
  • Oil fraction
  • Viscosity
  • Energy dissipation rate (from rpm & impeller diameter
  • Temperature (the change in HLD with temperature is not included)
  • For the NAC viscosity correction, Surf φ, MWt and A are required
  • For the min. surf. conc calculation, Surf MWt and A are required along with a chosen minimum drop radius

To build up your understanding of what is going on, you can choose to plot the settling time, the emulsion drop radius or the interfacial tension.

The app introduces a modification to the CCS model by considering collisions driven by settling or creaming, which are more appropriate for larger drops. If the creaming velocity is sufficiently high for drops to coalesce in timescales shorter than that from the Brownian motion (which is very large for large drops) then you see a "creaming" fall-off of the stability. The creaming velocity depends on the density difference between oil and water. If these are set equal then no creaming effects can take place. The app compares the characteristic collision time considering Brownian motion (Davies and Rideal) and settling/creaming. The lowest characteristic time is used as the pre-exponential of the activation energy term.

The app also considers the time to double the volume of a drop (similar to the coalescence of two drops of the same size) according to the Ostwald ripening mechanism. Again, the shortest time between the coalescence model and the Ostwald ripening is presented in the total calculation of stability. The Ostwald time depends on a Diffusion coefficient (assumed to be 1E-06 cm), a molar volume for the oil (assumed to be 200cc/mole), and the solubility of the oil, c in mole %, in the aqueous phase. This final value is the great unknown. It is certainly increased if the surfactant is inefficiently sitting in micelles able to transport the oil. At this level of precision it's probably best simply to see what sort of value of c is needed to cause a significant wing shift, then you can worry about whether those sorts of concentrations are likely to be realistic.

If either creaming or Ostwald effects create a tail then the appropriate term appears on the graph.

The NAC Viscosity Correction option allows you to show the likely effects of the viscosity increase that arises near HLD=0. As shown by Acosta (Sumit K. Kiran and Edgar J. Acosta, Predicting the Morphology and Viscosity of Microemulsions Using the HLD-NAC Model, Ind. Eng. Chem. Res. 2010, 49, 3424–3432), as the microemulsion drop size (not the emulsion drop size!) increases near to HLD=0, viscosity increases via a Doi and Edwards model for the interactions of ellipsoidal nanodrops. The equation breaks down close to HLD=0 so the viscosity is capped at 5* the original viscosity. This increase in viscosity is well-known experimentally though its effects on emulsion stability are less well-known and alarmed experimentalists might have smoothed out some "anomalous" data points.

You need a minimum concentration of surfactant in order to create an emulsion drop of a given radius. Because surface area to volume ratio increases with smaller drops, you need more surfactant for smaller drops. The calculated minimum surfactant concentration C (in g/l) for an oil fraction of φ, drop of radius r (in Å), with surfactant head area A (in Å) and MWt is given by:

C=φ*(3/r)*(1/A)*10^4/6.023*MWt

Given that a typical surfactant will be present at >1% (>10g/l) and that many emulsions have drops >1µm, the amount of surfactant is rarely limiting.

There is one more subtle effect. When the concentration of surfactant plus solubilized oil is less than ~0.4% then the osmotic pressure stabilisation ceases to be significant. As you go down to low surfactant concentrations the stability curve can show strange spikes when the diminishing solubilisation of the oil ceases to compensate for the low level of surfactant.