基于格子玻耳兹曼的等离子体射流中粉末热运动数值模拟

为了探索快速有效的数值计算方法用于研究等离子体射流中粉末群的热运动状态,本文采用格子玻耳兹曼方法计算等离子体射流,使用随机算法跟踪颗粒,模拟了等离子体射流中粉末的加热和运动。计算结果表明:颗粒温度的变化比速度快,颗粒高速区比高温区的范围大。离射流出口越近,沿射流横截面颗粒轴向平均速度和温度波动越大,能获得沿射流横截面方向温度速度较高且分布均匀的位置区间为距离射流出口150 mm~200 mm。实测的粒子温度与速度和计算结果吻合较好,验证了本模型的有效性。论文研究为等离子体加工工艺条件的制定提供了依据。     

消息与动态(九)

一种计算等离子喷涂颗粒吸收热量的方法计算等离子喷涂过程中颗粒吸收热量的方法一般直接引用能量和质量转移联立方程的解.可是被喷涂的颗粒材料在等离子射流中温度变化十分激烈,变化范围太宽,以致 ...     

消息与动态

一种计算等离子喷涂颗粒吸收热量的方法计算等离子喷涂过程中颗粒吸收热量的方法一般直接引用能量和质量转移联立方程的解.可是被喷涂的颗粒材料在等离子射流中温度变化十分激烈,变化范围太宽,以致     

消息与动态(九)

一种计算等离子喷涂颗粒吸收热量的方法计算等离子喷涂过程中颗粒吸收热量的方法一般直接引用能量和质量转移联立方程的解.可是被喷涂的颗粒材料在等离子射流中温度变化十分激烈,变化范围太宽,以致 ...     

水力喷砂射孔机理及力学特性分析

通过分析水力喷砂射孔机理,建立了磨料射流形成过程的动力学模型,推导出了液相的加速表达式以及磨料颗粒在喷嘴内的运动方程;基于塑性材料破坏的第四强度理论和岩石破坏的库伦-摩尔准则,建立了球形磨料颗粒对套管、岩石打击过程的力学模型,推导出了磨料颗粒的极限速度表达式;借助MATLAB软件,使用迭代法对上述力学模型进行了计算验证。研究结果表明:在喷嘴结构参数和颗粒直径不变时,液相速度越大、加速度越大,则磨料颗粒的加速也会越快;而磨料颗粒对材料的冲击作用力不仅与靶件本身的性能参数有关,还与颗粒的冲击速度和密度呈正相关。     

纳米W粉冲击烧结的分子动力学模拟

粉末冲击烧结是制备高品质W的一种有效方法，而分子动力学方法在尺度极小、过程迅速的数值模拟上有着独特的优势。因此运用分子动力学方法，结合W的嵌入原子势，对常温下的纳米W粉末的冲击烧结过程进行模拟，得到颗粒微观压实过程图、体系速度分布云图、p-U_p、T-U_p、T-p曲线以及径向分布函数。研究了不同颗粒速度及产生的射流对纳米W粉末冲击烧结影响，分析了微观冲击烧结机理。结果表明，低速冲击条件下（500 m/s以下），纳米颗粒无法压实。高速条件下（1 000 m/s及以上），颗粒能获得致密化很高的压实。颗粒间的相互挤压造成的高应力使颗粒表面的原子发生流动变形，原子向颗粒间空隙流动，形成压实。颗粒间产生的射流以及高速冲击导致的颗粒熔化，均促进烧结获得致密度更高的烧结体。     

气固射流加砂方式中颗粒速度场及浓度场的实验研究

将激光数字式粒子成像测速技术(DPIV)方法应用于砂尘环境实验中颗粒速度场及浓度场的实验研究,成功地对环境风洞中颗粒速度场及浓度分布进行了无干扰瞬态测量,实验结果较真实地反映了气固射流加砂方式中粒子速度场及浓度场的分布规律。文中所获得的规律性认识为我国自行研制大型砂尘实验设备,确定加砂/尘方法及实验段参数等提供了依据。     

逆壁射流中射流剪切层内湍流结构特性的实验研究

采用粒子图像测速技术对逆壁射流全流场进行了实验测量,射流与主流的速度比为8.89,基于射流圆管内径的雷诺数为对射流中心线上不同流向位置的脉动速度场统计分析发现:在占据主导地位.在驻点附近(尺度进行分析,在x/D=46～51总尺度向射流下游发展呈减小趋势.在x/D=35～41,参考点下游尺度大于上游尺度.在本征正交分解方法对湍流结构进行了定量分析,发现模态能量集中在低频,流场中能量最大的模态频率为f D/Uj=0.000 5,出现在再循环区.频率为产生了湍流结构,并且沿再循环区外围输运.高频结构的构型是类似的,均位于射流剪切层内,且频率越高,越接近射流出口,尺度越小.     

气固两相流场的湍流颗粒浓度理论模型

本文进行了气固两相流动颗粒湍流扩散现象的理论分析，提出了颗粒湍流扩散系数和气流弥散效应二个颗粒湍流模化新概念，在此基础上建立了气固两相流场湍流颗粒浓度模型。理论模型包括离心力和其它外加力场作用下颗粒运动和浓度分布的计算方法。运用湍流颗粒浓度模型，对直管气固两相流动、受限射流气固两相流动和９０°弯管气固两相流动等三种流动做了数值模拟，计算获得颗粒速度、颗粒浓度等主要流动参数。讨论了湍流颗粒浓度模型的适用性。     

水动力作用下低速射流碎化的数值模拟

为研究射流在水动力作用下的碎化特性,采用有限体积法对轴对称坐标下Navier-Stokes方程进行了求解,考虑重力和表面张力的影响,并通过Volume-of-Fluid法与Level-Set法成功捕捉到界面的不稳定发展、变形及射流碎化过程,分析了流场内部速度场和压力场分布,结果表明,射流碎化长度随Re/We“5数呈指数型增加,最后探讨了射流速度、直径及周围流体密度、粘性等参量对射流的碎化过程的影响规律.     

柔性微粒介电泳分离过程的多尺度模拟

介电泳分离是一种高效的微细颗粒分离技术,利用非均匀电场极化并操纵分离微流道中的颗粒. 柔性微粒在介电泳分离过程中同时受多种物理场、多相流和微粒变形等复杂因素的影响,仅用单一的计算方法对其进行模拟存在一定的难度,本文采用有限单元——格子玻尔兹曼耦合计算的方法处理这一难题.介观尺度的格子玻尔兹曼方法将流体看成由大量微小粒子组成,在离散格子上求解玻尔兹曼输运方程,易于处理多相流及大变形问题,特别适合模拟柔性颗粒在介电泳分离过程中的变形情况.另一方面,介电泳分离过程的模拟需求解流体、电场和微粒运动方程,计算量相当庞大,通过有限单元法求解介电泳力,提高计算效率.利用这种多尺度耦合计算方法,对一款现有的介电泳芯片分离过程进行了模拟.分析了微粒在电场作用下产生的介电泳力,揭示了介电泳力与电场变化率等因素之间的关系.对微粒运动轨迹及其变形的情况进行了研究,发现微粒的变形主要与流体剪切作用有关.这种多尺度耦合计算方法,为复杂微流体的计算提供了一种有效的解决方案.      

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…     

湍流两相流的脉动速度联合PDF输运方程

概率密度函数（PDF）的方法是构造两相湍流模型的一种重要的方法.构建气体-颗粒速度联合PDF输运方程的关键是颗粒所见气体微团速度的Langevin方程.首先由N-S方程出发,精确推导出颗粒所见气体微团脉动速度的Langevin方程,进而通过理论分析表明,对比通常采用的颗粒所见气体微团瞬时速度的Langevin方程而言,采用前者能有效地减少关联量的统计偏差.最后,给出了颗粒-气体脉动速度的联合PDF输运方程.     

Numerical modeling of electrostatic spray painting transfer processes in rotary bell cup for automotive painting

Here we implemented a 3D comprehensive Eulerian-Lagrangian model in order to investigate the electrostatic spray transfer processes in the high-speed rotary bell sprayer. This efficient algorithm contains spray dynamics, airflow, paint droplets tracking and an electrostatic effect to simulate atomization. The algorithm is implemented using the OpenFOAM package. A solver for the particle trajectory was used to illustrate the process of spray transport and also the interaction of the airflow and the particle that is solved by momentum coupling. Creating an initial condition of the particle approach has been proposed that is matched with practical applications. The fluid-dynamics is simulated by solving the unsteady 3D compressible Navier-Stokes equations. Unsteady flow is computed by using a Large eddy simulation (LES) turbulence approach, while the motion of the particles is simulated by tracking the droplet size distribution approach. The model correctly predicts that the bell cup spin forces the paint particles to fall off from the bell surface towards the high-velocity airflow. The present work illustrates a tentative benchmark and contains a systematic analysis of the recirculation zone length, the toroidal vortex, the overspray phenomena and the flowfield characteristics like mean velocity, pressure, turbulent kinetic energy and velocity fluctuation. The results indicate as dominant operating parameter the air-paint flow rate with voltage level deeply affecting the spray shape. A more uniform distribution of the coating is obtained by growing this high-velocity shaping airflow, although the values of the transfer efficiency (TE) are reduced. The distribution of the particle size is very sensitive to changes in the rotational speed. Experimental results obtained in this study put forward a clear link between the shaping air flow rate and the rotation frequency on the aerodynamics and also provide valuable insights to design modern ERBS. The paint spray distribution obtained in the present work is validated against coating experimental results with suitable accuracy.     

Experimental and mechanistic study of dispersed micrometer-sized particle resuspension in a square straight duct with rough walls

The resuspension of graphite dust is an important phenomenon in the release of radioactivity and the safety of nuclear reactors during severe accidents. In this study, a visualization experimental platform is constructed to study effects of particle size, flow velocity, and wall roughness on the resuspension characteristics of graphite particles. A statistical model of particle resuspension applicable to monolayer dispersed particles is developed based on the moment equilibrium of the particles and the flow field characteristics, as calculated by the large-eddy simulation framework. The results show that particle resuspension can be divided into short- and long-term resuspension stages. Most particle resuspension occurs during the short-term stage. With increases in flow velocity and particle diameter, the aerodynamic or adhesion force acting on the particles increases, and corresponding particle resuspension fraction increases. The influence of rough walls on particle resuspension is related to both the force on the particles and the arm ratio between the wall morphology and the particle diameter. A comparison with the experimental results demonstrates that the particle resuspension model developed in this study accurately predicts the impact of flow velocity, particle size, and wall roughness on particle resuspension.     

Simulation of turbulent dispersion using Weierstrass modes

A simple fractal model is proposed for the dispersion of passive scalars in an incompressible homogeneous turbulent flow field. The dispersion process is based on a three-dimensional velocity field which is assumed to be a linear superposition of Taylor-Green vortices with wave numbers and amplitudes as those in a Weierstrass function. The chosen velocity field satisfies the continuity equation as well as Kolmogorov's inertial range power law, and has a fractal dimension D between upper and lower length-scale bounds. A cloud of tracer particles is relased into the flow field and dispersed by the fluid motion, the trajectories are numerically integrated from the velocity field. The puffs which result from this process are used to examine some aspects of turbulent dispersion, through comparisons with integrated concentration wind-tunnel measurements. The agreement between numerical and experimental results indicates the significance of the proposed simulation model.     

Explosively driven particle fields imaged using a high speed framing camera and particle image velocimetry

A high speed framing camera and a particle image velocimetry instrument were used to determine the properties of explosively driven particle fields in early microsecond and later millisecond times. Test items were configured in a two inch long cylindrical shape with a half inch diameter core of organic explosive. The core was surrounded by a particle bed of aluminum or tungsten powder of a specific particle size distribution. Position data from the leading edge of the particle fronts for each charge was recorded with a high speed framing camera at early time and with a particle image velocimetry (PIV) instrument at later time to determine particle velocity. Using a PIV image, a velocity gradient along the length of the particle field was established by using the mean particle velocity value determined from three separate horizontal bands that transverse the particle field. The results showed slower particles at the beginning of the particle field closest to the source and faster ones at the end. Differences in particle dispersal, luminescence, and agglomeration were seen when changes in the initial particle size and material type were made. The aluminum powders showed extensive luminescence with agglomeration forming large particle structures while the tungsten powder showed little luminescence, agglomeration and no particle structures. Combining velocity data from the high speed framing camera and PIV, the average drag coefficient for each powder type was determined. The particle field velocities and drag coefficients at one meter showed good agreement with the numerical data produced from a computational fluid dynamics code that takes advantage of both Eulerian and Lagrangian solvers to track individual particles after a set post detonation time interval.     

Investigation of the acceleration of aluminum particles behind a shock wave using instantaneous Laser Doppler Velocimetry

The acceleration of aluminum particles with a 5μm diameter in the flow field behind an incident shock wave was investigated experimentally in a 10-m long and 70 mm inner diameter shock tube. By means of instantaneous Laser Doppler Velocimetry (LDV) the velocity of the particles was observed directly. The light scattered by the moving particles is Doppler shifted and sent to the laser Doppler velocimeter. The velocimeter essentially consists of a phase-stabilized Michelson interferometer used as a sensitive spectrometer. An electro-optical circuit ensures the phase stabilization that results in a voltage signal independent of the scattered light intensity and proportional to the mean velocity of the particles at the measurement point. Because of the very short response time (1μs) of the LDV system used here, the latter gives a continuous real-time signal of the particle acceleration. To avoid particle oxidation the particles were accelerated by a high-speed nitrogen gas flow. From the measured velocity the dimensionless drag coefficient was calculated. The drag coefficient is related to the fluid dynamic force exerted by the gas on the particles. The experimental data were compared to theoretical models from the literature. A significant deviation between the model and the experimental data was observed. This deviation is supposed to be induced by the shock wave, which hits the particles and breaks them into pieces of a smaller diameter. Further experiments will be carried out in the future to check the size distribution of the particles after the shock has gone past them.      

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

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

Experimental and numerical study of installing a dune model in a 90° bend in the horizontal-vertical pneumatic conveying system

In the pneumatic conveying process, particles move to the bend under the influence of inertia to form a particle rope, which will cause serious wear between the particles and the pipe wall, and then the dune model is designed and installed in the 90° bend to reduce energy consumption and wear in this study. Firstly, the minimum pressure drop velocity of particles transported by different size dune models was obtained through experimental study. Then the energy saving mechanism of the dune model is studied by CFD-DEM coupling. The experimental results show that the installation of the dune model reduces the minimum pressure drop velocity. The numerical simulation results show that the number of collisions between the particles and the tube wall in the vertical tube decreases after the installation of the dune model, which reduces the energy loss. Moreover, the increasing of tail size of the dune model is beneficial to the diffusion and acceleration of the particles in the vertical tube.