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

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

考虑颗粒旋转的颗粒动力学模拟提升管气固两相流动特性

流化床内颗粒旋转会影响颗粒相的流动特性,目前在流化床数值模拟中普遍采用的颗粒动力学模型却没有考虑颗粒的旋转效应。今运用基于颗粒动力学理论的欧拉-欧拉气固多相流模型,考虑颗粒旋转流动对颗粒碰撞能量交换和耗散的影响,提出了考虑颗粒旋转效应的颗粒动力学模型以及颗粒相守恒方程,数值模拟提升管内气体颗粒两相流动特性。计算结果表明提升管内中心区域为低浓度-高速的颗粒上升流动、壁面区域为高浓度-低速的颗粒下降流动。分析了颗粒粗糙度系数对颗粒相能量耗散、颗粒平动温度和黏度的影响。随着颗粒粗糙度系数的增加,颗粒碰撞能量耗散先逐渐增加后减小。颗粒平动温度和黏度的变化趋势是相反的,表明颗粒旋转产生摩擦将导致颗粒旋转脉动能量的改变,影响提升管内气体-颗粒两相宏观流动特性。     

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

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

圆管内含颗粒幂律流体的密相两相湍流

阐述了幂律流体控制方程的特点,研究了幂律流体的稠度系数和流动指数的变化对两相流动的影响.随着流动指数的增加,幂律流体在圆管中心附近的主流速度减小,同时颗粒相速度在圆管中心附近增大,而在管壁附近减小.随着稠度系数的增加,幂律流体和颗粒相的主流速度分布出现了与流动指数带来的影响相似的趋势.对带颗粒的幂律流体的两相流流动与液固两相流流动做了比较,幂律流体两相流的流体速度在管道中心附近的大部分区域比液固两相流的流体速度流动的速度大,而颗粒相的速度分布比较平坦.     

硝酸对AlOOH溶胶流变特性的影响

为控制氧化铝陶瓷膜制备中的涂膜工艺,以异丙醇铝为原料,硝酸为胶溶剂,制备了AlOOH溶胶。通过旋转流变仪、Zeta电位分析仪等表征手段,研究了硝酸用量对AlOOH溶胶流变特性的影响。结果表明,硝酸与异丙醇铝摩尔比R(n(HNO₃):n(Al(C₃H₇O₃)₃))为0.3~0.6时,AlOOH溶胶表现出假塑性流体特征。通过Herschel-Bulkley模型拟合流变测试数据得到的流动指数随着R的增大而减小,稠度系数随之增大;振幅扫描中储能模量(G′)与损耗模量(G″)均随着R增大而上升;频率扫描中G′与G″的交点随着R的增大而逐渐靠近低频。随着R的增大,AlOOH溶胶假塑性增强,稠度增大,弹性与黏性均上升,更易转变为凝胶。当R≥0.7时已经转化为凝胶,不再符合流体规律。Zeta电位的高低影响AlOOH溶胶的分散性,进而影响上述现象。     

鼓泡流化床内粉煤气化的数值模拟

将应用欧拉双流体模型对鼓泡气化炉内的气化过程进行研究。摒弃传统颗粒动理学理论中颗粒光滑无旋转的假设,引入颗粒的旋转运动,构建粗糙颗粒动理学理论来封闭双流体模型。基于燃烧理论建立粉煤热解、气化模型以及鼓泡床内气固之间以及气体和气体之间的传热、传质模型。采用该模型进行数值模拟计算,分析床内的气固反应过程,对比实验结果表明粗糙颗粒动理学理论适用于模拟鼓泡床气化炉内的反应。     

超细粉体对水泥净浆流变行为影响研究

采用截锥圆模法测定不同掺量超细粉体水泥复合浆液的流动度.采用ZNN-D6B型旋转黏度计研究超细玻璃粉和偏高岭土两种超细粉体对水泥净浆流变行为的影响,得到了剪切速率-剪切应力(γ-τ)曲线和剪切速率-表观粘度(γ-μa)曲线,并分别采用宾汉姆模型和赫-巴模型对γ-τ流变曲线进行拟合,得到不同掺量超细玻璃粉-水泥(GP-C)复合浆液和偏高岭土-水泥(MK-C)复合浆液的动切力、塑性粘度、稠度系数和流性指数等流变参数.结果表明:超细粉体的加入降低了复合浆液的流动度.随着掺量的增加,两种复合浆液的宾汉动切力τ0、塑性粘度η、赫-巴动切力τy均逐渐增大,MK-C复合浆液的稠度系数K和流性指数n逐渐减小,GP-C复合浆液的稠度系数K呈现增大-减小-增大的趋势,而流性指数n呈现减小-增大-减小的趋势.所有样本表观粘度μa都随着剪切速率的增大而减小,呈现剪切稀释现象.     

柴油低温粘温特性和流变性研究

考察了柴油加入低温流动改进剂前后的低温粘温特性和流变特性。结果表明,随着温度降低,柴油的表观粘度和粘流活化能增大,稠度系数增大．流动特性指数减小,柴油的流变特性越来越偏离牛顿流体,非牛顿性越来越强,柴油在低温下为假塑性非牛顿型流体。加入低温流动性改进剂T1804 后,柴油的低温表现粘度和粘流活化能显著降低,柴油的流动特性指数值增大,稠度系数值减小,与未加剂柴油相比,更接近于牛顿流体。     

考虑颗粒旋转的流化床流动和气化反应模拟

摒弃传统颗粒动力学模型中颗粒绝对光滑的假设,以粗糙颗粒为研究对象,同时考虑颗粒碰撞过程中的对心和切向分力建立了粗糙颗粒动力学模型,采用近似求解给出了相关本构关系式。结合粉煤气化反应模型模拟研究了鼓泡流化床内粉煤颗粒的流动-反应过程,获得了床内粗糙颗粒时均速度和浓度的径向分布。与光滑颗粒的计算结果相比,粗糙颗粒的脉动能量增大,床内不均匀特性进一步增强。同时得到的各气体组分的浓度分布与他人的实验结果相吻合。     

粗糙颗粒碰撞参数对鼓泡流化床内气固流动模拟的影响

结合粗糙颗粒动力学理论和双流体方法,数值模拟了碰撞参数对鼓泡流化床内稠密气固两相流动特性的影响. 结果表明,增大摩擦系数或减小法向弹性恢复系数会使床内颗粒分布更为不均,并增强床层膨胀及压力降脉动. 合理选取摩擦系数模拟得到时均气固流场分布,与实验吻合,罂粟籽颗粒的摩擦系数取0.3~0.6较合适. 法向弹性恢复系数改变不影响时均气固流场分布的基本形态,其取值敏感性不如摩擦系数;切向弹性恢复系数对鼓泡流化床动态特性及时均气固流场的影响相对较弱.     

A bubbling fluidization model using kinetic theory of rough spheres

Collisional motion of inelastic rough spheres is analyzed on the basis of the kinetic theory for flow of dense, slightly inelastic, slightly rough sphere with the consideration of gas–solid interactions. The fluctuation kinetic energy of particles is introduced to characterize the random motion of particles as a measure of the translational and rotational velocities fluctuations. The kinetic energy transport equation is proposed with the consideration of the redistribution of particle kinetic energy between the rotational and translational modes and kinetic energy dissipation by collisions. The solid pressure and viscosity are obtained in terms of the particle roughness and restitution coefficient. The partition of the random‐motion kinetic energy of inelastic rough particles between rotational and translational modes is shown to be strongly affected by the particle restitution coefficient and roughness. Hydrodynamics of gas–solid bubbling fluidized beds are numerically simulated on the basis of the kinetic theory for flow of rough spheres. Computed profiles of particles are in agreement with the experimental measurements in a bubbling fluidized bed. The effect of roughness on the distribution of energy dissipation, kinetic energy, and viscosity of particles is analyzed. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

Partial slip boundary conditions for collisional granular flows at flat frictional walls

We derive new boundary conditions (BCs) for collisional granular flows of spheres at flat frictional walls. A new theory is proposed for the solids stress tensor, translational and rotational energy dissipation rate per unit area and fluxes of translational and rotational fluctuation energy. In the theory we distinguish between sliding and sticking collisions and include particle rotation. The predictions are compared with literature results obtained from a discrete particle model evaluated at a given ratio of rotational to translational granular temperature. We find that the new theory is in better agreement with the observed stress ratios and heat fluxes than previous kinetic theory predictions. Finally, we carry out two fluid model simulations of a bubbling fluidized bed with the new BCs, and compare the simulation results with those obtained from discrete particle simulations. The comparison reveals that the new BCs are better capable of predicting solids axial velocity profiles, solids distribution near the walls and granular temperatures. © 2016 American Institute of Chemical Engineers AIChE J, 63: 1853–1871, 2017     

Multiphase flow with unequal granular temperatures

Multiphase kinetic theory has been generalized to include rotation of particles with unequal masses and diameters. Inelastic binary collisions of particles with normal and tangential restitution coefficients are considered. New expressions for number of binary collisions, viscosities and conductivities were developed. Collision integrals produced new expressions for energy dissipation involving tangential and normal restitution coefficients. Computed radial profiles of granular temperatures of 530 and 156 μm glass beads flowing in a two story riser matched the experimentally measured profiles of granular temperatures using a particle image velocity technique without rotation. The measured and the computed particle concentration profiles were nearly flat in the central portion of the riser without rotation. Computations show that rotation can alter these profiles. Sufficient particle spin can drive the particles to concentrate near the center of the tube.     

Boundary conditions for collisional granular flows of frictional and rotational particles at flat walls

Collisions between frictional particles and flat walls are determined using Coulomb friction and both tangential and normal restitution, and pseudothermal states of particles are described by both the translational and rotational granular temperatures. Then, new models for the stresses and the fluxes of fluctuation energy for the collisional granular flows at the walls are derived. These new models are tested and compared with the literature data and models. The ratio of rotational to translational granular temperatures is shown to be crucial on accurately predicting the shear stress and energy flux and is dependent on the normalized slip velocity as well as the collisional parameters. Using a theoretical but constant value for this ratio, predictions by the new models could still agree better with the literature data than those by the previous models. Finally, boundary conditions are developed to be used within the framework of kinetic theory of granular flow. © 2014 American Institute of Chemical Engineers AIChE J 60: 4065–4075, 2014     

Hydrodynamic modeling of particle rotation for segregation in bubbling gas-fluidized beds

A multi-fluid Eulerian model has been improved by incorporating particle rotation using kinetic theory for rapid granular flow of slightly frictional spheres. A simplified model was implemented without changing the current kinetic theory framework by introducing an effective coefficient of restitution to account for additional energy dissipation due to frictional collisions. Simulations without and with particle rotation were performed to study the bubble dynamics and bed expansion in a monodispersed bubbling gas-fluidized bed and the segregation phenomena in a bidispersed bubbling gas-fluidized bed. Results were compared between simulations without and with particle rotation and with corresponding experimental results. It was found that the multi-fluid model with particle rotation better captures the bubble dynamics and time-averaged bed behavior. The model predictions of segregation percentages agreed with experimental data in the fluidization regime where kinetic theory is valid to describe segregation and mixing.     

从质量流向漏斗流转变过程中的动力学分析

球床模块式高温气冷堆的堆芯是全陶瓷型包覆铀燃料制成的球形颗粒,与石墨颗粒混合堆积而成,堆芯颗粒流的流态取决于颗粒尺度的平移、旋转等动力学量,以及力链、涡旋等介尺度物理量。为了分析颗粒的平移、旋转等动力学量对颗粒流流态的影响。基于筒仓颗粒流的物理模型,首先开展了筒仓颗粒流流变过程的实验测量,并使用基于 Hertz-Mindlin和 RVD （relative velocity dependent）滚动摩擦接触模型的离散单元法 （distinct element method, DEM）,研究了锥形筒仓颗粒流流变过程中球形颗粒的动力学量。进一步,基于DEM计算结果进行分析,发现筒仓自上而下呈现出质量流向漏斗流过渡的混合流状态。在筒仓混合流的不同流型区域中,平移速度和旋转速度之间的相关性是相反的;颗粒间的相对切向运动较大的区域集中在漏斗流区域与边壁区域。了解筒仓流变过程中颗粒的动力学特征,有助于优化筒仓颗粒流动,并减少颗粒表面的磨损。     

循环流化床多组分颗粒气固两相流动模型和数值模拟

基于稠密气体分子运动论和颗粒动力学,考虑多组分颗粒中颗粒组分与颗粒组分、颗粒组分内颗粒之间的相互作用以及气体与颗粒之间的相互作用,提出多组分颗粒非等温颗粒气固两相流动模型.以颗粒压力、径向分布函数、黏度、颗粒碰撞耗散等耦合各颗粒组分间和颗粒间的相间作用.采用大涡模拟方法模拟气相湍流流动.提出了多组分颗粒的径向分布函数计算方法.对循环流化床上升管中双组分颗粒气固两相流动特性进行了数值模拟,模拟结果揭示了上升管中双组分颗粒气固两相流动的环-核流动结构,得到了平均颗粒粒径的轴向和径向分布规律,计算结果与文献中实验结果相吻合.