Microstructure evolution during impact on granular matter

Roles of packing fraction, microscopic friction, and projectile spin in cratering by impact

From small seeds falling from trees to asteroids colliding with planets and moons, the impact of projectiles onto granular targets occurs in nature at different scales. In this paper, we investigate open questions in the mechanics of granular cratering, in particular, the forces acting on the projectile and the roles of granular packing, grain-grain friction, and projectile spin. For that, we carried out discrete element method computations of the impact of solid projectiles on a cohesionless granular medium, where we varied the projectile and grain properties (diameter, density, friction, and packing fraction) for different available energies (within relatively small values). We found that a denser region forms below the projectile, pushing it back and causing its rebound by the end of its motion, and that solid friction affects considerably the crater morphology. Besides, we show that the penetration length increases with the initial spin of the projectile, and that differences in initial packing fractions can engender the diversity of scaling laws found in the literature. Finally, we propose an ad hoc scaling that collapsed our data for the penetration length and can perhaps unify existing correlations. Our results provide new insights into the formation of craters in granular matter.

Experimental investigations of crater formation as a result of high-velocity impacts on sand bed

The formation of craters is an important issue in the investigations of the surface of the earth and other planets. The aim of the study was to check whether the different textures of sand beds affect the size and dynamics of the formation of craters and ejection curtain after high-velocity impacts. The experiments were conducted using an aluminium impactor at two impact speeds (~700 and ~1300 m∙s^-1) and a sand bed composed of either a broad range of sizes (<2.0 mm) or any of the three fractions obtained from it (<0.5, 0.5–1, 1–2 mm). The diameters, depths, wall slope, and rim heights of the resulting craters were measured. The ejecta curtain was characterized by the inclination angle of walls, base diameter, and expansion velocity. The mass of the transferred material and the depth of the impactor penetration were also determined. Additionally, the results were used to calculate dimensionless parameters commonly considered in crater studies (π_V, π₂ and α). The texture of the sand most clearly influenced the diameters of the craters, its effect could also be seen in the case of the distance covered by the ejected material. This information appears to be relevant for future research, providing some rationale to help assess in which aspects of the phenomenon the texture may be important.

Sink versus tilt penetration into shaken dry granular matter: The role of the foundation

We study the behavior of cylindrical objects as they sink into a dry granular bed fluidized due to lateral oscillations. Somewhat unexpectedly, we have found that, within a large range of lateral shaking powers, cylinders with flat bottoms sink vertically, while those with a "foundation" consisting of a shallow ring attached to their bottom, tilt besides sinking. The latter scenario seems to dominate independently from the nature of the foundation when strong enough lateral vibrations are applied. We are able to explain the observed behavior by quasi-2D numerical simulations, which also demonstrate the influence of the intruder's aspect ratio. The vertical sink dynamics is explained with the help of a Newtonian equation of motion for the intruder. Our findings may shed light on the behavior of buildings and other manmade structures during earthquakes.

Study on Meso-Structure Evolution in Granular Matters Based on the Contact Loop Recognition and Determination Technique

On the mesoscopic scale, granular matter is tessellated into contact loops by a contact network. The stability of granular matter is highly dependent on the evolution of contact loops, including the number and area evolutions of contact loops with different geometric shapes (which can reflect the mechanical variables in the macroscale). For the features of numerous loops with complex geometry shapes in contact network images, a contact loop recognition and determination technique was developed in this study. Then, numerical biaxial compression tests were performed by the discrete element method (DEM) to investigate how the meso-structural indexes evolve along with the macro-mechanical indexes. The results show that the proposed Q-Y algorithm is effective in determining the geometric types of contact loops from contact network images. The evolution of contact loops is most active in the hardening stage, during which the number percentages of L₃ (loops with three sides) and L₆₊ (loops with six or more sides) show opposite evolution patterns. For the area percentage, only L₆₊ increases while others decrease. Considering the meso-structural indexes (number percentage and area percentage of loops) are sensitive to the change of macro-mechanical indexes (deviatoric stress, axial strain, and volumetric strain) in the hardening stage. Multivariate models were established to build a bridge between the meso-structure and the macro-mechanics.

Intrusion into Granular Media Beyond the Quasistatic Regime

The drag force exerted on an object intruding into granular media is typically assumed to arise from additive velocity and depth dependent contributions. We test this with intrusion experiments and molecular dynamics simulations at constant speed over four orders of magnitude, well beyond the quasistatic regime. For a vertical cylindrical rod we find velocity dependence only right after impact, followed by a crossover to a common, purely depth-dependent behavior for all intrusion speeds. The crossover is set by the timescale for material, forced to well up at impact, to subsequently settle under gravity. These results challenge current models of granular drag.

Collision-based understanding of the force law in granular impact dynamics

We study the stopping force felt by an intruder impacting onto a granular medium. Variations in the shape of the intruder can influence the penetration depth by changing the inertial drag. We find this observed correlation can be explained by associating the velocity-dependent inertial drag to the energy dissipation that occurs through intermittent collisions of force-chain-like clusters, the mean behavior of which can be statistically described. In consequence, the stopping force can be captured through a proposed collisional model with good accuracy, and the observed impact dynamics data can be reproduced quantitatively.

Penetration in bimodal, polydisperse granular material

We investigate the impact penetration of spheres into granular media which are compositions of two discrete size ranges, thus creating a polydisperse bimodal material. We examine the penetration depth as a function of the composition (volume fractions of the respective sizes) and impact speed. Penetration depths were found to vary between δ=0.5D_{0} and δ=7D_{0}, which, for mono-modal media only, could be correlated in terms of the total drop height, H=h+δ, as in previous studies, by incorporating correction factors for the packing fraction. Bimodal data can only be collapsed by deriving a critical packing fraction for each mass fraction. The data for the mixed grains exhibit a surprising lubricating effect, which was most significant when the finest grains [d_{s}∼O(30) μm] were added to the larger particles [d_{l}∼O(200-500) μm], with a size ratio, ε=d_{l}/d_{s}, larger than 3 and mass fractions over 25%, despite the increased packing fraction. We postulate that the small grains get between the large grains and reduce their intergrain friction, only when their mass fraction is sufficiently large to prevent them from simply rattling in the voids between the large particles. This is supported by our experimental observations of the largest lubrication effect produced by adding small glass beads to a bed of large sand particles with rough surfaces.

Steady flow dynamics during granular impact

We study experimentally and computationally the dynamics of granular flow during impacts where intruders strike a collection of disks from above. In the regime where granular force dynamics are much more rapid than the intruder motion, we find that the particle flow near the intruder is proportional to the instantaneous intruder speed; it is essentially constant when normalized by that speed. The granular flow is nearly divergence free and remains in balance with the intruder, despite the latter's rapid deceleration. Simulations indicate that this observation is insensitive to grain properties, which can be explained by the separation of time scales between intergrain force dynamics and intruder dynamics. Assuming there is a comparable separation of time scales, we expect that our results are applicable to a broad class of dynamic or transient granular flows. Our results suggest that descriptions of static-in-time granular flows might be extended or modified to describe these dynamic flows. Additionally, we find that accurate grain-grain interactions are not necessary to correctly capture the granular flow in this regime.

Conical wave propagation and diffraction in two-dimensional hexagonally packed granular lattices

Linear and nonlinear mechanisms for conical wave propagation in two-dimensional lattices are explored in the realm of phononic crystals. As a prototypical example, a statically compressed granular lattice of spherical particles arranged in a hexagonal packing configuration is analyzed. Upon identifying the dispersion relation of the underlying linear problem, the resulting diffraction properties are considered. Analysis both via a heuristic argument for the linear propagation of a wave packet and via asymptotic analysis leading to the derivation of a Dirac system suggests the occurrence of conical diffraction. This analysis is valid for strong precompression, i.e., near the linear regime. For weak precompression, conical wave propagation is still possible, but the resulting expanding circular wave front is of a nonoscillatory nature, resulting from the complex interplay among the discreteness, nonlinearity, and geometry of the packing. The transition between these two types of propagation is explored.

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.

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.

Topological analysis of tapped granular media using persistent homology

We use the first Betti number of a complex to analyze the morphological structure of granular samples in mechanical equilibrium. We investigate two-dimensional granular packings after a tapping process by means of both simulations and experiments. States with equal packing fraction obtained with different tapping intensities are distinguished after the introduction of a filtration parameter which determines the particles (nodes in the network) that are joined by an edge. This is accomplished by just using the position of the particles obtained experimentally and no other information about the possible contacts, or magnitude of forces.

Traveling waves in 2D hexagonal granular crystal lattices

This study describes the dynamic response of a two-dimensional hexagonal packing of uncompressed stainless steel spheres excited by localized impulsive loadings. The dynamics of the system are modeled using the Hertzian normal contact law. After the initial impact strikes the system, a characteristic wave structure emerges and continuously decays as it propagates through the lattice. Using an extension of the binary collision approximation for one-dimensional chains, we predict its decay rate, which compares well with numerical simulations and experimental data. While the hexagonal lattice does not support constant speed traveling waves, we provide scaling relations that characterize the directional power law decay of the wave velocity for various angles of impact. Lastly, we discuss the effects of weak disorder on the directional amplitude decay rates.

Molecular dynamic simulations of the elastic and inelastic surface scattering of nanoparticles

Molecular dynamics calculations have been undertaken to simulate the collision of a solid, rotating nanoparticle with a planar, two-dimensional surface at thermal velocities (linear and rotational) equivalent to 500 K. During the course of a collision, mechanisms have been introduced into the simulation process that allows for the dissipation of kinetic energy and for components of linear and angular velocity to couple. Although previous studies of particle-particle collisions have used a similar energy dissipation procedure, in these first calculations on particle-surface collisions, it is found that the mechanism actually facilitates the movement of particles across a surface. It is also shown that the direction of travel of particles on a surface is strongly influenced by their rotational motion.

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.

Quantifying stretching and rearrangement in epithelial sheet migration

Although understanding the collective migration of cells, such as that seen in epithelial sheets, is essential for understanding diseases such as metastatic cancer, this motion is not yet as well characterized as individual cell migration. Here we adapt quantitative metrics used to characterize the flow and deformation of soft matter to contrast different types of motion within a migrating sheet of cells. Using a Finite-Time Lyapunov Exponent (FTLE) analysis, we find that - in spite of large fluctuations - the flow field of an epithelial cell sheet is not chaotic. Stretching of a sheet of cells (i.e., positive FTLE) is localized at the leading edge of migration and increases when the cells are more highly stimulated. By decomposing the motion of the cells into affine and non-affine components using the metric D2min , we quantify local plastic rearrangements and describe the motion of a group of cells in a novel way. We find an increase in plastic rearrangements with increasing cell densities, whereas inanimate systems tend to exhibit less non-affine rearrangements with increasing density.

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.

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.