Researcher profile

Gary S. Grest

Gary S. Grest contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

preprint2026arXiv

Shape Elasticity in Colloidal Bent-Core Liquid Crystals

Curved particles have been shown to stabilize a range of states with unique order in dense suspensions of colloidal bent core liquid crystals. The shape of the colloidal rods encourages the formation of curved director fields. However, states of constant bend cannot uniformly fill either two or three dimensional Euclidean space and are therefore geometrically frustrated. As a result, curved rods are forced to couple their preference for bend with additional twist and splay deformations, giving rise to twist-bend and splay-bend states of nematic and smectic order. In this article, we study the effect of rod curvature on these diverse states of liquid crystalline order using molecular dynamics simulations of a bonded particle model of curved rods with tunable shape elasticity. Focusing on the case of intermediately curved rods, we find that curved rods go through a sequence of isotropic, nematic twist-bend and smectic splay-bend ordering as the density is increased from the dilute limit, in agreement with previous studies of rigid rods. As the rods become more elastic, the critical concentration separating these phases is shifted to higher density. Lastly, we find that flexibility weakens the first-order phase transition separating the isotropic and nematic twist-bend phases.

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.

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.

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.

preprint2022arXiv

Shear Is Not Always Simple: Rate-Dependent Effects of Flow Type on Granular Rheology

Despite there being an infinite variety of types of flow, most rheological studies focus on a single type such as simple shear. Using discrete element simulations, we explore bulk granular systems in a wide range of flow types at large strains and characterize invariants of the stress tensor for different inertial numbers and interparticle friction coefficients. We identify a strong dependence on the type of flow which grows with increasing inertial number or friction. Standard models of yielding, repurposed to describe the dependence of the stress on flow type in steady-state flow and at finite rates, are compared with data.

preprint2022arXiv

Super-stretchable Elastomer from Cross-linked Ring Polymers

The stretchability of polymeric materials is critical to many applications such as stretchable electronics and soft robotics, yet the stretchability of conventional cross-linked linear polymers is limited by the entanglements between polymer chains. We show using molecular dynamics simulations that cross-linked ring polymers are significantly more stretchable than cross-linked linear polymers. Compared to linear polymers, the entanglements between ring polymers do not act as effective cross-links. As a result, the stretchability of cross-linked ring polymers is determined by the maximum extension of polymer strands between cross-links, rather than between entanglements as in cross-linked linear polymers. The more compact conformation of ring polymers before deformation also contributes to the increase in stretchability.

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.

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.

preprint2019arXiv

Topological Linking Drives Anomalous Thickening of Ring Polymers In Weak Extensional Flows

Molecular dynamics simulations confirm recent extensional flow experiments showing ring polymer melts exhibit strong extension-rate thickening of the viscosity at Weissenberg numbers $Wi<<1$. Thickening coincides with the extreme elongation of a minority population of rings that grows with $Wi$. The large susceptibility of some rings to extend is due to a flow-driven formation of topological links that connect multiple rings into supramolecular chains. Links form spontaneously with a longer delay at lower $Wi$ and are pulled tight and stabilized by the flow. Once linked, these composite objects experience larger drag forces than individual rings, driving their strong elongation. The fraction of linked rings generated by flow depends non-monotonically on $Wi$, increasing to a maximum when $Wi\sim1$ before rapidly decreasing when the strain rate approaches the relaxation rate of the smallest ring loops $\sim 1/τ_e$.

preprint2010arXiv

End Grafted Polymer Nanoparticles in a Polymeric Matrix: Effect of Coverage and Curvature

It has recently been proposed that the miscibility of nanoparticles with a polymer matrix can be controlled by grafting polymer chains to the nanoparticle surface. As a first step to study this situation, we have used molecular dynamics simulations on a single nanoparticle of radius R ($4σ\le$R$\le 16σ$, where $σ$ is the diameter of a polymer monomer) grafted with chains of length 500 in a polymer melt of chains of length 1000. The grafting density $Σ$ was varied between $0.04$-$0.32$ chains/$σ^2$. To facilitate equilibration a Monte Carlo double-bridging algorithm is applied - new bonds are formed across a pair of chains, creating two new chains each substantially different from the original. For the long brush chains studied here, the structure of the brush assumes its large particle limit even for $R$ as small as 8$σ$, which is consistent with recent experimental findings. We study autophobic dewetting of the melt from the brush as a function of increasing $Σ$. Even these long brush and matrix chains of length $6$ and $12$ $N_e$, respectively, (the entanglement length is $N_e \sim 85$) give somewhat ambiguous results for the interfacial width, showing that studies of two or more nanoparticles are necessary to properly understand these miscibility issues. Entanglement between the brush and melt chains were identified using the primitive path analysis. We find that the number of entanglements between the brush and melt chains scale simply with the product of the local monomer densities of brush and melt chains.