壁面约束对柱状粒子在牛顿流体中沉降影响的研究

格子Boltzmann方法(LBM)直接数值模拟了在有壁面约束的流场中柱状粒子的沉降。结果说明粒子在壁面的影响下，会离开管壁向管中心移动，移动速度先是增大，然后随着与管壁距离的增加，逐渐减小，粒子中心最后并不一定停留在管中心。随着粒子与壁面距离的增加，粒子的沉降速度将增大，可见壁面对粒子沉降有阻碍作用，平行于粒子轴线的壁面对粒子沉降的影响比垂直于粒子轴线壁面的影响大。     

密度不同的颗粒在流体中的沉降特性

采用格子Boltzmann方法(LBM)对雷诺数范围为5≤Re≤12的双颗粒沉降进行了直接数值模拟研究,主要关注颗粒之间的密度差异k对其周期性振动的影响。根据雷诺数颗粒沉降特性可以分成三个阶段,即当Re较小时,颗粒的振动幅度随k的增大而减小,当Re较大时则正好相反,介于两者之间存在一个临界雷诺数,在此雷诺数附近,颗粒的沉降同时具有以上两种特征。同时还研究了两个颗粒的稳定沉降结构以及重颗粒摆脱轻颗粒的条件。     

双颗粒沉降的动力学分析

采用格子Boltzmann-虚拟区域方法对雷诺数范围为50≤Re≤200的双颗粒自由沉降进行了直接数值模拟。首先,研究发现双颗粒最终沉降的位置分别在通道的1/4和3/4位置附近,且与颗粒的初始间距、雷诺数以及通道宽度无关。其次,重点研究了颗粒沉降过程中所受的侧向力(与沉降方向垂直),首次揭示了侧向力的振动频率与雷诺数呈二次关系,且单颗粒的结果始终小于双颗粒的结果;研究还发现侧向力的振动频率与通道宽度近似呈幂律函数关系,且幂律指数与雷诺数有关,雷诺数越大,幂律指数的绝对值越小。最后还研究了雷诺数及通道宽度对侧向力振动振幅的影响。     

Numerical simulation of particle sedimentation in 3D rectangular channel

The 3D lattice Boltzmann method is used to simulate particle sedimentation in a rectangular channel. The results of single particle sedimentation indicate that the last position of the particle is along the center line of the channel regardless of the initial position, the particle diameter, and the particle Reynolds number. The wall effect on the terminal velocity is in good agreement with experimental results quantitatively. The drafting, kissing, and tumbling (DKT) process is reproduced and analyzed by simulating two-particle cluster sedimentation. The effects of the diameter ratio, initial position, and wall on the DKT process are investigated. When the two particles have equal diameter sediment in the rectangular channel, a periodical DKT process and the spiraling trajectory are found. The last equilibrium configuration is obtained from the simulation results. The interesting regular sedimentation phenomena are found when 49 particles fall down under gravity.     

The approximate solution of the self-similar problem for radial flow of non-newtonian fluids through porous media

In this paper, we study the approximate solution of the self-simikar problem for radial flow of non-Newtonian fluids through porous media. Assuming that the fluids obey the exponential function law, we obtain an exact solution for the exponent n=0 and compare it with the approximate solution in ref. [1]. For n>1 and n<1, we obtain respectively approximate solutions. Some exampls are presented.     

INTERACTIONS BETWEEN TWO SEDIMENTING PARTICLES WITH DIFFERENT SIZES

IntroductionParticulateflowsarecommonlyencounteredinvariousnaturalphenomenaandindustrialprocesses,suchasinthechemicalandpetroleumindustries,thesedimentationofparticlesisawayofseparatingparticlesfromfluid ,aswellasawayofseparatingparticleswithdifferentsettlingspeedsfromeachother[1].Animportantaspectoftheseflowsisparticlepairinteraction ,whichisfundamentalmechanismsthatenterstronglyintoallpracticalapplicationsofparticulateflows.Problemofthesedimentationoftwoequalparticleshasbeenpreviouslystudied…     

New formula for drag coefficient of cylindrical particles

The drag force on a cylindrical particle is calculated using lattice Boltzmann method.The results show that the drag coefficient of a particle with different orientation angles decreases with increasing Reynolds number.When the principal axis of the particle is parallel to flow,the drag coefficient is much larger than that of others and decreases fastest with increasing Reynolds number,which becomes more obvious with increasing particle aspect ratio.When the principal axis of the particle is inclined to flo...     

A direct‐forcing pressure‐based lattice Boltzmann method for solving fluid–particle interaction problems

A direct‐forcing immersed boundary‐lattice Boltzmann method (IB–LBM) is developed to simulate fluid–particle interaction problems. This method uses the pressure‐based LBM to solve the incompressible flow field and the immersed boundary method to handle the fluid–particle interactions. The pressure‐based LBM uses the pressure distribution functions instead of the density distribution functions as the independent dynamic variables. The main idea is to explicitly eliminate the compressible effect due to the density fluctuation. In the IB method, a direct‐forcing method is introduced to capture the particle motion. It directly computes an IB force density at each lattice grid from the differences between the pressure distribution functions obtained by the LBM and the equilibrium pressure distribution functions computed from the particle velocity. By applying this direct‐forcing method, the IB–LBM becomes a purely LBM version. Also, by applying the Gauss theorem, the formulas for computing the force and the torque acting on the particle from the flows are derived from the volume integrals over the particle volume instead of from the surface integrals over the particle surface. The order of accuracy of the IB–LBM is demonstrated on the errors of velocity field, wall stress, and gradients of velocity and pressure. As a demonstration of the efficiency and capabilities of the new method, sedimentation of a large number of spherical particles in an enclosure is simulated. Copyright © 2010 John Wiley & Sons, Ltd.     

等离子熔射粉末颗粒飞行过程格子Boltzmann法仿真

为考察等离子熔射过程中粉末的飞行过程,本文在已开发的正六边形7-b it格子Bo ltzm ann(LB)方法等离子射流的温度场和速度场的计算模型基础上,采用单个颗粒加速方程,建立了一个随机算法,实现了对粉末颗粒在射流场中运动过程的仿真;计算结果通过动画演示了粉末飞行的全过程,表明初始位置越靠近射流场出口中心的粉末颗粒加速越充分,并且在射流场一定的情况下,减小粉末颗粒直径可以提高粉末速度,但会降低粉末利用率。     

Integrating thermal-concentration smoothed profile with lattice Boltzmann methods for simulating sedimentation of nonisothermal circular particles

A thermal-concentration smoothed profile-lattice Boltzmann method is proposed to study the effect of the concentration field on the dynamic behavior of nonisothermal cylindrical particles during the sedimentation process. The velocity, temperature, and concentration equations are solved using the lattice Boltzmann method. Moreover, the smoothed profile method is employed to enforce the nonslip boundary condition as well as constant temperature and constant concentration boundary conditions at the particles surfaces. Moreover, the Boussinesq approximation is used to couple the velocities, temperatures, and concentrations fields. The proposed combined method is validated by comparing the present numerical results with those found in the literature, showing good consistency. Then, the effect of the concentration buoyancy on the behavior of nonisothermal particles is discussed. In addition, the effect of Prandtl, Schmidt, and thermal Grashof numbers on the settling process is investigated. The results show that, by adding the effect of concentration, the maximum settling velocity of hot particles is reduced more relative to the cold ones; accordingly, the cold particles are settled faster than the hot ones. Finally, the sedimentation of two particles in a container at high thermal Grashof is investigated. It is shown that, at high thermal Grashof, there is an intense competition between the buoyancy force and gravity for the hot particles. The buoyancy flow generated leads to the reversal of the drafting-kissing-tumbling motion of the hot particles, making the particles move upward.     

Inertial migration of sedimenting particles in a suspension flow through a Hele-Shaw cell

Within the framework of the model of two interpenetrating continua, a horizontal laminar dilute-suspension flow in a vertical Hele-Shaw cell is investigated. Using the method of matched asymptotic expansions, an asymptotic model of the transverse migration of sedimenting particles is constructed. The particle migration in the horizontal section of the cell is caused by an inertial lateral force induced by the particle sedimentation and the shear flow of the carrier phase. A characteristic longitudinal length scale is determined, on which the particles migrate across the slot through a distance of the order of the slot half-width. The evolution of the particle number concentration and velocity fields along the channel is studied using the full Lagrangian method. Depending on the particle inertia parameter, different particle migration regimes (with and without crossing of the channel central plane by the particles) are detected. A critical value of the particle inertia parameter corresponding to the change in migration regime is found analytically. The possibility of intersection of the particle trajectories and the formation of singularities in the particle number concentration is demonstrated.     

Numerical simulation of polygonal particles moving in incompressible viscous fluids

A two-dimensional coupled lattice Boltzmann immersed boundary discrete element method is introduced for the simulation of polygonal particles moving in incompressible viscous fluids. A collision model of polygonal particles is used in the discrete element method. Instead of a collision model of circular particles, the collision model used in our method can deal with particles of more complex shape and efficiently simulate the effects of shape on particle–particle and particle–wall interactions. For two particles falling under gravity, because of the edges and corners, different collision patterns for circular and polygonal particles are found in our simulations. The complex vortexes generated near the corners of polygonal particles affect the flow field and lead to a difference in particle motions between circular and polygonal particles. For multiple particles falling under gravity, the polygonal particles easily become stuck owing to their corners and edges, while circular particles slip along contact areas. The present method provides an efficient approach for understanding the effects of particle shape on the dynamics of non-circular particles in fluids.     

Lattice Boltzmann study of the interaction between a single solid particle and a thin liquid film

The isothermal single-component multi-phase lattice Boltzmann method(LBM) combined with the particle motion model is used to simulate the detailed process of liquid film rupture induced by a single spherical particle.The entire process of the liquid film rupture can be divided into two stages.In Stage 1,the particle contacts with the liquid film and moves into it due to the interfacial force and finally penetrates the liquid film.Then in Stage 2,the upper and lower liquid surfaces of the thin fi...     

悬浮颗粒运动的格子Boltzmann数值模拟

将固体颗粒的牛顿力学和格子Boltzmann方法相结合,研究不规则形状悬浮颗粒在流场中的运动。通过受力分析,精确求得其所受合力、合力矩、合力作用中心等。提出了跟随颗粒运动的动网格计算域技术和模拟悬浮颗粒转动运动的局部数组方法及Euler-Lagrange两套坐标技术。通过对椭圆颗粒运动的数值模拟和对照他人对矩形颗粒的研究,分析了其复杂运动规律,并提供了合理的物理解释。结果表明:运用格子Boltzmann方法和上述特殊技术可以得到与有限元方法相同的模拟精度,且具有计算速度快、对复杂形状边界处理方便灵活、程序简单及特别适合大规模并行计算等优点。     

Dynamics of two elliptical particles sedimentation in a vertical channel: chaotic state

The LB-DF/FD method derived from the Lattice Boltzmann Method and direct forcing/fictitious domain method is used to numerically investigate the dynamics and interaction of two elliptical particles settling in an infinitely long channel. One particle (EP0) is initially kept horizontal (major axis perpendicular to sedimentation) for all simulations while the other's (EP1) orientation is varied. It is found that if EP1strays away from horizontality, the particles undergo transitions from a steady state to reach a chaotic state. Furthermore, there are two distinct chaotic states for the particle motion when EP1 orientation is varied, in which a turning point is observed to distinguish the two states.     

Free vibration of layered cylindrical shells filled with fluid

The vibration of the layered cylindrical shells filled with a quiescent, incompressible, and inviscid fluid is analyzed. The governing equations of the cylindrical shells are derived by Love’s approximation. The solutions of the displacement functions are assumed in a separable form to obtain a system of coupled differential equations in terms of the displacement functions. The displacement functions are approximated by Bickley-type splines. A generalized eigenvalue problem is obtained and solved numerically for the frequency parameter and an associated eigenvector of the spline coefficients. Two layered shells with three different types of materials under clamped-clamped (C-C) and simply supported (S-S) boundary conditions are considered. The variations of the frequency parameter with respect to the relative layer thickness, the length-to-radius ratio, the length-to-thickness ratio, and the circumferential node number are analyzed.     

Lattice‐Boltzmann and finite element simulations of fluid flow in a SMRX Static Mixer Reactor

A detailed comparison between the finite element method (FEM) and the lattice‐Boltzmann method (LBM) is presented. As a realistic test case, three‐dimensional fluid flow simulations in an SMRX static mixer were performed. The SMRX static mixer is a piece of equipment with excellent mixing performance and it is used as a highly efficient chemical reactor for viscous systems like polymers. The complex geometry of this mixer makes such three‐dimensional simulations non‐trivial. An excellent agreement between the results of the two simulation methods was found. Furthermore, the numerical results for the pressure drop as a function of the flow rate were close to experimental measurements. Results show that the relatively simple LBM is a good alternative to traditional methods. Copyright © 1999 John Wiley & Sons, Ltd.     

Acoustic radiation force of a solid elastic sphere immersed in a cylindrical cavity filled with ideal fluid

An expression for the acoustic radiation force function on a solid elastic spherical particle placed in an infinite rigid cylindrical cavity filled with an ideal fluid is deduced when the incident wave is a plane progressive wave propagated along the cylindrical axis. The acoustic radiation force of the spherical particle with different materials was computed to validate the theory. The simulation results demonstrate that the acoustic radiation force changes demonstrably because of the influence of the reflective acoustic wave from the cylindrical cavity. The sharp resonance peaks, which result from the resonance of the fluid-filled cylindrical cavity, appear at the same positions in the acoustic radiation force curve for the spherical particle with different radii and materials. Relative radius, which is the ratio of the sphere radius and the cylindrical cavity radius, has more influence on acoustic radiation force. Moreover, the negative radiation forces, which are opposite to the progressive directions of the plane wave, are observed at certain frequencies.