Impact craters in loose granular media

Drag-force regimes in granular impact

We study the penetration dynamics of a projectile incident normally on a substrate comprising of smaller granular particles in three-dimensions using the discrete element method. Scaling of the penetration depth is consistent with experimental observations for small velocity impacts. Our studies are consistent with the observation that the normal or drag force experienced by the penetrating grain obeys the generalized Poncelet law, which has been extensively invoked in understanding the drag force in the recent experimental data. We find that the normal force experienced by the projectile consists of position and kinetic-energy-dependent pieces. Three different penetration regimes are identified in our studies for low-impact velocities. The first two regimes are observed immediately after the impact and in the early penetration stage, respectively, during which the drag force is seen to depend on the kinetic energy. The depth dependence of the drag force becomes significant in the third regime when the projectile is moving slowly and is partially immersed in the substrate. These regimes relate to the different configurations of the bed: the initial loose surface packed state, fluidized bed below the region of impact, and the state after the crater formation commences.

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.

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.

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.

Infinite penetration of a projectile into a granular medium

An object falling in a fluid reaches a terminal velocity when the drag force and its weight are balanced. Contrastingly, an object impacting into a granular medium rapidly dissipates all its energy and comes to rest always at a shallow depth. Here we study, experimentally and theoretically, the penetration dynamics of a projectile in a very long silo filled with expanded polystyrene particles. We discovered that, above a critical mass, the projectile reaches a terminal velocity and, therefore, an endless penetration.

Collapse of a rectangular well in a quasi-two-dimensional granular bed

We report on an experimental and numerical study of the collapse under gravity of a rectangular well in a quasi-two-dimensional granular bed. For comparison, we also perform experiments on the collapse of a single vertical step. Experiments are conducted in a vertical Hele-Shaw cell, which allows the flow to be recorded from the side using high-speed video. If the rectangular well is sufficiently narrow, the collapsing sidewalls collide at the center of the well and the dynamics of the collapse are dependant on the aspect ratio of the initial well. We follow the evolution of the free surface from the video images, and use particle image velocimetry to determine the subsurface velocity field. From these data, the potential and kinetic energy of the system are calculated. We observe two stages to the collapse flow: an initial gravity-dominated stage, during which the kinetic energy increases, and a later dissipation-dominated phase during which the kinetic energy decreases. We find that although both the width and depth of the depression that remains after the well has collapsed depend on the initial aspect ratio, the surface profiles are self-similar; that is, the shape of the final profile is independent of the aspect ratio of the initial well. We model the collapse of the well using a depth-averaged continuum model with basal friction and with a discrete element model. Both models give results which agree well with experiment. The discrete element model indicates that friction between the particles is the most important source of dissipation over the course of the collapse.

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.

Dynamics of crater formations in immersed granular materials

We report the formation of a crater at the free surface of an immersed granular bed, locally crossed by an ascending gas flow. In two dimensions, the crater consists of two piles which develop around the location of the gas emission. We observe that the typical size of the crater increases logarithmically with time, independently of the gas emission dynamics. We describe the related granular flows and give an account of the influence of the experimental parameters, especially of the grain size and of the gas flow.

Shape of impact craters in granular media

We present the results of experiments studying the shape of craters formed by the normal impact of a solid spherical projectile into a deep noncohesive granular bed at low energies. The resultant impact crater surfaces are accurately digitized using laser profilometry, allowing for the detailed investigation of the crater shape. We find that these impact craters are very nearly hyperbolic in profile. Crater radii and depths are dependent on impact energy, as well as the projectile density and size. The precise crater shape is a function of the crater aspect ratio. While the dimensions of the crater are highly dependent on the impact energy, we show that the energy required to excavate the crater is only a tiny fraction (0.1%-0.5%) of the kinetic energy of the projectile.

Dynamics of impact cratering in shallow sand layers

When a solid sphere impacts a shallow layer of sand deposited on a solid surface, a crater can be obtained. The dynamics of the opening of the crater can be followed accurately. During this opening, the radius of the crater can be conveniently modeled by an exponential saturation with a well-defined time constant. The crater then closes up partially once the opening phase is over as the sand avalanches down the slope of the crater. We here present a detailed study of the full dynamics of the crater formation as well as the dynamics of the corrola formed during this process. A simple model accounts for most of our observations.

Dynamics of shallow impact cratering

We present data for the time dependence of wooden spheres penetrating into a loose noncohesive packing of glass beads. The stopping time is a factor of 3 longer than the time d/v0 needed to travel the total penetration distance d at the impact speed v0. The acceleration decreases monotonically throughout the impact. These kinematics are modeled by a position- and velocity-dependent stopping force that is constrained to reproduce prior observations for the scaling of the penetration depth with the total drop distance.

Penetration depth for shallow impact cratering

We present data for the penetration of a variety of spheres, dropped from rest, into a loose noncohesive granular medium. We improve upon earlier work [J. S. Uehara, Phys. Rev. Lett. 90, 194301 (2003)] in three regards. First, we explore the behavior vs sphere diameter and density more systematically, by holding one of these parameters constant while varying the other. Second, we prepare the granular medium more reproducibly and, third, we measure the penetration depth more accurately. The new data support the previous conclusion that the penetration depth is proportional to the 1/2 power of sphere density, the 2/3 power of sphere diameter, and the 1/3 power of total drop distance.