Turbulence modulation in dilute particle-laden flow

A new particle source term to account for the effect of particles on the turbulence equations based on the Euler/Lagrange approach is introduced and compared with existing models and experimental data. Three different sizes of particles are considered to cover the range of large particles, where augmentation of the carrier phase turbulence is expected, and small particles, for which attenuation is expected. The new model is derived directly from the balance equations of fluid flow and represents a combination of the so-called standard and consistent approaches. The performance of the new model surpasses that of the standard and consistent models and it is able to predict both the suppression and enhancement of fluid turbulence for small and large particles.     

Turbulence production in flow over a wavy wall

Measurements of the spatial and time variation of two components of the velocity have been made over a sinusoidal solid wavy boundary with a height to length ratio of 2a/λ = 0.10 and with a dimensionless wave number of α⁺ = (2π/λ)(v/u ^?) = 0.02. For these conditions, both intermittent and time-mean flow reversals are observed near the troughs of the waves. Statistical quantities that are determined are the mean streamwise and normal velocities, the root-meansquare of the fluctuations of the streamwise and normal velocities, and the Reynolds shear stresses. Turbulence production is calculated from these measurements. The flow is characterized by an outer flow and by an inner flow extending to a distance of about α^?1 from the mean level of the surface. Turbulence production in the inner region is fundamentally different from flow over a flat surface in that it is mainly associated with a shear layer that separates from the back of the wave. Flow close to the surface is best described by an interaction between the shear layer and the wall, which produces a retarded zone and a boundary-layer with large wall shear stresses. Measurements of the outer flow compare favorably with measurements over a flat wall if velocities are made dimensionless by a friction velocity defined with a shear stress obtained by extrapolating measurements of the Reynolds stress to the mean levels of the surface (rather than from the drag on the wall).     

Turbulence statistics from optical whole-field measurements in particle-laden turbulence

Obtaining turbulence statistics in particle-laden flows using optical whole-field measurements is complicated due to the inevitable data loss. The effects of this data loss are first studied using synthetic data and it is shown that the interpolation of missing data leads to biased results for the turbulence spectrum and its derived quantities. It is also shown that the use of overlapping interrogation regions in images with a low image density can lead to biased results due to oversampling. The slotting method is introduced for the processing of particle image velocimetry (PIV) data fields with missing data. Next to this, it is extended to handle unstructured data. Using experimental data obtained by a dual-camera PIV/PTV (particle tracking velocimetry) system in particle-laden grid turbulence, the performance of the new approach is studied. Some preliminary two-phase results are presented to indicate the significant improvement in the statistics, as well as to demonstrate the unique capabilities of the system.An erratum to this article can be found at

---

Turbulence measurements in a flow generated by the collision of radially flowing wall jets

Early results of an experimental investigation of the abnormally high turbulence level and mixing layer growth rate characteristics found in the upwash regions of aircraft with vertical short takeoff and landing (V/STOL) flows in ground effect are presented. The upwash flow is formed from the collision of two opposing radially flowing wall jets. The wall jets are created in a unique way that allows the upwash to form without any interference due to the source jets. The objective of this work is to systematically characterize the development and structure of the flow. The upwash flow exhibits very large mixing rates compared to turbulent free or wall jet flows. A unique set of two component velocity profiles was taken in the upwash flow field. These measurements include several higher moment terms that appear in the turbulent kinetic energy equations, as well as length scales and intermittency determinations. Measurements were taken' along the axis connecting the two source jets as well as off this axis at six measurement stations above ground. The results provide detailed data on an important class of flows where none existed, and they are expected to significantly improve the computational empirical tools available for predicting V/STOL behavior near the ground.A version of this paper was presented at the 10th Symposium on Turbulence, University of Missouri-Rolla, September 22–24, 1986     

Laminar flow performance of a heated body in particle-laden water

The effects of small uniformly sized spherical particles seeded into the freestream flow of a water tunnel on the delayed transition of a heated laminar flow control body is examined experimentally. In separate trials, four different mean diameter particle seedings were added to the flow and the approach flow velocity was cycled from subcritical to supercritical conditions at three different body heating conditions. The transition Reynolds number based on the body arc length and the approach flow velocity decreases monotonically with increasing d/ ^*, where d is the particle diameter and ^* is the displacement thickness at a critical location. The location of initial turbulent spot formation defines the critical location, and, within the range of experimental conditions reported here, is independent of particle size, heating condition and the approach velocity. For the high unit Reynolds numbers considered (Re ^u 1.88 × 10⁷ per metre), there is no observed critical particle diameterbased Reynolds number threshold; all sizes of particles considered in the experiments (d = 37 to 218 m) have some effect on transition. In a second set of experiments, particles were injected into the laminar boundary layer from a small orifice located at the forward stagnation point. These injected particles have no observable effect on the laminar layer or transition, which suggests that the injected particles fail to produce wakes or vorticity within the laminar layer that may lead to turbulent spot production.Also with the Graduate Program in Acoustics, Penn State UniversityThis work has been supported by the Applied Research Laboratory of The Pennsylvania State University under contracts with the Office of Naval Research and the Naval Sea Systems Command. The authors are particularly indebted to Professor Ron Blackwelder and his colleagues for sharing their yet unpublished findings from particle-induced transition experiments being conducted at the University of Southern California.     

Turbulence modulation in two-phase jets

A simple model is presented which can be used to predict the modifying effect of a dispersed-phase on the turbulence structure of two-phase jets. It is based on Kolmogorov's concept of spectral energy transfer and takes into account the additional energy dissipation resulting from the inability of dispersed-phase particles to completely follow turbulent eddy fluctuations. According to the analysis presented, high-frequency eddies are attenuated preferentially and a reduction of the total rate of turbulent energy dissipation results. Turbulence intensities were also reduced. Good agreement between predictions and experimental findings were obtained.     

Transient growth in turbulent particle-laden channel flow

Linear transient growth of optimal perturbations in particle-laden turbulent channel flow is investigated in this work.The problem is formulated in the framework of a Eulerian-Eulerian approach,employing two-way coupling between fine particles and fluid flow.The model is first validated in laminar cases,after which the transient growth of coherent perturbations in turbulent channel flow is investigated,where the mean particle concentration distribution is obtained by direct numerical simulation.It is shown that the optimal small-scale structures for particles are streamwise streaks just below the optimal streamwise velocity streaks,as was previously found in numerical simulations of particle-laden channel flow.This indicates that the optimal growth of perturbations is a dominant mechanism for the distribution of particles in the near-wall region.The current study also considers the transient growth of small-and large-scale perturbations at relatively high Reynolds numbers,which reveals that the optimal large-scale structures for particles are in the near-wall region while the optimal large-scale structures for fluid enter the outer region.     

Lagrangian particle dispersion in turbulent flow over a wall mounted obstacle

Large-eddy simulations (LES) of particle-laden turbulent flows are presented in order to investigate the effects of particle response time on the dispersion patterns of a space developing flow with an obstruction, where solid particles are injected inside the wake of an obstacle [Vincont, J.Y., Simoens, S., Ayrault M., Wallace, J.M., 2000. Passive scalar dispersion in a turbulent boundary layer from a line source at the wall and downstream of an obstacle. J. Fluid Mech. 424, 127–167]. The numerical method is based on a fully explicit fractional step approach and finite-differences on Cartesian grids, using the immersed boundary method (IBM) to represent the existence of solid obstacles. Two different turbulence models have been tested, the classical Smagorinsky turbulence model and the filtered structure function model. The dispersed phase was modelled either by an Eulerian approach or a Lagrangian particle tracking scheme of solid particles with Stokes numbers in the range St = 0–25, assuming one-way coupling between the two phases. A very good agreement was observed between the Lagrangian and Eulerian approaches. The effect of particle size was found to significantly differentiate the dispersion pattern for the inhomogeneous flow over the obstacle. Although in homogeneous flows like particle-laden turbulent channels near-wall particle clustering increases monotonically with particle size, for the examined flow over an obstacle, preferential concentration effects were stronger only for an intermediate range of Stokes numbers.     

Turbulence control in wall jets

A new method for the control of mixing of a plane turbulent wall jet has been investigated. A thin wire, mounted in the vicinity of the wall-jet nozzle, changes the formation of the shear-layer structures in the early stages of the development of the wall jet. The wire is operated in two ways: (1) a still wire inhibits the natural shear layer roll-up and reduces the size of the turbulent structures and thereby the mixing; (2) a self-excited oscillating wire introduces large coherent structures and thereby enhances the mixing. The size of these structures does not depend on the shear-layer instability but rather on the wire frequency.Measurements of the mean and fluctuating velocities have been performed with hot-wire anemometry as well as measurements of the skin friction by means of Preston tubes, surface fences and wall hot-wires. The Reynolds number based on the slot width was Re_j=10000.     

Temperature statistics in a radiatively heated particle-laden turbulent square duct flow

Radiation absorption by preferentially concentrated particles in a turbulent square duct flow is studied experimentally. The particle-laden flow is exposed to near-infrared radiation, and the gas phase temperature statistics are measured along the wall bisector of the duct. It is found that the instantaneous temperature fluctuations are comparable to the overall mean temperature rise. The temperature statistics at the duct centerline and near the wall are qualitatively different. The former reflects preferential concentration in isotropic flows while the latter displays evidence of particle clustering into streamwise elongated streaks. Comparison of the experimental data to a simplified heat transfer model suggests that the Lagrangian evolution of particle clusters and voids, and turbulent mixing in the vicinity of particle clusters, are important. This work was motivated by particle solar receiver technology, but the findings are also relevant to systems where there is localized heat release or mass transfer from disperse particles or droplets. It shows that obtaining Lagrangian histories of particle trajectories is an important next step towards understanding thermal transport phenomena in particle-laden turbulent flows.     

Turbulence measurements in a swirling pipe flow

This paper reports laser-Doppler measurements of the mean flow and turbulence stresses in a swirling pipe flow. Experiments were carried out under well-controlled laboratory conditions in a refractive index-matched pipe flow facility. The results show pronounced asymmetry in mean and fluctuating quantities during the downstream decay of the swirl. Experimental data reveal that the swirl significantly modifies the anisotropy of turbulence and that it can induce explosive growth of the turbulent kinetic energy during its decay. Anisotropy invariant mapping of the turbulent stresses shows that the additional flow deformation imposed by initially strong swirling motion forces turbulence in the core region to tend towards the isotropic two-component state. When turbulence reaches this limiting state it induces rapid production of turbulent kinetic energy during the swirl decay.

---

Turbulence in multiphase flow

We establish in this paper the foundations of a two-field turbulent flow model that includes two turbulent fields. The case of dispersed particles in an incompressible carrier fluid is treated here, but the very presence of these two fields allows for the generalization of the model to the instability-induced turbulent mixing of two materials. This model describes both cases of turbulent mass diffusion and small drag regime, “wave-like” interpenetration of the two components. It also includes the damping of the turbulence due to the presence of the particles. In addition, a theoretical derivation of the drag-induced decay of the large-scale turbulence kinetic energy is proposed as another mechanism specific to turbulent multiphase flow.     

Direct numerical simulation of a particle-laden flow in a flat plate boundary layer

In this paper, a direct numerical simulation of particle-laden flow in a flat plate boundary layer is performed, using the Eulerian–Lagrangian point-particle approach. This is, as far as we know, the first simulation of a particle-laden spatially-developing turbulent boundary layer with two-way coupling. A local minimum of the particle number density is observed in the close vicinity of the wall. The present simulation results indicate that the inertial particles displace the quasi-streamwise vortices towards the wall, which, in turn, enhance the mean streamwise fluid velocity. As a result, the skin-friction coefficient is increased whereas the boundary layer integral thicknesses are reduced. The presence of particles augments the streamwise fluctuating velocity in the near-wall region but attenuates it in the outer layer. Nevertheless, the wall-normal and spanwise velocity fluctuations are significantly damped, and so is the Reynolds stress. In addition, the combined effect of a reduced energy production and an increased viscous dissipation leads to the attenuation of the turbulent kinetic energy.     

Stochastic simulations of particle-laden isotropic turbulent flow

Numerical simulations are performed of dispersion and polydispersity of particles in isotropic incompressible turbulence. The mass loading of the particles is assumed to be small; thus the effects of particles on turbulence is neglected (one-way coupling). A stochastic model is employed to simulate the carrier phase. The results of the simulations are compared with direct numerical simulation (DNS) data and theoretical results. The stochastic model predicts most of the trends as portrayed by DNS and theory. However, the continuity effect associated with the crossing trajectories effect is not captured. Also, the peaking in the variation of the particle asymptotic diffusivity coefficient with the particle time constant is not observed. For evaporating particles, the stochastic model predicts thinner probability density functions (pdfs) for the particle diameter as compared with DNS generated pdfs. The model is implemented to investigate the effects of gravity on evaporation. It is shown that the depletion rate increases with increase of the drift velocity at short and intermediate times, but an opposite trend is observed at long times. The standard deviation and skewness of the particle diameter indicate peak values in their variations with the drift velocity. Dispersion of evaporating particles decreases with respect to that of non-evaporating particles at small drift velocities; an opposite trend is observed at large drift velocities. The effects of the initial evaporation rate and the particle Schmidt number on the evaporation in the gravity environment are also studied.     

Experimental investigation of flow over a backward-facing step in proximity to a flexible wall

Turbulent flow between a flexible wall and a solid surface containing a backward-facing step (BFS) was investigated using digital particle image velocimetry and high-speed photography. Stationary sheet of paper under tension was positioned above the solid surface in proximity to the BFS. The incoming air flow emerged from a planar nozzle that was located in the solid wall upstream of the BFS. Flows corresponding to two values of the Reynolds number (3,000 and 3,600) based on the step height and the maximum flow velocity at the step location were characterized in terms of patterns of time-averaged velocity, out-of-plane vorticity, streamline topology, and turbulence statistics. In addition, paper sheet oscillation was characterized using high-speed photography. For the control case of a solid upper wall with the geometry that represented the time-averaged paper profile, hydrodynamic frequencies were characterized using unsteady pressure measurements. Frequencies of the natural vibration modes of the paper sheet were well separated from the hydrodynamic frequencies corresponding to the oscillations of the shear layer downstream of the BFS. As the inflow velocity increased, the paper sheet was pulled closer to the solid surface, which resulted in increased confinement of the incoming jet. The flow reattachment length calculated on the basis of time-averaged flow patterns increased with the increasing Reynolds number.     

Turbulent boundary layer heat transfer for a constant property particle-laden gas flow

Solution of a turbulent boundary layer for a constant property, particle-laden gas flow is obtained by a differential method. A dimensionless analysis shows importance of an interaction parameter in increasing heat flux. Boundary layer analysis is done in usual manner by transforming partial differential equations and solution is started at the leading edge by Runge-Kutta method. Velocity and temperature profiles at downstream planes for gas and particles are calculated by an implicit finite-difference iterative procedure, and numerical results are compared with available experimental data.     

Large Eddy Simulations of turbulent particle-laden channel flow

This paper scrutinises the Large Eddy Simulation (LES) approach to simulate the behaviour of inter-acting particles in a turbulent channel flow. A series of simulations that are fully (four-way), two-way and one-way coupled are performed in order to investigate the importance of the individual physical phenomena occurring in particle-laden flows. Moreover, the soft sphere and hard sphere models, which describe the interaction between colliding particles, are compared with each other and the drawbacks and advantages of each algorithm are discussed. Different models to describe the sub-grid scale stresses with LES are compared. Finally, simulations accounting for the rough walls of the channel are compared to simulations with smooth walls. The results of the simulations are discussed with the aid of the experimental data of Kussin J. and Sommerfeld M., 2002, Experimental studies on particle behaviour and turbulence modification in horizontal channel flow with different wall roughness, Exp. in Fluids, 33, pp. 143–159 of Reynolds number 42,000 based on the full channel height. The simulations are carried out in a three-dimensional domain of 0.175 m × 0.035 m  × 0.035 m where the direction of gravity is perpendicular to the flow. The simulation results demonstrate that rough walls and inter-particle collisions have an important effect in redistributing the particles across the channel, even for very dilute flows. A new roughness model is proposed which takes into account the fact that a collision in the soft sphere model is fully resolved and it is shown that the new model is in very good agreement with the available experimental data.