Hostname: page-component-586b7cd67f-t8hqh Total loading time: 0 Render date: 2024-11-26T08:41:26.308Z Has data issue: false hasContentIssue false

Disks aligned in a turbulent channel

Published online by Cambridge University Press:  28 April 2015

Greg A. Voth*
Affiliation:
Department of Physics, Wesleyan University, Middletown, CT 06459, USA
*
Email address for correspondence: [email protected]

Abstract

Anisotropic particles are suspended in a wide range of industrial, environmental and biological fluid flows. The orientations of these particles are sometimes randomized by turbulence, but often they are brought into preferential alignment by the fluid flow. In a recently published study, Challabotla, Zhao & Andersson (J. Fluid Mech., vol. 766, 2015, R2) performed the first numerical simulations of inertial disks in a turbulent channel flow. They find that disks can be made to preferentially align either parallel or perpendicular to the wall depending on the particle density. Particle shape also affects alignment, particularly for lower density particles, and the alignment of disks is quite different from the alignment of fibres.

Type
Focus on Fluids
Copyright
© 2015 Cambridge University Press 

1. Introduction

We see anisotropic particles carried by turbulent fluid flows all around us. The particles could be snowflakes or leaves blowing in the wind. They could be ice crystals in clouds that produce sundogs or seeds with shapes adapted by evolution to enhance their dispersal. Non-spherical particles also play important but less noticed roles in our lives. Paper is made using turbulent suspensions of wood fibres (Lundell, Söderberg & Alfredsson Reference Lundell, Söderberg and Alfredsson2011); the lowest trophic levels of the largest ecosystem on the planet are dominated by plankton of a wide variety of shapes (Guasto, Rusconi & Stocker Reference Guasto, Rusconi and Stocker2012; Byron et al. Reference Byron, Einarsson, Gustavsson, Voth, Mehlig and Variano2015); and the disk shape and deformability of red blood cells play a critical role in the human circulatory system (Chien Reference Chien1987).

A natural expectation might be that turbulence will randomize the orientation of these particles so that simple rotationally averaged models can describe light scattering by particles, rheological properties of a suspension or the rotation rates of particles. But even in isotropic turbulence, it is found that particles become aligned by fluid velocity gradients and this alignment results in particle rotation rates up to three times smaller than that predicted by rotationally averaged models (Shin & Koch Reference Shin and Koch2005; Parsa et al. Reference Parsa, Calzavarini, Toschi and Voth2012).

Particles occur in a wide variety of shapes, but many of them can be accurately modelled by effective ellipsoids (Bretherton Reference Bretherton1962). An axisymmetric ellipsoid, (or spheroid) can be specified by a single aspect ratio, ${\it\lambda}$ , defined as the ratio of the dimension along the symmetry axis to a perpendicular dimension: ${\it\lambda}=1$ is a sphere, ${\it\lambda}>1$ is a prolate spheroid or fibre and ${\it\lambda}<1$ is an oblate spheroid or disk. Fibres have received most of the attention from researchers working on anisotropic particles in turbulence. The reason for this prejudice against small aspect ratios is not entirely clear, but factors involve the industrial importance of fibre flows, the theoretical convenience of slender body theory, and the experimental convenience of imaging nearly one-dimensional objects.

Disks have recently climbed into the limelight as researchers realized that the rotation and alignment of disk-shaped particles provide a challenging test case for stochastic models of the velocity gradient tensor (Chevillard & Meneveau Reference Chevillard and Meneveau2013) as well as for analytic models of particle rotation rates (Gustavsson, Einarsson & Mehlig Reference Gustavsson, Einarsson and Mehlig2014). Experiments have recently been able to measure rotations of three-dimensional printed particles that rotate like disks and whose orientation can be accurately measured with multiple cameras (Marcus et al. Reference Marcus, Parsa, Kramel, Ni and Voth2014).

2. Overview

Recently, Challabotla, Zhao, & Andersson have presented a groundbreaking numerical study that explores the rotation and alignment of inertial disk-shaped particles in a turbulent channel flow (Challabotla, Zhao & Andersson Reference Challabotla, Zhao and Andersson2015). This has been a very fruitful system in the study of fibres in flows (Mortensen et al. Reference Mortensen, Andersson, Gillissen and Boersma2008; Marchioli, Fantoni & Soldati Reference Marchioli, Fantoni and Soldati2010), but the dynamics of disks in wall-bounded shear flow has not previously been explored. Challabotla et al. (Reference Challabotla, Zhao and Andersson2015) use a point particle model similar to those that have been widely used for spherical particles in turbulent flows (Balachandar & Eaton Reference Balachandar and Eaton2010). When spherical particles are sufficiently small, their equation of motion can be accurately approximated with a simple Stokes drag force. The equivalent model for ellipsoidal particles requires a resistance tensor to find the drag force vector from the particle orientation and its relative velocity with the fluid. Modelling the orientation requires integrating the torques on the particle produced by the fluid velocity gradients. Challabotla et al. (Reference Challabotla, Zhao and Andersson2015) provide a clear and compact guide for implementing this anisotropic Stokes drag model. This model does not include fluid inertia effects, but is a good approximation in many systems of interest where particles are small and their density is large compared with the fluid density.

Figure 1. Disks in a turbulent channel flow showing the preferential orientation near the wall. Data are the same as Challabotla et al. (Reference Challabotla, Zhao and Andersson2015, figure 1(b)). The mean velocity is in the $x$ (streamwise) direction, the flow has periodic boundary conditions in the $y$  direction, and no-slip boundary conditions at the upper and lower walls ( $z$  direction). The particles have aspect ratio  $1/3$ and Stokes number 30. Particle colour encodes the spanwise component of the particle angular velocity,  ${\it\omega}_{y}$ .

A snapshot of disks moving in the channel flow simulated by Challabotla et al. (Reference Challabotla, Zhao and Andersson2015) is shown in figure 1. The most striking feature is the preferential orientation of the disks near the walls where the disks’ symmetry axes align in the spanwise direction. In the centre of the channel, particles are more randomly oriented. This image is for heavy particles whose response time, ${\it\tau}_{p}$ , is much longer than the fluid viscous timescale, ${\it\tau}_{{\it\nu}}$ ; so the Stokes number is relatively large, $\mathit{St}={\it\tau}_{p}/{\it\tau}_{{\it\nu}}=30$ . At lower $\mathit{St}=1$ , they observe a very different alignment with the symmetry axis of the disks in the wall normal direction. None of their simulations find disks preferentially aligned in the streamwise direction which is the orientation chosen by fibres (Mortensen et al. Reference Mortensen, Andersson, Gillissen and Boersma2008; Marchioli et al. Reference Marchioli, Fantoni and Soldati2010).

The interaction of preferentially oriented particles with the mean and fluctuating velocity gradients produces fascinating rotation statistics. Heavy particles ( $\mathit{St}=30$ ) have a mean rotation that is almost equal to the mean fluid rotation, while light particles have a smaller rotation that is strongly dependent on particle shape. The fluctuating rotation rates show even more complex behaviour which reveals an important feature of anisotropic inertial particles: the rotational response time of particles can be quite different than their translational response time.

To understand these observations conceptually, the authors turn to the large body of work on the motion of anisotropic particles in simple shear flow where neutrally buoyant spheroidal particles undergo periodic motions called Jeffery orbits. The preferential alignment of the low-Stokes-number particles in the wall-normal direction can be understood simply by considering the Jeffery orbits of a disk in simple shear flow. Disk-shaped particles in orbits that bring their symmetry axis near the wall-normal direction will spend a long time in this orientation before tumbling, which leads to the observed orientation. The alignment of the high-Stokes-number particles is controlled by particle inertia. The alignment seen near the walls in figure 1 matches the steady state alignment observed in earlier work on disks in simple shear, where the symmetry axis aligns parallel to the flow vorticity. Turbulence competes with alignment by the mean shear so this preferential alignment diminishes toward the centre of the channel where the mean shear becomes small.

3. Future

These important results on the alignment and rotation of inertial disks in turbulent shear flow suggest several promising research directions. The authors observe that both rods and disks have a preferred alignment with a long axis in the streamwise direction. In isotropic turbulence, it has been found that rods and disks align with a long axis in the direction of maximum fluid stretching defined by the Cauchy–Green strain tensors (Marcus et al. Reference Marcus, Parsa, Kramel, Ni and Voth2014; Ni, Ouellette & Voth Reference Ni, Ouellette and Voth2014). Consideration of the mean and fluctuating stretching fields may provide additional insights into particle alignment.

Another interesting subject for future work is the instantaneous spatial distribution of particle orientations. For neutrally buoyant particles, the instantaneous orientation should be a nearly continuous field whose rich topology has begun to be explored (Szeri Reference Szeri1993; Wilkinson, Bezuglyy & Mehlig Reference Wilkinson, Bezuglyy and Mehlig2011). The addition of particle inertia introduces the possibility of caustics in the orientation field (Gustavsson et al. Reference Gustavsson, Einarsson and Mehlig2014; Siewert, Kunnen & Schröder Reference Siewert, Kunnen and Schröder2014b ), which can strongly affect collision probabilities.

Attempts to reproduce the gravity-free numerical results from Challabotla et al. (Reference Challabotla, Zhao and Andersson2015) in the laboratory will face the challenge that even quite small inertial disks show effects from gravitational sedimentation. But this raises opportunities as well. The behaviour of sedimenting anisotropic particles is significantly richer than the gravity-free case including phenomena such as sedimentation instabilities (Guazzelli & Hinch Reference Guazzelli and Hinch2011) and alignment of particles by inertial torques (Siewert et al. Reference Siewert, Kunnen, Meinke and Schröder2014a ).

Some key fundamental phenomena in the behaviour of anisotropic particles in turbulent flows have now been identified, but much of the parameter space remains to be explored. Some of the most pressing applications involve dense suspensions, sedimentation and deformability of anisotropic particles in turbulence. Our deepening understanding of foundational cases such as rigid disks in channel flow is essential for the development of predictive models of the complex turbulent suspensions all around us.

References

Balachandar, S. & Eaton, J. K. 2010 Turbulent dispersed multiphase flow. Annu. Rev. Fluid Mech. 42, 111133.CrossRefGoogle Scholar
Bretherton, F. P. 1962 The motion of rigid particles in a shear flow at low Reynolds number. J. Fluid Mech. 14 (2), 284304.CrossRefGoogle Scholar
Byron, M., Einarsson, J., Gustavsson, K., Voth, G., Mehlig, B. & Variano, E. 2015 Shape-dependence of particle rotation in isotropic turbulence. Phys. Fluids 27 (3), 035101.Google Scholar
Challabotla, N. R., Zhao, L. & Andersson, H. I. 2015 Orientation and rotation of inertial disk particles in wall turbulence. J. Fluid Mech. 766, R2.Google Scholar
Chevillard, L. & Meneveau, C. 2013 Orientation dynamics of small, triaxialellipsoidal particles in isotropic turbulence. J. Fluid Mech. 737, 571596.Google Scholar
Chien, S. 1987 Red-cell deformability and its relevance to blood flow. Annu. Rev. Physiol. 49, 177192.Google Scholar
Guasto, J. S., Rusconi, R. & Stocker, R. 2012 Fluid mechanics of planktonic microorganisms. Annu. Rev. Fluid Mech. 44 (1), 373400.Google Scholar
Guazzelli, É. & Hinch, J. 2011 Fluctuations and instability in sedimentation. Annu. Rev. Fluid Mech. 43, 97116.Google Scholar
Gustavsson, K., Einarsson, J. & Mehlig, B. 2014 Tumbling of small axisymmetric particles in random and turbulent flows. Phys. Rev. Lett. 112, 014501.CrossRefGoogle ScholarPubMed
Lundell, F., Söderberg, L. D. & Alfredsson, P. H. 2011 Fluid mechanics of papermaking. Annu. Rev. Fluid Mech. 43, 195217.Google Scholar
Marchioli, C., Fantoni, M. & Soldati, A. 2010 Orientation, distribution, and deposition of elongated, inertial fibers in turbulent channel flow. Phys. Fluids 22 (3), 033301.Google Scholar
Marcus, G. G., Parsa, S., Kramel, S., Ni, R. & Voth, G. A. 2014 Measurements of the solid-body rotation of anisotropic particles in 3D turbulence. New J. Phys. 16 (10), 102001.CrossRefGoogle Scholar
Mortensen, P. H., Andersson, H. I., Gillissen, J. J. J. & Boersma, B. J. 2008 Dynamics of prolate ellipsoidal particles in a turbulent channel flow. Phys. Fluids 20 (9), 093302.CrossRefGoogle Scholar
Ni, R., Ouellette, N. T. & Voth, G. A. 2014 Alignment of vorticity and rods with Lagrangian fluid stretching in turbulence. J. Fluid Mech. 743, R3.Google Scholar
Parsa, S., Calzavarini, E., Toschi, F. & Voth, G. A. 2012 Rotation rate of rods in turbulent fluid flow. Phys. Rev. Lett. 109 (13), 134501.Google Scholar
Shin, M. & Koch, D. L. 2005 Rotational and translational dispersion of fibres in isotropic turbulent flows. J. Fluid Mech. 540, 143173.Google Scholar
Siewert, C., Kunnen, R. P. J., Meinke, M. & Schröder, W. 2014a Orientation statistics and settling velocity of ellipsoids in decaying turbulence. Atmos. Res. 142, 4556.Google Scholar
Siewert, C., Kunnen, R. P. J. & Schröder, W. 2014b Collision rates of small ellipsoids settling in turbulence. J. Fluid Mech. 758, 686701.CrossRefGoogle Scholar
Szeri, A. J. 1993 Pattern-formation in recirculating-flows of suspensions of orientable particles. Phil. Trans. R. Soc. Lond. Ser. 345 (1677), 477506.Google Scholar
Wilkinson, M., Bezuglyy, V. & Mehlig, Bernhard 2011 Emergent order in rheoscopic swirls. J. Fluid Mech. 667, 158187.Google Scholar
Figure 0

Figure 1. Disks in a turbulent channel flow showing the preferential orientation near the wall. Data are the same as Challabotla et al. (2015, figure 1(b)). The mean velocity is in the $x$ (streamwise) direction, the flow has periodic boundary conditions in the $y$ direction, and no-slip boundary conditions at the upper and lower walls ($z$ direction). The particles have aspect ratio $1/3$ and Stokes number 30. Particle colour encodes the spanwise component of the particle angular velocity, ${\it\omega}_{y}$.