工业中粉体颗粒的荷电机理及数值模拟方法

工业过程中粉体颗粒不可避免地会相互摩擦碰撞而荷电. 荷电颗粒的存在可能会危害正常的工业生产过程, 也可能对工业过程起促进作用. 因此, 荷电粉体颗粒及其特性受到了广泛的关注, 但目前对粉体颗粒的荷电机理依然缺乏透彻的了解, 尤其是在气固两相流动中的粉体颗粒荷电现象. 事实上, 工业中存在的粉体颗粒的运动都受到流体的影响, 是典型的气固两相流系统, 流体对粉体颗粒的作用使粉体颗粒接触的荷电现象变得更为复杂, 因此从两相流动的观点来研究粉体颗粒荷电的物理本质就显得越来越重要. 本文介绍了工业过程中的几种不同类型的粉体颗粒荷电行为, 回顾了颗粒的荷电机理与描述颗粒荷电的数学模型. 对于工业过程中颗粒的荷电现象及颗粒在多相流体中的动力学行为, 介绍了研究颗粒受流体影响时荷电特性的数值模拟方法. 本文旨在对粉体颗粒的荷电机理、应用以及研究方法进行梳理与探讨, 为正确认识工业过程中粉体颗粒的荷电现象并加以控制利用提供理论借鉴.     

椭球颗粒搅拌运动及混合特性的数值模拟研究

为探讨在强制搅拌下同属性颗粒由分层到分布均匀状态的运动特征及规律, 本研究利用三维离散单元法模拟不同转速下U形罐体内等粒径椭球颗粒的混合过程. 从单颗粒随机运动轨迹、宏观颗粒流运动矢量图的角度分析颗粒混合过程的宏观混合规律及局部混合特征, 定量描述混合度与搅拌叶片旋转圈数的数学关系. 结果表明, 强制搅拌下同属性分层颗粒的混合是在对流混合及四个局部混合共同作用下实现的; 分层颗粒的混合度与搅拌轴的转速无关, 而与搅拌轴旋转圈数直接相关; 混合度与圈数的关系符合指数增长模型. 研究结果可为散体物料增混行业的设备改进及操作控制提供依据和参考.     

三维光滑粒子流体动力学并行计算程序

介绍了基于消息传递并行程序设计平台研制的三维光滑粒子流体动力学并行程序CSPH3D.包括计算格式、并行方案、并行程序逻辑,以及加快邻域粒子搜索的处理方法.对三维微喷射和斜侵彻的计算表明:CSPH3D程序可以较好地计算这类问题.并且程序具有较高的并行效率.对于粒子总数为1527402的微喷射算例和粒子总数为1454225的斜侵彻算例,使用100个处理器时,并行效率可以达到80%.     

Simulation of ultrasonic attenuation in theelastic mixing particle system

For the study of predicting ultrasonic attenuation of elastic, spherical mixing particles in the liquid-solid two-phase system, the Monte Carlo method(MCM) is introduced,serving as a probability and statistics technique to evaluate the inside ultrasonic events during the ultrasound propagation. On the basis of ultrasonic scattering and aborption, the continuous ultrasonic waves are represented as discrete and independent phonons. By recording the scattering events, tracing the trajectory of a moving phonon and calculating the number of phonons that finally reach the receiving transducer, the ultrasonic attenuation coefficient is obtained to be a frequence-dependent spectrum. Numerical investigations have been carried out to predict and compare the ultrasonic attenuation for a solid-liquid two-phase system with a single type particle. After verifing its feasibility, such a method is then appalied into mixing particle system, where the mixing iron particles and glass beads with various ratios are set as examples for the purpose of predicting ultrasonic attenuation for the monodisperse and polydisperse mixing particle systems. The results of MCM, the ECAH model, the Lloyd Berry(LB) model and the Waterman model match well when the particle volume concentration is lower than 10%,corresponding to iron particles and glass beads respectively. In the case of two-phase system with mixing particles, it is shown that as the particle volume concentration increases to 10%, the variation of the ultrasonic attenuation coefficient with mixing ratio yields a nonlinear tendency.The physical properties of particles can also influence ultrasonic attenuation significantly.     

基于拟颗粒的气固两相曳力模型研究

为了研究气固两相流动大涡模拟中合适的曳力计算模型,本文引入拟颗粒和拟颗粒表面能的概念,通过拟颗粒表面能与外界输入能量之间的平衡关系来确定拟颗粒的粒径。根据拟颗粒粒径,得到运算量较小且考虑颗粒团聚效应的曳力计算模型。应用本文的曳力计算模型对二维竖直槽道内稠密气固两相流动进行了大涡模拟,结果表明颗粒的浓度分布具有上稀下浓,壁面附近浓中心稀及颗粒聚集等特点。这与实验结果在定性上是一致的。对气相和颗粒相的瞬时速度场进行了分析,发现气相和颗粒相速度场分布的非对称性是形成颗粒浓度分布壁面附近浓中心稀的重要原因之一。     

高速流动PIV示踪粒子跟随性分析

通过理论推导提出了一种评价高速流动PIV示踪粒子随流能力的松弛特性分析模型,在法向Mach数大于1.4时具有良好的适用性.将新模型应用于试验测量,发展了高速流动PIV系统和示踪粒子布撒技术,验证了高速流动PIV的定量化测量能力.针对空间发展的二维超声速气固两相混合层,数值模拟了不同Stokes数和对流Mach数（M_c）下的粒子跟随性以及弥散和迁徙运动,结果表明:相同对流Mach数,粒径越小的示踪粒子跟随性越好,Stokes数在[1, 10]范围内的粒子有最大扩散距离.示踪粒子的直径大小决定其在超声速混合层大涡拟序结构中的分布特征,且粒径越小,气体与粒子的掺混越剧烈.相同粒径的粒子,对流Mach数越大跟随性越差.      

流速比对颗粒在液固两相平板混合层中沉降的影响

本文采用流动显示的方法对平板混合层中上下层流体速度比对固体颗粒在混合层中沉降的影响进行了研究。实验中分别采用粒径小于 ４０ ｕｍ,粒径 ９８～１０４ｕｍ,粒径 １５４～１６０ｕｍ的玻璃微珠以及环氧树脂作为固相颗粒,对这些颗粒在速度比分别为１：１．２、１：２和１：２．８的液相混合层中的运动进行了显示。结果表明混合层中大涡结构对固体颗粒的沉降具有迟滞作用,其作用程度取决于混合层中上下层流体速度比。速度比越大,颗粒的沉降越慢。     

薄滚筒内二元湿颗粒体系混合行为的离散单元模拟研究

采用所建立的四方程线性弹性-阻尼离散单元模型,同时考虑了液桥力的作用,对填充量为40%、含液量为3%的水平薄滚筒内S型(不同直径颗粒)二元湿颗粒体系混合过程进行了数值模拟,并与同等操作条件下不含液的干颗粒体系的混合行为进行了比较,分析了液体对颗粒体系混合行为的影响. 同时还将计算结果与文献中的实验结果进行了比较. 结果表明,由于湿颗粒间液桥力的牵引作用使不同性质的颗粒不易分离,使部分颗粒聚结成团,减弱了离析作用的影响,使得滚筒内湿颗粒的混合程度高于相同条件下的干颗粒体系,且接触力的分布较干颗粒体系更加均匀. 通过对混合过程的模拟,直观地反映了混合过程中颗粒的微观运动特性和内部的力学结构,为研究湿颗粒体系混合过程机理提供了依据和参考. 关键词： 湿颗粒体系 混合 分离 离散单元法     

气固两相混合层流场双向耦合的数值研究

采用双向耦合模型对有涡配对的二维气固两相混合展数值模拟,在考虑颗粒对流场反作用基础上进一步对颗粒间通过流体的相互作用进行分析。流场用拟谱方法求解,颗粒用颗粒轨道模型跟踪。结果发现,流场中大涡卷起和配对仍居主导地位;颗粒Ｓｔ数为Ｏ（０．１）～Ｏ（１）时,颗粒减弱了流场雷诺应力强度,加快涡量扩散; Ｓｔ数为Ｏ（１）时,颗粒分布极不均匀,主要集中在渴的边缘．     

复合管道内液-固两相流浆体水击压强计算

对于单一管道的单相液体水击现象计算已经比较成熟,但是很少有对于复合管道的液固两相流系统的浆体水击 计算方面的报道。根据实际工程设计需要,提出了具有复合管道的液-固两相流浆体水击压强的计算方法。同时对某电厂 液一固两相流输灰系统浆体水击压强进行了计算与分析,提出了防止浆体水击的措施。     

突扩圆管内液-固两相流固体颗粒运动特性的DPM数值模拟

采用考虑颗粒碰撞的欧拉一拉格朗日数值模拟方法(DPM),对水平突扩圆管中液固两相流固体颗粒的碰撞过程进行了数值计算.在模型中,对液相采用欧拉法建立控制方程,对离散颗粒采用拉格朗日方法模拟.采用硬球模型描述颗粒间的碰撞作用.计算结果表明,该模型可以真实地模拟液固两相流中固体颗粒运动的动态变化过程以及颗粒的非均匀分布特征,从单颗粒层次上提供颗粒的运动信息,这有助于深入研究液固两相流中固体颗粒的运动规律.     

混合弹性颗粒体系声衰减数值模拟

采用概率统计方法——蒙特卡罗方法,建立一种预测液固两相体系中混合弹性球形颗粒声衰减的理论模型。在单颗粒声散射和吸收的基础上,将连续超声波抽象离散化为大量独立的声子,追踪声散射过程,通过统计接收器探测声子数最终确定声衰减系数。采用数值方法对单一球形颗粒的液固两相体系中声衰减进行预测和比较,确定了该方法的可行性后将该方法推广到混合颗粒体系中,对玻璃微珠/铁粉构成的混合颗粒及多分散混合颗粒体系进行数值研究。结果表明:在体积浓度低于10%时,蒙特卡罗法预测得玻璃微珠或铁粉颗粒声衰减和ECAH,Lloyd和Berry,Waterman等建立的模型结果吻合。对于混合颗粒构成的两相体系,算例中,随着体积浓度增大到10%,声衰减系数随混合颗粒数目比的变化呈现出了非线性的变化,同时颗粒物性也会影响不同组分颗粒对声衰减的贡献,算例中铁粉颗粒比玻璃微珠对声衰减的影响更大。      

Critical Angle Refractometry for Simultaneous Measurement of Particles in Flow: Size and Relative Refractive Index

The principle of the optical technique critical angle refractometry, used to determine the size and refractive index of spherical particles (with relative refractive index below unity) in liquid flows, was investigated. This technique is based on the observation of the particle scattering pattern around the critical angle. Similarly to the recent technique developed for rainbow scattering pattern analysis for droplet temperature and size measurements, it is shown that the relative particle refractive index (m_r<1) and size can be determined from the position of the primary diffraction fringe and from the angular spacing between two fringes. Explicit equations for refractive index and particle size measurement were derived from the first-order term of the physical optics approximation. An experimental validation test and numerical computations based on the Lorenz-Mie theory were used to validate the principle of the proposed technique and to estimate its sensitivity, which was shown to be of the same order as that of the rainbow technique. This technique is considered to be useful for various applications in liquid multiphase flows where the particles size and material are to be characterized.     

Study of mass consistency LES/FDF techniques for chemically reacting flows

A hybrid large eddy simulation/filtered density function (LES/FDF) approach is used for studying chemically reacting flows with detailed chemistry. In particular, techniques utilised for ensuring a mass consistent coupling between LES and FDF are discussed. The purpose of these techniques is to maintain a correct spatial distribution of the computational particles representing specified amounts of fluid. A particular mass consistency technique due to Y.Z. Zhang and D.C. Haworth (A general mass consistency algorithm for hybrid particle/finite-volume PDF methods, J. Comput. Phys. 194 (2004), pp. 156–193) and their associated algorithms are implemented in a pressure-based computational fluid dynamics code suitable for the simulation of variable density flows, representative of those encountered in actual combustion applications. To assess the effectiveness of the referenced technique for enforcing LES/FDF mass consistency, two- and three-dimensional simulations of a temporal mixing layer using detailed and reduced chemistry mechanisms are carried out. The parametric analysis performed focuses on determining the influence on the level of mass consistency errors of parameters such as the initial number of particles per cell and the initial density ratio of the mixing layers. Particular emphasis is put on the computational burden that represents the use of such a mass consistency technique. The results show the suitability of this type of technique for ensuring the mass consistency required when utilising hybrid LES/FDF approaches. The level of agreement of the computed results with experimental data is also illustrated.     

Mixing and segregation of granular materials in chute flows

Mixing of granular solids is invariably accompanied by segregation, however, the fundamentals of the process are not well understood. We analyze density and size segregation in a chute flow of cohesionless spherical particles by means of computations and theory based on the transport equations for a mixture of nearly elastic particles. Computations for elastic particles (Monte Carlo simulations), nearly elastic particles, and inelastic, frictional particles (particle dynamics simulations) are carried out. General expressions for the segregation fluxes due to pressure gradients and temperature gradients are derived. Simplified equations are obtained for the limiting cases of low volume fractions (ideal gas limit) and equal sized particles. Theoretical predictions of equilibrium number density profiles are in good agreement with computations for mixtures of equal sized particles with different density for all solids volume fractions, and for mixtures of different sized particles at low volume fractions (nu<0.2), when the particles are elastic or nearly elastic. In the case of inelastic, frictional particles the theory gives reasonable predictions if an appropriate effective granular temperature is assumed. The relative importance of pressure diffusion and temperature diffusion for the cases considered is discussed. (c) 1999 American Institute of Physics.     

方腔内微细颗粒物运动特性的格子Boltzmann方法模拟

对Masselot and Chopard提出的模拟气固两相流动的格子Boltzmann-格子气(LBE-LGA)方法进行推广,能够反映两相间拖曳作用.利用该方法研究封闭方腔内的气固两相流运动特性,分析斯托克斯数St和模拟颗粒数目对颗粒群运动的影响,并与文献结果进行比较,表明LBE-LGA方法模拟颗粒运动是可行的.     

不同振动模式下颗粒分离行为的数值模拟

利用三维离散元法对垂直方向上的直线、圆和椭圆振动模式颗粒分离过程进行了数值模拟研究,对直线振动时上层大颗粒的波动及圆和椭圆振动时出现的聚集、循环等现象的形成机理进行了分析,并讨论了振动强度对各振动模式下颗粒分离形态的影响规律. 研究表明,综合运用空隙填充、侧面驱动的颗粒运动和能量非均匀分布三种机理,并结合颗粒群的速度矢量分布情况能够较好地解释各振动模式下的颗粒分离行为. 振动强度对圆和椭圆振动模式的分离形态具有显著的影响,并在振动强度约为3时,各种振动模式均具有良好的颗粒分离效果和稳定的颗粒运动状态. 关键词： 振动模式 颗粒分离 离散元法 数值模拟     

不同振动模式下颗粒分离行为的数值模拟

利用三维离散元法对垂直方向上的直线、圆和椭圆振动模式颗粒分离过程进行了数值模拟研究,对直线振动时上层大颗粒的波动及圆和椭圆振动时出现的聚集、循环等现象的形成机理进行了分析,并讨论了振动强度对各振动模式下颗粒分离形态的影响规律. 研究表明,综合运用空隙填充、侧面驱动的颗粒运动和能量非均匀分布三种机理,并结合颗粒群的速度矢量分布情况能够较好地解释各振动模式下的颗粒分离行为. 振动强度对圆和椭圆振动模式的分离形态具有显著的影响,并在振动强度约为3时,各种振动模式均具有良好的颗粒分离效果和稳定的颗粒运动状态.     

颗粒圆孔射流碰撞并行直接数值模拟算法研究

本文对空间模式发展的颗粒圆孔射流碰撞进行了并行直接数值模拟算法研究。气相采用可压缩的N-S方程直接求解。颗粒相采用Lagrangian方法跟踪实际的颗粒运动。利用并行求解算法,实现了颗粒穿越边界面的模拟以及高效颗粒碰撞算法。考虑了颗粒和流体的双相耦合以及颗粒之间的碰撞。在本文的计算条件下,颗粒的直径远小于网格的间距,平均的Kolmogorov尺度和网格的间距在一个量级。气相和颗粒相的应力与实验的对比研究表明,本文的颗粒并行程序是可信的。     

GPU accelerated simulations of bluff body flows using vortex particle methods

We present a GPU accelerated solver for simulations of bluff body flows in 2D using a remeshed vortex particle method and the vorticity formulation of the Brinkman penalization technique to enforce boundary conditions. The efficiency of the method relies on fast and accurate particle-grid interpolations on GPUs for the remeshing of the particles and the computation of the field operators. The GPU implementation uses OpenGL so as to perform efficient particle-grid operations and a CUFFT-based solver for the Poisson equation with unbounded boundary conditions. The accuracy and performance of the GPU simulations and their relative advantages/drawbacks over CPU based computations are reported in simulations of flows past an impulsively started circular cylinder from Reynolds numbers between 40 and 9500. The results indicate up to two orders of magnitude speed up of the GPU implementation over the respective CPU implementations. The accuracy of the GPU computations depends on the Re number of the flow. For Re up to 1000 there is little difference between GPU and CPU calculations but this agreement deteriorates (albeit remaining to within 5% in drag calculations) for higher Re numbers as the single precision of the GPU adversely affects the accuracy of the simulations.