No Arabic abstract
Granular intrusions, such as dynamic impact or wheel locomotion, are complex multiphase phenomena where the grains exhibit solid-like and fluid-like characteristics together with an ejected gas-like phase. Despite decades of modeling efforts, a unified description of the physics in such intrusions is as yet unknown. Here we show that a continuum model based on the simple notions of frictional flow and tension-free separation describes complex granular intrusions near free surfaces. This model captures dynamics in a variety of experiments including wheel locomotion, plate intrusions, and running legged robots. The model reveals that three effects (a static contribution and two dynamic ones) primarily give rise to intrusion forces in such scenarios. Identification of these effects enables the development of a further reduced-order technique (Dynamic Resistive Force Theory) for rapid modeling of granular locomotion of arbitrarily shaped intruders. The continuum-motivated strategy we propose for identifying physical mechanisms and corresponding reduced-order relations has potential use for a variety of other materials.
The drag force exerted on an object intruding into granular media can depend on the objects velocity as well as the depth penetrated. We report on intrusion experiments at constant speed over four orders in magnitude together with systematic molecular dynamics simulations well beyond the quasi-static regime. We find that velocity dependence crosses over to depth dependence at a characteristic time after initial impact. This crossover time scale, which depends on penetration speed, depth, gravity and intruder geometry, challenges current models that assume additive contributions to the drag.
We analyze the capabilities of various recently developed techniques, namely Resistive Force Theory (RFT) and continuum plasticity implemented with the Material Point Method (MPM), in capturing dynamics of wheel--dry granular media interactions. We compare results to more conventionally accepted methods of modeling wheel locomotion. While RFT is an empirical force model for arbitrarily-shaped bodies moving through granular media, MPM-based continuum modeling allows the simulation of full granular flow and stress fields. RFT allows for rapid evaluation of interaction forces on arbitrary shaped intruders based on a local surface stress formulation depending on depth, orientation, and movement of surface elements. We perform forced-slip experiments for three different wheel types and three different granular materials, and results are compared with RFT, continuum modeling, and a traditional terramechanics semi-empirical method. Results show that for the range of inputs considered, RFT can be reliably used to predict rigid wheel granular media interactions with accuracy exceeding that of traditional terramechanics methodology in several circumstances. Results also indicate that plasticity-based continuum modeling provides an accurate tool for wheel-soil interaction while providing more information to study the physical processes giving rise to resistive stresses in granular media.
Granular flow in a silo demonstrates multiple nonlocal rheological phenomena due to the finite size of grains. We solve the Nonlocal Granular Fluidity (NGF) continuum model in quasi-2D silo geometries and evaluate its ability to predict these nonlocal effects, including flow spreading and, importantly, clogging (arrest) when the opening is small enough. The model is augmented to include a free-separation criterion and is implemented numerically with an extension of the trans-phase granular flow solver described in arXiv:1411.5447, to produce full-field solutions. The implementation is validated against analytical results of the model in the inclined chute geometry, such as the solution for the $H_{mathrm{stop}}$ curve for size-dependent flow arrest, and the velocity profile as a function of layer height. We then implement the model in the silo geometry and vary the apparent grain size. The model predicts a jamming criterion when the opening competes with the scale of the mean grain size, which agrees with previous experimental studies, marking the first time to our knowledge that silo jamming has been achieved with a continuum model. For larger openings, the flow within the silo obtains a diffusive characteristic whose spread depends on the models nonlocal amplitude and the mean grain size. The numerical tests are controlled for grid effects and a comparison study of coarse vs refined numerical simulations shows agreement in the pressure field, the shape of the arch in a jammed silo configuration, and the velocity field in a flowing configuration.
Using high-speed video and magnetic resonance imaging (MRI) we study the motion of a large sphere in a vertically vibrated bed of smaller grains. As previously reported we find a non-monotonic density dependence of the rise and sink time of the large sphere. We find that this density dependence is solely due to air drag. We investigate in detail how the motion of the intruder sphere is influenced by size of the background particles, initial vertical position in the bed, ambient pressure and convection. We explain our results in the framework of a simple model and find quantitative agreement in key aspects with numerical simulations to the model equations.
We have extended the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) to support directional bonds and dynamic bonding. The framework supports stochastic formation of new bonds, breakage of existing bonds, and conversion between bond types. Bond formation can be controlled to limit the maximal functionality of a bead with respect to various bond types. Concomitant with the bond dynamics, angular and dihedral interactions are dynamically introduced between newly connected triplets and quartets of beads, where the interaction type is determined from the local pattern of bead and bond types. When breaking bonds, all angular and dihedral interactions involving broken bonds are removed. The framework allows chemical reactions to be modeled, and use it to simulate a simplistic, coarse-grained DNA model. The resulting DNA dynamics illustrate the power of the present framework.