Measurement of local size and velocity probability density distributions in two-phase suspension flows by Laser-Doppler technique

A technique is presented for the simultaneous measurement of the local number and velocity probability densities of a dilute two-phase suspension which has a distribution of particle sizes and a predominate direction of flow orientation such as in the cases of pipe and boundary-layer flows. It is shown that by a suitable scheme of discrimination on the amplitude as well as the residence time and frequency of the individual Laser-Doppler bursts, one can obtain the statistics on the size number density distribution and, for each size range, velocity distribution of the particulate phase together with the velocity probability distribution of the fluid phase.Results have been obtained for experiments conducted on a laminar uniform flow and a turbulent shear flow of a dilute glass particle-water suspension having a particle size distribution. Calibration needed for the scheme was accomplished by analyzing particle size and number density distribution data obtained from a Coulter particle sizing counter on a sample taken with an isokinetic probe.     

斜向自由剪切层流动特性的实验研究

本文用恒温热线风速仪对斜向自由剪切层流场进行了测量,得到了流场动量厚度、平均速度和二阶脉动速度相关统计量的变化规律;给出了不同位置处纵、展方向上速度信号随时间变化的波形、概率密度分布、自相关和互相关函数以及功率谱;表明斜向自由剪切层流场中有螺旋形涡存在的特征.     

Transport and deposition of colliding particles in turbulent channel flows

The purpose of this paper is twofold: (i) to present statistical models that describe particle–turbulence interactions as well as particle–particle collisions and (ii) to gain a better understanding of the effect of inter-particle collisions on transport, deposition, and preferential concentration of heavy particles in turbulent channel flows. The models presented are based on a kinetic equation for the probability density function of the particle velocity distribution in anisotropic turbulent flow. The model predictions compare reasonable well with numerical simulations and properly reproduce the crucial trends of computations.     

Modeling of particle motion in a turbulent flow with allowance for collisions

On the basis of the kinetic equation for the colliding-particle velocity probability density distribution in a turbulent flow, a model for calculating the dispersed phase motion is constructed for a broad range of variation of the number density and particle size.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 62–78, January–February, 1995.     

Subgrid continuum modeling of particle motion in a turbulent flow

An Eulerian continuum approach to modeling the motion of dispersed particles within the framework of the large-eddy simulation method is developed. The approach is based on a kinetic equation for the filtered probability density function for the particle velocity in a turbulent flow. Models for the subgrid turbulent stresses of the dispersed phase are presented.     

低雷诺数沟槽表面湍流/非湍流界面特性的实验研究

湍流/非湍流界面是流动中湍流和无旋流的边界,其相关研究在加深对湍流与无旋流之间的物质、动量和能量交换的理解有重要意义.本文采用时间解析的二维粒子图像测速技术,分别对零压梯度光滑、顺流向锯齿形沟槽表面平板在不同雷诺数下对湍流/非湍流界面的几何特征及动力学特性进行了实验研究.实验雷诺数为$Re_{\tau } =400\sim1000$.本文采用了湍动能准则对湍流/非湍流界面进行了识别,并分析界面高度分布、分形特征及界面附近的条件平均速度和涡量.结果表明在不同雷诺数下, 无论是光滑壁面还是沟槽壁面,界面平均高度在0.8 $\sim$ 0.9$\delta_{99} $附近. 对于沟槽壁面而言,减阻时对应的界面高度的概率密度分布与光滑壁面基本一致, 均遵循正态分布,而当阻力增大时, 界面高度分布偏离正态分布出现正的偏度. 在本实验情况下,界面分形维度、跨界面速度跳变均会随着雷诺数增大而增大. 此外,不同壁面情况下无量纲条件平均涡量在界面附近的分布相近,而界面附近无量纲速度梯度最大值近似为常数.      

EXPERIMENTAL STUDY ON PROPERTIES OF TURBULENT/NON-TURBULENT INTERFACE OVER RIBLETS SURFACES AT LOW REYNOLDS NUMBERS1)

湍流/非湍流界面是流动中湍流和无旋流的边界,其相关研究在加深对湍流与无旋流之间的物质、动量和能量交换的理解有重要意义.本文采用时间解析的二维粒子图像测速技术,分别对零压梯度光滑、顺流向锯齿形沟槽表面平板在不同雷诺数下对湍流/非湍流界面的几何特征及动力学特性进行了实验研究.实验雷诺数为$Re_{\tau } =400\sim1000$.本文采用了湍动能准则对湍流/非湍流界面进行了识别,并分析界面高度分布、分形特征及界面附近的条件平均速度和涡量.结果表明在不同雷诺数下, 无论是光滑壁面还是沟槽壁面,界面平均高度在0.8 $\sim$ 0.9$\delta_{99} $附近. 对于沟槽壁面而言,减阻时对应的界面高度的概率密度分布与光滑壁面基本一致, 均遵循正态分布,而当阻力增大时, 界面高度分布偏离正态分布出现正的偏度. 在本实验情况下,界面分形维度、跨界面速度跳变均会随着雷诺数增大而增大. 此外,不同壁面情况下无量纲条件平均涡量在界面附近的分布相近,而界面附近无量纲速度梯度最大值近似为常数.     

Statistical structure of the velocity field in cavitating flow around a 2D hydrofoil

Transformation of flow turbulence structure with cavitation occurrence, determination of the flow conditions favorable for nucleation of cavitation bubbles, influence of the statistical structure of turbulence on this process and the inverse effect of cavitation on the flow dynamics are challenging problems in modern fluid mechanics. The paper reports on the results of statistical processing of the velocity fields measured by a PIV technique in cavitating flow over a 2D symmetric hydrofoil for four flow conditions, starting from a cavitation-free regime and finishing by unsteady cloud cavitation. We analyze basic information on the statistical structure of velocity fluctuations in the form of histograms and Q-Q diagrams along with profiles of the mean velocity and turbulent kinetic energy. The research reveals that the flow turbulence pattern and distributions of turbulent fluctuations change significantly with the cavitation development. Under unsteady cloud cavitation conditions, the probability density function of the fluctuating velocity has a two-mode distribution, which indicates switching of two alternating flow conditions in a region above the hydrofoil aft part due to periodic passing of cavitation clouds. Behaviors of the mean and most probable velocities unexpectedly appear to be different with a monotonous increase of the incoming flow velocity. This finding must be caused by modification of the skewness coefficient of the fluctuating velocity.     

A Statistical Model for Predicting the Heat Transfer of Solid Particles in Turbulent Flows

The objective of this paper is twofold: (1) to present a statistical model of particle transport and heat transfer in turbulent flows and (2) to examine the performance of this model in various turbulent flows going from a simple flow to a more complicated one. This model is based on a kinetic equation for the probability density function of the particle velocity and temperature distributions in anisotropic turbulent flow. The model predictions compare reasonable well with numerical simulations and properly reproduce the crucial trends of computations performed in various turbulent flows.     

Analysis of high-order velocity moments in a strained channel flow

In the current study, model expressions for fifth-order velocity moments obtained from the truncated Gram-Charlier series expansions model for a turbulent flow field probability density function are validated using data from direct numerical simulation (DNS) of a planar turbulent flow in a strained channel. Simplicity of the model expressions, the lack of unknown coefficients, and their applicability to non-Gaussian turbulent flows make this approach attractive to use for closing turbulent models based on the Reynolds-averaged Navier-Stokes equations. The study confirms validity of the model expressions. It also shows that the imposed flow deformation improves an agreement between the model and DNS profiles for the fifth-order moments in the flow buffer zone including when the flow reverses its direction. The study reveals sensitivity of particularly odd velocity moments to the grid resolution. A new length scale is proposed as a criterion for the grid generation near the wall and in the other flow areas dominated by high mean velocity gradients when higher-order statistics have to be collected from DNS.     

Differential and algebraic models for the second moments of particle velocity and temperature fluctuations in turbulent flows

Differential and algebraic models are constructed for the dispersed-phase turbulent stresses and heat fluxes and for the mixed moments of the velocity and temperature fluctuations in the continuous and dispersed phases. The models are based on a kinetic equation for the joint probability density of the particle velocity and temperature in an anisotropic turbulent flow. The results are compared with the available direct numerical simulation (DNS) data.     

Statistical models for predicting the effect of bidisperse particle collisions on particle velocities and stresses in homogeneous anisotropic turbulent flows

The purpose of this paper is to present and compare two statistical models for predicting the effect of collisions on particle velocities and stresses in bidisperse turbulent flows. These models start from a kinetic equation for the probability density function (PDF) of the particle velocity distribution in a homogeneous anisotropic turbulent flow. The kinetic equation describes simultaneously particle–turbulence and particle–particle interactions. The paper is focused on deriving the collision terms in the governing equations of the PDF moments. One of the collision models is based on a Grad-like expansion for the PDF of the velocity distributions of two particles. The other model stems from a Grad-like expansion for the joint fluid–particle PDF. The validity of these models is explored by comparing with Lagrangian simulations of particle tracking in uniformly sheared and isotropic turbulent flows generated by LES. Notwithstanding the fact that the fluid turbulence may be isotropic, the particle velocity fluctuations are anisotropic due to the impact of gravitational settling. Comparisons of the model predictions and the numerical simulations show encouraging agreement.     

Particle collisions in a turbulent flow

An equation for the two-point probability density function of the two-particle the coordinate and velocity distribution is obtained. A closed system of equations for the first and second two-point moments of the velocity fluctuations of a pair of particles with allowance for the turbulent flow inhomogeneity is given. Boundary conditions for the equations of the particle concentration and the intensity of the relative random velocity during particle collision are obtained. A unified formula describing the interparticle collision process as a result of turbulent motion and the average relative particle velocity slip is obtained for the kernel of the coagulation equation. The effect of the average velocity slip of the particles and the carrier phase on the parameters of motion of the dispersed admixture and its coagulation is investigated on the basis of a two-point two-time velocity fluctuation autocorrelation function with two time and space scales representing the energy-bearing and small-scale motion of the fluid phase.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 104–116, March–April, 1996.     

Modeling of dynamics, heat transfer, and combustion in two-phase turbulent flows: 1. Isothermal flows

This paper presents a review of authors' collective works in the field of two-phase flow modeling done in the past few decades. The paper is aimed at the construction of mathematical models for simulation of particle-laden turbulent flows. A kinetic equation was obtained for the probability density function (PDF) of the particle velocity distribution in turbulent flows. The proposed kinetic equation describes both the interaction of particles with turbulent eddies of the carrier phase and particle-particle collisions. This PDF equation is used for the derivation of different schemes describing turbulent momentum transfer in the dispersed particle phase. The turbulent characteristics of the gaseous phase are calculated on the basis of the k - turbulence model with a modulation effect of particles on the turbulence.

The constructed models have been applied to the calculation of various two-phase gas-particle turbulent flows in jets and channels as well as particle deposition in tubes and separators. For validating the theoretical and numerical results, a wide range of comparisons with experimental data from Russian and foreign sources has been done.     

Singularity of Inertial Particle Concentration in the Viscous Sublayer of Wall-bounded Turbulent Flows

The exact kinetic equation for probability density function (PDF) of the velocity and the position of inertial particle transported by turbulent non-Gaussian fluid velocity fields in the viscous sublayer of wall-bounded turbulent flow is analyzed by the method of matched asymptotic expansions. It is shown that the particle concentration near the wall exhibits a power-law singularity giving rise to the phenomenon of particle accumulation. It is shown how the corresponding exponent depends upon the particle Stokes number. The result is in good agreement with previously published results of numerical simulations. A corresponding singularity is found for the standardized higher-order moments of particle velocity.     

紊动流场中悬浮颗粒分布的随机理论

通过分析固体颗粒在紊动流场中的随机运动,建立了二维流场中垂直于时均流动的方向上颗粒随机位移的概率密度分布函数所满足的方程。由该方程解出的分布函数在一定条件下即相当于颗粒浓度分布函数。运用这一方法研究了[1]、[2]中报道的壁面附近颗粒浓度降低的现象。     

Band-pass filtered velocity statistics in decaying turbulent box

For homogeneous isotropic turbulence study, the acquisition of band-pass filtered velocity increments (FVI) in a non-forced turbulent box is still a challenge both experimentally and numerically. Turbulence and associated physical processes, at a given instant, are permanently contaminated by a forcing process which can seldom be universal. The situation tends to be the origin of intermittency and the non-Gaussian probability density distribution for acceleration and velocity gradients. To reveal implied mechanism, grid turbulence is adapted to observe non-perturbed homogeneous isotropic turbulence. The velocity increments (VI) can be obtained following Comte-Bellot and Corrsin (GCBC) by means of two point-two time shifted velocity measurements. It is difficult to obtain decaying turbulence (DT) at large turbulent Reynolds number without pollution coming from walls. Nevertheless it is also significant to investigate DT in low Reynolds number regimes to determine non-polluted tendencies. The similarity of DT between particle image velocimetry (PIV) and hot wire anemometry measurements by GCBC are presented. Here we focus our tendency on VI and FVI probability density function (PDF) shapes in this letter. In conclusion, the tendency to Gaussian shape in inertial zone wavenumbers, demonstrates that there will be no intermittency if turbulent cascade is not perturbed.     

Relative particle velocity in a turbulent flow

On the basis of a statistical approach using a probability density function for the coordinates of two particles in a turbulent flow, the parameters of the relative particle motion are investigated. For the functions describing particle entrainment in the turbulence, rigorous results are obtained using a 3D turbulence spectrum. A method of calculating the particle relative-velocity rate with account for particle trajectory correlation is presented. The effects of particle inertia and velocity slip on the parameters of the relative particle motion are studied. Simple approximating formulas for calculating the relative particle motion in a turbulent flow are proposed. The calculation results are compared with the data of direct numerical simulation of stochastic particle trajectories in an isotropic turbulent field.     

Pipe flow of suspensions

This study shows that fully developed pipe flow of a particulate suspension is defined by four dimensionless parameters of particle-fluid interactions in addition to the Reynolds number. Effects accounted for include the Magnus effect due to fluid shear, electrostatic repulsion due to electric charges on the particles, and Brownian or turbulent diffusion. In the case of a laminar liquid-solid suspension, electrostatic effect is negligible but shear effect is prominent. Solution of the basic equations gives the density distribution of particles with a peak at the center (Einstein, Jeffery) or at other radii between the center and the pipe wall (Segré et al) depending on the magnitudes of the various flow parameters. In the case of a turbulent gas-solid suspension, the Magnus effect is significant only within the thickness of the laminar sublayer. However, charges induced on the particles by the impact of particles at the wall produce a higher density at the wall than at the center of the pipe. The velocity distribution of particles is characterized by a slip velocity at the wall and a lag in velocity in the core from the fluid phase. These results are verified by earlier measurements.