Simulation and Visualization of Flows Laden with Cylindrical Nanoparticles in a Mixing Layer

/e motion of cylindrical particles in a mixing layer is studied using the pseudospectral method and discrete particle model. /e effect of the Stokes number and particle aspect ratio on the mixing and orientation distribution of cylindrical particles is analyzed. /e results show that the rollup of mixing layer drives the particles to the edge of the vortex by centrifugal force. /e cylindrical particles with the small Stokes number almost follow fluid streamlines and are mixed thoroughly, while those with the large Stokes number, centrifugalized and accumulated at the edge of the vortex, are poorly mixed. /e mixing degree of particles becomes worse as the particle aspect ratio increases. /e cylindrical particles would change their orientation under two torques and rotate around their axis of revolution aligned to the vorticity direction when the shear rate is low, while aligning on the flow-gradient plane beyond a critical shear rate value. More particles are oriented with the flow direction, and this phenomenon becomes more obvious with the decrease of the Stokes number and particle aspect ratio.


Introduction
Mixing of nanoparticles in a flow was found in many applications and has been a subject of interest during the recent twenty years.
e mixing of nanoparticles affects the property of the flow and the behaviors of the final production.
erefore, understanding deeply such a flow is important to design a novel technology for nanoparticle manipulation.
Some efforts have been put into the numerical simulation and visualization of nanoparticles.Singh et al. [1] deposited In 2 O 3 and In 2 O 3 : Ag nanoparticle layers and indicated the presence of Ag 2 O and Ag in air-and vacuumannealed samples, respectively.Aruna et al. [2] reported the modifications in CO sensing of SnOx nanoparticle layers by using Pd nanoparticles; they found that the homogeneously mixed nanoparticle layers show capability between CO and ethanol as a manifestation of the dual conductance response.Hsiao et al. [3] used mixed PSS-DTMA Langmuir layers to incorporate with silver precursors from the subphase and transferred them onto mica substrates; they showed that enhancing the DTMA (+) in the mixed PSS-DTMA system would increase the hydrophobic property of the complexes.Finally, they inferred that the polyelectrolyte-surfactant template can offer a potential of designing structures of polyelectrolyte-nanoparticle materials.Liggieri et al. [4] utilized different techniques to measure the dilational viscoelasticity in a wide frequency range and used these data to provide qualitative and quantitative information about structural behaviors of complex mixed layers.Xie et al. [5] simulated numerically the impact coherent structure on the Brownian coagulation of particles in a mixing layer, and it is found that the number density of nanoparticles decreases gradually as the flow evolves, while the particle average volume increases.e impact of fluid advection on particle coagulation is small in the regions far away from the eddy structure.
e particle coagulation in the eddy core has a wave-like distribution.Yazhgur et al. [6] studied the adsorption films of silica nanoparticles modified by a cationic surfactant at the air-water interface.
ey divided the whole surfactant concentration range into four regions characterized by di erent surface rheological behaviors.Kim et al. [7] studied the optical, electrical, and morphological behaviors of the Ga 2 O 3 NP/SWNT layers by enhancing the thickness of SWNTs.Mur eld and Garrick [8] performed numerical simulations of nanoparticle nucleation in turbulent wakes and showed that nucleation initially occurs in the shear layers where molecular transport dominates and across the span of the wake.Fager and Garrick [9] presented the results of direct numerical simulation of zinc nanoparticle nucleation in a turbulent round jet; they indicated that particle nucleation occurs along the outer region of the jet, and the regions over which nucleation occurs increase signi cantly after collapse of the jet potential core.Guzman et al. [10] carried out measurements of the dilational viscoelastic modulus against the frequency by the oscillatory barrier method at di erent degrees of compression of the monolayer in order to deepen the understanding of the impacts of nanoparticles on the interfacial properties of biosystems.Orsi et al. [11] studied the interfacial dynamics of a 2D self-organized mixed layer consisted of silica nanoparticles interacting with phospholipid monolayers at the air-water interface and found a dynamical transition from the Brownian di usion to an arrested state.Garrick [12] showed that, in the proximal region of the jet, condensation is the dominant mechanism; however, once the jet potential core collapses and turbulent mixing begins, coagulation is the dominant mechanism.Kerli and Alver [13] investigated the mixture of ZnO and NiO e ect on the solar cell and observed that the solar cells made with ZnO have the highest performance with the e ciency of 0.542%.e researches available in the literature are mainly concerned with spherical particles as shown above.However, many particles in the practical usage are nonspherical among which the cylinder particle is a typical example.Compared with the spherical particle, the mixing of cylindrical particles in a ow is more complicated because the change of particle orientation is coupled with the translation motion.Understanding the property of cylindrical particles in the mixing layer is a great importance in the optimization of production processes.As far as we know, the motion and mixing behaviors of cylindrical nanoparticles in the mixing layer have not been reported in the literatures yet.erefore, the objective of the present paper is to explore the impact of the Stokes number and particle aspect ratio on the mixing property and orientation distributions of cylindrical particles in a mixing layer.

Mixing Layer Flow and Fluid Equation.
A mixing layer consisting of two streams and a cylindrical particle are shown in Figure 1 where velocities U 1 and U 2 are di erent.In the mixing layer, the hydrodynamic instabilities will happen, which leads to the ow to roll up and form coherent vortex and even pairing of the vortices.
Dilute cylindrical particle suspension is studied here, that is, the e ect of the cylindrical particle on the uid can be ignored; thus, the governing equations are the following continuity and momentum equations: where u and p are the uid velocity and pressure, respectively; ω is the vorticity; and Re is the Reynolds number.

Equation for Cylindrical Particles.
e size of a non-Brownian and rigid cylindrical particle is smaller than the characteristic scale of the ow eld.e motion of particles is caused by force F and torque T exerted by the uid: where m p , u p , and Ω are the particle mass, velocity, and angular velocity, respectively; I is the particle moment of inertia; f(s) is the force distribution along the particle axis; s is the distance from any point to the particle center; l is the half length of the particle; and p is the orientation vector of the particle.Equation ( 3) can be nondimensionalized with the scale of the ow eld as du p dt where a is the particle aspect ratio, St is the Stokes number (St ρ p r 2 U 0 /2μθ 0 , in which ρ p is the particle density, r is the particle radius, U 0 U 1 − U 2 , and θ 0 is the initial momentum thickness of the mixing layer), and ρ * is the particle-to-uid density ratio.(1) e velocity profile with hyperbolic-tangent form of the base flow: (2) e lower wavenumber disturbances with stream function: where ψ is the stream function; A 1 and A 2 are the amplitude of the fundamental and its subharmonic modes, respectively, and we set A 1 � 0.1 and A 2 � 0 and A 1 � 0.1 and A 2 � 0.06 for the process of vortex rollup and vortex pairing; Ɽ means the real part for a complex; φ 1 (y) and φ 2 (y) are the eigenfunctions for the fundamental and subharmonic modes, respectively; and α is the fundamental wavenumber and is set to be 0.4446 [14].

Boundary Condition.
e periodic boundary conditions in the streamwise (x) and transverse (y) directions are imposed by introducing an image flow far enough from the mixing layer center because all the perturbations vanish as y → ∞.
erefore, the Fourier spectral method can be exerted in both directions, and fast Fourier transformation can be used.e periodic L 1 and L 2 in the streamwise and transverse directions are, respectively, taken to be 2π/r and 28 for the vortex rollup, and both are doubled for the vortex pairing. 1) and ( 2) are numerically solved using the pseudospectral method.e Adams-Bashforth scheme is used for the nonlinear term, and the implicit Crank-Nicolson scheme is used for other terms.e time step is set to be 0.03.e collection points in the streamwise and transverse directions are taken to be 128 and 256, respectively.

Computational Approach. Equations (
Equations ( 4) and ( 5) are solved with the fourth-order Runge-Kutta method for 400 cylindrical particles.ese particles are distributed and oriented homogeneously in the upper layer of the flow as shown in Figure 1. e time step of the integration is small enough to guarantee a stable trajectory of the particles.

Visualization of Vortex Rolling Up and Vortex
Pairing.
e Kelvin-Helmholtz instability makes the mixing layer roll up and forms the large-scale vortex as shown in Figure 2, where the center of the vortex is called the center point and the middle point of two center points is called the saddle point.
e flow can be divided into the inhomogeneous region and the homogeneous region.e streamlines connecting two saddle points are defined as the demarcation lines; the region encompassed by the streamlines, that is, vortex core, is considered as the inhomogeneous region, and the other region is called the homogeneous region.
e process of vortex pairing is shown in Figure 3 where we can see that two vortices rotate each other and merge into one vortex finally.

Particle Distribution.
e mixing of cylindrical particles can be described by the particle spatial distributions.
e rollup of the mixing layer drives the particles to the edge of the vortex by the inertial centrifugal force, and at the same time it brings about an asymmetric shear rate distribution around the particle and generates a normal stress imbalance on the particle surface.
e cylindrical particle migration and mixing are dependent on the above two factors.
e mixing degree of the particles is negatively correlated with the inhomogeneous degree of particle distribution.
erefore, the inhomogeneous degree of particle distribution M n is used in order to quantitatively describe the mixing behavior: where (x pi , y pi ) and (x c , y c ) are the coordinates of the center of the ith particle and vortex center, respectively.A larger value of M n is corresponding to a poorer mixing.

Effect of the Stokes Number on the Mixing
Degree. Figure 4 shows the particle distribution at t � 100 for St � 0.02 and 1.
Stokes number (St) is directly related to the function of particle inertia.In the process of rollup of the mixing layer, the cylindrical particles for St � 0.02 almost follow fluid streamlines and distribute more homogeneously and are mixed thoroughly.On the contrary, the particles for St � 1.0 are centrifugalized from the vortex core to the edge by inertia and are poorly mixed.
Figure 5 shows the change of inhomogeneous degree with time at different Stokes numbers.We can see well mixed for St � 0.02.e mixing degree becomes worse with the increase of the Stokes number because more particles are centrifugalized to the edge of the vortex when St is larger.

Effect of the Particle Aspect Ratio on the Mixing
Degree. Figure 6 shows the change of inhomogeneous degree with time for different particle aspect ratios.e mixing degree becomes worse as the particle aspect ratio increases because the angular acceleration of the particles is inversely proportional to the particle aspect ratio.Particle angular acceleration and translational motion induced by the fluid are smaller for particles with larger particle aspect ratios so that the centrifugal force drives the particles towards the edge of the vortex.

Orientation Distribution of the Cylindrical Particle.
e cylindrical particles will change the orientation angles φ and θ (shown in Figure 1) under two torques.One torque makes the particle rotate around the vorticity axis, and another torque leads to the particle spin around the flow direction.For describing the orientation distribution, the orientation distribution function ψ is used: where θ is shown in Figure 1, N is the total particle number, and N i+1−i is the number of particles with the angles ranging from i to i + 1.

Effect of the Stokes Number on the Particle
Orientation. Figure 7 shows the orientation distribution function of particles for different St.We can see that, with the decrease of the Stokes number, more particles are oriented with the flow direction.is may be attributed that the angular acceleration of the particles is inversely proportional to the Stokes number, and the angular acceleration induced by the fluid stress is larger at smaller Stokes number so that the particles change their orientation more rapidly.

Effect of the Particle Aspect Ratio on the Particle
Orientation.e orientation distribution functions of cylindrical particles for different particle aspect ratios are shown in Figure 8. e particle would rotate around its axis of revolution aligned to the vorticity direction when the shear rate is low, while aligning on the flow-gradient plane beyond a critical shear rate value.From the figure, we can see that more particles are oriented on the flow-gradient plane (θ � 0 °) and more particles are oriented towards the flow direction with the decrease of the particle aspect ratio.

Conclusion
e effect of the Stokes number and particle aspect ratio on the mixing and orientation distribution of cylindrical particles in a mixing layer is numerically studied.e mixing of cylindrical particles can be described by the particle spatial distributions.
e rollup of the mixing layer drives the particles to the edge of the vortex by the inertial centrifugal force and at the same time brings about an asymmetric shear rate distribution around the particle.e cylindrical particles with the small Stokes number almost follow fluid streamlines and are mixed thoroughly.On the contrary, the particles with the large Stokes number are centrifugalized to the edge of the vortex and are poorly mixed.e mixing degree of particles becomes worse as the particle aspect ratio increases.Particle angular acceleration and translational motion are smaller for that particle with a larger particle aspect ratio so that the centrifugal force drives the particles towards the edge of the vortex.e cylindrical particles will change the orientation angle under the two torques with which the particles would rotate around its axis of revolution aligned to the vorticity direction when the shear rate is low, while aligning on the flow-gradient plane beyond a critical shear rate value.More particles are oriented with the flow direction, and this phenomenon becomes more obvious with the decrease of the Stokes number and particle aspect ratio.

Figure 6 :
Figure 6: Change of inhomogeneous degree with time at different a (Re � 200 and St � 0.1).