Blast shocks in quasi-two-dimensional supersonic granular flows

Organization of Embryonic Morphogenesis via Mechanical Information

Embryonic organizers establish gradients of diffusible signaling molecules to pattern the surrounding cells. Here, we elucidate an additional mechanism of embryonic organizers that is a secondary consequence of morphogen signaling. Using pharmacological and localized transgenic perturbations, 4D imaging of the zebrafish embryo, systematic analysis of cell motion, and computational modeling, we find that the vertebrate tail organizer orchestrates morphogenesis over distances beyond the range of morphogen signaling. The organizer regulates the rate and coherence of cell motion in the elongating embryo using mechanical information that is transmitted via relay between neighboring cells. This mechanism is similar to a pressure front in granular media and other jammed systems, but in the embryo the mechanical information emerges from self-propelled cell movement and not force transfer between cells. The propagation likely relies upon local biochemical signaling that affects cell contractility, cell adhesion, and/or cell polarity but is independent of transcription and translation.

Shock propagation in locally driven granular systems

We study shock propagation in a system of initially stationary hard spheres that is driven by a continuous injection of particles at the origin. The disturbance created by the injection of energy spreads radially outward through collisions between particles. Using scaling arguments, we determine the exponent characterizing the power-law growth of this disturbance in all dimensions. The scaling functions describing the various physical quantities are determined using large-scale event-driven simulations in two and three dimensions for both elastic and inelastic systems. The results are shown to describe well the data from two different experiments on granular systems that are similarly driven.

Blast Dynamics in a Dissipative Gas

The blast caused by an intense explosion has been extensively studied in conservative fluids, where the Taylor-von Neumann-Sedov hydrodynamic solution is a prototypical example of self-similarity driven by conservation laws. In dissipative media, however, energy conservation is violated, yet a distinctive self-similar solution appears. It hinges on the decoupling of random and coherent motion permitted by a broad class of dissipative mechanisms. This enforces a peculiar layered structure in the shock, for which we derive the full hydrodynamic solution, validated by a microscopic approach based on molecular dynamics simulations. We predict and evidence a succession of temporal regimes, as well as a long-time corrugation instability, also self-similar, which disrupts the blast boundary. These generic results may apply from astrophysical systems to granular gases, and invite further cross-fertilization between microscopic and hydrodynamic approaches of shock waves.

Self-organized shocks in the sedimentation of a granular gas

A granular gas in gravity heated from below develops a certain stationary density profile. When the heating is switched off, the granular gas collapses. We investigate the process of sedimentation using computational hydrodynamics, based on the Jenkins-Richman theory, and find that the process is significantly more complex than generally acknowledged. In particular, during its evolution, the system passes several stages which reveal distinct spatial regions of inertial (supersonic) and diffusive (subsonic) dynamics. During the supersonic stages, characterized by Mach>1, the system develops supersonic shocks which are followed by a steep front of the hydrodynamic fields of temperature and density, traveling upward.

Energy decay in three-dimensional freely cooling granular gas

The kinetic energy of a freely cooling granular gas decreases as a power law t(-θ) at large times t. Two theoretical conjectures exist for the exponent θ. One based on ballistic aggregation of compact spherical aggregates predicts θ=2d/(d+2) in d dimensions. The other based on Burgers equation describing anisotropic, extended clusters predicts θ=d/2 when 2≤d≤4. We do extensive simulations in three dimensions to find that while θ is as predicted by ballistic aggregation, the cluster statistics and velocity distribution differ from it. Thus, the freely cooling granular gas fits to neither the ballistic aggregation or a Burgers equation description.

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.

Unstable blast shocks in dilute granular flows

Shocks and blasts can be readily obtained in granular flows be they dense or dilute. Here, by examining the propagation of a blast shock in a dilute granular flow, we show that such a front is unstable with respect to transverse variations of the density of grains. This instability has a well-defined wavelength which depends on the density of the medium and has an amplitude which grows as an exponential of the distance traveled. These features can be understood using a simple model for the shock front, including dissipation which is inherent to granular flows. While this instability bears much resemblance to that anticipated in gases, it is distinct and has special features we discuss here.

Simulation of granular jets: is granular flow really a perfect fluid?

We perform three-dimensional simulations of the impact of a granular jet for both frictional and frictionless grains. Small shear stress observed in the experiment [X. Cheng et al., Phys. Rev. Lett. 99, 188001 (2007)] is reproduced through our simulation. However, the fluid state after the impact is far from a perfect fluid, and thus the similarity between granular jets and quark gluon plasma is superficial because the observed viscosity is finite and its value is consistent with the prediction of the kinetic theory.

Shock propagation in granular flow subjected to an external impact

We analyze a recent experiment [Boudet, Cassagne, and Kellay, Phys. Rev. Lett. 103, 224501 (2009)] in which the shock created by the impact of a steel ball on a flowing monolayer of glass beads is studied quantitatively. We argue that radial momentum is conserved in the process and hence show that in two dimensions the shock radius increases in time t as a power law t{1/3}. This is confirmed in event driven simulations of an inelastic hard sphere system. The experimental data are compared with the theoretical prediction and are shown to compare well at intermediate times. At long times the experimental data exhibit a crossover to a different scaling behavior. We attribute this to the problem becoming effectively three dimensional due to the accumulation of particles at the shock front and propose a simple hard sphere model that incorporates this effect. Simulations of this model capture the crossover seen in the experimental data.

Coarse-grained dynamics of the freely cooling granular gas in one dimension

We study the dynamics and structure of clusters in the inhomogeneous clustered regime of a freely cooling granular gas of point particles in one dimension. The coefficient of restitution is modeled as r(0)<1 or 1, depending on whether the relative speed is greater or smaller than a velocity scale δ. The effective fragmentation rate of a cluster is shown to rise sharply beyond a δ-dependent time scale. This crossover is coincident with the velocity fluctuations within a cluster becoming order δ. Beyond this crossover time, the cluster-size distribution develops a nontrivial power-law distribution, whose scaling properties are related to those of the velocity fluctuations. We argue that these underlying features are responsible for the recently observed nontrivial coarsening behavior in the one-dimensional freely cooling granular gas.