Shear-induced alignment and dynamics of elongated granular particles

Exploring noisy Jeffery orbits: A combined Fokker-Planck and Langevin analysis in two and three dimensions

The behavior of nonspherical particles in a shear flow is of significant practical and theoretical interest. These systems have been the object of numerous investigations since the pioneering work of Jeffery a century ago. His eponymous orbits describe the deterministic motion of an isolated, rodlike particle in a shear flow. Subsequently, the effect of adding noise was investigated. The theory has been applied to colloidal particles, macromolecules, anisometric granular particles, and most recently to microswimmers, for example, bacteria. We study the Jeffery orbits of elongated (uniaxial, prolate) particles subject to noise using Langevin simulations and a Fokker-Planck equation. We extend the analytical solution for infinitely thin needles (β=1) obtained by Doi and Edwards to particles with arbitrary shape factor (0≤β≤1) and validate the theory by comparing it with simulations. We examine the rotation of the particle around the vorticity axis and study the orientational order matrix. We use the latter to obtain scalar order parameters s and r describing nematic ordering and biaxiality from the orientational distribution function. The value of s (nematic ordering) increases monotonically with increasing Péclet number, while r (measure of biaxiality) displays a maximum value. From perturbation theory, we obtain simple expressions that provide accurate descriptions at low noise (or large Péclet numbers). We also examine the orientational distribution in the v-grad v plane and in the perpendicular direction. Finally, we present the solution of the Fokker-Planck equation for a strictly two-dimensional (2D) system. For the same noise amplitude, the average rotation speed of the particle in 3D is larger than in 2D.

Effect of Particle Shape on the Flow of an Hourglass

The flow rate of a granulate out of a cylindrical container is studied as a function of particle shape for flat and elongated ellipsoids experimentally and numerically. We find a nonmonotonic dependence of the flow rate on the grain aspect ratio a/b. Starting from spheres the flow rate grows and has two maxima around the aspect ratios of a/b≈0.6 (lentil-like ellipsoids) and a/b≈1.5 (ricelike ellipsoids) reaching a flow rate increase of about 15% for lentils compared to spheres. For even more anisometric shapes (a/b=0.25 and a/b=4) the flow rate drops. Our results reveal two contributing factors to the nonmonotonic nature of the flow rate: both the packing fraction and the particle velocity through the orifice are nonmonotonic functions of the grain shape. Thus, particles with slightly nonspherical shapes not only form a better packing in the silo but also move faster through the orifice than spheres. We also show that the resistance of the granulate against shearing increases with aspect ratio for both elongated and flat particles; thus change in the effective friction of the granulate due to changing particle shape does not coincide with the trend in the flow rate.

Shear flow of non-Brownian rod-sphere mixtures near jamming

We use the discrete element method, taking particle contact and hydrodynamic lubrication into account, to unveil the shear rheology of suspensions of frictionless non-Brownian rods in the dense packing fraction regime. We find that, analogously to the random close packing volume fraction, the shear-driven jamming point of this system varies in a nonmonotonic fashion as a function of the rod aspect ratio. The latter strongly influences how the addition of rodlike particles affects the rheological response of a suspension of frictionless non-Brownian spheres to an external shear flow. At fixed values of the total (rods plus spheres) packing fraction, the viscosity of the suspension is reduced by the addition of "short"(≤2) rods but is instead increased by the addition of "long"(≥2) rods. A mechanistic interpretation is provided in terms of packing and excluded-volume arguments.

Elongated particles discharged with a conveyor belt in a two-dimensional silo

The flow of elliptical particles out of a two-dimensional silo when extracted with a conveyor belt is analyzed experimentally. The conveyor belt-placed directly below the silo outlet-reduces the flow rate, increases the size of the stagnant zone, and it has a very strong influence on the relative velocity fluctuations as they strongly increase everywhere in the silo with decreasing belt speed. In other words, instead of slower but smooth flow, flow reduction by belt leads to intermittent flow. Interestingly, we show that this intermittency correlates with a strong reduction of the orientational order of the particles at the orifice region. Moreover, we observe that the average orientation of the grains passing through the outlet is modified when they are extracted with the belt, a feature that becomes more evident for large orifices.

Discharge of elongated grains in silos under rotational shear

The discharge of elongated particles from a silo with rotating bottom is investigated numerically. The introduction of a slight transverse shear reduces the flow rate Q by up to 70% compared with stationary bottom, but the flow rate shows a modest increase by further increasing the external shear. Focusing on the dependency of flow rate Q on orifice diameter D, the spheres and rods show two distinct trends. For rods, in the small-aperture limit Q seems to follow an exponential trend, deviating from the classical power-law dependence. These macroscopic observations are in good agreement with our earlier experimental findings [Phys. Rev. E 103, 062905 (2021)2470-004510.1103/PhysRevE.103.062905]. With the help of the coarse-graining methodology we obtain the spatial distribution of the macroscopic density, velocity, kinetic pressure, and orientation fields. This allows us detecting a transition from funnel to mass flow pattern caused by the external shear. Additionally, averaging these fields in the region of the orifice reveals that the strong initial decrease in Q is mostly attributed to changes in the flow velocity, while the weakly increasing trend at higher rotation rates is related to increasing packing fraction. Similar analysis of the grain orientation at the orifice suggests a correlation of the flow rate magnitude with the vertical orientation and the packing fraction at the orifice with the order of the grains. Lastly, the vertical profile of mean acceleration at the center of the silo denotes that the region where the acceleration is not negligible shrinks significantly due to the strong perturbation induced by the moving wall.

Dense shearing flows of soft, frictional cylinders

We perform discrete numerical simulations at a constant volume of dense, steady, homogeneous flows of true cylinders interacting via Hertzian contacts, with and without friction, in the absence of preferential alignment. We determine the critical values of the solid volume fraction and the average number of contacts per particle above which rate-independent components of the stresses develop, along with a sharp increase in the fluctuations of angular velocity. We show that kinetic theory, extended to account for a velocity correlation at solid volume fractions larger than 0.49, can quantitatively predict the measured fluctuations of translational velocity, at least for sufficiently rigid cylinders, for any value of the cylinder aspect ratio and friction investigated here. The measured pressure above and below the critical solid volume fraction is in agreement with a recent theory originally intended for spheres that conjugates extended kinetic theory, the finite duration of collisions between soft particles and the development of an elastic network of long-lasting contacts responsible for the rate-independency of the flows in the supercritical regime. Finally, we find that, for sufficiently rigid cylinders, the ratio of shear stress to pressure in the subcritical regime is a linear function of the ratio of the shear rate to a suitable measure of the fluctuations of translational velocity, in qualitative accordance with kinetic theory, with an intercept that increases with friction. A decrease in the particle stiffness gives rise to nonlinear effects that greatly diminishes the stress ratio.

Directional shear jamming of frictionless ellipses

In this work we study shear reversals of dense non-Brownian suspensions composed of cohesionless elliptical particles. By numerical simulations, we show that a new fragility appears for frictionless ellipses in the flowing states, where particles can flow indefinitely in one direction at applied shear stresses but shear jam in the other direction upon shear stress reversal. This new fragility, absent in the isotropic particle case, is linked to the directional order of the elongated particles at steady shear and its reorientation at shear stress reversal, which forces the suspensions to pass through a more disordered state with an increased number of contacts in which it might get arrested.

Shear-induced orientational ordering in an active glass former

Dense assemblies of self-propelled particles that can form solid-like states also known as active or living glasses are abundant around us, covering a broad range of length scales and timescales: from the cytoplasm to tissues, from bacterial biofilms to vehicular traffic jams, and from Janus colloids to animal herds. Being structurally disordered as well as strongly out of equilibrium, these systems show fascinating dynamical and mechanical properties. Using extensive molecular dynamics simulation and a number of distinct dynamical and mechanical order parameters, we differentiate three dynamical steady states in a sheared model active glassy system: 1) a disordered state, 2) a propulsion-induced ordered state, and 3) a shear-induced ordered state. We supplement these observations with an analytical theory based on an effective single-particle Fokker-Planck description to rationalize the existence of the shear-induced orientational ordering behavior in an active glassy system without explicit aligning interactions of, for example, Vicsek type. This ordering phenomenon occurs in the large persistence time limit and is made possible only by the applied steady shear. Using a Fokker-Planck description with parameters that can be measured independently, we make testable predictions for the joint distribution of single-particle position and orientation. These predictions match well with the joint distribution measured from direct numerical simulation. Our results are of relevance for experiments exploring the rheological response of dense active colloids and jammed active granular matter systems.

Numerical investigation on effect of particle aspect ratio on the dynamical behaviour of ellipsoidal particle flow

Flow of ellipsoidal particles in a modal shear cell was investigated at the microdynamic level based on discrete element method simulations. In a stress-controlled double-shear condition, the flow was studied by varying the aspect ratio of ellipsoidal particles and comparing with the flow of spherical particle assembly in terms of some key properties, including particle alignment, linear velocity, angular velocity, porosity, contact force and contact energy. It was found that particle elongation impacts the rotational displacement around the axis perpendicular to the shear direction, which causes that the ellipsoidal particles with higher elongation are more aligned with the direction of the shear velocity, with more uniform force network. This then affects other particle properties. The fluctuation of linear velocity and the angular velocity decreases with an increase in particle aspect ratio, although the particle elongation does not significantly affect the flow velocity gradient. There is a reduction in both normal and tangential forces per contact with an increase of particle elongation. Due to the variation of the particle alignment with elongation, the standard deviation of the contact energies increases and then reduces when an increase in particle aspect ratio occurs, and on contrary, the porosity has an opposite variation trend.

Shape or friction? Which of these characteristics drives the shear strength in granular systems?

The shape of the particles and local friction, separately, are known to strongly affect the macroscopic properties of an assembly of grains. But the combined effects of these two parameters still remain poorly described. By means of extensive two dimensional contact dynamics simulations, we perform a systematic analysis of the interplay between friction and shape on strength properties of granular systems. The shape of the particles is varied from disks to triangles, while the friction is varied from 0 to 0.7. We find that the macroscopic friction first increases with angularity, but it may decline (for low friction values), saturate (for intermediates friction values), or continue to increase (for large friction values) for the most angular shapes. In other words, the effect of the particle’s angularity on the shear strength depends on the level of sliding friction. In contrast, the effect of local friction on the shear strength does not depend on the specific properties of shape. The results presented here highlight the subtle coupling existing between shape and friction effects.

Translational and rotational velocities in shear-driven jamming of ellipsoidal particles

We study shear-driven jamming of ellipsoidal particles at zero temperature with a focus on the microscopic dynamics. We find that a change from spherical particles to ellipsoids with aspect ratio α=1.02 gives dramatic changes of the microscopic dynamics with much lower translational velocities and a new role for the rotations. Whereas the velocity difference at contacts-and thereby the dissipation-in collections of spheres is dominated by the translational velocities and reduced by the rotations, the same quantity is in collections of ellipsoids instead totally dominated by the rotational velocities. By also examining the effect of different aspect ratios we find that the examined quantities show either a peak or a change in slope at α≈1.2, which thus gives evidence for a crossover between different regions of low and high aspect ratio.

Orientational order in dense suspensions of elliptical particles in the non-Stokesian regime

Suspensions of neutrally buoyant elliptic particles are modeled in 2D using fully resolved simulations that provide two-way interaction between the particle and the fluid medium. Forces due to particle collisions are represented by a diffuse interface approach that allows the investigation of dense suspensions (up to 47% packing fraction). We focus on the role inertial forces play at low and high particle Reynolds numbers termed low Reynolds number and inertial regimes, respectively. The suspensions are characterized by the orientation distribution function (ODF) that reflects shear induced rotation of the particles at low Reynolds numbers, and nearly stationary (swaying) particles at high Reynolds numbers. In both cases, orientational ordering differs qualitatively from the behavior observed in the Stokesian-regime. The ODF becomes flatter with increasing packing fraction, as opposed to the sharpening previous work predicted in the Stokesian regime. The ODF at low particle concentrations differs significantly for the low Reynolds number and inertial regimes, whereas with increasing packing fraction convergence is observed. For dense suspensions, the particle-particle interactions dominate the particle motion.

Depletion forces in athermally sheared mixtures of frictionless disks and rods in two dimensions

We carry out numerical simulations to study the behavior of an athermal mixture of frictionless circular disks and elongated rods in two dimensions, under three different types of global linear deformation at a finite strain rate: (i) simple shearing, (ii) pure shearing, and (iii) isotropic compression. We find that the fluctuations induced by such deformations lead to depletion forces that cause rods to group in parallel oriented clusters for the cases of simple and pure shear, but not for isotropic compression. For simple shearing, we find that as the fraction of rods increases, this clustering increases, leading to an increase in the average rate of rotation of the rods, and a decrease in the magnitude of their nematic ordering.

Self-organization in brain tumors: How cell morphology and cell density influence glioma pattern formation

Modeling cancer cells is essential to better understand the dynamic nature of brain tumors and glioma cells, including their invasion of normal brain. Our goal is to study how the morphology of the glioma cell influences the formation of patterns of collective behavior such as flocks (cells moving in the same direction) or streams (cells moving in opposite direction) referred to as oncostream. We have observed experimentally that the presence of oncostreams correlates with tumor progression. We propose an original agent-based model that considers each cell as an ellipsoid. We show that stretching cells from round to ellipsoid increases stream formation. A systematic numerical investigation of the model was implemented in [Formula: see text]. We deduce a phase diagram identifying key regimes for the dynamics (e.g. formation of flocks, streams, scattering). Moreover, we study the effect of cellular density and show that, in contrast to classical models of flocking, increasing cellular density reduces the formation of flocks. We observe similar patterns in [Formula: see text] with the noticeable difference that stream formation is more ubiquitous compared to flock formation.

Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Spatial structure and correlations

We use numerical simulations to study the flow of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state simple shear applied at a fixed volume and a fixed finite strain rate γ[over ̇]. Energy dissipation is via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension with Newtonian rheology. We study the resulting spatial structure of the sheared system, and compute correlation functions of the velocity, the particle density, the nematic order parameter, and the particle angular velocity. Correlations of density, nematic order, and angular velocity are shown to be short ranged both below and above jamming. We compare a system of size-bidisperse particles with a system of size-monodisperse particles, and argue how differences in spatial order as the packing increases lead to differences in the global nematic order parameter. We consider the effect of shearing on initially well ordered configurations, and show that in many cases the shearing acts to destroy the order, leading to the same steady-state ensemble as found when starting from random initial configurations.

Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Particle rotations and orientational ordering

We use numerical simulations to study the flow of a bidisperse mixture of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state simple shear applied at a fixed volume and a fixed finite strain rate γ[over ̇]. Energy dissipation is via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension with Newtonian rheology. Considering a range of packing fractions ϕ and particle asphericities α at small γ[over ̇], we study the angular rotation θ[over ̇]_{i} and the nematic orientational ordering S_{2} of the particles induced by the shear flow, finding a nonmonotonic behavior as the packing ϕ is varied. We interpret this nonmonotonic behavior as a crossover from dilute systems at small ϕ, where single-particle-like behavior occurs, to dense systems at large ϕ, where the geometry of the dense packing dominates and a random Poisson-like process for particle rotations results. We also argue that the finite nematic ordering S_{2} is a consequence of the shearing serving as an ordering field, rather than a result of long-range cooperative behavior among the particles. We arrive at these conclusions by consideration of (i) the distribution of waiting times for a particle to rotate by π, (ii) the behavior of the system under pure, as compared to simple, shearing, (iii) the relaxation of the nematic order parameter S_{2} when perturbed away from the steady state, and (iv) by construction, a numerical mean-field model for the rotational motion of a particle. Our results also help to explain the singular behavior observed when taking the α→0 limit approaching circular disks.

Rheology in dense assemblies of spherocylinders: Frictional vs. frictionless

Using molecular dynamics simulations, we study the steady shear flow of dense assemblies of anisotropic spherocylindrical particles of varying aspect ratios. Comparing frictionless and frictional particles we discuss the specific role of frictional inter-particle forces for the rheological properties of the system. In the frictional system we evidence a shear-thickening regime, similar to that for spherical particles. Furthermore, friction suppresses the alignment of the spherocylinders along the flow direction. Finally, the jamming density in frictional systems is rather insensitive to variations in aspect ratio, quite contrary to what is known from frictionless systems.

Effects of particle shape mixture on strength and structure of sheared granular materials

Using bi-dimensional discrete element simulations, the shear strength and microstructure of granular mixtures composed of particles of different shapes are systematically analyzed as a function of the proportion of grains of a given number of sides and the combination of different shapes (species) in one sample. We varied the angularity of the particles by varying the number of sides of the polygons from 3 (triangles) up to 20 (icosagons) and disks. The samples analyzed were built keeping in mind the following cases: (1) increase of angularity and species starting from disks; (2) decrease of angularity and increase of species starting from triangles; (3) random angularity and increase of species starting from disks and from polygons. The results show that the shear strength vary monotonically with increasing numbers of species (it may increase or decrease), even in the random mixtures (case 3). At the micro-scale, the variation in shear strength as a function of the number of species is due to different mechanisms depending on the cases analyzed. It may result from the increase of both the geometrical and force anisotropies, from only a decrease of frictional anisotropy, or from compensation mechanisms involving geometrical and force anisotropies.

Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Stress, jamming, and contacts

We use numerical simulations to study the flow of a bidisperse mixture of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state shear strain applied at a fixed finite rate. Energy dissipation occurs via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension and resulting in a Newtonian rheology. We study the resulting pressure p and deviatoric shear stress σ of the interacting spherocylinders as a function of packing fraction ϕ, strain rate γ[over ̇], and a parameter α that measures the asphericity of the particles; α is varied to consider the range from nearly circular disks to elongated rods. We consider the direction of anisotropy of the stress tensor, the macroscopic friction μ=σ/p, and the divergence of the transport coefficient η_{p}=p/γ[over ̇] as ϕ is increased to the jamming transition ϕ_{J}. From a phenomenological analysis of Herschel-Bulkley rheology above jamming, we estimate ϕ_{J} as a function of asphericity α and show that the variation of ϕ_{J} with α is the main cause for differences in rheology as α is varied; when plotted as ϕ/ϕ_{J}, rheological curves for different α qualitatively agree. However, a detailed scaling analysis of the divergence of η_{p} for our most elongated particles suggests that the jamming transition of spherocylinders may be in a different universality class than that of circular disks. We also compute the number of contacts per particle Z in the system and show that the value at jamming Z_{J} is a nonmonotonic function of α that is always smaller than the isostatic value. We measure the probability distribution of contacts per unit surface length P(ϑ) at polar angle ϑ with respect to the spherocylinder spine and find that as α→0 this distribution seems to diverge at ϑ=π/2, giving a finite limiting probability for contacts on the vanishingly small flat sides of the spherocylinder. Finally, we consider the variation of the average contact force as a function of location on the particle surface.

Orientational Ordering in Athermally Sheared, Aspherical, Frictionless Particles

We numerically simulate the uniform athermal shearing of bidisperse, frictionless, two-dimensional spherocylinders and three-dimensional prolate ellipsoids. We focus on the orientational ordering of particles as an asphericity parameter α→0 and particles approach spherical. We find that the nematic order parameter S_{2} is nonmonotonic in the packing fraction ϕ and that, as α→0, S_{2} stays finite at jamming and above. The approach to spherical particles thus appears to be singular. We also find that sheared particles continue to rotate above jamming and that particle contacts preferentially lie along the narrowest width of the particles, even as α→0.

Anisotropic particles strengthen granular pillars under compression

We probe the effects of particle shape on the global and local behavior of a two-dimensional granular pillar, acting as a proxy for a disordered solid, under uniaxial compression. This geometry allows for direct measurement of global material response, as well as tracking of all individual particle trajectories. In general, drawing connections between local structure and local dynamics can be challenging in amorphous materials due to lower precision of atomic positions, so this study aims to elucidate such connections. We vary local interactions by using three different particle shapes: discrete circular grains (monomers), pairs of grains bonded together (dimers), and groups of three bonded in a triangle (trimers). We find that dimers substantially strengthen the pillar and the degree of this effect is determined by orientational order in the initial condition. In addition, while the three particle shapes form void regions at distinct rates, we find that anisotropies in the local amorphous structure remain robust through the definition of a metric that quantifies packing anisotropy. Finally, we highlight connections between local deformation rates and local structure.

Rheology of dense granular flows for elongated particles

We study the rheology of dense granular flows for frictionless spherocylinders by means of 3D numerical simulations. As in the case of spherical particles, the effective friction μ is an increasing function of the inertial number I, and we systematically investigate the dependence of μ on the particle aspect ratio Q, as well as that of the normal stress differences, the volume fraction, and the coordination number. We show in particular that the quasistatic friction coefficient is nonmonotonic with Q: from the spherical case Q=1, it first sharply increases, reaches a maximum around Q≃1.05, and then gently decreases until Q=3, passing its initial value at Q≃2. We provide a microscopic interpretation for this unexpected behavior through the analysis of the distribution of dissipative contacts around the particles: as compared to spheres, slightly elongated grains enhance contacts in their central cylindrical band, whereas at larger aspect ratios particles tend to align and dissipate by preferential contacts at their hemispherical caps.

Experimental Study of Ordering of Hard Cubes by Shearing

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.

Analyzing X-ray tomographies of granular packings

Starting from three-dimensional volume data of a granular packing, as, e.g., obtained by X-ray Computed Tomography, we discuss methods to first detect the individual particles in the sample and then analyze their properties. This analysis includes the pair correlation function, the volume and shape of the Voronoi cells, and the number and type of contacts formed between individual particles. We mainly focus on packings of monodisperse spheres, but we will also comment on other monoschematic particles such as ellipsoids and tetrahedra. This paper is accompanied by a package of free software containing all programs (including source code) and an example three-dimensional dataset which allows the reader to reproduce and modify all examples given.

Binary mixtures of disks and elongated particles: Texture and mechanical properties

We analyze the shear strength and microstructure of binary granular mixtures consisting of disks and elongated particles by varying systematically both the mixture ratio and degree of homogeneity (from homogeneous to fully segregated). The contact dynamics method is used for numerical simulations with rigid particles interacting by frictional contacts. A counterintuitive finding of this work is that the shear strength, packing fraction, and, at the microscopic scale, the fabric, force, and friction anisotropies of the contact network are all nearly independent of the degree of homogeneity. In other words, homogeneous mixtures have the same strength properties as segregated packings of the two particle shapes. In contrast, the shear strength increases with the proportion of elongated particles correlatively with the increase of the corresponding force and fabric anisotropies. By a detailed analysis of the contact network topology, we show that various contact types contribute differently to force transmission and friction mobilization.

Influence of Coalescence on the Anisotropic Mechanical and Electrical Properties of Nickel Powder/Polydimethylsiloxane Composites

Multifunctional polymer-based composites have been widely used in various research and industrial applications, such as flexible and stretchable electronics and sensors and sensor-integrated smart structures. This study investigates the influence of particle coalescence on the mechanical and electrical properties of spherical nickel powder (SNP)/polydimethylsiloxane (PDMS) composites in which SNP was aligned using an external magnetic field. With the increase of the volume fraction of the SNP, the aligned SNP/PDMS composites exhibited a higher tensile strength and a lower ultimate strain. In addition, the composites with aligned SNP showed a lower percolation threshold and a higher electrical conductivity compared with those with randomly dispersed SNP. However, when the concentration of the SNP reached a certain level (40 vol. %), the anisotropy of the effective material property became less noticeable than that of the lower concentration (20 vol. %) composites due to the change of the microstructure of the particles caused by the coalescence of the particles at a high concentration. This work may provide rational methods for the fabrication of aligned composites.

Metastable orientational order of colloidal discoids

The interplay between phase separation and kinetic arrest is important in supramolecular self-assembly, but their effects on emergent orientational order are not well understood when anisotropic building blocks are used. Contrary to the typical progression from disorder to order in isotropic systems, here we report that colloidal oblate discoids initially self-assemble into short, metastable strands with orientational order—regardless of the final structure. The model discoids are suspended in a refractive index and density-matched solvent. Then, we use confocal microscopy experiments and Monte Carlo simulations spanning a broad range of volume fractions and attraction strengths to show that disordered clusters form near coexistence boundaries, whereas oriented strands persist with strong attractions. We rationalize this unusual observation in light of the interaction anisotropy imparted by the discoids. These findings may guide self-assembly for anisotropic systems in which orientational order is desired, such as when tailored mechanical properties are sought.

The pathways available for self-assembly are affected by the shape anisotropy of the building blocks, but the details are still unclear. Here, Hsiao et al. show that colloidal discoids self-assemble into metastable states with orientational order when kinetic trapping is incorporated as a design principle.

Heaping, secondary flows and broken symmetry in flows of elongated granular particles

In this paper we report experiments where we shear granular rods in split-bottom geometries, and find that a significant heap of height of least 40% of the filling height can form at the particle surface. We show that heaping is caused by a significant secondary flow, absent for spherical particles. Flow reversal transiently reverses the secondary flow, leading to a quick collapse and slower regeneration of the heap. We present a symmetry argument and experimental data that show that the generation of the secondary flow is driven by a misalignment of the mean particle orientation with the streamlines of the flow. This general mechanism is expected to be important in all flows of sufficiently anisometric grains.

Particle alignment and clustering in sheared granular materials composed of platy particles

By means of molecular dynamics simulations, we investigate the texture and local ordering in sheared packings composed of cohesionless platy particles. The morphology of large packings of platy particles in quasistatic equilibrium is complex due to the combined effects of local nematic ordering of the particles and anisotropic orientations of contacts between particles. We find that particle alignment is strongly enhanced by the degree of platyness and leads to the formation of face-connected clusters of exponentially decaying size. Interestingly, due to dynamics in continuous shearing, this ordering phenomenon emerges even in systems composed of particles of very low platyness differing only slightly from spherical shape. The number of clusters is an increasing function of platyness. However, at high platyness the proportion of face-face interactions is too low to allow for their percolation throughout the system.

Evolution of shear zones in granular materials

The evolution of wide shear zones or shear bands was investigated experimentally and numerically for quasistatic dry granular flows in split bottom shear cells. We compare the behavior of materials consisting of beads, irregular grains, such as sand, and elongated particles. Shearing an initially random sample, the zone width was found to significantly decrease in the first stage of the process. The characteristic shear strain associated with this decrease is about unity and it is systematically increasing with shape anisotropy, i.e., when the grain shape changes from spherical to irregular (e.g., sand) and becomes elongated (pegs). The strongly decreasing tendency of the zone width is followed by a slight increase which is more pronounced for rodlike particles than for grains with smaller shape anisotropy (beads or irregular particles). The evolution of the zone width is connected to shear-induced packing density change and for nonspherical particles it also involves grain reorientation effects. The final zone width is significantly smaller for irregular grains than for spherical beads.

Effects of grain shape on packing and dilatancy of sheared granular materials

A granular material exposed to shear shows a variety of unique phenomena: Reynolds dilatancy, positional order and orientational order effects may compete in the shear zone. We study granular packing consisting of macroscopic prolate, oblate and spherical grains and compare their behaviour. X-ray tomography is used to determine the particle positions and orientations in a cylindrical split bottom shear cell. Packing densities and the arrangements of individual particles in the shear zone are evaluated. For anisometric particles, we observe the competition of two opposite effects. On the one hand, the sheared granules are dilated, on the other hand the particles reorient and align with respect to the streamlines. Even though aligned cylinders in principle may achieve higher packing densities, this alignment compensates for the effect of dilatancy only partially. The complex rearrangements lead to a depression of the surface above the well oriented region while neighbouring parts still show the effect of dilation in the form of heaps. For grains with isotropic shapes, the surface remains rather flat. Perfect monodisperse spheres crystallize in the shear zone, whereby positional order partially overcompensates dilatancy effects. However, even slight deviations from the ideal monodisperse sphere shape inhibit crystallization.

Packings of irregular polyhedral particles: strength, structure, and effects of angularity

We present a systematic numerical investigation of the shear strength and structure of granular packings composed of irregular polyhedral particles. The angularity of the particles is varied by increasing the number of faces from 8 (octahedronlike shape) to 596. We find that the shear strength increases with angularity up to a maximum value and saturates as the particles become more angular (below 46 faces). At the same time, the packing fraction increases to a peak value but declines for more angular particles. We analyze the connectivity and anisotropy of the microstructure by considering both the contacts and branch vectors joining particle centers. The increase of the shear strength with angularity is shown to be due to a net increase of the fabric and force anisotropies but at higher particle angularity a rapid falloff of the fabric anisotropy is compensated by an increase of force anisotropy, leading thus to the saturation of shear strength.

Granular gases of rod-shaped grains in microgravity

Granular gases are convenient model systems to investigate the statistical physics of nonequilibrium systems. In the literature, one finds numerous theoretical predictions, but only few experiments. We study a weakly excited dilute gas of rods, confined in a cuboid container in microgravity during a suborbital rocket flight. With respect to a gas of spherical grains at comparable filling fraction, the mean free path is considerably reduced. This guarantees a dominance of grain-grain collisions over grain-wall collisions. No clustering was observed, unlike in similar experiments with spherical grains. Rod positions and orientations were determined and tracked. Translational and rotational velocity distributions are non-Gaussian. Equipartition of kinetic energy between translations and rotations is violated.