Source author record

Karen E. Daniels

Karen E. Daniels 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

37works
14topics
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

37 published item(s)

preprint2022arXiv

Controlling rheology via boundary conditions in dense granular flows

Boundary shape, particularly roughness, strongly controls the amount of wall slip in dense granular flows. In this paper, we aim to quantify and understand which aspects of a dense granular flow are controlled by the boundary condition, and to incorporate these observations into a cooperative nonlocal model characterizing slow granular flows. To examine the influence of boundary properties, we perform experiments on a quasi-2D annular shear cell with a rotating inner wall and a fixed outer wall; the later is selected from among 6 walls with various roughness, local concavity, and compliance. We find that we can successfully capture the full flow profile using a single set of empirically determined model parameters, with only the wall slip velocity set by direct observation. Through the use of photoelastic particles, we observe how the internal stresses fluctuate more for rougher boundaries, corresponding to lower wall slip, and connect this observation to the propagation of nonlocal effects originating at the wall.

preprint2021arXiv

Betweenness centrality illuminates intermittent frictional dynamics

Dense granular systems subjected to an imposed shear stress undergo stick-slip dynamics with systematic patterns of dilation-compaction. During each stick phase, as the frictional strength builds up, the granular system dilates to accommodate shear strain, developing stronger force networks. During each slip event, when the stored energy is released, particles experience large rearrangements and the granular network can significantly change. Here, we use numerical simulations of 3D, sheared frictional packings to show that the mean betweenness centrality -- a property of network of interparticle connections -- follows consistent patterns during the stick-slip dynamics, showing sharp spikes at each slip event. We identify the source of this behavior as arising from the connectivity and contact arrangements of granular network during dilation-compaction cycles, and find that a lower potential for connection between particles leads to an increase of mean betweenness centrality in the system. Furthermore, we show that at high confinements, few particles lose contact during slip events, leading to a smaller change in granular connectivity and betweenness centrality.

preprint2021arXiv

Interfacial Tension Modulation of Liquid Metal via Electrochemical Oxidation

This progress report summarizes recent studies of electrochemical oxidation to modulate the interfacial tension of gallium-based alloys. These alloys, which are liquid at ambient conditions, have the largest interfacial tension of any liquid at room temperature. The ability to modulate the tension offers the possibility to create forces that change the shape and position of the metal. It has been known since the late 1800s that electrocapillarity-the use of potential to modulate the electric double layer on the surface of metals in electrolyte-lowers the interfacial tension of liquid metal. Yet, this phenomenon can only achieve modest changes in interfacial tension since it is limited to potential windows that avoid reactions. A recent discovery suggests that reactions driven by the electrochemical oxidation of gallium alloys cause the interfacial tension to decrease from ~500 mN/m at 0 V to ~0 mN/m at ~0.8 V, a change in tension that goes well beyond what is possible via conventional electrocapillarity or surfactants. The changes in tension are reversible; reductive potentials return the metal back to a state of high interfacial tension. This report aims to summarize key work and introduce beginners to this field by including electrochemistry basics while addressing misconceptions. We discuss applications that utilize modulations in interfacial tension of liquid metal and conclude with remaining opportunities and challenges that need further investigation.

preprint2020arXiv

Sponge-like rigid structures in frictional granular packings

We show how rigidity emerges in experiments of sheared frictional granular materials by using generalizations of two methods for identifying rigid structures. Both approaches, the force-based dynamical matrix and the topology-based rigidity percolation, agree with each other and identify similar rigid structures. As the system becomes jammed, at a contact number of $z=2.4\pm 0.1$, a rigid backbone interspersed with floppy, particle-filled holes of a broad range of sizes emerges, creating a sponge-like morphology. We also find that the pressure within rigid structures always exceeds the pressure outside the rigid structures, i.e. that the backbone is load-bearing. These findings shows that it is necessary to go beyond mean-field theory to capture the physics of frictional jamming and also suggests that mechanical stability arises through arch structures and hinges at the mesoscale.

preprint2019arXiv

Enlightening force chains: a review of photoelasticimetry in granular matter

A photoelastic material will reveal its internal stresses when observed through polarizing filters. This eye-catching property has enlightened our understanding of granular materials for over half a century, whether in the service of art, education, or scientific research. In this review article in honor of Robert Behringer, we highlight both his pioneering use of the method in physics research, and its reach into the public sphere through museum exhibits and outreach programs. We aim to provide clear protocols for artists, exhibit-designers, educators, and scientists to use in their own endeavors. It is our hope that this will build awareness about the ubiquitous presence of granular matter in our lives, enlighten its puzzling behavior, and promote conversations about its importance in environmental and industrial contexts. To aid in this endeavor, this paper also serves as a front door to a detailed wiki containing open, community-curated guidance on putting these methods into practice.

preprint2016arXiv

Capillary fracture of ultrasoft gels: heterogeneity and delayed nucleation

A droplet of surfactant spreading on an ultrasoft ($E \lesssim 100$ Pa) gel substrate will produce capillary fractures at the gel surface; these fractures originate at the contact-line and propagate outwards in a starburst pattern. There is an inherent variability in both the number of fractures formed and the time delay before fractures form. In the regime where single fractures form, we observe a Weibull-like distribution of delay times, consistent with a thermally-activated process. The shape parameter is close to 1 for softer gels (a Poisson process), and larger for stiffer gels (indicative of aging). For single fractures, the characteristic delay time is primarily set by the elastocapillary length of the system, calculated from the differential in surface tension between the droplet and the substrate, rather than the elastic modulus as for stiffer systems. For multiple fractures, all fractures appear simultaneously and long delay times are suppressed. The delay time distribution provides a new technique for probing the energy landscape and fracture toughness of ultrasoft materials.

preprint2016arXiv

Deformation of an Elastic Substrate Due to a Resting Sessile Droplet

On a sufficiently-soft substrate, a resting fluid droplet will cause significant deformation of the substrate. This deformation is driven by a combination of capillary forces at the contact line and the fluid pressure at the solid surface. These forces are balanced at the surface by the solid traction stress induced by the substrate deformation. Young's Law, which predicts the equilibrium contact angle of the droplet, also indicates an a priori radial force balance for rigid substrates, but not necessarily for soft substrates which deform under loading. It remains an open question whether the contact line transmits a non-zero radial force to the substrate surface in addition to the conventional vertical force. We present an analytic Fourier transform solution technique that includes general interfacial energy conditions which govern the contact angle of a 2D droplet. This includes evaluating the effect of gravity on the droplet shape in order to determine the correct fluid pressure at the substrate surface for larger droplets. Importantly, we find that in order to avoid a strain singularity at the contact line under a nonzero radial contact line force, it is necessary to include a previously-neglected radial traction boundary condition. To quantify the effects of the contact line and identify key quantities that will be experimentally-accessible for testing the model, we evaluate solutions for the substrate surface displacement field as a function of Poisson's ratio and zero/non-zero radial contact line forces.

preprint2016arXiv

Evolution of network architecture in a granular material under compression

As a granular material is compressed, the particles and forces within the system arrange to form complex heterogeneous structures. Force chains are a prime example and are thought to constrain bulk properties such as mechanical stability and acoustic transmission. However, characterizing the dynamic nature of mesoscale architectures in granular systems can be challenging. A growing body of work has shown that graph theoretic approaches may provide a useful foundation for tackling these problems. Here, we extend current approaches by utilizing multilayer networks as a framework for directly quantifying the evolution of mesoscale architecture in a compressed granular system. We examine a quasi-two-dimensional aggregate of photoelastic disks, subject to biaxial compression through a series of small, quasistatic steps. Treating particles as network nodes and inter-particle forces as network edges, we construct a multilayer network by linking together the series of static force networks that exist at each strain step. We then extract the inherent mesoscale structure from the system by using a generalization of community detection methods, and we define quantitative measures to characterize the reconfiguration and evolution of this structure throughout compression. By separately considering the network of normal and tangential forces, we find that they display different structural evolution. To test the sensitivity of the network model to particle properties, we examine whether the method can distinguish a subsystem of low-friction particles within a bath of higher-friction particles. We find that this can be done by considering the network of tangential forces. The results discussed throughout this study suggest that these novel network science techniques may provide a direct way to compare and classify data from systems under different external conditions or with different physical makeup.

preprint2016arXiv

High refractive index immersion liquid for super-resolution 3D imaging using sapphire-based aNAIL optics

Optically-transparent immersion liquids with refractive index (n ~ 1.77) to match sapphire-based aplanatic numerical aperture increasing lens (aNAIL) are necessary for achieving deep 3D imaging with high spatial resolution. We report that antimony tribromide (SbBr$_{3}$) salt dissolved in liquid diiodomethane (CH$_{2}$I$_{2}$) provides a new high refractive index immersion liquid for optics applications. The refractive index is tunable from n = 1.74 (pure) to n = 1.873 (saturated), by adjusting either salt concentration or temperature; this allows it to match (or even exceed) the refractive index of sapphire. Importantly, the solution gives excellent light transmittance in the ultraviolet to near-infrared range, an improvement over commercially-available immersion liquids. This refractive index matched immersion liquid formulation has enabled us to develop a sapphire-based aNAIL objective that has both high numerical aperture (NA = 1.17) and long working distance (WD = 12 mm). This opens up new possibilities for deep 3D imaging with high spatial resolution.

preprint2016arXiv

Nonaffine deformation under compression and decompression of a flow-stabilized solid

Understanding the particle-scale transition from elastic deformation to plastic flow is central to making predictions about the bulk material properties and response of disordered materials. To address this issue, we perform experiments on flow-stabilized solids composed of micron-scale spheres within a microfluidic channel, in a regime where particle inertia is negligible. Each solid heap exists within a stress gradient imposed by the flow, and we track the positions of particles in response to single impulses of fluid-driven compression or decompression. We find that the resulting deformation field is well-decomposed into an affine field, with a constant strain profile throughout the solid, and a non-affine field. The magnitude of this non-affine response decays with the distance from the free surface in the long-time limit, suggesting that the distance from jamming plays a significant role in controlling the length scale of plastic flow. Finally, we observe that compressive pulses create more rearrangements than decompressive pulses, an effect that we quantify using the $D^2_\mathrm{min}$ statistic for non-affine motion. Unexpectedly, the time scale for the compression response is shorter than for decompression at the same strain (but unequal pressure), providing insight into the coupling between deformation and cage-breaking.

preprint2016arXiv

Simulating Surfactant Spreading: Impact of a Physically Motivated Equation of State

For more than two decades, a single model for the spreading of a surfactant-driven thin liquid film has dominated the applied mathematics literature on the subject. Recently, through the use of fluorescently-tagged lipids, it has become possible to make direct, quantitative comparisons between experiments and models. These comparisons have revealed two important discrepancies between simulations and experiments: the spatial distribution of the surfactant layer, and the timescale over which spreading occurs. In this paper, we present numerical simulations that demonstrate the impact of the particular choice of the equation of state (EoS) relating the surfactant concentration to the surface tension. Previous choices of the model EoS have been an ad-hoc decreasing function. Here, we instead propose an empirically-motivated equation of state; this provides a route to resolving some discrepancies and raises new issues to be pursued in future experiments. In addition, we test the influence of the choice of initial conditions and values for the non-dimensional groups. We demonstrate that the choice of EoS improves the agreement in surfactant distribution morphology between simulations and experiments, and impacts the dynamics of the simulations. The relevant feature of the EoS, the gradient, has distinct regions for empirically motivated choices, which suggests that future work will need to consider more than one timescale. We observe that the non-dimensional number controlling the relative importance of gravitational vs. capillary forces has a larger impact on the dynamics than the other non-dimensional groups. Finally, we observe that the experimental approach of using a ring to contain the surfactant could affect the surfactant and fluid dynamics if it disrupts the intended initial surfactant distribution. However, the meniscus itself does not significantly affect the dynamics.

preprint2016arXiv

Spatio-temporal Patterns in Inclined Layer Convection

This paper reports on a theoretical analysis of the rich variety of spatio-temporal patterns observed recently in inclined layer convection at medium Prandtl number when varying the inclination angle $γ$ and the Rayleigh number $R$. The present numerical investigation of the inclined layer convection system is based on the standard Oberbeck-Boussinesq equations. The patterns are shown to originate from a complicated competition of buoyancy-driven and shear-flow driven pattern forming mechanisms. The former are expressed as \rm{longitudinal} convection rolls with their axes oriented parallel to the incline, the latter as perpendicular \rm{transverse} rolls. Along with conventional methods to study roll patterns and their stability, we employ direct numerical simulations in large spatial domains, comparable with the experimental ones. As a result, we determine the phase diagram of the characteristic complex 3D convection patterns above onset of convection in the $γ-R$ plane, and find that it compares very well with the experiments. In particular we demonstrate that interactions of specific Fourier modes, characterized by a resonant interaction of their wavevectors in the layer plane, are key to understanding the pattern morphologies.

preprint2015arXiv

Extraction of Force-Chain Network Architecture in Granular Materials Using Community Detection

Force chains form heterogeneous physical structures that can constrain the mechanical stability and acoustic transmission of granular media. However, despite their relevance for predicting bulk properties of materials, there is no agreement on a quantitative description of force chains. Consequently, it is difficult to compare the force-chain structures in different materials or experimental conditions. To address this challenge, we treat granular materials as spatially-embedded networks in which the nodes (particles) are connected by weighted edges that represent contact forces. We use techniques from community detection, which is a type of clustering, to find sets of closely connected particles. By using a geographical null model that is constrained by the particles' contact network, we extract chain-like structures that are reminiscent of force chains. We propose three diagnostics to measure these chain-like structures, and we demonstrate the utility of these diagnostics for identifying and characterizing classes of force-chain network architectures in various materials. To illustrate our methods, we describe how force-chain architecture depends on pressure for two very different types of packings: (1) ones derived from laboratory experiments and (2) ones derived from idealized, numerically-generated frictionless packings. By resolving individual force chains, we quantify statistical properties of force-chain shape and strength, which are potentially crucial diagnostics of bulk properties (including material stability). These methods facilitate quantitative comparisons between different particulate systems, regardless of whether they are measured experimentally or numerically.

preprint2015arXiv

Friction and Pressure-Dependence of Force Chain Communities in Granular Materials

Granular materials transmit stress via a network of force chains. Despite the importance of these chains to characterizing the stress state and dynamics of the system, there is no common framework for quantifying their their properties. Recently, attention has turned to the tools of network science as a promising route to such a description. In this paper, we apply community detection techniques to numerically-generated packings of spheres over a range of interparticle friction coefficients and confining pressures. In order to extract chain-like features, we use a modularity maximization with a recently-developed geographical null model \cite{Bassett2015}, and optimize the technique to detect branched structures by minimizing the normalized convex hull of the detected communities. We characterize the force chain communities by their size (number of particles), network strength (internal forces), and normalized convex hull ratio (sparseness). We find the that the first two exhibit an approximately linear correlation and are therefore largely redundant. For both pressure $P$ and interparticle friction $μ$, we observe crossovers in behavior. For $μ\lesssim 0.1$, the packings exhibit more sensitivity to pressure. In addition, we identify a crossover pressure where the frictional dependence switches from having more large/strong communities at low $μ$ vs. high $μ$. We explain these phenomena by comparison to the spatial distribution of communities along the vertical axis of the system. These results provide new tools for considering the mesoscale structure of a granular system and pave the way for reduced descriptions based on the force chain structure.

preprint2015arXiv

Stretching Rubber, Stretching Minds: a polymer physics lab for teaching entropy

Entropy is a difficult concept to teach using real-world examples. Unlike temperature, pressure, volume, or work, it's not a quantity which most students encounter in their day-to-day lives. Even the way entropy is often qualitatively described, as a measure of disorder, is incomplete and can be misleading. In an effort to address these obstacles, we have developed a simple laboratory activity, the stretching of an elastic rubber sheet, intended to give students hands-on experience with the concepts of entropy, temperature and work in both adiabatic and quasistatic processes. We present two versions of the apparatus: a double-lever system, which may be reproduced with relatively little cost, and a commercial materials testing system, which provides students experience with scientific instrumentation that is used in research.

preprint2014arXiv

Elastocapillary deformations on partially-wetting substrates: rival contact-line models

A partially-wetting liquid can deform the underlying elastic substrate upon which it rests. This situation requires the development of theoretical models to describe the wetting forces imparted by the drop onto the solid substrate, particularly those at the contact-line. We construct a general solution using a displacement potential function for the elastic deformations within a finite elastic substrate associated with these wetting forces, and compare the results for several different contact-line models. Our work incorporates internal contributions to the surface stress from both liquid/solid $Σ_{ls}$ and solid/gas $Σ_{sg}$ solid surface tensions (surface stress), which results in a non-standard boundary-value problem that we solve using a dual integral equation. We compare our results to relevant experiments and conclude that the generalization of solid surface tension $Σ_{ls} \neq Σ_{sg}$ is an essential feature in any model of partial-wetting. The comparisons also allow us to systematically eliminate some proposed contact-line models.

preprint2014arXiv

Nonlinear Elasticity of Flow-Stabilized Solids

Thermal fluctuations, geometric exclusion, and external driving all govern the mechanical response of dense particulate suspensions. Here, we measure the stress-strain response of quasi-two-dimensional flow-stabilized microsphere heaps in a regime in which all three effects are present using a microfluidic device. We observe that the elastic modulus and the mean interparticle separation of the heaps are tunable via the confining stress provided by the fluid flow. Furthermore, the measured stress-strain curves exhibit a universal nonlinear shape which can be predicted from a thermal van der Waals equation of state with excluded volume. This analysis indicates that many-body interactions contribute a significant fraction of the stress supported by the heap.

preprint2014arXiv

Self-healing Dynamics of Surfactant Coatings on Thin Viscous Films

We investigate the dynamics of an insoluble surfactant on the surface of a thin viscous fluid spreading inward to fill a surfactant-free region. During the initial stages of surfactant self-healing, Marangoni forces drive an axisymmetric ridge inward to coalesce into a growing central distension; this is unlike outward-spreading, in which the ridge decays. In later dynamics, the distension slowly decays and the surfactant concentration equilibrates. We present results from experiments in which we simultaneously measure the surfactant concentration (using fluorescently-tagged lipids) and the fluid height profile (via laser profilometry). We compare the results to simulations of a mathematical model using parameters from our experiments. For surfactant concentrations close to but below the critical monolayer concentration, we observe agreement between the height profiles in the numerical simulations and the experiment, but disagreement in the surfactant distribution. In experiments at lower concentrations, the surfactant spreading and formation of a Marangoni ridge are no longer present, and a persistent lipid-free region remains. This observation, which is not captured by the simulations, has undesirable implications for applications where uniform coverage is advantageous. Finally, we probe the generality of the effect, and find that distensions of similar size are produced independent of initial fluid thickness, size of initial clean region, and surfactant type.

preprint2014arXiv

The Statistical Physics of Athermal Materials

At the core of equilibrium statistical mechanics lies the notion of statistical ensembles: a collection of microstates, each occurring with a given a priori probability that depends only on a few macroscopic parameters such as temperature, pressure, volume, and energy. In this review article, we discuss recent advances in establishing statistical ensembles for athermal materials. The broad class of granular and particulate materials is immune from the effects of thermal fluctuations because the constituents are macroscopic. In addition, interactions between grains are frictional and dissipative, which invalidates the fundamental postulates of equilibrium statistical mechanics. However, granular materials exhibit distributions of microscopic quantities that are reproducible and often depend on only a few macroscopic parameters. We explore the history of statistical ensemble ideas in the context of granular materials, clarify the nature of such ensembles and their foundational principles, highlight advances in testing key ideas, and discuss applications of ensembles to analyze the collective behavior of granular materials.

preprint2013arXiv

Capillary fracture of soft gels

A liquid droplet resting on a soft gel substrate can deform that substrate to the point of material failure, whereby fractures develop on the gel surface that propagate outwards from the contact-line in a starburst pattern. In this paper, we characterize i) the initiation process in which the number of arms in the starburst is controlled by the ratio of surface tension contrast to the gel's elastic modulus and ii) the propagation dynamics showing that once fractures are initiated they propagate with a universal power law $L\propto t^{3/4}$. We develop a model for crack initiation by treating the gel as a linear elastic solid and computing the deformations within the substrate from the liquid/solid wetting forces. The elastic solution shows that both the location and magnitude of the wetting forces are critical in providing a quantitative prediction for the number of fractures and, hence, an interpretation of the initiation of capillary fractures. This solution also reveals that the depth of the gel is an important factor in the fracture process, as it can help mitigate large surface tractions; this finding is confirmed with experiments. We then develop a model for crack propagation by considering the transport of an inviscid fluid into the fracture tip of an incompressible material, and find that a simple energy-conservation argument can explain the observed material-independent power law. We compare predictions for both linear elastic and neo-Hookean solids finding that the latter better explains the observed exponent.

preprint2013arXiv

Correlations between electrical and mechanical signals during granular stick-slip events

Powders and grains exhibit unpredictable jamming-to-flow transitions that manifest themselves on geophysical scales in catastrophic slip events such as landslides and earthquakes, and on laboratory/industrial scales in profound processing difficulties. Over the past few years, insight into these transitions has been provided by new evidence that slip events may accompanied, or even preceded, by electrical effects. In the present work, we quantify the correlation between slip and the separation of electrical charges, using an archetypal granular material: photoelastic polymers. We measure a strong correlation between material displacement, acoustic emissions, and voltage. We find that the generation of voltage is associated with surface, rather than bulk properties of the granular materials. While voltage precursors are only occasionally observed in this system, there is some asymmetry in the cross-correlation between the slip and voltage signals that indicates differences between the pre-slip and post-slip dynamics.

preprint2013arXiv

Surfactant Spreading on a Thin Liquid Film: Reconciling Models and Experiments

The spreading dynamics of surfactant molecules on a thin fluid layer is of both fundamental and practical interest. A mathematical model formulated by Gaver and Grotberg 1990 describing the spreading of a single layer of insoluble surfactant has become widely accepted, and several experiments on axisymmetric spreading have confirmed its predictions for both the height profile of the free surface and the spreading exponent (the radius of the circular area covered by surfactant grows as t^{1/4}). However, these prior experiments have primarily utilized surfactant quantities exceeding (sometimes far exceeding) a monolayer. In this paper, we report that this regime is characterized by a mismatch between the timescales of the experiment and model, and additionally find that the spatial distribution of surfactant molecules differs substantially from the model prediction. For experiments performed in the monolayer regime for which the model was developed, the surfactant layer is observed to have a spreading exponent of less than 1/10, far below the predicted value, and the surfactant distribution is also in disagreement. These findings suggest that the model is inadequate for describing the spreading of insoluble surfactants on thin fluid layers.

preprint2012arXiv

Acoustic measurement of a granular density of modes

In glasses and other disordered materials, measurements of the vibrational density of states reveal that an excess number of long-wavelength (low-frequency) modes, as compared to the Debye scaling seen in crystalline materials, is associated with a loss of mechanical rigidity. In this paper, we present a novel technique for measuring the density of modes (DOM) in a real granular material, in which we mimic thermal excitations using white noise acoustic waves. The resulting vibrations are detected with piezoelectric sensors embedded inside a subset of the particles, from which we are able to compute the DOM via the spectrum of the velocity autocorrelation function, a technique previously applied in thermal systems. The velocity distribution for individual particles is observed to be Gaussian, but the ensemble distribution is non-Gaussian due to varying widths of the individual distributions. In spite of this deviation from a true thermal system, we find that the DOM exhibits several thermal-like features, including Debye scaling in a compressed hexagonally ordered packing, and an increase in low-frequency modes as the confining pressure is decreased. In disordered packings, we find that a characteristic frequency $f_c$ increases with pressure, but more weakly than has been observed in simulations of frictionless packings.

preprint2012arXiv

Equilibrating temperature-like variables in jammed granular subsystems

Although jammed granular systems are athermal, several thermodynamic-like descriptions have been proposed which make quantitative predictions about the distribution of volume and stress within a system and provide a corresponding temperature-like variable. We perform experiments with an apparatus designed to generate a large number of independent, jammed, two-dimensional configurations. Each configuration consists of a single layer of photoelastic disks supported by a gentle layer of air. New configurations are generated by alternately dilating and re-compacting the system through a series of boundary displacements. Within each configuration, a bath of particles surrounds a smaller subsystem of particles with a different inter-particle friction coefficient than the bath. The use of photoelastic particles permits us to find all particle positions as well as the vector forces at each inter-particle contact. By comparing the temperature-like quantities in both systems, we find compactivity (conjugate to the volume) does not equilibrate between the systems, while the angoricity (conjugate to the stress) does. Both independent components of the angoricity are linearly dependent on the hydrostatic pressure, in agreement with predictions of the stress ensemble.

preprint2012arXiv

Flow-Driven Formation of Solid-like Microsphere Heaps

We observe the formation of heaps of repulsive microspheres, created by flowing a colloidal microsphere suspension towards a flat-topped ridge placed within a quasi two-dimensional microfluidic channel. This configuration allows for both shear and normal forces on the microspheres in contact with the ridge. The heaps, which form against the ridge, are characterized by two distinct phases: a solid-like bulk phase in the interior, and a highly-fluctuating, liquid-like state which exists along its leading edge. We observe that heaps only form above a critical flow velocity, (v_c), and that they are destroyed by thermal rearrangements when the flow ceases. We monitor the dynamics of heap formation using fluorescence video microscopy, measuring the heap volume and the angle of repose in response to microsphere deposition and erosion processes. We find that the steady-state angle of repose, (θ_f), increases as a function of inflow velocity, (v_\infty), with a functional form (θ_f \propto \sqrt{v_\infty - v_c}).

preprint2012arXiv

Granular segregation in dense systems: the role of statistical mechanics and entropy

Granular segregation is ubiquitous in industrial, geological or daily-life context, but there is still no unifying theoretical approach. In this review, we examine two examples of granular segregation -- shallow rapid flows and rotating drums -- which suggest that the dynamics of systems at intermediate and high density might be amenable to a statistical mechanics approach.

preprint2012arXiv

The Influence of Network Topology on Sound Propagation in Granular Materials

Granular materials, whose features range from the particle scale to the force-chain scale to the bulk scale, are usually modeled as either particulate or continuum materials. In contrast with either of these approaches, network representations are natural for the simultaneous examination of microscopic, mesoscopic, and macroscopic features. In this paper, we treat granular materials as spatially-embedded networks in which the nodes (particles) are connected by weighted edges obtained from contact forces. We test a variety of network measures for their utility in helping to describe sound propagation in granular networks and find that network diagnostics can be used to probe particle-, curve-, domain-, and system-scale structures in granular media. In particular, diagnostics of meso-scale network structure are reproducible across experiments, are correlated with sound propagation in this medium, and can be used to identify potentially interesting size scales. We also demonstrate that the sensitivity of network diagnostics depends on the phase of sound propagation. In the injection phase, the signal propagates systemically, as indicated by correlations with the network diagnostic of global efficiency. In the scattering phase, however, the signal is better predicted by meso-scale community structure, suggesting that the acoustic signal scatters over local geographic neighborhoods. Collectively, our results demonstrate how the force network of a granular system is imprinted on transmitted waves.

preprint2012arXiv

Trajectory entanglement in dense granular materials

The particle-scale dynamics of granular materials have commonly been characterized by the self-diffusion coefficient $D$. However, this measure discards the collective and topological information known to be an important characteristic of particle trajectories in dense systems. Direct measurement of the entanglement of particle space-time trajectories can be obtained via the topological braid entropy $\Sbraid$, which has previously been used to quantify mixing efficiency in fluid systems. Here, we investigate the utility of $\Sbraid$ in characterizing the dynamics of a dense, driven granular material at packing densities near the static jamming point $ϕ_J$. From particle trajectories measured within a two-dimensional granular material, we typically observe that $\Sbraid$ is well-defined and extensive. However, for systems where $ϕ\gtrsim 0.79$, we find that $\Sbraid$ (like $D$) is not well-defined, signifying that these systems are not ergodic on the experimental timescale. Both $\Sbraid$ and $D$ decrease with either increasing packing density or confining pressure, independent of the applied boundary condition. The related braiding factor provides a means to identify multi-particle phenomena such as collective rearrangements. We discuss possible uses for this measure in characterizing granular systems.

preprint2011arXiv

Equipartition of rotational and translational energy in a dense granular gas

Experiments quantifying the rotational and translational motion of particles in a dense, driven, 2D granular gas floating on an air table reveal that kinetic energy is divided equally between the two translational and one rotational degrees of freedom. This equipartition persists when the particle properties, confining pressure, packing density, or spatial ordering are changed. While the translational velocity distributions are the same for both large and small particles, the angular velocity distributions scale with the particle radius. The probability distributions of all particle velocities have approximately exponential tails. Additionally, we find that the system can be described with a granular Boyle's Law with a van der Waals-like equation of state. These results demonstrate ways in which conventional statistical mechanics can unexpectedly apply to non-equilibrium systems.

preprint2010arXiv

Equilibration of granular subsystems

We experimentally investigate the steady states of two granular assemblies differing in their material properties and allowed to exchange volume with each other under external agitation in the vicinity of their jamming transition. We extract the statistics of various static and dynamic quantities, and uncover a materials-independent relationship between the average packing fraction and its fluctuations. This relationship defines an intensive parameter which decouples from the volume statistics, and remarkably takes the same value in both subsystems. We also observe that an effective diffusion coefficient also takes the same value in each subsystem, even as the structural relaxation time increases over several orders of magnitude. These observations provide strong constraints on the eventual establishment of a granular equation of state.

preprint2010arXiv

Fluorescent visualization of a spreading surfactant

The spreading of surfactants on thin films is an industrially and medically important phenomenon, but the dynamics are highly nonlinear and visualization of the surfactant dynamics has been a long-standing experimental challenge. We perform the first quantitative, spatiotemporally-resolved measurements of the spreading of an insoluble surfactant on a thin fluid layer. During the spreading process, we directly observe both the radial height profile of the spreading droplet and the spatial distribution of the fluorescently-tagged surfactant. We find that the leading edge of spreading circular layer of surfactant forms a Marangoni ridge in the underlying fluid, with a trough trailing the ridge as expected. However, several novel features are observed using the fluorescence technique, including a peak in the surfactant concentration which trails the leading edge, and a flat, monolayer-scale spreading film which differs from concentration profiles predicted by current models. Both the Marangoni ridge and surfactant leading edge can be described to spread as $R \propto t^δ$. We find spreading exponents, $δ_H \approx 0.30$ and $δ_Γ\approx 0.22$ for the ridge peak and surfactant leading edge, respectively, which are in good agreement with theoretical predictions of $δ= 1/4$. In addition, we observe that the surfactant leading edge initially leads the peak of the Marangoni ridge, with the peak later catching up to the leading edge.

preprint2010arXiv

Local origins of volume fraction fluctuations in dense granular materials

Fluctuations of the local volume fraction within granular materials have previously been observed to decrease as the system approaches jamming. We experimentally examine the role of boundary conditions and inter-particle friction $μ$ on this relationship for a dense granular material of bidisperse particles driven under either constant volume or constant pressure. Using a radical Voronoï tessellation, we find the variance of the local volume fraction $ϕ$ monotonically decreases as the system becomes more dense, independent of boundary condition and $μ$. We examine the universality and origins of this trend using experiments and the recent granocentric model \cite{Clusel-2009-GMR,Corwin-2010-MRP}, modified to draw particle locations from an arbitrary distribution ${\cal P}(s)$ of neighbor distances $s$. The mean and variance of the observed ${\cal P}(s)$ are described by a single length scale controlled by $\bar ϕ$. Through the granocentric model, we observe that diverse functional forms of ${\cal P}(s)$ all produce the trend of decreasing fluctuations, but only the experimentally-observed ${\cal P}(s)$ provides quantitative agreement with the measured $ϕ$ fluctuations. Thus, we find that both ${\cal P}(s)$ and ${\cal P}(ϕ)$ encode similar information about the ensemble of observed packings, and are connected to each other by the local granocentric model.

preprint2010arXiv

Scalar conservation laws with nonconstant coefficients with application to particle size segregation in granular flow

Granular materials will segregate by particle size when subjected to shear, as occurs, for example, in avalanches. The evolution of a bidisperse mixture of particles can be modeled by a nonlinear first order partial differential equation, provided the shear (or velocity) is a known function of position. While avalanche-driven shear is approximately uniform in depth, boundary-driven shear typically creates a shear band with a nonlinear velocity profile. In this paper, we measure a velocity profile from experimental data and solve initial value problems that mimic the segregation observed in the experiment, thereby verifying the value of the continuum model. To simplify the analysis, we consider only one-dimensional configurations, in which a layer of small particles is placed above a layer of large particles within an annular shear cell and is sheared for arbitrarily long times. We fit the measured velocity profile to both an exponential function of depth and a piecewise linear function which separates the shear band from the rest of the material. Each solution of the initial value problem is non-standard, involving curved characteristics in the exponential case, and a material interface with a jump in characteristic speed in the piecewise linear case.

preprint2010arXiv

Universal shapes formed by two interacting cracks

We investigate the origins of the widely-observed "en passant" crack pattern which forms through interactions between two approaching cracks. A rectangular elastic plate is notched on each long side and then subjected to quasistatic uniaxial strain from the short side. The two cracks propagate along approximately straight paths until they pass each other, after which they curve and release a lenticular fragment. We find that for materials with diverse mechanical properties, the shape of this fragment has an aspect ratio of 2:1, with the length scale set by the initial crack offset $s$ and the time scale set by the ratio of $s$ to the pulling velocity. The cracks have a universal square root shape which we understand using a simple geometric model of the crack-crack interaction.

preprint2009arXiv

Mixing and segregation rates in sheared granular materials

The size-segregation of granular materials, a process colloquially known as the Brazil Nut Effect, has generally been thought to proceed faster the greater the size difference of the particles. We experimentally investigate sheared bidisperse granular materials as a function of the size ratio of the two species, and find that the mixing rate at low confining pressure behaves as expected from percolation-based arguments. However, we also observe an anomalous effect for the re-segregation rates, wherein particles of both dissimilar and similar sizes segregate more slowly than intermediate particle size ratios. Combined with the fact that increasing the confining pressure significantly suppresses both mixing and segregation rates of particles of dissimilar size, we propose that the anomalous behavior may be attributed to a species-dependent distribution of forces within the system.

preprint2009arXiv

Shear-driven size segregation of granular materials: modeling and experiment

Granular materials segregate by size under shear, and the ability to quantitatively predict the time required to achieve complete segregation is a key test of our understanding of the segregation process. In this paper, we apply the Gray-Thornton model of segregation (developed for linear shear profiles) to a granular flow with an exponential profile, and evaluate its ability to describe the observed segregation dynamics. Our experiment is conducted in an annular Couette cell with a moving lower boundary. The granular material is initially prepared in an unstable configuration with a layer of small particles above a layer of large particles. Under shear, the sample mixes and then re-segregates so that the large particles are located in the top half of the system in the final state. During this segregation process, we measure the velocity profile and use the resulting exponential fit as input parameters to the model. To make a direct comparison between the continuum model and the observed segregation dynamics, we locally map the measured height of the experimental sample (which indicates the degree of segregation) to the local packing density. We observe that the model successfully captures the presence of a fast mixing process and relatively slower re-segregation process, but the model predicts a finite re-segregation time, while in the experiment re-segregation occurs only exponentially in time.

preprint2007arXiv

Competition and bistability of ordered undulations and undulation chaos in inclined layer convection

Experimental and theoretical investigations of undulation patterns in high-pressure, inclined layer gas convection at a Prandtl number near unity are reported. Particular focus is given to the competition between the spatiotemporal chaotic state of undulation chaos and stationary patterns of ordered undulations. In experiments a competition and bistability between the two states is observed, with ordered undulations most prevalent at higher Rayleigh number. The spectral pattern entropy, spatial correlation lengths, and defect statistics are used to characterize the competing states. The experiments are complemented by a theoretical analysis of the Oberbeck-Boussinesq equations. The stability region of the ordered undulation as a function of their wavevectors and the Rayleigh number is obtained with Galerkin techniques. In addition, direct numerical simulations are used to investigate the spatiotemporal dynamics. In the simulations both ordered undulations and undulation chaos were observed dependent on initial conditions. Experiment and theory are found to agree well.