The drag force on a spherical intruder in dense granular shear flows is studied using discrete element method simulations. Three regimes of the intruder dynamics are observed depending on the magnitude of the drag force (or the corresponding intruder velocity) and the flow inertial number: a fluctuation-dominated regime for small drag forces; a viscous regime for intermediate drag forces; and an inertial (cavity formation) regime for large drag forces. The transition from the viscous regime (linear force-velocity relation) to the inertial regime (quadratic force-velocity relation) depends further on the inertial number. Despite these distinct intruder dynamics, we find a quantitative similarity between the intruder drag in granular shear flows and the Stokesian drag on a sphere in a viscous fluid for intruder Reynolds numbers spanning five orders of magnitude. Beyond this first-order description, a modified Stokes drag model is developed that accounts for the secondary dependence of the drag coefficient on the inertial number and the intruder size and density ratios. When the drag model is coupled with a segregation force model for intruders in dense granular flows, it is possible to predict the velocity of gravity-driven segregation of an intruder particle in shear flow simulations.
more »
« less
Granular flows in drums of non-uniform widths
We study how channel width variations influence the dynamics of free-surface granular flows. For this purpose, we extend a continuum model framework to granular flows passing through channels that narrow or widen. Our theory uses a linearized approximation to an established dense granular flow rheology and a Coulomb friction law to model interaction between flow and sidewalls. We test the theoretical predictions using two novel 40 cm-diameter drums (convex and concave) filled halfway with 2 mm diameter particles rotated at rates in which the shear layer remains shallow and dense. We apply particle tracking velocimetry to enable quantitative comparisons between experimental data and theoretical predictions. We find that our experimental kinematics and energy profiles largely agree with the theoretical predictions. In general, flows through narrowing channels are faster and deeper than flows through widening channels. The influence of width variations grows with increasing flow speed, and the form of the rate dependence changes fundamentally as the regime changes from one in which kinetic energy is dissipated locally to one in which it is advected downstream. For both regimes, theoretical scaling analysis leads us to experimentally validated power laws, in which the exponent depends on the flow regime, and the multiplicative coefficient depends on channel geometry alone. Finally, we discuss how the differences between theoretical predictions and experimental data may be useful for improving our understanding of flows through non-uniform channels.
more »
« less
- Award ID(s):
- 2127476
- PAR ID:
- 10500028
- Publisher / Repository:
- Journal of Fluid Mechanics
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 954
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Active nematic liquid crystals have the remarkable ability to spontaneously deform and flow in the absence of any external driving force. While living materials with orientational order, such as the mitotic spindle, can self-assemble in quiescent active phases, reconstituted active systems often display chaotic, periodic, or circulating flows under confinement. Quiescent active nematics are, therefore, quite rare, despite the prediction from active hydrodynamic theory that confinement between two parallel plates can suppress flows. This spontaneous flow transition—named the active Fréedericksz transition by analogy with the conventional Fréedericksz transition in passive nematic liquid crystals under a magnetic field—has been a cornerstone of the field of active matter. Here, we report experimental evidence that confinement in spherical droplets can stabilize the otherwise chaotic dynamics of a 3D extensile active nematics, giving rise to a quiescent—yet still out-of-equilibrium—nematic liquid crystal. The active nematics spontaneously flow when confined in larger droplets. The composite nature of our model system composed of extensile bundles of microtubules and molecular motors dispersed in a passive colloidal liquid crystal allows us to demonstrate how the interplay of activity, nematic elasticity, and confinement impacts the spontaneous flow transition. The critical diameter increases when motor concentration decreases or nematic elasticity increases. Experiments and simulations also demonstrate that the critical confinement depends on the confining geometry, with the critical diameter in droplets being larger than the critical width in channels. Biochemical assays reveal that neither confinement nor nematic elasticity impacts the energy-consumption rate, confirming that the quiescent active phase is the stable out-of-equilibrium phase predicted theoretically. Further experiments in dense arrays of monodisperse droplets show that fluctuations in the droplet composition can smooth the flow transition close to the critical diameter. In conclusion, our work provides experimental validation of the active Fréedericksz transition in 3D active nematics, with potential applications in human health, ecology, and soft robotics. Published by the American Physical Society2024more » « less
-
Abstract Successful management of flooding and erosion hazards on floodplains depends on our ability to predict a river channel's shape and the lifespan during which it will continue to flow. Recent progress has improved our understanding of what sets the lifespan and width of single‐thread channels; the next challenge is to extend this knowledge to braided channels and their interwoven sub‐channels (threads). In this study, we investigate the lifespan and width of braided channel threads in a large experimental data set, coupled with particle‐image velocimetry‐derived measurements of riverbank erosion and accretion. We find that, unlike single‐thread channels, braided channels in the experiment do not exhibit an equilibrium between bank erosion and accretion. Instead, bank erosion outpaces lateral accretion, causing individual threads to widen and infill until they are abandoned. Thread lifespan is limited to the time it takes for threads to triple their width: tripling of the width yields enough bank material to aggrade more than half the channel depth, at which point flow is rerouted to a narrower thread. In consequence the width of active threads is limited to three times their initial width. Threshold channel theory accurately predicts the median thread width, which is roughly double the initial width and two‐thirds the limiting width. The results are consistent with existing field data and suggest that differential bank migration is sufficient to explain why braided channels show greater width variability and higher width‐to‐depth ratios than their single‐thread counterparts.more » « less
-
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.more » « less
-
Abstract. The width of valleys and channels affects the hydrology, ecology,and geomorphic functionality of drainage networks. In many studies, thewidth of valleys and/or channels (W) is estimated as a power-law function ofthe drainage area (A), W=kcAd. However, in fluvial systemsthat experience drainage reorganization, abrupt changes in drainage areadistribution can result in valley or channel widths that are disproportionalto their drainage areas. Such disproportionality may be more distinguishedin valleys than in channels due to a longer adjustment timescale forvalleys. Therefore, the valley width–area scaling in reorganized drainagesis expected to deviate from that of drainages that did not experiencereorganization. To explore the effect of reorganization on valley width–drainage areascaling, we studied 12 valley sections in the Negev desert, Israel,categorized into undisturbed, beheaded, and reversed valleys. We found thatthe values of the drainage area exponents, d, are lower in the beheadedvalleys relative to undisturbed valleys but remain positive. Reversedvalleys, in contrast, are characterized by negative d exponents, indicatingvalley narrowing with increasing drainage area. In the reversed category, wealso explored the independent effect of channel slope (S) through theequation W=kbAbSc, which yieldednegative and overall similar values for b and c. A detailed study in one reversed valley section shows that the valleynarrows downstream, whereas the channel widens, suggesting that, ashypothesized, the channel width adjusts faster to post-reorganizationdrainage area distribution. The adjusted narrow channel dictates the widthof formative flows in the reversed valley, which contrasts with the meaningfullywider formative flows of the beheaded valley across the divide. Thisdifference results in a step change in the unit stream power between thereversed and beheaded channels, potentially leading to a “width feedback”that promotes ongoing divide migration and reorganization. Our findings demonstrate that valley width–area scaling is a potential toolfor identifying landscapes influenced by drainage reorganization. Accountingfor reorganization-specific scaling can improve estimations of erosion ratedistributions in reorganized landscapes.more » « less
An official website of the United States government

