Inertial focusing and rotating characteristics of elliptical and rectangular particle pairs in channel flow

The lattice Boltzmann method is used to study the inertial focusing and rotating characteristics of two-dimensional elliptical particles and rectangular particles in channel flow. The results show that both elliptical particles and rectangular particles initially located on one side and two sides of channel centerline migrate first towards the equilibrium position. Then, the single-line particle train with an increasing spacing and the staggered particle train with stable spacing are formed. The axial spacing of the staggered particle pair increases with aspect ratio and Reynolds number increasing. The staggered elliptical or rectangular particle pairs form perpendicular orientation angles, which will be more obvious at larger aspect ratio and lower Reynolds number. The single-line particle trains with different shapes seldom form the perpendicular orientation angle.     

Lattice Boltzmann simulations of a dumbbell moving in a Poiseuille flow

In this paper, the lattice Boltzmann method is applied to simulate a dumbbell moving in a pressure-driven flow in a planar channel with the stress-integration method for the evaluation of hydrodynamic force acting on the cylinders. The simulation results show that the dumbbell also has the important feature of the Segr\'e--Silberberg effect like a particle in a Poiseuille flow. The dumbbell trajectories, orientations, the cylinders vertical velocities and angular velocities all reach their equilibrium values separately independent of their initial positions. It is also found that the dumbbell equilibrium positions depend on the flow Reynolds number, blockage ratio and elastic coefficient. This study is expected to be helpful to understand the dynamics of polymer solutions, polymer synthesis and reaction, etc.     

Effect of Rolling Massage on Particle Moving Behaviour in Blood Vessels

The rolling massage manipulation is a classic Chinese massage, which is expected to eliminate many diseases. Here the effect of the rolling massage on the particle moving property in the blood vessels under the rolling massage manipulation is studied by the lattice Boltzmann simulation. The simulation results show that the particle moving behaviour depends on the rolling velocity, the distance between particle position and rolling position. The average values, including particle translational velocity and angular velocity, increase as the rolling velocity increases almost linearly. The result is helpful to understand the mechanism of the massage and develop the rolling techniques.     

Describing Hydrodynamic Particle Removal from Surfaces Using the Particle Reynolds Number

The fundamental processes related to the removal of fine particles from surfaces in a hydrodynamic flow field are not adequately understood. A critical particle Reynolds number approach is proposed to assess these mechanisms for fine particles when surface roughness is small compared to particle diameter. At and above the critical particle Reynolds number, particle removal occurs, while below the critical value, particles remain attached to a surface. The system under consideration consists of glass particles adhering to a glass surface in laminar channel flow. Our results indicate rolling is the removal mechanism, which is in agreement with the literature. Theoretical results of the critical particle Reynolds number model for rolling removal are in general agreement with experimental data when particle size distribution, particle and surface roughness, and system Hamaker constant are taken into account.     

Nanoparticle electrostatic loss within corona needle charger during particle-charging process

A numerical investigation has been carried out to examine the electrostatic loss of nanoparticles in a corona needle charger. Two-dimensional flow field, electric field, particle charge, and particle trajectory were simulated to obtain the electrostatic deposition loss at different conditions. Simulation of particle trajectories shows that the number of charges per particle during the charging process depends on the particle diameter, radial position from the symmetry axis, applied voltage, Reynolds number, and axial distance along the charger. The numerical results of nanoparticle electrostatic loss agreed fairly well with available experimental data. The results reveal that the electrostatic loss of nanoparticles increases with increasing applied voltage and electrical mobility of particles; and with decreasing particle diameter and Reynolds number. A regression equation closely fitted the obtained numerical results for different conditions. The equation is useful for directly calculating the electrostatic loss of nanoparticles in the corona needle charger during particle-charging process.     

研究两相流中固粒对流体湍动特性影响的新方法

本文提出了计算两相流中固粒对流体湍动特性影响的一种新方法，得到不同情况下固粒对流体端动特性的影响。将该方法用于槽流湍流场的求解，说明该方法是可行的。     

颗粒增强湍流的新模型和气粒流动的湍流变动

用雷诺应力方程模型和极细的网格系对单个颗粒受湍流气体绕流进行了数值模拟,研究了改变颗粒直径和气体相对速度时颗粒增强气体湍流的规律.据此构造了颗粒尾涡增强气体湍流的新模型.将此子模型加入到两相流动模型中,对竖直和水平通道内气粒两相流动进行了数值模拟,和实验结果的对照表明,考虑颗粒尾涡增强气体湍流效应得到的气体湍流脉动速度的模拟结果比不考虑此效应的模拟结果好得多.     

两相混合层中颗粒运动的数值模拟

本文采用离散涡方法对平板混合层流动进行了数值模拟,得到了与实验完全定量符合的速度场。再用单向耦合方法模拟了混合层流场中颗粒的运动。分析了混合层流动中大尺度涡结构及Ｓｔｏｋｅｓ数对颗粒扩散的影响。与前人工作中所采用的每个时间步一个颗粒在固定的位置进入计算域的方法不同,本文中每个时间步有多个颗粒在入口处以随机的横向位置进入计算域。因此,在不需增加太多计算量的基础上,计算域中可以包含足够多的颗粒以获得较精确的统计结果。采用本文方法得到的颗粒速度场与实验结果定量符合得很好。     

基于拉格朗日跟踪法的波纹流道流体特性研究

采用正交实验方法考察了具有不同结构参数的三维周期波纹流道中的流体性能,并采用Webb评价方法对其进行性能评价。比较了不同波纹宽度的波纹流道的阻力因子e_f、传热因子e_Nu和能效因子η的值,结果表明三者都随Re的增大而增大,波纹宽度最小时能效因子η最大。流体在波纹流道中垂直于主流方向的横截面上产生二次流,随着Re增大,二次流增强,阻力增大,温度边界层减薄,温度等值线分布变得不均匀,传热增强。采用拉格朗日粒子跟踪技术分析了不同Re下,流体粒子在波纹流道内的运动轨迹,绘制了不同周期出口流体粒子的庞加莱截面图,结果表明流体粒子在波纹流道中被反复拉伸和折叠,增加了流体粒子的接触面积,提高混合效率,强化了传热。     

Investigation on the cavitation bubble collapse and the movement characteristics near spherical particles based on Weiss theorem

In this paper, the cavitation bubble dynamics near two spherical particles of the same size are investigated theoretically and experimentally. According to the Weiss theorem, the flow characteristics and the Kelvin impulse are obtained and supported by the sufficient experimental data. In terms of the initial bubble position, the bubble size and the distance between the two particles, the collapse morphology and the movement characteristics of the bubble are revealed in detail. The main findings include: (1) Based on a large number of experimental results, it is found that the Kelvin impulse theoretical model established in this paper can effectively predict the movement characteristics of the cavitation bubble near two particles of the same size. (2) When the initial bubble position is gradually away from the particles along the horizontal symmetry axis near two particles of the same size, the movement distance of the bubble centroid in the first period increases first and then decreases. (3) When the initial position of the bubble centroid is at the asymmetric position near the two particles, the movement direction of the bubble centroid is biased towards the particle closer to the bubble, but not towards the center of this particle.     

受碾区域内颗粒轴向流动特性的离散元模拟

为探讨受碾状态颗粒的稳定流动, 在碾辊轴与筛筒组成的受碾区域内, 建立了轴向运动的颗粒流离散元物理模型. 研究结果表明: 受碾区域内各颗粒沿轴向运动能力的差异造成了颗粒流密度不均匀; 颗粒与筛筒间的静摩擦系数影响颗粒轴向流动的形态、速率及集散程度, 受碾区域内单层颗粒的轴向均方偏差与流动时间的平方正相关, 属于“super”扩散; 整体分析受碾区域发现, 颗粒的轴向平均速度沿轴向坐标逐渐增大, 而颗粒的三轴合成平均速度沿轴向坐标逐渐降低; 受碾区域内各轴向位置处颗粒运动的剧烈程度不同, 沿轴向坐标颗粒的波动速度平方呈现先增大后降低而后又增大的趋势; 单颗粒的碰撞总能量损失能谱也表明了颗粒运动程度不同, 即轴向流动时在受碾区域的前半段碰撞剧烈, 能量损失多, 在后半段碰撞程度弱, 能量损失较少. 通过对受碾区域内颗粒流动的数值模拟分析, 明晰了颗粒在受碾条件下稳定流动特性, 有益于碾磨工业对产品品质控制及设备参数优化的研究.     

Particle transport behavior in air channel flow with multi-group Lagrangian tracking

The particle motions of dispersion and transport in air channel flow are investigated using a large eddy simulation(LES) and Lagrangian trajectory method. The mean and fluctuating velocities of the fluids and particles are obtained,and the results are in good agreement with the data in the literature. Particle clustering is observed in the near-wall and low-speed regions. To reveal the evolution process and mechanism of particle dispersion and transport in the turbulent boundary layer, a multi-group Lagrangian tracking is applied when the two-phase flow has become fully developed: the fluid fields are classified into four sub-regions based on the flow characteristics, and particles in the turbulent region are divided accordingly into four groups when the gas–particle flow is fully developed. The spatiotemporal transport of the four groups of particles is then tracked and analyzed. The detailed relationship between particle dispersion and turbulent motion is investigated and discussed.     

格子玻尔兹曼方法模拟弯流道中粒子的惯性迁移行为

本文发展了一个能够模拟微流场环境下粒子惯性迁移行为的三维耦合模型.该模型采用基于动理论的格子玻尔兹曼方法(LBM)描述流体流动,采用牛顿动力学模型描述粒子的平动和转动,采用基于LBM反弹格式的运动边界法实现流体与粒子模型的耦合.模拟了重力作用下粒子的沉降过程和Couette流条件下粒子的转动过程,通过将模拟结果与文献中的基准解进行对比定量验证了模型的可靠性.模拟了不同大小的球形粒子在环形流道中的迁移,成功复现了经典的流道截面二次流形成过程,分析了粒径大小对粒子在流道中平衡位置的影响机理.结果表明,粒子在弯流道中的平衡位置与粒径大小密切相关,小半径粒子的平衡位置靠近流道外侧而大半径粒子则靠近流道内侧.通过实验对模拟结果进行了定性验证.本模型为深入研究微流场环境下粒子的运动特性以及开发微流控粒子分选器件提供了参考依据.     

Density ratio effect on particle sedimentation in a vertical channel

A circular particle settling in a confined vertical channel with static fluid is investigated by lattice Boltzmann method, the effect of density ratio on the particle motion and flow pattern is discussed. It shows that the particle starting from an initial position off-center displays various moving modes accompanied by different flow patterns for 1.003 ≤ γ ≤ 5.0: when γ < 1.1, the particle is finally in a stable equilibrium near the centerline and settles uniformly with a symmetrical flow, but there are different modes before reaching the centerline at different density ratios, a monotonic approach appears at γ ≤ 1.005 and an transient overshoot appears at 1.005 < γ < 1.1; When γ = 1.1, a weak oscillation appears in the particle moving but the flow is almost symmetrical, it is in a critical state of the particle moving from stable equilibrium with a transient overshoot to oscillation; When 1.1 < γ ≤ 3.0, the particle oscillates regularly around the centerline with the uniform amplitude and frequency, and it turns into irregularity for γ > 3.0. Meanwhile the rotation of the particle depends on the distance from a wall; it rolls up the closer wall and stops at centerline.     

微通道内椭球颗粒惯性聚焦行为数值研究

基于浸没边界-格子Boltzmann方法，对方形截面微通道内椭球颗粒的惯性迁移与旋转动力学行为进行数值研究，发现微通道内椭球颗粒的惯性迁移存在两种主要的运动状态：①翻转状态，即椭球颗粒前进过程中长轴始终在中心对称平面内；②滚动状态，即椭球颗粒前进过程中长轴始终垂直于中心对称平面.研究表明：在低Re数（Re=10）下颗粒以两种状态随流体迁移至平衡位置；在较大Re数（50≤Re≤200）下最终均以翻转状态随流体迁移，随Re数增加，平衡位置先逼近壁面后远离壁面.通过对不同运动状态下椭球颗粒周围的微观流场进行分析，提示该微观流动在颗粒惯性聚焦行为特征中有重要影响，并从流体和颗粒的惯性角度对颗粒不同运动状态的转换机理给出解释.     

微细颗粒壁面沉积的数值研究

针对微细颗粒在壁面上的沉积特性,采用湍流雷诺应力模型结合颗粒拉格朗日轨道模型对无量纲弛豫时间τ~+为0.1～100量级的微细颗粒在壁面的沉积特性进行了研究,考查了流动方式和流速对颗粒沉积速率的影响。研究表明:在垂直向下流动和水平流动情况下,颗粒的沉积速率随着颗粒无量纲弛豫时间τ~+增加而增加,但是垂直向下流动中,当颗粒无量纲弛豫时间τ~+增大到一定值,颗粒的沉积速率基本保持定值。对于相同无量纲弛豫时间τ~+的颗粒,颗粒在壁面上的沉积速率随着摩擦速度u~*的增加而减少。     

Fouling dynamics in suspension flows

A particle suspension flowing in a channel in which fouling layers are allowed to form on the channel walls is investigated by numerical simulation. A two-dimensional phase diagram with at least four different behaviors is constructed. The fouling is modeled by attachment during collision with the deposits and by detachment caused by large enough hydrodynamic drag. For fixed total number of particles and small Reynolds numbers, the relevant parameters governing the fouling dynamics are the solid volume fraction of the suspension and the detachment drag force threshold. Below a critical curve in this 2D phase space only transient fouling takes place when the suspension is accelerated from rest by a pressure gradient. Above the fouling transition line, persistent fouling layers are formed via ballistic deposition for low and via homogeneous deposition for large solid volume fractions. Close to the fouling transition line, the flow path between the deposited layers meanders, while necking appears for increasing distance from the transition. Finally, another transition to a fully blocked flow path takes place. As determined by the estimated amount of deposited particles at saturation, both transitions seem to be discontinuous. Large fluctuations and long saturation times are typical of the dynamics of the system.     

Lattice Boltzmann Simulation of One Particle Migrating in a Pulsating Flow in Microvessel

A lattice Boltzmann model of two dimensions is used to simulate the movement of a single rigid particle suspended in a pulsating flow in micro vessel.
The particle is as big as a red blood cell, and the micro vessel is four times as wide as the diameter of the particle. It is found that Segré-Silberberg effect will not respond to the pulsation of the flow when the Reynolds number is relatively high. However, when the Reynolds number is low enough, Segré-Silberberg effect disappears. In the steady flow, different initial position leads to different equilibrium positions. In a pulsating flow, different frequencies of pulsation also cause different equilibrium positions. Particularly, when the frequency of pulsation is closed to the human heart rate, Segré-Silberberg effect presents again. The evolutions of velocity, rotation, and trajectory of the particle are investigated to find the dynamics of such abnormal phenomenon.     

Numerical Investigation of "Frog-Leap" Mechanisms of Three Particles Aligned Moving in an Inclined Channel Flow

Intrigued by our recent experimental work (H. Yamaguchi and X. D. Niu, J. Fluids Eng., 133 (2011), 041302), the present study numerically investigate the flow-structure interactions (FSI) of three rigid circular particles aligned moving in an inclined channel flow at intermediate Reynolds numbers by using a momentum-exchanged immersed boundary-lattice Boltzmann method. A "frog-leap" phenomenon observed in the experiment is successfully captured by the present simulation and flow characteristics and underlying FSI mechanisms of it are explored by examining the effects of the channel inclined angles and Reynolds numbers. It is found that the asymmetric difference of the vorticity distributions on the particle surface is the main cause of the "frog-leap" when particle moves in the boundary layer near the lower channel boundary.     

Lattice-Boltzmann Simulation of Particle Suspensions in Shear Flow

Inclusion of short-range particle–particle interactions for increased numerical stability in a lattice-Boltzmann code for particle-fluid suspensions, and handling of the particle phase for an effective implementation of the code for parallel computing, are discussed and formulated. In order to better understand the origin of the shear-thickening behavior observed in real suspensions, two simplified cases are considered with the code thus developed. A chain-like cluster of suspended particles is shown to increase the momentum transfer in a shear flow between channel walls, and thereby the effective viscosity of the suspension in comparison with random configurations of particles. A single suspended particle is also shown to increase the effective viscosity under shear flow of this simple suspension for particle Reynolds numbers above unity, due to inertial effects that change the flow configuration around the particle. These mechanisms are expected to carry over to large-scale particle-fluid suspensions.