Three-dimensional structure of the Lagrangian acceleration in turbulent flows

We report experimental results on the three-dimensional Lagrangian acceleration in highly turbulent flows. Tracer particles are tracked optically using four silicon strip detectors from high energy physics that provide high temporal and spatial resolution. The components of the acceleration are shown to be statistically dependent. The probability density function of the acceleration magnitude is comparable to a log-normal distribution. Assuming isotropy, a log-normal distribution of the magnitude can account for the observed dependency of the components. The time dynamics of the acceleration components is found to be typical of the dissipation scales, whereas the magnitude evolves over longer times, possibly close to the integral time scale.     

Finite-size particles in turbulent channel flow: quadrant analysis and acceleration statistics

The interaction between finite-size particles and turbulent channel flow in the absence of gravity is studied by direct numerical simulations (DNS). The study is motivated by DNS observations of a turbulent channel flow with high-density, pointwise particles, that cluster in regions of high streamwise root-mean-square (RMS) acceleration close to the wall, contrary to what is observed in homogeneous isotropic turbulence. The aim of the present study is to explore if this is still the case when size effects are taken into account in the DNS. Based on the analysis of the velocity and acceleration statistics, the present DNS shows that, close to the wall, particles with ρ_p/ρ_f ranging from 2 to 4 are surrounded by regions with low streamwise RMS velocity but high streamwise RMS acceleration. According to the normalised particle acceleration probability density functions (PDFs), size effects become important in the near-wall region. As particle inertia increases, the normalised PDFs of particle acceleration tend to a Gaussian distribution. The tails of the normalised PDFs of the fluid conditioned by the presence of particles are higher than that of the unconditioned fluid close to the wall.     

Statistics of Lagrangian velocities in turbulent flows

We present a generalized Fokker-Planck equation for the joint position-velocity probability distribution of a single fluid particle in a turbulent flow. Based on a simple estimate, the diffusion term is related to the two-point two-time Eulerian acceleration-acceleration correlation. Dimensional analysis yields a velocity increment probability distribution with normal scaling v approximately t(1/2). However, the statistics need not be Gaussian.     

Mean velocity profile in confined turbulent convection

In this Letter, highly resolved measurements of the horizontal velocity inside the boundary layer of turbulent Rayleigh-Bénard convection are reported. They were performed in a cylindrical box with an aspect ratio Gamma=1.13 which was filled with air with a Prandtl number Pr=0.7. The horizontal velocity was measured along the central axis close to the cooling plate in a range of Rayleigh numbers between Ra=10;{11} and Ra=10;{12} using a two-dimensional laser Doppler velocimeter. We demonstrate that the profile of the mean velocity strongly differs from that of classical shear flows like the Blasius shape of a laminar flat plate boundary layer or a turbulent logarithmic velocity profile with standard coefficients.     

Use of Lagrangian statistics for the analysis of the scale separation hypothesis in turbulent channel flow

Turbulence models often involve Reynolds averaging, with a closure providing the Reynolds stress tensor as function of mean velocity gradients, through a turbulence constitutive equation. The main limitation of this linear closure is that it rests on an analogy with kinetic theory. For this analogy to be valid there has to be a scale separation between the mean velocity variations and the turbulent Lagrangian free path whose mean value is the turbulent mixing length. The aim of this work is to better understand this hypothesis from a microscopic point of view. Therefore, fluid elements are tracked in a turbulent channel flow. The flow is resolved by direct numerical simulation (DNS). Statistics on particle trajectories ending on a certain distance y₀ from the wall are computed, leading to estimations of the turbulent mixing length scale and the Knudsen number. Comparing the computed values to the Knudsen number in the case of scale separation, we may know in which region of the flow and to what extent the turbulence constitutive equation is not verified. Finally, a new non-local formulation for predicting the Reynolds stress is proposed.     

Multifractal statistics of Lagrangian velocity and acceleration in turbulence

The statistical properties of velocity and acceleration fields along the trajectories of fluid particles transported by a fully developed turbulent flow are investigated by means of high resolution direct numerical simulations. We present results for Lagrangian velocity structure functions, the acceleration probability density function, and the acceleration variance conditioned on the instantaneous velocity. These are compared with predictions of the multifractal formalism, and its merits and limitations are discussed.     

Broadband sound generation by confined turbulent jets

Sound generation by confined stationary jets is of interest to the study of voice and speech production, among other applications. The generation of sound by low Mach number, confined, stationary circular jets was investigated. Experiments were performed using a quiet flow supply, muffler-terminated rigid uniform tubes, and acrylic orifice plates. A spectral decomposition method based on a linear source-filter model was used to decompose radiated nondimensional sound pressure spectra measured for various gas mixtures and mean flow velocities into the product of (1) a source spectral distribution function; (2) a function accounting for near field effects and radiation efficiency; and (3) an acoustic frequency response function. The acoustic frequency response function agreed, as expected, with the transfer function between the radiated acoustic pressure at one fixed location and the strength of an equivalent velocity source located at the orifice. The radiation efficiency function indicated a radiation efficiency of the order (kD)2 over the planar wave frequency range and (kD)4 at higher frequencies, where k is the wavenumber and D is the tube cross sectional dimension. This is consistent with theoretical predictions for the planar wave radiation efficiency of quadrupole sources in uniform rigid anechoic tubes. The effects of the Reynolds number, Re, on the source spectral distribution function were found to be insignificant over the range 20002.5. The influence of a reflective open tube termination on the source function spectral distribution was found to be insignificant, confirming the absence of a feedback mechanism.     

Lagrangian acceleration of passive tracers in statistically steady rotating turbulence

The statistical properties of the Lagrangian acceleration vector of passive tracers in statistically steady rotating turbulence is studied by particle tracking velocimetry. Direct effects of the background rotation are the suppression of high-acceleration events parallel to the (vertical) rotation axis, the enhancement of high-acceleration events for the horizontal acceleration, and the strong amplification of the autocorrelation of the acceleration component perpendicular to both the rotation vector Ω and local velocity vector u. The autocorrelation of the acceleration component in the plane set up by Ω and u is only mildly enhanced.     

Size scaling of turbulent transport in magnetically confined plasmas

Transport scaling with respect to device size in magnetically confined plasmas is critically examined for electrostatic ion-temperature-gradient turbulence using global gyrokinetic particle simulations. It is found, by varying device size normalized by ion gyroradius while keeping other dimensionless plasma parameters fixed, that fluctuation scale length is microscopic in the presence of zonal flows. The local transport coefficient exhibits a gradual transition from a Bohm-like scaling for device sizes corresponding to present-day experiments to a gyro-Bohm scaling for future larger devices.     

Lagrangian velocity statistics in turbulent flows: effects of dissipation

We use the multifractal formalism to describe the effects of dissipation on Lagrangian velocity statistics in turbulent flows. We analyze high Reynolds number experiments and direct numerical simulation data. We show that this approach reproduces the shape evolution of velocity increment probability density functions from Gaussian to stretched exponentials as the time lag decreases from integral to dissipative time scales. A quantitative understanding of the departure from scaling exhibited by the magnitude cumulants, early in the inertial range, is obtained with a free parameter function D(h) which plays the role of the singularity spectrum in the asymptotic limit of infinite Reynolds number. We observe that numerical and experimental data are accurately described by a unique quadratic D(h) spectrum which is found to extend from h(min) approximately 0.18 to h(max) approximately 1.     

Joint statistics of the Lagrangian acceleration and velocity in fully developed turbulence

We report experimental results on the joint statistics of the Lagrangian acceleration and velocity in highly turbulent flows. The acceleration was measured up to a microscale Reynolds number R(lambda)=690 using high speed silicon strip detectors from high energy physics. The acceleration variance was observed to be strongly dependent on the velocity, following a Heisenberg-Yaglom-like u(9/2) increase. However, the shape of the probability density functions of the acceleration component conditioned on the same component of the velocity when normalized by the acceleration variance was observed to be independent of velocity and to coincide with the unconditional probability density function of the acceleration components. This observation imposes a strong mathematical constraint on the possible functional form of the acceleration probability distribution function.     

Ultimate turbulent Taylor-Couette flow

The flow structure of strongly turbulent Taylor-Couette flow with Reynolds numbers up to Re(i)=2×10(6) of the inner cylinder is experimentally examined with high-speed particle image velocimetry (PIV). The wind Reynolds numbers Re(w) of the turbulent Taylor-vortex flow is found to scale as Re(w)∝Ta(1/2), exactly as predicted by Grossmann and Lohse [Phys. Fluids 23, 045108 (2011).] for the ultimate turbulence regime, in which the boundary layers are turbulent. The dimensionless angular velocity flux has an effective scaling of Nu(ω)∝Ta(0.38), also in correspondence with turbulence in the ultimate regime. The scaling of Nu(ω) is confirmed by local angular velocity flux measurements extracted from high-speed PIV measurements: though the flux shows huge fluctuations, its spatial and temporal average nicely agrees with the result from the global torque measurements.     

Turbulent vortex rings in weakly turbulent and in turbulent ambient flow

Turbulent vortex rings were investigated in weakly turbulent flow and in three different grid generated turbulent flows to clarify the reciprocal action of the vortex ring with defined external turbulence. Assuming self-similarity and turbulent viscosity as proportional to V₀D₀ the equations for the ring diameter D(t) and the velocity of propagation V(t) were derived. The time difference Δt between the virtual origins of 1/V(t) respectively D²(t) led to an invariant term. The equation of momentum is fulfilled. – Position and diameter of the vortex rings were determined till their decay by means of an optical system, which did not disturb the vortex rings. The experimental results in weakly turbulent ambient flow obtained by the author and by others confirm the theory very well. The ambient turbulence was nearly constant in the measuring region; its effect could be described by simply adding its viscosity to the vortex ring’s internal turbulent viscosity. The results could be represented in unified non-dimensional diagrams. Moreover, an explanation was found as to why the mean internal turbulent viscosity is constant.     

Layering instability in a confined suspension flow

We show [J. Fluid Mech. 592, 447 (2007)] that swapping (reversing) trajectories in confined suspension flows prevent collisions between particles approaching each other in adjacent streamlines. Here we demonstrate that by inducing layering this hydrodynamic mechanism changes the microstructure of suspensions in a confined Couette flow. Layers occur either in the near-wall regions or span the whole channel width, depending on the strength of the swapping-trajectory effect. While our theory focuses on dilute suspensions, we postulate that this new hydrodynamic mechanism controls the formation of a layered microstructure in a wide range of densities.     

Lagrangian investigation of local extinction,re-ignition and auto-ignition in turbulent flames

Lagrangian PDF investigations are performed of the Sandia piloted flame E and the Cabra H₂/N₂ lifted flame to help develop a deeper understanding of local extinction, re-ignition and auto-ignition in these flames, and of the PDF models' abilities to represent these phenomena. Lagrangian particle time series are extracted from the PDF model calculations and are analyzed. In the analysis of the results for flame E, the particle trajectories are divided into two groups: continuous burning and local extinction. For each group, the trajectories are further sub-divided based on the particles' origin: the fuel stream, the oxidizer stream, the pilot stream, and the intermediate region. The PDF calculations are performed using each of three commonly used models of molecular mixing, namely the EMST, IEM and modified Curl mixing models. The calculations with different mixing models reproduce the local extinction and re-ignition processes observed in flame E reasonably well. The particle behavior produced by the IEM and modified Curl models is different from that produced by the EMST model, i.e., the temperature drops prior to (and sometimes during) re-ignition. Two different re-ignition mechanisms are identified for flame E: auto-ignition and mixing-reaction. In the Cabra H₂/N₂ lifted flame, the particle trajectories are divided into different categories based on the particles' origin: the fuel stream, the oxidizer stream, and the intermediate region. The calculations reproduce the whole auto-ignition process reasonably well for the Cabra flame. Four stages of combustion in the Cabra flame are identified in the calculations by the different mixing models, i.e., pure mixing, auto-ignition, mixing-ignition, and fully burnt, although the individual particle behavior by the IEM and modified Curl models is different from that by the EMST model. The relative importance of mixing and reaction during re-ignition and auto-ignition are quantified for the IEM model.