Modulation of homogeneous turbulence seeded with finite size bubbles or particles

The dynamics of homogeneous, isotropic turbulence seeded with finite sized particles or bubbles is investigated in a series of numerical simulations, using the force-coupling method for the particle phase and low wavenumber forcing of the flow to sustain the turbulence. Results are given on the modulation of the turbulence due to massless bubbles, neutrally buoyant particles and inertial particles of specific density 1.4 at volumetric concentrations of 6%. Buoyancy forces due to gravity are excluded to emphasize finite size and inertial effects for the bubbles or particles and their interactions with the turbulence. Besides observing the classical entrapment of bubbles and the expulsion of inertial particles by vortex structures, we analyze the Lagrangian statistics for the velocity and acceleration of the dispersed phase. The turbulent fluctuations are damped at mid-range wavenumbers by the bubbles or particles while the small-scale kinetic energy is significantly enhanced. Unexpectedly, the modulation of turbulence depends only slightly on the dispersion characteristics (bubble entrapment in vortices or inertial sweeping of the solid particles) but is closely related to the stresslet component (finite size effect) of the flow disturbances. The pivoting wavenumber characterizing the transition from damped to enhanced energy content is shown to vary with the size of the bubbles or particles. The spectrum for the energy transfer by the particle phase is examined and the possibility of representing this, at large scales, through an additional effective viscosity is discussed.     

Measurement of inertial particle clustering and relative velocity statistics in isotropic turbulence using holographic imaging

We present the first measurements of relative velocity statistics of inertial particles in a homogeneous isotropic turbulent flow with three-dimensional holographic particle image velocimetry (holographic PIV). From the measurements we are able to obtain the radial relative velocity probability density function (PDF) conditioned on the interparticle separation distance, for distances on the order of the Kolmogorov length scale. Together with measurements of the three-dimensional radial distribution function (RDF) in our turbulence chamber, these statistics, in principle, can be used to determine interparticle collision rates via the formula derived by Sundaram and Collins (1997). In addition, we show temporal development of the RDF, which reveals the existence of an extended quasi-steady-state regime in our facility. Over this regime the measured two-particle statistics are compared to direct numerical simulations (DNS) with encouraging qualitative agreement. Statistics at the same Reynolds number but different Stokes numbers demonstrate the ability of the experiment to correctly capture the trends associated with particles of different inertia. Our results further indicate that even at moderate Stokes numbers turbulence may enhance collision rates significantly. Such experimental investigations may prove valuable in validating, guiding and refining numerical models of particle dynamics in turbulent flows.     

A hybrid Lagrangian–Eulerian particle‐level set method for numerical simulations of two‐fluid turbulent flows

A coupled Lagrangian interface‐tracking and Eulerian level set (LS) method is developed and implemented for numerical simulations of two‐fluid flows. In this method, the interface is identified based on the locations of notional particles and the geometrical information concerning the interface and fluid properties, such as density and viscosity, are obtained from the LS function. The LS function maintains a signed distance function without an auxiliary equation via the particle‐based Lagrangian re‐initialization technique. To assess the new hybrid method, numerical simulations of several ‘standard interface‐moving’ problems and two‐fluid laminar and turbulent flows are conducted. The numerical results are evaluated by monitoring the mass conservation, the turbulence energy spectral density function and the consistency between Eulerian and Lagrangian components. The results of our analysis indicate that the hybrid particle‐level set method can handle interfaces with complex shape change, and can accurately predict the interface values without any significant (unphysical) mass loss or gain, even in a turbulent flow. The results obtained for isotropic turbulence by the new particle‐level set method are validated by comparison with those obtained by the ‘zero Mach number’, variable‐density method. For the cases with small thermal/mass diffusivity, both methods are found to generate similar results. Analysis of the vorticity and energy equations indicates that the destabilization effect of turbulence and the stability effect of surface tension on the interface motion are strongly dependent on the density and viscosity ratios of the fluids. Copyright © 2007 John Wiley & Sons, Ltd.     

Filtered particle tracking in isotropic turbulence and stochastic modeling of subgrid-scale dispersion

A numerical study based on the Eulerian–Lagrangian formulation is performed for dispersed phase motion in a turbulent flow. The effect of spatial filtering, commonly employed in large-eddy simulations, and the role of the subgrid scale turbulence on the statistics of heavy particles, including preferential concentration, are studied through a priori analysis of DNS of particle-laden forced isotropic turbulence. In simulations where the subgrid scale kinetic energy attains 30–35% of the total we observe the impact of residual fluid motions on particles of a smaller inertia. It is shown that neglecting the influence of subgrid scale fluctuations has a significant effect on the preferential concentration of those particles. A stochastic Langevin model is proposed to reconstruct the residual (or subgrid scale) fluid velocity along particle trajectories. The computation results for a selection of particle inertia parameters are performed to appraise the model through comparisons of particle turbulent kinetic energy and the statistics of preferential concentrations.     

Turbulence-induced preferential concentration of solid particles in microgravity conditions

A box of near-isotropic, particle-laden turbulence was flown aboard NASA's reduced gravity aircraft in order to measure the turbulence-induced preferential concentration of solid particles in microgravity. Three particle sets of Stokes numbers based on the fluid Kolmogorov time scale of approximately 0.5, 5, and 50 were tested for relative amounts of preferential concentration. Eight fans in each corner of a Lexan box generated fluid turbulence. Particle concentrations were measured using an imaging system consisting of a camera viewing perpendicular to a white light sheet. Post-processing of video images found largest concentration gradients for the intermediate-sized particles of Stokes number 5, closely followed by the Stokes number 0.5 particles. The experimental results agreed well with the trends seen in direct numerical simulations. The quantitative effects of turbulence modulation by the presence of particles were not measured in the experiment, but were most likely present. Received: 10 February 2000 / Accepted: 9 November 2001     

Particle‐tracking within turbulent cylindrical electromagnetically‐driven flow

This paper discusses a simple model of electromagnetically stirred molten metal within a long cylinder, neglecting end effects. The radially‐dependent velocity profiles of the molten metal are calculated using the Prandtl mixing length model for turbulence. The paths of non‐conducting particles within the fluid are also computed numerically, under the influence of Stokes’ drag and a random force due to turbulence. The paths are calculated for a range of particle diameters and the general motion is observed. Copyright © 1999 John Wiley & Sons, Ltd.     

SPH modelling of depth‐limited turbulent open channel flows over rough boundaries

A numerical model based on the smoothed particle hydrodynamics method is developed to simulate depth‐limited turbulent open channel flows over hydraulically rough beds. The 2D Lagrangian form of the Navier–Stokes equations is solved, in which a drag‐based formulation is used based on an effective roughness zone near the bed to account for the roughness effect of bed spheres and an improved sub‐particle‐scale model is applied to account for the effect of turbulence. The sub‐particle‐scale model is constructed based on the mixing‐length assumption rather than the standard Smagorinsky approach to compute the eddy‐viscosity. A robust in/out‐flow boundary technique is also proposed to achieve stable uniform flow conditions at the inlet and outlet boundaries where the flow characteristics are unknown. The model is applied to simulate uniform open channel flows over a rough bed composed of regular spheres and validated by experimental velocity data. To investigate the influence of the bed roughness on different flow conditions, data from 12 experimental tests with different bed slopes and uniform water depths are simulated, and a good agreement has been observed between the model and experimental results of the streamwise velocity and turbulent shear stress. This shows that both the roughness effect and flow turbulence should be addressed in order to simulate the correct mechanisms of turbulent flow over a rough bed boundary and that the presented smoothed particle hydrodynamics model accomplishes this successfully. © 2016 The Authors International Journal for Numerical Methods in Fluids Published by John Wiley & Sons Ltd     

Buffeting in transonic flow prediction using time‐dependent turbulence model

In transonic flow conditions, the shock wave/turbulent boundary layer interaction and flow separations on wing upper surface induce flow instabilities, ‘buffet’, and then the buffeting (structure vibrations). This phenomenon can greatly influence the aerodynamic performance. These flow excitations are self‐sustained and lead to a surface effort due to pressure fluctuations. They can produce enough energy to excite the structure. The objective of the present work is to predict this unsteady phenomenon correctly by using unsteady Navier–Stokes‐averaged equations with a time‐dependent turbulence model based on the suitable (k–ε) turbulent eddy viscosity model. The model used is based on the turbulent viscosity concept where the turbulent viscosity coefficient (Cμ) is related to local deformation and rotation rates. To validate this model, flow over a flat plate at Mach number of 0.6 is first computed, then the flow around a NACA0012 airfoil. The comparison with the analytical and experimental results shows a good agreement. The ONERA OAT15A transonic airfoil was chosen to describe buffeting phenomena. Numerical simulations are done by using a Navier–Stokes SUPG (streamline upwind Petrov–Galerkin) finite‐element solver. Computational results show the ability of the present model to predict physical phenomena of the flow oscillations. The unsteady shock wave/boundary layer interaction is described. Copyright © 2005 John Wiley & Sons, Ltd.     

Hyperviscosity and Vorticity-Based Models for Subgrid Scale Modeling

Large-eddy simulations corresponding to the decaying isotropic turbulence experiment of Comte-Bellot and Corrsin are performed, using a pseudo-spectral code that incorporates four models: viscosity and hyperviscosity types, each implemented for both the subgrid scale stress tensor and the subgrid scale force. Two 1/T scalings are also considered for the viscosity amplitude. The dynamic procedure is extended to the four models and is tested. Results are obtained with and without this procedure and for both scalings. The main conclusions are: (a) the two viscosity models perform equally well; (b) the Kolmogorov scaling performs as well as the Smagorinsky scaling, yet it is computationally more efficient; (c) in the dynamic procedure, there is a fairly wide range of test to grid filter ratios which produces results insensitive to this ratio; and (d) the hyperviscosity models lead to energy decay curves that follow the experimental data as well as the usual viscosity models.     

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.     

Rotational dynamics of bottom-heavy rods in turbulence from experiments and numerical simulations

We successfully perform the three-dimensional tracking in a turbulent fluid flow of small axisymmetrical particles that are neutrally-buoyant and bottom-heavy, i.e., they have a non-homogeneous mass distribution along their symmetry axis. We experimentally show how a tiny mass inhomogeneity can affect the particle orientation along the preferred vertical direction and modify its tumbling rate. The experiment is complemented by a series of simulations based on realistic Navier–Stokes turbulence and on a point-like particle model that is capable to explore the full range of parameter space characterized by the gravitational torque stability number and by the particle aspect ratio. We propose a theoretical perturbative prediction valid in the high bottom-heaviness regime that agrees well with the observed preferential orientation and tumbling rate of the particles. We also show that the heavy-tail shape of the probability distribution function of the tumbling rate is weakly affected by the bottom-heaviness of the particles.     

各向同性湍流内颗粒碰撞率的直接模拟研究

对 Re_{\lambda } 约为51均匀各向同性湍流内 St_{k}(=\tau_{p}/\tau_{k}) 为 0 ~10.0 的 有限惯性颗粒的碰撞行为进行了直接数值模拟，以研究湍流对有限惯性 颗粒碰撞的影响. 结果表明，具有一定惯性颗粒的湍流碰撞率完全不同于零惯性的轻颗粒 (St_{k}=0) 和可忽略湍流作用的重颗粒 (St{k} \to \infty) , 其变化趋势极其复杂： 在Stk为 0~1.0 之间，颗粒的碰撞率随 St 的增加而近乎线性地剧烈增长，在 Stk≈1.0 3.0(对应的StE=τp/Te≈0.5)附近，颗粒碰撞率出现两个峰值，在Stk>3.0以后，颗粒的碰撞率随惯性增 大而逐渐趋向于重颗粒极限；在峰值处，有限惯性颗粒的平均碰撞率的峰值较轻颗粒增强了 30倍左右. 为进一步分析湍流作用下颗粒碰撞率的影响因素，分别使用可能发生碰撞 的颗粒对的径向分布函数和径向相对速度来量化颗粒的局部富集效应和湍流掺混效应，表明 St_{k} \approx 1.0 时局部富集效应最为强烈，使得颗粒的碰撞率出现第1个峰值； 湍流掺混效应则随着颗粒Stk的增大而渐近增大；局部富集和湍流掺混联合作用的结果， 使得颗粒碰撞率在 St_{k} \approx 3.0 附近出现另一个峰值.     

VELOCITY-ACCELERATION STRUCTURE FUNCTION IN TWO-DIMENSIONAL DECAYING TURBULENCE$^{\bf 1)}$

湍流场中二阶速度加速度结构函数 (velocity-acceleration structure function, VASF) 被认为与尺度间能量或者拟涡能的传递相关,其正负表明传递的方向. 三维湍流中,能量从大尺度向 小尺度传递,VASF 为负. 在二维湍流中,能量反向传递到大尺度,拟涡能正向传递到小尺度,因此理论上 VASF 无论在反向能量级串区还是在正向拟 涡能级串区均为正. 然而,相对于三维湍流中 VASF 的充分研究,二维湍流中 VASF 的正负性迄今尚无实验或数值模拟数据验证. 本文通过三维二维湍流中普适的公式推导,指出在空间非均匀湍流场中,VASF 除了尺度间传递,还受到非均匀项的影响. 一种常见的空间非均匀湍流场是在实验研究中常用的风洞或水洞中,湍流发生装置 (如栅格) 后的湍流. 该流场中,湍流强度随下游位置增大而逐渐衰减,这种衰减则带来空间上的非均匀性. 本文在基于竖直流动皂膜的二维衰减湍流场中,利用拉格朗日粒子追踪法测得在拟涡能级串区的 VASF,并分析各部分的影响. 结果表明,虽然尺度间传递项为正值,但由于衰减带来的非均匀项为负值,使得 VASF 的值为负,使之失去了表征拟涡能传递方向的意义. 因此,在类似风洞、水洞、水槽等衰减流场中对 VASF 的讨论不应忽略非均匀项. 最后对与 VASF 密切相关的弥散过程进行分析,发现后期弥散过程变慢是由于负的 VASF 导致.     

液固两相槽道流中湍流调制的数值研究

结合雷诺应力模型 (Reynolds Stress Model, RSM) 和混合模型 (Mixture Model) 对槽道湍流向下流动中的颗粒调制湍流问题进行了研究．该模型考虑了颗粒流的动能理论和颗粒对湍流的反馈作用．着重分析了颗粒对湍流的调制作用,以及颗粒参数变量（如颗粒密度和质量载荷）对湍流调制的影响．结果表明：(1)在颗粒抑制湍流的范围内,当颗粒密度小于载流体密度时,湍流强度的改变量与颗粒密度成反比;当颗粒密度大于载流体密度时,湍流强度的改变量与颗粒密度成正比;(2)在一定范围内,颗粒抑制湍流的能力随颗粒质量载荷增加而变强．     

The effect of two-way coupling and inter-particle collisions on turbulence modulation in a vertical channel flow

Turbulence modulation due to its interaction with dispersed solid particles in a downward fully developed channel flow was studied. The Eulerian framework was used for the gas-phase, whereas the Lagrangian approach was used for the particle-phase. The steady-state equations of conservation of mass and momentum were used for the gas-phase, and the effect of turbulence on the flow-field was included via the standard k–ε model. The particle equation of motion included the drag, the Saffman lift and the gravity forces. Turbulence dispersion effect on the particles was simulated as a continuous Gaussian random field. The effects of particles on the flow were modeled by appropriate source terms in the momentum, k and ε equations. Particle–particle collisions and particle–wall collisions were accounted for in these simulations. Gas-phase velocities and turbulence kinetic energy in the presence of 2–100% mass loadings of two particle classes (50 μm glass and 70 μm copper) were evaluated, and the results were compared with the available experimental data and earlier numerical results. The simulation results showed that when the inter-particle collisions were important and was included in the computational model, the fluid turbulence was attenuated. The level of turbulence attenuation increased with particle mass loading, particle Stokes number, and the distance from the wall. When the inter-particle collisions were negligible and/or was neglected in the model, the fluid turbulence was augmented for the range of particle sizes considered.     

On the nature of multitude scalings in decaying isotropic turbulence

We report multitude scaling laws for isotropic fully developed decaying turbulence through group theoretic method employing on the spectral equations both for modelling and without any modelling of nonlinear energy transfer. For modelling, besides the existence of classical power law scalings, an exponential decay of turbulent energy in time is obtained subject to exponentially decaying integral length scale at infinite Reynolds number limit. For the transfer without modelling, at finite Reynolds number, in addition to general power law decay of turbulence intensity with integral length scale growing as a square root of time, an exponential decay of energy in time is explored when integral length scale remains constant. Both the power and exponential decaying laws of energy agree to the theoretical results of George (1992), George and Wang (2009) and experimental results of fractal grid generated turbulence by Hurst and Vassilicos (2007). At infinite Reynolds number limit, a general power law scaling is obtained from which all classical scaling laws are recovered. Further, in this limit, turbulence exhibits a general exponential decaying law of energy with exponential decaying integral length scale depending on two scaling group parameters. The role of symmetry group parameters on turbulence dynamics is discussed in this study.     

Dispersion of particles in wall-bounded particle-laden turbulent flows with high wall permeability

A series of numerical simulations were performed to investigate the distribution and deposition properties of particles in turbulent flows bounded by permeable walls using the Large Eddy Simulation (LES) with a Lagrangian trajectory approach. The wall permeation speeds were taken from 10⁻⁴ to 10⁻² of the bulk velocity. The directions of the permeation speed were the same at both walls, and they were inward on one wall but outward on the other wall to reserve the fluid mass. Particles with Stokes number (respecting viscous time scale) around 0.1, 1 and 10 were released in the fully developed turbulent channel flow. The particle–particle interaction and the retroaction from particles to the fluid were neglected. The fluid-phase turbulence statistical properties and particle's transport characteristics by vortexes were then analyzed in details. If the wall permeation exists, the turbulence intensities will be depressed close to the outward permeable wall but increased near the inward permeable wall. Not influenced by the wall permeation, the suspended particles with St⁺ ∼O(1) tend to accumulate in the less vortical zones away from the wall, while those particles in the flow regions near the outward permeable wall will distribute disregarding of the vorticity. The turbulence structures near the outward permeable wall are found to exert promotional effects on the particle deposition rate, but such effects are different for particles with various Stokes number. A distribution tendency of streamwise streaks for the deposited particles is also found on the wall imposed by the high outward permeation speed and the clustering deposition pattern is more obvious with increasing particle Stokes number.     

Direct numerical simulation of a particle-laden low Reynolds number turbulent round jet

Direct numerical simulation method is used for the investigating of particle-laden turbulent flows in a spatially evolution of low Reynolds number axisymmetric jet, and the Eulerian–Lagrangian point-particle approach is employed in the simulation. The simulation uses an explicit coupling scheme between particles and the fluid, which considers two-way coupling between the particle and the fluid. The DNS results are compared well with experimental data with equal Reynolds number (R_e = 1700). Our objects are: (i) to investigate the correlation between the particle number density and the fluctuating of fluid streamwise velocity; (ii) to examine whether the three-dimensional vortex structures in the particle-laden jet are the same as that in the free-air jet and how the particles modulate the thee-dimensional vortex structures and turbulence properties with different Stokes number particles; (iii) to discover the particle circumferential dispersion with different Stokes number particles. Our findings: (i) all the particles, regardless of their particle size, tend to preferentially accumulate in the region with large-than-mean fluid streamwise velocity; (ii) the small Stokes number particles take an important part in the modulation of three-dimensional vortex structures, but for the intermediate and larger sized particles, this modulation effect seems not so apparent; (iii) the particle circumferential dispersion is more effective for the smaller and intermediate sized particles, especially for the intermediate sized particles.     

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.     

Lagrangian velocity and acceleration statistics of fluid and inertial particles measured in pipe flow with 3D particle tracking velocimetry

Three-dimensional particle tracking velocimetry (3D-PTV) has been applied to particle-laden pipe flow at Reynolds number 10,300, based on the bulk velocity and the pipe diameter. The volume fraction of the inertial particles was equal to 1.4 × 10⁻⁵. Lagrangian velocity and acceleration statistics were determined both for tracers and for inertial particles with Stokes number equal to 2.3, based on the particle relaxation time and the viscous time scale. The decay of Lagrangian velocity and acceleration correlation functions was measured both for the fluid and for the dispersed phase at various radial positions. The decay of Lagrangian velocity correlations is faster for inertial particles than for flow tracers, whereas the decay of Lagrangian acceleration correlations is about 25% slower for inertial particles than for flow tracers. Further differences between inertial and tracer particles are found in velocity fluctuations evaluated for both positive and negative time lags. The asymmetry in time of velocity cross-correlations is more pronounced for inertial particles. Quadrant analysis revealed another difference still near the wall: ejection and sweep events are less frequent for inertial particles than for tracers.