1 / 18

Measuring segregation of inertial particles in turbulent flows by a Full Lagrangian approach

Measuring segregation of inertial particles in turbulent flows by a Full Lagrangian approach. 4 th IMS Workshop on Clouds and Turbulence Institute for Mathematical Sciences Imperial College London 23-25 March 2009. E. Meneguz Ph.D. project: Rain in a box of turbulence

jeremy-caldwell
Download Presentation

Measuring segregation of inertial particles in turbulent flows by a Full Lagrangian approach

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Measuring segregation of inertial particles in turbulent flows by a Full Lagrangian approach 4th IMS Workshop on Clouds and Turbulence Institute for Mathematical Sciences Imperial College London 23-25 March 2009 E. Meneguz Ph.D. project: Rain in a box of turbulence Supervisor: M. W. Reeks PostDoc: R. H. A. IJzermans School of Mechanical and Systems Engineering UNIVERSITY OF NEWCASTLE

  2. Outline Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 2 • Introduction; • Equations of motion of inertial particles and other equations; • Compressibility of the Particle Velocity Field: • - MEPVF (Eulerian) - FLA (Lagrangian) • Method • flow field description - numerical simulations • Results: • Compressibility of the particle vel. field in both approaches - preliminary results in DNS of HIT - moments of the particle number density (comparison with theor. predictions) • Conclusions and future developments

  3. State of the art Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 3 • Importance of de-mixing of particles in turbulent flows: many environmental-industrial and statistical interest • Preferential concentration (Crowe et al, 1993 and Maxey (1987); Sundaram and Collins (1997) and Wang et al. (1998); • Recently studied from different viewpoints (Chen et al, 2006, Balkowsky et al, 2001, Sommerer & Ott (1993) and Wilkinson et al, 2005;2007) • Random Uncorrelated Motion - “sling effect” (Falkovich et al. 2002) or “crossing trajectories effect” (Wilkinson et al. 2005) • Fevrier et al. MEF as sum of two contributions: MEPVF and RUM; based on box-counting (EULERIAN); • Osiptsov’s method (Reeks 2004, Healy and Young 2005) LAGRANGIAN

  4. Equations of motion (Lagrangian frame) Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 4 • Particles/droplets: spherical, rigid, identical and heavy • Dilute system • Effect of gravity and Brownian motion not included particle position particle velocity fluid velocity at the position of the particle normalized particle relaxation number

  5. Compressibility of the PVF Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 5 • Flow field incompressible: • Preferential concentration: • PVF compressible • non-zero gradients in the particle number density • Continuity equation: • For sufficiently small Stokes number:

  6. MEPVF and compressibility of PVF Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 6 i-th cell • According to Fevrier et al. 2005: MEPVF = PVF + RUM • PVF = To be obtained from finite differences

  7. Lagrangian quantification of the compressibility Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 7 Evol. eqts. 2D EXAMPLE: Continuity equation and averaging over all particle trajectories:

  8. Model of synthetic turbulent flow Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 8 To study effect of RUM: Threshold value: St=0.25 2D carrier flow field (Babiano et al. 2000):

  9. Numerical methods (Lagrangian vs Eulerian) Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 9 • 10,000 inertial particles, uniformly distributed at t=0; • Periodic boundary conditions for particles; • Trajectories and equations for calculated by RK4 method; • Initial conditions: (Volume is initially a cube)‏. • 1,000,000 inertial particles, uniformly distributed at t=0; • Trajectories calculated by RK4 method; • Periodic boundary conditions for particles; • Divergence of PVF calculated using 2nd order finite difference scheme; • Numerical resolution varied between 102 and 602 cells.

  10. Compressibility in time Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 10 St=0.05 St=0.2 St=2 St=0.5

  11. Influence of numerical resolution Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 11 Resolution:102 cells ---- 202 cells ---- 602 cells ---- d<ln|J|>/dt ---- St=0.2 Lagrangian method corresponds to limiting case for infinitely fine grid in Eulerian box-counting method

  12. Compressibility in DNS of HIT Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 12 St=1 Picciotto et al. 2005 • Good agreement between Lagrangian and Eulerian method • Singularities seem to be detected better by Lagrangian method

  13. Moments of particle number density Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 13 • Along particle trajectory: particle number density n related to J by: • Particle averaged value of is related to spatially averaged value: • Any space-averaged momentis readily determined, if J is known for all particles in the sub-domain (equivalent to counting particles) Trivial limits:

  14. Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 14 Moments of particle number density St=0.2 St=0.5 • Particle number density is spatially strongly intermittent • Sudden peaks indicate singularities in particle velocity field

  15. Comparison with analytical estimate Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 15 If St is sufficiently small: Trivial limits: In agreement with Balkovsky et al (2001, PRL)

  16. Conclusions Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 16 • The Lagrangian and Eulerian method show good agreement in the compressibility of the PVF for a wide range of St; • Singularities in the PVF can be detected by Lagrangian method but not by Eulerian method, due to finite grid size in the latter; • Lagrangian method allows for determination of any space-averaged moment of the particle number density, in contrast with Eulerian which would have too limited spatial resolutions; • The determination of moments of the particle number density have shown very high spatial intermittency due to RUM. • For first time, numerical support for theory of • Balkovsky et al (2001, PRL): “is convex function of”.

  17. Further developments Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 17 • Open question: high intermittency in particle number density in DNS? • 3D DNS of stationary HIT for different St numbers • pdf methods for two particle dispersion at higher Re numbers

  18. Elena Meneguz ● 4th IMS Workshop on Clouds and Turbulence ● Imperial College London ● 24.03.09 slide 18 Thank you for your attention

More Related