Slow axial drift in three-dimensional granular tumbler flow

Slow axial drift in three-dimensional granular tumbler flow

Zafir Zaman Department of Chemical & Biological Engineering, Northwestern University, Evanston, Illinois 60208, USA    Umberto D’Ortona Laboratoire M2P2, UMR 7340 CNRS / Aix-Marseille Université, 13451 Marseille Cedex 20, France    Paul B. Umbanhowar Department of Mechanical Engineering, Northwestern University, Evanston, Illinois 60208, USA    Julio M. Ottino Department of Chemical & Biological Engineering, Northwestern University, Evanston, Illinois 60208, USA Department of Mechanical Engineering, Northwestern University, Evanston, Illinois 60208, USA The Northwestern Institute on Complex Systems (NICO), Northwestern University, Evanston, Illinois 60208, USA    Richard M. Lueptow Department of Mechanical Engineering, Northwestern University, Evanston, Illinois 60208, USA The Northwestern Institute on Complex Systems (NICO), Northwestern University, Evanston, Illinois 60208, USA
July 29, 2019

Models of monodisperse particle flow in partially filled three-dimensional tumblers often assume that flow along the axis of rotation is negligible. We test this assumption, for spherical and double cone tumblers, using experiments and discrete element method simulations. Cross sections through the particle bed of a spherical tumbler show that, after a few rotations, a colored band of particles initially perpendicular to the axis of rotation deforms: particles near the surface drift toward the pole, while particles deeper in the flowing layer drift toward the equator. Tracking of mm-sized surface particles in tumblers with diameters of 8–14 cm shows particle axial displacements of one to two particle diameters, corresponding to axial drift that is 1–3% of the tumbler diameter, per pass through the flowing layer. The surface axial drift in both double cone and spherical tumblers is zero at the equator, increases moving away from the equator, and then decreases near the poles. Comparing results for the two tumbler geometries shows that wall slope causes axial drift, while drift speed increases with equatorial diameter. The dependence of axial drift on axial position for each tumbler geometry is similar when both are normalized by their respective maximum values


I Introduction

Flow of granular media can be difficult to predict and challenging to model because of the inherent complexity of the collective motion of large numbers of particles. Yet, granular flows can sometimes be modeled using a continuum approximation, for example, as in the case of flow in a rotating tumbler Meier et al. (2007) and flow down an inclined plane or heap Midi (2004). The usual assumption in these geometries is that spanwise particle motion is primarily diffusive and averages to zero. However, there are situations where spanwise flow occurs. For example, in a partially-filled, long cylindrical tumbler rotating about its axis, endwall friction causes pathlines to curve near the endwalls Maneval et al. (2005); Chen et al. (2008). In the case of bidisperse particles, this out-of-plane flow can drive axial segregation and initiate axial band formation Zik et al. (1994); Hill and Kakalios (1995); Hill et al. (1997); Choo et al. (1997); Aranson and Tsimring (1999); Chen et al. (2010, 2011).

Figure 1: (a) Sketch of spherical tumbler flow modeled as a composite of two-dimensional slices. (b) Sketch of a two-dimensional slice model of circular tumbler flow in the continuous flow regime with a thin flowing layer near the surface and solid body rotation for material in the fixed bed. and denote the rotation speed and flowing layer velocity, respectively, and is the dynamic angle of repose of the free surface with respect to horizontal.
Figure 2: (Color online) Axial drift of particles demonstrated by band deformation in plane cross sections in a spherical tumbler ( cm) with fill fraction from experiment with  mm, (a) zero rotations and (b) three rotations (eight fixed bed passes), and from (c) DEM simulation with  mm after two rotations (three fixed bed passes). The colored band was initially located between  cm and  cm (measured axially from the tumbler equator) and spanned the tumbler in the direction. Drift occurs throughout the tumbler in the direction. White dotted lines on the 1-cm cross section highlight the similarity between band deformation in (b) experiment and (c) simulation, though drift per pass is greater in simulation.

Here we consider flow in partially-filled three-dimensional (3D) tumblers, such as spherical tumblers [Fig. 1(a)] Meier et al. (2007); Sturman et al. (2008); Juarez et al. (2010, 2012), rotating with angular velocity about a horizontal axis that intersects the tumbler at its “poles.” Our interest derives from the desire to extend the understanding of granular flow in quasi-two-dimensional (2D) cases Cantelaube and Bideau (1995); Clément et al. (1995); Hill et al. (1999); Khakhar et al. (2001); Hill et al. (2003); Jain et al. (2004); Midi (2004) to fully 3D flows. We consider the situation where the free surface is essentially flat and continuously flowing Meier et al. (2007). In this regime, the surface of the flowing layer maintains a dynamic angle of repose with respect to horizontal which depends on the frictional properties and diameter of the particles, and the rotational speed of the tumbler  Taberlet et al. (2003); du Pont et al. (2003); Doppler et al. (2007); Pignatel et al. (2012). A simple model of 3D flow in a spherical tumbler Meier et al. (2007); Juarez et al. (2010, 2012) assumes that flow in each plane perpendicular to the axis of rotation is essentially that of a two-dimensional circular slice of the appropriate diameter, as shown in Fig. 1(a). In this reduced quasi-2D geometry, particles enter the upstream end of the thin, rapidly flowing layer from the fixed bed (the region of particles in solid body rotation), flow downslope, and return to the fixed bed following the idealized streamlines shown in Fig. 1(b). The streamwise velocity profile in the flowing layer decreases approximately linearly with depth Khakhar et al. (1997); Jain et al. (2002); Ottino and Khakhar (2002), while the portion of the ‘solid body’ region nearest the flowing layer exhibits a much slower creeping motion Komatsu et al. (2001); Socie et al. (2005); Arndt et al. (2006).

Despite the attractiveness of its simplicity, there are indications that the 2D flow assumption is imperfect. Simulations of monodisperse flow in a partially-filled spherical tumbler Chen et al. (2009) indicate a slight out-of--plane curvature in the trajectories of surface particles. Asymmetries between the upstream and downstream portions of the curved trajectories manifest as axial drift. Unlike a partially-filled cylindrical tumbler Santomaso et al. (2004); Pohlman et al. (2009), the trajectory curvature cannot be directly attributed to frictional endwall effects. Rather, the curved particle trajectories appear to be related to the relative curvature of the tumbler walls to the surface of spherical particles as the paths of large particles curve more than the paths of small particles for the same tumbler diameter.

In this paper, axial drift is examined in experiment and simulation. Experimentally, particles are tracked on the surface and along the tumbler wall to measure the axial drift. Cross sections of the tumbler are also imaged to visualize the axial motion of colored particles within the bed. Discrete element method (DEM) simulations are performed to obtain particle trajectories and velocities to understand the axial drift in more detail. Similar axial drift occurs in both experiment and simulation over a range of tumbler and particle parameters, confirming the robust nature of the phenomenon.

Figure 3: Calculated axial variation of normalized solid body rotation time, , vs. normalized axial position, , for in a spherical tumbler. Curves extend from the equator to the maximum axial position of the free surface. The region between vertical dashed lines represents the initial position of the tracer band () in experiments and simulations.
Figure 4: (Color online) Band deformation in spherical tumblers is robust to variation in fill fraction as shown in -plane cross sections containing the axis of rotation from simulation. Fill fraction, rotation number and approximate number of passes through the flowing layer (in parenthesis) are indicated above each image.  mm,  cm, and  rpm.

Ii Band Deformation

The tumblers in these experiments were clear acrylic spheres rotated at constant angular velocity about a horizontal axis by a motor. In the band drift experiments, a -cm diameter tumbler was filled to a fill fraction (by volume) varying between and with  mm diameter soda-lime glass beads (SiLiglit Deco Beads, Sigmund Lindner GmbH, Germany) and rotated at rpm. A tracking band [Fig. 2(a)] was formed by filling the space between two thin partitions (at cm and cm) with light colored beads and filling the exterior with dark colored beads. To reduce electro-static interactions between particles and the wall, the inside of the tumbler was either wiped with an antistatic wipe (Staticide, ACL Inc., Chicago, IL) or treated with an antistatic spray (SP 610, Sprayon Inc., Cleveland, OH) prior to filling. To reduce inter-particle static charging when the relative humidity was below 50%, a small amount of deionized water (typically 2-5 L) was allowed to evaporate in the sealed tumbler to increase the relative humidity above 50%.

To view cross sections on planes perpendicular to the flow direction in experiment after tumbling, the grain bed was immobilized by pouring hot gelatin into the tumbler. Once the gelatin cooled to room temperature, the tumbler was placed in a freezer to set the gelatin. The tumbler was then cut in half along the meridional plane and photographed. Additional parallel planes behind the initial cut were exposed by carefully removing particles with a scraper. Two sets of images were obtained for each experiment from the halves created by the initial cut.

Before presenting the experimental results, we note that the time between surface particle passes through the flowing layer is always less than the tumbler rotation period for and, as a consequence, the number of particle passes (cycles) is always greater than the number of tumbler rotations. This occurs because the bed surface in any quasi-2D slice perpendicular to the rotation axis subtends an angle, measured with respect to the rotation axis (Fig. 3 inset). Ignoring the relatively much shorter time spent in the flowing layer (see Appendix A), . For all and thus varies with axial position since the surface of the particle bed does not extend to the poles. Figure 3 shows the predicted value of as a function of normalized axial position where is the tumbler radius. For fixed increases with while for fixed decreases (increases) with for (). At is independent of

Tumbler diameter (cm) 10, 14
Particle diameter (mm) 2
Particle density (kg m) 1308 Drake and Shreve (1986); Foerster et al. (1994); J. Schäfer et al. (1996)
Friction coefficient111For particle-particle and particle-wall interactions. 0.7
Restitution coefficient111For particle-particle and particle-wall interactions. 0.87
Fill fraction 20%-70%
Collision time (s)  Ristow (2000); Campbell (2002); Silbert et al. (2007)
Integration time step (s)  J. Schäfer et al. (1996)
Stiffness coefficient (N m)  J. Schäfer et al. (1996)
Damping coefficient (kg s) 0.206
Table 1: DEM simulation parameters
Figure 5: (Color online) Top view of a colored band in a spherical tumbler for increasing rotation from experiment (initial band position denoted by vertical lines). Rotation number and approximate number of passes through the flowing layer (in parenthesis) are indicated above each image. ,  mm,  cm, and  rpm.

Cross sections of the bed in Fig. 2(b) after three rotations for reveal axial drift of particles, although collisional diffusion has somewhat blurred the band boundaries. The deformation of the band at and beneath the free surface indicates the direction and degree of axial drift. The band shifts toward the pole at the top surface and at the bottom wall of the tumbler, while it shifts toward the equator in the mid-layer. White dashed lines for the cross-section 1 cm behind the meridional cut highlight the band deformation. Particles in the band pass through the flowing layer about eight times for after three tumbler rotations. Thus, axial drift is relatively small compared to the streamwise particle displacements. However, it is large enough to dominate over axial mixing due to collisional diffusion.

The colored band experiment is replicated using DEM simulations (see Appendix B for details and Table 1 for parameters). Figure 2(c) demonstrates that tracer band deformation in the simulation is similar to that in the experiment: particles drift toward the poles at the free surface and toward the equator lower in the flowing layer. (see Supplemental Material 111 See Supplemental Material for a video of the band deformation over time.) Particles near the surface and tumbler wall in the simulation drift further toward the pole than in the experiment, possibly due to differences in particle properties (i.e.,  Chen et al. (2009), , and ) between the experiment and the simulation. Additional experiments at and and simulations at (see Fig. 4) exhibit similar axial drift. Qualitatively similar drift for fill fractions equal to, less than, and greater than rules out potential mechanisms for drift related to axial variation of cycle time, as shown in Fig. 3.

Further insight into axial drift was obtained from images of the free surface distribution of the light colored particles initially in the tracking band (Fig. 5). After four rotations, the band is closer to the nearest pole and becomes wider and less distinct. After 25 rotations, beads from the band are closer to the equator. Collisional mixing of light and dark colored particles is substantial, but light colored particles remain almost exclusively in the hemisphere of the tumbler in which they started (left). To quantify this surface motion, the streamwise-averaged image intensity of the light colored beads was fit to a Gaussian. As determined by the peak of the Gaussian, centers of the light colored band are  cm for 0, 4, 10 and 25 rotations (respectively). Assuming that mean surface velocities are constant, the non-monotonic motion of the band center suggests that, in addition to axial flow, particles also move vertically through the layer, descending at the poles and ascending near the equator (see Sec. III.1).

Figure 6: Mean particle trajectories from simulation at two different initial depths viewed in the -plane exhibit net axial displacement per pass through the flowing layer. Bold (light) trajectories start (starting location denoted by asterisks) 0.2 cm (1.8 cm) from the tumbler wall in the fixed bed and traverse the flowing layer near the top (bottom) of the layer.  mm,  cm, , and  rpm.
Figure 7: (Color online) A representative mean particle trajectory for 25 tumbler rotations (45 flowing layer passes) from a spherical tumbler simulation with  cm viewed from (a) the side (-plane) and (b) the top (parallel to ); tumbler rotations (passes) are noted above each column. The trajectory is shown in its entirety (first column) and in 6.25 tumbler rotation segments. The free surface is indicated by the black curve in (a), and the start of each trajectory section is indicated by an asterisk (*). Straight portions of the trajectory in (b) correspond to solid body rotation, while curved portions correspond to flow. The full trajectory starts 6.8 mm from the tumbler wall at . It advances toward the pole through the top of the flowing layer and returns back to the equator through the bottom of the flowing layer. ,  mm, and  rpm.

Iii Drift Kinematics

iii.1 Mean trajectories

To better understand the axial drift, we examine average particle trajectories constructed by integrating the mean velocity field from simulation using a second order Runge-Kutta method. Figure 6 shows the effects of initial particle depth on pairs of trajectories at three axial locations viewed in the plane (looking down onto the free surface). The two trajectories in each pair start (indicated by a filled circle) in the fixed bed and the plane 0.2 cm and 1.8 cm away from the tumbler wall and enter the flowing layer near the free surface (bold curve) and near the bottom of the flowing layer (light curve), respectively. Near surface trajectories enter the flowing layer further upstream and exit it further downstream compared to deeper trajectories. More importantly, near surface trajectories move further toward the pole in the upstream half of the flowing layer than away from it in the downstream half of the flowing layer resulting in a net poleward motion as previously reported for segregating particles Chen et al. (2009). In contrast, deeper trajectories exhibit the opposite behavior and move closer to the equator each pass. Axial displacement is greatest at the upstream and downstream ends of the flowing layer, while the trajectory in the middle portion of the flowing layer is nearly straight.

Figure 8: Non-dimensional axial surface drift per pass vs. normalized axial position in a spherical tumbler with -mm particles. (a) from experiment (black points) and simulation (gray points) are similar. Mean axial drift with 95% confidence interval error bars for experiment (white circles) and for simulation (white squares) match within uncertainty (error bars for simulation are smaller than the symbols). Axial distance to the poles (dashed lines) bounds the axial displacement. and  rpm. (b) Mean axial surface drift in spherical tumbler simulations increases with . is the half width of the free surface.

To understand the global motion resulting from the axial drift per pass , Fig. 7 presents a nearly complete mean particle trajectory for 25 tumbler rotations and the same trajectory in shorter (6.25 rotation) increments. manifests as the spacing between the trajectory loops. In the top row, the view along the axis of rotation ( plane) shows the dynamic angle of repose of the particles as well as the depth of the particle trajectory. In the bottom row, the view is down the gravity vector with flow from top to bottom, so the trajectory is asymmetric with respect to the horizontal axis of rotation due to the repose angle of the tumbling particles.

From 0 to 6.25 tumbler rotations, the trajectory drifts toward the pole with each cycle through the flowing layer and fixed bed, and the axial displacement with each loop increases slightly approaching the pole. In the top view, trajectories are straight during solid body rotation in the fixed bed and curved in the flowing layer where the axial position changes. From 6.25 to 12.50 tumbler rotations, the trajectory is deep in the flowing layer and migrates back toward the equator, as is evident from the side view. From 12.50 to 18.75 tumbler rotations, the trajectory continues to move toward the equator, but with decreasing axial displacement, as it simultaneously approaches the free surface. From 18.75 to 25 tumbler rotations, the trajectory drifts slowly back toward the pole through the top of the flowing layer. For the entire trajectory, the axial displacement per pass is small near the equator and pole and larger in between. The axial drift behavior in Fig. 7 is consistent with that in Fig. 6, where poleward drift occurs when the particle is near the surface of the flowing layer, and drift toward the equator occurs when the particle is deep in the flowing layer. The trajectory circulates through the domain of the tumbler. The poleward axial flux of particles near the surface of the flowing layer is balanced by the axial flux of particles toward the equator near the bottom of the flowing layer.

iii.2 Surface Drift

The trajectory results reported in the previous section are derived from average velocity fields. Here we consider the axial drift per pass, of individual particles at the surface in experiments and simulations to better understand the stochastic nature of the drift process. is defined as the change in axial position of a surface particle in the fixed bed after a single pass through the flowing layer. In simulation, particles in the fixed bed whose centers remain within 1.5 mm (0.75) of the tumbler wall after passing through the flowing layer are considered surface particles. Figure 8(a) shows that measured from simulation (gray points) is stochastic with a nearly constant width of about 10 except near the poles where it is reduced due to the finite size of the tumbler which restricts the magnitude of poleward displacement (dashed lines). Axial displacements near the surface display a bias that depends on axial position. Poleward displacements outnumber equator directed displacements except near the pole () where particles are more likely to drift toward the equator. To quantify the mean drift, individual displacements are grouped into bins of width and the mean of each bin are computed (white squares in Fig. 8(a)). Error bars representing the 95% confidence interval are smaller than the symbols due to the large number of samples in each bin. The individual displacements are well fit by a cubic polynomial 222The weighted cubic polynomial was accepted as an appropriate fit for the mean axial drift using the F-test Kutner et al. (2005) for lack of fit at the 1% significance level. The individual displacements were weighted using sampling weights Lohr (2009) , where is the total number of data points and is the number of individual displacements in the bin to which the data point belongs. The sampling weight allows every axial position to be considered equally regardless of axial sampling variations.. The mean axial drift per cycle is zero at the equator (presumably due to symmetry of the tumbler about this plane), increases to a maximum of about at , and then decreases to reflecting the eventual reversal of the drift for particles near the pole (due to the geometric limits on poleward displacement as shown by the dashed lines in Fig. 8(a)).

Quantifying axial drift in experiment is more challenging because all particles cannot be easily tracked. Instead, a single tracer particle identical except for color was tracked at the clear wall of the tumbler in the fixed bed. A camera, mounted below the spherical tumbler, imaged fixed bed particles in contact with the wall at 30 fps. Since a single tracer appears at the wall infrequently, tumbler rotations were required to obtain a minimum of 20 sets of observations (pairs of tracer sightings within 6 passes) per bin of length . In comparison,  5000 observations per bin were obtained by tracking all particles at the wall in simulation. Moreover, in experiments, tracer displacements could only be obtained for due to the camera’s inability to image the steeply curved tumbler wall near the pole. However, within this restricted region, the distribution of individual axial displacements from experiment (black points in Fig. 8(a)) is similar to that from simulation (gray points) as is the mean axial drift (white circles calculated as in the simulation) whose error bounds overlap the mean axial drift from simulation (white squares). Since results from experiments and simulations are similar, simulations are again used to further probe axial drift at the surface.

The effects of varying tumbler diameter, speed, and fill fraction on the mean axial surface drift from simulations are shown in Fig. 8(b) for the right hemisphere of the tumbler. All simulations demonstrate the same qualitative behavior: drift is 0 at the equator and increases approximately linearly to a maximum around (except for , which occurs slightly further from the equator). For , maximum axial drift of the  cm simulation is 1.8 times that in the  cm simulation. Varying from 10 rpm to 30 rpm for both values of has no impact on the axial drift per pass. For  cm, the maximum axial drift is highest for , just slightly greater than that for but significantly greater than despite having the same free surface geometry as . The reduced drift at may be influenced by differences such as the greater bed depth near the pole and the existence of an unmixed core, which never enters the flowing layer Jain et al. (2002); Bonamy et al. (2002); Socie et al. (2005). In addition, the axial drift depends on particle properties. Increasing the restitution coefficient in simulations from 0.5 to 0.8 decreases the maximum axial drift by 30%, while increasing the friction coefficient from 0.5 to 0.9 decreases the maximum axial drift by 16%.

Figure 9: Top views of (a) spherical tumbler of diameter, , and (b) truncated double cone tumbler of width, equator diameter, , and pole diameter, . Rotation is about the -axis, which is also the axis of symmetry for the double cone tumbler. (Left column) Bed cross sections at fill fraction . (Right column) Sketches of the flowing layer surface (dashed curves) in the -plane for with flowing layer width, , along the axis of rotation. The flowing layer length, , varies axially for both spherical and conical tumblers. is the axial width of the tumbler while is the axial distance that the free surface extends in the tumbler, which may be less than the width of the tumbler depending of the fill level.
Figure 10: (Color Online) Mean particle trajectory for a double cone tumbler simulation: (a) side view ( plane) and (b) top view (parallel to ). Trajectory is shown for 35 rotations (63 flowing layer passes) in its entirety (first column) and in 8.75 tumbler rotation segments with  cm,  cm,  cm, and . Tumbler rotations are given above each column with flowing layer passes in parentheses. Similar to the spherical tumbler (Fig. 7), straight portions of the trajectory in (b) correspond to solid body rotation while curved portions correspond to flow Note4 (). The free surface is indicated by the black curves in (a). The trajectory starts 0.6 mm below the surface and directly below the axis of rotation at (asterisk), advances towards the pole through the top of the flowing layer, and returns to the equator through the bottom of the flowing layer. ,  mm, and  rpm.

Iv Tumbler geometry

While the experimental and simulation results for spherical tumbler flow reveal an axial drift that depends on axial position, they leave unresolved the mechanism of axial drift. Previous results Santomaso et al. (2004); Chen et al. (2008) indicate that axial drift is negligible in cylindrical tumblers (except near the endwalls where wall friction drives the flow), suggesting that axial drift depends on axial variation of the flowing layer length, or wall angle, measured relative to the rotation axis in the plane of the free surface. To explore the effects of wall geometry on axial drift, simulations were done in double cone tumblers because in this geometry varies with axial position while remains constant (Fig. 9).

We first note that the perimeter of the free surface in the -plane of spherical and double cone tumblers (Fig. 9) depends on . For the free surface perimeter (solid curves) is characterized by the tumbler diameter ( for spheres and and for double cones) and by the tumbler width ( for spheres and for double cones). However, for , the free surface perimeter (dashed curves in the right column of Fig. 9) does not necessarily span the width of the tumbler, resulting in an effective free surface width Regardless of , the free surface in a spherical tumbler retains the same shape (circular) unlike the double cone tumbler whose free surface is bounded by two hyperbolas that depend on  333Flow in double cones rotated about the -axis (in Fig. 9) has been explored Moakher et al. (2000); Brone and Muzzio (2000); Alexander et al. (2001). However, in this orientation, the flowing layer wall boundary changes shape from a circle to a double cone cross section during tumbling—an effect outside the scope of this study. We consider the double cone rotated about only the -axis (Fig. 9(b)) where the flowing layer retains the same shape during rotation.. Unless otherwise noted, for the double cone geometries considered, the fill level () was chosen such that the free surface extended just to the endwalls, corresponding to the outermost dashed curves for the case in Fig. 9(b). Additionally, the double cone tumbler simulation parameters are identical to those used for the spherical tumbler simulation (Table 1444For ease of simulation, the double cone walls were defined by a monolayer of fixed particles identical to the tumbler particles rather than a smooth surface as used for the spherical tumbler simulations..

In Fig. 10, particle trajectories calculated by integration of the mean velocity field in simulation show that axial drift in a double cone tumbler ( cm,  cm,  cm, and ) is qualitatively similar to that in a spherical tumbler (Fig. 7). A trajectory starting in the top of the flowing layer advances toward the pole and returns to the equator lower in the flowing layer 555The apparent axial drift in the fixed bed (during upward motion in the figure) for rotations 8.75–17.50 may be a consequence of two effects, neither of which represent actual flow in the fixed bed. First, the trajectory may be located at the very bottom of the flowing layer where the apparent reverse flow comes about because the total velocity is the sum of the very slow positive streamwise velocity of the flowing layer superimposed on the slightly faster negative velocity due to the overall rotation of the tumbler system. Second, it is difficult to spatially resolve the velocity field very near the boundary between the flowing layer and the fixed bed since discrete bins used for averaging the velocity field may locally overlap both the flowing layer and the fixed bed. Near the equator, it slowly returns to the top of the flowing layer and repeats the process. The persistence of drift in the double cone tumbler indicates that axially varying wall curvature is not required for axial drift to occur.

Figure 11: Non-dimensional mean axial surface drift of  mm particles in simulation vs. normalized axial position measured at the double conical tumbler wall by particle tracking as in Fig. 8 (,  mm,  rpm). Lines connecting data points emphasize three simulations that are compared directly for wall slope and equator diameter effects. Inset: Maximum dimensionless axial displacement increases linearly with . All simulations have  cm, since the free surface extends from endwall to endwall.

Comparing Fig. 10 for a double cone tumbler with the analogous results for the spherical tumbler in Fig. 7 raises the question of whether or determines the magnitude of axial drift. To address this question, we simulated axial drift in double cone tumblers with different . This was accomplished by holding and constant and selecting equator diameters, of 8, 10, 12, and 14 cm, which correspond to wall slope angles, , of , , , and , respectively. An additional simulation with  cm,  cm, and was performed to further investigate the role of and on the axial drift. For the simulations with ( cm,  cm) and ( cm,  cm), particles filled the tumbler above where the endwall meets the conical section.

Axial drift in the double cone tumbler simulations was measured with the same particle tracking procedure used for the spherical tumbler simulations (Fig. 8). is plotted in Fig. 11 as a function of for all the double cone simulations. Like the spherical tumbler simulations, is 0 at the equator and peaks away from the equator. However, the location of the peak, , occurs closer to the equator in the double cones () than in the sphere (). The maximum of , , increases linearly with (Fig. 11 inset) independent of . The linear scaling of with is expected to fail for small as the bed of particles becomes shallow. In addition, for the double cones with  cm and  cm (Fig. 11) is nearly the same as for the spheres (Fig. 8) with corresponding , suggesting that equator diameter controls in both geometries. Near the pole, the axial drift is more negative than in the spherical tumbler (except for the  cm case), likely due to the relative shallowness of the particle bed in this region compared to that of the spherical tumbler.

To further elucidate the role of in the drift, we focus on the simulations with ( cm,  cm), ( cm,  cm), and ( cm,  cm) (denoted by filled symbols in Fig. 11). The wall angle for ( cm,  cm) and ( cm, ) is the same, and both tumblers exhibit nearly the same rate of decrease of axial drift with axial position from the maximum to the endwall indicating that wall slope determines this behavior. As mentioned earlier, the maximum axial drift is controlled by . Thus, the maximum of axial drift depends on the length of the flowing layer, while the rate of change of with axial position from this maximum toward the pole depends on the wall slope. Moving from the maximum to the equator, decreases continuously to zero. is zero at the equator due to the reflection symmetry of the tumbler about the vertical midplane.

In Fig. 12, which compares axial drift in spherical and double cone tumblers, axial drift is normalized by and axial position is normalized by the axial location of the maximum, . This scaling reveals a separate master curve for the axial drift in each tumbler geometry. While the two curves are similar, the different functional dependence of flowing layer length on axial position results in differences between the two curves. Attempts were made to use a physically motivated scaling for the drift but these scalings did not collapse the data well.

Figure 12: Scaling of axial surface drift in double cone and spherical tumbler simulations (from data in Figs. 8(b) and 11) using normalized axial surface drift, , and normalized axial position, . A parabolic fit for double cone data and a cubic fit for sphere data are shown by solid curves. Plot symbols are as in Fig. 8(b) and Fig. 11, where black filled symbols and ’s are from spherical tumbler simulations and open and gray filled symbols are from double cone tumbler simulations.

V Conclusion

Spherical and double cone tumblers rotating about a horizontal axis exhibit weak three dimensional flow which takes the form of a slow axial drift. The axial drift is 1-2 particle diameters, corresponding to 1%-3% of the tumbler diameter, with each pass through the flowing layer. The drift is caused by asymmetries in the curvature of mean particle trajectories in the flowing layer and results in poleward drift near the free surface and equator directed drift deeper in the flowing layer. Drift is caused by the axial slope of the bounding walls, while the length of the flowing layer determines its magnitude. Axial drift is largely independent of wall friction as both smooth and rough walls yield similar results.

This research was motivated by the commonly used assumption in continuum models of granular flow that flow in any slice perpendicular to the rotation axis is independent of adjacent slices—in other words, the flow in a slice is two-dimensional. Implicitly, this assumption requires that the spanwise flow is negligible and any spanwise motion is primarily diffusive with zero mean. The results here demonstrate that this assumption has limitations, since axial displacements of 1%-3% of the tumbler diameter occur with each pass through the flowing layer. It is unclear if axial drift would be present for particles and tumbler sizes not considered by this study as other forces, such as cohesion and electrostatic forces, may play a larger role. In some circumstances, axial drift might be used to an advantage. For example, it may be possible to offset the surface axial drift away from endwalls observed in cylindrical tumblers with flat endwalls Santomaso et al. (2004); Chen et al. (2009) by using hemispherical (or other non-planar) endwalls which induce drift in the opposite direction.

While this work addresses aspects of axial drift in non-cylindrical tumblers, several details remain poorly understood. Specifically, the mechanism for axial drift in the flowing layer deserves further investigation. It is evident that, in a tumbler with axially sloped walls, poleward drift near the free surface and equator directed drift deeper in the flowing layer satisfy continuity. It is possible that the drift is initiated by particles in the free surface due to axial gradients in the streamwise velocity while particles deep in the flowing layer drift toward the equator to satisfy continuity. Particles do not cross the equatorial plane (except by collisional diffusion) in the spherical or double cone tumblers; this is clearly a consequence of the equator being a plane of symmetry with axial drift having opposite signs on opposite sides of the equator. However, it is not clear what would occur if the plane where the wall slope changes has a discontinuous and non-symmetric wall slope change, say for example, in the case of planes formed by joining a cylindrical tumbler with conical (or hemispherical) ends. It is also not apparent whether the drift occurs in the avalanching regime (at lower tumbler rotational speeds) and in the cataracting regime (curved free surface at higher tumbler rotational speeds). Additionally, it is not immediately clear how the axial variation in particle bed depth impacts the drift. Preliminary experiments at low fill fractions () have shown particle slip impacts particle trajectories in the flowing layer. These and other questions remain to be addressed.

This research was funded by NSF Grant CMMI-1000469. UDO thanks the Mésocentre d’Aix-Marseille Université for providing computing facilities. Travel for RML was supported by the Carnot STAR Institute. Wendy Chan provided useful insight regarding the statistical analysis.
Method (rpm) Passes
Experiment 1 20
Experiment 2 20
Experiment 3 20
Simulation 3 16
Simulation 5 12
Simulation 15 16
Simulation 30 27
Simulation111rough walls 30 69
Table 2: Normalized solid body rotation time, , and normalized flowing layer passage time, , obtained by tracking a  mm tracer particle in a spherical tumbler ( cm,  mm). is the standard deviation.

Appendix A Residence Times

The negligible flowing layer passage time assumption made in Sec. II was verified by experiments and simulations with  mm colored tracer bead in a bed of  mm beads. A larger tracer was used because it remained at the top of the flowing layer and at the tumbler wall when in solid body rotation, which allowed the time spent in solid body rotation to be measured. was used because should be independent of for this fill fraction. The relative flowing layer time, , and the solid body rotation residence time, , were measured for 20 tumbler rotations. Results in Table 2 indicate that is slightly longer than the predicted value of 0.5 for probably due to a small degree of internal slip and slight rearrangement Komatsu et al. (2001); Socie et al. (2005); Arndt et al. (2006) of particles in the fixed bed as the tumbler rotates. is about an order of magnitude less than . It increases slightly with rotation rate, which is probably a result of the thicker flowing layer at these high flow rates. From simulations, the flowing layer at the equator was 1.3 times thicker at 30 rpm than at 3 rpm, resulting in a twofold increase in flowing layer passage time. Flowing layer passage times can also be estimated from experiments in a quasi-2D circular tumbler by Jain et al. Jain et al. (2002) and from the model of Christov et al. Christov et al. (); both methods give values similar to those obtained in tracer experiments (Table 2).

Appendix B Simulation Method

For the DEM simulations, a standard linear-spring and viscous damper force model Cundall and Strack (1979); J. Schäfer et al. (1996); Ristow (2000); Chen et al. (2008) was used to calculate the normal force between two contacting particles: , where and are the particle overlap and the relative velocity () of contacting particles and respectively; is the unit vector in the direction between particles and ; is the reduced mass of the two particles; is the normal stiffness and is the normal damping, where is the collision time and is the restitution coefficient Ristow (2000); Chen et al. (2008). A standard tangential force model Cundall and Strack (1979); J. Schäfer et al. (1996) with elasticity was implemented: , where is the relative tangential velocity of two particles Rapaport (2002), is the tangential stiffness, and is the net tangential displacement after contact is first established at time, . The velocity-Verlet algorithm Ristow (2000); Allen and Tildesley (2002) was used to update the position, orientation, and linear and angular velocity of each particle. Tumbler walls were modeled as both smooth surfaces (smooth walls) and as compositions of bonded particles (rough walls). Both wall conditions have infinite mass for calculation of the collision force between the tumbling particles and the wall. Rough walls were used exclusively for double cone tumbler simulations. Most spherical tumbler simulations used smooth walls, though a few used rough walls for comparison—both wall types produced similar results. To characterize the mean velocity field, the spherical computational domain was divided into cubical bins of width . In most cases, this bin width adequately resolves the flowing layer, as the thickness of the flowing layer typically ranges from  Pignatel et al. (2012). Local flow properties were obtained by averaging values for all particles in each bin every 100 time steps for a total of 15 s of physical time.


Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description