Nuestros investigadores

Diego Maza Ozcoidi

Física y Matemática Aplicada
Facultad de Ciencias. Universidad de Navarra
Líneas de investigación
Medios Granulares y Materia Blanda, Cáos y Dinámica No Lineal, Fluidos complejos
Índice H
25, (WoS, 19/10/2020)

Publicaciones científicas más recientes (desde 2010)

Autores: Lopez, D.; Cruz, Raúl; et al.
ISSN 2470-0045  Vol. 102  Nº 1  2020  págs. 010902(R)
We report experimental evidence of clogging due to the spontaneous development of hanging arches when a granular sample composed of spherical particles flows down a narrow vertical pipe. These arches, akin to the ones responsible for silo clogging, can only be possible due to the role of frictional forces; otherwise they will be unstable. We find that, contrary to the silo case, the probability of clogging in vertical narrow tubes does not decrease monotonically with the ratio of the pipe-to-particle diameters. This behavior is related to the clogging prevention caused by the spontaneous ordering of particles apparent in certain aspect ratios. More importantly, by means of numerical simulations, we discover that the interparticle normal force distributions broaden in systems with higher probability of clogging. This feature, which has been proposed before as a distinctive feature of jamming in sheared granular samples, suggests that clogging and jamming are connected in pipe flow.
Autores: Feliciani, C., (Autor de correspondencia); Zuriguel, Iker; Garcimartín, Ángel; et al.
ISSN 2045-2322  Vol. 10  Nº 1  2020  págs. 15947
Although some experimental evidence showed that an obstacle placed in front of a door allows making people's evacuations faster, the efficacy of such a solution has been debated for over 15 years. Researchers are split between those who found the obstacle beneficial and those who could not find a significant difference without it. One of the reasons for the several conclusions lies in the variety of the experiments performed so far, both in terms of competitiveness among participants, geometrical configuration and number of participants. In this work, two unique datasets relative to evacuations with/without obstacle and comprising low and high competitiveness are analyzed using state-of-the-art definitions for crowd dynamics. In particular, the so-called congestion level is employed to measure the smoothness of collective motion. Results for extreme conditions show that, on the overall, the obstacle does not reduce density and congestion level and it could rather slightly increase it. From this perspective, the obstacle was found simply shifting the dangerous spots from the area in front of the exit to the regions between the obstacle and the wall. On the other side, it was however confirmed, that the obstacle can stabilize longitudinal crowd waves, thus reducing the risk of trampling, which could be as important (in terms of safety) as improving the evacuation time.
Autores: Darias, J. R.; Fernandez, M. E.; et al.
ISSN 0032-5910  Vol. 366  2020  págs. 488 - 496
Hoppers are one of the most popular devices implemented to allow precise flow control mechanism when dispensing granulate materials from silos and other containers. Despite its ubiquity in many industrial processes, the effect that hopper geometry has on the flow rate is still poorly understood. In this work, we study the influence of the hopper angle on the main two variables that determine the flow rate: the solid-fraction and the yelorities of the particles. To this end, we use a quasi-two-dimensional system which allows a precise characterization of the profiles of these variables at the orifice. Using these experimental results, we compute the flow-density vector and obtain the resulting expression for the volumetric flow rate. Finally, we compare this expression with an equation introduced back in 1961 by RL Brown. (C) 2020 Elsevier B.V. All rights reserved.
Autores: Maza, D; Zuriguel, Iker;
ISSN 0032-5910  Vol. 360  2020  págs. 104 - 111
By means of an experimental analysis we study the granular flow in a two-dimensional silo discharged through a conveyor belt placed below the outlet. The results exhibit a saturation of the flow rate, W, with the belt velocity, v(b). Moreover, we find a dependence of the flow rate and grains velocity on the outlet size D which differs from the purely gravitational regime. To explain it, we propose an analysis based on mass conservation arguments that agrees with the experimental data when v(b) is sufficiently low. For large values of this variable, it seems to be a smooth transition between a free discharge regime (for small D) and a belt extraction regime (for large D) where the proposed model is also valid. Our analysis provides a useful connection between the flow rate and the exit geometry, a feature that may be very useful from a practical point of view. (C) 2019 Elsevier B.V. All rights reserved.
Autores: Zuriguel, Iker, (Autor de correspondencia); Maza, D; et al.
ISSN 0925-7535  Vol. 121  2019  págs. 394 - 402
We report experimental measurements obtained during the evacuation of 180 soldiers through a narrow door. Several conditions are analyzed in the evacuation drills, such as the degree of competitiveness (from rush to shove) and the influence of an obstacle placed before the exit. From the data, we compute the flow rate through the door and the velocity and density fields, as well as a map of the local evacuation time. We also present novel results on the pressure that the individuals exert on the wall adjacent to the door. Our study challenges the idea that an obstacle could be beneficial for pedestrian evacuations because of a hypothetical alleviation of pressure at the door. At the same time, we discover a correlation among the largest pressure peaks and the development of clogging.
Autores: Maza, D, (Autor de correspondencia); Schins, J.;
ISSN 0960-0779  Vol. 119  2019  págs. 237 - 242
We present experimental data of the motion of a cylindrical slider interacting only by friction with a polished horizontal tray. The tray is harmonically shacked in the horizontal direction. Below a certain threshold of the driver acceleration, the slider permanently sticks to its substrate due to the static friction. Above that threshold, the observed slider dynamics is periodic (synchronous with the driver oscillation frequency) but not wholly harmonic: for driver accelerations little beyond the threshold, the slider velocity signal is quasi-triangular. A Markovian model shows that, with increasing driver acceleration, the slider motion increasingly tends to be harmonic again, though with a prominent phase difference respect to the driver.
Autores: Zuriguel, Iker, (Autor de correspondencia); Maza, D; Janda, A.; et al.
ISSN 1434-5021  Vol. 21  Nº 3  2019  págs. 47
In this paper we report experimental and numerical results on the velocity fluctuations of grains inside silos. Although simple models exist for the stationary and continuous approximation of the flow, the variability at the microscopic level (both ensemble averages and the fluctuations of individual particles around the average) reveal non-Gaussian statistics that resist a straightforward treatment. We also show that decreasing the orifice size causes an increase in the relative amplitude of the velocity fluctuations, meaning that the intermittency grows bigger.
Autores: Gella, D.; To, K.; et al.
ISSN 2470-0045  Vol. 99  Nº 3  2019  págs. 1 - 6
We present experimental results of the effect of the hopper angle on the clogging of grains discharged from a two-dimensional silo under gravity action. We observe that the probability of clogging can be reduced by three orders of magnitude by increasing the hopper angle. In addition, we find that for very large hopper angles, the avalanche size (s) grows with the outlet size (D) stepwise, in contrast to the case of a flat-bottom silo for which s grows smoothly with D. This surprising effect is originated from the static equilibrium requirement imposed by the hopper geometry to the arch that arrests the flow. The hopper angle sets the bounds of the possible angles of the vectors connecting consecutive beads in the arch. As a consequence, only a small and specific portion of the arches that jam a flat-bottom silo can survive in hoppers.
Autores: Garcimartín, Ángel, (Autor de correspondencia); Maza, D; Pastor, J. M. ; et al.
ISSN 1367-2630  Vol. 20  2018 
The placement of obstacles in front of doors is believed to be an effective strategy to increase the flow of pedestrians, hence improving the evacuation process. Since it was first suggested, this counter-intuitive feature is considered a hallmark of pedestrian flows through bottlenecks. Indeed, despite the little experimental evidence, the placement of an obstacle has been hailed as the panacea for solving evacuation problems. In this work, we challenge this idea and experimentally demonstrate that the pedestrians flow rate is not necessarily altered by the presence of an obstacle. This result-which is at odds with recent demonstrations on its suitability for the cases of granular media, sheep and mice- differs from the outcomes of most of existing numerical models, and warns about the risks of carelessly extrapolating animal behaviour to humans. Our experimental findings also reveal an unnoticed phenomenon in relation with the crowd movement in front of the exit: in competitive evacuations, an obstacle attenuates the development of collective transversal rushes, which are hazardous as they might cause falls.
Autores: Gella, D.; Zuriguel, Iker; Maza, D;
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.
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.
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: Madrid, M.; Maza, D;
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: Gella, D.; Maza, D; Zuriguel, Iker;
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.
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)
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: Gago, P. A.; Maza, D; Pugnaloni, L. A.;
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. ;
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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: Zuriguel, Iker; Maza, D;
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.
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.
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.
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.
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.
ISSN 0031-9007  Vol. 107  Nº 27  2011  págs. 278001
Autores: Zuriguel, Iker; Maza, D;
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; Cruz, Raúl; Maza, D; et al.
ISSN 1742-5468  2010  págs. P06020
The cooling dynamics of a 2D granular gas of elongated particles is analyzed. We perform simulations on the temporal evolution of soft particles, using a molecular dynamics algorithm. For weakly dissipative particles, we found a homogeneous cooling process where the overall translational kinetic energy decreases analogously to viscoelastic circular particles. In contrast, for strongly dissipative particles we observed an inhomogeneous cooling process where the diminishing of translational kinetic energy notably slows down. The rotational kinetic energy, however, always decays in agreement with Haff's prediction for the homogeneous cooling state of inelastic particles. We mainly found that the cooling kinetics of the system is controlled by the mechanisms that determine the local energy dissipation (collisions). However, we detected a strong influence of particle shape and inelasticity on the structure of the clusters which develop in the inhomogeneous cooling regimes. Our numerical outcomes suggest that strong dissipation and particle anisotropy induce the formation of ordered cluster structures that retards the relaxation to the final asymptotic regime.
Autores: Pugnaloni, Luis A.; Sánchez, Iván; Gago, Paula A.; et al.
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: Cruz, Raúl; Zuriguel, Iker; Maza, D; et al.
ISSN 1742-5468  2010  págs. P06025
We report experimental and theoretical results of the effect that particle shape has on the packing properties of granular materials. We have systematically measured the particle angular distribution, the cluster size distribution and the stress profiles of ensembles of faceted elongated particles deposited in a bidimensional box. Stress transmission through this granular system has been numerically simulated using a two-dimensional model of irregular particles. For grains of maximum symmetry (squares), the stress propagation localizes and forms chain-like forces analogous to those observed for granular materials composed of spheres. For thick layers of grains, a pressure saturation is observed for deposit depths beyond a characteristic length. This scenario correlates with packing morphology and can be understood in terms of stochastic models of aggregation and random multiplicative processes. As grains elongate and lose their symmetry, stress propagation is strongly affected. Lateral force transmission becomes less favored than vertical transfer, and hence, an increase in the pressure develops with depth, hindering force saturation.
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.
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: 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: 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: 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.



Física (F.Ciencias). 
Universidad de Navarra - Facultad de Ciencias.

Lenguajes de programación (MC2). 
Universidad de Navarra - Facultad de Ciencias.

Física II (F.Ciencias). 
Universidad de Navarra - Facultad de Ciencias.

Adquisición de datos (MC2). 
Universidad de Navarra - Facultad de Ciencias.



Desde hace 25 años he realizado importantes esfuerzos para desarrollar y consolidar la investigación experimental en sistemas no lineales en la UNAV así como en el sistema español de CyT. Hace casi dos décadas puse en marcha el Laboratorio de Medios Granulares de la Universidad, donde se llevan a cabo a un número importantes de experimentos que además de su interés básico tienen potenciales aplicaciones industriales. Este hecho nos ha permitido interaccionar con diferentes empresas (que han financiado en parte nuestros resultados). Los aspectos internacionales y de formación siempre han sido una prioridad en mis proyectos, razón por la cual he recibido estudiantes e interactuado con unidades de I+D+i de todo el mundo, participando en redes de formación tanto a nivel de postgrado como de grado. He realizado y realizo diversas actividades de gestión tanto a nivel docente como de investigación, incluyendo además las tareas de extensión universitaria colaborando con la Universidad Nacional de Educación a Distancia donde ejerzo como Prof. Tutor y dirijo trabajos fin de grado. Recientemente he combinado mis primeras experiencias experimentales relacionadas con instabilidades en fluidos con los conocimientos adquiridos de la dinámica granular para extender el alcance de nuestro laboratorio a fluidos complejos y no newtonianos.