Density patterns in two-dimensional hoppers

Pushing Droplet Through a Porous Medium

I use a mechanical model of a soft body to study the dynamics of an individual fluid droplet in a random, non-wettable porous medium. The model of droplet relies on the spring–mass system with pressure. I run hundreds of independent simulations. I average droplets trajectories and calculate the averaged tortuosity of the porous domain. Results show that porous media tortuosity increases with decreasing porosity, similar to single-phase fluid study, but the form of this relationship is different.

Density distribution of particles upon jamming after an avalanche in a 2D silo

We present a complete analysis of the density distribution of particles in a two dimensional silo after discharge. Simulations through a pseudo-dynamic algorithm are performed for filling and subsequent discharge of a plane silo. Particles are monosized hard disks deposited in the container and subjected to a tapping process for compaction. Then, a hole of a given size is open at the bottom of the silo and the discharge is triggered. After a clogging at the opening is produced, and equilibrium is restored, the final distribution of the remaining particles at the silo is analyzed by dividing the space into cells with different geometrical arrangements to visualize the way in which the density depression near the opening is propagated throughout the system. The different behavior as a function of the compaction degree is discussed.Received: 9 December 2014,  Accepted: 13 April 2015; Edited by: L. A. Pugnaloni; Reviewed by: F. Vivanco, Dpto. de Física, Universidad de Santiago de Chile, Chile; DOI: http://dx.doi.org/10.4279/PIP.070007Cite as: R. O. Uñac, J. L. Sales, M. V. Gargiulo, A. M. Vidales, Papers in Physics 7, 070007 (2015)This paper, by R. O. Uñac, J. L. Sales, M. V. Gargiulo, A. M. Vidales, is licensed under the Creative Commons Attribution License 3.0. 

Role of dipole elongation in orientationally ordered liquids

A recent study shows that the dipole elongation in the extended dipole model plays a significant role in the phase transitions of liquid crystal phases. In this paper, molecular dynamics (MD) simulations were performed for the dipole model with different distances between the two charges keeping the total dipole moment the same. The potential energy consists of the Lennard-Jones potential and the site-site electrostatic contribution of partial charges. Detailed analyses were made with respect to the average order parameters (P(1)) and (P(2)) as functions of density along with other thermodynamic properties. When the reduced dipole elongations are 0.16 and 0.32, respectively, it is shown that the chainlike structures in the low density regime, liquid phases with columnar and smectic orders, and solid phases are formed; the phase with nematic order is not present anymore. At 0.64, the phases with antiferroelectric order were favored. The transition is found at the reduced elongation 0.55. It shows that the phase transitions are quite sensitive to the molecular charge distribution; this simple system could exhibit rather rich phase behaviors, which represents a significant advance in identifying molecular and state parameters of the future ferroelectric liquids.

Model for density waves in gravity-driven granular flow in narrow pipes

A gravity-driven flow of grains through a narrow pipe in vacuum is studied by means of a one-dimensional model with two coefficients of restitution. Numerical simulations show clearly how density waves form when a strikingly simple criterion is fulfilled: that dissipation due to collisions between the grains and the walls of the pipe is greater per collision than that which stems from collisions between particles. Counterintuitively, the highest flow rate is observed when the number of grains per density wave grows large. We find strong indication that the number of grains per density wave always approaches a constant as the particle number tends to infinity, and that collapse to a single wave, which was often observed also in previous simulations, occurs because the number of grains is insufficient for multiple wave formation.

Process modeling in the pharmaceutical industry using the discrete element method

The discrete element method (DEM) is widely used to model a range of processes across many industries. This paper reviews current DEM models for several common pharmaceutical processes including material transport and storage, blending, granulation, milling, compression, and film coating. The studies described in this review yielded interesting results that provided insight into the effects of various material properties and operating conditions on pharmaceutical processes. Additionally, some basic elements common to most DEM models are overviewed. A discussion of some common model extensions such as nonspherical particle shapes, noncontact forces, and interstitial fluids is also presented. While these more complex systems have been the focus of many recent studies, considerable work must still be completed to gain a better understanding of how they can affect the processing behavior of bulk solids.

Coefficient of tangential restitution for viscoelastic spheres

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

Analytical approach to continuous and intermittent bottleneck flows

We propose a many-particle-inspired theory for granular outflows from a hopper and for the escape dynamics through a bottleneck based on a continuity equation in polar coordinates. If the inflow is below the maximum outflow, we find an asymptotic stationary solution. If the inflow is above this value, we observe queue formation, which can be described by a shock wave equation. We also address the experimental observation of intermittent outflows, taking into account the lack of space in the merging zone by a minimum function and coordination problems by a stochastic variable. This results in avalanches of different sizes even if friction, force networks, inelastic collapse, or delay-induced stop-and-go waves are not assumed. Our intermittent flows result from a random alternation between particle propagation and gap propagation. Erratic flows in congested merging zones of vehicle traffic may be explained in a similar way.

Friction effects and clogging in a cellular automaton model for pedestrian dynamics

We investigate the role of conflicts in pedestrian traffic, i.e., situations where two or more people try to enter the same space. Therefore a recently introduced cellular automaton model for pedestrian dynamics is extended by a friction parameter mu. This parameter controls the probability that the movement of all particles involved in a conflict is denied at one time step. It is shown that these conflicts are not an undesirable artifact of the parallel update scheme, but are important for a correct description of the dynamics. The friction parameter mu can be interpreted as a kind of an internal local pressure between the pedestrians which becomes important in regions of high density, occurring, e.g., in panic situations. We present simulations of the evacuation of a large room with one door. It is found that friction has not only quantitative effects, but can also lead to qualitative changes, e.g., of the dependence of the evacuation time on the system parameters. We also observe similarities to the flow of granular materials, e.g., arching effects.

Simulating dynamical features of escape panic

One of the most disastrous forms of collective human behaviour is the kind of crowd stampede induced by panic, often leading to fatalities as people are crushed or trampled. Sometimes this behaviour is triggered in life-threatening situations such as fires in crowded buildings; at other times, stampedes can arise during the rush for seats or seemingly without cause. Although engineers are finding ways to alleviate the scale of such disasters, their frequency seems to be increasing with the number and size of mass events. But systematic studies of panic behaviour and quantitative theories capable of predicting such crowd dynamics are rare. Here we use a model of pedestrian behaviour to investigate the mechanisms of (and preconditions for) panic and jamming by uncoordinated motion in crowds. Our simulations suggest practical ways to prevent dangerous crowd pressures. Moreover, we find an optimal strategy for escape from a smoke-filled room, involving a mixture of individualistic behaviour and collective 'herding' instinct.

Statistics of shock waves in a two-dimensional granular flow

We have investigated the dynamics of shock waves in a single layer of uniform balls in a small-angle two-dimensional funnel. When the funnel half-angle 0 degrees < or approximately beta < or approximately 2 degrees, the flow is intermittent and kinematic shock waves are observed to propagate against the flow. We have used fast video equipment and image analysis methods to study the statistics of the shock waves. It is found that their speed and frequency increase with the distance from the outlet. In particular, the shock speed scales as the ratio of the local funnel width to the width of the funnel outlet. Various kinds of interactions between shock waves are observed, including repulsion. New shock waves are only created at those sites where a close-packed triangular packing of the monodisperse balls fits across the funnel.