You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
138<br />
Dynamics of heavy particles near a helical vortex lament<br />
R. H.A. IJzermans ∗ , R. Hagmeijer ∗ and P.J. van Langen ∗<br />
The motion of small heavy particles near a helical vortex filament is investigated<br />
both numerically and analytically. We study the case of a helical vortex filament<br />
in unbounded space and the case of a helical vortex filament which is placed in a<br />
concentric cylindrical pipe. Potential flow is assumed which is not influenced by the<br />
presence of the heavy particles (one-way coupling). For both cases, the velocity field<br />
in a helical coordinate frame is described by a stream function 1 , 2 . In the equation<br />
of motion of the heavy particles, Stokes drag is taken into account.<br />
Numerical results show that heavy particles may accumulate in a fixed point in<br />
the helical coordinate frame, near an elliptic region in the stream function. In physical<br />
space this corresponds to a helically shaped equilibrium trajectory. The relation<br />
between the topology of the carrier flow field and the phenomenon of particle accumulation<br />
becomes clear from figure 1, where the positions of a group of initially uniformly<br />
distributed particles are plotted after a long time. The particle accumulation is proven<br />
analytically by a linear stability analysis. In addition, a full classification of possible<br />
flow field topologies is presented.<br />
These results can be useful in the design of gas-liquid separators, in which a helical<br />
vortex filament is used for the separation of small dust particles and/or very small<br />
droplets 3 .<br />
η<br />
∗ University of Twente, P/O box 217, 7500 AE, Enschede, Netherlands.<br />
1 Hardin, Phys. of Fluids 25(11), (1982).<br />
2 Alekseenko et al., J. Fluid Mech. 382, 185 (1999).<br />
3 Hagmeijer et al., Phys. Fluids 17, 056101 (2005).<br />
2<br />
1.5<br />
1<br />
0.5<br />
0<br />
-0.5<br />
-1<br />
-1.5<br />
-2<br />
-2 -1.5 -1 -0.5 0<br />
ξ<br />
0.5 1 1.5 2<br />
a) ɛ/a = 10 −3 , l/a = 0.5.<br />
η<br />
2<br />
1.5<br />
1<br />
0.5<br />
0<br />
-0.5<br />
-1<br />
-1.5<br />
-2<br />
-2 -1.5 -1 -0.5 0<br />
ξ<br />
0.5 1 1.5 2<br />
b) ɛ/a = 10 −2 , l/a = 1.<br />
η<br />
2<br />
1.5<br />
1<br />
0.5<br />
0<br />
-0.5<br />
-1<br />
-1.5<br />
-2<br />
-2 -1.5 -1 -0.5 0<br />
ξ<br />
0.5 1 1.5 2<br />
c) ɛ/a = 10 −2 , l/a =5.<br />
Figure 1: Positions of heavy particles (blue dots) after time tΓ/a 2 = 100, for three<br />
different values of the thickness of the vortex filament ɛ/a and the helix pitch 2πl/a;<br />
the tube radius R/a = 2 and the particle relaxation time τpΓ/a 2 = 0.4, whereΓ<br />
is the vortex strength and a the radius of the helix geometry. For comparison, the<br />
streamlines of passive tracers are plotted in the background (red lines).