Abstract
Predicting the mixing or segregation of granular of dense materials is difficult due to their athermal nature and lack of general theories. The goal of this research is to improve the design (diameter of the cylinder) and operation (rotational frequency of the screw conveyor) of the stirrer to attain good mixing. We characterize the mixing by some macroscopic measures and relate its behavior with the design parameters. Afterward, we compute some microscopic measurements that characterize the movement of the spheres and relate them with the operating conditions. The quality of mixing is found to depend on the rotation frequency of the screw conveyor and the gap between the latter and the enclosing cylinder. Besides, some relationships between micro- and macroscopic measures are drawn. This approach can provide guidance on the design or operation of devices to handle granular media when their microscopic behavior is previously known.
1 Introduction
Granular materials are justly famous for being difficult to mix. Because of their intrinsic out-of-equilibrium nature, they display anti-intuitive behaviors like the increase in order due to segregation (). Since granular materials do not naturally diffuse (their effective temperature is zero), an input of energy is required to start movement. This is usually done by vibrating, tapping, or rotating a recipient full of grains. Segregation under agitation (the Brazil nut effect) is one of the most studied problems in dense granular matter (; ; ; ). Granular particles also segregate under rotation in a cylindrical tumbler (; ; ; ; ; ; ) or a rotating drum (). Other surprising behaviors include the horizontal segregation under vertical vibrations (; ), the reversal of segregation (), and even pattern formation (). In all cases, the differences in the particles’ size, mass, or density are the subject of a mechanism that leads to the separation of the particles according to their nature. The segregation of a massive intruder can be caused by inertia () under tapping or by Archimedean sinking under vibration (); under similar conditions, a large intruder will be segregated due to steric effects (; ). In large-enough containers, convection currents are the leading mechanism (). Similar processes are at play in tumblers along with new processes like avalanching and waves. The cause can also have a more subtle origin like a difference in the friction coefficient (; . For dilute to moderately dense media, the kinetic theory can be used to explain segregation and mixing using a granular temperature proportional to the fluctuations of the velocities, as shown in ).
In addition to its scientific interest, segregation has important practical consequences (). In many production processes, it is necessary to mix particles, e.g., to produce cement (), create drugs (; ), and handle foodstuffs (). Procedures for sorting require the opposite (). Thus, understanding the causes and mechanisms of segregation will lead to optimized industrial processes. Regarding mixing, in the 1950s, it was thought that practice was the best source of knowledge on techniques and apparatuses to mix solids (). However, soon the necessity to conduct the basic theoretical work was recognized (; ). It was pointed out that the deduction of rheological properties on simple geometries would be helpful in understanding the behavior of particulated solids in complex mixers. This idea rooted in the field, and in the next years, studies that analyzed the diffusion of particles in mixers emerged (; ; ). These were complemented with the definition of phenomenological quantities, like the mixing index. This is just an adimensional quantity defined to take values between 0 and 1, the former corresponding to a segregated state and the latter to a well-mixed state. Several definitions have been put forward (; ; ; ; ; ). With the advancement of computational resources, the first numerical studies grew to full-scale simulations. The literature is abundant. DEM studies on apparatuses similar to the apparatus studied here can be found in ), ), ), ), ), ), and ). Researchers have also resorted to simulations using different methods, e.g., (; ).
In the aforementioned studies, one wishes to mix particles with different properties. However, a pertinent question arises: what happens when all particles are identical? In this case, we can still speak of mixing because we can put a label in each particle and follow its trajectory. It is possible to define measures of mixing that will tell us whether particles that were close at the initial instant remain close after a certain amount of time. The mixing of identical grains under vertical vibration was experimentally studied in ), concluding that it is a self-diffusive process. Computational studies of beds of identical grains under shear (; ) were also found to be diffusive processes. In ), two regimes were identified (flow and voidage) depending on the density of the bed. The influence of flow parameters on self-diffusion in a vibrated granular bed was studied by simulation and experiment in ). The crossover from a ballistic to diffusional regime was investigated experimentally in ) using image processing technology to measure the self-diffusion coefficients of hard disks in a vertical vibrated bed. Simulating granular convection () measured a self-diffusion coefficient and characterized the underlying physics as a phenomenon of mass diffusion induced by energy injection. In a different system, ) analyzed shearing experimentally and modeled it by a random walk, confirming the diffusive nature of the flow. More recently, ; ) measured the self-diffusion coefficients in discrete numerical simulations of steady, homogeneous, collisional shearing flows of nearly identical, frictional, and inelastic spheres. They considered homogeneous and inhomogeneous flows and modeled their results by a random walk, confirming the diffusional nature of the process.
In this work, we analyze the mixing of identical granular particles with the aim of optimizing mixing. Our motivation is a practical application: mixing heated particles to homogenize an applied temperature. The system chosen is a screw conveyor enclosed in a cylindrical case, as shown in Figure 1. The process starts by filling the mixer with ≃2000 balls. After they reach equilibrium at the bottom of the mixer, the screw conveyor starts rotating at a fixed velocity, causing the particles to move upward. The design of the screw conveyor is such that the particles accumulate close to the pole and fall by the sides. This pattern of movement should produce a homogeneous mixing and facilitate the even distribution of the applied temperature.
FIGURE 1
We run several simulations in which we change the rotational frequency of the screw conveyor, the diameter of the enclosing cylinder, and the size of the particles. In each case, the performance of the apparatus is assessed and analyzed as a function under the operating conditions. Next, we characterize the microscopic variables, such as the mean-squared displacement, and try to find a connection between the global performance of the mixer and the microscopic description of the medium. The aim is to predict whether a given granular material will mix well based only on known diffusive behaviors or to choose (or restrict the range of choices) the best granular material for given mixing application before having to make costly and time-consuming experiments.
2 Materials and methods
2.1 Numerical method
We perform simulations using the well-known discrete element method (DEM) (). LIGGGHTS software () is used to implement them. This is free and open-source software based on LAMMPS (), which incorporates many models specific to the simulation of granular materials. It is parallelized and optimized allowing to run fast simulations. The diameters of the particles are d = 12 mm for the reference case and 15 and 9 mm for the additional cases, which we also refer to as big and small particles. The diameter of the screw conveyor inside the stirrer is Ds = 128.23 mm and is fixed for all simulations. The diameter of the enclosing cylinder is variable. We used D = 20, 18, 17, and 16 cm for d = 12 mm, D = 28, 24, 20, and 18 cm for d = 15 mm, and D = 22, 20, 18, and 16 cm for d = 9 mm. To rationalize all these measures, we define the adimensional gap λ = (D − Ds)/2d, which gives the number of particles that fit between the border of the screw conveyor and the wall of the cylinder. The length of the cylinder is Lc = 55 cm and is fixed for all simulations. The horizontal coordinates are X and Z, while Y is the vertical coordinate.
We consider only contact forces. Given the size of the particles and the material (SiC), we do not expect that attractive forces develop. Two particles are in contact whenever their distance is less than the sum of their radii rij < 0.5 (di + dj). The contact model consists of two forces, one normal to the surfaces in contact, i.e., along the line that joins the centers of the particles, and a tangential force perpendicular to the first, which are given byIn the first equation, the spring force is proportional to the overlap of the particles. The second term is a damping term proportional to the relative velocity of the colliding particles. The parameter kn is the elastic constant, and γn is a viscoelastic damping constant whose role is to dissipate energy during the collision. Analogously, the force in the tangential direction depends on a restoring term, proportional to the sliding of the particles, and a damping term that dissipates energy. The damping constants are related to the restitution coefficient (Section 2.1). They control the amount of energy lost in a collision and thus influence the rate at which the medium settles in a static state. The tangential overlap ξtij increases, while the contact lasts but is truncated to fulfill the Coulomb criterion Ft ≤ μFn, where μ is the friction coefficient. For the simulation, we fix the values of the Young’s modulus Y, Poisson’s ratio ν, restitution coefficient e, and friction coefficient μ. These are related to the constants of the force model as follows:
The mixed variables are defined as follows:
For this research, all particles in a given simulation are identical, with radii as mentioned previously. The other parameters have values Y = 5 × 108 Pa, ν = 0.2, e = 0.9, and μ = 0.5, and the mass is fixed by setting the density ρ = 3600 kg/m3. These are typical magnitudes given by the manufacturers (; ). The value of gravity is set to g = 9.8 m/s2, and the integration time step is δt = 10–6 s. We have checked that this set of values gives rise to a realistic granular behavior; e.g., in Figure 2A, the average overlap between spheres is well below 1%, and increasing the Young’s modulus does not make a difference.
FIGURE 2
The collision of a particle with a wall is treated in the same way as particle–particle collisions, but assuming that the second particle has infinite mass and radius (flat wall limit). The properties of the second particle are identical to the particles of the granular medium.
2.2 Simulation protocol
A simulation starts with an empty stirrer which is filled by a rain of particles released from the top. The number of particles is fixed for design reasons to 2,060. During this phase, the screw conveyor is not moving. When all particles have been released, they are allowed to settle until their kinetic energy drops to J, at which point the bed can be considered at rest. Then, the screw conveyor starts to move at a given rotational speed. We have used three different frequencies ω = 36, 180, and 360 rpm. The simulations that run until 10 s have been simulated. As shown in Figure 2B, this provides ample time to reach a stationary state in all cases.
During a simulation, we store frames of the particles’ positions, velocities, forces, and other parameters at intervals of 0.2 and, sometimes, 0.1 s. These are later used for analysis.
3 Results
3.1 Macroscopic measures
The first result we focus on is the mixing index. As mentioned in Introduction, this adimensional quantity takes values between 0 (segregated state) and 1 (mixed state). Among the several definitions available in the literature, for this work, we choose the index proposed in ). This index is defined as follows:
Here, rij is the distance between particles i and j at time t = 0 and rt,ij is the distance at time t. The parameter R is a dimensionless quantity which can be chosen between 0 and 1. According to the results in ), we set it to R = 0.87. The quantity g is a scaling factor of dimensions length−1. We set its value to g = 1/Lc. As shown in Figure 2B, we see that this definition with these parameters allows distinguishing between the different scenarios during the run by changing the frequency of the screw conveyor and the parameter λ.
The mixing index starts at zero, as defined by Eqs 12, 13, and increases relatively fast up to a stationary value that remains fixed after a few seconds. The stationary value is higher for higher frequencies and exhibits lower values of λ. Not only the maximum Γ but the increasing speed also depends on the frequency and, to a lesser degree, on λ, with analogous trends. In Figures 2C, D, we can see this more clearly. The maximum Γmax increases approximately linearly with the frequency and λ. The mixing time τmix is calculated by fitting the data to the following equation:
As a function of the frequency, the mixing time decreases as the former increases. The decrease is not linear, being much faster for lower frequencies than for larger frequencies. Regarding λ, the mixing time does not change monotonously, and we observe that the curves in Figure 2D cross at some points.
An visual inspection of videos of the mixing process reveals two facts that explain those results. For the lowest frequency, the screw conveyor is not able to lift the spheres to a significant height. In this case, the bed remains rather dense and close to the bottom. When the frequency is increased, the particles acquire enough energy to climb up the screw conveyor. At this point, λ starts to matter. For relatively high values, the spheres fall from the screw conveyor since there is a wide gap available and, again, do not attain very high altitudes. Instead when λ is small, it is more difficult for the spheres to fall before reaching the upper part, and consequently, Γ increases.
We introduce the variable height by defining mobilized spheres: a sphere which reaches a height greater than 1/2Lc at some instant during the mixing process. In Figure 2E, we plot the percentage of mobilized spheres as a function of the maximum mixing index. It is worth noting that the values of Γmax = 0.7 can be reached with a small percent of mobilized particles. This means that most particles remain in the lower part of the stirrer, so the mixing is not actually effective as the beads do not increase and are exposed to the source of heat (not within the scope of this study). Most of these cases correspond to the lowest frequency, but cases with the highest frequency are also observed.
We can clarify the previous figure by adding the variable λ to the plot. As shown in Figure 2F, the size of the symbols is proportional to the adimensional gap λ. We can now observe that the low values of mixing index and mobilization are mostly associated with the large values of λ. The steep increase in mobilization at approximately Γmax = 0.7 corresponds to scenarios with low values of the adimensional gap. As mentioned previously, this is due to the fact that spheres fall easily if there is space available. So, in order to have a good mixing, we need a stirrer with high frequency and small gap, slightly above 1.
3.2 Microscopic variables
After studying the macroscopic measures that describe the mixing behavior, we turn to microscopic quantities that characterize the mobility of the spheres. In Figure 3A, we plot the mean-squared displacement of the spheres as a function of time for each coordinate and for the total displacement. We first note that in all cases, the displacement is superdiffusive, i.e., is a power law with an exponent larger than 1. Comparing with the reference line given, we note that it is actually larger than 1.5. This is not surprising since we are giving energy to the spheres. Before one second has elapsed, the MSD turns almost flat, which corresponds to the approaching of the maximum mixing index.
FIGURE 3
The magnitude of the total MSD increases with frequency, as shown in Figure 3B. As more energy is injected in the system, the spheres become more mobile. Note that the power law for the total MSD seems larger than 1.5. In addition to the quantitative differences, the qualitative behavior of this measure is quite similar in the three cases.
Plotted as a function of λ, as shown in Figure 3C, the total MSD increases in magnitude as the former decreases. This is consistent with the discussion in section 3.1. Since particles tend to fall from the screw conveyor in systems with large λ, the total displacement does not reach large values. So in this case, we see that a microscopic measure of the behavior of the medium seems correlated with the macroscopic behavior.
Finally, we show the velocity correlation function in Figure 3D. The velocities in each component decorrelate very fast, reaching Cv = 0 in split seconds. Afterward, the horizontal components alternate weak correlations and anticorrelations until seconds in which the correlation is zero. The vertical component behaves differently, showing only one oscillation of anticorrelation and remaining flat after 2 s. The extremely fast decorrelation at the beginning is due to the fact that spheres start in repose and are given a non-zero velocity suddenly as the screw conveyor starts to move at its prescribed frequency. In addition, the velocity of each sphere depends on whether it is in contact or not with the screw conveyor and the orientation of the contact. So in practice, the acquired velocities are very different. The oscillations seen in the horizontal components are due to the interplay between collisions (anticorrelation) and the effect of the screw conveyor that tends to displace the particles in the same direction (correlation).
To obtain more insights into these results, we fit the first part of the total mean-squared displacement (MSDr) to the following expression:
The value of the exponent, as shown in Figure 4A, Figures 4C, E, follows a descending trend with λ. The trend is noisy in some cases, but one can see that it is roughly linear for intermediate and large values of the adimensional gap. Instead, for values close to λ = 1, the trend flattens as it approaches an exponent equal to 2. The descending trend is understandable in terms of the role of λ discussed in previous sections. When λ decreases, the trend flattens as the exponent tends to the limiting value 2, which would indicate a ballistic regime. Large values of λ correspond to particles falling early from the screw conveyor, and these trajectories are very short and have a low weight in the computation of the MSD. This is the reason why we do not see a large value of a for low λ. Only the trajectories of spheres falling from a high altitude have an important enough weight to increase the exponent to near-ballistic values.
FIGURE 4
The natural logarithm of the prefactor is shown in Figures 4B, D, F. The linear trend indicates that b decays exponentially with λ. As could have been anticipated, the prefactor increases with the frequency of the screw conveyor. For all large values of λ obtained for d = 9 mm, there seems to be a flattening of the trend. This can be compared with Figure 2C, where the second part of the curve becomes more horizontal as the value of λ increases.
3.3 Correlation between micro- and macroscopic variables
Finally, we would like to know whether there exists some simple correlation between the macroscopic and microscopic variables studied. There are many possible combinations. In Figure 5A, we plot the maximum mixing index as a function of the natural logarithm of the prefactor of the total MSD for all scenarios. This relationship is linear, with some dispersion, indicating that the mixing increases with the logarithm of the mean-squared displacement, since we have seen that the prefactor sets the magnitude of the MSD. A relationship between Γmax and b should exist that could have been anticipated intuitively, but the precise nature of the relationship is not obvious.
FIGURE 5
The results shown in previous sections and intuition suggest that there should be a relationship between the maximum mixing index and the kinetic energy imparted to the spheres by the screw conveyor. In Figure 5B, we observe that Γmax increases linearly with the natural logarithm of the total kinetic energy. Hence, the relationship shown in Figure 5C indicates that the total kinetic energy and the mean-squared displacement are proportional (the plot is clearer on a logarithmic scale). The kinetic energy is measured in the steady state when the mixing index attains its maximum value. In this state, the total kinetic energy of the spheres has a well-defined average with small fluctuations.
In Figure 5D, we plot the maximum mixing index as a function of the decorrelation time τdecorr. The latter is defined from the velocity autocorrelation function for the X component, as the time needed for this function to reach the zero value for the first time . In this case, there is no simple relationship. The decorrelation time for ω = 360 and 180 rpm is in the range of 0.1–0.25 s for all the values of λ. Furthermore, there is no systematic dependence on the latter, and one can see large- and small-sized symbols covering the range of τdecorr values. The small symbols (low λ values) accumulate in the upper part of the plot (large Γmax) due to better mixing, as was discussed in reference to Figure 2F. The data points for ω = 360 tend to have lower values of the decorrelation time, but there is some mixing with the ω = 180 points and the dependence with frequency seems very weak in any case.
For the lowest frequency, we find a wide range of values of τdecorr, ranging from 0.45 to 1 s. For this frequency, the mixing degree increases with the decorrelation time, following a saturating curve. As in the case of the large frequencies, there is no systematic dependence on λ since we can find small and large symbols in very different values of the decorrelation time. From Figure 5D, one can infer that τdecorr depends on the frequency of the screw conveyor, strongly for low frequencies and very weakly for the larger frequencies. However, more data for intermediate frequencies would be needed to clarify the relationship between both quantities.
Finally, let us mention that plots of τdecorrversus frequency or even λ (not shown) show the same vertical organization. So in this case, we can conclude that τdecorr is not a useful variable to distinguish the different cases and discriminate the best scenario for optimizing mixing.
4 Summary and discussion
We have presented simulations to characterize and optimize the design and operation of a stirrer whose purpose is to mix a bed of identical granular spheres. We have first studied some macroscopic measures, such as the mixing index. We conclude that the mixing is improved by increasing the frequency of the screw conveyor and decreasing the adimensional gap λ. The first gives more energy to the spheres to move upward the screw conveyor, while the second prevents the spheres to fall by the side of the screw conveyor before reaching a high altitude. We have shown that the mixing index offers limited information though. Effectively, by defining what we have called mobilized spheres (those who reach above half the height of the stirrer), we show that high values of the mixing are compatible with a small mobilization. Thus, it is necessary to consider several measures to characterize mixing and attain a good design.
Next, we have studied some microscopic measures that characterize the movement of particles. The mean-squared displacement is superdiffusive in all cases due to the constant supply of energy. It behaves as a power law during the time in which the mixing index increases and afterward flattens. The value of the exponent changes linearly with the frequency of the screw conveyor and linearly with λ. In this case, the exponent tends to 2 as λ reaches 1; this case represents the limiting case in which the spheres reach a high altitude in the stirrer. The prefactor of the power law is similarly correlated with the frequency of the screw conveyor and λ. In the latter case, we observe a flattening of the curves for large λ which corresponds to the flattening of the MSD at large times.
Finally, we have seen that it is possible to find correlations between the microscopic and macroscopic measures. In particular, we find that the maximum mixing index increases with the logarithm of the prefactor of the mean-squared displacement. This type of relationship can be used to design or operate apparatuses to handle media whose microscopic properties are known. Although some correlation might be expected based on physical reasoning, the precise form it takes can be far from obvious. In that case, simulations can be critical, especially when dealing with dense granular media for which a general theory is still lacking.
The open question of whether these relationships can be fruitfully used in industrial design remains unanswered. More research is required to assess the degree of generalization to different systems or to define frameworks that can lead from microscopic properties to design guidelines.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
RA: conceptualization, data curation, formal analysis, investigation, methodology, supervision, visualization, writing–original draft, and writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work received funding from the Horizon Europe project PLASTICE (Grant Agreement ID: 101058540).
Acknowledgments
The author thanks his colleagues Carlos González and Ignacio Julián for helpful discussions concerning this work.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AokiK. M.AkiyamaT. (1996). Spontaneous wave pattern formation in vibrated granular materials. Phys. Rev. Lett.77, 4166–4169. 10.1103/physrevlett.77.4166
2
AokiK. M.AkiyamaT. (1998). Control parameter in granular convection. Phys. Rev. E58, 4629–4637. 10.1103/physreve.58.4629
3
ArtoniR.LarcherM.JenkinsJ. T.RichardP. (2021). Self-diffusion scalings in dense granular flows. Soft Matter17, 2596–2602. 10.1039/d0sm01846e
4
ArtoniR.LarcherM.JenkinsJ. T.RichardP. (2022). Self-diffusion in inhomogeneous granular shearing flows. Phys. Rev. E106, L032901. 10.1103/physreve.106.l032901
5
BaiL.ZhengQ.YuA. (2017). FEM simulation of particle flow and convective mixing in a cylindrical bladed mixer. Powder Technol.313, 175–183. 10.1016/j.powtec.2017.03.018
6
BarkerG. C.MehtaA. (1993). Transient phenomena, self-diffusion, and orientational effects in vibrated powders. Phys. Rev. E47, 184–188. 10.1103/physreve.47.184
7
BasinskasG.SakaiM. (2016). Numerical study of the mixing efficiency of a batch mixer using the discrete element method. Powder Technol.301, 815–829. 10.1016/j.powtec.2016.07.017
8
BourneJ. R. (1964). The mixing of powders, pastes and non-Newtonian fiuid. Trans. Inst. Chem. Eng.181, CE202.
9
BreuA. P. J.EnsnerH.-M.KruelleC. A.RehbergI. (2003). Reversing the Brazil-nut effect: competition between percolation and condensation. Phys. Rev. Lett.90, 014302. 10.1103/physrevlett.90.014302
10
BridgewaterJ. (1976). Fundamental powder mixing mechanisms. Powder Technol.15, 225–236. 10.1016/0032-5910(76)80051-4
11
CahnD. S.FuerstenauD. W. (1969). A probabilistic model of the diffusional mixing of particulate solids. Powder Technol.2, 215–222. 10.1016/0032-5910(69)80015-x
12
CahnD. S.FuerstenauD. W.HealyT. W.HoggR.RoseH. E. (1966). Diffusional mechanism of solid–solid mixing. Nature209, 494–496. 10.1038/209494a0
13
CampbellC. S. (1997). Self-diffusion in granular shear flows. J. Fluid Mech.348, 85–101. 10.1017/s0022112097006496
14
ChandratillekeG. R.DongK. J.ShenY. S. (2018). DEM study of the effect of blade-support spokes on mixing performance in a ribbon mixer. Powder Technol.326, 123–136. 10.1016/j.powtec.2017.12.055
15
CiamarraM. P.De ViziaM. D.FierroA.TarziaM.ConiglioA.NicodemiM. (2006). Granular species segregation under vertical tapping: effects of size, density, friction, and shaking amplitude. Phys. Rev. Lett.96, 058001. 10.1103/physrevlett.96.058001
16
CôtéP.AbatzoglouN. (2006). Powder and other divided solids mixing. scale-up and parametric study of a ribbon blender used in pharmaceutical powders mixing. Pharm. Dev. Technol.11, 29–45. 10.1080/10837450500463778
17
CullenP. J. (2009). Food mixing: principles and applications. (UK: John Wiley & Sons.
18
DippelS.LudingS. (1995). Simulation on size segregation: geometrical effects in the absence of convection. J. Phys. I Fr.5, 1527–1537. 10.1051/jp1:1995215
19
DuranJ.RajchenbachJ.ClémentE. E. (1993). Arching effect model for particle size segregation. Phys. Rev. Lett.70, 2431–2434. 10.1103/physrevlett.70.2431
20
FiedorS. J.OttinoJ. M. (2003). Dynamics of axial segregation and coarsening of dry granular materials and slurries in circular and square tubes. Phys. Rev. Lett.91, 244301. 10.1103/physrevlett.91.244301
21
GajjarP.JohnsonC. G.CarrJ.ChrispeelsK.GrayJ. M. N. T.WithersP. J. (2021). Size segregation of irregular granular materials captured by time-resolved 3d imaging. Sci. Rep.11, 8352. 10.1038/s41598-021-87280-1
22
GaoW.LiuL.LiaoZ.ChenS.ZangM.TanY. (2019). Discrete element analysis of the particle mixing performance in a ribbon mixer with a double u - shaped vessel. Granul. Matter21, 12–16. 10.1007/s10035-018-0864-4
23
García-ChamorroM.Gómez-GonzálezR.GarzóV. (2022). Kinetic theory of polydisperse granular mixtures: influence of the partial temperatures on transport properties. Entropy24, 1–40. 10.3390/e24060826
24
GillemotK. A.SomfaiE.BörzsönyiT. (2017). Shear-driven segregation of dry granular materials with different friction coefficients. Soft Matter13 (2), 415–420. 10.1039/C6SM01946C
25
HalidanM.ChandratillekeG. R.DongK. J.YuA. B. (2018). Mixing performance of ribbon mixers: effects of operational parameters. Powder Technol.325, 92–106. 10.1016/j.powtec.2017.11.009
26
HerseyJ. A. (1975). Ordered mixing: a new concept in powder mixing practice. Powder Technol.11, 41–44. 10.1016/0032-5910(75)80021-0
27
HillK. M.KakaliosJ. (1995). Reversible axial segregation of rotating granular media. Phys. Rev. E52, 4393–4400. 10.1103/physreve.52.4393
28
HuertaD. A.Ruiz-SuárezJ. C. (2004). Vibration-induced granular segregation: a phenomenon driven by three mechanisms. Phys. Rev. Lett.92, 114301. 10.1103/physrevlett.92.114301
29
JakličM.KočevarS.SrčičK.DreuR. (2015). Particle size-based segregation of pharmaceutical powders in a vertical chute with a closed bottom: an experimental evaluation. Powder Technol.278, 171–180. 10.1016/j.powtec.2015.03.021
30
JinX.ChandratillekeG. R.WangS.ShenY. (2022). DEM investigation of mixing indices in a ribbon mixer. Particuology60, 37–47. 10.1016/j.partic.2021.03.005
31
JingL.OttinoJ. M.LueptowR. M.UmbanhowarP. B. (2020). Rising and sinking intruders in dense granular flows. Phys. Rev. Res.2, 022069. (R). 10.1103/physrevresearch.2.022069
32
JuliánI.HerguidoJ.MenéndezM. (2016). Experimental and simulated solids mixing and bubbling behavior in a scaled two-section two-zone fluidized bedreactor. Chem. Eng. Sci.143, 240–255. 10.1016/j.ces.2015.12.025
33
KlossC.GonivaC.HagerA.AmbergerS.PirkerS. (2012). Models, algorithms and validation for opensource DEM and CFD-DEM. Prog. Comput. Fluid Dyn. Int. J.12, 140–152. 10.1504/pcfd.2012.047457
34
KnightJ. B.JaegerH. M.NagelS. R. (1993). Vibration-induced size separation in granular media: the convection connection. Phys. Rev. Lett.70, 3728–3731. 10.1103/physrevlett.70.3728
35
Kyocera (2023). Silcon carbide technical data. Available at: https://global.kyocera.com/prdct/fc/list/material/silicon_carbide/silicon_carbide.html.
36
LaceyP. M. C. (1954). Developments in the theory of particle mixing. J. Appl. Chem.4, 257–268. 10.1002/jctb.5010040504
37
LevanonM.RapaportD. C. (2001). Stratified horizontal flow in vertically vibrated granular layers. Phys. Rev. E64, 011304. 10.1103/physreve.64.011304
38
LiffmanK.MuniandyK.RhodesM.GutteridgeD.MetcalfeG. (2001). A segregation mechanism in a vertically shaken bed. Granul. Matter3, 205–214. 10.1007/s100350100093
39
MeierS. W.LueptowR. M.OttinoJ. M. (2007). A dynamical systems approach to mixing and segregation of granular materials in tumblers. Adv. Phys.56, 757–827. 10.1080/00018730701611677
40
Nahmad-MolinariY.Canul-ChayG.Ruiz-SuárezJ. (2003). Inertia in the Brazil nut problem. Phys. Rev. E68, 041301. 10.1103/physreve.68.041301
41
NajiL.StannariusR. (2009). Axial and radial segregation of granular mixtures in a rotating spherical container. Phys. Rev. E79, 031307. 10.1103/physreve.79.031307
42
NguyenD.RasmusonA.BjörnI. N.ThalbergK. (2014). CFD simulation of transient particle mixing in a high shear mixer. Powder Technol.258, 324–330. 10.1016/j.powtec.2014.03.041
43
OyamaY. (1939). Report no.5 (in Japanese). Bull. Inst. Phys. Chem. Res.18, 600.
44
OyamaY. (1940). Studies on mixing of solids. mixing of binary system of two sizes by ball mill motion. Sci. Pap. Inst. Phys. Chem. Res.37, 17.
45
Pinzón-MorenoD. D.RibeiroS.SaronC. (2020). Aggregates segregation in rubberized mortar. Structures26, 778–784. 10.1016/j.istruc.2020.05.009
46
PouxM.FayolleP.BertrandJ.BridouxD.BousquetJ. (1991). Powder mixing: some practical rules applied to agitated systems. Powder Technol.68, 213–234. 10.1016/0032-5910(91)80047-m
47
PrigozhinL.KalmanH. (1998). Radial mixing and segregation of a binary mixture in a rotating drum: model and experiment. Phys. Rev. E57, 2073–2080. 10.1103/physreve.57.2073
48
Qualitymaterial (2023). Porperties of silicon carbide. Available at: https://www.qualitymaterial.net/news_list85.html.
49
RapaportD. C. (2001). Mechanism for granular segregation. Phys. Rev. E64, 061304. 10.1103/physreve.64.061304
50
RapaportD. C. (2002). Simulational studies of axial granular segregation in a rotating cylinder. Phys. Rev. E65, 061306. 10.1103/physreve.65.061306
51
RapaportD. C. (2004). The art of molecular dynamics simulation. UK: Cambridge University Press.
52
RietzF.StannariusR. (2012). Oscillations, cessations, and circulation reversals of granular convection in a densely filled rotating container. Phys. Rev. Lett.108, 118001. 10.1103/physrevlett.108.118001
53
RosatoA.StrandburgK.PrinzF.SwendsenR. (1987). Why the Brazil nuts are on top. Phys. Rev. Lett.58, 1038–1041. 10.1103/PhysRevLett.58.1038
54
SavageS. B.DaiR. (1993). Studies of granular shear flows Wall slip velocities, ‘layering’ and self-diffusion. Mech. Mater.16, 225–238. 10.1016/0167-6636(93)90047-u
55
SchofieldC. (1976). The definition and assessment of mixture quality in mixtures of particulate solids. Powder Technol.15, 169–180. 10.1016/0032-5910(76)80045-9
56
SiiriäS.YliruusiJ. (2009). Determining a value for mixing: mixing degree. Powder Technol.196, 309–317. 10.1016/j.powtec.2009.08.009
57
SimonsT.Combarros-GarcíaM.GuptaP.TüzünU.ZiganS.SunJ.et al (2013). “Segregation and mixing of granular material in industrial processes,” in III International Conference on Particle-Based Methods — Fundamentals and Applications (PARTICLES 2013), Stuttgart, 18-20 September, 2013, 1.
58
SimonsT. A. H.BensmannS.ZiganS.FeiseH. J.ZetzenerH.KwadeA. (2016). Characterization of granular mixing in a helical ribbon blade blender. Powder Technol.293, 15–25. 10.1016/j.powtec.2015.11.041
59
StewartR. L.BridgwaterJ.ZhouY. C.YuA. B. (2001). Simulated and measured flow of granules in a bladed mixer—a detailed comparison. Chem. Eng. Sci.56, 5457–5471. 10.1016/s0009-2509(01)00190-7
60
ThompsonA. P.AktulgaH. M.BergerR.BolintineanuD. S.BrownW. M.CrozierP. S.et al (2022). Lammps - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comput. Phys. Commun.271, 108171. 10.1016/j.cpc.2021.108171
61
TsugenoY.SakaiS.YamazakiM.NishinomiyaT. (2021). DEM simulation for optimal design of powder mixing in a ribbon mixer. Adv. Powder Technol.32, 1735–1749. 10.1016/j.apt.2021.03.026
62
UtterB.BehringerR. P. (2004). Self-diffusion in dense granular shear flows. Phys. Rev. E69, 031308. 10.1103/physreve.69.031308
63
WeidenbaumS. S. (1958). Mixing of solids. Adv. Chem. Eng.2, 209. 10.1016/S0065-2377(08)60229-X
64
WildmanR. D.HuntleyJ. M.HansenJ. P. (1999). Self-diffusion of grains in a two-dimensional vibrofluidized bed. Phys. Rev. E60, 7066–7075. 10.1103/physreve.60.7066
65
YangS.-C.HsiauS.-S. (2001). Self-diffusion analysis in a vibrated granular bed. Adv. Powder Technol.12, 61–77. 10.1163/156855201744967
66
ZikO.StavansJ. (1991). Self-diffusion in granular flows. Europhys. Lett.6, 255–258. 10.1209/0295-5075/16/3/006
67
ZuriguelI.GrayJ. M. N. T.PeixinhoJ.MullinT. (2006). Pattern selection by a granular wave in a rotating drum. Phys. Rev. Lett.73, 061302. 10.1103/physreve.73.061302
Summary
Keywords
mixing, identical particles, simulation, microscopic variables, macroscopic measures
Citation
Arévalo R (2023) Connecting micro- and macroscopic variables in the design of a granular stirrer using simulations. Front. Soft Matter 3:1302367. doi: 10.3389/frsfm.2023.1302367
Received
26 September 2023
Accepted
09 November 2023
Published
27 November 2023
Volume
3 - 2023
Edited by
Lou Kondic, New Jersey Institute of Technology, United States
Reviewed by
Nagi Khalil, Rey Juan Carlos University, Spain
Ralf Stannarius, Otto von Guericke University Magdeburg, Germany
Updates
Copyright
© 2023 Arévalo.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Roberto Arévalo, rarevalo@fcirce.es
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.