双流体颗粒-壁面碰撞模型用于旋流流动

对旋流两相流动,用考虑壁面粗糙度的双流体颗粒-壁面碰撞模型,结合二阶矩两相湍流模型和不同颗粒相边界条件,进行了数值模拟。结果表明,由于考虑了各方向雷诺应力之间的相互转化,雷诺应力从平均运动中得到能量,以及壁面对颗粒运动的衰减作用等因素,该颗粒壁面碰撞模型给出的模拟结果与实验符合较好,而广泛应用的零梯度颗粒壁面边界条件给出的模拟结果最差。     

粗糙壁面湍流边界层流动和拟序结构的大涡模拟研究

采用LES方法和虚拟粗糙元模型,,模拟了Rer=180有粗糙壁面的充分发展槽道湍流流动.从流向平均速度、脉动速度均方根、雷诺剪切应力、脉动速度相关,以及瞬时近壁准流向涡结构等多个方面,对粗糙壁面湍流边界层流动和拟序结构进行了分析.结果表明,粗糙壁面处雷诺剪切应力、脉动速度均方根以及近壁条带平均距离和准流向涡尺度都大于光滑壁面,且随着粗糙元高度的增加这种趋势更加明显.本工作为进一步研究颗粒在近壁区域的弥散规律奠定了基础.     

循环流化床颗粒相湍流结构研究

从不同层次上研究了循环流化床颗粒相的湍流结构。ＰＤＡ测量发现颗粒在壁面附近的运动表现出显著的成团特征和湍流度显著降低,小波分析法得到的直观图形象地揭示了颗粒速度脉动的自相似、分叉等具有混沌特征的微观结构,且在不同的尺度上颗粒速度脉动表现出各向异性的特征。     

含气量对激光空泡在刚性壁面空蚀的影响

采用强脉冲激光器设计液体环境下刚性壁面空蚀实验平台,改变液体中含气量,利用高速相机观察不同含气量条件下激光空泡在壁面附近的脉动过程,并对刚性壁面造成的空蚀结果进行了观测。实验研究发现,随着液体中相对空气含量的提高,激光空泡脉动的最大尺寸增大,空泡的膨胀运动变剧烈,溃灭运动速度降低,空泡的溃灭强度降低,从而影响到溃灭冲击波和壁面微射流对刚性壁面的冲击速度,减弱了壁面空蚀,而液体中含气量的提高能够降低激光空泡对刚性壁面的空蚀程度。     

横向粗糙元壁面槽道湍流的大涡模拟研究

采用大涡模拟和浸没边界法相结合对不同高度和不同间距横向粗糙元壁面槽道湍流进行了模拟,得到了光滑壁面和粗糙壁面湍流的流向平均速度分布,雷诺剪切应力,脉动速度均方根和近壁区拟序结构。结果发现横向粗糙元降低了流向平均速度,增大了流动阻力,粗糙壁面湍流的雷诺剪切应力大于光滑壁面。粗糙元降低了流向脉动速度,增强了展向和法向脉动速度。粗糙元高度越高,对湍流流动影响越大,而粗糙元间距对湍流统计特性的影响不大。粗糙壁面仍然存在着和光滑壁面类似的条带结构。     

温度场内可吸入颗粒物运动特性的实验研究

对温度场内可吸入颗粒物的运动特性进行了实验研究,使用PDA测量了场内可吸入颗粒物的速度和浓度分布等参数,研究了实验段入口速度和温度等操作参数对温度场内颗粒相运动特性的影响.结果表明,在近冷壁区边界层内,PM2.5的轴向速度和脉动速度变化很大,颗粒因热泳力及扩散作用会向冷壁面运动并产生沉积.在实验研究的基础上,提出了计算温度场内PM2.5沉积效率的经验公式.     

圆突扩两相流中气相脉动速度与速度梯度的关系

采用PDA对三种不同颗粒质量载荷的竖直向下圆管突扩流动进行了测量,发现在剪切较强的突扩发展区中气相脉动速度与速度梯度存在线性关系.采用强剪切流假设,进一步对管流中气相湍流的理论推导得到了与实验结果一致的结论.理论分析表明,颗粒浓度较小时,颗粒对气相湍流作用较弱,单相流中气相脉动速度与速度梯度的线性关系得以保持.同时推测,当增大颗粒浓度,颗粒对气相湍流的非线性作用增大,气相脉动速度与速度梯度的线性关系消失.对经典文献中的实验数据进行分析,验证了以上推论.     

高炉渣熔融液滴撞击壁面行为特性实验研究

针对高温熔渣粒化技术的开发,本文通过可视化实验研究了高炉渣熔融液滴撞击不锈钢壁面的动态行为特性。结果表明:随液滴雷诺数增大,撞击过程液滴形态演变模式依次由铺展-回缩转变为铺展-回缩破碎和铺展-破碎-凝固;增大壁面粗糙度可减弱与壁面换热,抑制液滴铺展;减小壁面粗糙度促进液滴发生破碎;液膜回缩过程普遍出现回卷现象,壁面倾角越大,液膜回卷和液滴滚动现象越显著,且液滴铺展面积越大,在壁面停留时间越长;减小液滴雷诺数、减小壁面粗糙度并采用垂直粒化仓壁面有利于壁面防黏结。     

圆锥料仓内颗粒周期性脉动特征研究

料仓卸料过程中的颗粒脉动会引起料仓振动甚至导致结构失效.为了明晰颗粒脉动特征,本文进行了椭球颗粒在不同半锥角深仓的模拟卸料实验,将料仓圆筒部划分为4个固定区域以研究区域内颗粒的运动信息,分析了料仓圆筒部颗粒系统的运动特征.研究结果表明:整个卸料过程颗粒速度始终处于波动变化中,卸料前期表现为大振幅、周期性的剧烈脉动,卸料后期平均速度的变化则是小振幅无规律的波动;剧烈脉动时段各区域的颗粒层平均受力的变化规律与颗粒速度脉动特征相似,越接近储料顶端颗粒脉动振幅越大,表现出更规律的周期性脉动,相邻颗粒层间的脉动波形相似且周期相同,剧烈脉动过程中顶面颗粒呈周期性的自由落体运动,该时段内顶层颗粒每一次的自由落体运动都会引起该范围内颗粒间接触力消失;料仓半锥角越小时剧烈脉动频率越高、振幅越大且脉动持续时间也越长,卸料速度越稳定,且颗粒速度不会出现带有上升趋势的波动.研究结果可为卸料设备的安全设计提供参考.     

三维槽道两相流颗粒运动的大涡模拟

利用基于动态粘性亚网格模型的大涡模拟方法,分别模拟了Reτ=180的三维槽道湍流内的三种颗粒,以及Reτ分别为180和395时的塑料颗粒的运动状况。其中气相控制方程的求解采用分步投影方法,泊松方程采用基于FFT的快速求解,颗粒的拉氏控制方程采用二阶龙格-库塔积分．模拟统计得到了流体和颗粒的平均速度、脉动速度等,不同的颗粒在流场中(特别是近壁区)表现出来了不同的行为,并再现了颗粒在壁面附近局部富集的现象．随Reτ提高,颗粒各方向上的速度脉动都有所增强．     

Large Eddy Simulation of Turbulent Channel Flow with 3D Roughness Using a Roughness Element Model

Large eddy simulation of turbulent channel flow with dense and small 3D roughness elements is carried out using a roughness element model Profiles of mean Reynolds stress, mean velocity and rms velocity as well as turbulent structures near the wall are obtained. The shear stress in the rough wall is larger than that in the smooth wall side and the rough wall has a larger influence on the channel flow. Profiles of mean streamwise velocity near the wall have logarithmic velocity distributions for both smooth and roughness walls, while there is a velocity decrease for the rough wall due to larger fractional drag. All the three components of rms velocities in the rough wall region are larger than that in the smooth wall region, and the roughness elements on the wall increase turbulent intensity in all directions. The s~reak spacing and average diameter of near wall quasi-s~reamwise vortices increase with the presence of roughness elements on the wall and it is shown that the rough wall induces complex and strong streamwise vortices. Results of dense and small 3D roughness elements in both turbulent statistics and structure, obtained with a relatively simple method, are found to be comparable to related experiments.     

Comparison of influence of fluidization time on electrostatic charge build-up in the bubbling vs. slugging flow regimes in gas–solid fluidized beds

Electrostatic charge generation poses significant problems in some commercial gas–solid fluidized bed reactors such as those in gas-phase polyolefin production. Understanding the contributing factors to charge generation is important in determining the charge generation mechanisms, leading to the development of methods to reduce or prevent this phenomenon. This work focused on determining the effect of fluidization time on particle charging and the amount of particle adhesion on the fluidization column wall in both the bubbling and slugging flow regimes. The charging effect was investigated for particles in three regions of the fluidized bed: elutriated fines, bulk particles inside the bed, and particles adhered to the column wall. The particles size distribution, mass and charge were measured for all three regions. Fluidization was carried out with polyethylene resins from an industrial reactor; times of 15, 30, 60, 120, 180, and 360 min were evaluated. Increased fluidization time decreased the amount of particles mass collected in the bulk region and increased those adhered to the column wall during the velocities tested in the bubbling flow regime. Whereas the quantity of particles in each region was not affected by fluidization time for the velocities examined in the slugging flow regime. Bipolar charging was observed with relatively smaller particles becoming predominately positively charged and larger particles becoming predominately negatively charged. Each region of the bed affected the magnitude of net q/m, with elutriated fines having the largest magnitude, followed by those adhered to the column wall, and finally those in the bulk of the bed. Charge saturation was attained for fluidization times greater than 60 min for particles in the bulk and along the column wall for all gas velocities. However, extended fluidization times were required with the entrained fines in bubbling flow; whereas charge saturation of fines in slugging flow occurred shortly after the onset of fluidization. Mean particle diameter for each measurement region was not impacted by the fluidization time for any of the gas velocities tested. The bed hydrodynamics was found to definitely have an impact on the particle–wall fouling where the particle layer continued to develop on the inner column wall as fluidization time increased for those velocities in the bubbling regime while comparatively less impact on particle layer growth was observed in the slugging flow regime. In addition, the bubbling flow regime resulted in particle layers formed on the column wall to be longer and thinner whereas those formed in the slugging flow regime were shorter and thicker.     

Turbulent plane wall jets over smooth and rough surfaces

Large-eddy simulations were carried out to study the effects of surface roughness on a plane wall-jet using the Lagrangian dynamic eddy-viscosity subgrid-scale model, at Re = 7500 (based on the jet bulk velocity and height). Results over both smooth and rough surfaces were validated by experimental data at the same Reynolds number. As the jet is injected into the still environment, large-scale rollers are generated in the shear layer between the high-momentum fluid of the jet and the surrounding and are convected downstream with the flow. To understand the extent to which the outer-layer structures modify the flow in the inner layer and the extent to which the effect of roughness spreads away from the wall, both instantaneous and mean flow fields were investigated. The results revealed that, for the Reynolds number and roughness height considered in this study, the effect of roughness is mostly confined to the near-wall region of the wall jet. There is no structural difference between the outer layer of the wall jet over the smooth and rough surfaces. Roughness does not affect the size of the outer-layer structures or the scaling of the profiles of Reynolds stresses in the outer layer. However, in the inner layer, roughness redistributes stresses from streamwise to wall-normal and spanwise directions toward isotropy. Contours of joint probability-density function of the streamwise and wall-normal velocity fluctuations at the bottom of the logarithmic region match those of the turbulent boundary layer at the same height; while the traces of the outer-layer structure were detected at the top of the logarithmic region, indicating that they do not affect the flow very close to the wall, but still modify a major portion of the inner layer. This modification must be taken into consideration when the inner layer of a wall jet is compared with the conventional turbulent boundary layer.     

Laser anemometry measurements in a gas-solid suspension flow

A laser anemometer has been used to measure local solid particle velocities and turbulence intensities, mean particle velocity and particle wall velocity for a fully developed gas-solid suspension flowing in a vertical pipe. The technique, which does not disturb the flow and requires no calibration, can be used successfully for internal pipe flow measurements where the solids to air mass loading ratio is below about two.     

Effect of surface roughness element on near wall turbulence with zero-pressure gradient

In the present paper,we present an investigation on the effect of roughness elements onto near-wall kinematics of a zeropressure-gradient turbulent boundary layer.An array of spanwisely-aligned cylindrical roughness elements was attached to the wall surface to regulate the near-wall low-speed streaky structures.With both qualitative visualization and quantitative measurement,we found that the regularization only occurs in the region below the height of the roughness elements.Statistical analysis on the probability distribution of the streak spanwise spacing showed that the mean spanwise streak spacing is dominated by the roughness elements;however,the latter's effect is in competition with the intrinsic streak generation mechanisms of smooth wall turbulence.Below the top of the roughness elements,local streamwise turbulent fluctuation intensity can be reduced by about 10%.We used POD analysis to depict such regularization effect in terms of near-wall structure modulation.We further found that if the spanwise spacing of roughness elements increased to be larger than the mean streak spacing in the smooth wall turbulence,there is no streak-regularization effect in the buffer region,so that the near-wall streamwise turbulent fluctuation intensity doesn't reduce.     

Surface roughness effects on the turbulence statistics in a low Reynolds number channel flow

A high-resolution particle image velocimetry was used to characterize a low Reynolds number turbulent flow in a channel. Experiments were conducted over a sand grain-coated surface of large relative roughness, and the results were compared with measurements over a smooth surface. The roughness perturbation significantly modified the outer layer. Even though the streamwise Reynolds stress shows less sensitivity in the outer layer to the boundary condition, significant enhancements were observed in the wall-normal Reynolds stress and the Reynolds shear stress. These modifications were considered as footprints of the larger-scale eddies transporting intense wall-normal motions away from the rough wall. A quadrant decomposition shows that strong and more frequent ejections are responsible for the larger values of the mean Reynolds shear stress over the rough wall. The results also indicate that spanwise vortex cores with mean vorticity of the same sign as the mean shear are the dominant smaller-scale vortical structures over the smooth and rough walls. A linear stochastic estimation-based analysis shows that the average larger-scale structure associated with these vortices is a shear layer that strongly connects the outer layer flow to the near-wall flow. A proper orthogonal decomposition of the flow suggests that the large-scale eddy is more energetic for the rough wall, and contributes more significantly to the resolved turbulent kinetic energy and the Reynolds shear stress than the smooth wall.     

Recycling inflow method for simulations of spatially evolving turbulent boundary layers over rough surfaces

The technique by Lund et al. to generate turbulent inflow for simulations of developing boundary layers over smooth flat plates is extended to the case of surfaces with roughness elements. In the Lund et al. method, turbulent velocities on a sampling plane are rescaled and recycled back to the inlet as inflow boundary condition. To rescale mean and fluctuating velocities, appropriate length scales need be identified and for smooth surfaces, the viscous scale l_ν = ν/u_τ (where ν is the kinematic viscosity and u_τ is the friction velocity) is employed for the inner layer. Different from smooth surfaces, in rough wall boundary layers the length scale of the inner layer, i.e. the roughness sub-layer scale l_d, must be determined by the geometric details of the surface roughness elements and the flow around them. In the proposed approach, it is determined by diagnosing dispersive stresses that quantify the spatial inhomogeneity caused by the roughness elements in the flow. The scale l_d is used for rescaling in the inner layer, and the boundary layer thickness δ is used in the outer region. Both parts are then combined for recycling using a blending function. Unlike the blending function proposed by Lund et al. which transitions from the inner layer to the outer layer at approximately 0.2δ, here the location of blending is shifted upwards to enable simulations of very rough surfaces in which the roughness length may exceed the height of 0.2δ assumed in the traditional method. The extended rescaling–recycling method is tested in large eddy simulation of flow over surfaces with various types of roughness element shapes.     

Lattice Boltzmann Simulations of Particle-Particle Interaction in Steady Poiseuille Flow

The moving behaviour of two- and three-particles in a pressure-driven flow is studied by the lattice Boltzmann simulation in two dimensions. The time-dependent values, including particles＇ radial positions, translational velocities, angular velocities, and the x-directional distance between the particles are analysed extensively. The effect of flow Reynolds number on particle motion is also investigated numerically. The simulation results show that the leading particle equilibrium position is closer to the channel centre while the trailing particle equilibrium position is closer to the channel wall. If Reynolds number Re is less than 85.30, the larger flow Reynolds number results in the smaller x-directional equilibrium distance, otherwise the x-directional distance increases almost linearly with the increase of time and the particles separate finally. The simulation results are helpful to understand the particle-particle interaction in suspensions with swarms of particles.     

用格子Boltzmann模型模拟可压缩完全气体流动

采用一种新的格子Boltzmann模型模拟超音速流动。在这种模型中,粒子的速度不受限制,可以取得很广。而平衡分布函数的支集却相对集中,使模型得以简化。粒子速度的这种自适应特性允许流体以较高的马赫数流动。通过引入粒子的势能使得该模型适用于具有任意比热比的完全气体。利用Chapman-Enskog方法,从BGK型Boltzmann方程推导出Navier-Stokes方程。在六边形网格上模拟了马赫数为3的前台阶绕流,得到了合理的结果。