Granulation and suspension rheology: a unified treatment

Mixing a small amount of liquid into a powder can give rise to dry-looking granules; increasing the amount of liquid eventually produces a flowing suspension. We perform experiments on these phenomena using Spheriglass, an industrially-realistic model powder. Drawing on recent advances in understanding friction-induced shear thickening and jamming in suspensions, we offer a unified description of granulation and suspension rheology. A 'liquid incorporation phase diagram' explains the existence of permanent and transient granules and the increase of granule size with liquid content. Our results point to rheology-based design principles for industrial granulation.

solid-content dispersions are the desired end point for other applications, e.g., ceramic pastes, construction and medical cements, and molten chocolate 5 .If insufficient liquid is added in these systems, flow fails -they fracture and jam, e.g., near constrictions.
Applied research into these two areas has hitherto proceeded separately: where the interest of one community ends ('overwet'), the attention of the other begins.However, in terms of the amount of liquid incorporated, the preparation of granules and high-solid dispersions form a continuum, as is visually apparent in published snapshots taken from the process of 'conching' chocolate 5 (see also Supplementary Information 6 ) and mixing concrete 7 .Thus, it is of interest to enquire whether, on some fundamental level, these two areas of industrial practice may be amenable to a single, unified description, with insights from each enriching the understanding of the other.Indeed, it has been suggested that the shear-induced jamming may be related to granulation 8,9 : the fracturing of an 'underwet dispersion' into granules is the system's response to the impossibility of homogeneous flow.
We report an experimental study of the entire range of phenomena as progressively more liquid is incorporated into a model industrial powder of ∼ 10 µm glass particles.We observe successively the formation of two different kinds of granules, which, at the point of 'overwetting', merge into a high-solid-content dispersion, Fig. 1.We characterise the rheology of these dispersions, which shear thicken, Fig. 2. Interpreting our data in the light of recent advances in suspension rheology [10][11][12][13][14][15] , we construct a 'liquid incorporation phase diagram', Fig. 3 'tie line construction' satisfying the 'lever rule' explains the dependence of granule size on solid content, Fig. 4. Our results confirm theoretical proposals 8,9 that link suspension rheology and wet granulation.
Suspensions of glass spheres in a mixture of glycerol and water at various solid volume fractions, φ, were mixed using different methods to apply 'high' and 'low' stresses (see Methods for details).Note that φ = V solid /(V solid + V liquid ), so that in a dry powder φ = 1.
The result of high and low stress mixing at 0.5 ≤ φ ≤ 0.75 is shown in Fig. 1, where we have color coded different regimes.Samples with φ = 0.75 and 0.70 produced solid granules upon highstress mixing, which remained solid but grew in size after subsequent low-stress mixing.Samples with φ = 0.65 and 0.60 also produced solid granules after high-stress mixing, but these merged to give flowing suspensions after low-stress mixing.Agitating the flowing suspensions under high shear regenerated the solid granules.Samples with φ = 0.55 and 0.50 always produced flowing suspensions.
We next characterised the rheology of the flowing suspensions produced in the range 0.3 φ 0.65.These samples were prepared from granules made in the high shear mixer at φ = 0.65, which were then vortex mixed to a flowing suspension and diluted to the required volume fraction.
These suspensions shear thicken, Fig. 2 (inset): the viscosity increases from a lower to a higher Newtonian plateau as stress increases, with the onset stress for thickening, σ ≈ 1 Pa, being approximately independent of φ.The viscosity at a fixed stress, η(σ; φ), diverges at some φ J (σ).
φ ≥ φ rcp : samples never flow.Inset: η r as a function of σ for various φ., = points plotted using the same symbols in the main graph.
These two volume fractions demarcate three regimes in Fig. 2. At φ < φ m (red) the system flows at all applied stresses, transitioning from a low-to a high-viscosity state as σ increases beyond σ .At φ m ≤ φ < φ rcp (yellow) the system flows at σ < σ , but jams into a solid-like state 9, 17 at higher stresses.Finally, at φ > φ rcp (green), there is inadequate liquid to disperse all the particles, and there is no meaningful 'suspension rheology'.Significantly, the three types of behavior shown in Fig. 1 occur in precisely the three concentration regimes in Fig. 2, pointing to a connection between rheology and liquid incorporation via the physics of jamming 8,9 .
To elucidate this connection, we turn to the measured jamming point as a function of applied stress, φ J (σ), which can be fitted well by a stretched exponential 6 .The inverse function σ J (φ), In order to situate our observations, Fig. 1, on Fig. 3, we need to estimate typical stresses encountered in our protocol.The main mechanical action in our high-shear mixer occurs in the h ≈ 1 mm gap between the tip of the blade (speed v ≈ 3 m s −1 ) and the bowl, at a shear rate of ∼ v/h ∼ 3000 s −1 .Multiplying by the appropriate viscosity turns this into a stress.Below φ m , we use the high-stress plateau value; above φ m we use the high-stress plateau value of the most concentrated sample below φ m as a lower bound.In our vortex mixer, the centrifugal stress acting on a particle is ≈ 0.2 Pa σ 6 .These crude estimates locate the observations shown in Fig. 1 on the jamming phase diagram to within order of magnitude on the stress axis, Fig. 3.It is clear that the jamming phase boundary demarcates granules ( ) and flowing suspensions ( ) in our experiments.Granulation is indeed directly related to jamming.
Cates and others have adumbrated a physical mechanism for this connection 8,9 .A suspension droplet at φ subjected to stress σ > σ J (φ) will be in a jammed, solid state.Trying to shear a jammed state at constant γ generates very large stresses, which we assume exceed σ cap ∼ Σ/a, a capillary stress scale (where Σ is the interfacial tension and a = d/2 is the particle radius).Particles will therefore protrude from the droplet-air interface.If σ cap > σ J (φ), then upon removal of the applied stress, the negatively curved interfaces between protruding particles continue to exert a high enough stress on the droplet to trap it in a jammed, solid state -it is a granule, which is matt because of protruding particles.Note, in passing, that there has long been experimental evidence for a critical Σ/a ratio for granulation 18 .
Such granules are 'fragile' 19 : each offers solid-like resistance only to the particular config-uration of stresses that jammed it in the first place.Even very low stresses in other directions will unjam the state 20 .Thus, a granule with φ < φ rcp subjected to random stresses below σ J (φ) will unjam and remain fluid, because the external stresses are now insufficient to re-jam it.A collection of such droplets will merge into a bulk, flowing suspension.
As in a conventional phase diagram, a 'lever rule' applies to Fig. 3, in both cases reflecting mass conservation.Consider a mixture with average composition φ = φ 0 (point A) consisting initially of dry powder at φ = 1 (point C) and liquid sprayed into the mixer at φ = 0.As mixing progresses, particles are incorporated into droplets (all d) until their composition reach point B, which for σ σ will be at ≈ φ m .Here the droplets jam solid and will not incorporate any more particles, so that some unincorporated dry powder remains.Mass conservation then predicts that the proportion of granules to dry powder is given by the lever-rule construction, viz., AC/AB.
Similarly, granulating a mixture starting at state point D gives granules at φ rcp (point E) coexisting with powder (point F) in proportion DF/DE.
We observed only matt granules but no loose, dry powder up to φ 0.85.The most parsimonious assumption is to postulate granules in which the outmost layer of particles protrude, which can be treated as a layer of dry powder.Based on the lever rule (= mass conservation), we can predict the radius of such granules.Working in more general terms for later use, the volume of solid and liquid in a granule of radius R with a dry powder shell of thickness t s , which in general may be dependent on stress, and φ = 1 (dry = no liquid) and a jammed core (φ = φ J ) are respectively.Since φ = V solid / (V solid + V liquid ), we find The predictions for t s = d and σ J = φ rcp and φ m are shown in Fig. 4 (dot-dash).In both cases, R decreases with φ, starting from a divergence at φ J .The latter is qualitatively consistent with previous reports of diverging granule size 'at 100% liquid saturation' 22,23 .
To test Eq. 4 quantitatively, we measured the volume-weighted mean granule size, R(φ), prepared using a smaller version of the high-shear mixer used for the experiments reported so far (to conserve material) 6 .Our results for granulating at high and low stresses are plotted in Fig. 4 (points).Note that experiments at high stress below φ rcp proved impractical: granules fluidised into a high viscosity suspension in the low-stress regions of our instrument, which required a higher torque to flow than this mixer could generate.Nevertheless, R does indeed increases with decreasing φ, with clear evidence for a divergence for the low-σ data at ≈ φ rcp .The predictions of Eq. 4 capture the form of the experimental data, but are numerically two orders of magnitude too small.
To investigate the cause of this discrepancy we cut open a granule 6 , revealing a wet core and 0.55 0.65 0.75 0.85 φ Figure 4: The volume-weighted mean granule size, R, in units of the mean primary particle radius, a, as a function of the average system composition, φ; red: low stress, black: high stress.Lines: calculated R/a for shell thickness t s = d in Eq. 4 (dot-dashed, multiplied by 100 for comparison to data), and fitting Eq. 4 to return t s = 74d and 56d for high and low stress data sets respectively (dotted).
an essentially dry shell that was many particles thick2 .Such core-shell granule have been observed before 24,25 .Presumably there is a very small amount of liquid in this dry shell forming capillary bridges to bind the particles together 2, 26 ; but the liquid content is too small to visualise, justifying the assumption of φ = 1 in the shell we used in arriving at Eq. 4. The assumption that the jammed core is at φ J is justified by X-ray tomography 6 , which detects 2 % vol.air inside granules.
We now fit Eq. 4 to our data, with the shell thickness t s (σ) as a stress-dependent parameter, Fig. 4 (dotted).This returns shell thicknesses of 74d and 54d for high-and low-σ respectively, possibly suggesting that high-stress mixing makes available somewhat more liquid for capillary bridges to build a thicker dry shell.Fitting our two data sets using a single shell thickness returns t s = 60d, while measurement of a single granule 6 gave ≈ 45d.(Multiple measurements failed because granules typically disintegrated during sectioning.)Performing a similar analysis using a polydisperse granule distribution 6 does not materially change these conclusions.
To summarise, we have shown that the phenomenology of liquid incorporation into powders with repulsive interactions is dominated by shear thickening, which, in extremis, leads to jamming.
Thus, the main control parameter is the stress to which the system is subjected during mixing, because there is a volume-fraction-dependent stress, σ J (φ), above which frictional contacts result in jamming.Mixing below σ J (φ) yields flowing suspensions.Mixing above σ J (φ) yields wet, jammed homogeneous cores coexisting with dry powder.We observed these coexisting states as core and shell respectively in heterogeneous granules.Conservation of mass leads to the prediction of divergent granule sizes as φ → φ J from above, as observed here and in previous work.Quantitatively, a stress-dependent shell size gives good fits to the observed granule size as a function of φ.
Mixing stress has not been identified as a key variable before in granulation: a well-known 'granule regime map' 4 uses the 'maximum pore saturation' and a 'deformation number' as control parameters.The former roughly plays the role of φ −1 , while the latter measures kinetic energy density inside the mixer.These variables control kinetic processes, and affect, e.g., the thickness of the dry 'shell' on core-shell granules.
Our incorporation phase diagram, Fig. 3, suggests novel design principles.For example, the formulation space in which granulation can occur could be tuned by varying φ rcp and φ m through modification of e.g., the particle size or shape distributions.Separately, while the outside of a core-shell granule will always be at φ = 1, the density of the core can be 'tuned' by moving along σ J (φ) by using different mixing stresses.Furthermore, advances in shear-thickening rheology tell us how to modify the σ J (φ) curve itself: e.g., φ m can be tuned by changing the interparticle friction coefficient 10,11 , e.g., by using a variety of surface additives 27 .A full understanding, of course, has to await future work that brings together our incorporation phase diagram and the kinetic regime map 4 .Note that the kinetic energy density axis in the latter has the dimensions of stress, which is an axis in the former.A deeper connection seems probable.
Finally, many systems exhibit more complex rheology such as shear thinning when additional constraints, e.g., adhesion, are considered 28 .The impact of such complexity on the incorporation of liquid into powders is not immediately clear and should be probed in future work.

Methods
We used soda-lime Potters Spheriglass (A-Type 5000, dried at 120 • C for 3 hours), which consists mostly of polydisperse hard spheres (HS; mean diameter d = 7.2 µm, polydispersity = 147%), but with some irregular shards 6 .The powder was dispersed into a 9:1 by volume glycerol-water mixture (viscosity η s = 0.336 Pa s) using a two-stepped protocol at various solid volume fractions φ = V solid /(V solid + V liquid ) (so that in a dry powder φ = 1) calculated from mass fractions 6 .High stress mixing using a bespoke mixer is followed by lower stress mixing using an Ika Vortex Genius 3 mixer.
Our bespoke device consisted of a high-torque overhead mixer (Ika Eurostar Power Control-Visc) driving an aluminium impeller with three equi-spaced blades with a 45 • rake at 500 rpm.

Figure 1 :
Figure 1: The result of mixing at different volume fractions and stresses.(a) Output from the high-shear mixer, and (b) and subsequent low stress vortex mixing.Tube diameter = 20 mm.Colors correspond to the φ regimes in Fig. 2. φ < φ m : samples are flowing, liquid suspensions in both cases.φ m < φ < φ rcp : high-shear mixing produces solid-like granules, which melt to a flowing liquid suspension upon vortex mixing.φ > φ rcp : solid granules from the high-shear mixing remain solid upon vortex mixing.