The role of deformability in determining the structural and mechanical properties of bubbles and emulsions

Elasticity and rheology of auxetic granular metamaterials

The flowing, jamming, and avalanche behavior of granular materials is satisfyingly universal and vexingly hard to tune: A granular flow is typically intermittent and will irremediably jam if too confined. Here, we show that granular metamaterials made from particles with a negative Poisson's ratio yield more easily and flow more smoothly than ordinary granular materials. We first create a collection of auxetic grains based on a re-entrant mechanism and show that each grain exhibits a negative Poisson's ratio regardless of the direction of compression. Interestingly, we find that the elastic and yielding properties are governed by the high compressibility of granular metamaterials: At a given confinement, they exhibit lower shear modulus, lower yield stress, and more frequent, smaller avalanches than materials made from ordinary grains. We further demonstrate that granular metamaterials promote flow in more complex confined geometries, such as intruder and hopper geometries, even when the packing contains only a fraction of auxetic grains. Moreover, auxetic granular metamaterials exhibit enhanced impact absorption. Our findings blur the boundary between complex fluids and metamaterials and could help in scenarios that involve process, transport, and reconfiguration of granular materials.

A soft departure from jamming: the compaction of deformable granular matter under high pressures

The high-pressure compaction of three dimensional granular packings is simulated using a bonded particle model (BPM) to capture linear elastic deformation. In the model, grains are represented by a collection of point particles connected by bonds. A simple multibody interaction is introduced to control Poisson's ratio and the arrangement of particles on the surface of a grain is varied to model both high- and low-frictional grains. At low pressures, the growth in packing fraction and coordination number follow the expected behavior near jamming and exhibit friction dependence. As the pressure increases, deviations from the low-pressure power-law scaling emerge after the packing fraction grows by approximately 0.1 and results from simulations with different friction coefficients converge. These results are compared to predictions from traditional discrete element method simulations which, depending on the definition of packing fraction and coordination number, may only differ by a factor of two. As grains deform under compaction, the average volumetric strain and asphericity, a measure of the change in the shape of grains, are found to grow as power laws and depend heavily on the Poisson's ratio of the constituent solid. Larger Poisson's ratios are associated with less volumetric strain and more asphericity and the apparent power-law exponent of the asphericity may vary. The elastic properties of the packed grains are also calculated as a function of packing fraction. In particular, we find the Poisson's ratio near jamming is 1/2 but decreases to around 1/4 before rising again as systems densify.

Jamming on convex deformable surfaces

Jamming is a fundamental transition that governs the behavior of particulate media, including sand, foams and dense suspensions. Upon compression, such media change from freely flowing to a disordered, marginally stable solid that exhibits non-Hookean elasticity. While the jamming process is well established for fixed geometries, the nature and dynamics of jamming for a diverse class of soft materials and deformable substrates, including emulsions and biological matter, remains unknown. Here we propose a new scenario, metric jamming, where rigidification occurs on a surface that has been deformed from its ground state. Unlike classical jamming processes that exhibit discrete mechanical transitions, surprisingly we find that metric jammed states possess mechanical properties continuously tunable between those of classically jammed and conventional elastic media. The compact and curved geometry significantly alters the vibrational spectra of the structures relative to jamming in flat Euclidean space, and metric jammed systems also possess new types of vibrational mode that couple particle and shape degrees of freedom. Our work provides a theoretical framework that unifies our understanding of solidification processes that take place on deformable media and lays the groundwork to exploit jamming for the control and stabilization of shape in self-assembly processes.

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 ∼ (w/σ_avg - k)^β, where σ_avg is the average particle diameter, kσ_avg is an offset where Q ∼ 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 λ → 0 limit and d - 3/2 in the λ → ∞ 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 ∼ k_max in the hard-particle limit, where k_max ∼ 3.5 is larger for λ → ∞ compared to that for λ → 0. Finally, we show that the simulations of hopper flows of deformable particles in the λ → ∞ limit recapitulate the experimental results for quasi-2D hopper flows of oil droplets in water.

Structural Measures as Guides to Ultrastable States in Overjammed Packings

Jammed, disordered packings of given sets of particles possess a multitude of equilibrium states with different mechanical properties. Identifying and constructing desired states, e.g., of superior stability, is a complex task. Here, we show that in two-dimensional particle packings the energy of all metastable states (inherent structures) is reliably classified by simple scalar measures of local steric packing. These structural measures are insensitive to the particle interaction potential and so robust that they can be used to guide a modified swap algorithm that anneals polydisperse packings toward low-energy metastable states exceptionally fast. The low-energy states are extraordinarily stable against applied shear, so that the approach also efficiently identifies ultrastable packings.

The structural, vibrational, and mechanical properties of jammed packings of deformable particles in three dimensions

We investigate the structural, vibrational, and mechanical properties of jammed packings of deformable particles with shape degrees of freedom in three dimensions (3D). Each 3D deformable particle is modeled as a surface-triangulated polyhedron, with spherical vertices whose positions are determined by a shape-energy function with terms that constrain the particle surface area, volume, and curvature, and prevent interparticle overlap. We show that jammed packings of deformable particles without bending energy possess low-frequency, quartic vibrational modes, whose number decreases with increasing asphericity and matches the number of missing contacts relative to the isostatic value. In contrast, jammed packings of deformable particles with non-zero bending energy are isostatic in 3D, with no quartic modes. We find that the contributions to the eigenmodes of the dynamical matrix from the shape degrees of freedom are significant over the full range of frequency and shape parameters for particles with zero bending energy. We further show that the ensemble-averaged shear modulus 〈G〉 scales with pressure P as 〈G〉 ∼ Pβ, with β ≈ 0.75 for jammed packings of deformable particles with zero bending energy. In contrast, β ≈ 0.5 for packings of deformable particles with non-zero bending energy, which matches the value for jammed packings of soft, spherical particles with fixed shape. These studies underscore the importance of incorporating particle deformability and shape change when modeling the properties of jammed soft materials.

Embryonic Tissues as Active Foams

The physical state of embryonic tissues emerges from non-equilibrium, collective interactions among constituent cells. Cellular jamming, rigidity transitions and characteristics of glassy dynamics have all been observed in multicellular systems, but it is unclear how cells control these emergent tissue states and transitions, including tissue fluidization. Combining computational and experimental methods, here we show that tissue fluidization in posterior zebrafish tissues is controlled by the stochastic dynamics of tensions at cell-cell contacts. We develop a computational framework that connects cell behavior to embryonic tissue dynamics, accounting for the presence of extracellular spaces, complex cell shapes and cortical tension dynamics. We predict that tissues are maximally rigid at the structural transition between confluent and non-confluent states, with actively-generated tension fluctuations controlling stress relaxation and tissue fluidization. By directly measuring strain and stress relaxation, as well as the dynamics of cell rearrangements, in elongating posterior zebrafish tissues, we show that tension fluctuations drive active cell rearrangements that fluidize the tissue. These results highlight a key role of non-equilibrium tension dynamics in developmental processes.

Viscoelasticity of non-colloidal hydrogel particle suspensions at the liquid-solid transition

Suspensions of soft particles transition from a viscous fluid to a soft material upon increases in phase volume. The criteria defining the transition to this jammed state are difficult to define due to the porous and deformable nature of soft particles. Here, we characterise the rheology of aqueous suspensions of industrially relevant non-colloidal, polydisperse, frictional agarose microgels and evaluate shear and viscoelastic behaviour across a range of phase volumes from the dilute regime to the highly concentrated regime. In order to model the viscoelastic response of suspensions without free fitting parameters, the random close packing volume fraction (φrcp) and the particle modulus are determined, respectively, from particle size distribution measurements and direct measurements of reduced elastic modulus of individual particles (Erp) using Atomic Force Microscopy. It is found that at φrcp, previously shown to correspond to divergence of the viscosity, also corresponds to the suspension transition from a viscous to viscoelastic fluid. However, the transition to a jammed solid-like state (φj) occurs at phase volumes exceeding this value (i.e. φj > φrcp). The suspension modulus and its sudden growth at φj are well-predicted by the Evans and Lips model that incorporates the Erp of the hydrogel particles. This rheological behaviour showing a dual transition is reminiscent of two families of systems: (i) colloidal suspensions and (ii) frictional-adhesive non-colloidal suspensions. However, it does not strictly follow either case. We propose that the width of the transition region is dictated by frictional contact, particle size distribution and particle modulus, and plan to further probe this in future work.

Adhesion as a trigger of droplet polarization in flowing emulsions

Tissues are subjected to large external forces and undergo global deformations during morphogenesis. We use synthetic analogues of tissues to study the impact of cell-cell adhesion on the response of cohesive cellular assemblies under such stresses. In particular, we use biomimetic emulsions in which the droplets are functionalized in order to exhibit specific droplet-droplet adhesion. We flow these emulsions in microfluidic constrictions and study their response to this forced deformation via confocal microscopy. We find that the distributions of avalanche sizes are conserved between repulsive and adhesive droplets. However, adhesion locally impairs the rupture of droplet-droplet contacts, which in turn pulls on the rearranging droplets. As a result, adhesive droplets are a lot more deformed along the axis of elongation in the constriction. This finding could shed light on the origin of polarization processes during morphogenesis.

Poisson-bracket formulation of the dynamics of fluids of deformable particles

Using the Poisson-bracket method, we derive continuum equations for a fluid of deformable particles in two dimensions. Particle shape is quantified in terms of two continuum fields: an anisotropy density field that captures the deformations of individual particles from regular shapes and a shape tensor density field that quantifies both particle elongation and nematic alignment of elongated shapes. We explicitly consider the example of a dense biological tissue as described by the Vertex model energy, where cell shape has been proposed as a structural order parameter for a liquid-solid transition. The hydrodynamic model of biological tissue proposed here captures the coupling of cell shape to flow and provides a starting point for modeling the rheology of dense tissue.

Structural characterization and statistical properties of jammed soft ellipsoid packing

The jamming transition and jammed packing structures of hydrogel soft ellipsoids are studied using magnetic resonance imaging techniques. As the packing fraction increases, the fluctuation of local free volume decreases and the fluctuation of particle deformation increases. Effective thermodynamic quantities are obtained by characterizing these fluctuations using k-gamma distributions based on an underlying statistical model for granular materials. Surprisingly, the two granular temperatures measuring the relative fluctuations of both free volume and particle deformation remain basically unchanged as the packing fraction increases. The total configurational entropy is also approximately constant for packing with different packing fractions. The significantly different behaviors of these effective thermodynamic quantities compared with hard sphere systems are further attributed to a statistically affine structural transformation of the packing structures along with particle deformations when the packing fraction changes.

Size distribution dependence of collective relaxation dynamics in a two-dimensional wet foam

Foams can be ubiquitously observed in nature and in industrial products. Despite the relevance of their properties to deformation, fluidity, and collapse, all of which are essential for applications, there are few experimental studies of collective relaxation dynamics in a wet foam. Here, we directly observe how the relaxation dynamics changes with increasing liquid fraction in both monodisperse and polydisperse two-dimensional foams. As we increase the liquid fraction, we quantitatively characterize the slowing-down of the relaxation, and the increase of the correlation length. We also find two different relaxation modes which depend on the size distribution of the bubbles. It suggests that the bubbles which are simply near to each other play an important role in large rearrangements, not just those in direct contact. Finally, we confirm the generality of our experimental findings by a numerical simulation for the relaxation process of wet foams.

Scaling law for the kinetics of water imbibition in polydisperse foams

We investigated the kinetics of water imbibition in polydisperse foams. We used a Hele-Shaw cell, and horizontal imbibition was observed for a timescale of up to 10³ s in which the gravity effect was negligible. While several papers have reported kinetics for imbibition in foams, imbibition kinetics in polydisperse foams and its variations in longer timescales are not well understood. The tip position of imbibition was proportional to the square root of time in the initial stage of imbibition, but it showed plateauing in the late stage of imbibition. We evaluated the proportional constant A in the initial stage of imbibition as a kinetic constant for the time-dependent increase in the tip position, which showed a clear dependency on the initial and final water volume fractions in the foams. Conversely, the mean initial radius of the curvature and the channel length in the Plateau borders did not show any clear correlations with A, although both valuables are frequently used in modeling for liquid imbibition in foams. On the basis of the t^1/2 dependence, the correlation of A with the water volume fraction and the increase in the water volume fraction during imbibition, we proposed a simple equation to describe the tip position over the entire period of imbibition. We used them to scale all of the experimental data, which showed good agreement with the theoretical line. This clearly showed that the water volume fraction in the foams during imbibition was the key factor to quantitatively describe the rate of water imbibition. Features in the kinetics of imbibition were discussed.

Horizontal imbibition of water in foams is scaled well by a simple mathematical expression that considers t^1/2 dependence and changes in volume fraction of water in foams.

Strain localization and failure of disordered particle rafts with tunable ductility during tensile deformation

Quasi-static tensile experiments were performed for a model disordered solid consisting of a two-dimensional raft of polydisperse floating granular particles with capillary attractions. The ductility is tuned by controlling the capillary interaction range, which varies with the particle size. During the tensile tests, after an initial period of elastic deformation, strain localization occurs and leads to the formation of a shear band at which the pillar later fails. In this process, small particles with long-ranged interactions can endure large plastic deformation without forming significant voids, while large particles with short-range interactions fail dramatically by fracturing at small deformation. Particle-level structure was measured, and the strain-localized region was found to have higher structural anisotropy than the bulk. Local interactions between anisotropic sites and particle rearrangements were the main mechanisms driving strain localization and the subsequent failure, and significant differences of such interactions exist between ductile and brittle behaviors.

Compaction Model for Highly Deformable Particle Assemblies

The compaction behavior of deformable grain assemblies beyond jamming remains bewildering, and existing models that seek to find the relationship between the confining pressure P and solid fraction ϕ end up settling for empirical strategies or fitting parameters. Using a coupled discrete-finite element method, we analyze assemblies of highly deformable frictional grains under compression. We show that the solid fraction evolves nonlinearly from the jamming point and asymptotically tends to unity. Based on the micromechanical definition of the granular stress tensor, we develop a theoretical model, free from ad hoc parameters, correctly mapping the evolution of ϕ with P. Our approach unveils the fundamental features of the compaction process arising from the joint evolution of grain connectivity and the behavior of single representative grains. This theoretical framework also allows us to deduce a bulk modulus equation showing an excellent agreement with our numerical data.

Depletion attraction impairs the plasticity of emulsions flowing in a constriction

We study the elasto-plastic behavior of dense attractive emulsions under a mechanical perturbation. The attraction is introduced through non-specific depletion interactions between the droplets and is controlled by changing the concentration of surfactant micelles in the continuous phase. We find that such attractive forces are not sufficient to induce any measurable modification on the scalings between the local packing fraction and the deformation of the droplets. However, when the emulsions are flowed through 2D microfluidic constrictions, we uncover a measurable effect of attraction on their elasto-plastic response. Indeed, we measure higher levels of deformation inside the constriction for attractive droplets. In addition, we show that these measurements correlate with droplet rearrangements that are spatially delayed in the constriction for higher attraction forces.

Pressure Dependent Shear Response of Jammed Packings of Frictionless Spherical Particles

The mechanical response of packings of purely repulsive, spherical particles to athermal, quasistatic simple shear near jamming onset is highly nonlinear. Previous studies have shown that, at small pressure p, the ensemble-averaged static shear modulus ⟨G-G_{0}⟩ scales with p^{α}, where α≈1, but above a characteristic pressure p^{**}, ⟨G-G_{0}⟩∼p^{β}, where β≈0.5. However, we find that the shear modulus G^{i} for an individual packing typically decreases linearly with p along a geometrical family where the contact network does not change. We resolve this discrepancy by showing that, while the shear modulus does decrease linearly within geometrical families, ⟨G⟩ also depends on a contribution from discontinuous jumps in ⟨G⟩ that occur at the transitions between geometrical families. For p>p^{**}, geometrical-family and rearrangement contributions to ⟨G⟩ are of opposite signs and remain comparable for all system sizes. ⟨G⟩ can be described by a scaling function that smoothly transitions between two power-law exponents α and β. We also demonstrate the phenomenon of compression unjamming, where a jammed packing unjams via isotropic compression.