Search

The Discrete Element Method is widely employed for simulating granular flows, but conventional integration techniques may produce unphysical results for simulations with static friction when particle size ratios exceed $R \approx 3$. These inaccuracies arise because some variables in the velocity-Verlet algorithm are calculated at the half-timestep, while others are computed at the full timestep. To correct this, we develop an improved velocity-Verlet integration algorithm to ensure physically accurate outcomes up to the largest size ratios examined ($R=100$). The implementation of this improved integration method within the LAMMPS framework is detailed, and its effectiveness is validated through a simple three-particle test case and a more general example of granular flow in mixtures with large size-ratios, for which we provide general guidelines for selecting simulation parameters and accurately modeling inelasticity in large particle size-ratio simulations.

Predicting particle segregation has remained challenging due to the lack of a general model for the segregation velocity that is applicable across a range of granular flow geometries. Here, a segregation velocity model for dense granular flows is developed by exploiting momentum balance and recent advances in particle-scale modelling of the segregation driving and drag forces over a wide range of particle concentrations, size and density ratios, and flow conditions. This model is shown to correctly predict particle segregation velocity in a diverse set of idealized and natural granular flow geometries simulated using the discrete element method. When incorporated in the well-established advection-diffusion-segregation formulation, the model has the potential to accurately capture segregation phenomena in many relevant industrial application and geophysical settings.

Lift and drag forces on moving intruders in granular materials are of fundamental interest. While the drag force on an intruder in granular flow has been studied, the few studies characterizing the lift force explore a relatively limited range of parameters. Here we use discrete element method (DEM) simulations to measure the lift force, $F_\mathrm{L}$, on a spherical intruder in a uniformly sheared bed of smaller spheres for a range of intruder slip velocities, $u_\mathrm{s}$, relative to the unperturbed flow. In what at first appears as a puzzling result, $F_\mathrm{L}$ in granular shear flow acts in the opposite direction to the Saffman lift force on a sphere in a sheared fluid at low $u_\mathrm{s}$, reaches a maximum value, and then decreases, eventually reversing direction and becoming comparable to $F_\mathrm{L}$ for a fluid. This non-monotonic response holds over a range of flow conditions, and the $F_\mathrm{L}$ versus $u_\mathrm{s}$ data can be collapsed by scaling both quantities using the particle sizes, shear rate, and overburden pressure. Analogous fluid simulations demonstrate that the flow field around the intruder particle is similar in the granular and fluid cases. However, the shear stress acting on the intruder in a granular shear flow is much less than that in a fluid shear flow. This difference, combined with a void region behind the intruder in granular flow, which alters the pressure and shear stress on the trailing side of the intruder, significantly changes the lift-force inducing stresses acting on the intruder between the granular and fluid cases., Comment: 27 pages, 17 figures

Static granular packings play a central role in numerous industrial applications and natural settings. In these situations, fluid or fine particle flow through a bed of static particles is heavily influenced by the narrowest passage connecting the pores of the packing, commonly referred to as pore throats or constrictions. Existing studies predominantly assume monodisperse rigid particles, but this is an oversimplification of the problem. In this work, we illustrate the connection between pore throat size, polydispersity, and particle deformation. Simple analytical expressions are provided to link these properties of the packing, followed by examples from Discrete Element Method (DEM) simulations of fine particle percolation demonstrating the impact of polydispersity and particle deformation. Our intent is to emphasize the substantial impact of polydispersity and particle deformation on constriction size, underscoring the importance of accounting for these effects in particle transport in granular packings.

Legged robot locomotion is hindered by a mismatch between applications where legs can outperform wheels or treads, most of which feature deformable substrates, and existing tools for planning and control, most of which assume flat, rigid substrates. In this study we focus on the ramifications of plastic terrain deformation on the hop-to-hop energy dynamics of a spring-legged monopedal hopping robot animated by a switched-compliance energy injection controller. From this deliberately simple robot-terrain template, we derive a hop-to-hop energy return map, and we use physical experiments and simulations to validate the hop-to-hop energy map for a real robot hopping on a real deformable substrate. The dynamical properties (fixed points, eigenvalues, basins of attraction) of this map provide insights into efficient, responsive, and robust locomotion on deformable terrain. Specifically, we identify constant-fixed-point surfaces in a controller parameter space that suggest it is possible to tune control parameters for efficiency or responsiveness while targeting a desired gait energy level. We also identify conditions under which fixed points of the energy map are globally stable, and we further characterize the basins of attraction of fixed points when these conditions are not satisfied. We conclude by discussing the implications of this hop-to-hop energy map for planning, control, and estimation for efficient, agile, and robust legged locomotion on deformable terrain., Comment: 20 pages, 13 figures, submitted to IEEE Transactions on Robotics

We study the percolation velocity, $v_p,$ of a fine spherical particle in a sheared large-particle bed under gravity using discrete element method simulations for large-to-fine particle diameter ratios, $R=d/d_f,$ below and above the free-sifting threshold, $R_t\approx6.5.$ For $RR_t$, $v_p$ is constant at low $\dot\gamma$ but decreases toward zero at higher shear rates due to fine-particle excitation.

Particle segregation in dense flowing size-disperse granular mixtures is driven by gravity and shear, but predicting the associated segregation force due to both effects has remained an unresolved challenge. Here, a model of the combined gravity- and kinematics-induced segregation force on a single intruder particle is integrated with a model of the concentration dependence of the gravity-induced segregation force. The result is a general model of the net particle segregation force in flowing size-bidisperse granular mixtures. Using discrete element method simulations for comparison, the model correctly predicts the segregation force for a variety of mixture concentrations and flow conditions in both idealized and natural shear flows.

In dense flowing bidisperse particle mixtures varying in size or density alone, smaller particles sink (driven by percolation) and lighter particles rise (driven by buoyancy). But when the particle species differ from each other in both size and density, percolation and buoyancy can either enhance (large/light and small/heavy) or oppose (large/heavy and small/light) each other. In the latter case, a local equilibrium condition can exist in which the two segregation mechanisms balance and particles remain mixed: this allows the design of minimally-segregating mixtures by specifying particle size ratio, density ratio, and mixture concentration. Using experimentally validated DEM simulations, we show that mixtures specified by the methodology remain relatively well-mixed in the thin rapid surface flows characteristic of heaps and tumblers commonly used in industry. Furthermore, minimally-segregating particle mixtures prepared in a fully segregated state in a tumbler mix over time and eventually reach a state of nearly uniform species concentration., Comment: 19 pages, 17 figures, submitted to AIChE Journal

We present a continuum approach to model segregation of size-bidisperse granular materials in unsteady bounded heap flow as a prototype for modeling segregation in other time varying flows. In experiments, a periodically modulated feed rate produces stratified segregation like that which occurs due to intermittent avalanching, except with greater layer-uniformity and higher average feed rates. Using an advection-diffusion-segregation equation and characterizing transient changes in deposition and erosion after a feed rate change, we model stratification for varying feed rates and periods. Feed rate modulation in heap flows can create well-segregated layers, which effectively mix the deposited material normal to the free surface at lengths greater than the combined layer-thickness. This mitigates the strong streamwise segregation that would otherwise occur at larger particle-size ratios and equivalent steady feed rates and can significantly reduce concentration variation during hopper discharge. Coupling segregation, deposition and erosion is challenging but has many potential applications.