Spontaneous emergence of counterclockwise vortex motion in assemblies of pedestrians roaming within an enclosure.

The emergence of coherent vortices has been observed in a wide variety of many-body systems such as animal flocks, bacteria, colloids, vibrated granular materials or human crowds. Here, we experimentally demonstrate that pedestrians roaming within an enclosure also form vortex-like patterns which, intriguingly, only rotate counterclockwise. By implementing simple numerical simulations, we evidence that the development of swirls in many-particle systems can be described as a phase transition in which both the density of agents and their dissipative interactions with the boundaries play a determinant role. Also, for the specific case of pedestrians, we show that the preference of right-handed people (the majority in our experiments) to turn leftwards when facing a wall is the symmetry breaking mechanism needed to trigger the global counterclockwise rotation observed.

Visual analysis of density and velocity profiles in dense 3D granular gases.

Granular multiparticle ensembles are of interest from fundamental statistical viewpoints as well as for the understanding of collective processes in industry and in nature. Extraction of physical data from optical observations of three-dimensional (3D) granular ensembles poses considerable problems. Particle-based tracking is possible only at low volume fractions, not in clusters. We apply shadow-based and feature-tracking methods to analyze the dynamics of granular gases in a container with vibrating side walls under microgravity. In order to validate the reliability of these optical analysis methods, we perform numerical simulations of ensembles similar to the experiment. The simulation output is graphically rendered to mimic the experimentally obtained images. We validate the output of the optical analysis methods on the basis of this ground truth information. This approach provides insight in two interconnected problems: the confirmation of the accuracy of the simulations and the test of the applicability of the visual analysis. The proposed approach can be used for further investigations of dynamical properties of such media, including the granular Leidenfrost effect, granular cooling, and gas-clustering transitions.

The role of initial speed in projectile impacts into light granular media.

Projectile impact into a light granular material composed of expanded polypropylene (EPP) particles is investigated systematically with various impact velocities. Experimentally, the trajectory of an intruder moving inside the granular material is monitored with a recently developed non-invasive microwave radar system. Numerically, discrete element simulations together with coarse-graining techniques are employed to address both dynamics of the intruder and response of the granular bed. Our experimental and numerical results of the intruder dynamics agree with each other quantitatively and are in congruent with existing phenomenological model on granular drag. Stepping further, we explore the 'microscopic' origin of granular drag through characterizing the response of granular bed, including density, velocity and kinetic stress fields at the mean-field level. In addition, we find that the dynamics of cavity collapse behind the intruder changes significantly when increasing the initial speed . Moreover, the kinetic pressure ahead of the intruder decays exponentially in the co-moving system of the intruder. Its scaling gives rise to a characteristic length scale, which is in the order of intruder size. This finding is in perfect agreement with the long-scale inertial dissipation type that we find in all cases.

Clogging transition of many-particle systems flowing through bottlenecks.

When a large set of discrete bodies passes through a bottleneck, the flow may become intermittent due to the development of clogs that obstruct the constriction. Clogging is observed, for instance, in colloidal suspensions, granular materials and crowd swarming, where consequences may be dramatic. Despite its ubiquity, a general framework embracing research in such a wide variety of scenarios is still lacking. We show that in systems of very different nature and scale -including sheep herds, pedestrian crowds, assemblies of grains, and colloids- the probability distribution of time lapses between the passages of consecutive bodies exhibits a power-law tail with an exponent that depends on the system condition. Consequently, we identify the transition to clogging in terms of the divergence of the average time lapse. Such a unified description allows us to put forward a qualitative clogging state diagram whose most conspicuous feature is the presence of a length scale qualitatively related to the presence of a finite size orifice. This approach helps to understand paradoxical phenomena, such as the faster-is-slower effect predicted for pedestrians evacuating a room and might become a starting point for researchers working in a wide variety of situations where clogging represents a hindrance.

Brittle-to-ductile transition in a fiber bundle with strong heterogeneity.

We analyze the failure process of a two-component system with widely different fracture strength in the framework of a fiber bundle model with localized load sharing. A fraction 0≤α≤1 of the bundle is strong and it is represented by unbreakable fibers, while fibers of the weak component have randomly distributed failure strength. Computer simulations revealed that there exists a critical composition α(c) which separates two qualitatively different behaviors: Below the critical point, the failure of the bundle is brittle, characterized by an abrupt damage growth within the breakable part of the system. Above α(c), however, the macroscopic response becomes ductile, providing stability during the entire breaking process. The transition occurs at an astonishingly low fraction of strong fibers which can have importance for applications. We show that in the ductile phase, the size distribution of breaking bursts has a power law functional form with an exponent µ=2 followed by an exponential cutoff. In the brittle phase, the power law also prevails but with a higher exponent µ=9/2. The transition between the two phases shows analogies to continuous phase transitions. Analyzing the microstructure of the damage, it was found that at the beginning of the fracture process cracks nucleate randomly, while later on growth and coalescence of cracks dominate, which give rise to power law distributed crack sizes.

Role of particle shape on the stress propagation in granular packings.

We present an experimental and numerical study on the influence that particle aspect ratio has on the mechanical and structural properties of granular packings. For grains with maximal symmetry (squares), the stress propagation in the packing localizes forming chainlike forces analogous to the ones observed for spherical grains. This scenario can be understood in terms of stochastic models of aggregation and random multiplicative processes. As the grains elongate, the stress propagation is strongly affected. The interparticle normal force distribution tends toward a Gaussian, and, correspondingly, the force chains spread leading to a more uniform stress distribution reminiscent of the hydrostatic profiles known for standard liquids.

Phase transitions of semiflexible hard-sphere chain liquids.

We present a density-functional theory for describing liquid-crystalline phase transitions in a fluid of semiflexible hard-sphere chain molecules based on the Onsager second-virial approximation for the free energy. Key ingredients of this theory are the generation of semiflexible chain conformations and calculation of the pair excluded volume and excluded area using a single-chain Monte Carlo enumeration method. First, we investigate the isotropic-nematic phase transition. Next, the theory is extended to account for a smectic-A phase by a bifurcation analysis around the nematic solution. The perturbation is calculated using a lowest-order Fourier series representation.

Smectic ordering of homogeneous semiflexible polymers.

A self-consistent-field theory for fluids of homogeneous wormlike polymers exhibiting a one-dimensional spatial variation is presented. We have extended the treatment of excluded-volume effects by adding an effective interaction term which describes the excluded volume between wormlike cylindrical segments and terminal (or end) segments of the polymer molecules. This enables us to find a smectic-A phase in the case of homogeneous semiflexible polymers. Using this framework, we have investigated the occurrence of smectic-A, nematic, and isotropic phases in the second-virial (Onsager) approximation. Phase diagrams are calculated for systems characterized by different rigidities (i.e., persistence lengths). For the case of infinitely rigid molecules, the nematic-smectic transition appears to be mostly second order. Systems of semiflexible molecules exhibit mainly a first-order smectic-nematic transition, and their isotropic-nematic-smectic triple points are accessed for different rigidity values. The nematic-smectic transition line is in good agreement with previous analytical calculations, which were also performed assuming the second-virial approximation. However, the values of the volume fraction at the nematic-smectic transition are large compared with computer simulation results, indicating limitations of the second-virial approximation.

Time evolution of damage under variable ranges of load transfer.

We study the time evolution of damage in a fiber bundle model in which the range of interaction of fibers varies through an adjustable stress transfer function recently introduced. We find that the lifetime of the material exhibits a crossover from mean-field to short-range behavior as in the static case. Numerical calculations showed that the value at which the transition takes place depends on the system's disorder. Finally, we have performed a microscopic analysis of the failure process. Our results confirm that the growth dynamics of the largest crack is radically different in the two limiting regimes of load transfer during the first stages of breaking.

Scaling laws of creep rupture of fiber bundles.

We study the creep rupture of fiber composites in the framework of fiber bundle models. Two fiber bundle models are introduced based on different microscopic mechanisms responsible for the macroscopic creep behavior. Analytical and numerical calculations show that above a critical load the deformation of the creeping system monotonically increases in time resulting in global failure at a finite time t(f), while below the critical load the system suffers only partial failure and the deformation tends to a constant value giving rise to an infinite lifetime. It is found that approaching the critical load from below and above the creeping system is characterized by universal power laws when the fibers have long-range interaction. The lifetime of the composite above the critical point has a universal dependence on the system size.

Evolution of percolating force chains in compressed granular media.

The evolution of effective force chains percolating through a compressed granular system is investigated. We compressed an ensemble of spherical particles monitoring the macroscopic constitutive behavior and the acoustic signals emitted by microscopic rearrangements of particles. We applied the continuous damage model of fiber bundles to describe the evolution of the array of force chains. The model provides a nonlinear constitutive behavior in good quantitative agreement with the experimental results. The amplitude distribution of acoustic signals was found experimentally to follow a power law of an exponent delta=1.15+/-0.05, which is in good agreement with the analytic solution of the model.

Fracture model with variable range of interaction.

We introduce a fiber bundle model where the interaction among fibers is modeled by an adjustable stress-transfer function that can interpolate between the two limiting cases of load redistribution, i.e., the global and the local load sharing schemes. By varying the range of interaction, several features of the model are numerically studied and a crossover from mean-field to short-range behavior is obtained. The properties of the two regimes and the emergence of the crossover in between are explored by numerically studying the dependence of the ultimate strength of the material on the system size, the distribution of avalanches of breakings, and of the cluster sizes of broken fibers. Finally, we analyze the moments of the cluster size distributions to accurately determine the value at which the crossover is observed.

Creep rupture of viscoelastic fiber bundles.

We study the creep rupture of bundles of viscoelastic fibers occurring under uniaxial constant tensile loading. A fiber bundle model is introduced that combines the viscoelastic constitutive behavior and the strain controlled breaking of fibers. Analytical and numerical calculations showed that above a critical external load the deformation of the system monotonically increases in time resulting in global failure at a finite time t(f), while below the critical load the deformation tends to a constant value giving rise to an infinite lifetime. Our studies revealed that the nature of the transition between the two regimes, i.e., the behavior of t(f) at the critical load sigma(c), strongly depends on the range of load sharing: for global load sharing t(f) has a power law divergence at sigma(c) with a universal exponent of 0.5, however, for local load sharing the transition becomes abrupt: at the critical load t(f) jumps to a finite value, analogous to second- and first-order phase transitions, respectively. The acoustic response of the bundle during creep is also studied.