采用格子-Boltzmann方法模拟周期性边界计算域内的颗粒流化。计算采用的流化系统Archimedes数为1 432，对应于颗粒终端Reynolds数为30。研究模拟颗粒浓度为25% ，颗粒-流体密度比为2~1 000时，密度比对流体-颗粒流动不稳定性的影响。密度比的范围对应由液固到气固的两相流动。颗粒与颗粒之间的碰撞采用弹性碰撞。研究获得颗粒平均速度、速度方差、偏度及峰度随密度比变化的规律。结合结构因子的分析，因密度比变化使颗粒-流体流动由稳定转变为不稳定的过程中颗粒速度特性变化与聚团形成的关系被确定，也确定了不稳定流动产生时所对应的密度比范围。
A lattice-Boltzmann method is used to simulate the hydrodynamic property of particles in a periodic domain. An Archimedes number of 1 432 is used for the system corresponding to the terminal Reynolds number of 30. Simulations are carried out when the average solid volume fraction is 25[WTB4]%[WTBZ], and the density ratio ranges from 2 to 1 000 (from liquid-particle density ratio to gas-particle density ratio). Investigation of the effect of density ratio on stability of the fluid-particle flow system is carried out. The collision among particles is consideved to be elastic. The variations regulations of mean particle velocity, variance, skewness and kurtosis are obtained and analyzed. By combining with the structure factor analysis, the relationships between particle velocity properties and dynamic clusters are determined at different density ratios when the stable-unstable transition occurs, and the density ratio range is also determined.
[1] Derksen J J, Sundaresan S. Direct numerical simulations of dense suspensions:wave instabilities in liquid-fluidized beds[J]. Journal of Fluid Mechanics, 2007, 587:303-336.
[2] Wen C Y, Yu Y H. Mechanics of fluidization[J]. American Institute of Chemical Engineers Symposium Series, 1966,62:100-111.
[3] Beetstra R, Van der Hoef M A, Kuipers J A M. A lattice-boltzmann simulation study of the drag coefficient of clusters of spheres[J]. Computers & Fluids, 2006, 35:966-970.
[4] Yang N, Wang W, Li J H, CFD Simulation of concurrent-up gas-solid flow in circulation fluidized beds with structure-dependent drag coefficient[J]. Chemical Engineering Journal, 2003, 96:71-80.
[5] Liu G D, Wang P, Wang S, et al. Numerical simulation of flow behavior of liquid and particles in liquid-solid risers with multi scale interfacial drag method[J]. Advanced Powder Technology, 2013, 24(2):537-548.
[6] Hill R J, Koch D L, Ladd A J C. The first effects of fluid inertia on flows in ordered and random arrays of spheres[J]. Journal of Fluid Mechanics, 2001, 449:213-241.
[7] Mitrano P P, Zenk J R, Benyahia S, et al. Kinetic-theory predictions of clustering instabilities in Granular flows:beyond the small-Knudsen-number regime[J]. Journal of Fluid Mechanics, 2014, 738 R2:1-12.
[8] William D F, Mitrano P P, Li X Q, et al. Validation of a new kinetic theory based two-fluid model for monodisperse gas-solid particulate flows[C]//Japan-US Chemical Engineering, 2015.
[9] 郭照立, 郑楚光. 格子Boltzmann方法的原理及应用[M]. 北京:科学出版社, 2009.
[10] Ladd A J C. Numerical simulations of particulate suspensions via a discretized Boltzmann equation. Part I:Theoretical foundation[J]. Journal of Fluid Mechanics, 1994, 271:285-310.
[11] Ladd A J C, Verberg R. Lattice-Boltzmann simulations of particle-fluid suspensions[J]. Journal of Statistical Physics, 2001, 104:1191-1251.
[12] Yin X L, Koch D L. Hindered settling velocity and microstructure in suspensions of solid spheres with moderate Reynolds numbers[J]. Physics of Fluids, 2007, 19:093302, 1-15.
