scholarly journals Statistical Model of Turbulent Dispersion Recapitulated

Fluids ◽  
2021 ◽  
Vol 6 (5) ◽  
pp. 190
Author(s):  
J. J. H. Brouwers

A comprehensive summary and update is given of Brouwers’ statistical model that was developed during the previous decade. The presented recapitulated model is valid for general inhomogeneous anisotropic velocity statistics that are typical of turbulence. It succeeds and improves the semiempirical and heuristic models developed during the previous century. The model is based on a Langevin and diffusion equation of which the derivation involves (i) the application of general principles of physics and stochastic theory; (ii) the application of the theory of turbulence at large Reynolds numbers, including the Lagrangian versions of the Kolmogorov limits; and (iii) the systematic expansion in powers of the inverse of the universal Lagrangian Kolmogorov constant C0, C0 about 6. The model is unique in the collected Langevin and diffusion models of physics and chemistry. Presented results include generally applicable expressions for turbulent diffusion coefficients that can be directly implemented in numerical codes of computational fluid mechanics used in environmental and industrial engineering praxis. This facilitates the more accurate and reliable prediction of the distribution of the mean concentration of passive or almost passive admixture such as smoke, aerosols, bacteria, and viruses in turbulent flow, which are all issues of great societal interest.

2010 ◽  
Vol 659 ◽  
pp. 405-419 ◽  
Author(s):  
K. YEO ◽  
B.-G. KIM ◽  
C. LEE

The behaviour of fluid-particle acceleration in near-wall turbulent flows is investigated in numerically simulated turbulent channel flows at low to moderate Reynolds numbers, Reτ = 180~600). The acceleration is decomposed into pressure-gradient (irrotational) and viscous contributions (solenoidal acceleration) and the statistics of each component are analysed. In near-wall turbulent flows, the probability density function of acceleration is strongly dependent on the distance from the wall. Unexpectedly, the intermittency of acceleration is strongest in the viscous sublayer, where the acceleration flatness factor of O(100) is observed. It is shown that the centripetal acceleration around coherent vortical structures is an important source of the acceleration intermittency. We found sheet-like structures of strong solenoidal accelerations near the wall, which are associated with the background shear modified by the interaction between a streamwise vortex and the wall. We found that the acceleration Kolmogorov constant is a linear function of y+ in the log layer. The Reynolds number dependence of the acceleration statistics is investigated.


1978 ◽  
Vol 88 (4) ◽  
pp. 685-709 ◽  
Author(s):  
U. Schumann ◽  
G. S. Patterson

The spectral method of Orszag & Patterson has been extended to calculate the static pressure fluctuations in incompressible homogeneous decaying turbulence at Reynolds numbers Reλ [lsim ] 35. In real space 323 points are treated. Several cases starting from different isotropic initial conditions have been studied. Some departure from isotropy exists owing to the small number of modes at small wavenumbers. Root-mean-square pressure fluctuations, pressure gradients and integral length scales have been evaluated. The results agree rather well with predictions based on velocity statistics and on the assumption of normality. The normality assumption has been tested extensively for the simulated fields and found to be approximately valid as far as fourth-order velocity correlations are concerned. In addition, a model for the dissipation tensor has been proposed. The application of the present method to the study of the return of axisymmetric turbulence to isotropy is described in the companion paper.


Author(s):  
K. M. Akyuzlu ◽  
J. Farkas

An experimental study was conducted to observe/visualize, the formation of circulation patterns inside a square cavity due to the movement of a lid at constant velocity. Lid driven cavity flow is one of the benchmark studies used in the verification/improvement of CFD codes for internal flow applications/predictions. Previous work on this topic is primarily focused on improving the steady state predictions of the CFD codes using different numerical schemes and algorithms. Furthermore, almost all of the studies reported in computational fluid mechanics literature relates to steady state predictions of lid or shear driven flows. Experimental work that is reported in these studies is limited in scope and number. This paper reports on the measurements we made using Particle Image Velocimeter (PIV) technique to determine the flow field as it develops from stagnation to steady state inside a square cavity driven by a lid. For this purpose, we employed a 2-D PIV system, which uses a double-cavity, Nd:Yag laser to illuminate the test cavity. Experiments were conducted using water as the working fluid inside a square cavity that is one inch (25.4 mm) high and one inch wide. The depth of the cavity is five inches (127 mm) to ensure two-dimensional circulations patterns. Hollow glass sphere particles with 10 microns in diameter were used as seeding of the working fluid, water. Experiments were repeated for different lid velocities corresponding to lid Reynolds numbers (laminar to beginning of transition of turbulence.) Velocity fields were captured during the development of the circulations patters each being unique for the time of the measurement and value of the lid velocity. The center of the circulation pattern and its path inside the cavity is constructed from the captured images as steady state is attained. Also, the strength of the circulation (as manifested by the increase in the diameter of the circulation) is determined at different times for different Reynolds numbers.


1998 ◽  
Vol 357 ◽  
pp. 167-198 ◽  
Author(s):  
BURKHARD M. O. HEPPE

The relative velocity of two fluid particles in homogeneous and stationary turbulence is considered. Looking for reduced dynamics of turbulent dispersion, we apply the nonlinear Mori–Zwanzig projector method to the Navier–Stokes equations. The projector method decomposes the Lagrangian acceleration into a conditionally averaged part and a random force. The result is an exact generalized Langevin equation for the Lagrangian velocity differences accounting for the exact equation of the Eulerian probability density. From the generalized Langevin equation, we obtain a stochastic model of relative dispersion by stochastic estimation of conditional averages and by assuming the random force to be Gaussian white noise. This new approach to dispersion modelling generalizes and unifies stochastic models based on the well-mixed condition and the moments approximation. Furthermore, we incorporate viscous effects in a systematic way. At a moderate Reynolds number, the model agrees qualitatively with direct numerical simulations showing highly non-Gaussian separation and velocity statistics for particle pairs initially close together. At very large Reynolds numbers, the mean-square separation obeys a Richardson law with coefficient of the order of 0.1.


Author(s):  
A. Gharakhani ◽  
J. Sitaraman ◽  
M. J. Stock

A grid-free method for the simulation of incompressible flow over complex 3-D objects is presented. The algorithm is based on the Lagrangian Vortex Element and the Vorticity Redistribution methods to account for vorticity dynamics and diffusion, respectively. The wall velocity boundary conditions are imposed by solving a Petrov-Galerkin discretization of the Fredholm boundary integral equation of the second kind for the vector-valued vortex sheets at the wall. Computations are accelerated using a MPI-parallel adaptive Fast Multipole Method (oct-tree code). In this paper, the computational algorithm is described briefly, and the first in a series of benchmark results presented using the example of flow over a sphere at low Reynolds numbers. Preliminary diagnostics of the present simulations show very good agreement with previously published data.


2017 ◽  
Vol 826 ◽  
Author(s):  
Christian E. Willert ◽  
Julio Soria ◽  
Michel Stanislas ◽  
Joachim Klinner ◽  
Omid Amili ◽  
...  

This paper reports on near-wall two-component–two-dimensional (2C–2D) particle image velocimetry (PIV) measurements of a turbulent pipe flow at shear Reynolds numbers up to $Re_{\unicode[STIX]{x1D70F}}=40\,000$ acquired in the CICLoPE facility of the University of Bologna. The 111.5 m long pipe of 900 mm diameter offers a well-established turbulent flow with viscous length scales ranging from $85~\unicode[STIX]{x03BC}\text{m}$ at $Re_{\unicode[STIX]{x1D70F}}=5000$ down to $11~\unicode[STIX]{x03BC}\text{m}$ at $Re_{\unicode[STIX]{x1D70F}}=40\,000$. These length scales can be resolved with a high-speed PIV camera at image magnification near unity. Statistically converged velocity profiles were determined using multiple sequences of up to 70 000 PIV recordings acquired at sampling rates of 100 Hz up to 10 kHz. Analysis of the velocity statistics shows a well-resolved inner peak of the streamwise velocity fluctuations that grows with increasing Reynolds number and an outer peak that develops and moves away from the inner peak with increasing Reynolds number.


2016 ◽  
Vol 796 ◽  
pp. 617-658 ◽  
Author(s):  
Peter J. Ireland ◽  
Andrew D. Bragg ◽  
Lance R. Collins

In this study, we analyse the statistics of both individual inertial particles and inertial particle pairs in direct numerical simulations of homogeneous isotropic turbulence in the absence of gravity. The effect of the Taylor microscale Reynolds number, $R_{{\it\lambda}}$, on the particle statistics is examined over the largest range to date (from $R_{{\it\lambda}}=88$ to 597), at small, intermediate and large Kolmogorov-scale Stokes numbers $St$. We first explore the effect of preferential sampling on the single-particle statistics and find that low-$St$ inertial particles are ejected from both vortex tubes and vortex sheets (the latter becoming increasingly prevalent at higher Reynolds numbers) and preferentially accumulate in regions of irrotational dissipation. We use this understanding of preferential sampling to provide a physical explanation for many of the trends in the particle velocity gradients, kinetic energies and accelerations at low $St$, which are well represented by the model of Chun et al. (J. Fluid Mech., vol. 536, 2005, pp. 219–251). As $St$ increases, inertial filtering effects become more important, causing the particle kinetic energies and accelerations to decrease. The effect of inertial filtering on the particle kinetic energies and accelerations diminishes with increasing Reynolds number and is well captured by the models of Abrahamson (Chem. Engng Sci., vol. 30, 1975, pp. 1371–1379) and Zaichik & Alipchenkov (Intl J. Multiphase Flow, vol. 34 (9), 2008, pp. 865–868), respectively. We then consider particle-pair statistics, and focus our attention on the relative velocities and radial distribution functions (RDFs) of the particles, with the aim of understanding the underlying physical mechanisms contributing to particle collisions. The relative velocity statistics indicate that preferential sampling effects are important for $St\lesssim 0.1$ and that path-history/non-local effects become increasingly important for $St\gtrsim 0.2$. While higher-order relative velocity statistics are influenced by the increased intermittency of the turbulence at high Reynolds numbers, the lower-order relative velocity statistics are only weakly sensitive to changes in Reynolds number at low $St$. The Reynolds-number trends in these quantities at intermediate and large $St$ are explained based on the influence of the available flow scales on the path-history and inertial filtering effects. We find that the RDFs peak near $St$ of order unity, that they exhibit power-law scaling for low and intermediate $St$ and that they are largely independent of Reynolds number for low and intermediate $St$. We use the model of Zaichik & Alipchenkov (New J. Phys., vol. 11, 2009, 103018) to explain the physical mechanisms responsible for these trends, and find that this model is able to capture the quantitative behaviour of the RDFs extremely well when direct numerical simulation data for the structure functions are specified, in agreement with Bragg & Collins (New J. Phys., vol. 16, 2014a, 055013). We also observe that at large $St$, changes in the RDF are related to changes in the scaling exponents of the relative velocity variances. The particle collision kernel closely matches that computed by Rosa et al. (New J. Phys., vol. 15, 2013, 045032) and is found to be largely insensitive to the flow Reynolds number. This suggests that relatively low-Reynolds-number simulations may be able to capture much of the relevant physics of droplet collisions and growth in the adiabatic cores of atmospheric clouds.


2015 ◽  
Vol 774 ◽  
pp. 5-36 ◽  
Author(s):  
Man Mohan Rai

The near and very near wake of a flat plate with a circular trailing edge, with vigorous vortex shedding, is investigated with data from direct numerical simulations (DNS). Computations were performed for four different combinations of the Reynolds numbers based on plate thickness ($D$) and momentum thickness near the trailing edge (${\it\theta}$). Unlike the case of the cylinder, these Reynolds numbers are independent parameters for the flat plate. The objectives of the study are twofold, to investigate the entrainment process when the separating boundary layers are turbulent and to better understand the instability of the detached shear layers (DSLs). A visualization of the entrainment process, the effect of changing the ratio ${\it\theta}/D$ on entrainment and wake-velocity statistics, and a way of understanding entrainment in a phase-averaged sense via distributions of the turbulent transport rate are provided here. The discussion on shear-layer instability focuses on the role of log-layer eddies in the destabilization process, the effect of high-speed streaks in the turbulent boundary layer in the vicinity of the trailing edge on shear-layer vortex generation rates, and a relationship between the prevalence of shear-layer vortex generation and shedding phase that is a result of an interaction between the shedding process and the shear-layer instability mechanism. A power-law relationship between the ratio of shear-layer and shedding frequencies and the Reynolds numbers mentioned above is obtained. A discussion of the relative magnitudes of the exponents is provided. A second power-law relationship between shed-vortex strength and these two Reynolds numbers is also proposed.


2010 ◽  
Vol 645 ◽  
pp. 497-528 ◽  
Author(s):  
J. BEC ◽  
L. BIFERALE ◽  
A. S. LANOTTE ◽  
A. SCAGLIARINI ◽  
F. TOSCHI

The relative dispersion of pairs of inertial point particles in incompressible, homogeneous and isotropic three-dimensional turbulence is studied by means of direct numerical simulations at two values of the Taylor-scale Reynolds number Reλ ~ 200 and Reλ ~ 400, corresponding to resolutions of 5123 and 20483 grid points, respectively. The evolution of both heavy and light particle pairs is analysed by varying the particle Stokes number and the fluid-to-particle density ratio. For particles much heavier than the fluid, the range of available Stokes numbers is St ∈ [0.1 : 70], while for light particles the Stokes numbers span the range St ∈ [0.1 : 3] and the density ratio is varied up to the limit of vanishing particle density. For heavy particles, it is found that turbulent dispersion is schematically governed by two temporal regimes. The first is dominated by the presence, at large Stokes numbers, of small-scale caustics in the particle velocity statistics, and it lasts until heavy particle velocities have relaxed towards the underlying flow velocities. At such large scales, a second regime starts where heavy particles separate as tracers' particles would do. As a consequence, at increasing inertia, a larger transient stage is observed, and the Richardson diffusion of simple tracers is recovered only at large times and large scales. These features also arise from a statistical closure of the equation of motion for heavy particle separation that is proposed and is supported by the numerical results. In the case of light particles with high density ratio, strong small-scale clustering leads to a considerable fraction of pairs that do not separate at all, although the mean separation increases with time. This effect strongly alters the shape of the probability density function of light particle separations.


Sign in / Sign up

Export Citation Format

Share Document