Simultaneous PIV and PTV measurements of wind and sand particle velocities

Wind-blown sand is a typical example of two-phase particle-laden flows. Owing to lack of simultaneous measured data of the wind and wind-blown sand, interactions between them have not yet been fully understood. In this study, natural sand of 100–125 μm taken from Taklimakan Desert was tested at the freestream wind speed of 8.3 m/s in an atmospheric boundary layer wind tunnel. The captured flow images containing both saltating sand and small wind tracer particles, were separated by using a digital phase mask technique. The 2-D PIV (particle imaging velocimetry) and PTV (particle tracking velocimetry) techniques were employed to extract simultaneously the wind velocity field and the velocity field of dispersed sand particles, respectively. Comparison of the mean streamwise wind velocity profile and the turbulence statistics with and without sand transportation reveal a significant influence of sand movement on the wind field, especially in the dense saltating sand layer (y/δ < 0.1). The ensemble-averaged streamwise velocity profile of sand particles was also evaluated to investigate the velocity lag between the sand and the wind. This study would be helpful in improving the understanding of interactions between the wind and the wind-blown sand.     

Dispersion de particules solides en mouvement de saltation dans un écoulement turbulentDispersion of solid saltating particles in a turbulent flow

The solid particle dispersion in saltating motion is studied in an homogeneous turbulence and in a turbulent boundary layer. The fluid velocity along the particle trajectory is estimated using a continuous stochastic differential equation in which the correlation integral time takes into account gravity and inertia effects. As far as the boundary layer is concerned, the aerodynamic entrainment of particles and the rebound are modelised as random variables with Gaussian probability density functions. Compared with experimental results, the numerical results show good agreement for dispersion, although velocity fluctuations are slightly under evaluated. To cite this article: C. Aguirre et al., C. R. Mecanique 332 (2004).     

Effect of thermal diffusion and electrostatic force on evolution of wind-blown sand flow

A theoretical model is suggested to mathematically describe the effect of thermal diffusion from a sand-bed on evolution of a wind-blown sand flow.An upward wind field is engendered by the thermal diffusion and the coupling interaction among the horizontal and upward wind flow,saltating grains,and a kind of electrostatic force exerted on the grains are considered in this theoretical model.The numerical results show that the effect of the thermal diffusion on the evolution process of wind-blown grain flow is quite obvious and very similar to the effect of the electrostatic force on the evolution.Not only the time for the entire system to reach a steady state(called the duration time),the transport rate of grains,the mass-flux profiles and the trajectory of saltating grains are affected by the thermal diffusion and the electrostatic force exerted on saltating grains, but also the wind profiles and the temperature profiles at the steady state are affected by the wind-blown sand flow.     

风沙两相流跃移层中沙粒相的速度分布

从单个跃移沙粒在气流中的运动方程出发导出了风沙两相流中沙粒相速度分布函数的Boltzmann方程，对风沙流研究中几种不同的分布函数及其相应的统计平均值等基本概念给出了严密的数学定义，指出了不同分布函数之间的区别和联系，在略去铅垂方向空气阻力的情况下，给出了沙粒相速度分布函数沿铅垂方向的边缘分布，作为风沙流中跃移理论的主要基础之一。利用结果对前人在风沙流研究中发现的某些重要规律和现象进行了解释。     

含尘离心风机中边界层对叶轮侵蚀的影响

本文采用二维可压缩紊流边界层计算模型,计算了离心风机内考虑边界层影响后的粒子轨迹和粒子与叶片表面的碰撞角、碰撞速度、磨损率沿叶片(相对弧长)的变化规律,计算结果表明:对于小直径的粒子,无粘流动假设是不合适的,计及边界层的影响,对于准确预测风机内的粒子轨迹和磨损是必不可少的。     

Experimental study and numerical simulation of the characteristics of a percussive gas–solid separator

The presence of solid particles in the flow of hypersonic wind tunnels damages the appearance of the experiment models in the wind tunnel and influences the accuracy of experimental results. The design of a highly efficient gas–solid separator was therefore undertaken. Particle trajectory imaging methods were used to measure trajectories under different conditions. The flow field and particle movement characteristics for different head angles (HAs) and separation tooth angles (STAs), inlet velocities, and the exhaust gas outlet pressures in the separator, were calculated using simulations based on the discrete phase model. The particle separation efficiency, pressure loss, and flow loss resulting from different structural parameters were also studied. In line with experimental observations, the characteristic angle of particle movements in the separator and the separation efficiency of the separator were found to increase with decreasing HA and with increasing STA. Separation efficiency improves with increasing inlet velocity and with increasing negative pressure of the exhaust gas outlet; however, the corresponding pressure loss and the flow rate of the waste gas also increased.     

Near-wake flow structure downwind of a wind turbine in a turbulent boundary layer

Wind turbines operate in the surface layer of the atmospheric boundary layer, where they are subjected to strong wind shear and relatively high turbulence levels. These incoming boundary layer flow characteristics are expected to affect the structure of wind turbine wakes. The near-wake region is characterized by a complex coupled vortex system (including helicoidal tip vortices), unsteadiness and strong turbulence heterogeneity. Limited information about the spatial distribution of turbulence in the near wake, the vortex behavior and their influence on the downwind development of the far wake hinders our capability to predict wind turbine power production and fatigue loads in wind farms. This calls for a better understanding of the spatial distribution of the 3D flow and coherent turbulence structures in the near wake. Systematic wind-tunnel experiments were designed and carried out to characterize the structure of the near-wake flow downwind of a model wind turbine placed in a neutral boundary layer flow. A horizontal-axis, three-blade wind turbine model, with a rotor diameter of 13 cm and the hub height at 10.5 cm, occupied the lowest one-third of the boundary layer. High-resolution particle image velocimetry (PIV) was used to measure velocities in multiple vertical stream-wise planes (x–z) and vertical span-wise planes (y–z). In particular, we identified localized regions of strong vorticity and swirling strength, which are the signature of helicoidal tip vortices. These vortices are most pronounced at the top-tip level and persist up to a distance of two to three rotor diameters downwind. The measurements also reveal strong flow rotation and a highly non-axisymmetric distribution of the mean flow and turbulence structure in the near wake. The results provide new insight into the physical mechanisms that govern the development of the near wake of a wind turbine immersed in a neutral boundary layer. They also serve as important data for the development and validation of numerical models.     

压力梯度对边界层两相流动稳定性的影响

分析了有压力梯度的边界层两相流动稳定性,推导出类似于Saffman理论的修正的稳定性方程,数值计算采用高精度的谱方法。结果说明,压力梯度对边界层两相流动稳定性有显著的影响,顺压梯度增强流动稳定性,而逆压梯度则促进流动失稳。在不同的压力梯度和浓度下,Stokes数对流动稳定性的影响是一致的,存在一个临界Stokes数,小Stokes数促进流动失稳,而大Stokes数则提高临界雷诺数,抑制流动失稳的最佳Stokes数为10的量级。     

Quantitative visualization of compressible turbulent shear flows using condensate-enhanced Rayleigh scattering

This paper describes several flow visualization experiments carried out in Mach 3 and Mach 8 turbulent shear flows. The experimental technique was based on laser scattering from particles of H₂O or CO₂ condensate that form in the wind tunnel nozzle expansion process. The condensate particles vaporize extremely rapidly on entering the relatively hot fluid within a turbulent structure, so that a sharp vaporization interface marks the outer edge of the rotational shear layer fluid. Calculations indicate that the observed thin interface corresponds to a particle size of 10 nm or less, which is consistent with optical measurements, and that particles of this size track the fluid motions well. Further, calculations and experiments show that the freestream concentration of condensate required for flow visualization has only a small effect on the wind tunnel pressure distribution. Statistics based on the image data were compared to corresponding results from probe measurements and agreement was obtained in statistical measures of speed, scale, and orientation of the large-scale structures in the shear layer turbulence. The condensate-enhanced Rayleigh scattering technique is judged to be a useful tool for quantitative studies of shear layer structure, particularly for identifying the instantaneous boundary layer edge and for extracting comparative information on the large-scale structures represented there.     

Dynamic wind loads and wake characteristics of a wind turbine model in an atmospheric boundary layer wind

An experimental study was conducted to characterize the dynamic wind loads and evolution of the unsteady vortex and turbulent flow structures in the near wake of a horizontal axis wind turbine model placed in an atmospheric boundary layer wind tunnel. In addition to measuring dynamic wind loads (i.e., aerodynamic forces and bending moments) acting on the wind turbine model by using a high-sensitive force-moment sensor unit, a high-resolution digital particle image velocimetry (PIV) system was used to achieve flow field measurements to quantify the characteristics of the turbulent vortex flow in the near wake of the wind turbine model. Besides conducting “free-run” PIV measurements to determine the ensemble-averaged statistics of the flow quantities such as mean velocity, Reynolds stress, and turbulence kinetic energy (TKE) distributions in the wake flow, “phase-locked” PIV measurements were also performed to elucidate further details about evolution of the unsteady vortex structures in the wake flow in relation to the position of the rotating turbine blades. The effects of the tip-speed-ratio of the wind turbine model on the dynamic wind loads and wake flow characteristics were quantified in the terms of the variations of the aerodynamic thrust and bending moment coefficients of the wind turbine model, the evolution of the helical tip vortices and the unsteady vortices shedding from the blade roots and turbine nacelle, the deceleration of the incoming airflows after passing the rotation disk of the turbine blades, the TKE and Reynolds stress distributions in the near wake of the wind turbine model. The detailed flow field measurements were correlated with the dynamic wind load measurements to elucidate underlying physics in order to gain further insight into the characteristics of the dynamic wind loads and turbulent vortex flows in the wakes of wind turbines for the optimal design of the wind turbines operating in atmospheric boundary layer winds.     

Analysis of sand particles’ lift-off and incident velocities in wind-blown sand flux

In the research of windblown sand movement,the lift-off and incident velocities of saltating sand particles play a significant role in bridging the spatial and temporal scales from single sand particle’s motion to windblown sand flux.In this paper,we achieved wind tunnel measurements of the movement of sand particles near sand bed through improving the wind tunnel experimental scheme of paticle image velocimetry(PIV) and data processing method.And then the influence of observation height on the probability distributions of lift-off and incident velocities of sand particles was analyzed.The results demonstrate that the observation height has no obvious influence on the distribution pattern of the lift-off and incident velocities of sand particles,i.e.,the probability distribution of horizontal and vertical velocities of lift-off and incident sand particles follow a Gaussian distribution and a negative exponential distribution,respectively.However,it influences the center of the Gaussian distribution,the decay constant and the amplitude of the negative exponential distribution.     

HYDRODYNAMIC RESISTANCE EFFECT OF FLUID LAYER BETWEEN TWO IMMERSED APPROACHING PARTICLES

1. Introduction The mechanisms of impact and rebound of solid parti- cles in particulate flow systems are of interest over a wide range of application areas such as fluidized beds, pneu- matic transport, filtration processes, erosion and pollution control of suspended particles. In many cases, the colli- sions of particles against themselves and against walls may affect the properties of the mixture. Efforts have been made to describe the fundamental mechanics of particle collisions. The conta…     

可侵蚀地表上方含尘大气运动的数值模型

在稀相气固两相流的双流体模型框架下，导出含尘大气运动的基本方程，其中的源项考虑相间双向耦合作用；通过对基本方程进行无量纲化，求得控制气相和弥散相动力学行为的相似准则，包括弗洛德数、罗斯贝数、颗粒雷诺数、颗粒惯性参数、颗粒剪切横移参数、颗粒旋转横移参数和沙尘质量载荷率等. 作为一个模型问题，研究可侵蚀地表上方充分发展湍流大气边界层流动对土壤颗粒的气动卷扬作用，其中计及气动阻力、萨弗曼升力和重力等因素的影响. 为了克服由于沙尘跃移运动轨迹交叉导致流动参数多值性所造成的困难，引进拉格朗日坐标下弥散相连续方程. 对两种不同风速和4种不同粒径条件下的沙尘运动进行数值模拟，给出相应情况中沙尘运动特性和密度分布剖面并讨论风速和粒径等参数的影响，还细致分析了含尘大气边界层中沙尘与气流之间的能量传递过程. 所得的方程、准则和方法可以为土壤风蚀和沙尘暴等自然灾害的预报提供理论基础.     

垂直旋转流场中单颗粒运动状态判别及受力分析

为研究单颗粒在旋转流场中的运动状态及受力情况,以毫米级球形颗粒为例,利用旋转流场颗粒运动装置,通过使用摄像机记录颗粒在流场中的运动轨迹以获取其运动参数,分析了不同转速和颗粒直径条件下颗粒的运动轨迹,拟合得到了颗粒运动状态判别公式以及颗粒运动轨迹公式,分析了颗粒在旋转流场中的受力情况。结果表明,颗粒在旋转流场平衡状态下运动状态主要分为两类,一类是未离开壁面保持静止,另一类是离开壁面保持稳定周向运动;颗粒进行周向运动的轨迹为椭圆形,并且圆心随着转速的增大靠近旋转中心,而随着粒径的增大靠近壁面;颗粒在旋转流场的运动过程中主要受到离心力和旋转科式力作用。     

Determination of real-time predictors of the wind turbine wake meandering

The present work proposes an experimental methodology to characterize the unsteady properties of a wind turbine wake, called meandering, and particularly its ability to follow the large-scale motions induced by large turbulent eddies contained in the approach flow. The measurements were made in an atmospheric boundary layer wind tunnel. The wind turbine model is based on the actuator disc concept. One part of the work has been dedicated to the development of a methodology for horizontal wake tracking by mean of a transverse hot wire rake, whose dynamic response is adequate for spectral analysis. Spectral coherence analysis shows that the horizontal position of the wake correlates well with the upstream transverse velocity, especially for wavelength larger than three times the diameter of the disc but less so for smaller scales. Therefore, it is concluded that the wake is actually a rather passive tracer of the large surrounding turbulent structures. The influence of the rotor size and downstream distance on the wake meandering is studied. The fluctuations of the lateral force and the yawing torque affecting the wind turbine model are also measured and correlated with the wake meandering. Two approach flow configurations are then tested: an undisturbed incoming flow (modelled atmospheric boundary layer) and a disturbed incoming flow, with a wind turbine model located upstream. Results showed that the meandering process is amplified by the presence of the upstream wake. It is shown that the coherence between the lateral force fluctuations and the horizontal wake position is significant up to length scales larger than twice the wind turbine model diameter. This leads to the conclusion that the lateral force is a better candidate than the upstream transverse velocity to predict in real time the meandering process, for either undisturbed (wake free) or disturbed incoming atmospheric flows.     

Suspended load and bed-load transport of particle-laden gravity currents: the role of particle–bed interaction

The development of particle-enriched regions (bed-load) at the base of particle-laden gravity currents has been widely observed, yet the controls and relative partitioning of material into the bed-load is poorly understood. We examine particle-laden gravity currents whose initial mixture (particle and fluid) density is greater than the ambient fluid, but whose interstitial fluid density is less than the ambient fluid (such as occurs in pyroclastic flows produced during volcanic eruptions or when sediment-enriched river discharge enters the ocean, generating hyperpycnal turbidity currents). A multifluid numerical approach is employed to assess suspended load and bed-load transport in particle-laden gravity currents under varying boundary conditions. Particle-laden flows that traverse denser fluid (such as pyroclastic flows crossing water) have leaky boundaries that provide the conceptual framework to study suspended load in isolation from bed-load transport. We develop leaky and saltation boundary conditions to study the influence of flow substrate on the development of bed-load. Flows with saltating boundaries develop particle–enriched basal layers (bed-load) where momentum transfer is primarily a result of particle–particle collisions. The grain size distribution is more homogeneous in the bed-load and the saltation boundaries increase the run-out distance and residence time of particles in the flow by as much as 25% over leaky boundary conditions. Transport over a leaky substrate removes particles that reach the bottom boundary and only the suspended load remains. Particle transport to the boundary is proportional to the settling velocity of particles, and flow dilution results in shear and buoyancy instabilities at the upper interface of these flows. These instabilities entrain ambient fluid, and the continued dilution ultimately results in these currents becoming less dense than the ambient fluid. A unifying concept is energy dissipation due to particle–boundary interaction: leaky boundaries dissipate energy more efficiently at the boundary than their saltating counterparts and have smaller run-out distance.

---

     

A direct‐forcing pressure‐based lattice Boltzmann method for solving fluid–particle interaction problems

A direct‐forcing immersed boundary‐lattice Boltzmann method (IB–LBM) is developed to simulate fluid–particle interaction problems. This method uses the pressure‐based LBM to solve the incompressible flow field and the immersed boundary method to handle the fluid–particle interactions. The pressure‐based LBM uses the pressure distribution functions instead of the density distribution functions as the independent dynamic variables. The main idea is to explicitly eliminate the compressible effect due to the density fluctuation. In the IB method, a direct‐forcing method is introduced to capture the particle motion. It directly computes an IB force density at each lattice grid from the differences between the pressure distribution functions obtained by the LBM and the equilibrium pressure distribution functions computed from the particle velocity. By applying this direct‐forcing method, the IB–LBM becomes a purely LBM version. Also, by applying the Gauss theorem, the formulas for computing the force and the torque acting on the particle from the flows are derived from the volume integrals over the particle volume instead of from the surface integrals over the particle surface. The order of accuracy of the IB–LBM is demonstrated on the errors of velocity field, wall stress, and gradients of velocity and pressure. As a demonstration of the efficiency and capabilities of the new method, sedimentation of a large number of spherical particles in an enclosure is simulated. Copyright © 2010 John Wiley & Sons, Ltd.     

湍流边界层中固体颗粒运动的数值模拟

根据Ｌａｇｒａｎｇｅ颗粒运动微分方程及不可压缩湍流边界层中流体的壁面速度分布规律，数值求解了颗粒在湍流边界层中的运动，考虑了Ｓａｆｆｍａｎ升为对颗粒运动的影响，壁面对运动阻力的影响，给出了固体颗粒沉积边壁，在边界层外缘上所需的最小速度和最小入射角，计算结果还表明边界层对固体颗粒撞击边壁的速度和入射角有较大影响，从数值结果可可以发现一个重要现象。     

CFD-DEM simulation of spouting of corn-shaped particles

Three dimensionally coupled computational fluid dynamics (CFD) and discrete element method (DEM) were used to investigate the flow of corn-shaped particles in a cylindrical spouted bed with a conical base. The particle motion was modeled by the DEM, and the gas motion by the k-? two-equation turbulent model. A two-way coupling numerical iterative scheme was used to incorporate the effects of gas–particle interactions in terms of momentum exchange. The corn-shaped particles were constructed by a multi-sphere method. Drag force, contact force, Saffman lift force, Magnus lift force, and gravitational force acting on each individual particle were considered in establishing the mathematical modeling. Calculations were carried out in a cylindrical spouted bed with an inside diameter of 200 mm, a height of 700 mm, and a conical base of 60°. Comparison of simulations with experiments showed the availability of the multi-sphere method in simulating spouting action with corn-shaped particles, but it depended strongly on the number and the arrangement of the spherical elements. Gas–solid flow patterns, pressure drop, particle velocity and particle concentration at various spouting gas velocity were discussed. The results showed that particle velocity reaches a maximum at the axis and then decreases gradually along the radial direction in the whole bed. Particle concentration increases along the radial direction in the spout region but decreases in the fountain region, while it is nearly constant in the annulus region. Increasing spouting gas velocity leads to larger pressure drop, remarkably increased speed of particle moving upward or downward, but decreased particle concentration.     

双柱状颗粒在线性剪切流场中的动力学分析

为研究柱状颗粒在线性剪切流场中的运动状态和受力情况,本文以颗粒长径比为2,颗粒之间的初始距离ΔS_Py=4D为例,基于直接力浸入边界法数值模拟了双柱状颗粒在三维线性剪切流场中的运动过程。根据模拟结果分析了柱状颗粒周围流场参数分布,在考虑壁面对颗粒的影响和颗粒之间相互影响的条件下,研究了颗粒的受力和运动的变化,探索了流体曳力导致柱状颗粒迁移和转动的规律。研究结果表明,双柱状颗粒在线性剪切流场中易向速度大的流体区域运动;前后两颗粒运动状态和轨迹不同,颗粒之间距离较近时,曳力会产生较大的波动;只有当颗粒在壁面附近时,滞后颗粒才能追上领先颗粒,两颗粒发生牵引、翻滚和分离过程。