• Volume/Page
  • Keyword
  • DOI
  • Citation
  • Advanced
   
 
 
 

Flickr Twitter iResearch App Facebook

Your access to this publication is provided through the subscription of Staats U Unibibl Bremen. Library Welcome Message Help

Library Welcome Message Help

close

Attention, institutional online journal administrators: if the institution name printed in this screen message is incorrect, you can request a change by sending an email to librarynames@aip.org.

Please include your institutional account number, the account administrator's email address and the preferred online display name.

FULL-TEXT OPTIONS:

Phys. Fluids 17, 073301 (2005); http://dx.doi.org/10.1063/1.1940367 (11 pages)

Clustering and collisions of heavy particles in random smooth flows

J. Bec1, A. Celani2, M. Cencini3, and S. Musacchio4

1Dipartimento di Fisica, Università di Roma “La Sapienza,” Piazzale Aldo Moro, 2 I-00185 Roma, Italy and Département Cassiopée, Observatoire de la Côte d’Azur, Boîte Postale 4229, 06304 Nice Cedex 4, France
2Département Cassiopée, Observatoire de la Côte d’Azur, Boîte Postale 4229, 06304 Nice Cedex 4, France and CNRS, INLN, 1361 Route des Lucioles, 06560 Valbonne, France
3Istituto dei Sistemi Complessi ISC-CNR, Via dei Taurini, 19 I-00185 Roma, Italy and SMC-INFM Dipartimento di Fisica, Universit à di Roma “La Sapienza”, Piazzale Aldo Moro, 2 I-00185 Roma, Italy
4Dipartimento di Fisica, Universit à di Roma “La Sapienza,” Piazzale Aldo Moro, 2 I-00185 Roma, Italy and SMC-INFM Dipartimento di Fisica, Universit à di Roma “La Sapienza” Piazzale Aldo Moro, 2 I-00185 Roma, Italy

View MapView Map

(Received 8 July 2004; accepted 4 May 2005; published online 24 June 2005)

Finite-size impurities suspended in incompressible flows distribute inhomogeneously, leading to a drastic enhancement of collisions. A description of the dynamics in the full position-velocity phase space is essential to understand the underlying mechanisms, especially for polydisperse suspensions. These issues are studied here for particles much heavier than the fluid by means of a Lagrangian approach. It is shown that inertia enhances collision rates through two effects: correlation among particle positions induced by the carrier flow and uncorrelation between velocities due to their finite size. A phenomenological model yields an estimate of collision rates for particle pairs with different sizes. This approach is supported by numerical simulations in random flows.

© 2005 American Institute of Physics

Article Outline

  1. INTRODUCTION
  2. DYNAMICS AND STATISTICS OF DILUTE SUSPENSIONS
    1. Lagrangian statistics
  3. MODEL FLOWS
  4. LOCAL DYNAMICS OF MONODISPERSE SUSPENSIONS
  5. AN EXTENSION TO POLYDISPERSE SUSPENSIONS
  6. PHENOMENOLOGICAL MODEL FOR THE COLLISION KERNEL
  7. CONCLUDING REMARKS

KEYWORDS and PACS

PACS

ARTICLE DATA

PUBLICATION DATA

ISSN

1070-6631 (print)  
1089-7666 (online)

    References

  1. J. K. Eaton and J. R. Fessler, "Preferential concentrations of particles by turbulence," Int. J. Multiphase Flow 20, 169 (1994).
  2. M. B. Pinsky and A. P. Khain, "Turbulence effects on droplet growth and size distribution in clouds—a review," J. Aerosol Sci. 28, 1177 (1997). [ISI]
  3. G. Falkovich, A. Fouxon, and M. G. Stepanov, "Acceleration of rain initiation by cloud turbulence," Nature (London) 419, 151 (2002). [ISI] [MEDLINE]
  4. R. A. Shaw, "Particle-turbulence interactions in atmospheric clouds," Annu. Rev. Fluid Mech. 35, 183 (2003). [ISI]
  5. B. J. Rothschild and T. R. Osborn, "Small-scale turbulence and plankton contact rates," J. Plankton Res. 10, 465 (1988). [ISI]
  6. S. Sundby and P. Fossum, "Feeding conditions of arcto-nowegian cod larvae compared with the Rothschild-Osborn theory on small-scale turbulence and plankton contact rates," J. Plankton Res. 12, 1153 (2000). [ISI]
  7. J. Mann, S. Ott, H. L. Pécseli, and J. Trulsen, "Predator-prey enconunters in turbulent waters," Phys. Rev. E 65, 026304 (2002).
  8. A. E. Motter, Y.-C. Lai, and C. Grebogi, "Reactive dynamics of inertial particles in nonhyperbolic chaotic flows," Phys. Rev. E 68, 056307 (2003). [ISI]
  9. T. Nishikawa, Z. Toroczkai, and C. Grebogi, "Advective coalescence in chaotic flows," Phys. Rev. Lett. 87, 038301 (2001). [ISI] [MEDLINE]
  10. P. G. Saffman and J. S. Turner, "On the collision of drops in turbulent clouds," J. Fluid Mech. 1, 16 (1956). [ISI]
  11. J. Abrahamson, "Collision rates of small particles in a vigorously turbulent fluid," Chem. Eng. Sci. 30, 1371 (1975).
  12. J. J. E. Williams and R. I. Crane, "Particle collision rate in turbulent flow," Int. J. Multiphase Flow 9 421, (1983).
  13. F. E. Kruis and K. A. Kusters, "The collision rate of particles in turbulent media," J. Aerosol Sci. 27, S263 (1996).
  14. L. I. Zaichik, O. Simonin, and V. Alipchenkov, "Two statistical models for predicting collision rates of inertial particles in homogeneous isotropic turbulence," Phys. Fluids 15, 2995 (2003)PHFLE6000015000010002995000001.
  15. S. Sundaram and L. R. Collins, "Collision statistics in an isotropic, particle-laden turbulent suspensions," J. Fluid Mech. 335, 75 (1997).
  16. W. C. Reade and L. R. Collins, "Effect of preferential concentration on turbulent collision rates," Phys. Fluids 12, 2530 (2000)PHFLE6000012000010002530000001.
  17. Y. Zhou, A. S. Wexler, and L.-P. Wang, "Modelling turbulent collision of bidisperse inertial particles," J. Fluid Mech. 433, 77 (2001). [ISI]
  18. H. Sigurgeirsson and A. M. Stuart, "A model for preferential concentration," Phys. Fluids 14, 4352 (2002)PHFLE6000014000012004352000001. [ISI]
  19. J.-P. Eckmann and D. Ruelle, "Ergodic theory of chaos and strange attractors," Rev. Mod. Phys. 57, 617 (1985).
  20. L.-P. Wang, A. S. Wexler, and Y. Zhou, "On the collision rate of small particles in isotropic turbulence. Part I. Zero-inertia case," Phys. Fluids 10, 266 (1998)PHFLE6000010000001000266000001. [ISI]
  21. Y. Zhou, L.-P. Wang, and A. S. Wexler, "On the collision rate of small particles in isotropic turbulence. Part II. Finite-inertia case," Phys. Fluids 10, 1206 (1998)PHFLE6000010000005001206000001.
  22. M. R. Maxey and J. Riley, "Equation of motion of a small rigid sphere in a nonuniform flow," Phys. Fluids 26, 883 (1983)PFLDAS000026000004000883000001.
  23. M. R. Maxey, "The gravitational settling of aerosol particles in homogeneous turbulence and random flow fields," J. Fluid Mech. 174, 441 (1987).
  24. P. Grassberger, "Generalized dimensions of strange attractors," Phys. Lett. 97, 227 (1983).
  25. H. G. E. Hentschel and I. Procaccia, "The infinite number of generalized dimensions of fractals and strange attractors," Physica D 8, 435 (1983).
  26. K. Gawedzki and M. Vergassola, "Phase transition in the passive scalar advection," Physica D 138, 63 (2000). [Inspec] [ISI]
  27. B. Mehlig and M. Wilkinson, "Coagulation by random velocity fields as a Kramers problem," Phys. Rev. Lett. 92, 250602 (2004). [ISI] [MEDLINE]
  28. G. Wurm, J. Blum, and J. E. Colwell, "Aerodynamical sticking of dust aggregates," Phys. Rev. E 64, 046301 (2001). [ISI] [MEDLINE]
  29. W. C. Reade and L. R. Collins, "A numerical study of the particle size distribution of an aerosol undergoing turbulent coagulation," J. Fluid Mech. 415, 45 (2000).
  30. This requires to assume ergodicity for the particle dynamics, which is usually the case.
  31. H. Sigurgeirsson and A. M. Stuart, "Inertial particles in a random field," Stochastics Dyn. 2, 295 (2002).
  32. E. Balkovsky, G. Falkovich, and A. Fouxon, "Intermittent distribution of inertial particles in turbulent flows," Phys. Rev. Lett. 86, 2790 (2001). [ISI] [MEDLINE]
  33. L. I. Zaichik and V. Alipchenkov, "Pair dispersion and preferential concentration of particles in isotropic turbulence," Phys. Fluids 15, 1776 (2003)PHFLE6000015000006001776000001. [ISI]
  34. J. Bec, "Fractal clustering of inertial particles in random flows," Phys. Fluids 15, L81 (2003)PHFLE6000015000011000L81000001.
  35. G. Falkovich and A. Pumir, "Intermittent distribution of heavy particles in a turbulent flow," Phys. Fluids 16, L47 (2004)PHFLE6000016000007000L47000001. [ISI]
  36. L.-P. Wang, A. S. Wexler, and Y. Zhou, "Statistical mechanics description and modeling of turbulent collision of inertial particles," J. Fluid Mech. 415, 117 (2000).
  37. S. E. Elghobashi and T. W. Abou-Arab, "A two-equation turbulence model for two-phase flows," Phys. Fluids 26, 931 (1983)PFLDAS000026000004000931000001.
  38. C.T. Crowe, M. Sommerfeld, Y. Tsuji, and C. Crowe, Multiphase Flows with Droplets and Particles (CRC, Boca Raton, FL, 1997).
  39. G. Boffetta, F. De Lillo, and A. Gamba, "Large scale inhomogeneity of inertial particles in turbulent flows," Phys. Fluids 16, L20 (2004)PHFLE6000016000004000L20000001. [ISI]
  40. B. I. Shraiman and E. Siggia, "Scalar turbulence," Nature (London) 405, 639 (2000). [MEDLINE]
  41. G. Falkovich, K. Gawedzki, and M. Vergassola, "Particles and fields in fluid turbulence," Rev. Mod. Phys. 73, 913 (2001).


Figures (click on thumbnails to view enlargements)

FIG.1
(Color online). (a) Scaling exponents μ of the density correlations and γ of the rates of approach as a function of the Stokes number for the time-correlated Gaussian random flow. Time averages are performed over (5×106)τf. (b) Same as in (a) for the shear flow with T = L = 1, U = math. Time averages are over (2×109)T. In both cases the exponents are measured as the mean logarithmic derivative of K over two to three decades in r.

Export Figure to PowerPoint

FIG.2
Modulus of the particle velocities as a function of their positions for two different values of the Stokes number (a) S = 10−3 and (b) S = 1. At small S the surface identified by the particles velocities is very close to that of the modulus of the Eulerian velocity field, meaning that at fixed spatial position the distribution of the particle velocity is sharply peaked on the fluid one. For larger Stokes numbers, the attractor folds in the velocity direction allowing particles to be very close with very different velocities.

Export Figure to PowerPoint

FIG.3
(Color online) Snapshot of the positions of N = 4×105 particles associated to two different values of the Stokes number, S = 0.4 (black) and S = 0.5 (gray, red online) for the random flow ( 9 ). The upper-left inset shows a zoom illustrating the effects induced by the difference in Stokes number (see text for details).

Export Figure to PowerPoint

FIG.4
Rescaled cumulative probability P1,2(r) that two particles are at a distance smaller than r (lower curve) and cumulative approaching rate K1,2(r) (upper curve) as a function of r/r*. Results are shown for S = 0.18 for different values of θ as in the label. Note the collapse of the different curves; the dashed lines indicate the scaling behavior below and above r*. It is worth also noticing that for r<r* both P and K scale as r2 while for r>r* they scale as rμ(S)+1 and rγ(S)+1, respectively. As one can see for such small values of S γμ+1. Data refer to simulations in the random shear flow ( 10 ) with T = 1, L = 1, U = math. The curve for P1,2 has been shifted down for plotting purposes.

Export Figure to PowerPoint

FIG.5
Sketch of the different regions in the (a1,a2) plane corresponding to different contributions of κ1,2 to the effective collision rates.

Export Figure to PowerPoint

FIG.6
(Color online). Typical functional shape (bold line) of the collision kernel, represented here as a function of a2 through a cut with a1 = 0.02 fixed. We chose here γ = 1.8 and the numerical factors and constants to fit the order of magnitude obtained in our numerical experiments. The effective kernel takes one of the different functional forms ( 23 , 25 , 26 ) (represented as dotted lines), depending on whether (a1,a2) is in the region A, B, or C.

Export Figure to PowerPoint

FIG.7
Effective interparticle collision rate Q(a1,a2) obtained numerically by considering ghost collisions in the case of the time-correlated random flow. The Stokes number S and the particle radius a are related by S = 2ρpa2/(9ρfL2) with the choice ρp/ρf = 4.5×103. To obtain these rates, interparticle statistics are computed for 100 different values of the Stokes number.

Export Figure to PowerPoint

FIG.8
(Color online). Collision rate Q(a1,a2) obtained with the same settings as for Fig. 7 and represented here for equal-size particles, i.e., a1 = a2 (solid curve, black online) and as a function of a2 for two fixed values of a1: a1 = 1.58×10−2 (dashed curve, blue online), (dotted curve, red online) a1 = 3.54×10−2.

Export Figure to PowerPoint

FIG.9
(Color online). Solid curve (black online): collision rate Q(a,a) for equal-size particles, normalized by that obtained when neglecting particle inertia, i.e., the Saffman and Turner result equation ( 20 ). Dashed curve (blue online): collision rate obtained when neglecting correlations between the velocity difference and the density (i.e., assuming that the first is just proportional to the separation between the particles), normalized in the same way. These curves were obtained numerically with the same setting as for Fig. 7; the Saffman–Turner kernel was evaluated by following tracer particles advected by the same flow.

Export Figure to PowerPoint



Close
Google Calendar
ADVERTISEMENT

close