Cohesion-driven mixing and segregation of dry granular media

Granular segregation is a common, yet still puzzling, phenomenon encountered in many natural and engineering processes. Here, we experimentally investigate the effect of particles cohesion on segregation in dry monodisperse and bidisperse systems using a rotating drum mixer. Chemical silanization, glass surface functionalization via a Silane coupling agent, is used to produce cohesive dry glass particles. The cohesive force between the particles is controlled by varying the reaction duration of the silanization process, and is measured using an in-house device specifically designed for this study. The effects of the cohesive force on flow and segregation are then explored and discussed. For monosized particulate systems, while cohesionless particles perfectly mix when tumbled, highly cohesive particles segregate. For bidisperse mixtures of particles, an adequate cohesion-tuning reduces segregation and enhances mixing. Based on these results, a simple scheme is proposed to describe the system’s mixing behaviour with important implications for the control of segregation or mixing in particulate industrial processes.

Mixing of granular materials is important in several particulate processes such as concrete preparation, chemical formulation and pharmaceutical engineering 1,2 . Under the presence of shear, granular materials in these processes often segregate owing to differences in particle properties such as the size and density, causing them to sink or rise within the bed and to eventually accumulate, leading to a variation of the bulk composition and density 3,4 . This segregation behavior can be a difficult problem in some processes that may critically degrade the quality of the final product (e.g., solid dosage forms preparation and food processes), and, in some other cases, can be used for benefit (e.g., recycling and separation processes). Thus, it is important to understand particulate segregation mechanism to efficiently prevent, trigger, enhance or reduce it.
Among the several size segregation mechanisms 5 , the kinetic sieving is the dominant mechanism causing segregation in gravity-driven free-surface flows 6,7 , where the small particles move effectively, shifting downwards, as they are more likely to move into the voids between larger particles. This mechanism performs well when the volume fractions of large and small particles are equal 8 . Using numerical simulations, Hong et al. 9 considered the interplay between size and mass in a binary mixture of spherical particles, and predicted the existence of the reversed Brazil Nut effect (RBN). Khakhar and co-workers 10 showed that the segregation rate of dry particles is independent of the filling level of the rotating drum. Alonso et al. 11 investigated the segregation of particles due to differences in size and density in a two-dimensional horizontal rotating cylinder, and proposed an expression of segregation index to predict if a large particle tends to float or sink in a bed of fine particles. More recently, Hongyi et al. 12 applied unsteady flows to strongly segregating granular materials in an attempt to control the segregation pattern and enhance mixing.
A typical approach to prevent size segregation is by making particulate components sizes as close as possible 13 , or by using cohesive particles 14 . However, cohesion between particles is difficult to introduce and control. As a consequence, several studies regard particles as non cohesive [15][16][17] , but cohesion is always present in most industrial particulate processes, and thus has to be taken into account. In some other studies, cohesion was introduced using various viscous liquids with different capillary forces 14,[18][19][20][21] . Jarray et al. 19 investigated the effect of liquid induced cohesion on granular flow in a rotating drum and showed that capillary cohesion increases the angle of the slope and decreases the granular temperature at the free flowing surface. Shi et al. 22 performed experiments on various size fractions of cohesive limestone powders and showed that cohesion dominates the flow behaviour for fine particles. Li and McCarthy 14 showed that segregation in granular systems can be controlled by adding moisture. Roy et al. 23 studied the effect of wet cohesion on the granular flow and proposed a rheological model that predicts various features such as the shear thinning behavior. Chou et al. 24 demonstrated that the segregation index in wet granular materials decreases with an increase of the angle of the slope in a rotating drum, regardless of the volume or viscosity of the added liquid. Apart from studies on segregation in wet systems with cohesion due to capillary forces, there are less experimental investigations focusing on dry cohesive systems, especially on studies where the cohesive force can be adjusted.
One way to modify the glass surface and obtain dry cohesive particles is the so-called silanization. The original purpose of chemical silanization is to change the hydrophobicity of glass surface 25 . Here, we use extensive silanization to modify the contact cohesion of millimetric size glass particles in a controlled way, with the goal to link this micro-cohesion change to the macroscopic flow and mixing behavior. In industry and within the scientific community, several geometries are used for particulates mixing 1,4,26,27 . In this study, we employ a rotating drum to study the flow and segregation of cohesive and non-cohesive granular materials. This apparatus is extensively used in experiments and as a model system due to its simple geometry comparing to other mixing devices.
To our knowledge, this is the first paper that experimentally investigates the effect of adjustable cohesive forces on the segregation of dry granular systems in a rotary shearing device. We characterize the cohesive particles properties using different experimental approaches, which include Atomic Force Microscopy (AFM) for surface scanning of the particles, bulk flowability in the drum, and cohesion force measurement using an in-house setup. The effects of silanization reaction duration on the surface properties of the particles are presented, and the cohesive and adhesive forces between particles are measured and discussed. Then, cohesive particles are mixed with non-cohesive particles in a rotating drum mixer, and the segregation is investigated under different combinations of cohesive forces and particle sizes, concluding that the cohesion-dependency of segregation may be manipulated to mitigate or enhance particulate mixing.

Results and Discussions
characterization of the surface and cohesive properties of the particles. Typically, to investigate the effect of cohesion on particulate segregation and flow, researchers use capillary forces [28][29][30] . When particles are in contact with each other in the presence of a small amount of liquid, the interstitial liquid forms a liquid bridge between the particle and the cohesion is dominated by capillary forces (see Fig. 1(b)). This approach has several limitations. For instance, the liquid distribution, and hence capillary force distribution, within the particulate system is not uniform especially under dynamic conditions 31,32 . Additionally, tuning the capillary force is not straightforward as it depends on both liquid and particulate properties. To avoid these complications, we use silanization where cohesiveness is only due to immediate chemical bonding with no formation of capillary bridges. In this case, cohesion between particles is only upon contact. If the particles do not touch, the cohesion force is equal to zero due to the very short range of the attractive force. This can be seen in Fig. 1(c,d), where silanized particles are perfectly dry and capillary bridges are absent, making their cohesive behaviour as close as possible to that of real dry cohesive powders. Figure 1(e) shows the cohesive (F AA ) and adhesive forces (F AB ) between two particles of radius 0.85 mm as a function of silanization reaction time. Here, "A" particles are cohesive and "B" particles are non-cohesive. The cohesion and adhesive forces increase with silanization reaction duration. The adhesive forces between non-cohesive and cohesive particles are order of half of the cohesive force between two cohesive particles. Highly cohesive particles of about 0.11 mN are obtained when the silanization reaction duration is longest (i.e., all the Heptane in the Silane solution evaporates during the silanization process). This increased cohesive force can be attributed to the formation of more covalent bonds (-Si-O-Si-) as the reaction progresses. Figure 2 shows microscopic and AFM measurements of dry and silane treated glass particles. Cohesive silanized particles seem to have smoother surfaces with less irregularities than non-silanized particles. Less irregularities infers a more uniform coating, and hence is an important factor in determining the cohesive force between particles. Since the thickness of the Silane coating is in the range of 30 nm, it is assumed that the main elastic properties of the individual particles are still the same after silanization. In addition to direct cohesive measurements, the cohesive force can be correlated with the angle of the flow (i.e., dynamic angle of repose) at the bulk level in a rotating drum. We show in Fig. 3(a) the smoothed angle of the flow θ in a rotating drum as a function of time for non-cohesive and cohesive particles. The particles are of radius 0.85 mm and the rotation speed of the drum is 25 rpm (i.e., Froude number Fr = 0.21). As the drum rotates, particles are lifted to the upper part of the bed and the angle of the slope increases until it reaches a maximum, then avalanches start to occur. For dry particles, successive avalanches happen with large amplitude followed by a relatively steady flow with small avalanche amplitude variations. For highly cohesive particles, successive avalanches remain as the drum rotates. This is because cohesive forces dominate the flow and particles become closely packed and start to flow as a bulk. For medium cohesive particles, the amplitude of the avalanches is the lowest indicating a liquid-like flow behavior. When the cohesion is low, the dynamic angle of repose is close to that of the non-cohesive case. Figure 3(b) shows the dynamic angle of repose averaged over the last 20 seconds versus the cohesion force. The dynamic angle of repose increases with the cohesive force, qualitatively confirming the results obtained by the microbalance setup in Fig. 2. By analogy to the flow of wet particles in the presence of capillary forces 18,19 , we can infer that the increase of the dynamic angle of repose is because the cohesive forces prevent the rolling and cascading of individual particles in favor of bulk sliding. These experiments show that the cohesive force between particles can be adjusted using the silanization reaction duration, and can also be characterized at the bulk level using the angle of the flow in a rotating drum.  conduct a set of experiments with different particle sizes and different cohesive forces as shown in Table 1. In all cases, the drum contains a mixture of 50-50% by weight (w/w) (i.e., 50% of particles "A" and 50% of particles "B") glass particles with the same density of 2500 kg/m 3 .
We plot in Fig. 4(a,b) the mixing index for monosized particles, of radius r A = 0.85 mm, as a function of time for cohesionless (case 1) and cohesive particles (case 4), respectively. The mixing index curves are smoothed using bisquare weighting smoothing to reduce the fluctuation of the mixing index. The raw values obtained directly from image post-processing, before smoothing, are shown in dots in Fig. 4 and in the Supplementary Information file, Appendix A. A Mixing Index (MI) close to 1 means the system is mixed, while an MI value close to 0 indicates a segregated system. The mixing index in Fig. 4(a) increases at first and then, after few rotation of the drum, remains relatively constant at around MI = 0.9, indicating good mixing of the system. However, as shown in Fig. 4(b), when the transparent particles are cohesive, a lower MI value around 0.6 is obtained and segregation occurs between monosized black and transparent particles. This can be explained by the clustering of the transparent particles by cohesion, which pushes them outside and keeps the black particles within the core region of the bed, where the flow is quasistatic. By performing Discrete Element Method (DEM) simulations of a mixture of cohesive and non-cohesive monosized particles in a rotating drum, Yazdani and Hashemabadi 33 observed similar behaviour, where the granular system tended to segregate and cohesive particles moved towards the outer layer of the rotating granular media.  www.nature.com/scientificreports www.nature.com/scientificreports/ We show in Fig. 4(c) the mixing index (MI) versus time for a bidesperse system composed of non-cohesive transparent particles of radius r A = 0.85 mm mixed with non-cohesive red particles of radius r B = 1.25 mm (case 5). In contrast with the previous case, when the system is composed of bidisperse particles, after few rotations of the drum, the mixing index reaches a steady state at about 0.5, indicating segregation, where larger particles are found in the upper region of the bed and near the wall of the drum, and smaller particles are concentrated in the core of the bed surrounded by the large particles. When the same experiment is performed with cohesive particles, the results are different. Figure 4(d) shows the mixing index (MI) of case 8, where particles of radius r A = 0.85 mm are cohesive and mixed with non-cohesive larger particles of radius r B = 1.25 mm. The value of the mixing index in this case is higher than for case 5, indicating better mixing with fewer large particles found in the outer region and more of them in the core of the bed. This demonstrates the ability of cohesion to improve dry particulate mixing. The cohesive force clumps the small particles together, decreasing their probability of moving into gaps between larger particles, which is similar to the explanation given by Li and McCarthy 14 for the case of cohesion due to moisture between the particles.
We report in Fig. 5 three plots of the mixing index measurements for all the cases described in Table 1 as a function of time. In Fig. 5(a), we show cases 1 to 4, which correspond to systems where the particles have the same size. In all four cases, the mixing index reaches a steady state after approximately 7 revolutions of the drum (i.e., 20 seconds). The final mixing index value decreases as the cohesive force increases, confirming that cohesion reduces mixing. We also notice that the higher the cohesiveness, the more time is needed by the system to reach the steady state. Figure 5(b) shows the mixing index for a system composed of particles of radii r A = 0.85 mm and r B = 1.25 mm for different cohesive forces (cases 5 to 8). Again, after few rotations of the drum, the mixing index curves reach a plateau. When the cohesion is low (case 6), MI reaches slightly higher values than in case 5, where particles are not cohesive. As the cohesive force increases, segregation reduces, indicating that cohesion plays a key role in improving the mixing of particulate systems through the formation of clusters. We infer that cohesion makes small cohesive particles to clump together impeding them to pass through the gaps between the large ones. By performing DEM simulations, Aarons et al. 34 arrived to the same statement. They found that the tendency for the small particles to form agglomerates and their ability to move between the big particles are the main mechanisms that control cohesive particulate mixing. Since the cohesive small particles ability to mitigate segregation is dependant on the size of the large particles, we carried out the same experiments but with larger non-cohesive particles of radius r = 2 mm for different cohesive forces (cases 9 to 12). The mixing index for these cases is shown in Fig. 5(c). Similarly to Fig. 5(a,b), a plateau is observed after 20 seconds of drum rotation. For the non-cohesive system (case 9), the value of the mixing index at the plateau is about 0.48, slightly lower than the value obtained in case 5 (i.e., non-cohesive combination of 1.25 mm and 0.85 mm), which is in accordance with the findings of previous work 24, [35][36][37] , stating that a greater difference in size of cohesionless particles increases the degree of www.nature.com/scientificreports www.nature.com/scientificreports/ segregation. After 7 rotations of the drum, low and medium cohesive particles (cases 10 and 11, respectively) reach almost the same level of mixing as the cohesionless case (i.e., case 9), indicating that the cohesion is not strong enough to prevent 2 mm particles segregation. Only highly cohesive particles (case 12) are able to reduce 2 mm particles segregation. All these results show that the mobility of the non-cohesive particles depends strongly on the cohesive force between the particles, which, more importantly, can be tuned to control segregation.
Next, we present diagrams to visually identify zones where systems tend to mix or segregate. We plot in Fig. 6(a) the mixing index resulting from cohesion, Δ = − − MI MI MI coh n on coh , averaged over the last 10 seconds of the duration of the experiment, as a function of the cohesive work change (i.e., surface energy change per particle) when the system goes from a segregated to a mixed state Δ = − W W W m s , with MI coh and MI non-coh are the averaged mixing index for the cohesive case and the non-cohesive case, respectively, and W m and W s are the cohesive work per particle of the perfectly mixed system and the segregated system, respectively. When the mixture is bidisperse (circle symbols in Fig. 6(a)), the mixing degree increases with ΔW, and the cohesion of the particles provides work (i.e., surface energy) that reduces the segregation of the system. ΔW = 0 mN/m means either W s = W m , or all particles are non-cohesive. In case W s = W m , the adhesive work exerted by the cohesive small particles on the large particles is counterbalanced by the cohesive work between the small particles. Either way, the system will segregate by size. This suggests that, as long as ΔW is sufficiently large, the mixing of the bidisperse system is enhanced. This is inline with the work of Chaudhuri et al. 38 , who examined computationally the effect of adhesive forces in a binary system in the absence of cohesive forces. They found that adhesion favors the mixing process. However, we expect that further increasing ΔW will reach a point where the system becomes perfectly mixed. Beyond this point, cohesion will start reducing the mixing of the system, because relatively large clumps of small particles, larger than the large individual particles, will form and migrate to the outer region of the bed by size segregation. This was emphasized by Sarkar and Wassgren 39 , who stated that further increasing the cohesion forces may start hindering the mixing. On the contrary, for monodisperse particles (square symbols in Fig. 6(a)), the mixing decreases with ΔW and cohesion induces segregation. For ΔW = 0 mN/m -i.e., the segregation and mixing work are balanced, the monodisperse system is well mixed.
To quantify the importance of cohesion in comparison to the weight of the particles, we plot in Fig. 6(b) the mixing index due to cohesion, Δ = − − MI MI MI coh n on coh , as a function of the granular Bond number Bo g . The granular Bond number refers to the ratio of the maximum cohesive force, max F F ( , ) AA AB , to the effective weight of the particles:  Figure 6(b) shows that for bidisperse systems (circle symbols), the mixing increases with Bo g . When Bo g > 1, the maximum cohesive force is greater than the effective weight of the particles, and mixing occurs. When Bo g < 1, the interparticle cohesive forces are weak comparing to the effective weight, and thus cannot prevent segregation. We expect that when the bond number reaches a value higher than a critical Bond number Bo g c , the bidisperse system will start segregating because the size of the clustered particles will exceed that of the large particles. We define Bo g c as the bond number where the maximum cohesive force is equal to the weight of the large particles, which gives 3 (see Supplementary file for more details, Appendix B). Unlike the bidisperse case, for monosized particles (square symbols in Fig. 6), segregation is enhanced when Bo g increases due to cohesion. Here, Bo g < 1 indicates mixing. This is because cohesive forces are not strong enough to counterbalance the equally weighted particles of the perfectly mixed system. This counterbalancing occurs when Bo g > 1, causing the system to start segregating due to small glass particles agglomeration. Notice that in this case, = Bo 1 g c , which also indicates that cohesion will reduce mixing. This is supported by the numerical simulations of Sarkar and Wassgren 39 , who showed that cohesive interactions between monosized particles promote agglomerate formation, which decreases mixing. Similarly, Halidan et al. 40 showed numerically that monosized particles mixing is better at moderate and low interparticle cohesion because the cohesive bonds between particles are easy to break. www.nature.com/scientificreports www.nature.com/scientificreports/ Comparing the results given by the cohesive work change approach (Fig. 6(a)) and those given by the Bond number approach (Fig. 6(b)), both give similar results, but the former seems to be more accurate in terms of predicting the mixing degree of the system.

Summary
Using controlled chemical silanization, it was possible to generate glass particles surfaces with adjustable cohesion. The cohesive force between particles was measured using an in-house experimental device and also characterized on the bulk level using the dynamic angle of repose in a rotating drum. We determined experimentally that an increase in cohesive forces increases the dynamic angle of repose, and the smooth flow of cohesionless particles transforms into an irregular flow of a bulk of clumped particles. Then, we investigated the effect of cohesive forces on the degree of segregation in the rotary drum. While strong enough cohesion of small particles improves mixing for bidisperse systems, it causes segregation for monosized particles. Therefore, one should be able to switch between size segregation or mixing of monosized particles simply by tuning on or off the cohesive forces between the particles. For bidisperse systems, avoiding segregation requires a balance between size segregation and cohesive clustering of the small particles. Our experiments illustrate the importance of cohesion for the controlled minimization of particle segregation during large-scale processing of powdery materials, capturing the mixing dependence on particle size and adhesive/cohesive forces between particles. The developed state diagrams can be extended to other size ratios, and may serve as a useful tool for controlling particle mixing in other devices or processes (e.g., slow chute flow, slowly vibrating bed) where shearing is not so strongly localized. Additionally, it would be of great interest to investigate whether similar mixing behaviour is obtained when the large particles are cohesive. Also, it would be important to extend our analysis to other weight fractions of large and fine particles.

Methods
Glass surface pretreatment: Silanization process. Silanization is based on the adsorption, self-assembly and covalent binding of Silane molecules onto the surface of glass particles. Silanes are coupling agents that can interact with both organic and inorganic materials. The silanization process has been used by many scientists for the surface hydrophobization of glass and other materials 19 . It was also applied to fine silica particles as rheological additives for adhesives, resins and paints [41][42][43] .
Chemical compounds used for silanization are: silanization solution 5% (V/V, 5% in volume of Dimethyldichlorosilane in Heptane, Selectophore), Hydrochloric acid (HCl, 0.1 mol), Acetone and Ethanol. The procedure for making the glass particles cohesive is as follows: First, to ensure that the surface of Silica glass particles is free of contamination, they were cleaned for at least one hour by immersion into freshly prepared HCl solution under agitation using a rotor-stator homogenizer. Then, they were rinsed thoroughly with deionized water and oven dried for 3 hours at 60 °C. Afterwards, the freshly cleaned samples of glass particles were immersed in the silanization solution under low agitation speed at room temperature for various durations of time (30, 60 and 90 min) to obtain samples of particles with different cohesive forces. During this process, the inorganic functional groups of the Silane reacts with the different OH groups formed after cleaning with HCl and forms Si-OH groups. Finally, the treated glass particles were allowed to air-dry under a fume hood for 24 hours, forming a nanometric thin coating around the surface of the particle. The cohesive force of the formed coating mainly depends on the chemical silanization reaction duration and the Silane concentration.
cohesion force measurement. The cohesive force between the particles was measured using an in-house setup that we specifically designed for this study at the Laboratory of Physics and Physical Chemistry of Foods (Wageningen University, NL). The setup consists of a micro-balance (Sartorius) and a micro-positioner as shown in Fig. 7. The micro-positioner was driven by a DC electric motor with an adjustable speed varying from 0.05 to 1 mm/min. The micro-positioner was mounted above the micro-balance. The whole setup was installed on an www.nature.com/scientificreports www.nature.com/scientificreports/ optical table to reduce mechanical vibrations. One particle was fixed on the balance by a double-sided tape while the other particle was glued to a flexible thin rod connected to the micropositioner.
In order to stick the upper particle to the flexible rod without contaminating the whole particle surface with the glue, the tip of the rod was first wetted with about 0.05 micro-liters of instant liquid resin. After about 20 seconds, the liquid resin became more viscous, and then the upper side of the particle was precisely put in contact with the glue. After another 30 seconds, the particle was firmly connected to the rod. In order to measure the adhesive force, first the upper and lower particles were aligned as shown in Fig. 7. The weight in the balance is set to zero, subtracting the weight of one glass particle. Then, the top particle was moved downward with a speed of 1 mm/min. Upon the first mechanical contact between the particles, the balance shows a positive weight. Then, the stage is moved upward slowly (0.05 mm/min). When the two particles separate, the maximum absolute value obtained is recorded as the cohesive pull-off force between the particles. The stronger the particle-particle cohesive force, the stronger the snap-in force. This measurement has been repeated (more than 3 times) for several particles of different sizes and cohesive forces. The ideal accuracy of the micro-balance would be about 10 −8 mN. However, airflow around the setup and mechanical vibrations reduce the accuracy of our measurements to about 10 −5 mN.

Surface characterization: AfM measurement. The AFM observations were performed with a
Dimension Fast Scan atomic force microscope (Bruker's ScanAsyst and PeakForce Tapping AFM, Bruker Corporation, UK), located in the NanoLab NL, University of Twente. All AFM images were collected using the tapping mode with a silicon nitride probe and a scanning area of 2 × 2 μm. The glass particle samples were glued onto a double side silica tape on top of a glass substrate before the tapping measurements. the drum apparatus. The drum is made by a cylinder of R = 60.5 mm inner radius, and L = 22 mm length, held between two circular Plexiglass (PMMA) plates of 5 mm thickness to allow optical access. For a quasi-two-dimensional rotating drum, Jain et al. 3 argued that the length of the drum in the axial direction should be larger than 6.4 × r, with r the average radius of the particles, to neglect the front and back wall friction on the flow characteristics. The drum length used in this study is indeed larger than 6.4 × r of all the particles used. The drum was placed on a horizontal rotating axis driven by a variable-speed motor. Images of the rotating drum were recorded using a MotionBLITZ EoSens high speed camera working at 120 fps.
Experiments were performed in the cascading regime, at a rotation speed of the drum ω = 25 rpm, corresponding to a Froude number Fr = 0.21. This number represents the ratio of centrifugal to gravitational acceleration: 2 where ω is the rotation speed of the drum, R its inner radius and g the acceleration due to gravity. Experiments are conducted using a selected set of Borosilicate glass particles of density 2500 kg/m 3 . The drum is less than half filled with the same mass of particles (i.e., 125 g of particles). Characteristics of the drum and the glass particles are summarized in Table 2.
image post-processing and mixing index computation. The images, acquired using the high-speed camera, were post-processed to obtain the dynamic angle of repose and the mixing index in the rotating drum. The dynamic angle of repose was obtained using the particle tracking package Trackmate within the FIJI ImageJ distribution 44 . First, we removed the background in the image sequence and adjusted the light to enhance particle detection. Particles outlines (i.e., spots) that stand out from the background were segmented and identified based on the difference of Gaussians distributions 45 with an estimated particle diameter of 8 pixels for a particle radius of r = 0.85 mm. The detected spots were converted into a table of positions and visualized using ParaView 46 . The dynamic angle of repose was computed by linear regression of the positions of the particles in the free surface of the bed using an in-house python code. For the computation of the mixing index, particles with different colors were used in the experiments to identify cohesive, non-cohesive, big and small particles during image post-processing using variable tolerances in FIJI ImageJ software. Then, binary images with different particle species were obtained and converted into continuum density fields using a method similar to Coarse Graining method (CG) 47 . The main difference between the two methods is the application of the CG to the pixel data rather than the particle positions. This fully removes the need for a particle size and position detection algorithm, leading to a significant reduction in analysis time. To obtain the density field, we assumed the mass to be equal to 1, which gives a CG density function of:

Properties Value
Drum, R × L (mm) 60.5 × 22 Drum rotation speed ω (rpm) 25 Glass particles radii r (mm) 0.85, 1.25 and 2 particle density ρ (kg/m 3 ) 2500 Filling level of the drum β (%) 35 with N A and N B , respectively, the number of particles "A" and "B" in the bulk. Another way of calculating the number of contacts per particles was proposed by Chandratilleke et al. 52 by using the coordination number, who also showed that this number is useful to identify the mechanisms of mixing and demixing in particulate systems. The free binding energy, or work of adhesion, done to separate unit areas of two surfaces or media from contact to infinity in a vacuum can be related to the adhesion force F AB between two particles by 53 : AA AA A F AB and F AA denote, respectively, the adhesion and cohesion forces existing at a single interparticle contact, r A and r B are the radii of particles "A" and "B" respectively, and W AB refers to the energies per unit area of the "A"-"B" interface. In a bulk of a mixed particles, each reference particle "A" will share Z AB adhesive bonds and Z AA cohesive bonds. Thus, the average bonding work per particle of the perfectly mixed system (see Fig. 8 For a completely segregated system, the number of adhesive bonds is very low comparing to the number of cohesive bonds (see Fig. 8(e)), and therefore the former can be neglected. We can then write the average bonding work of the segregated system as follows:

s A A AA A BB BB B
Here, for a segregated system (Fig. 8(e)), Z AA = Z BB = Z = 6. The difference of work ΔW = W m − W s gives the average cohesive work change (i.e., energy change per contact area) per particle, when the system goes from a segregated to a mixed state. Notice that for a monodisperse system, if F AB = F AA , ΔW is equal to zero.