Scaling and dynamics of sphere and disk impact into granular media

Penetration of spherical projectiles into wet granular media

We measure experimentally the penetration depth d of spherical particles into a water-saturated granular medium made of much smaller sand-sized grains. We vary the density, size R, and velocity U of the impacting spheres, and the size δ of the grains in the granular medium. We consider velocities between 7 and 107 m/s, a range not previously addressed, but relevant for impacts produced by volcanic eruptions. We find that d∝R(1/3)δ(1/3)U(2/3). The scaling with velocity is similar to that identified in previous, low-velocity collisions, but it also depends on the size of the grains in the granular medium. We develop a model, consistent with the observed scaling, in which the energy dissipation is dominated by the work required to rearrange grains along a network of force chains in the granular medium.

Granular dynamics during impact

We study the impact of a projectile onto a bed of 3 mm grains immersed in an index-matched fluid. We vary the amount of prestrain on the sample, strengthening the force chains within the system. We find this affects only the prefactor of the linear depth-dependent term in the stopping force. We propose a simple model to account for the strain dependence of this term, owing to increased pressure in the pile. Interestingly, we find that the presence of the fluid does not affect the impact dynamics, suggesting that dynamic friction is not a factor. Using a laser sheet scanning technique to visualize internal grain motion, we measure the trajectory of each grain throughout an impact. Microscopically, our results indicate that weaker initial force chains result in more irreversible, plastic rearrangements, suggesting static friction between grains does play a substantial role in the energy dissipation.

High-velocity drag friction in granular media near the jamming point

Drag friction that acts on a disk in a two-dimensional granular medium is studied at high packing fractions. We concentrate on a high-velocity region, in which the dynamic component of the force, obtained as an average of a strongly fluctuating force, clearly scales with velocity squared. We find that the total force composed of dynamic and static components, as well as its fluctuation, diverges with practically the same exponent as the packing fraction approaches the jamming point. To explain the critical behavior, we propose a simple theory equipped with a diverging length scale, which agrees well with the data and elucidates physical pictures for the divergence.

Experimental investigation into the impact of a liquid droplet onto a granular bed using three-dimensional, time-resolved, particle tracking

An experimental investigation into the interaction that occurs between an impacting water droplet and a granular bed of loose graded sand has been carried out. High-speed imaging, three-dimensional time-resolved particle tracking, and photogrammetric surface profiling have been used to examine individual impact events. The focus of the study is the quantification and trajectory analysis of the particles ejected from the sand bed, along with measurement of the change in bed morphology. The results from the experiments have detailed two distinct mechanisms of particle ejection: the ejection of water-encapsulated particles from the edge of the wetted region and the ejection of dry sand from the periphery of the impact crater. That the process occurs by these two distinct mechanisms has hitherto been unobserved. Presented in the paper are distributions of the particle ejection velocities, angles, and transport distances for both mechanisms. The ejected water-encapsulated particles, which are few in number, are characterized by low ejection angles and high ejection velocities, leading to large transport distances; the ejected dry particles, which are much greater in number, are characterized by high ejection angles and low velocities, leading to lower transport distances. From the particle ejection data, the momentum of the individual ballistic sand particles has been calculated; it was found that only 2% of the water-droplet momentum at impact is transferred to the ballistic sand particles. In addition to the particle tracking, surface profiling of the granular bed postimpact has provided detailed information on its morphology; these data have demonstrated the consistent nature of the craters produced by the impact and suggest that particle agglomerations released from their edges make up about twice the number of particles involved in ballistic ejection. It is estimated that, overall, about 4% of the water-droplet momentum is taken up in particle movement.

Collisional model for granular impact dynamics

When an intruder strikes a granular material from above, the grains exert a stopping force which decelerates and stops the intruder. Many previous studies have used a macroscopic force law, including a drag force which is quadratic in velocity, to characterize the decelerating force on the intruder. However, the microscopic origins of the force-law terms are still a subject of debate. Here, drawing from previous experiments with photoelastic particles, we present a model which describes the velocity-squared force in terms of repeated collisions with clusters of grains. From our high speed photoelastic data, we infer that "clusters" correspond to segments of the strong force network that are excited by the advancing intruder. The model predicts a scaling relation for the velocity-squared drag force that accounts for the intruder shape. Additionally, we show that the collisional model predicts an instability to rotations, which depends on the intruder shape. To test this model, we perform a comprehensive experimental study of the dynamics of two-dimensional granular impacts on beds of photoelastic disks, with different profiles for the leading edge of the intruder. We particularly focus on a simple and useful case for testing shape effects by using triangular-nosed intruders. We show that the collisional model effectively captures the dynamics of intruder deceleration and rotation; i.e., these two dynamical effects can be described as two different manifestations of the same grain-scale physical processes.

Depth-dependent resistance of granular media to vertical penetration

We measure the quasistatic friction force acting on intruders moving downwards into a granular medium. By utilizing different intruder geometries, we demonstrate that the force acts locally normal to the intruder surface. By altering the hydrostatic loading of grain contacts by a sub-fluidizing airflow through the bed, we demonstrate that the relevant frictional contacts are loaded by gravity rather than by the motion of the intruder itself. Lastly, by measuring the final penetration depth versus airspeed and using an earlier result for inertial drag, we demonstrate that the same quasistatic friction force acts during impact. Altogether this force is set by a friction coefficient, hydrostatic pressure, projectile size and shape, and a dimensionless proportionality constant. The latter is the same in nearly all experiments, and is surprisingly greater than one.

Lift and drag in intruders moving through hydrostatic granular media at high speeds

Recently, experiments showed that forces on intruders dragged horizontally through dense, hydrostatic granular packings mainly depend on the local surface orientation and can be seen as the sum of the forces exerted on small surface elements. In order to understand such forces more deeply, we perform a two-dimensional soft-sphere molecular dynamics simulation, on a similar setup, of an intruder dragged through a 50-50 bi-disperse granular packing, with diameters 0.30 and 0.34 cm. We measure, for both circular and half-circle shapes, the forces parallel (drag) and perpendicular (lift) to the drag direction as functions of the drag speed, with V=10.3-309 cm/s, and intruder depths, with D=3.75-37.5 cm. The drag forces on an intruder monotonically increase with V and D, and are larger for the circle. However, the lift force does not depend monotonically on V and D, and this relationship is affected by the shape of the intruder. The vertical force was negative for the half-circle, but for a small range of V and D, we measure positive lift. We find no sign change for the lift on the circle, which is always positive. The explanation for the nonmonotonic dependence is related to the decrease in contacts on the intruder as V increases. This is qualitatively similar to supersonic flow detachment from an obstacle. The detachment picture is supported by simulation measurements of the velocity field around the intruder and force profiles measured on its surface.

Ground fluidization promotes rapid running of a lightweight robot

We study the locomotor mechanics of a small, lightweight robot (DynaRoACH, 10 cm, 25 g) which can move on a granular substrate of 3 mm diameter glass particles at speeds up to 5 body length/s, approaching the performance of certain desert-dwelling animals. To reveal how the robot achieves this performance, we used high-speed imaging to capture its kinematics, and developed a numerical multi-body simulation of the robot coupled to an experimentally validated simulation of the granular medium. Average speeds measured in experiment and simulation agreed well, and increased nonlinearly with stride frequency, reflecting a change in propulsion mode. At low frequencies, the robot used a quasi-static “rotary walking” mode, in which the substrate yielded as legs penetrated and then solidified once vertical force balance was achieved. At high frequencies the robot propelled itself using the speed-dependent fluid-like inertial response of the material. The simulation allows variation of parameters which are inconvenient to modify in experiment, and thus gives insight into how substrate and robot properties change performance. Our study reveals how lightweight animals can achieve high performance on granular substrates; such insights can advance the design and control of robots in deformable terrains.

Penetration of projectiles into granular targets

Energetic collisions of subatomic particles with fixed or moving targets have been very valuable to penetrate into the mysteries of nature. But the mysteries are quite intriguing when projectiles and targets are macroscopically immense. We know that countless debris wandering in space impacted (and still do) large asteroids, moons and planets; and that millions of craters on their surfaces are traces of such collisions. By classifying and studying the morphology of such craters, geologists and astrophysicists obtain important clues to understand the origin and evolution of the Solar System. This review surveys knowledge about crater phenomena in the planetary science context, avoiding detailed descriptions already found in excellent papers on the subject. Then, it examines the most important results reported in the literature related to impact and penetration phenomena in granular targets obtained by doing simple experiments. The main goal is to discern whether both schools, one that takes into account the right ingredients (planetary bodies and very high energies) but cannot physically reproduce the collisions, and the other that easily carries out the collisions but uses laboratory ingredients (small projectiles and low energies), can arrive at a synergistic intersection point.

Drag force scaling for penetration into granular media

Impact dynamics is measured for spherical and cylindrical projectiles of many different densities dropped onto a variety non-cohesive granular media. The results are analyzed in terms of the material-dependent scaling of the inertial and frictional drag contributions to the total stopping force. The inertial drag force scales similar to that in fluids, except that it depends on the internal friction coefficient. The frictional drag force scales as the square-root of the density of granular medium and projectile, and hence cannot be explained by the combination of granular hydrostatic pressure and Coulomb friction law. The combined results provide an explanation for the previously observed penetration depth scaling.

Experimental velocity fields and forces for a cylinder penetrating into a granular medium

We present here a detailed granular flow characterization together with force measurements for the quasi-bidimensional situation of a horizontal cylinder penetrating vertically at a constant velocity in dry granular matter between two parallel glass walls. In the velocity range studied here, the drag force on the cylinder does not depend on the velocity V(0) and is mainly proportional to the cylinder diameter d. While the force on the cylinder increases with its penetration depth, the granular velocity profile around the cylinder is found to be stationary with fluctuations around a mean value leading to the granular temperature profile. Both mean velocity profile and temperature profile exhibit strong localization near the cylinder. The mean flow perturbation induced by the cylinder decreases exponentially away from the cylinder on a characteristic length λ that is mainly governed by the cylinder diameter for a large enough cylinder/grain size ratio d/d(g): λ~d/4+2d(g). The granular temperature exhibits a constant plateau value T(0) in a thin layer close to the cylinder of extension δ(T(0))~λ/2 and decays exponentially far away with a characteristic length λ(T) of a few grain diameters (λ(T)~3d(g)). The granular temperature plateau T(0) that scales as V(0)(2)d(g)/d is created by the flow itself from the balance between the "granular heat" production by the shear rate V(0)/λ over δ(T(0)) close to the cylinder and the granular dissipation far away.

Particle scale dynamics in granular impact

We perform an experimental study of granular impact, where intruders strike 2D beds of photoelastic disks from above. High-speed video captures the intruder dynamics and the local granular force response, allowing investigation of grain-scale mechanisms in this process. We observe rich acoustic behavior at the leading edge of the intruder, strongly fluctuating in space and time, and we show that this acoustic activity controls the intruder deceleration, including large force fluctuations at short time scales. The average intruder dynamics match previous studies using empirical force laws, suggesting a new microscopic picture, where acoustic energy is carried away and dissipated.

Influence of particle rotation on the oblique penetration in granular media

The rotation of a particle has significant influence on the dynamic response of a granular bed subjected to the oblique impact of a spherical projectile. Based on the discrete element method, the dynamical behavior of two-dimensional granular media impacted obliquely by rotating particles has been examined in this work, especially for the influence of rotational angular velocity on its penetration depth. The simulations show that the incident rotational velocity will not only act on the characteristics of velocity distribution of the bed after impacting, but also influence the trajectory of the projectile qualitatively. For low angular velocities, particle rotations will significantly increase the vertical penetration depths, while the different directions of rotation will exert opposite effects on the horizontal penetration depths. In addition, the influence of particle rotation on its penetration depth will be enhanced with increasing angular velocity, but such effect will reach an asymptotic plateau for sufficiently large angular velocities. This indicates that the angular velocity has an obvious criticality. Furthermore, the variation of critical angular velocity may be linear with the impact velocity and square with the impacting angle, approximately. Finally, the influence of the initial particle rotation on the scaling law of penetration depth is also considered, and we find that the linear scaling with impact velocity is still applicable for most impact conditions.

Using computational and mechanical models to study animal locomotion

Recent advances in computational methods have made realistic large-scale simulations of animal locomotion possible. This has resulted in numerous mathematical and computational studies of animal movement through fluids and over substrates with the purpose of better understanding organisms' performance and improving the design of vehicles moving through air and water and on land. This work has also motivated the development of improved numerical methods and modeling techniques for animal locomotion that is characterized by the interactions of fluids, substrates, and structures. Despite the large body of recent work in this area, the application of mathematical and numerical methods to improve our understanding of organisms in the context of their environment and physiology has remained relatively unexplored. Nature has evolved a wide variety of fascinating mechanisms of locomotion that exploit the properties of complex materials and fluids, but only recently are the mathematical, computational, and robotic tools available to rigorously compare the relative advantages and disadvantages of different methods of locomotion in variable environments. Similarly, advances in computational physiology have only recently allowed investigators to explore how changes at the molecular, cellular, and tissue levels might lead to changes in performance at the organismal level. In this article, we highlight recent examples of how computational, mathematical, and experimental tools can be combined to ultimately answer the questions posed in one of the grand challenges in organismal biology: "Integrating living and physical systems."

Sphere impact and penetration into wet sand

We present experimental results for the penetration of a solid sphere when released onto wet sand. We show, by measuring the final penetration depth, that the cohesion induced by the water can result in either a deeper or shallower penetration for a given release height compared to dry granular material. Thus the presence of water can either lubricate or stiffen the granular material. By assuming the shear rate is proportional to the impact velocity and using the depth-averaged stopping force in calculating the shear stress, we derive effective viscosities for the wet granular materials.

Multi-functional foot use during running in the zebra-tailed lizard (Callisaurus draconoides)

A diversity of animals that run on solid, level, flat, non-slip surfaces appear to bounce on their legs; elastic elements in the limbs can store and return energy during each step. The mechanics and energetics of running in natural terrain, particularly on surfaces that can yield and flow under stress, is less understood. The zebra-tailed lizard (Callisaurus draconoides), a small desert generalist with a large, elongate, tendinous hind foot, runs rapidly across a variety of natural substrates. We use high speed video to obtain detailed three-dimensional running kinematics on solid and granular surfaces to reveal how leg, foot, and substrate mechanics contribute to its high locomotor performance. Running at ~10 body length/s (~1 m/s), the center of mass oscillates like a spring-mass system on both substrates, with only 15% reduction in stride length on the granular surface. On the solid surface, a strut-spring model of the hind limb reveals that the hind foot saves about 40% of the mechanical work needed per step, significant for the lizard's small size. On the granular surface, a penetration force model and hypothesized subsurface foot rotation indicates that the hind foot paddles through fluidized granular medium, and that the energy lost during irreversible deformation of the substrate does not differ from the reduction in the mechanical energy of the center of mass. The upper hind leg muscles must perform three times as much mechanical work on the granular surface as on the solid surface to compensate for the greater energy lost within the foot and to the substrate.

Microstructure evolution during impact on granular matter

We study the impact of an intruder on a dense granular material. The process of impact and interaction between the intruder and the granular particles is modeled using discrete element simulations in two spatial dimensions. In the first part of the paper we discuss how the intruder's dynamics depends on (1) the intruder's properties, including its size, shape and composition, (2) the properties of the grains, including friction, polydispersity, structural order, and elasticity, and (3) the properties of the system, including its size and gravitational field. It is found that polydispersity and related structural order, and frictional properties of the granular particles, play a crucial role in determining impact dynamics. In the second part of the paper we consider the response of the granular system itself. We discuss the force networks that develop, including their topological evolution. The influence of friction and structural order on force propagation, including the transition from hyperbolic-like to elastic-like behavior is discussed, as well as the affine and nonaffine components of the grain dynamics. Several broad observations include the following: tangential forces between granular particles are found to play a crucial role in determining impact dynamics; both force networks and particle dynamics are correlated with the dynamics of the intruder itself.

The scaling and dynamics of a projectile obliquely impacting a granular medium

In this paper, the dynamics of a spherical projectile obliquely impacting into a two-dimensional granular bed is numerically investigated using the discrete element method. The influences of projectile's initial velocities and impacting angles are mainly considered. Numerical results show that the relationship between the final penetration depth and the initial impact velocity is very similar to that in the vertical-impact case. However, the dependence of the stopping time on the impact velocity of the projectile exhibits critical characteristics at different impact angles: the stopping time approximately increases linearly with the impact velocity for small impact angles but decreases in an exponential form for larger impact angles, which demonstrates the existence of two different regimes at low and high impact angles. When the impact angle is regarded as a parametric variable, a phenomenological force model at large impact angles is eventually proposed based on the simulation results, which can accurately describe the nature of the resistance force exerted on the projectile by the granular medium at different impact angels during the whole oblique-impact process. The degenerate model agrees well with the existing experimental results in the vertical-impact cases.

Dense granular flow around a penetrating object: experiment and hydrodynamic model

We present in this Letter experimental results on the bidimensional flow field around a cylinder penetrating into dense granular matter, together with drag force measurements. A hydrodynamic model based on extended kinetic theory for dense granular flow reproduces well the flow localization close to the cylinder and the corresponding scalings of the drag force, which is found to not depend on velocity, but linearly on the pressure and on the cylinder diameter and weakly on the grain size. Such a regime is found to be valid at a low enough "granular" Reynolds number.

Undulatory swimming in sand: experimental and simulation studies of a robotic sandfish

A previous study of a sand-swimming lizard, the sandfish, revealed that it swims within granular media at speeds up to 0.4 body-lengths/cycle using body undulations (approximately a single period sinusoidal traveling wave) without limb use. Inspired by the organism, we develop a numerical model of a robot swimming in a simulated granular medium to guide the design of a physical device. Both in simulation and experiment the robot swims limblessly subsurface at speeds up to 0.3 body-lengths/cycle and, like the animal, increases its speed by increasing its oscillation frequency. The performance of the robot measured in terms of its wave efficiency η, the ratio of its forward speed to wave speed, is 0.34 ± 0.02 , within 8% of the simulation prediction. Both in simulation and experiment, η increases with increasing particle—particle friction but decreases with increasing body—particle friction. On a flat, rigid surface the robot fails to move forward, as expected, due to the frictional isotropy between the interacting surfaces. However, the surface and subsurface performance of the robot on low friction particles are comparable. Our work provides a validated simulation tool and the design of a robot that can move on or within yielding terrestrial substrates.

Drag induced lift in granular media

Laboratory experiments and numerical simulation reveal that a submerged intruder dragged horizontally at a constant velocity within a granular medium experiences a lift force whose sign and magnitude depend on the intruder shape. Comparing the stress on a flat plate at varied inclination angle with the local surface stress on the intruders at regions with the same orientation demonstrates that intruder lift forces are well approximated as the sum of contributions from flat-plate elements. The plate stress is deduced from the force balance on the flowing media near the plate.

Granular impact and the critical packing state

Impact dynamics during collisions of spheres with granular media reveal a pronounced and nontrivial dependence on volume fraction ϕ. Postimpact crater morphology identifies the critical packing state ϕcps, where sheared grains neither dilate nor consolidate, and indicates an associated change in spatial response. Current phenomenological models fail to capture the observed impact force for most ϕ; only near ϕcps is force separable into additive terms linear in depth and quadratic in velocity. At fixed depth the quadratic drag coefficient decreases (increases) with depth for ϕ<ϕcps (ϕ>ϕcps). At fixed low velocity, depth dependence of force shows a Janssen-type exponential response with a length scale that decreases with increasing ϕ and is nearly constant for ϕ>ϕcps.

Morphology scaling of drop impact onto a granular layer

We investigate the impact of a free-falling water drop onto a granular layer. First, we constructed a phase diagram of crater shapes with two control parameters, impact speed and grain size. A low-speed impact makes a deeper cylindrical crater in a fluffy granular target. After high-speed impacts, we observed a convex bump higher than the initial surface level instead of a crater. The inner ring can be also observed in a medium impact speed regime. Quantitatively, we found a scaling law for a crater radius with a dimensionless number consisting of impact speed and density ratio between the bulk granular layer and water drop. This scaling demonstrates that the water drop deformation is crucial to understanding the crater morphology.

Dynamics of grain ejection by sphere impact on a granular bed

The dynamics of grain ejection consecutive to a sphere impacting a granular material is investigated experimentally and the variations of the characteristics of grain ejection with the control parameters are quantitatively studied. The time evolution of the corona formed by the ejected grains is reported, mainly in terms of its diameter and height, and favorably compared with a simple ballistic model. A key characteristic of the granular corona is that the angle formed by its edge with the horizontal granular surface remains constant during the ejection process, which again can be reproduced by the ballistic model. The number and the kinetic energy of the ejected grains are evaluated and allow for the calculation of an effective restitution coefficient characterizing the complex collision process between the impacting sphere and the fine granular target. The effective restitution coefficient is found to be constant when varying the control parameters.

From the Cover: Sensitive dependence of the motion of a legged robot on granular media

Legged locomotion on flowing ground (e.g., granular media) is unlike locomotion on hard ground because feet experience both solid- and fluid-like forces during surface penetration. Recent bioinspired legged robots display speed relative to body size on hard ground comparable with high-performing organisms like cockroaches but suffer significant performance loss on flowing materials like sand. In laboratory experiments, we study the performance (speed) of a small (2.3 kg) 6-legged robot, SandBot, as it runs on a bed of granular media (1-mm poppy seeds). For an alternating tripod gait on the granular bed, standard gait control parameters achieve speeds at best 2 orders of magnitude smaller than the 2 body lengths/s (approximately 60 cm/s) for motion on hard ground. However, empirical adjustment of these control parameters away from the hard ground settings restores good performance, yielding top speeds of 30 cm/s. Robot speed depends sensitively on the packing fraction phi and the limb frequency omega, and a dramatic transition from rotary walking to slow swimming occurs when phi becomes small enough and/or omega large enough. We propose a kinematic model of the rotary walking mode based on generic features of penetration and slip of a curved limb in granular media. The model captures the dependence of robot speed on limb frequency and the transition between walking and swimming modes but highlights the need for a deeper understanding of the physics of granular media.

Projectile interactions in granular impact cratering

We present evidence for the interactions between a ball and the container boundaries, as well as between two balls, that are mediated by the granular medium during impact cratering. The presence of the bottom boundary affects the final penetration depth only for low drop heights with shallow filling, in which case, surprisingly, the penetration becomes deeper. By contrast the presence of the sidewall causes less penetration and also an effective repulsion. Repulsion is also found for two balls dropped side by side.

Influence of confinement on granular penetration by impact

We study experimentally the influence of confinement on the penetration depth of impacting spheres into a granular medium contained in a finite cylindrical vessel. The presence of close lateral walls reduces the penetration depth, and the characteristic distance for these lateral wall effects is found to be of the order of one sphere diameter. The influence of the bottom wall is found to have a much shorter range.