Search for: All records

Capillary droplets form due to surface tension when two immiscible fluids are mixed. We describe the motion of gravity-driven capillary droplets flowing through narrow constrictions and obstacle arrays in both simulations and experiments. Our new capillary deformable particle model recapitulates the shape and velocity of single oil droplets in water as they pass through narrow constrictions in microfluidic chambers. Using this experimentally validated model, we simulate the flow and clogging of single capillary droplets in narrow channels and obstacle arrays and find several important results. First, the capillary droplet speed profile is nonmonotonic as the droplet exits the narrow orifice, and we can tune the droplet properties so that the speed overshoots the terminal speed far from the constriction. Second, in obstacle arrays, we find that extremely deformable droplets can wrap around obstacles, which leads to decreased average droplet speed in the continuous flow regime and increased probability for clogging in the regime where permanent clogs form. Third, the wrapping mechanism causes the clogging probability in obstacle arrays to become nonmonotonic with surface tension Γ. At large Γ, the droplets are nearly rigid and the clogging probability is large since the droplets can not squeeze through the gaps between obstacles. With decreasing Γ, the clogging probability decreases as the droplets become more deformable. However, in the small-Γ limit, the clogging probability increases since the droplets are extremely deformable and wrap around the obstacles. The results from these studies are important for developing a predictive understanding of capillary droplet flows through complex and confined geometries. 

Distributions of strictly positive numbers are common and can be characterized by standard statistical measures such as mean, standard deviation, and skewness. We demonstrate that for these distributions the skewnessD₃is bounded from below by a function of the coefficient of variation (CoV) δ as D₃ > δ − 1/δ. The results are extended to any distribution that is bounded with minimum value x_min and/or bounded with maximum value x_max. We build on the results to provide bounds for kurtosis D₄, and conjecture analogous bounds exists for higher statistical moments. 

We study particle-scale motion in sheared highly polydisperse amorphous materials, in which the largest particles are as much as ten times the size of the smallest. We find strikingly different behavior from the more commonly studied amorphous systems with low polydispersity. In particular, an analysis of the nonaffine motion of particles reveals qualitative differences between large and small particles: The smaller particles have dramatically more nonaffine motion, which is induced by the presence of the large particles. We characterize how the nonaffine motion changes from the low- to high-polydispersity regimes. We further demonstrate a quantitative way to distinguish between “large” and “small” particles in systems with broad distributions of particle sizes. A macroscopic consequence of the nonaffine motion is a decrease in the energy dissipation rate for highly polydisperse samples, which is due both to a geometric consequence of the changing jamming conditions for higher polydispersity and to the changing character of nonaffine motion. 

Abstract Magnetorheological fluids (MRF) are smart materials of increasing interest due to their great versatility in mechanical and mechatronic systems. As main rheological features, MRFs must present low viscosity in the absence of magnetic field (0.1–1.0 Pa.s) and high yield stress (50–100 kPa) when magnetized, in order to optimize the magnetorheological effect. Such properties, in turn, are directly influenced by the composition, volume fraction, size, and size distribution (polydispersity) of the particles, the latter being an important piece in the improvement of these main properties. In this context, the present work aims to analyze, through experiments and simulations, the influence of polydispersity on the maximum packing fraction, on the yield stress under field (on-state) and on the plastic viscosity in the absence of field (off-state) of concentrated MRF (φ= 48.5 vol.%). Three blends of carbonyl iron powder (CIP) in polyalphaolefin oil were prepared. These blends have the same mode, but different polydispersity indexes (α), ranging from 0.46 to 1.44. Separate simulations show that the random close packing fraction increases from about 68% to 80% as the polydispersity indexes increase over this range. The on-state yield stress, in turn, is raised from 30 ± 0.5 kPa to 42 ± 2 kPa (B≈ 0.57 T) and the off-state plastic viscosity, is reduced from 4.8 Pa.s to 0.5 Pa.s. Widening the size distributions, as is well known in the literature, increases packing efficiency and reduces the viscosity of concentrated dispersions, but beyond that, it proved to be a viable way to increase the magnetorheological effect of concentrated MRF. The Brouwers model, which considers the void fraction in suspensions of particles with lognormal distribution, was proposed as a possible hypothesis to explain the increase in yield stress under magnetic field. 

We have studied shear deformation of binary Lennard-Jones glasses to investigate the extent to which the transient part of the stress strain curves is invariant when the thermodynamic state point is varied along an isomorph. Shear deformations were carried out on glass samples of varying stability, determined by cooling rate, and at varying strain rates, at state points deep in the glass. Density changes up to and exceeding a factor of two were made. We investigated several different methods for generating isomorphs but none of the previously developed methods could generate sufficiently precise isomorphs given the large density changes and nonequilibrium situation. Instead, the temperatures for these higher densities were chosen to give state points isomorphic to the starting state point by requiring the steady-state flow stress for isomorphic state points to be invariant in reduced units. In contrast to the steady-state flow stress, we find that the peak stress on the stress strain curve is not invariant. The peak stress decreases by a few percent for each ten percent increase in density, although the differences decrease with increasing density. Analysis of strain profiles and nonaffine motion during the transient phase suggests that the root of the changes in peak stress is a varying tendency to form shear bands, with the largest tendency occurring at the lowest densities. We suggest that this reflects the effective steepness of the potential; a higher effective steepness gives a greater tendency to form shear bands. 

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. 

We conduct molecular dynamics simulations of a bidisperse Kob–Andersen (KA) glass former, modified to add in additional polydispersity. The original KA system is known to exhibit dynamical heterogeneity. Prior work defined propensity, the mean motion of a particle averaged over simulations reconstructing the initial positions of all particles but with randomized velocities. The existence of propensity shows that structure and dynamics are connected. In this paper, we study systems which mimic problems that would be encountered in measuring propensity in a colloidal glass former, where particles are polydisperse (they have slight size variations). We mimic polydispersity by altering the bidisperse KA system into a quartet consisting of particles both slightly larger and slightly smaller than the parent particles in the original bidisperse system. We then introduce errors into the reconstruction of the initial positions that mimic mistakes one might make in a colloidal experiment. The mistakes degrade the propensity measurement, in some cases nearly completely; one no longer has an iso-configurational ensemble in any useful sense. Our results show that a polydisperse sample is suitable for propensity measurements provided one avoids reconstruction mistakes.