Penetration depth for shallow impact cratering

We present data for the penetration of a variety of spheres, dropped from rest, into a level non-cohesive granular medium. We improve upon our earlier work [Uehara {\it et al.} Phys. Rev. Lett. {\bf 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 our 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

Numerical simulations of liquid-gas-solid three-phase flows in microgravity

Three-phase liquid-gas-solid flows under microgravity condition are studied. An Eulerian- Lagrangian computational model was developed and used in the simulations. In this approach, the liquid flow was modeled by a volume-averaged system of governing equations, whereas motions of particles and bubbles were evaluated using the Lagrangian trajectory analysis procedure. It was assumed that the bubbles remained spherical, and their shape variations were neglected. The bubble-liquid, particle-liquid and bubbl- particle two-way interactions were accounted for in the analysis. The discrete phase equations used included drag, lift, buoyancy, and virtual mass forces. Particle-particle interactions and bubble-bubble interactions were accounted for by the hard sphere model. Bubble coalescence was also included in the model. The transient flow characteristics of the three-phase flow were studied; and the effects of gravity, inlet bubble size and g-jitter acceleration on variation of flow characteristics were discussed. The low gravity simulations showed that most bubbles are aggregated in the inlet region. Also, under microgravity condition, bubble transient time is much longer than that in normal gravity. As a result, the Sauter mean bubble diameter, which is proportional to the transient time of the bubble, becomes rather large, reaching to more than 9 mm. The bubble plume in microgravity exhibits a plug type flow behavior. After the bubble plume reaches the free surface, particle volume fraction increases along the height of the column. The particles are mainly located outside the bubble plume, with very few particles being retained in the plume. In contrast to the normal gravity condition, the three phases in the column are poorly mixed under microgravity conditions. The velocities of the three phases were also found to be of the same order. Bubble size significantly affects the characteristics of the three-phase flows under microgravity conditions. For the same inlet bubble number density, the flow with larger bubbles evolves faster. The simulation results showed that the effect of g-jitter acceleration on the gas-liquid-particle three phase flows is small

Using the fractional interaction law to model the impact dynamics in arbitrary form of multiparticle collisions

Using the molecular dynamics method, we examine a discrete deterministic model for the motion of spherical particles in three-dimensional space. The model takes into account multiparticle collisions in arbitrary forms. Using fractional calculus we proposed an expression for the repulsive force, which is the so called fractional interaction law. We then illustrate and discuss how to control (correlate) the energy dissipation and the collisional time for an individual article within multiparticle collisions. In the multiparticle collisions we included the friction mechanism needed for the transition from coupled torsion-sliding friction through rolling friction to static friction. Analysing simple simulations we found that in the strong repulsive state binary collisions dominate. However, within multiparticle collisions weak repulsion is observed to be much stronger. The presented numerical results can be used to realistically model the impact dynamics of an individual particle in a group of colliding particles.Comment: 17 pages, 8 figures, 1 table; In review process of Physical Review

Compaction and dilation rate dependence of stresses in gas-fluidized beds

A particle dynamics-based hybrid model, consisting of monodisperse spherical solid particles and volume-averaged gas hydrodynamics, is used to study traveling planar waves (one-dimensional traveling waves) of voids formed in gas-fluidized beds of narrow cross sectional areas. Through ensemble-averaging in a co-traveling frame, we compute solid phase continuum variables (local volume fraction, average velocity, stress tensor, and granular temperature) across the waves, and examine the relations among them. We probe the consistency between such computationally obtained relations and constitutive models in the kinetic theory for granular materials which are widely used in the two-fluid modeling approach to fluidized beds. We demonstrate that solid phase continuum variables exhibit appreciable ``path dependence'', which is not captured by the commonly used kinetic theory-based models. We show that this path dependence is associated with the large rates of dilation and compaction that occur in the wave. We also examine the relations among solid phase continuum variables in beds of cohesive particles, which yield the same path dependence. Our results both for beds of cohesive and non-cohesive particles suggest that path-dependent constitutive models need to be developed.Comment: accepted for publication in Physics of Fluids (Burnett-order effect analysis added

Assessing a Hydrodynamic Description for Instabilities in Highly Dissipative, Freely Cooling Granular Gases

An intriguing phenomenon displayed by granular flows and predicted by kinetic-theory-based models is the instability known as particle "clustering," which refers to the tendency of dissipative grains to form transient, loose regions of relatively high concentration. In this work, we assess a modified-Sonine approximation recently proposed [Garz\'o et al., Physica A 376, 94 (2007)] for a granular gas via an examination of system stability. In particular, we determine the critical length scale associated with the onset of two types of instabilities -vortices and clusters- via stability analyses of the Navier-Stokes-order hydrodynamic equations by using the expressions of the transport coefficients obtained from both the standard and the modified-Sonine approximations. We examine the impact of both Sonine approximations over a range of solids fraction \phi <0.2 for small restitution coefficients e=0.25--0.4, where the standard and modified theories exhibit discrepancies. The theoretical predictions for the critical length scales are compared to molecular dynamics (MD) simulations, of which a small percentage were not considered due to inelastic collapse. Results show excellent quantitative agreement between MD and the modified-Sonine theory, while the standard theory loses accuracy for this highly dissipative parameter space. The modified theory also remedies a (highdissipation) qualitative mismatch between the standard theory and MD for the instability that forms more readily. Furthermore, the evolution of cluster size is briefly examined via MD, indicating that domain-size clusters may remain stable or halve in size, depending on system parameters.Comment: 4 figures; to be published in Phys. Rev.

Which effective viscosity?

Magmas undergoing shear are prime examples of flows that involve the transport of solids and gases by a separate (silicate melt) carrier phase. Such flows are called multiphase, and have attracted much attention due to their important range of engineering applications. Where the volume fraction of the dispersed phase (crystals) is large, the influence of particles on the fluid motion becomes significant and must be taken into account in any explanation of the bulk behaviour of the mixture. For congested magma deforming well in excess of the dilute limit (particle concentrations >40% by volume), sudden changes in the effective or relative viscosity can be expected. The picture is complicated further by the fact that the melt phase is temperature- and shear-rate-dependent. In the absence of a constitutive law for the flow of congested magma under an applied force, it is far from clear which of the many hundreds of empirical formulae devised to predict the rheology of suspensions as the particle fraction increases with time are best suited. Some of the more commonly used expressions in geology and engineering are reviewed with an aim to home in on those variables key to an improved understanding of magma rheology. These include a temperature, compositional and shear-rate dependency of viscosity of the melt phase with the shear-rate dependency of the crystal (particle) packing arrangement. Building on previous formulations, a new expression for the effective (relative) viscosity of magma is proposed that gives users the option to define a packing fraction range as a function of shear stress. Comparison is drawn between processes (segregation, clustering, jamming), common in industrial slurries, and structures seen preserved in igneous rocks. An equivalence is made such that congested magma, viewed in purely mechanical terms as a high-temperature slurry, is an inherently non-equilibrium material where flow at large Péclet numbers may result in shear thinning and spontaneous development of layering

Hydrodynamic Description of Granular Convection

We present a hydrodynamic model that captures the essence of granular dynamics in a vibrating bed. We carry out the linear stability analysis and uncover the instability mechanism that leads to the appearance of the convective rolls via a supercritical bifurcation of a bouncing solution. We also explicitly determine the onset of convection as a function of control parameters and confirm our picture by numerical simulations of the continuum equations.Comment: 14 pages, RevTex 11pages + 3 pages figures (Type csh

Self-Structuring of Granular material under Capillary Bulldozing

An experimental observation of the structuring of a granular suspension under the progress of a gas/liquid meniscus in a narrow tube is reported here. The granular material is moved and compactifies as a growing accumulation front. The frictional interaction with the confining walls increases until the pore capillary entry pressure is reached. The gas then penetrates the clogged granular packing and a further accumulation front is formed at the far side of the plug. This cyclic process continues until the gas/liquid interface reaches the tube’s outlet, leaving a trail of plugs in the tube. Such 1D pattern formation belongs to a larger family of patterning dynamics observed in 2D Hele-Shaw geometry. The cylindrical geometry considered here provides an ideal case for a theoretical modelling for forced granular matter oscillating between a long frictional phase and a sudden viscous fluidization

Gas turbulence modulation in a two-fluid model for gas-solid flows

Recent rapid progress in the theoretical and experimental study of turbulence modulation has led to greater understanding of the physics of particle-gas turbulence interactions. A new two-fluid model incorporating these advances for relatively dilute gas-solid flows containing high-inertia particles is established. The effect of aerodynamic forces upon the particulate stresses is considered in this kinetic theory-based model, and the influence of the particles on the turbulent gas is addressed: the work associated with drag forces contributes to the gas turbulent energy, and the space occupied by particles restricts the turbulent length scale. The interparticle length scale, which is usually ignored, has been incorporated into a new model for determining the turbulent length scale. This model also considers the transport effect on the turbulent length scale. Simulation results for fully developed steady flows in vertical pipes are compared with a wide range of published experimental data and, generally, good agreement is shown. This comprehensive and validated model accounts for many of the interphase interactions that have been shown to be important