Autores: Gella, D.; Zuriguel, Iker; Maza, D;

Revista: PHYSICAL REVIEW LETTERS

ISSN 0031-9007
Vol. 121
Nº 13
2018
págs. 138001

Based on the implementation of a novel silo discharge procedure, we are able to control the grains velocities regardless of the outlet size. This allows isolating the geometrical and kinematic contributions to the clogging process. We find that, for a given outlet size, reducing the grains velocities to extremely low values leads to a clogging probability increment of almost two orders of magnitude, hence revealing the importance of particle kinematics in the silo clogging process. Then, we explore the contribution of both variables, outlet size and grains velocity, and we find that our results agree with an already known exponential expression that relates clogging probability with outlet size. We propose a modification of such expression revealing that only two parameters are necessary to fit all the data: one is related with the geometry of the problem, and the other with the grains kinematics.

Autores: Aumaître, S.; Behringer, R. P.; Cazaubiel, A.; et al.

Revista: REVIEW OF SCIENTIFIC INSTRUMENTS

ISSN 0034-6748
Vol. 89
Nº 075103
2018
págs. 1 - 10

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.

Autores: Asencio, K.; Acevedo, M.; Zuriguel, Iker; et al.

Revista: PHYSICAL REVIEW LETTERS

ISSN 0031-9007
Vol. 119
Nº 22
2017
págs. 228002

We experimentally analyze the compaction dynamics of an ensemble of cubic particles submitted to a novel type of excitation. Instead of the standard tapping procedure used in granular materials we apply alternative twists to the cylindrical container. Under this agitation, the development of shear forces among the different layers of cubes leads to particle alignment. As a result, the packing fraction grows monotonically with the number of twists. If the intensity of the excitations is sufficiently large, an ordered final state is reached where the volume fraction is the densest possible compatible with the boundary condition. This ordered final state resembles the tetratic or cubatic phases observed in colloids.

Autores: Gella, D.; Maza, D; Zuriguel, Iker;

Revista: PHYSICAL REVIEW E

ISSN 2470-0045
Vol. 95
Nº 5
2017
págs. 052904

We experimentally analyze the effect that particle size has on the mass flow rate of a quasi two-dimensional silo discharged by gravity. In a previous work, Janda et al. [Phys. Rev. Lett. 108, 248001 (2012)PRLTAO0031-900710.1103/PhysRevLett.108.248001] introduced a new expression for the mass flow rate based on a detailed experimental analysis of the flow for 1-mm diameter beads. Here, we aim to extend these results by using particles of larger sizes and a variable that was not explicitly included in the proposed expression. We show that the velocity and density profiles at the outlet are self-similar and scale with the outlet size with the same functionalities as in the case of 1-mm particles. Nevertheless, some discrepancies are evidenced in the values of the fitting parameters. In particular, we observe that larger particles lead to higher velocities and lower packing fractions at the orifice. Intriguingly, both magnitudes seem to compensate giving rise to very similar flow rates. In order to shed light on the origin of this behavior we have computed fields of a solid fraction, velocity, and a kinetic-stress like variable in the region above the orifice.

Autores: Gella, D.; Maza, D; Zuriguel, Iker; et al.

Revista: PHYSICAL REVIEW FLUIDS

ISSN 2469-990X
Vol. 2
Nº 8
2017
págs. 084304

We demonstrate experimentally that clogging in a silo correlates with some features of the particle velocities in the outlet proximities. This finding, that links the formation of clogs with a kinematic property of the system, is obtained by looking at the effect that the position of the lateral walls of the silo has on the flow and clogging behavior. Surprisingly, the avalanche size depends nonmonotonically on the distance of the outlet from the lateral walls. Apart from evidencing the relevance of a parameter that has been traditionally overlooked in bottleneck flow, this nonmonotonicity supposes a benchmark with which to explore the correlation of clogging probability with different variables within the system. Among these, we find that the velocity of the particles above the outlet and their fluctuations seem to be behind the nonmonotonicity in the avalanche size versus wall distance curve.

Autores: Maza, D; Cruz, Raúl, (Autor de correspondencia)

Revista: COMPUTATIONAL PARTICLE MECHANICS

ISSN 2196-4378
Vol. 4
Nº 4
2017
págs. 419 - 427

Very recently, we have examined experimentally and numerically the micro-mechanical details of monodisperse particle flows through an orifice placed at the bottom of a silo (Rubio-Largo et al. in Phys Rev Lett 114:238002, 2015). Our findings disentangled the paradoxical ideas associated to the free-fall arch concept, which has historically served to justify the dependence of the flow rate on the outlet size. In this work, we generalize those findings examining large-scale polydisperse particle flows in silos. In the range of studied apertures, both velocity and density profiles at the aperture are self-similar, and the obtained scaling functions confirm that the relevant scale of the problem is the size of the aperture. Moreover, we find that the contact stress monotonically decreases when the particles approach the exit and vanish at the outlet. The behavior of this magnitude is practically independent of the size of the orifice. However, the total and partial kinetic stress profiles suggest that the outlet size controls the propagation of the velocity fluctuations inside the silo. Examining this magnitude, we conclusively argue that indeed there is a well-defined transition region where the particle flow changes its nature. The general trend of the partial kinetic pressure profiles and the location of the transition region results the same for all particle types. We find that the partial kinetic stress is larger for bigger particles. However, the small particles carry a higher fraction of kinetic stress respect to their concentration, which suggest that the small particles have larger velocity fluctuations than the large ones and showing lower strength of correlation with the global flow. Our outcomes explain why the free-fall arch picture has served to describe the polydisperse flow rate in the discharge of silos.

Autores: Madrid, M.; Maza, D;

Revista: PHYSICAL REVIEW E

ISSN 2470-0045
Vol. 96
Nº 2
2017
págs. 022904

We present numerical and experimental results on the mass flow rate during the discharge of three-dimensional silos filled with a bidisperse mixture of grains of different sizes. We analyzed the influence of the ratio between coarse and fine particles on the profile of volume fraction and velocity across the orifice. By using numerical simulations, we have shown that the velocity profile has the same shape as that in the monodisperse case and is insensitive to the composition of the mixture. On the contrary, the volume fraction profile is strongly affected by the composition of the mixture. Assuming that an effective particle size can be introduced to characterize the mixture, we have shown that previous expression for the mass flow rate of monodisperse particles can be used for binary mixtures. A comparison with Beverloo's correlation is also presented.

Autores: Gago, P. A.; Maza, D; Pugnaloni, L. A.;

Revista: PAPERS IN PHYSICS

ISSN 1852-4249
Vol. 8
2016
págs. 080001

Static granular packs have been studied in the last three decades in the frame of a modified equilibrium statistical mechanics that assumes ergodicity as a basic postulate. The canonical example on which this framework is tested consists in the series of static configurations visited by a granular column subjected to taps. By analyzing the response of a realistic model of grains, we demonstrate that volume and stress variables visit different regions of the phase space at low tap intensities in different realizations of the experiment. We show that the tap intensity beyond which sampling by tapping becomes ergodic coincides with the forcing necessary to break all particle-particle contacts during each tap. These results imply that the well-known "reversible" branch of tapped granular columns is only valid at relatively high tap intensities.

Autores: Gago, P. A. ; Maza, D; Pugnaloni, L. A. ;

Revista: PHYSICAL REVIEW E

ISSN 1539-3755
Vol. 91
Nº 3
2015
págs. 032207

We investigate the steady-state packing fraction phi and force moment tensor Sigma of quasi-two-dimensional granular columns subjected to tapping. Systems of different height h and width L are considered. We find that phi and Sigma which describe the macroscopic state of the system, are insensitive to L for L > 50d (with d the grain diameter). However, results for granular columns of different heights cannot be conciliated. This suggests that comparison between results of different laboratories on this type of experiments can be done only for systems of same height. We show that a parameter epsilon = 1 + (A omega)(2)/(2gh), with A and omega the amplitude and frequency of the tap and g the acceleration of gravity, can be defined to characterize the tap intensity. This parameter is based on the effective flight of the granular bed, which takes into account the h dependency. When phi is plotted as a function of epsilon, the data collapses for systems of different h. However, this parameter alone is unable to determine the steady state to be reached since different Sigma can be observed for a given epsilon if different column heights are considered.

Autores: Zuriguel, Iker; Garcimartín, Ángel; et al.

Revista: GRANULAR MATTER

ISSN 1434-5021
Vol. 17
Nº 5
2015
págs. 545 - 551

We report experimental results for pipe flow of granular materials discharged through vertical narrow tubes by means of a conveyor belt placed at the bottom. When the diameter of the tube is not much larger than the particle size, the system clogs due to the development of hanging arches that are able to support the weight of the grains above them. We find that the time it takes to develop a stable clog decays exponentially, which is compatible with a clogging probability that remains constant during the discharge. From this, and making an analogy with the discharge of silos, we introduce the avalanche size, measured in terms of the number of discharged tubes before the system clogs. The mean avalanche size is found to increase as the tube diameter is enlarged, the velocity of the conveyor belt grows, and the tube tilt deviates from the vertical. © 2015, Springer-Verlag Berlin Heidelberg.

Autores: Lind, P. G. ; Maza, D; et al.

Revista: COMPUTATIONAL PARTICLE MECHANICS

ISSN 2196-4378
Vol. 2
Nº 2
2015
págs. 127 - 138

We present a hybrid GPU¿CPU implementation of an accurate discrete element model for a system of ellipsoids. The ellipsoids have three translational degrees of freedom, their rotational motion being described through quaternions and the contact interaction between two ellipsoids is described by a force which accounts for the elastic and dissipative interactions. Further we combine the exact derivation of contact points between ellipsoids (Wang et al. in Computing 72(1¿2):235¿246, 2004) with the advantages of the GPU-NVIDIA parallelization strategy (Owens et al. in Comput Graph Forum 26:80¿113, 2007). This novelty makes the analytical algorithm computationally feasible when dealing with several thousands of particles. As a benchmark, we simulate a granular gas of frictionless ellipsoids identifying a classical homogeneous cooling state for ellipsoids. For low dissipative systems, the behavior of the granular temperature indicates that the cooling dynamics is governed by the elongation of the ellipsoids and the restitution coefficient. Our outcomes comply with the statistical mechanical laws and the results are in agreement with previous findings for hard ellipsoids (Bereolos et al. in J Chem Phys 99:6087, 1993; Villemot and Talbot in Granul Matter 14:91¿97, 2012). Additionally, new insight is provided namely suggesting that the mean field description of the cooling dynamics of elongated particles is conditioned by the particle shape and the degree of energy equipartition.

Autores: Janda, A. ; Maza, D; et al.

Revista: PHYSICAL REVIEW LETTERS

ISSN 0031-9007
Vol. 114
Nº 23
2015
págs. 238002

Several theoretical predictions of the mass flow rate of granular media discharged from a silo are based
on the spontaneous development of a free-fall arch region, the existence of which is still controversial. In
this Letter, we study experimentally and numerically the particle flow through an orifice placed at the
bottom of 2D and 3D silos. The implementation of a coarse-grained technique allows a thorough
description of all the kinetic and micromechanical properties of the particle flow in the outlet proximities.
Though the free-fall arch does not exist as traditionally understood¿a region above which particles have
negligible velocity and below which particles fall solely under gravity action¿we discover that the kinetic
pressure displays a well-defined transition in a position that scales with the outlet size. This universal
scaling explains why the free-fall arch picture has served as an approximation to describe the flow rate in
the discharge of silos.

Autores: Pastor, José Martín; Garcimartín, Ángel; Gago, P. A. ; et al.

Revista: PHYSICAL REVIEW E

ISSN 2470-0045
Vol. 92
Nº 6
2015
págs. 062817

The ¿faster-is-slower¿ (FIS) effect was first predicted by computer simulations of the egress of pedestrians through a narrow exit [D. Helbing, I. J. Farkas, and T. Vicsek, Nature (London) 407, 487 (2000)]. FIS refers to the finding that, under certain conditions, an excess of the individuals' vigor in the attempt to exit causes a decrease in the flow rate. In general, this effect is identified by the appearance of a minimum when plotting the total evacuation time of a crowd as a function of the pedestrian desired velocity. Here, we experimentally show that the FIS effect indeed occurs in three different systems of discrete particles flowing through a constriction: (a) humans evacuating a room, (b) a herd of sheep entering a barn, and (c) grains flowing out a 2D hopper over a vibrated incline. This finding suggests that FIS is a universal phenomenon for active matter passing through a narrowing.

Autores: Ardanza-Trevijano, Sergio; Zuriguel, Iker; Arévalo, R.; et al.

Revista: PHYSICAL REVIEW E

ISSN 1539-3755
Vol. 89
2014
págs. 052212

We use the first Betti number of a complex to analyze the morphological structure of granular samples in mechanical equilibrium. We investigate two-dimensional granular packings after a tapping process by means of both simulations and experiments. States with equal packing fraction obtained with different tapping intensities are distinguished after the introduction of a filtration parameter which determines the particles (nodes in the network) that are joined by an edge. This is accomplished by just using the position of the particles obtained experimentally and no other information about the possible contacts, or magnitude of forces.

Autores: Acevedo, M. ; Zuriguel, Iker; Maza, D; et al.

Revista: GRANULAR MATTER

ISSN 1434-5021
Vol. 16
Nº 4
2014
págs. 411 - 420

We present experimental and numerical results for particle alignment and stress distribution in packings of faceted particles deposited in a small-scale bi-dimensional silo. First, we experimentally characterize the deposits' morphology in terms of the particles' aspect ratio and feeding rate. Then we use the experimental results to validate our discrete element method (DEM) based on spheropolygons. After achieving excellent agreement, we use contact forces and fabric provided by the simulations to calculate the coarse-grained stress tensor. For low feeding rates, square particles display a strong tendency to align downwards, i.e., with a diagonal parallel to gravity. This morphology leads to stress transmission towards the walls, implying a quick development of pressure saturation, in agreement with the Janssen effect. When the feed rate is increased, both the disorder and the number of horizontal squares in the silo increase, hindering the Janssen effect. Conversely, for elongated particles the feed rate has a weak effect on the final deposit properties. Indeed, we always observe highly ordered structures of horizontal rods where the stress is transmitted mainly in the vertical direction.

Autores: Arévalo, R.; Zuriguel, Iker; Maza, D; et al.

Revista: PHYSICAL REVIEW E

ISSN 1539-3755
Vol. 89
Nº 4
2014
págs. 042205

We present numerical results of the effect that the driving force has on the clogging probability of inert particles passing through a bottleneck. When the driving force is increased by four orders of magnitude, the mean avalanche size remains almost unaltered (increases 1.6 times) while the flow rate and the avalanche duration display strong dependence on this magnitude. This indicates that in order to characterize the ability of a system to clog, the right variable to consider is the number of particles that pass through the outlet. The weak dependence of this magnitude on the driving force is explained in terms of the average kinetic energy of the flowing grains that has to be dissipated in order to get an arch stabilized.

Autores: Zuriguel, Iker; Parisi, D. R.; Cruz, Raúl; et al.

Revista: SCIENTIFIC REPORTS

ISSN 2045-2322
Vol. 4
2014
págs. 7324

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.

Autores: Arévalo, R.; Pugnaloni, L. A.; Maza, D; et al.

Revista: PHYSICAL REVIEW E

ISSN 1539-3755
Vol. 87
Nº 2
2013
págs. 022203

We analyze the contact network of simulated two-dimensional granular packings in different states of mechanical equilibrium obtained by tapping. We show that topological descriptors of the contact network allow one to distinguish steady states of the same mean density obtained with different tap intensities. These equal-density states were recently proven to be distinguishable through the mean force moment tensor. In contrast, geometrical descriptors, such as radial distribution functions, bond order parameters, and Voronoi cell distributions, can hardly discriminate among these states. We find that small-order loops of contacts-the polygons of the network-are especially sensitive probes for the contact structure. DOI: 10.1103/PhysRevE.87.022203

Autores: Arévalo, R.; Pugnaloni, L. A.; Maza, D; et al.

Revista: PHILOSOPHICAL MAGAZINE

ISSN 1478-6435
Vol. 93
Nº 31 - 33
2013
págs. 4078 - 4089

We characterize the structure of simulated two-dimensional granular packings using concepts from complex networks theory. The packings are generated by a simulated tapping protocol, which allows us to obtain states in mechanical equilibrium in a wide range of densities. We show that our characterization method is able to discriminate non-equivalent states that have the same density. We do this by examining differences in the topological structure of the contact network of the packings. In particular, we find that the polygons of the network are specially sensitive probes for the contact structure. Additionally, we compare the network properties obtained in two different scenarios: the tapped and a compressed system.

Autores: Acevedo, Manuel Francisco; Cruz, Raúl; Zuriguel, Iker; et al.

Revista: PHYSICAL REVIEW E

ISSN 1539-3755
Vol. 87
Nº 1
2013
págs. 012202

In a previous paper [Hidalgo et al., Phys. Rev. Lett. 103, 118001 (2009)] it was shown that square particles deposited in a silo tend to align with a diagonal parallel to the gravity, giving rise to a deposit with very particular properties. Here we explore, both experimentally and numerically, the effect on these properties of the filling mechanism. In particular, we modify the volume fraction of the initial configuration from which the grains are deposited. Starting from a very dilute case, increasing the volume fraction results in an enhancement of the disorder in the final deposit characterized by a decrease of the final packing fraction and a reduction of the number of particles oriented with their diagonal in the direction of gravity. However, for very high initial volume fractions, the final packing fraction increases again. This result implies that two deposits with the same final packing fraction can be obtained from very different initial conditions. The structural properties of such deposits are analyzed, revealing that, although the final volume fraction is the same, their micromechanical properties notably differ.

Autores: Maza, D; et al.

Revista: REVISTA DE EDIFICACION

ISSN 0213-8948
Nº 41 - 42
2013
págs. 102 - 107

Los medios granulados ¿materia inerte compuesta por sólidos divididos¿ al pasar por una abertura que sea sólo un poco mayor
al tamaño de las partículas se pueden atascar. Se ha estudiado como disminuye la probabilidad de que este sistema se atasque
gracias a la colocación de un obstáculo delante de la salida. El resultado obtenido demuestra que se disminuye la probabilidad
de atasco dependiendo dónde esté colocado el obstáculo, siendo una posición óptima donde la distancia entre el obstáculo y
la salida es similar al tamaño de la abertura. Estos resultados pueden ser interesantes para el diseño óptimo de una salida de
evacuación.

Autores: Altshuler, E.; Pastor, José Martín; Garcimartín, Ángel; et al.

Revista: PLOS ONE

ISSN 1932-6203
Vol. 8
Nº 8
2013
págs. e67838

While ¿vibrational noise¿ induced by rotating components of machinery is a common problem constantly faced by engineers, the controlled conversion of translational into rotational motion or vice-versa is a desirable goal in many scenarios ranging from internal combustion engines to ultrasonic motors. In this work, we describe the underlying physics after isolating a single degree of freedom, focusing on devices that convert a vibration along the vertical axis into a rotation around this axis. A typical Vibrot (as we label these devices) consists of a rigid body with three or more cantilevered elastic legs attached to its bottom at an angle. We show that these legs are capable of transforming vibration into rotation by a ¿ratchet effect¿, which is caused by the anisotropic stick-slip-flight motion of the leg tips against the ground. Drawing an analogy with the Froude number used to classify the locomotion dynamics of legged animals, we discuss the walking regime of these robots. We are able to control the rotation frequency of the Vibrot by manipulating the shaking amplitude, frequency or waveform. Furthermore, we have been able to excite Vibrots with acoustic waves, which allows speculating about the possibility of reducing the size of the devices so they can perform tasks into the human body, excited by ultrasound waves from the outside.

Autores: Arévalo, R.; Pugnaloni, L. A.; Maza, D; et al.

Revista: PHILOSOPHICAL MAGAZINE

ISSN 1478-6435
Vol. 93
Nº 31 - 33
2013
págs. 4078 - 4089

We characterize the structure of simulated two-dimensional granular packings using concepts from complex networks theory. The packings are generated by a simulated tapping protocol, which allows us to obtain states in mechanical equilibrium in a wide range of densities. We show that our characterization method is able to discriminate non-equivalent states that have the same density. We do this by examining differences in the topological structure of the contact network of the packings. In particular, we find that the polygons of the network are specially sensitive probes for the contact structure. Additionally, we compare the network properties obtained in two different scenarios: the tapped and a compressed system.

Autores: Zuriguel, Iker; Maza, D;

Revista: PHYSICAL REVIEW LETTERS

ISSN 0031-9007
Vol. 108
Nº 24
2012
págs. 248001

"Beverloo's law" is considered as the standard expression to estimate the flow rate of particles through apertures. This relation was obtained by simple dimensional analysis and includes empirical parameters whose physical meaning is poorly justified. In this Letter, we study the density and velocity profiles in the flow of particles through an aperture. We find that, for the whole range of apertures studied, both profiles are self-similar. Hence, by means of the functionality obtained for them the mass flow rate is calculated. The comparison of this expression with the Beverloo's one reveals some differences which are crucial to understanding the mechanism that governs the flow of particles through orifices.

Autores: Garcimartín, Ángel; et al.

Revista: PHYSICAL REVIEW E

ISSN 1539-3755
Vol. 86
Nº 3
2012
págs. 031306

In a recent paper [Zuriguel et al., Phys. Rev. Lett. 107, 278001 (2011)] it has been shown that the presence of an obstacle above the outlet can significatively reduce the clogging probability of granular matter pouring from a silo. The amount of this reduction strongly depends on the obstacle position. In this work, we present new measurements to analyze different outlet sizes, extending foregoing results and revealing that the effect of the obstacle is enhanced as the outlet size is increased. In addition, the effect of the obstacle position on the flow rate properties and in the geometrical features of arches is studied. These results reinforce previous evidence of the pressure reduction induced by the obstacle. In addition, it is shown how the mean avalanche size and the average flow rate are not necessarily linked. On the other hand, a close relationship is suggested between the mean avalanche size and the flow rate fluctuations.

Autores: Pugnaloni, L. A. ; Damas, J. ; Zuriguel, Iker; et al.

Revista: PAPERS IN PHYSICS

ISSN 1852-4249
Vol. 3
Nº 0
2011
págs. 030004

We prepare static granular beds under gravity in different stationary states by tapping the system with pulsed excitations of controlled amplitude and duration. The macroscopic state---defined by the ensemble of static configurations explored by the system tap after tap---for a given tap intensity and duration is studied in terms of volume, V, and force moment tensor, \Sigma. In a previous paper [Pugnaloni et al., Phys. Rev. E 82, 050301(R) (2010)], we reported evidence supporting that such macroscopic states cannot be fully described by using only V or \Sigma, apart from the number of particles N. In this work, we present an analysis of the fluctuations of these variables that indicates that V and \Sigma may be sufficient to define the macroscopic states. Moreover, we show that only one of the invariants of \Sigma is necessary, since each component of \Sigma falls onto a master curve when plotted as a function of \rm{Tr}(\Sigma). This implies that these granular assemblies have a common shape for the stress tensor, even though it does not correspond to the hydrostatic type. Although most results are obtained by molecular dynamics simulations, we present supporting experimental results.

Autores: Garcimartín, Ángel; Zuriguel, Iker; et al.

Revista: PHYSICAL REVIEW E

ISSN 1539-3755
Vol. 84
Nº 3
2011
págs. 1 - 8

We present experimental data corresponding to a two-dimensional dense granular flow, namely, the gravity-driven discharge of grains from a small opening in a silo. We study the local velocity field at the scale of single grains at different places with the help of particle-tracking techniques. From these data, the velocity profiles can be obtained and the validity of some long-standing approaches can be assessed. Moreover, the fluctuations of the velocities are taken into consideration to characterize the features of the advective motion (due to the gravity force) and the diffusive motion, which shows nontrivial behavior.

Autores: Kanzaki, T. ; Acevedo, Manuel Francisco; Zuriguel, Iker; et al.

Revista: EUROPEAN PHYSICAL JOURNAL E

ISSN 1292-8941
Vol. 34
Nº 12
2011
págs. 133

We present experimental and numerical results of the effect that a partial discharge has on the morphological and micro-mechanical properties of non-spherical, convex particles in a silo. The comparison of the particle orientation after filling the silo and its subsequent partial discharge reveals important shear-induced orientation, which affects stress propagation. For elongated particles, the flow induces an increase in the packing disorder which leads to a reduction of the vertical stress propagation developed during the deposit generated prior to the partial discharge. For square particles, the flow favors particle alignment with the lateral walls promoting a behavior opposite to the one of the elongated particles: vertical force transmission, parallel to gravity, is induced. Hence, for elongated particles the flow developed during the partial discharge of the silo leads to force saturation with depth whereas for squares the flow induces hindering of the force saturation observed during the silo filling.

Autores: Zuriguel, Iker; Garcimartín, Ángel; et al.

Revista: PHYSICAL REVIEW LETTERS

ISSN 0031-9007
Vol. 107
Nº 27
2011
págs. 278001

Autores: Zuriguel, Iker; Maza, D;

Revista: PHYSICAL REVIEW E

ISSN 1539-3755
Vol. 81
Nº 4
2010
págs. 41302

The jamming transition of an isotropically compressed granular packing is studied by means of molecular dynamics simulations. The system is shown to undergo a critical transition which is analyzed by looking at the topological structure of the force network. At the critical packing fraction there is a sudden growth of the number of polygons in the network. Above the critical packing fraction the number of triangles keeps growing while the number of the rest of polygons is weakly reduced. Then, we prove that in the jammed regime, there is a linear relationship between the number of triangles and the coordination number. Furthermore, the presence of these minimal structures is revealed to be connected with the evolution of some important topological properties, suggesting its importance to understand the physical properties of the packing and the onset of rigidity during the compression.

Autores: Kanzaki, T; Hidalgo, RC; Maza, D; et al.

Revista: JOURNAL OF STATISTICAL MECHANICS-THEORY AND EXPERIMENT

ISSN 1742-5468
2010
págs. P06020

Autores: Pugnaloni, Luis A.; Sánchez, Iván; Gago, Paula A.; et al.

Revista: PHYSICAL REVIEW E

ISSN 1539-3755
Vol. 82
Nº 5
2010
págs. 050301

We analyze, experimentally and numerically, the steady states, obtained by tapping, of a two-dimensional granular layer. Contrary to the usual assumption, we show that the reversible (steady state branch) of the density-acceleration curve is nonmonotonous. Accordingly, steady states with the same mean volume can be reached by tapping the system with very different intensities. Simulations of dissipative frictional disks show that equal volume steady states have different values of the force moment tensor. Additionally, we find that steady states of equal stress can be obtained by changing the duration of the taps; however, these states present distinct mean volumes. These results confirm previous speculations that the volume and the force moment tensor are both needed to describe univocally equilibrium states in static granular assemblies.

Autores: Hidalgo, RC; Zuriguel, Iker; Maza, D; et al.

Revista: JOURNAL OF STATISTICAL MECHANICS-THEORY AND EXPERIMENT

ISSN 1742-5468
2010
págs. P06025

Autores: Zuriguel, Iker; Ardanza-Trevijano, Sergio; et al.

Revista: International Journal of Bifurcation and Chaos

ISSN 0218-1274
Vol. 20
Nº 3
2010
págs. 897 - 903

The existence of small order loops of contacts is presented as an intrinsic characteristic of force granular networks. Based on molecular dynamics simulations, it is proposed that the presence of these small order loops - and in particular third order loops of contacts - is important to understand the transition from fluid-like to solid-like behavior of granular packings. In addition, we show a close relationship between the development of third order loops and the small forces of the granular packing in the sense that almost all third order loops allocate a force component smaller than the average.

Autores: Sánchez, I; Gutiérrez, G; Zuriguel, Iker; et al.

Revista: PHYSICAL REVIEW E

ISSN 1539-3755
Vol. 81
Nº 6
2010
págs. 062301

We present an experimental study of the displacement of a light intruder immersed in a vibrated granular bed. Using high speed video we resolve the motion, during one cycle of oscillation, of a cylindrical object inside a Plexiglas box partially filled with grains. We report experimental evidence that, in the absence of convection, at least two forces are behind the intruder's motion: an air drag force-due to the airflow through the granular bed-and a buoyancy force produced by an air-mediated granular fluid.

Autores: Zuriguel, Iker; Janda, A.; Arévalo, R.; et al.

Libro: EPJ Web of Conferences 140, 01002 (2017). Powders & Grains 2017

Vol. 140
2017
págs. 01002

When a group of discrete particles pass through a narrowing, the flow may become arrested due to the
development of structures that span over the size of the aperture. Then, it is said that the system is clogged. Here,
we will discuss about the existence of a phase diagram for the clogged state that has been recently proposed,
arguing on its usefulness to describe different systems of discrete bodies ranging from granular materials, to
colloidal suspensions and live beings. This diagram is built based on the value of a flowing parameter which
characterizes the intermittent flow observed in all these discrete systems provided that there is an external or
internal energy supply. Such requirement, which is necessary to destabilize the clogging arches, is absent in a
standard static silo, which is therefore examined as a particular case. This view will help to understand some a
priori inconsistencies concerning the role of driving force in the clogging process that have been found in the
last years.

Autores: Asencio, K.V.; Zuriguel, Iker; Maza, D;

Libro: EPJ Web of Conferences 140, 02010 (2017). Powders & Grains 2017

Vol. 140
2017
págs. 02010

Nowadays, the common method to pack granular materials is to tap the ensemble against the gravity. Despite the apparent simplicity of that method, the asymptotic states reached by the tapped systems have strongly dependences on parameters like the shape of the tapping pulse, the container geometry or the ratio between lateral and axial dimensions. Beyond the restrictions imposed by the system boundaries, the particle shape (like rods or tetrahedrons) plays a central role in the evolution and the final state of the ensemble. In this work, we introduce an unconventional method for compacting granular ensembles by applying a sequence of alternating counterrotating pulses or ¿twists¿. By using spherical particles we analyze the efficiency of this method to achieve highly packed configurations.

Autores: López, D.; Zuriguel, Iker; Maza, D;

Libro: EPJ Web of Conferences 140, 03033 (2017). Powders & Grains 2017

Vol. 140
2017
págs. 03033

We report an experimental study on the flow of spherical particles through a vertical pipe discharged at constant velocity by means of a conveyor belt placed at the bottom. For a pipe diameter 3.67 times the diameter of the particles, we observe the development of hanging arches that stop the flow as they are able to support the weight of the particles above them. We find that the distribution of times that it takes until a stable clog develops, decays exponentially. This is compatible with a clogging probability that remains constant during the discharge. We also observe that the probability of clogging along the pipe decreases with the height, i.e. most of the clogs are developed near the bottom. This spatial dependence may be attributed to different pressure values within the pipe which might also be related to a spontaneous development of an helical structure of the grains inside the pipe.

Autores: Gella, D; Maza, D; Zuriguel, Iker;

Libro: EPJ Web of Conferences 140, 03021 (2017). Powders & Grains 2017

Vol. 140
2017
págs. 03021

Recently Janda et al. [Phys. Rev. Lett. 108, 248001 (2012)] reported an experimental study where it was measured the velocity and volume fraction fields of 1 mm diameter stainless steel beads in the exit of a two-dimensional silo. In that work, they proposed a new expression to predict the flow of granular media in silos which does not explicitly include the particle size as a parameter. Here, we study if effectively, there is not such influence of the particle size in the flux equations as well as investigate any possible effect in the velocity and volume fraction fields. To this end, we have performed high speed motion measurements of these magnitudes in a two-dimensional silo filled with 4 mm diameter beads of stainless steel, the same material than the previous works. A developed tracking program has been implemented to obtain at the same time both, the velocity and volume fraction. The final objective of this work has been to extend and generalize the theoretical framework of Janda et al. for all sizes of particles. We have found that the obtained functionalities are the same than in the 1 mm case, but the exponents and other fitting parameters are different.

Autores: Janda, A. ; Zuriguel, Iker; et al.

Libro: Particle-based methods V : fundamentals and applications

2015
págs. 71 - 80

Autores: Cruz, Raúl; Acevedo, Manuel Francisco; Zuriguel, Iker; et al.

Libro: Powders and grains 2013 : Proceedings of the 7th International Conference on Micromechanics of Granular Media

2013
págs. 895 - 898

The effect of the filling mechanism on the packing of faceted particles with different aspect ratios has been examined. We have experimentally measured the particle angular distribution and the packing fraction of ensembles of faceted particles deposited in a bidimensional box. The granular system has been numerically simulated using a two-dimensional model of faceted particles. We found that increasing the feeding rate results in an enhancement of the disorder in the final deposit and, consequently, in a reduction of the number of particles oriented in their most stable configuration. In this regime, the final packing fraction monotonically decreases as the feeding rate increases. The correlations between the final packing morphology and the stress transmission were examined by describing the micromechanical properties of the deposits. For the case of elongated particles, increasing the feeding rate leads to an enhancement of the stress transmission towards the sides of the box. On the contrary, for the case of square particles, increasing the feeding rate promotes vertical transmission of the stress.

Autores: Zuriguel, Iker; Bienzobas, J.; et al.

Libro: Powders and grains 2013 : Proceedings of the 7th International Conference on Micromechanics of Granular Media

2013
págs. 710 - 713

Transport of material through pipes or channels in mines or gravel quarries seems to be a simple and economic
form of conveying blasted ore between different levels. Despite the apparent advantages of moving the material by means of
the gravity force, there exists an important problem that makes the applicability of this method more difficult: the election of
the pipe diameter to prevent clogging of the stones. It was R. Kvapil in the sixties who extended the ideas of granular flows in
silos to underground mining. Nevertheless, after his pioneering works there are only a few manuscripts focused on this topic,
and many questions remain unsolved. In this work, we present experimental results about the flow of particles (gravel) driven
by gravity through tilted tubes. The amount of material discharged between clogs shows that the probability of clogging can
be estimated by the same procedures introduced for silos. Finally, by changing the ratio between the tube diameter and the
typical particle size, we discuss about the existence or not of a critical size beyond which clogging is not possible.

Autores: Zuriguel, Iker; et al.

Libro: Powders and grains 2013 : Proceedings of the 7th International Conference on Micromechanics of Granular Media

2013
págs. 698 - 701

We present an experimental study of the effect that an obstacle above the outlet of a silo has on the clogging
probability. Both, the size of the orice and the obstacle position are varied for a chosen obstacle size and shape. If the
position of the obstacle is properly selected the clogging probability can be importantly reduced. Indeed, as the outlet size
is increased ¿ and we approach the critical size above which there is not clogging ¿ the obstacle effect is enhanced. For the
largest outlet size studied, the clogging probability is reduced by a factor of more than one hundred. We will show, using
numerical simulations, that the physical parameter behind the reduction of the silo clogging seems to be the decrease of the
vertical pressure at the outlet proximities.

Autores: Maza, D; Rubio-Largo, S. M.; et al.

Libro: Powders and grains 2013 : Proceedings of the 7th International Conference on Micromechanics of Granular Media

2013
págs. 674 - 677

The role of density and velocity profiles in the flow of particles through apertures has been recently put on evidence
in a two-dimensional experiment (Phys. Rev. Lett. 108, 248001). For the whole range of apertures studied, both velocity and
density profiles are selfsimilar and the obtained scaling functions allow to derive the relevant scales of the problem. Indeed, by
means of the functionality obtained for these profiles, an exact expression for the mass flow rate was proposed. Such expression
showed a perfect agreement with the experiential data. In this work, we generalize this study to the three dimensional case.We
perform numerical simulations of a 3D silo in which the velocity and volume fraction profiles are determined. Both profiles
shows that the scaling obtained for 2D can be generalized to the 3D case. Finally, the scaling of the mass flow rate with the
outlet radius is discussed.