竖直振动管中颗粒毛细效应的离散元模拟

将一根细管插入填充有颗粒的静止容器中并对管施加竖直振动,颗粒将在管内发生上升运动,并最终稳定在一定高度,这一现象与液体毛细效应类似,被称为颗粒毛细效应.为探究颗粒毛细效应过程中伴随的颗粒尺度动力学行为及机理,基于离散元方法建立颗粒运动模型,对颗粒毛细效应动力学过程和特性开展数值模拟研究.模拟再现了文献中实验得到的颗粒毛细效应全过程,给出了管内颗粒柱高度随时间的演变规律,结果表明,受到颗粒系统参数的影响,本模拟条件下颗粒毛细效应过程呈现单周期上升、倍周期上升和倍周期稳定三个阶段,在倍周期上升阶段颗粒柱上升速度逐渐减小,平缓过渡到稳定阶段.在此基础上,分析了管内颗粒速度场和填充率分布随时间的演变特性,揭示了颗粒毛细效应过程中由容器传输到管内的颗粒的占比分布.研究发现,管内不同高度位置颗粒的运动并不同步,随着管的振动,管内出现速度波,速度波的传播引起管内颗粒出现膨胀和压缩交替的情况,从而管内颗粒填充率随时间发生周期性波动;在上升阶段,越接近管壁由容器传输到管内的颗粒占比越大,在稳定阶段,管内上层颗粒的对流引起容器传输到管内的颗粒占比发生反转.     

颗粒毛细效应影响因素的离散元分析

颗粒毛细效应是指将一根细管插入填充有颗粒物质的容器中并对管施加竖直振动时颗粒在管内上升并最终达到一个稳定的高度的现象,该现象为颗粒物料的逆重力输运提供了一种潜在的技术途径.为探究颗粒毛细效应的影响因素,采用离散元方法,模拟再现了颗粒毛细效应过程,展示了不同管径下颗粒竖直方向速度演变特性,考察了不同容器宽度和振动条件下颗粒最终毛细上升高度随管径的演变规律.结果表明,在容器宽度与粒径比为40、管振幅与粒径比为14.33、管振动频率为12 Hz情况下,管径与粒径比D/d=3.33时,管内颗粒堵塞严重,使得颗粒上升缓慢,并造成颗粒柱中断; D/d=8.33时,起初毛细上升高度增加迅速,随后毛细上升高度的增大逐渐减缓,管内颗粒在管径方向几乎不存在速度梯度; D/d=15时,随着颗粒毛细上升高度的增大,管内颗粒柱分离为速度截然不同的两层,上层颗粒在管径方向几乎不存在速度梯度,而下层颗粒存在明显的速度梯度.研究还发现,在毛细效应能够发生的管径范围内,存在一个对应于颗粒最终毛细上升高度最大值的临界管径,当管径小于临界管径时,颗粒最终毛细上升高度随管径的增大而增大,当管径大于临界管径时,颗粒最终毛细上升高度随管径的增大而趋于减小;增大容器宽度,临界管径有所增大;增大振幅、适当提高频率能够有效促进临界管径的增大.     

颗粒毛细效应影响因素的离散元分析

颗粒毛细效应是指将一根细管插入填充有颗粒物质的容器中并对管施加竖直振动时颗粒在管内上升并最终达到一个稳定的高度的现象, 该现象为颗粒物料的逆重力输运提供了一种潜在的技术途径. 为探究颗粒毛细效应的影响因素, 采用离散元方法, 模拟再现了颗粒毛细效应过程,展示了不同管径下颗粒竖直方向速度演变特性, 考察了不同容器宽度和振动条件下颗粒最终毛细上升高度随管径的演变规律. 结果表明, 在容器宽度与粒径比为40、管振幅与粒径比为14.33、管振动频率为12 Hz情况下, 管径与粒径比$D/d = 3.33$时, 管内颗粒堵塞严重, 使得颗粒上升缓慢,并造成颗粒柱中断; $D/d = 8.33$时, 起初毛细上升高度增加迅速, 随后毛细上升高度的增大逐渐减缓, 管内颗粒在管径方向几乎不存在速度梯度; $D/d =15$时, 随着颗粒毛细上升高度的增大, 管内颗粒柱分离为速度截然不同的两层, 上层颗粒在管径方向几乎不存在速度梯度, 而下层颗粒存在明显的速度梯度.研究还发现, 在毛细效应能够发生的管径范围内, 存在一个对应于颗粒最终毛细上升高度最大值的临界管径, 当管径小于临界管径时, 颗粒最终毛细上升高度随管径的增大而增大, 当管径大于临界管径时, 颗粒最终毛细上升高度随管径的增大而趋于减小; 增大容器宽度,临界管径有所增大; 增大振幅、适当提高频率能够有效促进临界管径的增大.      

竖直振动颗粒床对流机制的颗粒尺度实验研究

为了得到竖直振动颗粒床形成对流的运动模式及形成机制,本文通过高速摄影技术对竖直振动颗粒床进行了实验研究。实验发现,随着振动加速度的增加,对流环覆盖的粒子层数和强度明显增加。通过分析颗粒速度矢量图的演化,可以获得控制对流运动的各种应力波在在粒子床中传播的信息,发现应力波的强度和持续时间与振动加速度密切相关。通过实验发现,对流运动发端于重力波面上粒子从侧壁向粒子床中心的不可逆跃迁,这种横向对流的强度与重力波面的曲率密切相关,而持续的时间随粒子床振动周期的变化而变化。     

基于LBM的铝微滴斜柱沉积水平偏移研究

金属微滴沉积制造技术采用逐点堆砌方式成型, 为斜柱沉积提供无支撑制造方式, 具有高度灵活性. 本文针对铝液滴斜柱连续沉积过程, 建立格子玻尔兹曼模型进行数值模拟, 研究液滴在凝固表面上的水平偏移运动. 根据表面能充放过程, 沉积运动被划分为下落、快速扩张、慢速扩张、回弹4个阶段, 其受力状态由表面能、重力势能、动能和黏性耗散趋势得到. 液滴内部流动在扩张阶段中表现为滑动状态, 而在回弹阶段中表现为滚动状态. 液滴偏移运动的加速阶段主要发生在扩张阶段, 而偏移距离则在回弹阶段中产生. 扩张阶段的受力状态表明偏移运动的主要推动力为重力和毛细力. 随着液滴轴线距离的增大, 扩张阶段中的加速段时间缩短、速度峰值提高, 使水平偏移距离呈先增大后减小的趋势, 这种阶段化特征源于加速段时长和速度极大值的竞争关系. 不同沉积高度和固液浸润度下, 偏移距离均保持相同的演化趋势. 在相同的轴线距离下, 偏移距离随固液浸润度的增大而减小, 随沉积高度的增大而减小. 通过拟合水平偏移距离演化规律、优化扫描步距, 能够实现斜柱的均匀沉积, 并使倾角与理论结果一致.      

受限颗粒体对制动系统非线性振动的影响

制动系统在工作时,往往受到沙粒、尘土以及磨屑等受限颗粒体的影响,这些受限颗粒体在摩擦副中的高度分布具有较强的随机性,一定程度诱发了制动系统的非线性振动. 本文中基于制动片切向振动模型,引入了新的受限颗粒体摩擦模型,提出了用波动系数来描述受限颗粒体高度分布随机性的强弱. 发现在特定参数下,当此系数为0时,制动片切向振动为周期运动;但是当此系数不为0时,制动片切向振动呈现拟周期或混沌运动,此时的切向振动分岔特性图的稳定轨道也会出现数量或分布的变化,甚至表现出混沌特性. 同一时变信号内,受限颗粒体引发制动片切向非线性振动包括发散、收敛以及拟周期运动等多种形式.      

探头在准二维颗粒介质中慢速运行的受阻实验研究

测量了圆锥探头缓慢压入准二维颗粒介质过程中所受阻力随深度的变化,发现阻力曲线在不同区域呈现不同的增长规律,曲线全过程存在两次增长规律的转变.本实验条件下观测到的现象表明:当填充颗粒总高度较高时,阻力曲线存在两个拐点,前部拐点是由颗粒回填所致,拐点深度约65mm(与探头尺寸有关),且不随颗粒填充高度变化,后部拐点是由于底部边界影响所致;后部拐点出现的深度Z2随着颗粒填充高度Zmax的增加线性增加,容器底部边界开始影响探头受力时探头距离底部的距离H也会随着颗粒填充高度Zmax的增加而增加,Z2-Zmax和H-Zmax关系拟合直线的斜率均约为0.5.     

台阶地形爆破振动放大与衰减效应研究

基于台阶地形爆破振动数值模拟与边坡振动监测实验，研究台阶地形爆破振动速度在传播过程中高程放大效应的产生及变化规律。结果表明，单个台阶坡顶质点的振动速度放大效应是在距爆源一定距离、达到一定高差的条件下产生的；坡顶质点振动速度放大倍数并不随台阶高度的增加而单调增加，在台阶高度超过某一临界值后，放大倍数随台阶高度的增加而减小。台阶高程对爆破振动速度既有放大作用，同时也随高度的增加产生衰减作用。根据模拟计算及现场观测数据分析结果，给出了台阶地形爆破振动速度预测模型，该模型为类似边坡工程的爆破地震波传播规律研究提供一定的参考。     

高速破片撞击充液容器拖拽阶段气腔特性研究

为研究液压水锤效应拖拽阶段的气腔特性,利用数值模拟与实验相结合的方法对破片撞击充液容器的过程进行研究,并分析了破片撞击速度和液体介质对液压水锤效应拖拽阶段气腔的影响。结果表明:破片撞击充液容器时,在液体中形成的气腔形状为圆锥形,其最大直径和长度随破片的运动逐渐增大,气腔长径比最终趋于一稳定变化区域,约在3.8~3.9之间;气腔最大直径随着破片撞击速度的增大而增大;柴油介质中形成气腔的最大直径和长径比变化规律与水介质中形成的相同,气腔长径比最终在4.25左右浮动,柴油介质中形成气腔的最大直径和长径比均大于水介质中形成的。     

具有间歇性缺陷的混合流体行进波对流斑图

本文通过流体力学基本方程组的数值模拟,探讨了具有中等Soret效应的混合流体行进波斑图的动力学特性.当分离比Ψ=-0.3时,首次发现一种没有源缺陷的左右相对传播的CPW(Counter propagating waves)状态向行进波状态的过渡形式,并且在r=1.50-1.60的范围内,行进波对流斑图中存在着间歇性缺陷结构.这种缺陷出现的周期随瑞利数r增大而增加.在缺陷出现的周期内,对流振幅也以行进波的周期在周期的变化,对流振幅的振动次数或行进波的周围数也随相对瑞利数r增大而增加.当r增加到1.65时,行进波对流斑图中的缺陷结构消失.由于缺陷引起的对流振幅的周期性变化也随之消失,而以行进波的周期在整个时间段上周期的振动.     

不同尺度分布散粒材料砂堆形成过程的二维离散元模拟

用作者开发的离散元程序,模拟不同尺寸分布的砂堆形成过程.把散体颗粒简化为圆形颗粒,模拟过程分三步:首先利用参考网格生成颗粒的松散堆积结构;为了避免颗粒下落的冲击作用对砂堆安息角的影响,先模拟颗粒在重力作用下在圆柱容器内的自由下落与堆积,直至堆积达到稳定;最后,移除容器,只保留一个底部边界,模拟颗粒体系的散落过程,直至形成一个稳定的砂堆.模拟结果表明,在其他参数保持相同的情况下,随着颗粒尺寸的减小,砂堆的安患角逐步减小并趋向于一常值.对模拟中的两组颗粒体系进行相同条件下的砂堆形成实验,结果表明,模拟与实验所得安息角大体相当.     

Interaction of a shock wave with a loose dusty bulk layer

The interaction of a planar shock wave with a loose dusty bulk layer has been investigated both experimentally and numerically. Experiments were conducted in a shock tube. The incident shock wave velocity and particle diameters were measured with the use of pressure transducers and a Malvern particle sizer, respectively. The flow fields, induced by shock waves, of both gas and granular phase were visualized by means of shadowgraphs and pulsed X-ray radiography with trace particles added. In addition, a two-phase model for granular flow presented by Gidaspow is introduced and is extended to describe such a complex phenomenon. Based on the kinetic theory, such a two-phase model has the advantage of being able to clarify many physical concepts, like particulate viscosity, granular conductivity and solid pressure, and deduce the correlative constitutive equations of the solid phase. The AUSM scheme was employed for the numerical calculation. The flow field behind the shock wave was displayed numerically and agrees well with our corresponding experimental results.      

Nonlinear dynamics of excited plate immersed in granular matter

The interaction between granular matter and the elastic body is a complex issue due to the complex properties of granular matter. An experiment involving a sinusoidally excited plate buried in glass bead particles contained in a box is conducted. The motion behavior of the plate is observed and recorded by the strain gauge. The amplitude–frequency and phase–frequency curves are recorded to study the natural property of the plate in granular matter. In this experiment, jump phenomena are found in both the amplitude–frequency and phase–frequency planes in circumstances with smaller particle sizes, lower buried depths, and larger amplitudes of the excitation force. Otherwise, the period-doubling bifurcation, especially 3T, is found with the increase in the excitation force. These bifurcations usually occur in specific buried depth and excitation frequency band and require smaller particle sizes. The experiments with random-shaped particles exhibit no-jump phenomenon, but period-doubling bifurcation and chaos. These phenomena are sensitive to parameters and closely related to the varying process of the excitation frequency and force. Reasonable mechanisms are summarized qualitatively through some of our recent researches in this paper.     

颗粒物质在竖直管中下行蠕动流动的实验测量

颗粒流蠕动行为是颗粒物质在竖直管中流动时常见的一种流动现象,其产生机理较复杂。为此,本文在在内径为150mm、高为5000mm的竖直管实验装置上,以FCC催化剂为固体颗粒物料,采用PV6型颗粒速度测量仪,测量不同颗粒流率下竖直管中的颗粒下行蠕动流动速度以及颗粒固含率,系统地考察了颗粒物质在竖直管中下行流动时的蠕动流动特性及产生机理。实验结果表明,颗粒物质在竖直管中下行流动时的流动行为可划分为两种形式。在颗粒流率较小时,颗粒物质下行速度呈现脉冲式变化,有速度停滞,可称之为蠕动I型流动。随着颗粒流率的增加,颗粒下行速度停顿消失,但仍是起伏变化,为蠕动II型流动。当颗粒流率增大到一定值后,颗粒物质下行蠕动行为消失,转变为流化流动。颗粒物质下行的蠕动行为是出口区颗粒成拱与崩塌、颗粒与器壁滑动摩擦和颗粒力链作用的综合反映。     

DEM analysis of the effect of electrostatic interaction on particle mixing for carrier-based dry powder inhaler formulations

Particle interactions play a significant role in controlling the performance of dry powder inhalers (DPIs), which mainly arise through van der Waals potentials, electrostatic interactions, and capillary forces. Our aim is to investigate the influence of electrostatic charge on the performance of DPIs as a basis for improving the formulation of the particle ingredients. The mixing process of carrier and active pharmaceutical ingredient (API) particles in a vibrating container is investigated using a discrete element method (DEM). The number of API particles attaching to the carrier particle (i.e., contact number) increases with increasing charge and decreases with increasing container size. The contact number decreases with increasing vibrational velocity amplitude and frequency. Moreover, a mechanism governed by the electrostatic force is proposed for the mixing process. This mechanism is different from that previously proposed for the mixing process governed by van der Waals forces, indicating that long-range and short-range adhesive forces can result in different mixing behaviours.     

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.     

FLOW OF A TRAIN OF DEFORMABLE FLUID PARTICLES IN A TUBE

I.IntroductionWangandSkalakll]studiedforsolidspheresintubes.HymanandSkalak12--31studiedthemovementofacompactfluidspheres(non-deformable)intubes.Theyonlystudiedcompactdropswhoseradiusesarenotbiggerthan0.7.Pozrikidisl4]studiedthenlotiollofdeformabletluidparticlesintubes.Inthispaper,westudiedtwo-phasenowregime,consistingofaperiodictrainofequallyspacedliquidslugsmovingaxisymmetricallyinacontinuousliquidphasewhichisimmisciblewiththeslugphase(Fig.l(a)).Wereportnumericalsolutionsforthehydrodynamic…     

Polysiloxane coatings on microspheres based on Multiphase Flow with Interface Exchange-Discrete Element Modelling

Polysiloxane coatings on yttria stabilized zirconia (YSZ) microspheres of 500 μm were simulated in a spout fluidized bed coater using Multiphase Flow with Interface Exchange-Discrete Element Modelling (MFiX-DEM). Two different coater configurations were developed to study the influence of gas velocity and its distribution on particle dynamics. The presence of the Wurster tube not only enhances the distribution but also increases the overall residence time of the particles. Investigations were also carried out with different Wurster tube positions (normal, 10% and 20% lowered from its initial position). Among these, 20% lowered Wurster tube position demonstrated the most effective coating process. The effects of gas inlet pressure on the average gas velocity and the distribution of particles were analyzed. More than 97% of the particles can be retained. The derived results, including average gas velocity, particle retention percentage, and distribution of particles with gas velocity, are being used to guide the experimental work in obtaining defect-free coatings for YSZ microspheres.     

Prediction of segregation characterization based on granular velocity and concentration in rotating drum

In a binary granular system composed of two types of particles with different granule sizes and the same density, particle sorting occurs easily during the flow process. The segregation pattern structure is mainly affected by the granular velocity and granular concentration in the flow layer. This paper reports on the experimental velocity and concentration measurement results for spherical particles in a quasi-two-dimensional rotating drum. The relationship between the granular velocity along the depth direction of the flow layer and granular concentration was established to characterize structures with different degrees of segregation. The corresponding relationships between the granular velocity and concentration and the segregation pattern were further analyzed to improve the theoretical models of segregation (convection–diffusion model and continuous flow model) and provide a reference for granular segregation control in the production process.