粗糙颗粒动理学及流化床内气固流动的数值模拟

基于气体分子动理学和颗粒动理学理论,考虑颗粒旋转流动对颗粒碰撞能量交换和耗散的影响,建立粗糙颗粒动理学。采用Chapman-Enskog颗粒速度分布函数,提出了颗粒相应力、热通量和颗粒碰撞能量耗散计算模型。采用欧拉-欧拉气固双相流模型,数值模拟鼓泡流化床内气体-颗粒两相流动特性。模拟结果得到了床内颗粒相速度和脉动速度分布,与Yuu等实验结果相吻合。分析不同的切向弹性恢复系数对颗粒相拟总温的变化规律,结果表明在低颗粒浓度时颗粒拟总温随切向弹性恢复系数而增加。     

考虑多组分颗粒运动各向异性的鼓泡流化床生物质气化数值模拟

考虑鼓泡流化床生物质气化过程多组分颗粒运动特点,建立多组分颗粒速度脉动二阶矩模型,结合化学反应动力学方法描述鼓泡流化床内生物质气化过程。模拟的气体组分结果与采用原始颗粒动理学模型的模拟结果进行了比较,并给出了两种粒径碳颗粒的浓度与温度瞬时分布。分析了两种碳颗粒的速度时均径向分布及速度脉动二阶矩时均径向分布,两种碳颗粒的速度分布一致,说明不同粒径碳颗粒混合充分。粒径较大的碳颗粒速度脉动二阶矩在轴向与径向上均较大,粒径的增加使得颗粒速度脉动增强。模拟统计了计算域内两种颗粒的速度脉动各向异性随颗粒浓度变化关系,各向异性随颗粒浓度的增加逐渐减弱,粒径较大的碳颗粒在计算域内的各向异性平均效果不如粒径较小的碳颗粒明显。     

大涡模拟方法模拟鼓泡床内气固两相流动

采用大涡模拟（LES）方法模拟气相湍流,颗粒动理学方法考虑颗粒相碰撞产生的动量和能量传递和耗散,采用颗粒相大涡模拟方法（LESp）模拟颗粒脉动导致的能量耗散,同时考虑介观尺度对颗粒相压力的影响,建立了气体-颗粒LES-θ-LESp双流体模型,研究鼓泡流化床内气固两相流动的特性。数值模拟与文献实测颗粒速度和实测颗粒浓度结果具有相同的变化趋势。     

颗粒旋转对鼓泡流化床内气固两相流动特性影响的数值模拟

运用考虑颗粒自旋转流动对颗粒碰撞能量交换和耗散影响的颗粒动理学方法,建立鼓泡流化床气固两相Euler-Euler双流体模型,数值模拟流化床内气体颗粒两相流动特性。分析表明,颗粒平动温度与旋转温度之比是法向和切向颗粒弹性恢复系数和摩擦系数的函数。与不考虑颗粒旋转效应计算结果相比,考虑颗粒旋转效应后床内较容易形成气泡,颗粒自旋转运动将导致床内非均匀结构更明显。并且床层平均空隙率和床层膨胀高度增加,床中心区域颗粒轴向速度提高,床内颗粒平动温度下降。考虑颗粒旋转效应后预测的颗粒轴向速度和颗粒脉动速度与文献实验结果基本吻合。考虑颗粒旋转效应后获得的气泡直径更接近于前人经验关联式。     

亚格子湍动能模型及循环流化床气固湍流特性分析

以气相大涡-颗粒相二阶矩双流体模型为框架,基于单相流亚格子湍动能推导方法,考虑固相影响推导气相亚格子湍动能方程,建立了适用于气固两相流动的气相亚格子湍动能模型;同时考虑气相亚格子湍动能与颗粒相速度脉动二阶矩之间的脉动能量传递,补充了气固相间脉动能量作用模型。模拟了循环流化床内气固两相湍流流动过程,模拟结果与实验数据吻合较好,并较未考虑湍流模型的模拟结果更接近实验值。比较了不同亚格子湍流模型对颗粒运动的影响,与Smagorinsky亚格子涡黏模型相比,亚格子湍动能模型能够更好地模拟两相流的湍流特性。分析了气体表观速度对湍流作用的影响。研究表明,随着气体表观速度的增加,气相亚格子湍动能和亚格子能量耗散逐渐增加,径向分布的非均匀性增强。     

基于气泡和聚团的结构模型及其应用

综合考虑鼓泡流化床内气泡及聚团对床内细颗粒流动的影响,建立基于气泡和聚团的结构曳力模型及结构参数模型,同时,借助计算流体软件预测细颗粒在鼓泡床中流动状态。首先,基于细颗粒在鼓泡流化床的流动状态,在介观尺度上将床内气固流动结构划分为3个子结构,即气泡相、相间相及乳化相(聚团相)。然后,综合考虑细颗粒鼓泡流化床中气泡和聚团对气固流动的影响,根据力平衡、质量守恒建立基于气泡和聚团的结构参数模型及结构曳力模型。通过对结构参数模型封闭求解,得到11个结构参数值(f_b,U_b,d_b,U_(gb),f_i,U_(gi),ε_i,f_e,U_(ge),ε_e,d_c)。对结构参数计算结果进行分析,结构参数模型能够很好反映床内流动情况及床内各结构参数之间的关系,并能有效地预测颗粒聚团直径。此外借用非均匀因子,耦合结构曳力模型及结构参数模型到欧拉双流体模型对气固在床内流动行为进行数值模拟。模拟结果表明,使用基于气泡和聚团的结构曳力模型能够较好地预测细颗粒在鼓泡流化床中的流动行为。在对模拟结果中颗粒径向浓度比较时,可以发现,相对比基于气泡模型的结构曳力模型,使用基于气泡和聚团的结构曳力模型的模拟结果与实验结果更一致。     

B类颗粒在鼓泡流化床中流动特性的数值模拟

采用欧拉双流体模型对鼓泡流化床中气-固两相流动行为进行数值模拟。模拟结果表明,采用结构曳力模型能够较好地预测B类颗粒在鼓泡流化床中的流动行为。对比初始流态化颗粒浓度图和完全流态化颗粒浓度分布图,可以发现结构曳力模型能够较好地展现鼓泡流化床中气泡的运动特性。当比较不同曳力模型下的模拟结果时,结构曳力模型比传统曳力模型能够更好地预测颗粒的径向分布。     

鼓泡床内气泡-液体两相湍流代数应力模型的数值模拟

现有的气泡 -液体两相流动的数值模拟中 ,或者不考虑湍流 ,或者仅仅考虑液体湍流 ,但是直接模拟和PIV测量结果都表明气泡由于尾迹的作用有强烈的湍流脉动 .本文首次推导和封闭了同时模拟气泡湍流脉动和液体湍流脉动的二阶矩输运方程两相湍流模型 ,并在此基础上建立了代数应力气泡 -液体两相湍流模型 .用代数应力模型模拟了二维矩形断面鼓泡床内气泡 -液体两相流动 .预报结果给出了气泡和液体两相速度场、两相Reynolds应力及湍动能分布和气泡体积分数分布 .模拟结果与PIV测量结果符合很好 ,表明了模型的合理性 .研究结果表明 ,原先静止的液体在气泡因浮力而产生的上升运动的作用下产生回流流动 ,而气泡则只有上升运动 .气泡速度始终大于液体速度 .在床内气泡湍流脉动确实始终很强烈 .液体则由于气泡的作用以及自身速度梯度产生的双重作用而发生湍流脉动 .气泡的脉动显著地大于液体的脉动 .两相湍流脉动都是各向异性的 ,而且气泡湍流脉动的各向异性比液体的更强烈     

鼓泡床取热器流动特性的数值模拟

应用计算颗粒流体力学(CPFD)方法,采用点源注射的进气方式对气固鼓泡流化床取热器内流动特性进行数值模拟。考察不同气速下床层膨胀高度、轴径向时均固含率分布、颗粒轴向速度分布及床层颗粒内循环流率的分布。模拟结果与实验数据吻合较好,表明该模型可以用于描述鼓泡流化床取热器内气固两相的流动规律。     

考虑颗粒间碰撞的稠密气/固流动二阶矩两相湍流模型

将统一二阶矩两相湍流模型和颗粒动力学理论结合,推导并封闭了稠密两相流考虑颗粒间碰撞的统一二阶矩两相湍流模型.该模型用颗粒动力学理论模拟颗粒之间的碰撞,用各向异性的统一二阶矩模型考虑气相和颗粒相的湍流脉动,并用输运方程描述气固两相湍流之间的相互作用.最后用该模型对狭窄槽道内的气粒两相流动进行了模拟,模拟所得的颗粒水平方向和垂直方向的雷诺应力和实验结果吻合良好.结果表明,考虑颗粒间碰撞之后,颗粒水平方向雷诺应力的预报得到了明显改进.     

A second-order moment method of dense gas-solid flow for bubbling fluidization

A gas-solid two-fluid model with the second-order moment method is presented to close the set of equations applied to fluidization. With the kinetic theory of granular flow, transport equations for the velocity moments are derived for the particle phase. Closure equations for the third-order moments of velocity and for the fluid-particle velocity correlation are presented. The former is based on a modified model with the contribution of the increase of the binary collision probability, and the latter uses an algebraic model proposed by Koch and Sangani [1999. Particle pressure and marginal stability limits for a homogeneous monodisperse gas-fluidized bed: kinetic theory and numerical simulations. Journal of Fluid Mechanics 400, 229-263]. Boundary conditions for the set of equations describing flow of particles proposed by Strumendo and Canu [2002. Method of moments for the dilute granular flow of inelastic spheres. Physical Review E 66, 041304/1-041304/20] are modified with the consideration of the momentum exchange by collisions between the wall and particles. Flow behavior of gas and particles is performed by means of gas-solid two-fluid model with the second-order moment model of particles in the bubbling fluidized bed. The distributions of velocity and moments of particles are predicted in the bubbling fluidized bed. Predictions are compared with experimental data measured by Muller et al. [2008. Granular temperature: comparison of magnetic resonance measurements with discrete element model simulations. Powder Technology 184, 241-253] and Yuu et al. [2000. Numerical simulation of air and particle motions in bubbling fluidized bed of small particles. Powder Technology 110, 158-168]. in the bubbling fluidized beds. The simulated second-order moment in the vertical direction is 1.1-2.5 [Muller, C.R., Holland, D.J., Sedeman, A.J., Scott, S.A., Dennis, J.S., Gladden, L.F., 2008. Granular temperature: comparison of magnetic resonance measurements with discrete element model simulations. Powder Technology 184, 241-253] and 1.1-4.0 [Yuu, S., Umekage, T., Johno, Y., 2000. Numerical simulation of air and particle motions in bubbling fluidized bed of small particles. Powder Technology 110, 158-168] times larger than that in the lateral direction because of higher velocity fluctuations for particles in the bubble fluidized bed. The bubblelike Reynolds normal stresses per unit bulk density used by Gidaspow et al. [2004. Hydrodynamics of fluidization using kinetic theory: an emerging paradigm 2002 Flour-Daniel lecture. Powder Technology 148, 123-141.] are computed from the simulated hydrodynamic velocities. The predictions are in agreement with experimental second-order moments measured by Muller et al. [2008. Granular temperature: comparison of magnetic resonance measurements with discrete element model simulations. Powder Technology 184, 241-253] and fluctuating velocity of particles measured by Yuu et al. [2000. Numerical simulation of air and particle motions in bubbling fluidized bed of small particles. Powder Technology 110, 158-168].     

Flow of gas and particles in a bubbling fluidized bed with a filtered two-fluid model

Numerical simulations of gas-particles flow in a bubble fluidized bed with two large eddy simulations of gas and solid phases are presented. For gas phase and solid phase, the sub-grid scale model for the viscosity is based on the Smagorinsky form. The sub-grid model for the particle pressure proposed by Igci et al. (2008) is modified by replacing the minimum fluidization velocity. The collisional interaction of particles is considered by the kinetic theory of granular flow. Flow behavior of gas and particles is performed by means of these two sub-grid scale models. The subgrid closure for the particle phase viscosity and pressure led to a qualitative change in the simulation results. Predictions are compared with experimental data measured by Yuu et al. (2000) and Taghipour et al. (2005) in the bubbling fluidized beds. The distributions of concentration and velocity of particles are predicted in the bubbling fluidized bed. The predicted filtered particle phase pressure increases and the filtered particle phase viscosity decreases with the increase of particle concentration. The qualitative importance of the model constant c_s of particles is demonstrated.     

水平管加压密相煤粉气力输送数值模拟

针对加压密相气力输送,对现有的颗粒静摩擦力模型进行适当修正,并将其与颗粒动理学理论相结合,建立了可以描述加压密相气力输送的气固湍流流动状况的多相流模型。该模型充分考虑了颗粒间碰撞和摩擦力作用,以及气相和颗粒团湍流脉动之间的相互作用。采用该模型对水平管内加压密相气力输送进行了三维数值模拟研究,模拟得到了气相和固相的速度、浓度和湍流强度分布,以及压降梯度的变化规律,再现了颗粒沉积层的形成和运动的动态过程。并进行了加压密相煤粉气力输送试验研究,预测的压降梯度与试验测量结果相符合。     

Numerical simulation of bubble and particles motions in a bubbling fluidized bed using direct simulation Monte-Carlo method

Flow behavior of bubbles and particles in a bubbling fluidized bed were numerically computed using Euler-Lagrange approach. Particle collision was simulated by means of the direct simulation Monte-Carlo (DSMC) method and hard-sphere model. The computed velocities and fluctuations of particles were in agreement with experimental data of Yuu et al. [S. Yuu, H. Nishikawa, T. Umekage, Numerical simulation of air and particle motions in group-B particle turbulent fluidized bed, Powder Technol. 118 (2001) 32-44]. The distributions of velocity, concentration, granular temperature and collision frequency of particles in a bubbling fluidized bed were analyzed. The wavelet multi-resolution analysis was used to investigate flow behavior of bubbles and particles. The bubble frequency of random-like bubble fluctuation was determined from the wavelet multi-resolution analysis over a time-frequency plane.     

气固流化床内射流穿透深度的CFD模拟及其实验验证

在经典的Gidaspow无黏性双流体模型中考虑离散颗粒对流体和固体动量守恒方程的影响后,建立了一个具有模拟大规模流化床内气固两相流体动力学特性潜在优势的简化数学模型。在CFX4.4商业化软件平台上通过增加用户自定义子程序考察了二维气固流化床(高2.00 m、宽0.30 m)内射流气速、喷嘴尺寸、环隙气速和静床高度对射流穿透深度的影响,并以树脂颗粒(粒径670 μm、密度1474 kg·m^-3)为研究对象在厚度为0.025 m的矩形床内进行了对比实验。结果表明,选取空隙率为0.8的等高线作为射流边界比较合适;射流穿透深度随射流气速或射流喷口尺寸的增加而增大;射流周围环隙气速由0变到最小流化速度时,射流穿透深度随环隙气速增加而增大,在最小流化速度时达到最大值,然后随环隙气速增加单调减小,当环隙气速大于2.5倍最小流化速度时,射流穿透深度减小程度变缓;在相同射流气速下射流穿透深度随着静床高度的增加而减小,静床高度对射流穿透深度的影响随着射流气速增加呈现扩大的趋势。     

A second‐order moment method applied to gas–solid risers

Second‐order moment method of particles is proposed on the basis of the kinetic theory of granular flow. Closure equations for the third‐order velocity moments are presented to account for the increase of the probability of collisions of particles on the basis of the elementary kinetic theory and order of magnitude analysis. The boundary conditions for the set of equations describing flow of particles are proposed with the consideration of the momentum exchange by collisions between the wall and the particles. The distributions of velocity, concentration and moments of particles are predicted. Simulated results are compared with experimental data measured by Tartan and Gidaspow and Bhusarapu et al. in risers, and Tsuji et al. in a vertical pipe. The effects of the closure equations for the third‐order velocity moments and the fluid‐particle velocity correlation tensor on flow behavior of particles are analyzed. © 2012 American Institute of Chemical Engineers AIChE J, 2012     

CFD simulations of bubbling beds of rough spheres

Hydrodynamics of three-dimensional gas-solid bubbling fluidized beds are numerically analyzed. The particle-particle interactions are simulated from the kinetic theory for flow of dense, slightly inelastic, slightly rough sphere proposed by Lun [1991. Kinetic theory for granular flow of dense, slightly inelastic, slightly rough sphere. Journal of Fluid Mechanics 233, 539-559] to account for rough sphere binary collisions and the frictional stress model proposed by Johnson et al. [1990. Frictional-collisional equations of motion for particulate flows and their application to chutes. Journal of Fluid Mechanics 210, 501-535] to consider the frictional contact forces between particles. The present model is evaluated by measured particle distributions and velocities of Yuu et al. [2001. Numerical simulation of air and particle motions in group-B particle turbulent fluidized bed. Powder Technology 118, 32-44] and experimental bed expansion of Taghipour et al. [2005. Experimental and computational study of gas-solid fluidized bed hydrodynamics. Chemical Engineering Science 60, 6857-6867]. Our computed results indicated that the present model gives better agreement with experimental data than the results from original kinetic theory for frictionless slightly inelastic sphere of Ding and Gidaspow [1990. A bubbling fluidization model using kinetic theory of granular flow. A.I.Ch.E. Journal 36, 523-538] with and without solid friction stress model.     

Hydrodynamic simulations of gas–solid spouted bed with a draft tube

Flow behavior of particles in a two-dimensional spouted bed with a draft tube is studied using a continuous kinetic-friction stresses model. The kinetic stress of particles is predicted from kinetic theory of granular flow, while the friction stress is computed from a model proposed by Johnson et al. (1990). A stitching function is used to smooth from the rapid shearing viscous regime to the slow shearing plastic regime. The distributions of concentration and velocities of particles are predicted in the spouted bed with a draft tube. Simulated results compare with the vertical velocity of particles (Zhao et al., 2008) measured and in the spout bed with draft plates and solid circulation rate (Ishikura et al., 2003) measured in the spouted bed with a draft tube. The impact of the friction stress of particles on the spout, annulus, fountain and entrancement regions is analyzed in gas–solid spouted bed with a draft tube. Numerical results show that the gas flow rate through the annulus increases with the increase of the entrainment zone. The solids circulation rate decreases with the decrease of inlet gas velocity and the height of the entrainment zone. The effect of spouting gas velocity on distributions of concentration, velocity and particle circulation is discussed.     

Numerical predictions of flow behavior and cluster size of particles in riser with particle rotation model and cluster-based approach

Flow behavior of particles in a circulating fluidized bed (CFB) riser is numerically simulated using a two-fluid model incorporating with the kinetic theory for particle rotation and friction stress models. The particle rotations resulting from slightly friction particle-particle collisions was considered by introducing an effective coefficient of restitution based on the kinetic theory for granular flow derived by Jenkins and Zhang [2002. Kinetic theory for identical, frictional, nearly elastic spheres. Physics of Fluids 14, 1228-1235]. The normal friction stress model proposed by Johnson et al. [1990. Frictional-collisional equations of motion for particles flows and their application to chutes. Journal of Fluid Mechanics 210, 501-535] and a modified frictional shear viscosity model proposed by Syamlal et al. [1993. MFIX Documentation and Theory Guide, DOE/METC94/1004, NTIS/DE94000087] were used as the particle frictional stress model. The drag force between gas and particle phases was modified with cluster-based approach (CBA). The flow behavior of particles and the cluster size in a riser of Wei et al. [1998. Profiles of particle velocity and solids fraction in a high-density riser. Powder Technology 100, 183-189] and Issangya et al. [2000. Further measurements of flow dynamics in a high-density circulating fluidized bed riser. Powder Technology 111, 104-113] experiments are predicted. Effects of the rotation and friction stress models on the computed results are analyzed. It is concluded that particle rotations reduce the cluster size and alter the particle flows and distributions through more particle fluctuation energy dissipations. Effects of frictional stress on flow behavior and cluster size are not significant because the particle phase in the CFB riser is not dense enough to take into account for the particle-particle contact interactions.