液固两相流冲洗油管道的冲蚀磨损特性数值模拟及分析

以液固两相流冲洗油管道为研究对象,采用Realizableκ-ε湍流模型、随机轨道模型,结合液固两相流冲蚀磨损试验,建立修正的冲蚀磨损数理模型,数值预测典型工况下冲洗油管道内速度、压力、冲蚀磨损率等流动参数分布情况,分析了冲蚀磨损的形成机制.研究结果表明:受曲率半径影响,冲洗油管道冲蚀磨损速率随曲率半径的增加而减小;由于颗粒惯性及管内二次流影响,弯头中间区域外侧壁面和出口直管段内侧面磨损较为严重,三通管件的最大冲蚀磨损率位于孔口处,数值预测结果与失效解剖案例吻合.本文建立的冲蚀磨损定量预测方法,适用于压力管道的风险评定及寿命预测.     

离心泵叶轮内密相液固两相湍流的数值模拟

利用作者建立的描述密相液固两相湍流的 KET模型和推导的基本控制方程组 ,在处理壁面边界条件时考虑了颗粒和叶片的相互碰撞作用 ,对离心泵叶轮内密相液固两相流动进行了数值模拟 ,得到了泵叶轮内两相流动的一些规律 ,为液固两相流泵的设计提供了一定的理论依据。     

90°弯管内气固两相湍流流动的数值研究

利用两相湍流KET模型对90°弯管内气固两相湍流流动进行了数值模拟,得到了弯管内两相流动的一些规律,并提出用颗粒动理学压力来定性表征弯管内磨损严重部位,为管道抗磨损设计提供了一定的理论依据。     

稀相气固两相流90°圆截面弯管流动特性

针对工业中常见的圆截面90°弯管中的稀相气固两相流的流动特征难于观测,以及由此造成的对其流动特性缺乏认识的问题,以天津大学的气固装置具备的实验条件为基础,以压力测量值作为模型有效性的评价依据,建立了适用于圆截面管道的计算流体力学(Computational Fluid Dynamics,CFD)模型.针对FLUENT中通过定义反射系数设置颗粒-壁面碰撞关系不能反映实际过程的问题,建立了考虑颗粒与壁面接触过程中的不同运动状态的颗粒-壁面碰撞模型,并引入计算.以此为基础,对不同固气质量比、流量和管道曲率情况下的流动特性进行了研究.发现:两相流体在流经90°弯管后,速度分布的恢复长度与固气质量比、管道曲率半径等有密切关系,曲率半径为2的管道整体压损最小,过小或过大均会导致较大压损,与Mason弯管磨损试验的结论规律相近.     

固-液两相流黑水管道冲蚀磨损的数值模拟研究

煤气化黑水处理系统管道由于其流体介质高含固体颗粒和腐蚀性介质,且工作在高温、高压差环境中,极易受到冲蚀磨损和腐蚀的耦合作用而失效,影响其服役寿命. 采用计算流体力学(CFD)方法数值模拟研究了煤气化黑水处理系统固-液两相流管道的冲蚀磨损行为和机理,以及流体介质速度和固体颗粒粒径对管道冲蚀磨损的影响规律,并分析了盲通管和涡室结构对弯管冲蚀磨损行为的优化改善效果. 研究结果显示,煤气化黑水处理系统管线的冲蚀高危区主要分布在弯管外拱和变径管等结构突变区域;管道冲蚀磨损行为与其内部流体的运动和颗粒冲击特性有关;管道的冲蚀率均随着流体速度的增加而加剧,而粒径对弯管和变径管冲蚀率的影响并非单调关系,这与颗粒受力作用有关;弯管优化分析显示,涡室结构可以降低弯管的最大冲蚀率,减缓弯管的冲蚀磨损.      

大弯径比挠性弯管内气固两相流冲蚀特性研究

以90°大弯径比挠性弯管为研究对象,基于离散相模型和General冲蚀模型,利用Fluent数值仿真软件,开展不同工况下气相进口速度、夹带固体颗粒质量流量、弯管肋条及弯径比等影响因素对弯管内气固两相流冲蚀特性的数值模拟研究。研究结果表明：弯管纵切面气固两相流场速度分布存在分离现象;弯管的最大冲蚀速率随气相进口速度的增长呈指数关系,随颗粒质量流量的增长呈线性关系,随肋条数目的增长先减小后略微增大最后趋于稳定;大弯径比弯管内部流场二次流强度小且肋条的存在使弯管的耐冲蚀性能得到一定的提升;弯管冲蚀最严重区域发生在弯管沿流动方向偏转26°左右夹角的区域且呈现出V型冲蚀形貌,这种现象是由于弯管处独特的颗粒运动轨迹决定的;在弯管沿流动方向60°～80°夹角范围内出现沟槽状冲蚀形貌,它是由于颗粒在弯管内的二次及多次碰撞。     

一种减轻固粒对壁面冲蚀磨损的新方法

针对带固粒的近壁流场，提出了在壁面上开凿纵向浅沟槽以减轻固粒对壁面冲蚀磨损的方法．采用改进的κ-ε湍流模式和气-固两相双向耦合模型，计算了流场中颗粒的速度、轨迹以及固粒对壁面的冲蚀磨损，并进行了相应的试验对比．结果表明：在壁面上开凿纵向沟槽能减轻固粒对壁面的磨损；在一定的沟槽高度下，当沟槽宽度同沟槽问距离相等时，减轻磨损的效果最明显；在一定颗粒尺寸范围内，大颗粒导致磨损质量损失增大，当颗粒尺寸超过某一临界值后，颗粒尺寸对磨损的影响较小．进口处颗粒运动方向对壁面磨损具有直接影响．     

稠密气固两相流各向异性颗粒相矩方法

基于气体分子动力学和颗粒动理学方法,考虑颗粒速度脉动各向异性,建立颗粒相二阶矩模型.应用初等输运理论,对三阶关联项进行模化和封闭.考虑颗粒与壁面之间的能量传递和交换,建立颗粒相边界条件模型.采用Koch等计算方法模拟气固脉动速度关联矩.考虑气体-颗粒间相互作用,建立稠密气体-颗粒流动模型.数值模拟提升管内气固两相流动特性,模拟结果表明提升管内颗粒相湍流脉动具有明显的各向异性.预测颗粒速度、浓度和颗粒脉动速度二阶矩与Tartan等实测结果相吻合.模拟结果表明轴向颗粒速度脉动强度约为平均颗粒相脉动强度的1.5倍,轴向颗粒脉动能大约是径向颗粒脉动能3.0倍.     

非球形颗粒两相流的数值模拟研究进展

非球形颗粒两相流是多相流的重要研究方向之一, 常见于自然界及工业生产过程中. 不同于球形颗粒, 由于非球形颗粒形状的各向异性, 除了颗粒平动行为, 还需要考虑颗粒的转动与取向行为, 颗粒的取向与转动行为会影响颗粒所受的力和力矩. 为了准确模拟非球形颗粒的运动行为, 目前非球形颗粒两相流的数值模拟研究主要基于欧拉?拉格朗日的求解框架展开, 常见的非球形颗粒两相流数值模拟方法主要包括点颗粒法与全分辨颗粒法. 本文将对这两类方法进行介绍, 同时会全面介绍非球形颗粒两相流研究的基础理论模型, 并系统总结非球形颗粒在简单基本流和复杂湍流中的研究进展, 包括对于非球形颗粒在湍流中的取向与转动行为机理, 以及颗粒对湍流减阻调制作用的研究. 最后, 本文提出了非球形颗粒两相流研究存在的问题及未来研究方向.      

Y型分支管内气固两相流动的数值模拟

本文采用Euler-Lagrange两相流研究方法,固相采用离散相(DPM)模型,对三种不同夹角的Y型分支管内气固两相流动进行了数值模拟.气相湍流采用Realizable k-ε模型,固相的湍流耗散采用随机轨道模型.模拟结果较好地预估了颗粒在分叉处的流动形态、颗粒在分支管内的运动轨迹,以及重新实现颗粒相流场均匀分布所需的距离.同时发现,随机轨道模型较适于评估分支与主管夹角较大情况下,固体颗粒在分叉处的运动.将分支管内固体颗粒质量分配的数值模拟结果与试验结果比较,发现两者较吻合,相对误差较小.     

固体颗粒对沟槽湍流边界层影响的实验研究

本文采用粒子图像测速技术(particles image velocimetry, PIV)研究固体颗粒对放置在平板湍流边界层中的平壁和沟槽壁面减阻效果的影响. 实验对清水和加入粒径为155 μm聚苯乙烯颗粒的流法向二维速度场信息进行采集, 对不同工况下的平均速度剖面、雷诺应力和湍流度等统计量进行对比, 分析流体在边界层中的行为. 运用空间局部平均结构函数提取了不同工况湍流边界层喷射?扫掠行为的空间拓扑结构并进行比较. 结果发现, 在不同的壁面条件下, 粒子加入后的对数律区中无量纲速度均略大于清水组, 雷诺切应力有所降低, 湍流度有所减弱. 对于不同流场速度下的沟槽而言, 颗粒的加入均降低了壁面附近的阻力, 而颗粒单独作用于光滑壁面的减阻效果并不明显. 加入粒子后的相干结构数目有所增加, 法向脉动速度下降. 沟槽壁面附近的相干结构数目有所增加, 法向脉动速度在自由来流速度较大时有所上升, 在速度较小时有所下降. 这表明不同减阻状况下的沟槽均能将大涡破碎成更多的涡, 并且粒子的加入强化了这种破碎作用.      

The simulation of turbulent particle‐laden channel flow by the Lattice Boltzmann method

We perform direct numerical simulation of three‐dimensional turbulent flows in a rectangular channel, with a lattice Boltzmann method, efficiently implemented on heavily parallel general purpose graphical processor units. After validating the method for a single fluid, for standard boundary layer problems, we study changes in mean and turbulent properties of particle‐laden flows, as a function of particle size and concentration. The problem of physical interest for this application is the effect of water droplets on the turbulent properties of a high‐speed air flow, near a solid surface. To do so, we use a Lagrangian tracking approach for a large number of rigid spherical point particles, whose motion is forced by drag forces caused by the fluid flow; particle effects on the latter are in turn represented by distributed volume forces in the lattice Boltzmann method. Results suggest that, while mean flow properties are only slightly affected, unless a very large concentration of particles is used, the turbulent vortices present near the boundary are significantly damped and broken down by the turbulent motion of the heavy particles, and both turbulent Reynolds stresses and the production of turbulent kinetic energy are decreased because of the particle effects. We also find that the streamwise component of turbulent velocity fluctuations is increased, while the spanwise and wall‐normal components are decreased, as compared with the single fluid channel case. Additionally, the streamwise velocity of the carrier (air) phase is slightly reduced in the logarithmic boundary layer near the solid walls. Copyright © 2015 John Wiley & Sons, Ltd.     

含砂流对弯管段冲蚀的机理分析

管道的弯管段、闸阀等处常常面临含砂流动的冲蚀磨损而导致管壁减薄或泄漏。为研究弯管段冲蚀磨损的主要机制,进行了加入5%石英砂的固液两相管流循环试验。试验按砂粒粒径分组,主要观测布置在弯管段的一组碳钢贴片在2~3小时后的壁面冲蚀特征。用失重法计算各贴片的冲蚀速率,基于扫描电镜(SEM)研究了贴片表面的损伤形貌。冲蚀后的贴片表面既有犁沟、挤压凹坑及划痕,亦有圆形蚀坑。能谱分析(EDS)显示有一定量的铁氧化物,且不同位置处的贴片失重率有正有负,由此可见导致贴片表面破坏的机制既包括切削、挤压和二次冲击等含沙流流动的机械作用,也伴随着腐蚀作用。这两种机制在弯管各区有不同的体现,且受砂粒粒径的影响。用数值模拟的两相流冲蚀弯管段流场压力和颗粒分布特征做了验证。结果表明弯管段的壁面损伤机制在外拱壁面以颗粒冲击和切削为主,在内拱璧面则以氧化腐蚀为主。     

Using digital holographic microscopy for simultaneous measurements of 3D near wall velocity and wall shear stress in a turbulent boundary layer

A digital holographic microscope is used to simultaneously measure the instantaneous 3D flow structure in the inner part of a turbulent boundary layer over a smooth wall, and the spatial distribution of wall shear stresses. The measurements are performed in a fully developed turbulent channel flow within square duct, at a moderately high Reynolds number. The sample volume size is 90 × 145 × 90 wall units, and the spatial resolution of the measurements is 3–8 wall units in streamwise and spanwise directions and one wall unit in the wall-normal direction. The paper describes the data acquisition and analysis procedures, including the particle tracking method and associated method for matching of particle pairs. The uncertainty in velocity is estimated to be better than 1 mm/s, less than 0.05% of the free stream velocity, by comparing the statistics of the normalized velocity divergence to divergence obtained by randomly adding an error of 1 mm/s to the data. Spatial distributions of wall shear stresses are approximated with the least square fit of velocity measurements in the viscous sublayer. Mean flow profiles and statistics of velocity fluctuations agree very well with expectations. Joint probability density distributions of instantaneous spanwise and streamwise wall shear stresses demonstrate the significance of near-wall coherent structures. The near wall 3D flow structures are classified into three groups, the first containing a pair of counter-rotating, quasi streamwise vortices and high streak-like shear stresses; the second group is characterized by multiple streamwise vortices and little variations in wall stress; and the third group has no buffer layer structures.     

Coherent structures and riblets

This work deals with the effect of the riblets on the coherent structures near the wall. The emphasis is put on the genesis of the quasi-streamwise vortices in the presence of the riblets. The quasi-streamwise vortices regenerate by the tilting of wall normal vorticity induced by prevailing structures. This requires a mechanism which leads to a temporal streamwise dependence near the elongated flow structures and to a subsequent formation of new wall normal vorticity. It is suggested here that the action of existing quasi-streamwise vortices on the sidewalls of wall normal vorticity may create a local, streamwise dependent spanwise velocity and therefore, a secondary wall normal vorticity field. A preliminary analysis of the set-up and the time and space development of this secondary three-dimensional flow associated with the regeneration mechanism, is given. An attempt is made, in order to explain the drag reduction performed by the riblets through an intermittent model, based on the protrusion height. Logical estimates of the amount of drag reduction are obtained. The differences between the mechanism suggested here and those based on forced control experiments are also discussed.     

Investigation of particle-laden turbulent pipe flow at high-Reynolds-number using particle image/tracking velocimetry (PIV/PTV)

An experimental investigation of a high Reynolds number flow (Re = 320 000) of a dilute liquid-solid mixture (<1% by volume) was conducted. The turbulent motion of both the liquid phase (water) and particles (0.5, 1, and 2 mm glass beads) was evaluated in an upward pipe flow using a particle image/tracking velocimetry (PIV/PTV) technique. Results show that the Eulerian mean axial velocity of the glass beads is lower than that of the liquid phase in the central region but higher in the near-wall region. Moreover, the presence of the coarse particles has a negligible effect on the turbulence intensity of the liquid phase. Particles show higher streamwise and radial fluctuations than the liquid-phase at the tested conditions. The profiles of particle concentration across the pipe radius show almost constant concentration in the core of the pipe with a decrease towards the near wall region for 0.5 and 1 mm particles. For the 2 mm particles, a nearly linear concentration gradient from centre to the pipe wall is observed. The results presented here provide new information concerning the effect of a dispersed particulate phase on the turbulence modulation of the liquid carrier phase, especially at high Reynolds numbers. The present study also demonstrates how correlations developed to determine if particles cause turbulence attenuation/augmentation are not applicable for solid-liquid flows at high Reynolds numbers. Finally, the importance of particle-fluid slip velocity on fluid phase turbulence modulation is illustrated.     

DEM investigation on conveying of non-spherical particles in a screw conveyor

Screw conveyors are extensively used in modern industry such as metallurgy, architecture and pharmaceutical due to their high-efficiency in the transportation of granular materials. And substantial efforts have been devoted to the study of the screw conveyors. Numerical method is an effective way to study screw conveyor. However, previous studies have mainly focused in the regime of spherical particles while the in-depth investigations for non-spherical particles that should be the most encountered in practical applications are still limited. In view of the above situations, discrete element method (DEM), which has been widely accepted in simulating the discrete systems, is utilized to investigate the conveying process of non-spherical particles in a horizontal screw conveyor, with particles being modeled by super-ellipsoids. In addition, a wear model called SIEM (Shear Impact Energy Model) is incorporated into DEM to predict the wear of screw conveyor. The DEM simulation results demonstrate that the particle shape is influential for the flow behaviors of particles and the wear of conveyor. The conveying performance evaluated quantitatively of both mass flow rate and power consumption is subsequently obtained to investigate the effect of sphericity of particle with different operation parameters. Moreover, particle collision frequency and collision energy consumption are acquired to investigate the possible particle breakage between particles and screw blade. The comparisons between particle–particle collision and particle–wall collision reveal that particles with large shape index have more possibility to be damaged in particle–wall impingement.     

纵向涡对近壁湍流影响的激光测试

用激光测试方法详细地研究了单涡、上洗双涡、下洗双涡对近壁湍流的影响．在涡的下洗区，流体流速高，湍流度低，流体流向壁面；在涡的上洗区，流体流速低，湍流度高，流体远离壁面     

Numerical simulation of tangentially injected turbulent swirling flow in a divergent tube

This paper studies the properties of turbulent swirling decaying flow induced by tangential inlets in a divergent pipe using the realizable k–ε turbulence model and discusses the effects of the injector pressure and injection position. The results of transient solutions show that both the recirculation zone near the wall in upstream of the injectors and the vortex breakdown in downstream of the injectors increase in size during the whole period. A nearly axisymmetric conical breakdown is formed and its internal structure consists of two asymmetric spiral‐like vortices rotating in opposite directions. The stagnation point shifts slowly toward the pipe outlet over time. The maxima of the three velocity components, which are located near the wall, decrease gradually with streamwise direction. It can also be inferred that Mach number approaches 1.0 near the injector outlets. The velocities increase with the increasing injector pressure. However, its increasing trend is not significant. With the increase of the injection position, vortex breakdown moves in downstream direction and the pitch along the streamwise direction increases. Copyright © 2008 John Wiley & Sons, Ltd.