1
|
Lee S, Lee C. Identification of a particle collision as a finite-time blowup in turbulence. Sci Rep 2023; 13:181. [PMID: 36604443 PMCID: PMC9814516 DOI: 10.1038/s41598-022-27305-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2022] [Accepted: 12/29/2022] [Indexed: 01/06/2023] Open
Abstract
We propose an Eulerian approach to investigate the motion of particles in turbulence under the assumption that the motion of particles remains smooth in space and time until a collision between particles occurs. When the first collision happens, particle velocity loses [Formula: see text] continuity, resulting in a finite-time blowup. The corresponding singularities in particle velocity gradient, particle number density, and particle vorticity for various Stokes numbers and gravity factors are numerically investigated for the first time in a simple two-dimensional Taylor-Green vortex flow, two-dimensional decaying turbulence, and three-dimensional isotropic turbulence. In addition to the critical Stokes number above which a collision begins to occur, the flow condition leading to collision is revealed; particles tend to collide in very thin shear layer constructed by two parallel same-signed vortical structures when Stokes number is above the critical one.
Collapse
Affiliation(s)
- Seulgi Lee
- grid.15444.300000 0004 0470 5454Department of Mechanical Engineering, Yonsei University, Seoul, 03722 South Korea
| | - Changhoon Lee
- grid.15444.300000 0004 0470 5454Department of Mechanical Engineering & School of Mathematics and Computing, Yonsei University, Seoul, 03722 South Korea
| |
Collapse
|
2
|
Nicolleau FCGA, Farhan M, Nowakowski AF. Effect of the energy-spectrum law on clustering patterns for inertial particles subjected to gravity in kinematic simulation. Phys Rev E 2016; 94:043109. [PMID: 27841627 DOI: 10.1103/physreve.94.043109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/29/2016] [Indexed: 06/06/2023]
Abstract
We study the clustering of inertial particles using a periodic kinematic simulation. Particle clustering is observed for different pairs of Stokes number and Froude number and different spectral power laws (1.4≤p≤2.1). The main focus is to identify and then quantify the effect of p on the clustering attractor-by attractor we mean the set of points in the physical space where the particles settle when time tends to infinity. It is observed that spectral power laws can have a dramatic effect on the attractor shape. In particular, we observed an attractor type which was not present in previous studies for Kolmogorov spectra (p=5/3).
Collapse
Affiliation(s)
- F C G A Nicolleau
- Sheffield Fluid Mechanics Group, Department of Mechanical Engineering, The University of Sheffield, Sheffield, United Kingdom and Sorbonne Universités, Université Pierre et Marie Curie, Paris 6 and Institut Jean le Rond dAlembert, CNRS UMR 7190, Paris, France
| | - M Farhan
- Department of Mechanical Engineering of the University of Engineering and Technology Lahore, Pakistan and Sheffield Fluid Mechanics Group, Department of Mechanical Engineering, The University of Sheffield, Sheffield, United Kingdom
| | - A F Nowakowski
- Sheffield Fluid Mechanics Group, Department of Mechanical Engineering, The University of Sheffield, Sheffield, United Kingdom
| |
Collapse
|
3
|
Choi JI, Park Y, Kwon O, Lee C. Interparticle collision mechanism in turbulence. Phys Rev E 2016; 93:013112. [PMID: 26871160 DOI: 10.1103/physreve.93.013112] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2015] [Indexed: 11/06/2022]
Abstract
Direct numerical simulations of particle-laden homogeneous isotropic turbulence are performed to investigate interparticle collisions in a wide range of Stokes numbers. Dynamics of the particles are described by Stokes drag including particle-particle interactions via hard-sphere collisions, while fluid turbulence is solved using a pseudospectral method. Particular emphasis is placed on interparticle-collision-based conditional statistics of rotation and dissipation rates of the fluid experienced by heavy particles, which provide essential information on the collision process. We also investigate the collision statistics of collision time interval and angle. Based on a Lamb vortex model for a vortex structure, we claim that collision events occur in the edge region for vortical structures in the intermediate-Stokes-number regime, suggesting that the sling effect enhances collision as well as clustering.
Collapse
Affiliation(s)
- Jung-Il Choi
- Department of Computational Science and Engineering, Yonsei University, Seoul 120-749, South Korea
| | - Yongnam Park
- Department of Mechanical Engineering, Yonsei University, Seoul 120-749, South Korea
| | - Ohjoon Kwon
- Department of Mechanical Engineering, Yonsei University, Seoul 120-749, South Korea
| | - Changhoon Lee
- Department of Computational Science and Engineering, Yonsei University, Seoul 120-749, South Korea.,Department of Mechanical Engineering, Yonsei University, Seoul 120-749, South Korea
| |
Collapse
|
4
|
Perrin VE, Jonker HJJ. Relative velocity distribution of inertial particles in turbulence: A numerical study. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:043022. [PMID: 26565347 DOI: 10.1103/physreve.92.043022] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/12/2014] [Indexed: 06/05/2023]
Abstract
The distribution of relative velocities between particles provides invaluable information on the rates and characteristics of particle collisions. We show that the theoretical model of Gustavsson and Mehlig [K. Gustavsson and B. Mehlig, J. Turbul. 15, 34 (2014)], within its anticipated limits of validity, can predict the joint probability density function of relative velocities and separations of identical inertial particles in isotropic turbulent flows with remarkable accuracy. We also quantify the validity range of the model. The model matches two limits (or two types) of relative motion between particles: one where pair diffusion dominates (i.e., large coherence between particle motion) and one where caustics dominate (i.e., large velocity differences between particles at small separations). By using direct numerical simulation combined with Lagrangian particle tracking, we assess the model prediction in homogeneous and isotropic turbulence. We demonstrate that, when sufficient caustics are present at a given separation and the particle response time is significantly smaller than the integral time scales of the flow, the distribution exhibits the same universal power-law form dictated by the correlation dimension as predicted by the model of Gustavsson and Mehlig. In agreement with the model, no strong dependency on the Taylor-based Reynolds number is observed.
Collapse
|
5
|
Fouxon I, Park Y, Harduf R, Lee C. Inhomogeneous distribution of water droplets in cloud turbulence. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:033001. [PMID: 26465550 DOI: 10.1103/physreve.92.033001] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/23/2014] [Indexed: 06/05/2023]
Abstract
We consider sedimentation of small particles in the turbulent flow where fluid accelerations are much smaller than acceleration of gravity g. The particles are dragged by the flow by linear friction force. We demonstrate that the pair-correlation function of particles' concentration diverges with decreasing separation as a power law with negative exponent. This manifests fractal distribution of particles in space. We find that the exponent is proportional to ratio of integral of energy spectrum of turbulence times the wave number over g. The proportionality coefficient is a universal number independent of particle size. We derive the spectrum of Lyapunov exponents that describes the evolution of small patches of particles. It is demonstrated that particles separate dominantly in the horizontal plane. This provides a theory for the recently observed vertical columns formed by the particles. We confirm the predictions by direct numerical simulations of Navier-Stokes turbulence. The predictions include conditions that hold for water droplets in warm clouds thus providing a tool for the prediction of rain formation.
Collapse
Affiliation(s)
- Itzhak Fouxon
- Department of Computational Science and Engineering, Yonsei University, Seoul 120-749, South Korea
| | - Yongnam Park
- Department of Mechanical Engineering, Yonsei University, Seoul 120-749, South Korea
| | - Roei Harduf
- Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, Ramat Aviv, Tel Aviv 69978, Israel
| | - Changhoon Lee
- Department of Computational Science and Engineering, Yonsei University, Seoul 120-749, South Korea
- Department of Mechanical Engineering, Yonsei University, Seoul 120-749, South Korea
| |
Collapse
|
6
|
Fouxon I, Leshansky A. Phytoplankton's motion in turbulent ocean. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:013017. [PMID: 26274279 DOI: 10.1103/physreve.92.013017] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/10/2015] [Indexed: 06/04/2023]
Abstract
We study the influence of turbulence on upward motion of phytoplankton. Interaction with the flow is described by the Pedley-Kessler model considering spherical microorganisms. We find a range of parameters when the upward drift is only weakly perturbed or when turbulence completely randomizes the drift direction. When the perturbation is small, the drift is either determined by the local vorticity or is Gaussian. We find a range of parameters where the phytoplankton interaction with the flow can be described consistently as diffusion of orientation in effective potential. By solving the corresponding Fokker-Planck equation we find exponential steady-state distribution of phytoplankton's propulsion orientation. We further identify the range of parameters where phytoplankton's drift velocity with respect to the flow is determined uniquely by its position. In this case, one can describe phytoplankton's motion by a smooth flow and phytoplankton concentrates on fractal. We find fractal dimensions and demonstrate that phytoplankton forms vertical stripes in space with a nonisotropic pair-correlation function of concentration increased in the vertical direction. The probability density function of the distance between two particles obeys power law with the negative exponent given by the ratio of integrals of the turbulent energy spectrum. We find the regime of strong clustering where the exponent is of order one so that turbulence increases the rate of collisions by a large factor. The predictions hold for Navier-Stokes turbulence and stand for testing.
Collapse
Affiliation(s)
- Itzhak Fouxon
- Department of Chemical Engineering, Technion, Haifa 32000, Israel
- Department of Computational Science and Engineering, Yonsei University, Seoul 120-749, South Korea
| | | |
Collapse
|
7
|
Farhan M, Nicolleau FCGA, Nowakowski AF. Effect of gravity on clustering patterns and inertial particle attractors in kinematic simulations. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:043021. [PMID: 25974594 DOI: 10.1103/physreve.91.043021] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/28/2014] [Indexed: 06/04/2023]
Abstract
In this paper, we study the clustering of inertial particles using a periodic kinematic simulation. The systematic Lagrangian tracking of particles makes it possible to identify the particles' clustering patterns for different values of particle inertia and drift velocity. The different cases are characterized by different pairs of Stokes number (St) and Froude number (Fr). For the present study, 0≤St≤1 and 0.4≤Fr≤1.4. The main focus is to identify and then quantify the clustering attractor-when it exists-that is the set of points in the physical space where the particles settle when time goes to infinity. Depending on the gravity effect and inertia values, the Lagrangian attractor can have different dimensions, varying from the initial three-dimensional space to two-dimensional layers and one-dimensional attractors that can be shifted from a horizontal to a vertical position.
Collapse
Affiliation(s)
- M Farhan
- Sheffield Fluid Mechanics Group, Department of Mechanical Engineering, The University of Sheffield, Sheffield, United Kingdom
| | - F C G A Nicolleau
- Sheffield Fluid Mechanics Group, Department of Mechanical Engineering, The University of Sheffield, Sheffield, United Kingdom
| | - A F Nowakowski
- Sheffield Fluid Mechanics Group, Department of Mechanical Engineering, The University of Sheffield, Sheffield, United Kingdom
| |
Collapse
|