skip to main content

Explore Research Products in the NSF-PAR

Title: Collision rate of bidisperse, hydrodynamically interacting spheres settling in a turbulent flow
The collisions in a dilute polydisperse suspension of sub-Kolmogorov spheres with negligible inertia settling in a turbulent flow and interacting through hydrodynamics including continuum breakdown on close approach are studied. A statistically significant decrease in ideal collision rate without gravity is resolved via a Lagrangian stochastic velocity-gradient model at Taylor microscale Reynolds number larger than those accessible by current direct numerical simulation capabilities. This arises from the difference between the mean inward velocity and the root-mean-square particle relative velocity. Differential sedimentation, comparable to the turbulent shear relative velocity, but minimally influencing the sampling of the velocity gradient, diminishes the Reynolds number dependence and enhances the ideal collision rate i.e. the rate without interactions. The collision rate is retarded by hydrodynamic interactions between sphere pairs and is governed by non-continuum lubrication as well as full continuum hydrodynamic interactions at larger separations. The collision efficiency (ratio of actual to ideal collision rate) depends on the relative strength of differential sedimentation and turbulent shear, the size ratio of the interacting spheres and the Knudsen number (defined as the ratio of the mean-free path of the gas to the mean radius of the interacting spheres). We develop an analytical approximation to concisely report computed results across the parameter space. This accurate closed form expression could be a critical component in computing the evolution of the size distribution in applications such as water droplets in clouds or commercially valuable products in industrial aggregators.  more » « less
Award ID(s):
1803156
NSF-PAR ID:
10288861
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
912
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Journal of Fluid Mechanics)
    null (Ed.)
    Collisions in a dilute polydisperse suspension of spheres of negligible inertia interacting through non-continuum hydrodynamics and settling in a slow uniaxial compressional flow are studied. The ideal collision rate is evaluated as a function of the relative strength of gravity and uniaxial compressional flow and it deviates significantly from a linear superposition of these driving terms. This non-trivial behaviour is exacerbated by interparticle interactions based on uniformly valid non-continuum hydrodynamics, that capture non-continuum lubrication at small separations and full continuum hydrodynamic interactions at larger separations, retarding collisions driven purely by sedimentation significantly more than those driven purely by the linear flow. While the ideal collision rate is weakly dependent on the orientation of gravity with the axis of compression, the rate including hydrodynamic interactions varies by more than $100\,\%$ with orientation. This dramatic shift can be attributed to complex trajectories driven by interparticle interactions that prevent particle pairs from colliding or enable a circuitous path to collision. These and other important features of the collision process are studied in detail using trajectory analysis at near unity and significantly smaller than unity size ratios of the interacting spheres. For each case analysis is carried for a large range of relative strengths and orientations of gravity to the uniaxial compressional flow, and Knudsen numbers (ratio of mean free path of the media to mean radius). 
    more » « less
  2. ; ( , Journal of Rheology)
    The rheology of suspensions of rings (tori) rotating in an unbounded low Reynolds number simple shear flow is calculated using numerical simulations at dilute particle number densities ( n ≪ 1 ). Suspensions of non-Brownian rings are studied by computing pair interactions that include hydrodynamic interactions modeled using slender body theory and particle collisions modeled using a short-range repulsive force. Particle contact and hydrodynamic interactions were found to have comparable influences on the steady-state Jeffery orbit distribution. The average tilt of the ring away from the flow-vorticity plane increased during pairwise interactions compared to the tilt associated with Jeffery rotation and the steady-state orbit distribution. Particle stresses associated with the increased tilt during the interaction were found to be comparable to the stresses induced directly by particle contact forces and the hydrodynamic velocity disturbances of other particles. The hydrodynamic diffusivity coefficients in the gradient and vorticity directions were also obtained and were found to be two orders of magnitude larger than the corresponding values in fiber suspensions at the same particle concentrations. Rotary Brownian dynamics simulations of isolated Brownian rings were used to understand the shear rate dependence of suspension rheology. The orbit distribution observed in the regime of weak Brownian motion, P e ≫ ϕ T − 3, was surprisingly similar to that obtained from pairwise interaction calculations of non-Brownian rings. Here, the Peclet number P e is the ratio of the shear rate and the rotary diffusivity of the particle and ϕ T is the effective inverse-aspect ratio of the particle (approximately equal to 2 π times the inverse of its non-dimensional Jeffery time period). Thus, the rheology results obtained from pairwise interactions should retain accuracy even for weakly Brownian rings ( n ≪ 1 and ϕ T − 3 ≪ P e ). 
    more » « less
  3. ; ; ; ( , Journal of Fluid Mechanics)
    High Reynolds number wall-bounded turbulent flows subject to buoyancy forces are fraught with complex dynamics originating from the interplay between shear generation of turbulence ( $S$ ) and its production or destruction by density gradients ( $B$ ). For horizontal walls, $S$ augments the energy budget of the streamwise fluctuations, while $B$ influences the energy contained in the vertical fluctuations. Yet, return to isotropy remains a tendency of such flows where pressure–strain interaction redistributes turbulent energy among all three velocity components and thus limits, but cannot fully eliminate, the anisotropy of the velocity fluctuations. A reduced model of this energy redistribution in the inertial (logarithmic) sublayer, with no tuneable constants, is introduced and tested against large eddy and direct numerical simulations under both stable ( $B<0$ ) and unstable ( $B>0$ ) conditions. The model links key transitions in turbulence statistics with flux Richardson number (at $Ri_{f}=-B/S\approx$ $-2$ , $-1$ and $-0.5$ ) to shifts in the direction of energy redistribution. Furthermore, when coupled to a linear Rotta-type closure, an extended version of the model can predict individual variance components, as well as the degree of turbulence anisotropy. The extended model indicates a regime transition under stable conditions when $Ri_{f}$ approaches $Ri_{f,max}\approx +0.21$ . Buoyant destruction $B$ increases with increasing stabilizing density gradients when $Ri_{f} more » « less
  4. ; ; ; ( , Journal of fluids engineering)
    null (Ed.)
    The Marangoni propulsion of spheres and elliptical disks floating on the air–water interface were studied to understand the effect of particle shape on its motion and its stability at moderate Reynolds numbers. Self-propulsion of the Marangoni surfer was achieved by coating half of the spheres and the elliptical disks with either a solution of soap or isopropyl alcohol (IPA). The presence of the soap or IPA resulted in a surface tension gradient across the particles which propelled the particles in the direction of increasing surface tension. Beyond a critical velocity, a transition was observed from a straight-line motion to a rotational motion. These vortices were observed to shed above a critical Reynolds number resulting in an unbalanced torque that caused the particles to rotate. Increasing the aspect ratio between the major and minor axes of the elliptical disks was found to decrease their stability and greatly enhance their rate of rotation. This was especially true for elliptical disks traveling in a direction parallel to their major axis. The interactions between the particles and the wall of a Petri dish were also studied. Repulsive, concave curvature was found to decrease stability and enhance rotational motion, while attractive, convex curvature was shown to stabilize the straight-line motion of the spheres. For the neutrally buoyant elliptical disks, the presence of the bounding wall was found to greatly stabilize the straight-line motion of the elliptical disks when they were traveling in a direction parallel to their minor axis. 
    more » « less
  5. ; ( , Journal of Fluids Engineering)
    null (Ed.)
    Abstract The performance of different versions of the discrete random walk models in turbulent flows with nonuniform normal root-mean-square (RMS) velocity fluctuations and turbulence time scales were carefully investigated. The OpenFOAM v2−f low Reynolds number turbulence model was used for evaluating the fully developed streamwise velocity and the wall-normal RMS velocity fluctuations profiles in a turbulent channel flow. The results were then used in an in-house matlab particle tracking code, including the drag and Brownian forces, and the trajectories of randomly injected point-particles with diameters ranging from 10 nm to 30 μm were evaluated under the one-way coupling assumption. The distributions and deposition velocities of fluid-tracer and finite-size particles were evaluated using the conventional-discrete random walk (DRW) model, the modified-DRW model including the velocity gradient drift correction, and the new improved-DRW model including the velocity and time gradient drift terms. It was shown that the conventional-DRW model leads to superfluous migration of fluid-point particles toward the wall and erroneous particle deposition rate. The concentration profiles of tracer particles obtained by using the modified-DRW model still are not uniform. However, it was shown that the new improved-DRW model with the velocity and time scale drift corrections leads to uniform distributions for fluid-point particles and reasonable concentration profiles for finite-size heavy particles. In addition, good agreement was found between the estimated deposition velocities of different size particles by the new improved-DRW model with the available data. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email