BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Leonardo E. Silbert

Leonardo E. Silbert appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechcond-mat.mtrl-scicond-mat.dis-nnphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

18works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2026arXiv

Reverse segregation and self-organization in inclined chute flows of bidisperse granular mixtures

In the usual segregation scenario for stable inclined chute flows of bidisperse mixtures of fine and coarse spherical particles, coarse particles rise toward the free surface, forming a coarse-rich region atop the flowing pile. Beyond a threshold coarse-to-fine diameter ratio of approximately 4, conversely, the weight of the coarse particles exceeds the segregation driving forces, causing individual coarse particles to sink within the pile and producing a reversed segregation state. However, an understanding of the collective evolution of the pile structure is still lacking when the particle diameter ratio exceeds 4 {\textit{and}} the coarse particle mass fraction is appreciable. To explore this broadly bidisperse limit, we perform discrete element method simulations considering mean particle diameter ratios of up to 8 and coarse particle mass fractions spanning 0.1 to 0.9. The steady-state flow profiles reveal several intriguing behaviors that depend on the diameter ratio and mass fraction. These include a previously identified transition from usual to reverse segregation and a newfound tendency to self-organize into alternating coarse- and fine-rich particle layers stacked along the shear gradient direction, with layer thickness dictated by the coarse particle diameter. A fuller understanding of segregation at this scale could pave the way for enhanced mixing or demixing techniques at the commercial scale.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Flow and Arrest in Stressed Granular Materials

Flowing granular materials often abruptly arrest if not driven by sufficient applied stresses. Such abrupt cessation of motion can be economically expensive in industrial materials handling and processing, and is significantly consequential in intermittent geophysical phenomena such as landslides and earthquakes. Using discrete element simulations, we calculate states of steady flow and arrest for granular materials under the conditions of constant applied pressure and shear stress, which are also most relevant in practice. Here the material can dilate or compact, and flow or arrest, in response to the applied stress. Our simulations highlight that under external stress, the intrinsic response of granular materials is characterized by uniquely-defined steady states of flow or arrest, which are highly sensitive to interparticle friction. While the flowing states can be equivalently characterized by volume fraction, coordination number or internal stress ratio, to characterize the states of shear arrest, one needs to also consider the structural anisotropy in the contact network. We highlight the role of dilation in the flow-arrest transition, and discuss our findings in the context of rheological transitions in granular materials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fluctuations and power-law scaling of dry, frictionless granular rheology near the hard-particle limit

The flow of frictionless granular particles is studied with stress-controlled discrete element modeling simulations for systems varying in size from 300 to 100,000 particles. The volume fraction and shear stress ratio $μ$ are relatively insensitive to system size fo a wide range of inertial numbers $I$. Second-order effects in strain rate, such as second normal stress differences, require large system sizes to accurately extract meaningful results, notably a non-monotonic dependence in the first normal stress difference with strain rate. The first-order rheological response represented by the $μ(I)$ relationship works well at describing the lower-order aspects of the rheology, except near the quasi-static limit of these stress-controlled flows. The pressure is varied over five decades, and a pressure dependence of the coordination number is observed, which is not captured by the inertial number. Large fluctuations observed for small systems $N\le$ 1,000 near the quasi-static limit can lead to arrest of flow resulting in challenges to fitting the data to rheological relationships. The inertial number is also insufficient for capturing the pressure-dependent behavior of property fluctuations. Fluctuations in the flow and microstructural properties are measured in both the quasi-static and inertial regimes, including shear stress, pressure, strain rate, normal stress differences, volume fraction, coordination number and contact fabric anisotropy. The fluctuations in flow properties scale self-similarly with pressure and system size. A transition in the scaling of fluctuations of stress properties and contact fabric anisotropy are measured and proposed as a quantitative identification of the transition from inertial to quasi-static flow.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Large-scale frictionless jamming with power-law particle size distributions

Due to significant computational expense, discrete element method simulations of jammed packings of size-dispersed spheres with size ratios greater than 1:10 have remained elusive, limiting the correspondence between simulations and real-world granular materials with large size dispersity. Invoking a recently developed neighbor binning algorithm, we generate mechanically-stable jammed packings of frictionless spheres with power-law size distributions containing up to nearly four million particles with size ratios up to 1:100. By systematically varying the width and exponent of the underlying power laws, we analyze the role of particle size distributions on the structure of jammed packings. The densest packings are obtained for size distributions that balance the relative abundance of large-large/intermediate and small-small particle contacts. Although the proportion of rattler particles and mean coordination number strongly depend on the size distribution, the mean coordination of non-rattler particles attains the frictionless isostatic value of six in all cases. The size distribution of non-rattler particles that participate in the load-bearing network exhibits no dependence on the width of the total particle size distribution beyond a critical particle size for low-magnitude exponent power laws. This signifies that only particles with sizes greater than the critical particle size contribute to the mechanical stability. However, for high-magnitude exponent power laws, all particle sizes participate in the mechanical stability of the packing.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Evolution of Internal Granular Structure at the Flow-Arrest Transition

The evolution of the internal granular structure in shear-arrested and shear-flowing states of granular materials is characterized using fabric tensors as descriptors of the internal contact and force networks. When a dilute system of frictional grains is subjected to a constant pressure and shear stress, the bulk stress ratio is well-predicted from the anisotropy of its contact and force networks during transient flow prior to steady shear flow or shear arrest. Although the onset of shear arrest is a stochastic process, the fabric tensors upon arrest are distributed around nearly equal contributions of force and contact network anisotropy to the bulk stress ratio. The distribution becomes seemingly narrower with increasing system size. The anisotropy of the contact network in shear-arrested states is reminiscent of the fabric anisotropy observed in shear-jammed packings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Granular packings with sliding, rolling and twisting friction

Intuition tells us that a rolling or spinning sphere will eventually stop due to the presence of friction and other dissipative interactions. The resistance to rolling and spinning/twisting torque that stops a sphere also changes the microstructure of a granular packing of frictional spheres by increasing the number of constraints on the degrees of freedom of motion. We perform discrete element modeling simulations to construct sphere packings implementing a range of frictional constraints under a pressure-controlled protocol. Mechanically stable packings are achievable at volume fractions and average coordination numbers as low as 0.53 and 2.5, respectively, when the particles experience high resistance to sliding, rolling and twisting. Only when the particle model includes rolling and twisting friction, were experimental volume fractions reproduced.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cut-off nonlinearities in the low-temperature vibrations of glasses and crystals

We present a computer simulation study of glassy and crystalline states using the standard Lennard-Jones interaction potential that is truncated at a finite cut-off distance, as is typical of many computer simulations. We demonstrate that the discontinuity at the cut-off distance in the first derivative of the potential (corresponding to the interparticle force) leads to the appearance of cut-off nonlinearities. These cut-off nonlinearities persist into the very-low-temperature regime thereby affecting low-temperature thermal vibrations, which leads to a breakdown of the harmonic approximation for many eigen modes, particularly for low-frequency vibrational modes. Furthermore, while expansion nonlinearities which are due to higher order terms in the Taylor expansion of the interaction potential are usually ignored at low temperatures and show up as the temperature increases, cut-off nonlinearities can become most significant at the lowest temperatures. Anharmonic effects readily show up in the elastic moduli which not only depend on the eigen frequencies, but are crucially sensitive to the eigen vectors of the normal modes. Whereas, those observables that rely mainly on static structural information or just the eigen frequencies, such as the vibrational density of states, total potential energy, and specific heat, show negligible dependence on the presence of the cut-off. Similar aspects of nonlinear behavior have recently been reported in model granular materials, where the constituent particles interact through finite-range, purely-repulsive potentials. These nonlinearities have been ascribed to the nature of the sudden cut-off at contact in the force-law, thus we demonstrate that cut-off nonlinearities emerge as a general feature of ordered and disordered solid state systems interacting through truncated potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Elastic Moduli and Vibrational Modes in Jammed Particulate Packings

When we elastically impose a homogeneous, affine deformation on amorphous solids, they also undergo an inhomogeneous, non-affine deformation, which can have a crucial impact on the overall elastic response. To correctly understand the elastic modulus $M$, it is therefore necessary to take into account not only the affine modulus $M_A$, but also the non-affine modulus $M_N$ that arises from the non-affine deformation. In the present work, we study the bulk ($M=K$) and shear ($M=G$) moduli in static jammed particulate packings over a range of packing fractions $φ$. One novelty of this work is to elucidate the contribution of each vibrational mode to the non-affine $M_N$ through a modal decomposition of the displacement and force fields. In the vicinity of the (un)jamming transition, $φ_{c}$, the vibrational density of states, $g(ω)$, shows a plateau in the intermediate frequency regime above a characteristic frequency $ω^\ast$. We illustrate that this unusual feature apparent in $g(ω)$ is reflected in the behavior of $M_N$: As $φ\rightarrow φ_c$, where $ω^\ast \rightarrow 0$, those modes for $ω< ω^\ast$ contribute less and less, while contributions from those for $ω> ω^\ast$ approach a constant value which results in $M_N$ to approach a critical value $M_{Nc}$, as $M_N-M_{Nc} \sim ω^\ast$. At $φ_c$ itself, the bulk modulus attains a finite value $K_c=K_{Ac}-K_{Nc} > 0$, such that $K_{Nc}$ has a value that remains below $K_{Ac}$. In contrast, for the critical shear modulus $G_c$, $G_{Nc}$ and $G_{Ac}$ approach the same value so that the total value becomes exactly zero, $G_c = G_{Ac}-G_{Nc} =0$. We explore what features of the configurational and vibrational properties cause such the distinction between $K$ and $G$, allowing us to validate analytical expressions for their critical values.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spatial Distributions of Local Elastic Moduli Near the Jamming Transition

Recent progress on studies of the nanoscale mechanical responses in disordered systems has highlighted a strong degree of heterogeneity in the elastic moduli. In this contribution, using computer simulations, we study the elastic heterogeneities in athermal amorphous solids, composed of isotropic, static, sphere packings, near the jamming transition. We employ techniques, based on linear response methods, that are amenable to experimentation. We find that the local elastic moduli are randomly distributed in space and are described by Gaussian probability distributions, thereby lacking any significant spatial correlations, that persists all the way down to the transition point. However, the shear modulus fluctuations grow as the jamming threshold is approached, which is characterized by a new power-law scaling. Through this diverging behavior we are able to identify a characteristic length scale, associated with shear modulus heterogeneities, that distinguishes between bulk and local elastic responses.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Critical-like Features of Stress Response in Frictional Packings

The mechanical response of static, unconfined, overcompressed face centred cubic, granular arrays is studied using large-scale, discrete element method simulations. Specifically, the stress response due to the application of a localised force perturbation - the Green function technique - is obtained in granular packings generated over several orders of magnitude in both the particle friction coefficient and the applied forcing. We observe crossover behaviour in the mechanical state of the system characterised by the changing nature of the resulting stress response. The transition between anisotropic and isotropic stress response exhibits critical-like features through the identification of a diverging length scale that distinguishes the spatial extent of anisotropic regions from those that display isotropic behaviour. A multidimensional phase diagram is constructed that parameterises the response of the system due to changing friction and force perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Diffusion in jammed particle packs

Using random walk simulations we explore diffusive transport through monodisperse sphere packings over a range of packing fractions, $ϕ$, in the vicinity of the jamming transition at $ϕ_{c}$. Various diffusion properties are computed over several orders of magnitude in both time and packing pressure. Two well-separated regimes of normal, "Fickian" diffusion, where the mean squared displacement is linear in time, are observed. The first corresponds to diffusion inside individual spheres, while the latter is the long-time bulk diffusion. The intermediate anomalous diffusion regime and the long-time value of the diffusion coefficient are both shown to be controlled by particle contacts, which in turn depend on proximity to $ϕ_{c}$. The time required to recover normal diffusion $t^*$ scales as $(ϕ-ϕ_c)^{-0.5}$ and the long-time diffusivity $D_\infty \sim (ϕ- ϕ_c)^{0.5}$, or $D_\infty \sim 1/t^*$. It is shown that the distribution of mean first passage times associated with the escape of random walkers between neighboring particles controls both $t^*$ and $D_\infty$ in the limit $ϕ\rightarrow ϕ_c$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Jamming of frictional spheres and random loose packing

The role of friction coefficient, $μ$, on the jamming properties of disordered, particle packings is studied using computer simulations. Compressed, soft-sphere packings are brought towards the jamming transition - the point where a packing loses mechanical stability - by decreasing the packing fraction. The values of the packing fraction at the jamming transition, $ϕ^μ_{c}$, gradually decrease from the random close packing point for zero friction, to a value coincident with random loose packing as the friction coefficient is increased over several orders of magnitude. This is accompanied by a decrease in the coordination number at the jamming transition, $z^μ_{c}$, which varies from approximately six to four with increasing friction. Universal power law scaling is observed in the pressure and coordination number as a function of distance from the generalised, friction-dependent jamming point. Various power laws are also reported between the $ϕ^μ_{\rm c}$, $z^μ_{\rm c}$, and $μ$. Dependence on preparation history of the packings is also investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stress Response in Confined Arrays of Frictional and Frictionless Particles

Stress transmission inside three dimensional granular packings is investigated using computer simulations. Localized force perturbation techniques are implemented for frictionless and frictional shallow, ordered, granular arrays confined by solid boundaries for a range of system sizes. Stress response profiles for frictional packings agree well with the predictions for the semi-infinite half plane of classical isotropic elasticity theory down to boxes of linear dimensions of approximately forty particle diameters and over several orders of magnitude in the applied force. The response profiles for frictionless packings exhibit a transitional regime to strongly anisotropic features with increasing box size. The differences between the nature of the stress response are shown to be characterized by very different particle displacement fields.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Tuning Jammed Frictionless Disk Packings from Isostatic to Hyperstatic

We perform extensive computational studies of two-dimensional static bidisperse disk packings using two distinct packing-generation protocols. The first involves thermally quenching equilibrated liquid configurations to zero temperature over a range of thermal quench rates $r$ and initial packing fractions followed by compression and decompression in small steps to reach packing fractions $ϕ_J$ at jamming onset. For the second, we seed the system with initial configurations that promote micro- and macrophase-separated packings followed by compression and decompression to $ϕ_J$. We find that amorphous, isostatic packings exist over a finite range of packing fractions from $ϕ_{\rm min} \le ϕ_J \le ϕ_{\rm max}$ in the large-system limit, with $ϕ_{\rm max} \approx 0.853$. In agreement with previous calculations, we obtain $ϕ_{\rm min} \approx 0.84$ for $r > r^*$, where $r^*$ is the rate above which $ϕ_J$ is insensitive to rate. We further compare the structural and mechanical properties of isostatic versus hyperstatic packings. The structural characterizations include the contact number, bond orientational order, and mixing ratios of the large and small particles. We find that the isostatic packings are positionally and compositionally disordered, whereas bond-orientational and compositional order increase with contact number for hyperstatic packings. In addition, we calculate the static shear modulus and normal mode frequencies of the static packings to understand the extent to which the mechanical properties of amorphous, isostatic packings are different from partially ordered packings. We find that the mechanical properties of the packings change continuously as the contact number increases from isostatic to hyperstatic.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Long wavelength structural anomalies in jammed systems

The structural properties of static, jammed packings of monodisperse spheres in the vicinity of the jamming transition are investigated using large-scale computer simulations. At small wavenumber $k$, we argue that the anomalous behavior in the static structure factor, $S(k) \sim k$, is consequential of an excess of low-frequency, collective excitations seen in the vibrational spectrum. This anomalous feature becomes more pronounced closest to the jamming transition, such that $S(0) \to 0$ at the transition point. We introduce an appropriate dispersion relation that accounts for these phenomena that leads us to relate these structural features to characteristic length scales associated with the low-frequency vibrational modes of these systems. When the particles are frictional, this anomalous behavior is suppressed providing yet more evidence that jamming transitions of frictional spheres lie at lower packing fractions that that for frictionless spheres. These results suggest that the mechanical properties of jammed and glassy media may therefore be inferred from measurements of both the static and dynamical structure factors.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Stress response inside perturbed particle assemblies

The effect of structural disorder on the stress response inside three dimensional particle assemblies is studied using computer simulations of frictionless sphere packings. Upon applying a localised, perturbative force within the packings, the resulting {\it Green's} function response is mapped inside the different assemblies, thus providing an explicit view as to how the imposed perturbation is transmitted through the packing. In weakly disordered arrays, the resulting transmission of forces is of the double-peak variety, but with peak widths scaling linearly with distance from the source of the perturbation. This behaviour is consistent with an anisotropic elasticity response profile. Increasing the disorder distorts the response function until a single-peak response is obtained for fully disordered packings consistent with an isotropic description.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Normal Modes in Model Jammed Systems in Three Dimensions

Vibrational spectra and normal modes of mechanically stable particle packings in three dimensions are analyzed over a range of compressions, from near the jamming transition, where the packings lose their rigidity, to far above it. At high frequency, the normal modes are localized at all compressions. At low frequency, the nature of the modes depends somewhat on compression. At large compressions, far from the transition, the lowest-frequency normal modes have some plane-wave character, though less than one would expect for a crystalline or isotropic solid. At low compressions near the jamming transition, the lowest-frequency modes are neither plane-wave-like nor localized. We characterize these differences, highlighting the unusual dispersion behavior that emerges for marginally jammed solids.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Force heterogeneities in particle assemblies: From order to disorder

The effect of increasing structural disorder on the distribution of contact forces P(f), inside three dimensional particle assemblies is systematically studied using computer simulations of model granular packings. Starting from a face-centred cubic array, where all contact forces are identical, an increasing number of defects is introduced into the assembly, after which the system is then allowed to relax into a new mechanically stable state. Three distinct protocols for imposing disorder are compared. A quantitative measure of the disorder is obtained from distributions of the coordination number and three-particle contact angle. The distribution of normal contact forces show dramatic qualitative changes with increasing disorder. In the regime where the disorder is relatively weak, the pressure and the lowest normal mode frequency scale approximately linearly in the coordination number, with distance from the crystalline state. These results for P(f) are discussed in the context of jamming phenomena in glassy and granular materials.

0 citations0 reviewsNo code linkedOpen access