Nuestros investigadores

Diego Maza Ozcoidi

Departamento
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
24, (WoS, 20/06/2018)

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

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: 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.

ACTIVIDAD DOCENTE

   

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

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