Orientational order and alignment of elongated particles induced by shear

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.

SonoPrint: Acoustically Assisted Volumetric 3D Printing for Composites

Advances in additive manufacturing in composites have transformed aerospace, medical devices, tissue engineering, and electronics. A key aspect of enhancing properties of 3D-printed objects involves fine-tuning the material by embedding and orienting reinforcement within the structure. Existing methods for orienting these reinforcements are limited by pattern types, alignment, and particle characteristics. Acoustics offers a versatile method to control the particles independent of their size, geometry, and charge, enabling intricate pattern formations. However, integrating acoustics into 3D printing has been challenging due to the scattering of the acoustic field between polymerized layers and unpolymerized resin, resulting in unwanted patterns. To address this challenge, SonoPrint, an innovative acoustically assisted volumetric 3D printer is developed that enables simultaneous reinforcement patterning and printing of the entire structure. SonoPrint generates mechanically tunable composite geometries by embedding reinforcement particles, such as microscopic glass, metal, and polystyrene, within the fabricated structure. This printer employs a standing wave field to create targeted particle motifs-including parallel lines, radial lines, circles, rhombuses, hexagons, and polygons-directly in the photosensitive resin, completing the print in just a few minutes. SonoPrint enhances structural properties and promises to advance volumetric printing, unlocking applications in tissue engineering, biohybrid robots, and composite fabrication.

Evident structural anisotropies arising from near-zero particle asphericity in granular spherocylinder packings

With magnetic resonance imaging experiments, we study packings of granular spherocylinders with merely 2% asphericity. Evident structural anisotropies across all length scales are identified. Most interestingly, the global nematic order decreases with increasing packing fraction, while the local contact anisotropy shows an opposing trend. We attribute this counterintuitive phenomenon to a competition between gravity-driven ordering aided by frictional contacts and a geometric frustration effect at the marginally jammed state. It is also surprising to notice that such slight particle asphericity can trigger non-negligible correlations between contact-level and mesoscale structures, manifested in drastically different nonaffine structural rearrangements upon compaction from that of granular spheres. These observations can help improve statistical mechanical models for the orientational order transformation of nonspherical granular particle packings, which involves complex interplays between particle shape, frictional contacts, and external force field.

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.

Orientation of graphene nanosheets in suspension under an electric field: theoretical model and molecular dynamic simulations

Orientation regulation of nanoparticles in a suspension by an electric field is a powerful tool to tune its mechanical, thermal, optical, electrical properties etc. However, how molecular modification can affect the orientation of two-dimensional nanoparticles is still unclear. In this paper, the influence of molecular modification on the orientation of graphene nanosheets (GNS) in water was investigated through theoretical analyses and molecular dynamics (MD) simulations. Firstly, a new orientation angle model was proposed, which considers hydration effects, dipole moments and resistance torque. Then, MD simulations were conducted to investigate the effects of position, direction, type, and number of functional groups on the orientation of GNS. The trend observed in MD simulations is consistent with the proposed theoretical model. The results reveal that, under the combined influence of the dipole moment and hydration effects, the modification with hydrophilic functional groups can reduce the orientation angle from 21.31° to 8.34°, while the modification with hydrophobic functional groups increases it to 26.43°. Among the hydrophilic functional groups, orientation of hydroxylated GNS is the best. With an increase in the number of hydroxyl groups, orientation angle is decreased from 12.61° to 8.34°. This work can provide valuable guidance for the design of high-performance suspensions and composites, such as thermal smart materials with adjustable thermal conductivity and intelligent devices with tailored capabilities.

Dynamic phases and combing effects for elongated particles moving over quenched disorder

We consider a two-dimensional system of elongated particles driven over a landscape containing randomly placed pinning sites. For varied pinning site density, external drive magnitude, and particle elongation, we find a wide variety of dynamic phases, including random structures, stripe or combed phases with nematic order, and clogged states. The different regimes can be identified by examining nematic ordering, cluster size, number of pinned particles, and transverse diffusion. In some regimes we find that the pinning can enhance the particle alignment, producing a nonmonotonic signature in the nematic ordering with a maximum at a particular combination of pinning density and drive. The optimal nematic occurs when a sufficient number of particles can be pinned, generating a local shear and leading to what we call a combing effect. At high drives, the combing effect is reduced when the number of pinned particles decreases. For stronger pinning, the particles form a heterogeneous clustered or clogged state that depins into a fluctuating state with high diffusion.

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.

Orientational-induced strain hardening of axisymmetric grains

The rheological response of oriented axisymmetric grains has additional degrees of complexity associated with their microstructure orientation. These additional kinematic degrees of freedom that give rise to complex transient macroscale rheological responses are not well understood. In this Letter, we study the rheology of axisymmetric grains subjected to transient flow. We identify strong coupling between the microstructure rearrangement and strain hardening which, under certain conditions, can yield jamming. We identify the critical conditions corresponding to jamming and the dependency on the shape of the grains. It is shown that this is a particular form of jamming that is directional in nature, since unjamming occurs if the shear direction is reversed.

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.

Structural evolution of granular cubes packing during shear-induced ordering

Packings of granular particles may transform into ordered structures under external agitation, which is a special type of out-of-equilibrium self-assembly. Here, evolution of the internal packing structures of granular cubes under cyclic rotating shearing has been analyzed using magnetic resonance imaging techniques. Various order parameters, different types of contacts and clusters composed of face-contacting cubes, as well as the free volume regions in which each cube can move freely have been analyzed systematically to quantify the ordering process and the underlying mechanism of this granular self-assembly. The compaction process is featured by a first rapid formation of orientationally ordered local structures with faceted contacts, followed by further densification driven by free-volume maximization with an almost saturated degree of order. The ordered structures are strongly anisotropic with contacting ordered layers in the vertical direction while remaining liquid-like in the horizontal directions. Therefore, the constraint of mechanical stability for granular packings and the thermodynamic principle of entropy maximization are both effective in this system, which we propose can be reconciled by considering different depths of supercooling associated with various degrees of freedom.

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.

Discharge of elongated grains from silo with rotating bottom

We study the flow of elongated grains (wooden pegs of length L=20 mm with circular cross section of diameter d_{c}=6 and 8 mm) from a silo with a rotating bottom and a circular orifice of diameter D. In the small orifice range (D/d<5) clogs are mostly broken by the rotating base, and the flow is intermittent with avalanches and temporary clogs. Here d≡(3/2d_{c}^{2}L)^{1/3} is the effective grain diameter. Unlike for spherical grains, for rods the flow rate W clearly deviates from the power law dependence W∝(D-kd)^{2.5} at lower orifice sizes in the intermittent regime, where W is measured in between temporary clogs only. Instead, below about D/d<3 an exponential dependence W∝e^{κD} is detected. Here k and κ are constants of order unity. Even more importantly, rotating the silo base leads to a strong-more than 50%-decrease of the flow rate, which otherwise does not depend significantly on the value of ω in the continuous flow regime. In the intermittent regime, W(ω) appears to follow a nonmonotonic trend, although with considerable noise. A simple picture, in terms of the switching from funnel flow to mass flow and the alignment of the pegs due to rotation, is proposed to explain the observed difference between spherical and elongated grains. We also observe shear-induced orientational ordering of the pegs at the bottom such that their long axes in average are oriented at a small angle 〈θ〉≈15^{∘} to the motion of the bottom.

Measuring the order parameter of vertically aligned nanorod assemblies

Vertically aligned nanorod assemblies are of great interest both for fundamental studies of anisotropic physical properties arising from the structures and for the development of functional devices utilizing such anisotropic characteristics. Simultaneous measurement of the homeotropic order parameter (Shomeo) of assemblies in dynamic states can allow further optimization of the assembly process and the device performance. Although many techniques (e.g. birefringence measurement, SAXS analysis, and high-resolution microscopy) have been proposed to characterise Shomeo, these do not yet meet the essential criteria such as for rapid, in situ and non-destructive analyses. Here, we propose a novel approach employing a unique photoluminescence behaviour of lanthanide-doped crystalline nanorods, of which the emission spectrum contains the detailed information on the structure of the assembly. We demonstrate a rapid in situ determination of Shomeo of Eu3+-doped NaYF4 nanorods of which the orientation is controlled under an external electric field. The method does not require the consideration of polarization and can be performed using a conventional fluorescence microscopy setup. This new methodology would provide a more in-depth examination of various assembled nanostructures and the collective dynamics of their building blocks.

Evolution of shear zones in granular packings under pressure

Stress transmission in realistic granular media often occurs under external load and in the presence of boundary slip. We investigate shear localization in a split-bottom Couette cell with smooth walls subject to a confining pressure experimentally and by means of numerical simulations. We demonstrate how the characteristics of the shear zone, such as its center position and width, evolve as the confining pressure and wall slip modify the local effective friction coefficient of the material. For increasing applied pressure, the shear zone evolves toward the center of the cylinder and grows wider and the angular velocity reduces compared to the driving rate of the bottom disk. Moreover, the presence of slip promotes the transition from open shear zones at the top surface to closed shear zones inside the bulk. We also systematically vary the ratio of the effective friction near the bottom plate and in the bulk in simulations and observe the resulting impact on the surface flow profile. Besides the boundary conditions and external load, material properties such as grain size are also known to influence the effective friction coefficient. However, our numerical results reveal that the center position and width of the shear zone are insignificantly affected by the choice of the grain size as far as it remains small compared to the radius of the rotating bottom disk.

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.

Packings of frictionless spherocylinders

We present simulation results on the properties of packings of frictionless spherocylindrical particles. Starting from a random distribution of particles in space, a packing is produced by minimizing the potential energy of interparticle contacts until a force-equilibrated state is reached. For different particle aspect ratios α=10⋯40, we calculate contacts z, pressure as well as bulk and shear modulus. Most important is the fraction f_{0}(α) of spherocylinders with contacts at both ends, as it governs the jamming threshold z_{c}(α)=8+2f_{0}(α). These results highlight the important role of the axial "sliding" degree of freedom of a spherocylinder, which is a zero-energy mode but only if no end contacts are present.

Discrete Element Method Investigation of Binary Granular Flows with Different Particle Shapes

The effects of particle shape differences on binary mixture shear flows are investigated using the Discrete Element Method (DEM). The binary mixtures consist of frictionless rods and disks, which have the same volume but significantly different shapes. In the shear flows, stacking structures of rods and disks are formed. The effects of the composition of the mixture on the stacking are examined. It is found that the number fraction of stacking particles is smaller for the mixtures than for the monodisperse rods and disks. For binary mixtures with small particle shape differences, the mixture stresses are bounded by the stresses of the two corresponding monodisperse systems. However, for binary mixtures with large particle shape differences, the stresses of the mixtures can be larger than the stresses of the monodisperse systems at large solid volume fractions because larger differences in particle shape cause geometrical interference in packing, leading to stronger particle–particle interactions in the flow. The stresses in dense binary mixtures are found to be exponential functions of the order parameter, which is a measure of particle alignment. Based on the simulation results, an empirical expression for the bulk friction coefficient (ratio of the shear stress to normal stress) for dense binary flows is proposed by accounting for the effects of particle alignment and solid volume fraction.

Role of particle orientational order during shear thickening in suspensions of colloidal rods

Rheology of dense anisotropic colloidal suspensions often exhibits unsteady flow at constant imposed shear stress and/or shear rate. Using simultaneous high-resolution confocal microscopy and rheology, we find that the temporal behavior arises due to a strong coupling between shear flow and particle orientation. At smaller applied stresses, the orientation of rods fluctuates around the flow direction. A transition to an intermittent disordered state is observed at higher stresses when the angle between the flow and the rod orientation reaches a critical value. This disordered state is associated with transient drop in shear rate and an increase in viscosity. Simultaneous visualization of boundary stresses and orientation shows that the disordered regions lead to heterogeneous stresses and positive normal forces at the boundary, indicating the formation of systems spanning disordered particle contact networks.

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.

Extracellular matrix anisotropy is determined by TFAP2C-dependent regulation of cell collisions

The isotropic or anisotropic organization of biological extracellular matrices has important consequences for tissue function. We study emergent anisotropy using fibroblasts that generate varying degrees of matrix alignment from uniform starting conditions. This reveals that the early migratory paths of fibroblasts are correlated with subsequent matrix organization. Combined experimentation and adaptation of Vicsek modelling demonstrates that the reorientation of cells relative to each other following collision plays a role in generating matrix anisotropy. We term this behaviour 'cell collision guidance'. The transcription factor TFAP2C regulates cell collision guidance in part by controlling the expression of RND3. RND3 localizes to cell-cell collision zones where it downregulates actomyosin activity. Cell collision guidance fails without this mechanism in place, leading to isotropic matrix generation. The cross-referencing of alignment and TFAP2C gene expression signatures against existing datasets enables the identification and validation of several classes of pharmacological agents that disrupt matrix anisotropy.

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.

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.

Nanocombing Effect Leads to Nanowire-Based, in-Plane, Uniaxial Thin Films

The fabrication of thin films comprising ordered nanowire assemblies with emerging, precisely defined properties and adjustable functionalities enables highly integrated technologies in the fields of microelectronics and micro system technology, as well as for efficient power generation, storage, and utilization. Shear force, theoretically, is deemed the most promising method for obtaining in-plane, uniaxial thin films comprising nanowires. The success depends largely on the assembly process, and uniform structural control throughout multiple length scales can be achieved only if a rational strategy is executed. Here, we report that in shearing processes dopants such as lyotropic cellulose nanorods can give rise to the uniaxial alignment of V₂O₅· nH₂O nanowires. Our systematic study indicates that this finding, namely, the nanocombing effect, can be a general principle for the continuous production of uniaxial thin films comprising densely packed nanowires varying in chemical composition and aspect ratios. Conversion of the V₂O₅· nH₂O constituents via in situ oxidative polymerization leads to in-plane, uniaxial polyaniline (PANI) thin films with anisotropic electric and optical properties, which are otherwise difficult to fabricate due to the poor processability of PANI. The uniaxial PANI thin films can be utilized to fabricate flexible gas sensors for distinguishing various analytes, including similar homologues such as methanol and ethanol. We expect the methodology to be applied to a broad spectrum of synthetic and biogenic nanowires for the integration of their collective properties in high-performance electronic devices.

Translational and Rotational Dynamical Heterogeneities in Granular Systems

We use x-ray tomography to investigate the translational and rotational dynamical heterogeneities of a three dimensional hard ellipsoid granular packing driven by oscillatory shear. We find that particles which translate quickly form clusters with a size distribution given by a power law with an exponent that is independent of the strain amplitude. Identical behavior is found for particles that are translating slowly, rotating quickly, or rotating slowly. The geometrical properties of these four different types of clusters are the same as those of random clusters. Different cluster types are considerably correlated or anticorrelated, indicating a significant coupling between translational and rotational degrees of freedom. Surprisingly, these clusters are formed already at time scales that are much shorter than the α-relaxation time, in stark contrast to the behavior found in glass-forming systems.

Flow of anisometric particles in a quasi-two-dimensional hopper

The stationary flow field in a quasi-two-dimensional hopper is investigated experimentally. The behavior of materials consisting of beads and elongated particles with different aspect ratio is compared. We show, that while the vertical velocity in the flowing region can be fitted with a Gaussian function for beads, in the case of elongated grains the flowing channel is narrower and is bordered with sharper velocity gradient. For this case, we quantify deviations from the Gaussian velocity profile. Relative velocity fluctuations are considerably larger and slower for elongated grains.

Kinematic Model of Transient Shape-Induced Anisotropy in Dense Granular Flow

Nonspherical particles are ubiquitous in nature and industry, yet previous theoretical models of granular media are mostly limited to systems of spherical particles. The problem is that in systems of nonspherical anisotropic particles, dynamic particle alignment critically affects their mechanical response. To study the tendency of such particles to align, we propose a simple kinematic model that relates the flow to the evolution of particle alignment with respect to each other. The validity of the proposed model is supported by comparison with particle-based simulations for various particle shapes ranging from elongated rice-like (prolate) to flattened lentil-like (oblate) particles. The model shows good agreement with the simulations for both steady-state and transient responses, and advances the development of comprehensive constitutive models for shape-anisotropic particles.

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.

Dynamic X-ray radiography reveals particle size and shape orientation fields during granular flow

When granular materials flow, the constituent particles segregate by size and align by shape. The impacts of these changes in fabric on the flow itself are not well understood, and thus novel non-invasive means are needed to observe the interior of the material. Here, we propose a new experimental technique using dynamic X-ray radiography to make such measurements possible. The technique is based on Fourier transformation to extract spatiotemporal fields of internal particle size and shape orientation distributions during flow, in addition to complementary measurements of velocity fields through image correlation. We show X-ray radiography captures the bulk flow properties, in contrast to optical methods which typically measure flow within boundary layers, as these are adjacent to any walls. Our results reveal the rich dynamic alignment of particles with respect to streamlines in the bulk during silo discharge, the understanding of which is critical to preventing destructive instabilities and undesirable clogging. The ideas developed in this paper are directly applicable to many other open questions in granular and soft matter systems, such as the evolution of size and shape distributions in foams and biological materials.

Collisional dissipation rate in shearing flows of granular liquid crystals

We make use of discrete-element-method numerical simulations of inelastic frictionless cylinders in simple shearing at different length-to-diameter ratios and solid volume fractions to analyze the rate of collisional dissipation of the fluctuation kinetic energy. We show that the nonspherical geometry of the particles is responsible, by inducing rotation, for increasing the dissipation rate of the fluctuation kinetic energy with respect to that for frictionless spheres. We also suggest that the partial alignment of the cylinders induced by shearing concurs with the particle inelasticity in generating correlation in the velocity fluctuations and thus affecting the collisional dissipation rate as the solid volume fraction increases. Finally, we propose simple phenomenological modifications to the expression of the collisional dissipation rate of kinetic theory of granular gases to take into account our findings.

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.