Steady-state droplet size in montmorillonite stabilised emulsions.

The formation of hexadecane-in-water emulsions stabilised by montmorillonite platelets was studied. In this system the platelets form a monolayer around the droplets and the droplet size decreases with increasing platelet volume fraction. However, the number of platelets present exceeds that required for monolayer coverage. The kinetics of emulsification were investigated and coalescence of droplets during turbulent mixing was found to continue even after the droplets had reached their ultimate size. Non-spherical droplets, resulting from arrested coalescence, were not observed suggesting that particles may be desorbing from the interface during the turbulent flow. A kinetic model based on a competition between droplet break-up and coalescence, mediated by particle adsorption and desorption, reproduces experimental trends in droplet diameter. The model can be used to predict the most eﬃcient formulation to minimise droplet diameters for given materials and mixing conditions and sheds light on the processes occurring during emulsification in this system.


Introduction
Solid stabilised emulsions (also known as Pickering emulsions 1,2 ) are alternatives to surfactant or polymer stabilised emulsions. 3 Adsorption of nano-or microparticles to the interface between two immiscible fluids is energetically favourable when the energies of the interfaces created by adsorption of a particle are lower than those removed in the process. Typical conditions (hydrocarbon oil and particles that are wet by both phases) produce adsorption energies many times the thermal energy resulting in droplets that are very stable to coalescence under quiescent conditions. Particle stabilised emulsions therefore offer advantages over surfactant stabilised emulsions in that they can form a strongly adsorbed solid barrier to coalescence and the particles used are typically non-hazardous.
Many different particles have been found to possess the correct wetting characteristics to stabilise oil/water interfaces. Silica, 4 latex, 5 metal oxide, 6 clay 7,8 and microgel 9,10 particle stabilised emulsions have all been reported. Upon a simple geometric balance between the volume (V) of dispersed phase present and the total area (A) that the particles can occupy on a curved interface, after consideration of packing constraints, a typical droplet diameter of d ¼ 6V A can in principle be achieved if enough interfacial area can be generated and all particles are able to adsorb. Indeed many systems have been reported to achieve this high surface area. [11][12][13] It is also commonplace for there to be well defined trends between droplet size and processing or compositional parameters where this idealised surface coverage is not achieved. Sturzenegger et al. reported the effect of the competition between droplet break-up and coalescence in the formation of zinc oxide nanoand microparticle stabilised water-in-toluene emulsions formed using a rotor-stator mixer. 6 They showed that an increase in particle concentration resulted in a fall in average droplet diameter but there was always an excess of un-adsorbed particles in the continuous phase. The droplet size was therefore not controlled by the maximum area that the particles could cover, but rather by the speed at which the newly formed interface could be stabilised. Cui et al. reported a similar trend of size control with particle concentration for montmorillonite platelet stabilised polydecene-in-water emulsions (also using a rotor stator mixer) where a simple balance of geometry showed that the particles were in a 10Â excess for the size of droplets produced. 8 Ashby and Binks reported that large amounts of LAPONITE s platelets were required to stabilise toluene-in-water emulsions formed using a rotor-stator mixer. 14 Similarly Lagaly et al., 15,16 Garcia et al. 17 and Binks et al. 18 described the use of smectite and hectorite clays as dual emulsion stabilisers and rheology modifiers due to their presence in the continuous phase after emulsification. Additionally Binks and Kirkland showed the presence of particles in the continuous phase of oil-in-water and water-in-oil emulsions after emulsification by scanning electron microscopy. 19 There has been some discussion of the individual processes of break-up 20 and coalescence 21 of particle laden droplets and much discussion of the competition between these processes 6,22,23 but little attempt to quantify the interplay between all processes occurring during turbulent mixing of particle stabilised emulsions. A model where droplets break-up and are coated by particles (and also coalesce somewhat) until the particle surface is covered implies that, given enough time, the minimum droplet size allowed by the mixing conditions would be reached. In the present system and those discussed above this is not found to be the case, therefore there may be some as yet unconsidered factor preventing full surface coverage. Previously reported factors that may affect particle adsorption at the interface include electrostatic barriers, 24,25 mass transport of particles 6,26 and particle aggregation. 14,22 Upon consideration of experimental findings and the energies involved (see ESI, † for calculations) we find that in conditions pertaining to our experiments it is likely that particles may be forced both onto and off the interface during high shear mixing, despite desorption not occurring under quiescent conditions. In a recent paper Phipps and Gittins reach the same conclusions. 27 Moreover, it has recently been shown that desorption of particles from the oil/water interface can occur in high energy environments. 28 We show that montmorillonite platelet stabilised hexadecane-inwater emulsions consist of droplets covered by a monolayer of particles, that droplet size falls with increasing platelet volume fraction and decreasing oil volume fraction and that droplet coalescence continues once the droplets have reached their ultimate size. A kinetic model is constructed and it is shown that experimental trends are best reproduced when the model allows particle desorption during high shear mixing. The success of this model gives the first suggestion that a competition between particle adsorption and desorption is kinetically important for some systems and may be an indication of the position nanometric montmorillonite platelets on the continuum of emulsion stabilisers ranging from small molecule surfactants to microscopic particles.

Materials and methods
Wyoming montmorillonite (SWy-2) was purchased from the Clay Minerals Society source clays repository at Purdue University. The composition of SWy-2 is (Si 7.94 A l0.06 ) (Al 2.88 Fe 0.5 Mg 0.62 )-O 20 (OH) 4 Na 0.68 29 and it has a cation exchange capacity of 84 mequiv./100 g. 30 Berol R648 is an alkyl polyglycol ether ammonium methyl chloride and was kindly donated by Akzo Nobel. Reagent grade NaCl was purchased from Fischer, n-hexadecane (99%) was purchased from Sigma Aldrich. Montmorillonite dispersions were prepared as previously reported. 31 45 g L À1 powdered clay was added to deionised water and stirred for 24 hours. The suspension was then dialysed against aqueous NaCl (1 M) for 1 week changing the solution everyday, to remove unwanted ions such as Ca 2+ , and then dialysed against deionised water, changing water every day, until the conductivity of the dialysate was below 5 mS cm À1 . The dialysed suspension was then diluted to 1 wt% and left for 24 hours for large impurities such as sand and quartz to settle out. The supernatant was concentrated up by evaporation and used as stock. Material parameters relevant to the kinetic model are shown in Table 1.
The emulsions were prepared by first diluting stock montmorillonite suspension with appropriate NaCl and Berol solutions. NaCl concentration was fixed at 0.01 M and Berol added at 25 mg g À1 montmorillonite as this composition was previously found to produce stable Pickering emulsions. 8 Hexadecane was then added to the suspensions at different oil (f o ) and platelet (f p ) volume fractions (calculated from densities shown in Table 1), hand-shaken to pre-emulsify and finally mixed under high shear using an Ultra Turrax IKA T18 with a S18-10G dispersing element at 22 000 rpm for 5 min. Total sample volumes were 10 mL.
Unadsorbed platelets were tested for their emulsifying ability by preparing a dilute emulsion, allowing the droplets to cream, isolating some of the free particles and using them to prepare another emulsion. It was found that particles could stabilise emulsions regardless of their history suggesting that all platelets present possess equal emulsifying ability.
Emulsions were imaged using an Olympus BX51 optical microscope and Pixelink PL-B62CU colour CCD camera using bright field optics with a light balancing daylight filter (BF), differential interference contrast (DIC) and polarising (POM) optics. Droplet size histograms were produced by counting a minimum of 100 droplets from each sample using Linear Intercept software (TU Darmstadt). The surface weighted droplet diameter D 32 is used throughout because the relationship between average droplet diameter and population surface area is used in the analysis of results.
Small angle X-ray scattering (SAXS) measurements were carried out using a SAXSLAB Ganesha 300 XL+ with a Genix 3D Cu-source and Pilatus 300k detector running at 20 Hz. Data was reduced and water background subtracted using SAXSgui software (JJ X-Ray Systems ApS).

Structure of a montmorillonite stabilised emulsion
The organisation of platelets at the oil/water interface was determined by microscopy and small angle X-ray scattering. Fig. 1 shows DIC and POM images of a typical emulsion. A distinctive Maltese cross pattern is observed in POM around the droplet edge where the brightness increases from darkness at the angle of the polariser or analyser to maximum at 451 to the pr p 2 direction of the polariser or analyser (as has been observed before 8 ). Such a pattern was common across all emulsions prepared. This suggests that the platelets are organised parallel to the droplet edge which is consistent with an adsorption energy argument as this configuration would lead to the maximum reduction in oil/water interfacial area per platelet. Control images for a surfactant stabilised emulsion showing no Maltese cross pattern can be found in the ESI. † The possibility of particle stacking was tested using SAXS. Fig. 2 shows scattering of a suspension of Berol R648 treated montmorillonite platelets at 0.01 M NaCl and the corresponding hexadecane-in-water emulsion. Both datasets show a scaling of I(Q) p Q À2 at low Q characteristic of two dimensional plate-like objects. At high Q the emulsion signal begins to deviate from the power law. At wider angles (see ESI †) the emulsion sample shows a peak at 1.46 Å À1 which is too large to be due to platelet stacking and is, in fact, a result of hydrocarbon chain to chain separation. 35 The absence of any correlations in the measured Q range in addition to the patterns observed in polarising optical microscopy show that the platelets form an adsorbed monolayer at the oil-water interface and are fully exfoliated in suspension.

Trends in droplet size
The effect of changing emulsion composition on the final droplet populations was investigated by preparing emulsions at a range of f o and f p values. Fig. 3 shows DIC images of emulsions prepared at f o = 0.47 as a function of f p and Fig. 4 shows associated droplet size histograms. These figures illustrate the effect of increasing the platelet concentration (f p ) on the droplet populations. At lower particle volume fractions the droplet size histograms have a bimodal shape with a significant population of larger droplets resulting from coalescence. As f p increases, bare interface is coated with platelets more quickly, inhibiting coalescence, and the fraction of large droplets decreases, reducing the average droplet size. At lower oil volume fractions the change in this balance is even more notable, with the size distributions being dominated by the smaller diameter population associated with droplet break-up. See ESI, † for histograms of emulsions covering a wider range of emulsion compositions. Fig. 5 summarises the observed relationships between droplet size and emulsion composition for the full range of tested f p and f o (points correspond to experimental data). Clear trends in droplet size with increasing f p and also decreasing f o are observed. This can be understood as an increase in particle attachment rate as f p    Despite such systematic scaling of droplet size with f p , upon a balance of surface area between the platelets present and the droplets produced it is found that of order 10-100Â the droplet surface area is required to produce drops of a given size. Taking the total surface area occupied by platelets at the interface per unit of sample volume as A p = f p /t and total droplet surface area per unit of sample volume as A o = 3f o /% r, where t is the platelet thickness (see Table 1) and % r is the average droplet radius taken as D 32 /2, an area ratio A p /A o can be calculated. Fig. 6 shows the area ratio for all experimental data. It is shown that as f p is increased to push droplet sizes down the emulsification process becomes less efficient and suggests that there may be a kinetic limitation to the full utilisation of platelets in stabilisation of droplet interface.

Coalescence of droplets under high shear
Following the finding that a large excess of platelet surface area is required to drive droplet size down (see Fig. 5 and 6) the emulsification kinetics were investigated. Multiple phenomena that may act to impede particle adsorption at droplet interface have been reported. One is that there is a high energetic barrier to adsorption. This is unlikely to be a factor here as the main driving force for this is electrostatic repulsion and in our emulsions (following previous work 8 ) a high electrolyte concentration was used to screen platelet-platelet and platelet-droplet electrostatic interactions. Additionally it has been reported that for non-spherical particles the adsorption barrier is reduced as particles can make contact with the interface with their smallest dimension, reducing any interaction forces greatly. 36 A second potential reason for the required excess of particles is that they do not adsorb as single particles. This is unlikely here as Fig. 1 and 2 show that the particles are adsorbed as a single monolayer and although the montmorillonite platelets have been reported to form particle clusters in suspension 31,37 it has also been shown in high shear rotational rheometry that the particles behave hydrodynamically as single disk-like particles. 38,39 Finally there is the possibility that particles may desorb during high shear mixing. Upon comparison of turbulent kinetic energy in a typical laboratory scale mixer and the adsorption energy of a montmorillonite platelet it is apparent that they may desorb as a result of such energetic hydrodynamic conditions. One consequence of this would be that the droplet size would reach a steady state, if particles are both adsorbing and desorbing then droplets should also continue to coalesce and break-up as long as mixing continues. A previously prepared emulsion that is subjected to high shear mixing again would therefore undergo further particle adsorption/ desorption and droplet break-up/coalescence but the steady state droplet size would not change. Fig. 7 shows the result of an experiment where two sets of two emulsions were prepared. Each set contained emulsions of identical composition except that in one the dispersed phase was dyed with Sudan II. These two droplet populations were first combined (bottom left of Fig. 7) and then subjected to the same high shear mixing as when they were initially emulsified (bottom right of Fig. 7).
Upon simple combination the dyed and undyed droplets do not mix. Upon high shear mixing the dye is redistributed across the two populations suggesting that further coalescence and break-up of the droplets takes place during high shear mixing resulting in mixing of the oil phases.  If coalescence were to cease then it would be expected that, in the absence of any significant impediment to adsorption, droplets should either be as small as can be produced by the mixer if particles are in excess or as small as the particles can cover if they are not. Neither of these cases are observed in the present system further suggesting that there must be some kinetic limitation preventing particles remaining at the interface. The next section explores a kinetic model describing the competition between droplet break-up and coalescence being mediated by a competition between particle adsorption and desorption during turbulent flow. Predictions of this model are then compared to experimental data in an attempt to shed light on the mechanisms at play during high shear emulsification.

Emulsification in turbulent flow
High shear mixers, such as the Ultra Turrax rotor-stator used in this study, produce high fluid velocities in small volumes. The type of flow can be characterised by the ratio of inertial to viscous forces quantified by the Reynolds number Re ¼ rvL Z , where r, v and Z are the density, velocity and dynamic viscosity of the fluid respectively and L is a characteristic length scale (in this case 7.5 mm which is the diameter of the rotor). In high shear mixers Re is typically 410 3 indicating turbulent flow where the local fluid velocity differs from the time average and turbulent eddies are present. 40 If it is assumed that the turbulence is isotropic (which is valid at small length scales such as that of emulsion droplets) then the theory of Kolmogorov can be used to determine kinetic parameters. 40 The key quantity in analysing a turbulent flow regime is the energy dissipation rate (e, with units W kg À1 ) which has been estimated for toothed rotor-stator mixers 41,42 to be of order 10 5 W kg À1 . From this the Kolmogorov microscale l can be defined which is the smallest size of eddies present (see ESI †). When droplets are larger than l break-up is due to hydrodynamic interactions between the droplet and the eddies. The resulting local pressure fluctuations around a droplet, caused by variations in local flow velocity, that are larger than the Laplace pressure will cause the droplet to rupture; this is referred to as break-up in the turbulent inertial regime. When droplets are smaller than l the local flow near the drop is likely to be laminar. The stress acting on the droplet is therefore a function of the fluid viscosity and it is this stress that causes the drops to deform and eventually rupture. This viscosity dependence is referred to as the turbulent viscous regime. Due to the inverse relationship between Laplace pressure and droplet radius a critical length scale for both regimes results where the kinetic energy of the flow is no longer sufficient to break up the droplets and this can be used as an upper limit to the lowest attainable droplet size. 40 The two regimes are can be identified using a droplet Reynolds number (Re d ) where the droplet diameter is used in place of L; Re d 4 1 indicates turbulent inertial and Re d o 1 indicates turbulent viscous. The parameters for the system studied here are shown in Table 2. The bulk and droplet Reynolds numbers (Re and Re d ) are larger than 10 3 and 1 respectively therefore break-up can be described by the turbulent inertial regime. 40

Kinetics of emulsification
Understanding the emulsification in turbulent flow is key to understanding processing of these systems for industrial use and thus there have been many detailed theoretical studies into the kinetics of the processes involved and the resulting droplet size distributions. [43][44][45][46][47][48][49][50][51] Here we take only the minimum number of arguments required to realistically represent our system and only focus on the ultimate average droplet size. Such an approach allows for the straightforward construction of models  representing multiple phenomena for comparison with experimental results. Two parameters were followed; n d the droplet number density and A c the area of droplet surface covered with particles per unit volume. The droplet number density increases upon a break-up event and decreases upon a coalescence event and can be related to the average droplet size by eqn (1): where f o is the dispersed phase volume fraction. A c increases upon an adsorption event and decreases upon a coalescence event. As A c approaches the total droplet surface area (A total d = 4p% r 2 n d ) the probability of two regions of bare droplet interface coming into contact falls and therefore the coalescence rate drops.
One interpretation of such quantities (when there is an excess of particles) is that the collisions between particles and droplets will occur frequently so platelets adsorb, increasing A c and decreasing the coalescence rate, until A c approaches A total d , at which point coalescence will stop and the final droplet size will be reached. This suggests that the droplet size will always tend to the smallest that can be produced by the mixer and that, given sufficient emulsification time, coalescence will stop entirely. Neither of these phenomena are observed in our system. The kinetic model considered here includes particle desorption in the flow field therefore allowing A c o A total d even when the platelets are in excess. Upon consideration of the turbulent kinetic and platelet adsorption energies involved this would be a reasonable assumption (see ESI, † for calculations).
Combining the relevant processes and their effects on n d and A c gives rise to the following set of equations: Often R detach is ignored, however considering the scaling of particle desorption energy being with the square of the particle size (see calculations in ESI †) a continuum from small molecules surfactants to large solid particles is anticipated. Low molecular weight surfactants can displaced by thermal fluctuations, small particles can be displaced by the energy produced by a typical laboratory scale mixer and large particles are so strongly adsorbed that even this high power mixing cannot displace them. Following the predictions of particle coagulation models [43][44][45]51 we arrive at expressions for the rates R coalesce and R attach given by eqn (4) and (5). The platelets are treated as single particles as SAXS data (Fig. 2) shows that there is no particle stacking and previous rheological studies of montmorillonite suspensions at high shear rates in laminar flow find that they behave hydrodynamically as single disk-like particles. 38,39 Following Delichatsios and Probstein 45 R coalesce and R attach can be expressed as: Here A collision is the collision cross section between colliding bodies. For the purposes of this model the platelets are assumed to occupy an effective swept out spherical volume with radius equal to that of the platelet (r p in Table 1). The efficiency parameter a is included to represent an energetic barrier to platelet adsorption. 24,25,52 In practice, when fitting the model to experimental data, it appeared that a and k detach were correlated adjustable parameters. In what follows, a was set at 0.5 and the detachment rate constant used as an adjustable parameter. Interesting further work would include defining an expression for k detach that illustrates the mechanism of desorption, allowing a and k detach to be decoupled. Since the exact mechanisms of desorption are currently not fully understood this was not attempted here. n p is the particle number density, a p is the area a single particle occupies at the interface and f = A c /A total d is the fraction of droplet interface covered with platelets. u r is the relative velocity between the two colliding bodies. In isotropic turbulence when the colliding bodies are larger than l this is given by u r E 1.37(ed) 1/3 where d = 2% r. 45 The potential for bridging (where two droplets share interfacial particles) was ignored here as due to the extremely small platelet thickness (B1 nm) it follows that for droplets to be close enough to share a platelet the it is likely that the fluid menisci from the droplets will come into contact triggering a coalescence event. 21,53 The break-up of emulsion droplets has been studied for both laminar 54 and turbulent flow 48 and recent simulations have shown the effects of particles at the droplet surface on break-up. 20 The model for droplet break-up in turbulent flow used here follows the work of Coulaloglou and Tavlarides who applied an activation energy approach to represent the transfer of kinetic energy from turbulent eddies to a droplet upon a collision and this has been shown to agree with experimental data. 26,46 The break-up rate is defined as the fraction of droplets with energy high enough to break multiplied by the reciprocal of the time taken for a droplet to break: where E surface is proportional to g ow d 2 and E kinetic is proportional to r d d 3 u r 2 . g ow is the interfacial tension between the oil and water, r d is the density of the dispersed phase and u r is the relative velocity between two points of fluid as defined above. t break can be defined as the time taken for two points of fluid to move the distance of the droplet diameter away from each other and is proportional to d 2/3 e 1/3 . In applying these definitions to eqn (6) the proportionality constants are grouped into C 1 and C 2 in eqn (7) and in the following calculations are allowed to vary to within an order of magnitude of unity, which was found previously to be a suitable approximation. 46 The individual rates are expressed in terms of known parameters as shown in Table 1. These can be substituted into eqn (2) and (3) giving a set of equations that can be solved for n d and A c and, given f o and using eqn (1), the evolution of average droplet size over time can be evaluated.
R coalesce = 3.45pr 7/3 e 1/3 n d 2 (1 À f ) R attach = a1.37p(r + r p ) 7/3 e 1/3 n d n p a p (1 À f ) (9) The model was solved numerically and predicted a fall of droplet diameter to a steady state for parameters corresponding to experimental samples. k detach , C 1 and C 2 were varied to fit the model to the entire ensemble of experimental droplet diameters using a least squares approach giving values of 1.7 Â 10 5 s À1 , 1.2 and 1.2 respectively. Fig. 5 shows a comparison between experimental droplet diameters (points) and the steady state droplet diameters predicted by the model (solid lines). Predictions of the model for identical conditions where k detach = 0 are also shown (dotted lines).
In reality it is likely that k detach would vary with system parameters such as e and the size of the stabilising particles. Neither of these parameters vary in the present experiments, therefore k detach is fixed for all compositions and represents the number of particles desorbing per second in the given flow conditions. At high f p , R attach is initially fast which rapidly increases f and decreases R coalesce . This results in a droplet size that is determined by break-up processes and therefore constant with varying f p .
This model predicts the behaviour in the high shear region (the volume swept out by the rotor 55 ) and not the entire sample volume. The time any given fluid element spends in the high shear region must therefore be compared to the characteristic time scales of droplet break-up, particle attachment and particle detachment. The residence time in the high shear region can be estimated as: Where L is a characteristic dimension and e is the energy dissipation rate. Using L = 7.5 mm (the rotor diameter) and e = 6.45 Â 10 5 m 2 s À3 , t res E 450 ms.
The typical break-up time for a single droplet can be calculated by taking the inverse of eqn (7) when applied to a single droplet which accounts for the time the droplet takes to deform and the probability that it will collide with a turbulent eddy of sufficient energy to cause rupture. Inputting values as above the break-up time scales are estimated to be of order 10 ms.
The typical time for a particle desorption event can be estimated as simply the inverse of the estimated desorption rate constant giving a time scale of 6 ms. Since all of the systems studied reach a steady state the typical steady state attachment time scale must also be approximately 6 ms.
It follows that all of the processes described by the model occur quickly enough that a fluid element entering the mixing region will remain there for sufficient time for droplets to break-up and particles to attach and detach. This means that the model should require no adjustment to account for the probability of processes occurring in the mixing volume.
At low f p the effects of a non-zero R detach become apparent. Fig. 5 shows a steep increase of droplet size with decreasing f p for both experimental data and the model where k detach is nonzero. R attach is lower at low f p resulting in a lower f and therefore higher R coalesce . This shift in the competition between break-up and coalescence results in a larger steady state droplet size at low f p and is in agreement with data. Interestingly it is k detach that determines the magnitude of the deviation from the minimum droplet size. Fig. 8 shows the steady state droplet diameters predicted by the model for f o = 0.4 and varying f p at different k detach and shows that the magnitude of the desorption rate constant determines the f p at which the droplet size hits the high f p plateau. Fig. 6 shows the prediction of the model in the form of an area ratio between the platelets present in the system and steady state droplet size. There is good agreement with experimental data showing that the largest droplet surface area stabilised per total particle surface area occurs at lower f p . Interestingly the model area ratio appears to go through a minimum for f o = 0.36 and 0.47. This suggests that there exist optimum compositions for the most efficient use of stabiliser (largest droplet surface area stabilised per total surface area of particles), however given the small number of experimental data points in this region this is only a speculative conclusion.
The model reproduces the general experimental trends well (as shown in Fig. 5) however is still some quantitative discrepancy between model and data. Such agreement is still encouraging as this model combines a small number of arguments and produces a realistic representation of the experimental system and is a first step towards incorporating particle desorption during mixing into the general understanding of solid stabilised emulsions where it is appropriate. Possible advancements to the model include explicit treatment of adsorption barriers, more detailed treatment of coalescence including film thinning and particle ejection 21 and study into the precise mechanisms of particle desorption under high shear mixing.
The comparison of the model and the experimental data show that to gain real insight into the mechanisms of formation of solid stabilised emulsions many factors must be considered. Here we show evidence that particle desorption during shearing may be an important factor for montmorillonite stabilised emulsions.

Conclusions
It has been shown that a kinetic model based on a balance between droplet break-up and coalescence mediated by particle adsorption and desorption can reproduce experimental trends in droplet diameter of emulsions formed and stabilised in turbulent flow. Under circumstances where these processes occur on similar time-scales a steady state determines the ultimate droplet size as opposed to a regime where the droplet surface area is determined by the particle surface area. 11 This explains the need for a large excess of stabilising particles to drive the droplet size down and is in line with experimental work in this (Fig. 5) and previous studies. 6,8,19 Further work could see this model applied to other regimes such as ultrasound emulsification, which can produce higher e, 56 and emulsions with different dispersed phases. Also refinements to the model to truly capture the approach of a platelet to the interface accounting for orientation, electrostatic repulsion and an investigation into the mechanisms of particle desorption would make interesting further work.
Additionally this could be combined with previous work on rheology control using montmorillonite platelets 16,18,31 to produce a system where the minority component acts as a dual emulsion stabiliser and rheology modifier. This improvement in the understanding of platelet stabilised droplets in turbulent flow can be used to optimise the design of tightly controlled and efficient emulsion-based formulations for example Fig. 5 and 6 shows that the smallest droplet sizes are produced at high f p or low f o but the lowest area ratios at low f p and high f o . It is only by combining insight gained from studies on many different systems that the intricacies of Pickering emulsification will be fully uncovered.