scholarly journals Reduced particle settling speed in turbulence

2016 ◽  
Vol 808 ◽  
pp. 153-167 ◽  
Author(s):  
Walter Fornari ◽  
Francesco Picano ◽  
Gaetano Sardina ◽  
Luca Brandt

We study the settling of finite-size rigid spheres in sustained homogeneous isotropic turbulence (HIT) by direct numerical simulations using an immersed boundary method to account for the dispersed solid phase. We study semi-dilute suspensions at different Galileo numbers, $Ga$. The Galileo number is the ratio between buoyancy and viscous forces, and is here varied via the solid-to-fluid density ratio $\unicode[STIX]{x1D70C}_{p}/\unicode[STIX]{x1D70C}_{f}$. The focus is on particles that are slightly heavier than the fluid. We find that in HIT, the mean settling speed is less than that in quiescent fluid; in particular, it reduces by 6 %–60 % with respect to the terminal velocity of an isolated sphere in quiescent fluid as the ratio between the latter and the turbulent velocity fluctuations $u^{\prime }$ is decreased. Analysing the fluid–particle relative motion, we find that the mean settling speed is progressively reduced while reducing $\unicode[STIX]{x1D70C}_{p}/\unicode[STIX]{x1D70C}_{f}$ due to the increase of the vertical drag induced by the particle cross-flow velocity. Unsteady effects contribute to the mean overall drag by about 6 %–10 %. The probability density functions of particle velocities and accelerations reveal that these are closely related to the features of the turbulent flow. The particle mean-square displacement in the settling direction is found to be similar for all $Ga$ if time is scaled by $(2a)/u^{\prime }$ (where $2a$ is the particle diameter and $u^{\prime }$ is the turbulence velocity root mean square).

2016 ◽  
Vol 788 ◽  
pp. 640-669 ◽  
Author(s):  
Walter Fornari ◽  
Francesco Picano ◽  
Luca Brandt

Sedimentation of a dispersed solid phase is widely encountered in applications and environmental flows, yet little is known about the behaviour of finite-size particles in homogeneous isotropic turbulence. To fill this gap, we perform direct numerical simulations of sedimentation in quiescent and turbulent environments using an immersed boundary method to account for the dispersed rigid spherical particles. The solid volume fractions considered are ${\it\phi}=0.5{-}1\,\%$, while the solid to fluid density ratio ${\it\rho}_{p}/{\it\rho}_{f}=1.02$. The particle radius is chosen to be approximately six Kolmogorov length scales. The results show that the mean settling velocity is lower in an already turbulent flow than in a quiescent fluid. The reductions with respect to a single particle in quiescent fluid are approximately 12 % and 14 % for the two volume fractions investigated. The probability density function of the particle velocity is almost Gaussian in a turbulent flow, whereas it displays large positive tails in quiescent fluid. These tails are associated with the intermittent fast sedimentation of particle pairs in drafting–kissing–tumbling motions. The particle lateral dispersion is higher in a turbulent flow, whereas the vertical one is, surprisingly, of comparable magnitude as a consequence of the highly intermittent behaviour observed in the quiescent fluid. Using the concept of mean relative velocity we estimate the mean drag coefficient from empirical formulae and show that non-stationary effects, related to vortex shedding, explain the increased reduction in mean settling velocity in a turbulent environment.


Fluids ◽  
2021 ◽  
Vol 6 (3) ◽  
pp. 111
Author(s):  
Leonid M. Ivanov ◽  
Collins A. Collins ◽  
Tetyana Margolina

Using discrete wavelets, a novel technique is developed to estimate turbulent diffusion coefficients and power exponents from single Lagrangian particle trajectories. The technique differs from the classical approach (Davis (1991)’s technique) because averaging over a statistical ensemble of the mean square displacement (<X2>) is replaced by averaging along a single Lagrangian trajectory X(t) = {X(t), Y(t)}. Metzler et al. (2014) have demonstrated that for an ergodic (for example, normal diffusion) flow, the mean square displacement is <X2> = limT→∞τX2(T,s), where τX2 (T, s) = 1/(T − s) ∫0T−s(X(t+Δt) − X(t))2 dt, T and s are observational and lag times but for weak non-ergodic (such as super-diffusion and sub-diffusion) flows <X2> = limT→∞≪τX2(T,s)≫, where ≪…≫ is some additional averaging. Numerical calculations for surface drifters in the Black Sea and isobaric RAFOS floats deployed at mid depths in the California Current system demonstrated that the reconstructed diffusion coefficients were smaller than those calculated by Davis (1991)’s technique. This difference is caused by the choice of the Lagrangian mean. The technique proposed here is applied to the analysis of Lagrangian motions in the Black Sea (horizontal diffusion coefficients varied from 105 to 106 cm2/s) and for the sub-diffusion of two RAFOS floats in the California Current system where power exponents varied from 0.65 to 0.72. RAFOS float motions were found to be strongly non-ergodic and non-Gaussian.


1991 ◽  
Vol 46 (7) ◽  
pp. 616-620 ◽  
Author(s):  
Junko Habasaki

MD simulation has been performed to learn the microscopic mechanism of diffusion of ions in the Li2SiO3 system. The motion of lithium ions can be explained by the trapping model, where lithium is trapped in the polyhedron and moves with fluctuation of the coordination number. The mean square displacement of lithium was found to correlate well with the net changes in coordination number.


1994 ◽  
Vol 08 (24) ◽  
pp. 3411-3422 ◽  
Author(s):  
W. SCHOMMERS

The effect of premelting is of particular interest in connection with the theory of melting. In this paper, we discuss the structural and dynamical properties of the surfaces of semi-infinite crystals as well as of nano-clusters, which show the effect of premelting. The investigations are based on molecular-dynamics calculations: different models are used for the systematic study of the effect of premelting. In particular, the behaviour of the following functions have been studied: pair correlation function, generalized phonon density of states, and the mean-square displacement as a function of time. The calculations have been done for krypton since for this substance a reliable interaction potential is available.


2018 ◽  
Vol 32 (19) ◽  
pp. 1850210
Author(s):  
Chun-Yang Wang ◽  
Zhao-Peng Sun ◽  
Ming Qin ◽  
Yu-Qing Xu ◽  
Shu-Qin Lv ◽  
...  

We report, in this paper, a recent study on the dynamical mechanism of Brownian particles diffusing in the fractional damping environment, where several important quantities such as the mean square displacement (MSD) and mean square velocity are calculated for dynamical analysis. A particular type of backward motion is found in the diffusion process. The reason of it is analyzed intrinsically by comparing with the diffusion in various dissipative environments. Results show that the diffusion in the fractional damping environment obeys the Langevin dynamics which is quite different form what is expected.


1977 ◽  
Vol 44 (3) ◽  
pp. 487-491 ◽  
Author(s):  
S. F. Masri ◽  
F. Udwadia

The transient mean-square displacement, slope, and relative motion of a viscously damped shear beam subjected to correlated random boundary excitation is presented. The effects of various system parameters including the spectral characteristics of the excitation, the delay time between the beam support motion, and the beam damping have been investigated. Marked amplifications in the mean-square response are shown to occur for certain dimensionless time delays.


2019 ◽  
Vol 2019 ◽  
pp. 1-5
Author(s):  
Long Shi

In this work, a generalization of continuous time random walk is considered, where the waiting times among the subsequent jumps are power-law correlated with kernel function M(t)=tρ(ρ>-1). In a continuum limit, the correlated continuous time random walk converges in distribution a subordinated process. The mean square displacement of the proposed process is computed, which is of the form 〈x2(t)〉∝tH=t1/(1+ρ+1/α). The anomy exponent H varies from α to α/(1+α) when -1<ρ<0 and from α/(1+α) to 0 when ρ>0. The generalized diffusion equation of the process is also derived, which has a unified form for the above two cases.


Author(s):  
Ali Khalili Golmankhaneh ◽  
Saleh Ashrafi ◽  
Dumitru Baleanu ◽  
Arran Fernandez

AbstractIn this paper, we have investigated the Langevin and Brownian equations on fractal time sets using Fα-calculus and shown that the mean square displacement is not varied linearly with time. We have also generalized the classical method of deriving the Fokker–Planck equation in order to obtain the Fokker–Planck equation on fractal time sets.


Sign in / Sign up

Export Citation Format

Share Document