Source author record

Jeremy B. Lechman

Jeremy B. Lechman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

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.

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.

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.

preprint2020arXiv

Heterogeneous partition of cellular blood-borne nanoparticles through microvascular bifurcations

Blood flowing through microvascular bifurcations has been an active research topic for many decades, while the partitioning pattern of nanoscale solutes in the blood remains relatively unexplored. Here, we demonstrate a multiscale computational framework for direct numerical simulation of the nanoparticle (NP) partitioning through physiologically-relevant vascular bifurcations in the presence of red blood cells (RBCs). The computational framework is established by embedding a newly-developed particulate suspension inflow/outflow boundary condition into a multiscale blood flow solver. The computational framework is verified by recovering a tubular blood flow without a bifurcation and validated against the experimental measurement of an intravital bifurcation flow. The classic Zweifach-Fung (ZF) effect is shown to be well captured by the method. Moreover, we observe that NPs exhibit a ZF-like heterogeneous partition in response to the heterogeneous partition of the RBC phase. The NP partitioning prioritizes the high-flow-rate daughter branch except for extreme (large or small) suspension flow partition ratios under which the complete phase separation tends to occur. By analyzing the flow field and the particle trajectories, we show that the ZF-like heterogeneity in NP partition can be explained by the RBC-entrainment effect caused by the deviation of the flow separatrix preceded by the tank-treading of RBCs near the bifurcation junction. The recovery of homogeneity in the NP partition under extreme flow partition ratios is due to the plasma skimming of NPs in the cell-free layer. These findings, based on the multiscale computational framework, provide biophysical insights to the heterogeneous distribution of NPs in microvascular beds that are observed pathophysiologically.

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

preprint2011arXiv

Evaporation of Lennard-Jones Fluids

Evaporation and condensation at a liquid/vapor interface are ubiquitous interphase mass and energy transfer phenomena that are still not well understood. We have carried out large scale molecular dynamics simulations of Lennard-Jones (LJ) fluids composed of monomers, dimers, or trimers to investigate these processes with molecular detail. For LJ monomers in contact with a vacuum, the evaporation rate is found to be very high with significant evaporative cooling and an accompanying density gradient in the liquid domain near the liquid/vapor interface. Increasing the chain length to just dimers significantly reduces the evaporation rate. We confirm that mechanical equilibrium plays a key role in determining the evaporation rate and the density and temperature profiles across the liquid/vapor interface. The velocity distributions of evaporated molecules and the evaporation and condensation coefficients are measured and compared to the predictions of an existing model based on kinetic theory of gases. Our results indicate that for both monatomic and polyatomic molecules, the evaporation and condensation coefficients are equal when systems are not far from equilibrium and smaller than one, and decrease with increasing temperature. For the same reduced temperature $T/T_c$, where $T_c$ is the critical temperature, these two coefficients are higher for LJ dimers and trimers than for monomers, in contrast to the traditional viewpoint that they are close to unity for monatomic molecules and decrease for polyatomic molecules. Furthermore, data for the two coefficients collapse onto a master curve when plotted against a translational length ratio between the liquid and vapor phase.

preprint2005arXiv

Three-dimensional shear in granular flow

The evolution of granular shear flow is investigated as a function of height in a split-bottom Couette cell. Using particle tracking, magnetic-resonance imaging, and large-scale simulations we find a transition in the nature of the shear as a characteristic height $H^*$ is exceeded. Below $H^*$ there is a central stationary core; above $H^*$ we observe the onset of additional axial shear associated with torsional failure. Radial and axial shear profiles are qualitatively different: the radial extent is wide and increases with height while the axial width remains narrow and fixed.