Impact in granular matter: Force at the base of a container made with one movable wall

In geotechnics as well as in planetary science, it is important to find a means by which to protect a base from impacts of micrometeoroids. In the moon, for example, covering a moon base with regolith, and housing such regolith by movable bounding walls, could work as a stress-leaking shield. Using a numerical model, by performing impacts on a granular material housed in a rectangular container made with one movable sidewall, it is found that such wall mobility serves as a good means for controlling the maximum force exerted at the container's base. We show that the force exerted at the container's base decreases as the movable wall decreases in mass, and it follows a Janssen-like trend. Moreover, by making use of a dynamically defined redirecting coefficient K(X), proposed by Windows-Yule et al. [Phys. Rev. E 100, 022902 (2019)2470-004510.1103/PhysRevE.100.022902], which depends on the container's width X, we propose a model for predicting the maxima measured at the container's base. The model depends on the projectile and granulate properties, and the container's geometry.

Janssen effect in dynamic particulate systems

The Janssen model of stress redistribution within laterally bounded particulate assemblies is a longstanding and valuable theoretical framework, widely used in the design of industrial systems. However, the model relies on the assumption of a static packing of particles and has never been tested in a truly dynamic regime nor for a constraining system whose geometry is dynamically altered. In this paper, we explore the pressure distributions of granular beds housed within a container possessing a laterally mobile sidewall, allowing the depth, height, and cross-sectional areas of the systems studied to be dynamically altered, thus, inducing particle rearrangements and flow in the particulate system constrained thereby. We demonstrate that the systems studied can be successfully described by the Janssen model across a wide range of system expansion rates, including those for which liquidlike flow is clearly observed and propose an extension to the model allowing for an improved characterization of constrained dynamic systems.

An instrument for studying granular media in low-gravity environment

A new experimental facility has been designed and constructed to study driven granular media in a low-gravity environment. This versatile instrument, fully automatized, with a modular design based on several interchangeable experimental cells, allows us to investigate research topics ranging from dilute to dense regimes of granular media such as granular gas, segregation, convection, sound propagation, jamming, and rheology-all without the disturbance by gravitational stresses active on Earth. Here, we present the main parameters, protocols, and performance characteristics of the instrument. The current scientific objectives are then briefly described and, as a proof of concept, some first selected results obtained in low gravity during parabolic flight campaigns are presented.

Subharmonic instability of a self-organized granular jet

Downhill flows of granular matter colliding in the lowest point of a valley, may induce a self-organized jet. By means of a quasi two-dimensional experiment where fine grained sand flows in a vertically sinusoidally agitated cylinder, we show that the emergent jet, that is, a sheet of ejecta, does not follow the frequency of agitation but reveals subharmonic response. The order of the subharmonics is a complex function of the parameters of driving.

Movers and shakers: granular damping in microgravity

The response of an oscillating granular damper to an initial perturbation is studied using experiments performed in microgravity and granular dynamics simulations. High-speed video and image processing techniques are used to extract experimental data. An inelastic hard sphere model is developed to perform simulations and the results are in excellent agreement with the experiments. In line with previous work, a linear decay of the amplitude is observed. Although this behavior is typical for a friction-damped oscillator, through simulation it is shown that this effect is still present even when friction forces are absent. A simple expression is developed which predicts the optimal damping conditions for a given amplitude and is independent of the oscillation frequency and particle inelasticities.

Coefficient of tangential restitution for viscoelastic spheres

The collision of frictional granular particles may be described by an interaction force whose normal component is that of viscoelastic spheres while the tangential part is described by the model by Cundall and Strack (Géotechnique 29, 47 (1979)) being the most popular tangential collision model in Molecular Dynamics simulations. Albeit being a rather complicated model, governed by 5 phenomenological parameters and 2 independent initial conditions, we find that it is described by 3 independent parameters only. Surprisingly, in a wide range of parameters the corresponding coefficient of tangential restitution, epsilont, is well described by the simple Coulomb law with a cut-off at epsilont = 0. A more complex behavior of the coefficient of restitution as a function on the normal and tangential components of the impact velocity, gn and gt, including negative values of epsilont, is found only for very small ratio gt/gn. For the analysis presented here we neglect dissipation of the interaction in normal direction.

Translations and rotations are correlated in granular gases

In a granular gas of rough particles the axis of rotation is shown to be correlated with the translational velocity of the particles. The average relative orientation of angular and linear velocities depends on the parameters which characterize the dissipative nature of the collision. We derive a simple theory for these correlations and validate it with numerical simulations for a wide range of coefficients of normal and tangential restitution. The limit of smooth spheres is shown to be singular: even an arbitrarily small roughness of the particles gives rise to orientational correlations.

Correction algorithm for finite sample statistics

Assume in a sample of size M one finds M(i) representatives of species i with i = 1..N*. The normalized frequency pi* triple bond Mi/M, based on the finite sample, may deviate considerably from the true probabilities p(i). We propose a method to infer rank-ordered true probabilities r(i) from measured frequencies M(i). We show that the rank-ordered probabilities provide important informations on the system, e.g., the true number of species, the Shannon- and the Renyi-entropies.

Scaling properties of granular materials

Given an assembly of viscoelastic spheres with certain material properties, we raise the question how the macroscopic properties of the assembly will change if all lengths of the system, i.e. radii, container size etc., are scaled by a constant. The result leads to a method to scale down experiments to lab size.

Extremal collision sequences of particles on a line: optimal transmission of kinetic energy

The transmission of kinetic energy through chains of inelastically colliding spheres is investigated for the case of constant coefficient of restitution epsilon=const and impact-velocity-dependent coefficient epsilon(v) for viscoelastic particles. We derive a theory for the optimal distribution of particle masses which maximize the energy transfer along the chain and check it numerically. We found that for epsilon=const, the mass distribution is a monotonous function which does not depend on the value of epsilon. In contrast, for epsilon(v) the mass distribution reveals a pronounced maximum, depending on the particle properties and on the chain length. The system investigated demonstrates that even for small and simple systems, the velocity dependence of the coefficient of restitution may lead to new effects with respect to the same systems under the simplifying approximation epsilon=const.

Coefficient of restitution of colliding viscoelastic spheres

We perform a dimension analysis for colliding viscoelastic spheres to show that the coefficient of normal restitution epsilon depends on the impact velocity g as epsilon=1-gamma(1)g(1/5)+gamma(2)g(2/5)-/+..., in accordance with recent findings. We develop a simple theory to find explicit expressions for coefficients gamma(1) and gamma(2). Using these and few next expansion coefficients for epsilon(g) we construct a Padé approximation for this function which may be used for a wide range of impact velocities where the concept of the viscoelastic collision is valid. The obtained expression reproduces quite accurately the existing experimental dependence epsilon(g) for ice particles.

The hypercube structure of the genetic code explains conservative and non-conservative aminoacid substitutions in vivo and in vitro

A representation of the genetic code as a six-dimensional Boolean hypercube is described. This structure is the result of the hierarchical order of the interaction energies of the bases in codon-anticodon recognition. In this paper it is applied to study molecular evolution in vivo and in vitro. In the first case we compared aligned positions in homologous protein sequences and found two different behaviors: (a) There are sites in which the different amino acids may be explained by one or two 'attractor nodes' (coding for the dominating amino acid(s)) and their one-bit neighbors in the codon hypercube; and (b) There are sites in which the amino acids correspond to codons located in closed paths in the hypercube. In the second case we studied the 'Sexual PCR'1 experiment described by Stemmer [Stemmer (1994)] and found that the success of this combination of usual PCR and recombination is in part due to the Gray code structure of the genetic code.