skip to main content


Title: Effect of Surface Topography on Particle Deposition from Liquid Suspensions in Channel Flow
A Eulerian—Lagrangian model has been developed to simulate particle attachment to surfaces with arc-shaped ribs in a two-dimensional channel flow at low Reynolds numbers. Numerical simulation has been performed to improve the quantitative understanding of how rib geometries enhance shear rates and particle-surface interact for various particle sizes and flow velocities. The enhanced shear rate is attributed to the wavy flows that develop over the ribbed surface and the weak vortices that form between adjacent ribs. Varying pitch-to-height ratio can alter the amplitude of the wavy flow and the angle of attack of the fluid on the ribs. In the presence of these two competing factors, the rib geometry with a pitch-to-height ratio of two demonstrates the greatest shear rate and the lowest fraction of particle attachment. However, the ribbed surfaces have negligible effects on small particles at low velocities. A force analysis identifies a threshold shear rate to reduce particle attachment. The simulated particle distributions over the ribbed surfaces are highly non-uniform for larger particles at higher velocities. The understanding of the effect of surface topography on particle attachment will benefit the design of surface textures for mitigating particulate fouling in a wide range of applications.  more » « less
Award ID(s):
1705538
NSF-PAR ID:
10319164
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Fluids
Volume:
5
Issue:
1
ISSN:
2311-5521
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Spiral ribs are among the most common morphological features in mollusk shells and previous studies have shown them to have functional significance with expected evolutionary consequences. Many previous studies, however, have treated these features as potentially analogous across taxa, without examining whether they may have important constructional dissimilarities. Mollusk shells are made of multiple layers of calcite or aragonite which may exhibit different microstructure or microstructure orientations which may in turn impact their mechanical properties. In this study, five specimens of marine mollusks with spiral ribs, including three turritellid gastropods and two bivalves, were examined under SEM to examine microstructure of ribbed region in comparison of non-ribbed region. SEM imaging revealed differences in the number and thickness of distinct microstructural layers of each shell and allowed comparisons to be made between the ribbed and non-ribbed region of each specimen, providing a greater understanding of how these ribs were constructed during shell deposition. Ribs in all specimens are formed through the thickening of single or multiple crossed-lamellar layers, but differences in rib ultrastructures were found among species and different ribs of same species, showing great diversity and complexity of constructional mechanisms. This diversity in rib construction might indicate heterology in the development of shell sculpture, especially mechanisms for differences in concurrently deposited rib strength. This is especially notable for turritellids where the pattern of onset of spiral ornamentation is phylogenetically informative, suggesting homology of rib identity. Further study will be conducted on turritellid gastropods in different lineages to explore the taxonomic meaning of different rib constructional mechanisms. 
    more » « less
  2. Using a numerical model, we analyse the effects of shape on both the orientation and transport of anisotropic particles in wavy flows. The particles are idealized as prolate and oblate spheroids, and we consider the regime of small Stokes and particle Reynolds numbers. We find that the particles preferentially align into the shear plane with a mean orientation that is solely a function of their aspect ratio. This alignment, however, differs from the Jeffery orbits that occur in the residual shear flow (that is, the Stokes drift velocity field) in the absence of waves. Since the drag on an anisotropic particle depends on its alignment with the flow, this preferred orientation determines the effective drag on the particles, which in turn impacts their net downstream transport. We also find that the rate of alignment of the particles is not constant and depends strongly on their initial orientation; thus, variations in initial particle orientation result in dispersion of anisotropic-particle plumes. We show that this dispersion is a function of the particle’s eccentricity and the ratio of the settling and wave time scales. Due to this preferential alignment, we find that a plume of anisotropic particles in waves is on average transported farther but dispersed less than it would be if the particles were randomly oriented. Our results demonstrate that accurate prediction of the transport of anisotropic particles in wavy environments, such as microplastic particles in the ocean, requires the consideration of these preferential alignment effects. 
    more » « less
  3. Dense suspensions of particles in viscous liquid often demonstrate the striking phenomenon of abrupt shear thickening, where their viscosity increases strongly with increase of the imposed stress or shear rate. In this work, discrete-particle simulations accounting for short-range hydrodynamic, repulsive, and contact forces are performed to simulate flow of shear thickening bidisperse suspensions, with the packing parameters of large-to-small particle radius ratio δ = 3 and large particle fraction ζ = 0.15, 0.50, and 0.85. The simulations are carried out for volume fractions 0.54 ≤ ϕ ≤ 0.60 and a wide range of shear stresses. The repulsive forces, of magnitude F R , model the effects of surface charge and electric double-layer overlap, and result in shear thinning at small stress, with shear thickening beginning at stresses σ ∼ F R a −2 . A crossover scaling analysis used to describe systems with more than one thermodynamic critical point has recently been shown to successfully describe the experimentally-observed shear thickening behavior in suspensions. The scaling theory is tested here on simulated shear thickening data of the bidisperse mixtures, and also on nearly monodisperse suspensions with δ = 1.4 and ζ = 0.50. Presenting the viscosity in terms of a universal crossover scaling function between the frictionless and frictional maximum packing fractions collapses the viscosity for most of the suspensions studied. Two scaling regimes having different exponents are observed. The scaling analysis shows that the second normal stress difference N 2 and the particle pressure Π also collapse on their respective curves, with the latter featuring a different exponent from the viscosity and normal stress difference. The influence of the fraction of frictional contacts, one of the parameters of the scaling analysis, and its dependence on the packing parameters are also presented. 
    more » « less
  4. Understanding the orientation dynamics of anisotropic colloidal particles is important for suspension rheology and particle self-assembly. However, even for the simplest case of dilute suspensions in shear flow, the orientation dynamics of non-spherical Brownian particles are poorly understood. Here we analytically calculate the time-dependent orientation distributions for non-spherical axisymmetric particles confined to rotate in the flow–gradient plane, in the limit of small but non-zero Brownian diffusivity. For continuous shear, despite the complicated dynamics arising from the particle rotations, we find a coordinate change that maps the orientation dynamics to a diffusion equation with a remarkably simple ratio of the enhanced rotary diffusivity to the zero shear diffusion: $D_{eff}^{r}/D_{0}^{r}=(3/8)(p-1/p)^{2}+1$ , where $p$ is the particle aspect ratio. For oscillatory shear, the enhanced diffusion becomes orientation dependent and drastically alters the long-time orientation distributions. We describe a general method for solving the time-dependent oscillatory shear distributions and finding the effective diffusion constant. As an illustration, we use this method to solve for the diffusion and distributions in the case of triangle-wave oscillatory shear and find that they depend strongly on the strain amplitude and particle aspect ratio. These results provide new insight into the time-dependent rheology of suspensions of anisotropic particles. For continuous shear, we find two distinct diffusive time scales in the rheology that scale separately with aspect ratio $p$ , as $1/D_{0}^{r}p^{4}$ and as $1/D_{0}^{r}p^{2}$ for $p\gg 1$ . For oscillatory shear flows, the intrinsic viscosity oscillates with the strain amplitude. Finally, we show the relevance of our results to real suspensions in which particles can rotate freely. Collectively, the interplay between shear-induced rotations and diffusion has rich structure and strong effects: for a particle with aspect ratio 10, the oscillatory shear intrinsic viscosity varies by a factor of ${\approx}2$ and the rotational diffusion by a factor of ${\approx}40$ . 
    more » « less
  5. 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