Source author record

C. S. O'Hern

C. S. O'Hern 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

6works
3topics
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

6 published item(s)

preprint2022arXiv

Hopper flows of deformable particles

Numerous experimental and computational studies show that continuous hopper flows of granular materials obey the Beverloo equation that relates the volume flow rate $Q$ and the orifice width $w$: $Q \sim (w/σ_{\rm avg}-k)^β$, where $σ_{\rm avg}$ is the average particle diameter, $kσ_{\rm avg}$ is an offset where $Q\sim 0$, the power-law scaling exponent $β=d-1/2$, and $d$ is the spatial dimension. Recent studies of hopper flows of deformable particles in different background fluids suggest that the particle stiffness and dissipation mechanism can also strongly affect the power-law scaling exponent $β$. We carry out computational studies of hopper flows of deformable particles with both kinetic friction and background fluid dissipation in two and three dimensions. We show that the exponent $β$ varies continuously with the ratio of the viscous drag to the kinetic friction coefficient, $λ=ζ/μ$. $β= d-1/2$ in the $λ\rightarrow 0$ limit and $d-3/2$ in the $λ\rightarrow \infty$ limit, with a midpoint $λ_c$ that depends on the hopper opening angle $θ_w$. We also characterize the spatial structure of the flows and associate changes in spatial structure of the hopper flows to changes in the exponent $β$. The offset $k$ increases with particle stiffness until $k \sim k_{\rm max}$ in the hard-particle limit, where $k_{\rm max} \sim 3.5$ is larger for $λ\rightarrow \infty$ compared to that for $λ\rightarrow 0$. Finally, we show that the simulations of hopper flows of deformable particles in the $λ\rightarrow \infty$ limit recapitulate the experimental results for quasi-2D hopper flows of oil droplets in water.

preprint2016arXiv

Protocol dependence of the jamming transition

We propose a theoretical framework for predicting the protocol dependence of the jamming transition for frictionless spherical particles that interact via purely repulsive contact forces. We study isostatic jammed disk packings obtained via two protocols: isotropic compression and simple shear. We show that for frictionless systems, all jammed packings can be obtained via either protocol. However, the probability to obtain a particular jammed packing depends on the packing-generation protocol. We predict the average shear strain required to induce jamming in initially unjammed packings from the measured probability to jam at packing fraction $ϕ$ from isotropic compression. We compare our predictions to results from numerical simulations of jamming and find quantitative agreement. We also show that the packing fraction range, over which strain-induced jamming occurs, tends to zero in the large system limit for frictionless packings with overdamped dynamics.

preprint2015arXiv

Experiments demonstrate that the null space of the rigidity matrix determines grain motion during vibration-induced compaction

Using a previously developed experimental method to reduce friction in mechanically stable packings of disks, we find that frictional packings form tree-like structures of geometrical families that lie on reduced dimensional manifolds in configuration space. Each branch of the tree begins at a point in configuration space with an isostatic number of contacts and spreads out to sequentially higher dimensional manifolds as the number of contacts are reduced. We find that gravitational deposition of disks produces an initially under-coordinated packing stabilized by friction on a high-dimensional manifold. Using short vibration bursts to reduce friction, we compact the system through many stable configurations with increasing contact number and decreasing dimensionality until the system reaches an isostatic frictionless state. We find that this progression can be understood as the system moving through the null-space of the rigidity matrix defined by the interparticle contact network in the direction of the gravitational force. We suggest that this formalism can also be used to explain the evolution of frictional packings under other forcing conditions.

preprint2015arXiv

On the origin of multi-component bulk metallic glasses: Atomic size mismatches and de-mixing

The critical cooling rate $\mathcal{R}_c$, below which liquids crystallize upon cooling, characterizes the glass-forming ability (GFA) of the system. While pure metals are typically poor glass formers with $\mathcal {R}_c>10^{12}\, {\rm K/s}$, specific multi-component alloys can form bulk metallic glasses (BMGs) even at cooling rates below $\mathcal {R}\sim 1\, {\rm K/s}$. Conventional wisdom asserts that metal alloys with three or more components are better glass formers (with smaller ${\cal R}_c$) than binary alloys. However, there is currently no theoretical framework that provides quantitative predictions for $\mathcal{R}_c$ for multi-component alloys. We perform simulations of ternary hard-sphere systems, which have been shown to be accurate models for the glass-forming ability of BMGs, to understand the roles of geometric frustration and demixing in determining $\mathcal {R}_c$. Specifically, we compress ternary hard sphere mixtures into jammed packings and measure the critical compression rate, below which the system crystallizes, as a function of the diameter ratios $σ_B/σ_A$ and $σ_C/σ_A$ and number fractions $x_A$, $x_B$, and $x_C$. We find two distinct regimes for the GFA in parameter space for ternary hard spheres. When the diameter ratios are close to $1$, such that the largest ($A$) and smallest ($C$) species are well-mixed, the GFA of ternary systems is no better than that of the optimal binary glass former. However, when $σ_C/σ_A \lesssim 0.8$ is below the demixing threshold for binary systems, adding a third component $B$ with $σ_C < σ_B < σ_A$ increases the GFA of the system by preventing demixing of $A$ and $C$. Analysis of the available data from experimental studies indicates that most ternary BMGs are below the binary demixing threshold with $σ_C/σ_A < 0.8$.

preprint2011arXiv

Calculations of the Structure of Basin Volumes for Mechanically Stable Packings

There are a finite number of distinct mechanically stable (MS) packings in model granular systems composed of frictionless spherical grains. For typical packing-generation protocols employed in experimental and numerical studies, the probabilities with which the MS packings occur are highly nonuniform and depend strongly on parameters in the protocol. Despite intense work, it is extremely difficult to predict {\it a priori} the MS packing probabilities, or even which MS packings will be the most versus the least probable. We describe a novel computational method for calculating the MS packing probabilities by directly measuring the volume of the MS packing `basin of attraction', which we define as the collection of initial points in configuration space at {\it zero packing fraction} that map to a given MS packing by following a particular dynamics in the density landscape. We show that there is a small core region with volume $V^c_n$ surrounding each MS packing $n$ in configuration space in which all initial conditions map to a given MS packing. However, we find that the MS packing probabilities are very weakly correlated with core volumes. Instead, MS packing probabilities obtained using initially dilute configurations are determined by complex geometric features of the basin of attraction that are distant from the MS packing.

preprint1999arXiv

Sliding Phases in XY-Models, Crystals, and Cationic Lipid-DNA Complexes

We predict the existence of a totally new class of phases in weakly coupled, three-dimensional stacks of two-dimensional (2D) XY-models. These ``sliding phases'' behave essentially like decoupled, independent 2D XY-models with precisely zero free energy cost associated with rotating spins in one layer relative to those in neighboring layers. As a result, the two-point spin correlation function decays algebraically with in-plane separation. Our results, which contradict past studies because we include higher-gradient couplings between layers, also apply to crystals and may explain recently observed behavior in cationic lipid-DNA complexes.