No Arabic abstract
The orientation dynamics of small anisotropic tracer particles in turbulent flows is studied using direct numerical simulation (DNS) and results are compared with Lagrangian stochastic models. Generalizing earlier analysis for axisymmetric ellipsoidal particles (Parsa et al. 2012), we measure the orientation statistics and rotation rates of general, triaxial ellipsoidal tracer particles using Lagrangian tracking in DNS of isotropic turbulence. Triaxial ellipsoids that are very long in one direction, very thin in another, and of intermediate size in the third direction exhibit reduced rotation rates that are similar to those of rods in the ellipsoids longest direction, while exhibiting increased rotation rates that are similar to those of axisymmetric discs in the thinnest direction. DNS results differ significantly from the case when the particle orientations are assumed to be statistically independent from the velocity gradient tensor. They are also different from predictions of a Gaussian process for the velocity gradient tensor, which does not provide realistic preferred vorticity-strain-rate tensor alignments. DNS results are also compared with a stochastic model for the velocity gradient tensor based on the recent fluid deformation approximation (RFDA). Unlike the Gaussian model, the stochastic model accurately predicts the reduction in rotation rate in the longest direction of triaxial ellipsoids since this direction aligns with the flows vorticity, with its rotation perpendicular to the vorticity being reduced. For disc-like particles, or in directions perpendicular to the longest direction in triaxial particles, the model predicts {noticeably} smaller rotation rates than those observed in DNS, a behavior that can be understood based on the probability of vorticity orientation with the most contracting strain-rate eigen-direction in the model.
We investigate the dynamics of cohesive particles in homogeneous isotropic turbulence, based on one-way coupled simulations that include Stokes drag, lubrication, cohesive and direct contact forces. We observe a transient flocculation phase characterized by a growing average floc size, followed by a statistically steady equilibrium phase. We analyze the temporal evolution of floc size and shape due to aggregation, breakage, and deformation. Larger turbulent shear and weaker cohesive forces yield elongated flocs that are smaller in size. Flocculation proceeds most rapidly when the fluid and particle time scales are balanced and a suitably defined Stokes number is textit{O}(1). During the transient stage, cohesive forces of intermediate strength produce flocs of the largest size, as they are strong enough to cause aggregation, but not so strong as to pull the floc into a compact shape. Small Stokes numbers and weak turbulence delay the onset of the equilibrium stage. During equilibrium, stronger cohesive forces yield flocs of larger size. The equilibrium floc size distribution exhibits a preferred size that depends on the cohesive number. We observe that flocs are generally elongated by turbulent stresses before breakage. Flocs of size close to the Kolmogorov length scale preferentially align themselves with the intermediate strain direction and the vorticity vector. Flocs of smaller size tend to align themselves with the extensional strain direction. More generally, flocs are aligned with the strongest Lagrangian stretching direction. The Kolmogorov scale is seen to limit floc growth. We propose a new flocculation model with a variable fractal dimension that predicts the temporal evolution of the floc size and shape.
Small scale characteristics of turbulence such as velocity gradients and vorticity fluctuate rapidly in magnitude and oscillate in sign. Much work exists on the characterization of magnitude variations, but far less on sign oscillations. While averages performed on large scales tend to zero because of the oscillatory character, those performed on increasingly smaller scales will vary with the averaging scale in some characteristic way. This characteristic variation at high Reynolds numbers is captured by the so-called cancellation exponent, which measures how local averages tend to cancel out as the averaging scale increases, in space or time. Past experimental work suggests that the exponents in turbulence depend on whether one considers quantities in full three-dimensional space or uses their one- or two-dimensional cuts. We compute cancellation exponents of vorticity and longitudinal as well as transverse velocity gradients in isotropic turbulence at Taylor-scale Reynolds number up to 1300 on $8192^3$ grids. The 2D cuts yield the same exponents as those for full 3D, while the 1D cuts yield smaller numbers, suggesting that the results in higher dimensions are more reliable. We make the case that the presence of vortical filaments in isotropic turbulence leads to this conclusion. This effect is particularly conspicuous in magnetohydrodynamic turbulence, where an increased degree of spatial coherence develops along the imposed magnetic field.
We present an experimental study on the settling velocity of dense sub-Kolmogorov particles in active-grid-generated turbulence in a wind tunnel. Using phase Doppler interferometry, we observe that the modifications of the settling velocity of inertial particles, under homogeneous isotropic turbulence and dilute conditions $phi_vleq O(10)^{-5}$, is controlled by the Taylor-based Reynolds number $Re_lambda$ of the carrier flow. On the contrary, we did not find a strong influence of the ratio between the fluid and gravity accelerations (i.e., $gammasim(eta/tau_eta^2)/g$) on the particle settling behavior. Remarkably, our results suggest that the hindering of the settling velocity (i.e. the measured particle settling velocity is smaller than its respective one in still fluid conditions) experienced by the particles increases with the value of $Re_lambda$, reversing settling enhancement found under intermediate $Re_lambda$ conditions. This observation applies to all particle sizes investigated, and it is consistent with previous experimental data in the literature. At the highest $Re_lambda$ studied, $Re_lambda>600$, the particle enhancement regime ceases to exist. Our data also show that for moderate Rouse numbers, the difference between the measured particle settling velocity and its velocity in still fluid conditions scales linearly with Rouse, when this difference is normalized by the carrier phase rms fluctuations, i.e., $(V_p-V_T)/usim -Ro$.
We analyze the vector nulls of velocity, Lagrangian acceleration, and vorticity, coming from direct numerical simulations of forced homogeneous isotropic turbulence at $Re_lambda in [40-610]$. We show that the clustering of velocity nulls is much stronger than those of acceleration and vorticity nulls. These acceleration and vorticity nulls, however, are denser than the velocity nulls. We study the scaling of clusters of these null points with $Re_lambda$ and with characteristic turbulence lengthscales. We also analyze datasets of point inertial particles with Stokes numbers $St = 0.5$, 3, and 6, at $Re_lambda = 240$. Inertial particles display preferential concentration with a degree of clustering that resembles some properties of the clustering of the Lagrangian acceleration nulls, in agreement with the proposed sweep-stick mechanism of clustering formation.
Ice crystals settling through a turbulent cloud are rotated by turbulent velocity gradients. In the same way, turbulence affects the orientation of aggregates of organic matter settling in the ocean. In fact most solid particles encountered in Nature are not spherical, and their orientation affects their settling speed, as well as collision rates between particles. Therefore it is important to understand the distribution of orientations of non-spherical particles settling in turbulence. Here we study the angular dynamics of small prolate spheroids settling in homogeneous isotropic turbulence. We consider a limit of the problem where the fluid torque due to convective inertia dominates, so that rods settle essentially horizontally. Turbulence causes the orientation of the settling particles to fluctuate, and we calculate their orientation distribution for prolate spheroids with arbitrary aspect ratios for large settling number Sv (a dimensionless measure of the settling speed), assuming small Stokes number St (a dimensionless measure of particle inertia). This overdamped theory predicts that the orientation distribution is very narrow at large Sv, with a variance proportional to ${rm Sv}^{-4}$. By considering the role of particle inertia, we analyse the limitations of the overdamped theory, and determine its range of applicability. Our predictions are in excellent agreement with numerical simulations of simplified models of turbulent flows. Finally we contrast our results with those of an alternative theory predicting that the orientation variance scales as ${rm Sv}^{-2}$ at large Sv.