大颗粒气固流化床内两相流动的CFD模拟
大颗粒流化床上升管内两相流动的数值模拟
0 前 言
液固两相流动进行数值研究是两相流研 究的一个热 点. 研究颗粒团聚造成的多重尺度上的非均匀结
构 对于优化循 环流化床蒸发器 的设计 、调控 都有重要 的意义 .
对液 固两相流在垂直管束 中流动的均匀性 、稳定性 的研究多数研究者以实验 为基础 “ ,采用统计和 随机理论方法研 究床层 中颗粒 的随机运动规律 ,但对管 内颗粒 团的运动规律研究很少 , 尤其是低液速 范围.由于颗粒之 间的相互作用 ,颗粒浓度变化较大,颗粒 空隙率分布很不均匀 .Sg 和 I e e er e lrr bbg 在稀悬浮液垂直管流 中观察到 中性 悬浮的颗粒球移动并聚集成一个小的环状 区域 , 这种现象在矩形截面 的垂直流动中也 存在 ,通常的观察认 为径 向平衡位置在 = .的地方 .Ka n 认为在换热管中保持 0 6 l e r 流化状态的固体颗粒 由于和管 中的液体相互作用 , 提高 了整个管 内流体的能量 , 同时颗粒对 管壁具有刮 磨作用 . skr 等 同对 固液两相流在竖直管 中颗粒碰撞进行 了二维模拟 ,但进一步深入的研 究还未 见 A aua 报道 ,尤其是对于较大颗粒 ( =1 m .m )在垂直 向上管流中的流化特性 、速度分布 、压力及 相互影 响 0
Abs r c CFD o lo uliha efo c e lz d b hec m b n to e - l rm o l ndkie i o y f ta t m de fm tp s w anber a ie y t o i a in ofEulrEu e de n tct r o l a he muli a efo tph s w, Ba e n t s he re ,m od lngw a s d t i u ae alq d-ol fo l s d o he e t o is e i su e o sm l t iui s i l w, ditiuton e o iy d srb i sofv l ct a d c r ce itco wo p s o w a a c ae nave tc l pe Th e ulss w ha l t r xiti pea n ha a trsi ft ha ef w sc lult d i ri a l pi . er s t ho t tcuse se s nt pi nd he
基于CFD模拟的喷粉塔内气固两相流行为分析
基于CFD模拟的喷粉塔内气固两相流行为分析喷粉塔是一种用于固体颗粒物的喷涂、喷粉等工艺的设备。
在塔内,固体颗粒物与气体同时存在并发生相互作用,形成气固两相流。
为了更好地理解和优化喷粉塔的工艺条件,研究人员常常使用计算流体力学(CFD)模拟来分析塔内气固两相流行为。
CFD模拟是一种基于数学模型和计算方法的流体力学分析工具,通过离散化和数值计算等技术手段,可以模拟流体和固体在空间和时间上的运动和相互作用。
在喷粉塔的研究中,CFD模拟可以帮助研究人员了解塔内气固两相流的分布、速度、浓度等重要参数,从而优化塔内流体的运动和固体颗粒物的沉积。
首先,在进行CFD模拟前,需要建立一个合适的数学模型来描述塔内气固两相流的行为。
常见的模型包括 Euler-Euler 双流模型和 Euler-Lagrange 两流模型。
Euler-Euler 双流模型将气体和固体颗粒物视为两个相互作用的连续流体,通过求解两个连续体的质量守恒、动量守恒和能量守恒方程来模拟流动行为。
Euler-Lagrange 两流模型则将气体视为连续流体,而固体颗粒物视为离散的颗粒,通过追踪颗粒在流场中的轨迹来模拟其运动和沉积过程。
接下来,进行CFD模拟时,需要确定适当的边界条件和初始条件。
边界条件包括入口和出口条件,用于描述流体和颗粒物的初始速度、浓度以及其他相关参数。
初始条件则是指在模拟开始时,流场和颗粒物分布的初始状态。
通过合理选择和设置这些条件,可以更准确地模拟塔内气固两相流的行为。
在模拟过程中,需要选择合适的数值计算方法和离散化技术来求解流体和颗粒物的守恒方程。
常用的方法包括有限体积法、有限元法和拉格朗日法等。
这些方法可以根据具体问题的特点和要求,选择合适的离散化格式和求解算法,提高模拟结果的准确性和计算效率。
进行CFD模拟时,需要考虑塔内气固两相流的复杂性和多物理场的相互作用。
气固两相流的特点包括固体颗粒物的沉积、悬浮、弥散、聚集等过程,以及气体的扩散、对流等传输机制。
基于CFD的循环流化床旋风分离器数值模拟
基于CFD的循环流化床旋风分离器数值模拟
循环流化床旋风分离器是广泛应用于化学反应、热处理和废气处理等领域的重要设备。
该设备将固体物料以流化床形式进行循环,通过旋风分离器分离气固两相,实现了在化工、环保工程等领域的重要应用。
为了在设计和改进该设备时准确预测流动和分离的特性,数值模拟成为一种有效的方法。
基于CFD的循环流化床旋风分离器数值模拟已经成为一种快速、准确预测流动和分离过程的工具。
CFD数值模拟的过程首先需要准确建立数值模型,进行网格划分和边界条件设定。
然后利用计算机运算能力,通过数值模拟解决数学模型得到流动和分离的过程。
最后,通过数值模拟的结果,可以解决实际设计和改进的问题。
在这个过程中,需要确定正确的物理模型、适当的边界条件和网格密度以及求解方法等。
对于循环流化床旋风分离器的数值模拟,需要考虑以下几个因素。
首先,需要考虑气固两相之间的相互作用。
流化床内气固两相的运动应当由玻意耳数学模型来描述,同时应当将压降等现象考虑在内。
其次,流化床内的颗粒运动也应当考虑,因此某些模拟需采用多相流模型。
此外,旋风分离器中气相的旋转和离心效应都应当加以考虑。
基于CFD的循环流化床旋风分离器数值模拟具有许多优点。
其能够快速准确地预测流动和分离过程,可以帮助研究人员了
解设备内部的运动规律,进而优化设计、提高设备效率。
此外,该方法还可以帮助研究人员验证实验数据,并减少试验成本和时间。
因此CFD模型应用于该设备的设计和改进中,可以提
高效率,节约成本,同时优化产品质量和制造n优化生产成果。
外场作用下流化床中气固两相流动数值模拟
外场作用下流化床中气固两相流动数值模拟气固两相流化床已广泛应用于能源、化工、制药、石油等领域。
为了改善流化床的流化质量,通常采用对颗粒表面进行改性或者加入外能量场的方法,消除流化过程中出现的气固混合不均匀、扬析、沟流、颗粒损失等现象。
目前常用的外能量场有振动场、磁场、声场、电场等。
随着计算机性能的提高,离散元方法(DEM)在稠密气固两相流动数值模拟中得到广泛应用。
研究结果较好地复现了实际颗粒流化过程,预测了颗粒流动机理。
本文将对振动场、磁场和声场作为外加能量场的流化床内气固两相流动特性进行数值模拟,从宏观运动和受力分析角度研究外场对气固流动的影响。
采用Euler-Lagrange方法模拟气体和颗粒流动,颗粒碰撞采用软球模型。
同时考虑外场对颗粒受力的影响,建立不同外场作用下颗粒运动模型。
采用FORTRAN语言,自行编写计算程序。
为减小数值模拟运算量,在颗粒搜索方式上采用了定区域升序搜索,以提高运算速度。
通过上述模拟方法对外场作用下的流动现象进行复现,讨论了不同参数对气固流化特性的影响。
对于振动辅助气固流化床,考虑床体振动引起布风板所在的计算网格中心位置变化对空隙率和气体压力计算的影响,建立了振动辅助气固流化床的Euler-DEM计算模型,数值模拟研究床体竖直振动(整床振动)流化床中气体-颗粒流动过程。
研究振动幅值和振动频率对颗粒速度、浓度分布等的影响,分析振动能量从布风板传入气固两相流体的传播机理。
数值模拟发现,布风板振动导致布风板表面形成周期的低颗粒浓度区,振动空隙的出现促使床层内大气泡生成。
沿床高形成了受振动空隙影响的近布风板低颗粒浓度区域、床层中部高浓度区域和床层表面的过渡区域。
随着振动幅值和振动频率增加,平均颗粒浓度、颗粒速度、曳力径向分布都趋于均匀。
随布风板振动床层气体压力和气体压降均呈现周期振荡,由快速傅立叶变换(FFT)得到的气体压力波传播速度随振动频率增加而增大。
布风板产生的振动能量主要通过:(1)在布风板加速运动周期中布风板与颗粒之间的非弹性碰撞作用;(2)布风板减速运动周期中由气体压力波传递给床内气体-颗粒两相流体。
提升管反应器进料段气-固两相流的CFD模拟及结构优化
中图 分 类 号 :TE 2 64
G AS S LI FLo W —o D CFD I ULATI N SM o AND TRUCTURAL S oPTI I M ZATI N o
oF FEED NJ I ECTI NG ZoNE S oF Rl ER REACToR
pa tce p r il n a tce wa lc l e c d.The c lb a i n r s t ho d t tby usng t r i l a tce a d p r il— l ou d ber du e a i r to e ulss we ha i he
i i e t ie e c o .Ac o di g t heCFD i u a i e u t nsd he rs rr a t r c r n o t sm l ton r s ls,an o i ie t uc ur ft ie ptm z d s r t eo hers r r a t rwa c ie o r du e c a k ng o a a ys . Af e t uc ur o i z d o i e r a t he e c o s a qu r d t e c r c i f c t l t t r s r t e ptmi e f rs r e c or t
也大大降低 。
关
键
词: TMP ;提 升 管 反 应 器 ;气 一 固两 相 流 ;C D;催 化 剂 磨 损 F
循环流化床锅炉炉膛内气固两相流的数值模拟
循环流化床锅炉炉膛内气固两相流的数值模拟第 41卷第 3期 2020 年 5月锅炉技术BOIL ER TECHNOLO GYVol. 41, No. 3May. ,2020收稿日期 :2020 205221简介 :王建军 (19712 , 男 , 博士 , 副教授 , 主要从事流态化、多相流分离的研究。
文章编号 : CN3121508(2020 0520021206循环流化床锅炉炉膛内气固两相流的数值模拟王建军 1, 李东芳 2, 姬广勤 1, 金有海 1(1. 中国石油大学 (华东机电工程学院 , 山东东营 257061; 2. 海洋石油工程股份 ,河北塘沽 300451关键词 :循环流化床锅炉 ; 双流体模型 ; 气固两相流 ; 数值模拟摘要 :利用 CFD 软件 Fluent , ( 流的宏观流动特性进行了数值模拟。
准确性。
通过定性与定量分析 , , 核” 流动结构及颗粒轴向速度中心处向上 , , 沿轴向炉膛中下部区域及沿同时 , 操作条件对颗粒轴向速度的影响都表现为中心区域颗粒向边壁处的气固两相流动规律还有待于进一步研究。
中图分类号 : T K 227. 1文献标识码 : A0前言目前 , 对于循环流化床内的气固两相流主要集中在对循环流化床反应器[1-2]及鼓泡床 [3-4]的研究。
循环流化床锅炉炉膛内和循环流化床反应器内的气固两相流动特性有一定的差别 , 不仅体现在燃烧室的高径比 , 循环系统中采用的颗粒循环流率 , 床料的特性 , 而且循环流化床锅炉有二次风的加入 , 对循环流化床锅炉内气固两相流的研究并不多 [5-6]。
本文以欧拉双流体模型和颗粒动力学理论为基础采用 CFD 软件 Fluent 研究对循环流化床锅炉炉膛内气固两相流动特性的影响进行数值模拟。
1计算模型及数值方法1. 1几何模型及计算条件图 1为整个循环流化床锅炉循环系统几何模型及网格模型 , 模型按照工业装置 12∶ 1缩小得到。
气固湍动流化床结构曳力模型的建立及其cfd模拟
气固湍动流化床结构曳力模型的建立及其cfd模拟湍动流化床反应器因其气固接触充分、热质传递效率高、处理量大等优点,已广泛应用于化学和石油工业等诸多领域。
目前,湍动流化床的研究多集中在流动行为方面,传质行为的研究相对较少,传质行为的计算流体力学(putational fluid dynamics, cfd)模拟更是鲜有报道。
针对湍动流化床稀密两相均为半连续相的特性,本论文将其气固流动结构分为拟离散的气穴相和拟离散的聚团相,并使用cc,cd,csc,csd,cf和csf六个浓度参数描述湍动流化床的传质过程,从而根据质量守恒等原理建立基于结构的气固传质模型,最后通过cfd模拟验证气固传质模型的可靠性,并对湍动流化床的气固传质行为进行分析研究。
提出将湍动流化床的流动结构分为拟均匀的稀密两相,根据两相间质量守恒建立基于湍动流化床非均匀结构的气固传质模型的新思路。
首先,通过推导获得每一相的传质速率方程,将各相传质速率进行加和即可得总体传质速率,再结合平均传质速率定义,即可获得基于结构的湍动流化床传质系数的表达式,用于描述非均匀结构对湍动流化床气固传质的影响;其次,结合传质平衡原理、传质与反应的平衡关系等可得到一维和二维的组分输送方程,并实现传质方程所需六个浓度参数的封闭求解。
气固流动参数由基于结构的湍动流化床曳力模型求解,传质模型的数值模拟由商业软件fluent实现,以甲烷燃烧和臭氧分解实验结果对模拟结果进行校验,结果表明模拟结果与实验数据吻合较好,证明该传质模型具有较高的准确性。
通过模拟研究还发现,虽然甲烷浓度在催化剂浓度较高的位置会因化学反应消耗而降低,但组分流动以及扩散的影响使甲烷浓度与催化剂浓度之间不存在绝对的大小对应关系。
随着气速的增加,由于组分流入速率增加和停留时间变短,臭氧浓度有所增加,但因为气固传质得到强化,反应消耗的臭氧总量是增加的。
此外,稀密两相间的组分交换过程是臭氧分解反应的控制步骤。
本工作利用臭氧分解作为模型反应,通过以上的湍动流化床传质模型分析了传质模拟结果与相间传质系数求解关联式、气速的关系。
流化床内颗粒流体两相流的CFD模拟
万方数据万方数据万方数据万方数据万方数据第9期张锴等:流化床内颗粒流体两相流的CFD模拟时难以获得颗粒的真实堆积率,因此研究者们需要假设最大颗粒堆积率,如洪若瑜等[49’56巧71采用o.55,Chen等№143取o.60,Lettieri等[45]选O.62。
3.1液固体系在O.5m(高)×0.1m(宽)的二维流化床考察了液(IDl=1000kg・m一,产l一1.o×10-3Pa・s)固(佛=3000kg・m~,或一2.5×10-3m)体系内网格尺度、时间步长和收敛判据对床层固含率分布特性的影响。
结果表明:(1)从整体来看网格数目和时间步长对床层固含率分布的影响不大,但是从局部放大图可以发现,当网格数目(10×50和15×75)较少时,平衡时垂直方向上的固含率出现振荡,且10×50网格的振荡幅度大于15×75的网格,而网格数目(20×100和30×150)较多时,床层固含率趋于均匀分布特征;(2)通过对0.01、O.005、0.001、O.O005s和O.o001s时间步长的模拟表明,o.001s时间步长给出了更适宜的模拟结果;(3)收敛判据取10一、10-6和10_。
,所得模拟结果几乎完全一致,详细结果见文献[58]。
3.2气固体系首先采用摄像法考察了图2所示中心孔口为O.010m的2.Om(高)×O.3m(宽)拟二维流化床内射流形成及发展过程、射流穿透深度和射流频率。
实验以常温和常压下的空气为流化介质,GeldartB类物料的玻璃珠(佛=2550kg・m一,矾一250~300肛m,“mf一0.07m・s-1)为固体。
通过对射流气速为7.07m・s。
1的1200张图像进图2实验装置流程示意图Fig.2{khematicdiagramofexperimentalapparatus行逐帧分析,发现当时间为o.025s时射流已经形成并开始逐渐长大,到o.150s时,该射流在分布器上方脱落形成气泡,并有新的射流产生。
流化床干燥设备中物料粘流行为的数值模拟
流化床干燥设备中物料粘流行为的数值模拟流化床干燥设备是一种广泛应用于化工、食品、医药等领域的重要设备。
在流化床干燥过程中,物料粘流问题一直是制约设备性能和干燥效果的关键因素之一。
因此,对物料粘流行为进行准确的数值模拟具有重要的理论指导和工程应用价值。
各种不同性质的物料在流化床干燥过程中,由于其颗粒间作用力和干燥气体的流动引起的剪切力的存在,往往会产生粘流现象。
这种粘流行为严重影响了物料的流动性和干燥速度,甚至会导致设备淤积、堵塞、结块等问题。
因此,了解和预测物料在流化床中的粘流行为对于设计高效、稳定的流化床干燥设备至关重要。
数值模拟是研究流化床干燥设备中物料粘流行为的一种重要方法。
通过建立适当的物料颗粒力学模型和流体动力学模型,可以对物料在流化床中的运动和颗粒间相互作用进行模拟和预测。
在数值模拟中,常用的方法包括离散元法(DEM)、计算流体力学(CFD)、耦合DEM-CFD方法等。
离散元法(DEM)是一种基于颗粒间相互作用的方法,它通过考虑颗粒之间的碰撞和力学相互作用,模拟物料颗粒在流化床中的运动和变形。
DEM方法可以考虑物料颗粒间的接触力、弹性力、摩擦力等,从而准确预测颗粒运动和粘流行为。
通过DEM方法,可以对不同粒径、形状、密度等特性的物料在流化床中的流动性进行评估和优化。
计算流体力学(CFD)是研究流体流动和传热过程的一种数值模拟方法。
在流化床干燥设备中,CFD方法可以用来模拟干燥气体的流动和物料颗粒的运动。
通过建立流化床干燥设备的几何模型和物理模型,采用适当的流体动力学方程和边界条件,可以对物料颗粒的粘流行为进行准确预测。
CFD方法还可以考虑不同的操作参数对粘流行为的影响,优化流化床干燥设备的设计和运行参数。
耦合DEM-CFD方法是将DEM和CFD方法相结合的一种数值模拟方法,通过对物料颗粒和干燥气体的耦合模拟,可以更加准确地预测物料的粘流行为。
在耦合DEM-CFD方法中,DEM模拟可以用来预测物料颗粒的运动和相互作用,CFD模拟可以用来模拟干燥气体的流动和传热。
气固两相流动的数值模拟与建模
气固两相流动的数值模拟与建模气固两相流动是指在管道或设备中,同时存在气体和固体颗粒的流动现象。
这种流动在许多行业中都很常见,例如化工、能源、环境保护等领域。
通过数值模拟与建模,可以更好地理解和预测气固两相流动的特性,提高流动过程的效率和安全性。
在进行气固两相流动的数值模拟时,首先需要进行流体性质的建模。
气固两相流动中,气体和固体颗粒的物理性质和运动行为是不同的,因此需要对两相流动中的气相和固相进行单独建模。
对于气相,常用的模型有Navier-Stokes 方程和连续介质假设,通过这些模型可以描述气体在流动中的速度、压力和密度等特性。
对于固相颗粒,通常采用离散相模型,这个模型假设颗粒之间互相不作用,并体现出颗粒的运动和排列状态。
通过对气相和固相的建模,可以建立气固两相流动的数值模型。
数值模拟中最常用的方法之一是计算流体力学(CFD)方法。
CFD是通过离散化的数学方程和计算方法,对流场进行求解的一种方法。
在气固两相流动的数值模拟中,CFD方法可以用来解决气体和颗粒的速度、压力、浓度和能量等方程。
通过CFD方法,可以得到气固两相流动的速度和压力分布、颗粒浓度分布等参数,从而有效地描述了流动的特性。
除了CFD方法外,还可以采用粒子流体动力学(SPH)方法进行气固两相流动的数值模拟。
SPH方法是一种基于颗粒的数值计算方法,通过模拟颗粒的运动和相互作用,得到流场的分布和特性。
在气固两相流动中,SPH方法可以考虑颗粒之间的碰撞、沉积和湍流扩散等现象,从而更加准确地描述气固两相流动的特性。
数值模拟与建模的目的是为了更好地理解和预测气固两相流动的行为,以便优化流动过程的设计和操作。
通过数值模拟,可以得到气固两相流动中关键参数的分布规律,进而优化设备的结构和工艺参数。
例如,在化工领域中,通过数值模拟可以优化固体颗粒的输送设备,减小颗粒的堵塞和磨损程度,提高流动过程的效率和稳定性。
在能源领域中,数值模拟能够预测煤粉燃烧过程中的颗粒分布和燃烧效率,从而优化燃烧设备的设计和操作。
气固流化床反应器颗粒运动规律研究及计算流体力学模拟
气固流化床反应器颗粒运动规律研究及计算流体力学模拟结合我国“贫油、少气、多煤”的化石能源结构以及社会、经济、生态可持续发展的要求,应着力发展利用率高、环境污染小的新一代煤化工技术。
作为煤化工关键技术的费托合成煤制油和煤经甲醇制烯烃技术均大量使用了流化床反应器。
气、固流态化涉及到复杂的流动结构,其动态难以尽察;再加之应用范围扩大、装置大型化、过程强化等原因,使已有的流态化知识略显匮乏,仍需要大量深入的研究工作,以完善、优化反应器的设计和控制。
本文针对湍动流化床、循环流化床提升管、环形汽提器中的颗粒运动规律进行了大量的实验研究,并建立了计算流体力学(CFD)模型,对湍动流化床中的固体浓度、颗粒速度等进行仿真模拟。
建立了高4.8 m、内径0.15 m的湍动床大型冷模装置。
利用光纤浓度探针(PC6M)和激光多普勒测速仪(LDV)分别对床层中的固体浓度和颗粒速度进行测量。
研究表明:固体浓度呈现上稀、下浓的轴向分布和中心稀、壁面浓的径向分布。
增大表观气速或静床高,径向浓度梯度变大。
颗粒沿轴向先加速后减速,静床高越高,加速区间越长。
颗粒整体下行的环形边壁层沿床层高度逐渐变窄。
随着表观气速、静床高的增大或粒径的减小,中心高、边壁低的颗粒速度径向分布更为陡峭,环形边壁层增厚。
根据操作条件、颗粒尺寸以及测量位置的影响,建立用于预测气、固湍动床无因次边壁层厚的经验关联式,其计算值与测定值相吻合。
利用PC6M和PV6D对循环流化床提升管中不同颗粒(Glass beads I,Glass beads II,White fused alumina和SAPO-34)的固体浓度和颗粒速度进行测量。
研究了表观气速、固体循环率、颗粒物性以及轴、径向测量位置对固体浓度和颗粒运动速度的影响。
截面平均固体浓度呈现上稀、下浓的轴向分布以及中心稀、壁面浓的径向分布。
增大固体循环率、颗粒粒径、颗粒密度或减小表观气速,固体浓度升高,轴向发展愈趋缓慢,径向浓度梯度增大。
流化床干燥设备中气固两相流动的研究进展
流化床干燥设备中气固两相流动的研究进展气固两相流动在流化床干燥设备中起着重要的作用。
流化床干燥设备是一种常用的工业设备,广泛应用于化工、冶金、食品、医药等领域。
在流化床干燥设备中进行气固两相流动的研究,对于提高设备的干燥效率、降低能耗、改善产品质量具有重要意义。
气固两相流动是流化床干燥设备中的核心过程之一。
在流化床干燥设备中,气态的干燥介质通过喷嘴或旋转臂进入设备,与固体颗粒进行充分的接触和混合,使固体颗粒从内部吸收水分,实现干燥的目的。
在这个过程中,气态介质的流动特性有着重要的影响。
研究表明,气固两相流动的研究进展主要体现在以下几个方面。
首先,研究人员对气固两相流动的基本特性进行了深入的研究。
他们通过实验和数值模拟等方法,探究了气态介质在流化床中的流速分布、压力分布和温度分布等相关性质。
这些研究结果为进一步优化流化床干燥设备的操作参数提供了理论基础。
其次,研究者对气固两相流动的传质过程进行了详细的研究。
在流化床干燥设备中,气态介质通常扮演着传递热量和质量的角色。
研究人员通过实验和数值模拟等手段,研究了传热传质机制、传质速率和传质系数等关键参数,并提出了一些改进的方法来提高干燥效率。
此外,研究人员还关注了气固两相流动的流态转换过程。
在干燥过程中,气固两相流体可能处于不同的流态,如床层流态、气泡流态和涡流态等。
研究者通过实验和数值模拟,揭示了不同流态之间的转换机制和条件,并提出了一些改进的方法来优化设备的操作效率。
此外,研究者还在气固两相流动的实际应用中取得了一些进展。
例如,在干燥剂选择、干燥介质的循环利用和设备结构优化等方面,研究者提出了一些创新性的方法和理念,以提高设备的干燥效率和产品的质量。
综上所述,流化床干燥设备中气固两相流动的研究进展涉及了流动特性、传质过程、流态转换及相关应用等多个方面。
这些研究成果为改进流化床干燥设备的操作参数、提高干燥效率和优化产品质量提供了重要的理论和实践指导。
未来,需要进一步加强理论研究与实际应用的结合,以推动流化床干燥设备的进一步发展和创新。
大颗粒气固流化床内两相流动的CFD模拟
大颗粒气固流化床内两相流动的CFD模拟摘要:采用欧拉双流体模型和颗粒动力学方法,数值模拟了大颗粒流化床在不同密度、布风装置及曳力模型情况下的气固两相流动,考察了大颗粒流化床流化和流动特点,颗粒体积分率分布,床层压力瞬时变化,床层碰撞比,以及颗粒速度径向和空隙率轴向分布规律.研究结果表明,与直型布风板流化床比较,凹型布风板流化床内的气泡产生快,颗粒横向运动能力强;随着颗粒密度的增大,其在凹型布风板流化床边壁处的速度比中心位置处减小的快;比较3种曳力模型,发现其模拟的轴向空隙率分布和床层压力存在较大差异,且与床层膨胀比实验关联式相比,3种模型预测的值比实验关联式要大一些.通过研究,3个曳力模型中Gidaspow模型相对适用于大颗粒气固流化床的数值模拟.关键词:流化床;欧拉双流体模型;并行计算;大颗粒近年来,随着流态化技术的发展,大颗粒流化床在煤粉流态化燃烧和水泥熟料流态化煅烧等领域的应用也越来越广泛.由于流化床内两相流动情况复杂,使得人们对气固两相间的作用、固相应力本构方程的建立、两相湍流的认识以及多种因素的相对控制和协调的理解等变得很困难[】].实际上大多数流化床反应器都是根据经验设计的,大颗粒流化床的设计更是如此.文献[2]在研究颗粒的粒度及颗粒的表观密度等对流化特性影响后,将颗粒分成了A(30~100 tma)、B(100~600 tLm)、C(一般情况下粒度小于20 tLm)、D(600 Fm以上)4类_3].依据此分类,粒度在600肿以上的颗粒称为过粗颗粒.然而由于颗粒的表观密度与气体密度之差不同,本文所用颗粒直径为855 可能为B类(鼓泡颗粒),也有可能为D类(喷动用颗粒).其中,D类颗粒流化时极易产生大气泡或节涌,使实验难以操作,然而数值模拟可以克服这一困难,而且D类颗粒粒度在1.5 rain以下时,是完全可以流化的[3].文献[4]用粒径为3 mm的颗粒进行了模拟与实验,研究了气体进口速度和温度对床内含湿量、颗粒温度等的影响,得出模拟与实验的结果大体是一致的.文献[5]研究了表观气速、床内有无管道及布风方式对大颗粒流动的影响.模拟和试验的结果都表明,布风方式对颗粒体积分率及速度径向分布有着很大的影响,而且不论有无管道,某些布风方式都有助于气固形成环核流动结构.文献[6]通过改变颗粒粒径(从o.25 mm到1 mm)、密度、进口气速等参数后进行了模拟,结果表明:颗粒的粒径和进口气速对颗粒滑移速度的影响较大;合适的进口气速对减少能耗起着很重要的作用.本文借助CFD软件FLUENT对大颗粒气固流化床进行了模拟计算.对比并分析了不同密度颗粒、曳力模型及布风装置对流化床流动特性的影响.有些曳力模型采用皿F(用户自定义函数)实现.通过这些研究,从数值计算的角度揭示出了一些大颗粒的流化及流动特性.1 控制方程及曳力系数模型1.1 流体控制方程由于气固间没有质量交换,且升力、附加质量力等对流化床的影响很小,故气固两相流动所遵循的连续方程和动量方程可以简化成如下形式:动量方程1.2 曳力系数模型颗粒在流场中受到的作用力包括曳力、重力、浮力和其他作用力(如Basset力、Magnus力和Saff.man力等).若忽略其他力的作用,则可认为气固间作用主要为曳力作[1].Syamla1.0BriencAmstoopour~。
流化床干燥设备中气固两相流动的实验研究
流化床干燥设备中气固两相流动的实验研究随着科技进步和工业发展的推动,物料及生产过程的改进已成为追求高效和优质生产的必然要求。
在许多工业领域,干燥过程是一个重要的步骤,它对于物料的质量和工艺效率具有至关重要的影响。
干燥技术的不断发展和创新有助于提高物料的干燥速度和质量,并降低能源消耗。
其中,流化床干燥设备作为一种常见而重要的干燥设备,已经受到广泛关注。
在流化床干燥设备中,气固两相流动的研究对于提高干燥效果和优化设备运行至关重要。
气固两相流动是流化床干燥设备中一种常见的流动形式。
在流化床干燥设备中,物料被干燥的同时,气体以一定的速度通过床层,与物料颗粒发生气固两相流动。
这种气固两相流动的特点是流体具有较高的流动速度,可以实现物料的混合、传质和传热,从而提高干燥速度和质量。
实验研究是深入了解和优化气固两相流动的重要方法之一。
通过实验研究,可以获得气固两相流动的基本参数和特性,为流化床干燥设备的设计和优化提供依据。
在进行实验研究时,需要选择适当的实验装置和测量方法,以获取准确和可靠的实验数据。
在进行气固两相流动的实验研究时,首先需要确定实验的目标和研究内容。
比如,在流化床干燥设备中,可以研究气体的流速对气固两相流动的影响,物料的入口湿度对干燥效果的影响等。
然后,需要选择合适的实验装置和测量方法。
常用的实验装置包括流化床干燥设备、压力传感器、温度传感器等。
通过这些实验装置,可以获取气体和颗粒物料的压力、温度、湿度等参数,进一步分析和研究气固两相流动的特性。
在实验过程中,需要注意一些实验技巧和操作细节。
首先,需要选择合适的物料和颗粒粒度,以及合适的气体流速。
物料的选择应符合实际应用的要求,颗粒粒度的选择应考虑物料的性质和干燥速度的需求。
气体流速的选择应根据实验研究的目标和要求确定。
其次,需要进行实验前的准备工作,包括清洁实验装置、校准传感器等。
在实验过程中,需要及时记录和监测实验数据,保证实验数据的准确性和可靠性。
流化床中的CFD-DEM耦合模拟
Ii 指的是颗粒惯性矩 i 指的是颗粒角速度
Fd ,i 为流体作用于当前运动颗粒的曳力
Fp,i 为颗粒所受到流体作用力的压力梯度力
Fc,ij 为颗粒i与颗粒j之间的接触作用力;颗粒和壁面接触时,将壁面视为一个具有无限大直径的静
止颗粒,以此来计算两者之间的作用力和力矩。
三、CFD-DEM数学模型
颗粒相运动方程 Fd ,i
四、计算模型及边界条件
四、计算模型及边界条件
床体底部采用均匀速度入口,而其出口则釆用压力出口边界条件。对管壁和左右壁面,则 采用无滑移壁面边界条件来进行计算。对温度边界条件的设定,入流气体温度设定为293K, 埋管壁面温度固定为373K,其余边界条件则设定为Neumann边界条件。
在计算过程中,对固相和气相运动的时间步分别设定为1X10^-7s和1X10^-5s。对在不同流 化风速下进行的每个计算算例,均进行了15s的计算。但对床体内时均特性的统计,则在 10s后计算数据上进行。
流体作用于当前运动颗粒的曳力
gs 为曳力系数
Vp 为颗粒体积 为当前网格内流体速度矢量
vp 为颗粒的速度矢量 p 为当前计算网格内的固相浓度
f 为空隙率 p 1 f
三、CFD-DEM数学模型
颗粒相运动方程 f
颗粒切割法,即按照一定的数值算法将单个颗粒分为许多个颗粒元,根据每个颗粒元 的位置探测其所在的计算网格,得到该颗粒在该网格内的体积份额,以此来计算整个 流场网格内的空隙率。假定在计算过程中将单个颗粒分割为个部分,分别根据每个部 分的质心得到其所在计算网格。对计算网格,其内的空隙率计算为:
二、CFD-DEM耦合思路
气相运动的计算釆用欧拉方法在计算网格单元尺度上进行求解 (CFD方法);固相运动则采用离散单元法(DEM方法)在颗粒尺 度进行跟踪,即拉格朗日法。
流化床气化炉气固两相流三维数值模拟
i n i t i l a l f u i d i z a t i o n p r o c e s s ,a n a l y z i n g a x i l, a r a d i l s a o l i d h o l d u p a n d v e l o c i t y i f e l d d i s t i r b u t i o n o f p a r t i c l e a n d
流 化 床 技 术 被 广 泛 应 用 于 燃 料 的燃 烧 、 热解 、
度 矢量 图 , 结论 表明 3 D数 值 模 拟 对 于 预 测 流 体 动 量、 质量 与能 量 的传递 更为 准确 。本 文应 用 C F D软
气化等工业过程 , 流化床内气 固两相 的流动特性 已 成 为人 们 研 究 的重 点 ] 。流 化 床 反 应 器 气 固 两
摘 要: 为 了研究流 化床 内气 固两相的流动特性 , 阐述 了应用 C F D
双流化床全场气固流动的数值模拟研究
8双流化床全场气固流动的数值模拟研究双流化床全场气固流动的数值模拟研究Numerica I Simu I ation of Fu I I-fie I d Gas—soI id FI ow in DuaI FI uidized Bed罗健威金保昇胡华军黄亚继(东南大学能源热转换及其过程测控教育部重点实验室,能源与环境学院,江苏南京210096)摘要:以自行设计的双流化床装置为研究对象。
以欧拉双流体模型为基础,采用FLUENT软件对双流化床装置的全场气固流动特性进行数值模拟探究。
分析了燃烧炉的流化风量对其内部颗粒分布与运行情况的影响。
认为随着燃烧炉底部流化风量的增加,燃烧炉内整体颗粒浓度逐渐降低,靠近出口位置的颗粒上行速度增加。
关键词:双流化床;FLUENT曰数值模拟;流化风量Abstract:A se I f-designed dua I-f Iuidized bed unit is taken as the research object in this paper.Based on Eu I e r two-f I u id mode I,FLUENT software is used to simu I a te the fu I I-fie I d gas-so I id flow characteristics of a dua I f I u idized bed device.The inf I u ence of fluidization airf I o w rate on partic Ie distribution and operation of combustion furnace is ana I y zed.It is considered that the overa I I partic I e concentration in the combustion furnace decreases gradua I I y and the upward ve I o city of the partic I e s near the exit increases with the increase of fluidized air rate.Keywords:dua I f I u idized bed,FLUENT,numericaI simu I a tion,f I u idization airf I o w rate双流化床技术由于其双床的反应区域分离,既可以保证反应气氛不同,又可以通过循环床料来实现双床间能量交互的特点,被广泛应用于固体废弃物热解、CO2捕集等方面[1]遥但双流化床的种类繁多,内部气固流动特性复杂,对双流化床装置的设计与运行都造成了不小的困难。
流化床内气固两相流模拟参数敏感性分析
流化床内气固两相流模拟参数敏感性分析郑建祥;丛云龙【摘要】利用CFD-DEM耦合模型,数值模拟了竖直微型流化床内粘性颗粒的流动过程,研究分析了对数值模拟准确性具有重要影响的气相壁面条件、网格尺寸大小和杨氏模量等敏感性参数.研究结果表明,气相壁面条件对固相流场的影响可以忽略.推荐使用2.5倍~4倍粒径大小的欧拉网格,以使计算网格内可以容纳足够数量的颗粒,进而满足数值模拟统计平均的要求.临界流化速度随颗粒杨氏模量的增加而明显减小;当采用较大的杨氏模量时,流化床内出现明显的环核结构,此时更为贴近对比文献中的实验现象.【期刊名称】《东北电力大学学报》【年(卷),期】2019(039)004【总页数】8页(P41-48)【关键词】细颗粒;CFD-DEM;粘性;数值模拟;杨氏模量【作者】郑建祥;丛云龙【作者单位】东北电力大学能源与动力工程学院,吉林吉林132012;东北电力大学能源与动力工程学院,吉林吉林132012【正文语种】中文【中图分类】TQ021.1在对流化床内粘性细颗粒[1]流动的数值模拟中,细颗粒的粒径范围一般为30μm~100 μm,此粒径范围内的颗粒主要受粘性力[2~3]与重力的影响.流场中颗粒间的相互作用会导致碰撞现象的发生,在范德华力与静电力的作用下,发生碰撞的颗粒极易粘附在一起,通常来说,粘性力是一个比流场作用力更强的因素.因为基于欧拉-拉格朗日框架的计算流体动力学(Computational Fluid Dynamics,CFD)与离散单元法(Discrete Element Method,DEM)耦合模型可以对流化床内颗粒的流动碰撞过程[4~5]进行准确描述,因而该耦合模型受到愈加广泛的关注. 在利用CFD-DEM耦合模型对流化床内的气固流动行为进行数值模拟的过程中,需要设定众多对模拟结果准确性具有重要影响的模拟参数.模拟中的近似物性选择需要经验和技巧,同时气相流场和固相颗粒的参数敏感性测试尤为重要,对模拟参数的敏感性分析也一直是研究重点.Prata等[6]在2012年对流化床内颗粒的聚团现象进行了参数敏感性分析,确定了模拟中需要主要控制的参数,得到改变重要参数值的模拟结果分析.Abboud等[7]在2018年研究了工业流化床反应器中13个独立参数敏感性对模拟结果的影响.本文主要对CFD模拟中气相流场的壁面边界条件、网格尺寸参数和DEM模拟中固相颗粒的杨氏模量参数进行了模拟分析.1 数学模型在本文针对流化床内粘性细颗粒流动过程的数值模拟中,对气固两相采用欧拉-拉格朗日方法进行求解.1.1 气相计算模型气相被处理为连续相,使用时间平均的Navier Stokes(N-S)方程描述其流动行为. 质量守恒方程为+[(1-ε)ug]=0,(1)公式中:ε为固相体积分数;ug为气相速度矢量,m·s-1;t为时间,s.动量守恒方程为ug)=-pg+τg+Φd+ρg(1-ε)g,(2)公式中:ρg为气相密度,m·s-2;pg为气相压力,N;τg为气相的应力张量,N;Φd为颗粒与气体之间的相互作用项,N.1.2 固相计算模型颗粒相采用离散单元法进行求解,采用牛顿第二定律求解颗粒速度,通过对运动方程进行积分更新颗粒位置,其运动是由颗粒的重力、接触力、以范德华力为主带有静电作用的粘性力以及气固流动相间作用力控制的,颗粒的运动通过牛顿方程来求解,其质量和动量守恒方程分别为(3)(4)公式中:vi为颗粒i的移动速度,m·s-1;ωi为颗粒i的转动速度,rad·s-1;mi为颗粒i的质量,kg;Ii为颗粒i的转动惯量; fc,ij为颗粒j作用在颗粒i上的接触力,N; fc,iw为壁面与颗粒i间的接触力,N; fv,ik为颗粒k与颗粒i间的粘性力,N; fv,iw为壁面与颗粒i间的粘性力,N; fpf,i为颗粒i与气流间的相间作用力,N;Mt,ij为颗粒j对颗粒i的扭矩,N·m;Mt,iw为壁面对颗粒i 的扭矩,N·m.颗粒j对颗粒i的扭矩Mt,ij包括两个部分:一部分为切向力导致的颗粒i旋转的扭矩;另一部分为滚动摩擦扭矩,而这一部分在本文模拟中不予考虑.1.2.1 接触力模型fc,ij为颗粒j作用在颗粒i上的接触力,颗粒间的接触力表现为法向接触力和切向接触力,表达式为(5)公式中:等式右四项分别为颗粒j对颗粒i的法向力、颗粒j对颗粒i的法向阻尼力、颗粒j对颗粒i的切向力和颗粒j对颗粒i的切向阻尼力.在本文中,由于颗粒间的范德华力作用以及在接触模型中考虑颗粒间的重叠量和颗粒的表面能,因而在选择接触模型时选用Hertz-Mindlin with JKR接触模型,在JKR模型中颗粒所受的法向力与颗粒的接触半径和表面能间满足(6)公式中:ri为颗粒i的接触半径,m;R*为颗粒的有效半径,m;Y*为颗粒的有效杨氏模量系数,MPa;γi为颗粒i的表面能,J.颗粒的表面能定义为在考虑颗粒间的表面能时认为当两个颗粒在碰撞的过程中,一些能量储存在两个碰撞颗粒的弹性变形中,还有一些储存在接触区域内形成的粘性键中,另一些在周围流体中消散.所以在离同一材料的两个粘合刚性平面表面的单位面积所需要的功为两平面表面能γ的两倍,表面能可以由表面功计算得到[9]. 颗粒i与颗粒j间法向阻尼力的表达式为(7)公式中:ηn为法向阻尼系数,其大小与颗粒的恢复系数E相关,其表达式为(8)颗粒i和颗粒j间切向接触力的表达式为(9)公式中:ηt为切向阻尼系数,其大小与法向阻尼系数相近,本论文中认为ηt≡ηn;δt为颗粒在碰撞过程中的切向位移,m;μs为滑动摩擦系数;G*为颗粒的有效剪切模量系数,MPa.有效剪切模量系数的表达式为(10)在本文中不考虑切向的耗散.1.2.2 气固相间作用力颗粒与气流间的相互作用力表示为fpf,i=-Vpg+Vτg+fd,i,(11)公式中:V为颗粒的体积,m3;τg为气相偏应力张量,N·m; fd,i为颗粒i由于气体流动所受的力,N.本文采用Wen-Yu曳力模型[8],其表达式为(12)公式中:βgs为曳力系数;εg为气相体积分数;CD为阻力系数.(13)公式中:Re为雷诺系数.雷诺系数的表达式为Re=2Rsεgρgvg-vi/μg,(14)公式中:μg为气相动力粘度,Pa·s.1.2.3 范德华力计算模型对处于接触状态的弹性光滑颗粒,采用Hertz弹性模型[8].在研究中一般用黏性力表示颗粒间的黏结性.细颗粒所受粘性力主要包括静电力与范德华力,相比较而言,范德华力为起主要作用的黏性力.范德华力的基础公式Lennard-Jones关系式为(15)通过公式的相关变形,可以得到本文中颗粒i与颗粒k之间的范德华力表达式为(16)颗粒i与壁面之间的范德华力表达式为(17)2 模拟参数敏感性分析2.1 模型及数值模拟方法图1 模型尺寸示意图(单位:mm)本文数值模拟所采用的物理模型,如图1所示.XYZ轴方向的尺寸分别为1.5 mm、3.2 mm和150 mm,计算网格为正六面体网格.气相从模型底部进入流化床,并采用速度入口的边界条件;模型上部为气相的压力出口,其余设为无滑移壁面.采用结构化网格离散计算域.在模拟中,采用降速法模拟不同气相流速下的床内气固流场特性,模拟参数如表1所示.每一个气相流速稳定流化0.2 s~0.3 s,在进行CFD-DEM耦合计算前,首先在整个计算域内随机生成颗粒并自由沉降.自由沉降结束后通入气相,待气固相流场稳定后开始计时.整个求解过程采用有限体积法对气相控制方程进行离散求解,采用SIMPLE算法迭代计算从而对压力和速度进行解耦.非稳态项采用二阶隐式,动量方程采用二阶迎风格式离散.求解过程中气相流场计算时间步长为1.39×10-5s,颗粒计算时间步长为1.39×10-7s,即CFD时间步长亦为DEM时间步长的100倍.表1 CFD-DEM模拟参数表参数数值参数数值气相(空气)颗粒-颗粒静摩擦系数0.21粘度1.8 × 10-5Pa·s颗粒-颗粒静摩擦系数0.21密度1.205kg/m3颗粒-颗粒滚动摩擦系数0.01表观流速0.02~0.001m/s壁面固相密度7800粒径69μm杨氏模量7e10密度2500kg/m3泊松比0.3球形度1.0颗粒-壁面法向/切向恢复系数1剪切模量1.2e7N/m颗粒-壁面静摩擦系数0泊松比0.27颗粒-壁面滚动摩擦系数0颗粒-颗粒表面能0.000456833J/m2颗粒-壁面表面能0颗粒-颗粒法向/切向恢复系数0.972.2 壁面边界条件固体颗粒的壁面边界条件比较简单,但研究者对气相边界条件的选择却持有异议.在不同的壁面条件下,表观气速1.9 cm/s、静床高1.5 mm时,在床高1.0 mm处的时均固相Z向速度沿着Y向分布如图2所示,时均固相空隙率沿着Y向分布如图3所示(X轴数值进行归一化处理,即将网格中心Y坐标除以计算域Y向宽度).由图2、图3可知,对于微型流化床内气固两相流动模拟来说,壁面边界条件对固相颗粒的影响较小.这主要由于床内颗粒-颗粒相互作用以碰撞接触力和粘性力为主,颗粒间作用力远大于气固相间作用力导致的.同时使得壁面附近区域颗粒流场的瞬态变化不敏感.综上所述,下文的数值模拟中气相壁面边界条件选定为无滑移边界条件.2.3 网格尺寸在微型流化床气固流动特性的CFD-DEM模拟中,网格尺度划分不仅影响气相流场的精度,同时也影响固相的颗粒时均速度、时均空隙率等相关参数的计算和统计.在CFD-DEM计算中,通常只有当颗粒质心位于某个网格内部时才会考虑颗粒对网格空隙率的影响,这就造成空隙率随时间的阶梯型不连续变化.对于细网格的数值模拟来说,这样的不合理性更趋于严重.与此同时,局部空隙率通常采用体积加权平均法计算,当网格中颗粒分布表现出强烈不均匀时,密相区的局部空隙率理论上应低于稀相区的局部空隙率,但体积加权法计算的局部空隙率很难合理地反映密相区颗粒的聚集效应,并且会导致较大程度的曳力偏差.所以,CFD-DEM数值模拟方法对网格有强烈依赖性,要求网格既不能太粗以求解流场细节,又不能过小来保证容纳的颗粒数量满足统计平均的要求.基于CFD-DEM方法对网格的要求,本文按照前人常用的倍数变化关系来验证网格无关性.静床高1.5 cm、表观气速2 cm/s时沿着床高度方向的时均空隙率数值模拟结果如图4所示.由图4可知,网格尺寸小于2倍粒径时,统计平均的结果和文献[10]模拟结果偏差较大;而当网格边长为2.5 dp~4 dp时,数值模拟预测结果与文献[10]数值模拟结果吻合良好.因此,应用CFD-DEM耦合模拟时需要注意网格尺度,即使有相关经验也需要对网格进行预先评估.根据模拟数据,推荐的网格大小应是2.5 dp~4dp,这个结论与文献[10]中所使用得网格尺寸(2.6 dp)一致. 图2 表观气速1.9cm/s、静床高1.5mm时,在床高1.0mm处时均固相Z向速度沿着Y向分布图3 表观气速1.9cm/s、静床高1.5mm时,在床高1.0mm处时均固相空隙率沿着Y向分布图4 静床高1.5 cm、表观气速2 cm/s时沿着床轴向高度方向的时均空隙率分布曲线2.4 杨氏模量对固相颗粒物性的重要参数,也要进行敏感性分析研究,其中,固相颗粒的杨氏模量决定了计算的时间步长,是影响CFD-DEM耦合数值模拟计算量和顺利进行的最重要参数之一.在物体的弹性范围内,颗粒所受的应力与应变成正比,比值为颗粒的杨氏模量,大小仅取决于颗粒本身的物理性质.杨氏模量的大小说明了材料的刚性,颗粒的杨氏模量越大,也就越不容易发生形变.本文选取了3个量级的杨氏模量(73 MPa、100 MPa 、73 000 MPa)进行参数敏感性分析,其中73 000 MPa是颗粒实际杨氏模量.随着流速降低,不同流速下的典型固相速度云图如图5所示.由图5可知,当气相表观流速大于15 cm/s时,两数值模拟结果均显示流化床处于鼓泡流型状态;当气相表观流速小于5 cm/s时,两数值模拟结果显示流化床处于固定床状态;当气相表观流速小于介于5 cm/s~15 cm/s时,杨氏模量为73 MPa的数值模拟结果显示颗粒运动受到抑制,而当杨氏模量为73 000 MPa的数值模拟结果显示部分颗粒运动具有较高的运动速度.这表明较高的杨氏模量降低了从流化床转化为固定床的速度(即临界流化速度).另外,较高的杨氏模量也导致颗粒-颗粒之间法向力增大,令处于固定床状态时的床层静床高增加.图6和图7给出了表观气速1.9 cm/s、静床高1.5 mm时,不同杨氏模量下的粘性颗粒和非粘性颗粒流化床的逆流态化曲线,其中横坐标为气相表观流速、纵坐标为归一化的压降(床层压降除以流化态床层时均压降).理想状态下,随着气速的增大,归一化的压降由1逐渐变为0;归一化的压降开始小于1.0时的流速对应的流速即为临界流化速度,意味着床层开始由完全流化状态开始转变为固定床.本文设定归一化压降降低为0.98,流化床即转化为固定床.由图6可知,对非粘性颗粒流化床而言,杨氏模量对逆流态化曲线几乎无影响,气相表观流速-压降关系随着杨氏模量的变化而未表现出明显变化.但对于粘性颗粒流化床来说,杨氏模量显著影响了逆流态化曲线,临界流化速度随颗粒杨氏模量的增加而明显减小.其中,当杨氏模量分别为10 MPa、100 MPa和73 GPa时,数值模拟获得的临界流化速度如表2所示,与文献[10]结果吻合较好,这表明本文建立的数学模型能较好的模拟粘性颗粒的微型流化床内气固两相流动特性.另外,由图6和图7可知,由于颗粒间的粘性力,粘性颗粒流化床的临界流化速度明显高于非粘性颗粒流化床的临界流化速度.图5 不同杨氏模量时的典型固相速度云图图6 非粘性颗粒的逆流态化曲线图7 粘性颗粒的逆流态化曲线图8 表观气速1.9cm/s、静床高1.5mm时,在床高1.0mm的时均固相Z向速度沿Y向分布图9 表观气速1.9cm/s、静床高1.5mm时,在床高1.0mm的时均固相空隙率沿Y向分布表3 粘性颗粒临界流化速度(单位:cm/s)对比杨氏模量本文结果文献[9]结果10MPa1.271.18100MPa0.920.8273GPa0.750.64气相表观流速1.9 m/s、静床高1.5cm时,在床高1.0 mm处的时均固相速度沿着Y向分布的模拟结果如图8所示,时均空隙率沿着Y向分布的模拟结果如图9所示.由图8和图9可知,杨氏模量显著影响了沿着径向的空隙率分布和固相速度分布,且在壁面处和床层中心处表现出不同的差异.这种现象的出现主要是因为壁面区域颗粒-壁面间作用很剧烈,导致颗粒运动受到约束,频繁的碰撞就体现出了杨氏模量的影响;但在床层中心区域,数值模拟结果几乎没有受到影响.由图8可知,随着杨氏模量的增大,壁面附近颗粒纵向移动速度越小;由图2可知,随着杨氏模量的增大,壁面出颗粒的空隙率也越小,因此都可以表现出流化床内出现的环核结构越来越明显,这与文献[9]中出现的实验现象更加接近.因此,针对粘性细颗粒,应在计算机资源允许前提下采用较大的杨氏模量值进行模拟.3 结论采用前文建立的CFD-DEM强耦合模型,对颗粒-颗粒间存在粘性力的微型流化床进行了数值模拟,并进行了参数敏感性分析,获得了如下结论:(1)气相壁面条件对固相流场的影响可以忽略;(2)推荐使用的欧拉网格大小应是2.5 dp~4 dp,从而保证其可容纳足够数量的颗粒以满足数值模拟统计平均的要求;(3)杨氏模量显著影响了固相速度和空隙率分布,且临界流化速度随颗粒杨氏模量的增加而明显减小;(4)采用较大的杨氏模量,流化床内环核结构越明显,更加接近对比文献中描述的实验现象.参考文献【相关文献】[1] 白静,方红卫,何国建,等.细颗粒泥沙净冲刷和输移的大涡模拟研究[J].力学学报,2017,49(1):65-74.[2] 季佳圆.范德华力对Lennard-Jones体黏弹性的影响[J].化工学报,2018,69(8):40-46.[3] 王芳,徐怡,于恒修,等.细颗粒对气固流化床中静电行为的影响[J].化工学报,2008,59(2):341-347.[4] M.Saidi,H.B.Tabrizi,J.R.Grace,et al.Hydrodynamic investigation of gas-solid flow in rectangular spout-fluid bed using CFD-DEM modeling[J].Powder Technology,2015,284:355-364.[5] C.Duan,S.Cheng,L.Wu,et al.CFD-DEM simulation of fluid-solid flow of a tapered column separation bed[J].International Journal of Mining Science & Technology,2015,25(5):855-859.[6] A.S.Prata,A.Maudhuit,L.Boillereaux,et al.Development of a control system to anticipate agglomeration in fluidised bed coating[J].Powder Technology,2012,224:168-174.[7] A.W.Abboud,D.P.Guillen.Sensitivity study of a full-scale industrial spray-Injected fluidized bed reactor[J].Powder Technology,2018,334:36-52.[8] 郑建祥,王京阳,许帅.聚并器为颗粒聚团数值模拟[J].东北电力大学学报,2017,37(4):52-55.[9] 黄剑锋,余涛,陈江义.刚体碰撞的Hertz接触力模型比较分析[J].机械设计与制造,2017(8):28-30.[10] Peiyuan Liu,Q.Casey,LaMarche,et al.Fine-particle defluidization:interacton between cohesion,yong’s modulus and static bed height[J].Chemical Engineering Science,2016,145(1):266-278.。
新型气固环流反应器内颗粒流动的CFD模拟
新型气固环流反应器内颗粒流动的CFD模拟孟振亮;刘梦溪;李飞;王维;卢春喜【摘要】The behavior of particle flow in a new type of gas-solid air loop reactor (GSALR), which acted as a particle mixer, was numerically simulated by multi-scale computational fluid dynamics (CFD) with the structure-dependent EMMS drag model. The suitability of the drag model was verified by agreement of the simulated results and the experimental datain time averages of solid holdup and particle velocity. The particle upward velocity increased but the mean solid holdup on beds decreased when the superficial gas velocity increased in the draft tube. Several mixing regionsof particle crossflow and mixed flow in GSALR improved the efficiency of radial particle mixing. In the groove region, the distributions of solid holdup and particle velocity were more uniform, while a particle concentrating area was formed at circular overlap due to particle flow from the groove. In the region of 0≤L≤0.058 m and 0<r/R<0.3, the solid holdup was increased and particle flow was evidently enhanced along the radial direction.%采用基于结构的EMMS曳力模型,对一种新型气固环流反应器中的颗粒流动特性进行数值模拟。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
上海理工大学学报第32卷 第4期J.University of Shanghai for Science and TechnologyVol.32 No.4 2010文章编号:1007-6735(2010)04-0333-07收稿日期:2009-11-02基金项目:上海市浦江人才计划资助项目(07pj14072);上海市重点学科建设资助项目(J50501) 作者简介:晁东海(1985-),男,硕士研究生. E ma il:x yguo@大颗粒气固流化床内两相流动的CFD 模拟晁东海, 郭雪岩(上海理工大学能源与动力工程学院,上海 200093)摘要:采用欧拉双流体模型和颗粒动力学方法,数值模拟了大颗粒流化床在不同密度、布风装置及曳力模型情况下的气固两相流动,考察了大颗粒流化床流化和流动特点,颗粒体积分率分布,床层压力瞬时变化,床层碰撞比,以及颗粒速度径向和空隙率轴向分布规律.研究结果表明,与直型布风板流化床比较,凹型布风板流化床内的气泡产生快,颗粒横向运动能力强;随着颗粒密度的增大,其在凹型布风板流化床边壁处的速度比中心位置处减小的快;比较3种曳力模型,发现其模拟的轴向空隙率分布和床层压力存在较大差异,且与床层膨胀比实验关联式相比,3种模型预测的值比实验关联式要大一些.通过研究,3个曳力模型中Gidaspow 模型相对适用于大颗粒气固流化床的数值模拟.关键词:流化床;欧拉双流体模型;并行计算;大颗粒中图分类号:TQ 051.13 文献标志码:ACFD simulation on two phase flow in gas solidfluidized beds with coarse granulesCH AO Dong hai, GUO Xue yan(School of En er gy an d Pow er En gin eering ,Un iversit y of S han ghai for S cienceand Technology ,Shanghai 200093,Chin a )Abstract:Eulerian pseudo fluid model combined with the granule kinetics model,by integrating them in a CFD code(Fluent 6.3)was used to numerically simulate the gas solid flow patterns in fluidized beds of coarse granules.Different conditions including particle density,distributor types and drag models were taken into account for paramter study.The dependance of characteristics of fluidization and flow patterns,as well as the influences of phase fraction distribution,instantaneous pressure,radial particle velocity,expansion ratio and axial voidage distribution,on the parameters were thoroughly investigated.Simulation results show that two phase flow characteristics in the bed with a concave distributor is rather different from that in the bed with a flat distributor.For ex ample,bubbles will occur so oner and more particles move laterally in the concave distributor bed.It is also found that for larger solid gas density ratio,particle velocity profile near the wall becomes much flatter in the concave distributor bed.A comparison among the Syamlal O Brien,Gidaspow and Arastoopour models illustrates that the predicted axial voidage distributions and pressure drops by the three models are very different.Numerical prediction based on all the three drag models un derestimates the bed expansion ratio,comparing with the published experimental correlation.It can be concluded that numerical results based on Gidaspow drag model are of the least deviation in the上海理工大学学报2010年第32卷simulation of gas solid fluidized beds with coarse granules.Ke y words:flu idized bed;eu ler ia n pseudo fluid two fluid m odel;pa r allel com pu ta tion;coar se gr an u les近年来,随着流态化技术的发展,大颗粒流化床在煤粉流态化燃烧和水泥熟料流态化煅烧等领域的应用也越来越广泛.由于流化床内两相流动情况复杂,使得人们对气固两相间的作用、固相应力本构方程的建立、两相湍流的认识以及多种因素的相对控制和协调的理解等变得很困难[1].实际上大多数流化床反应器都是根据经验设计的,大颗粒流化床的设计更是如此.文献[2]在研究颗粒的粒度及颗粒的表观密度等对流化特性影响后,将颗粒分成了A(30~100 m)、B(100~600 m)、C(一般情况下粒度小于20 m)、D(600 m以上)4类[3].依据此分类,粒度在600 m以上的颗粒称为过粗颗粒.然而由于颗粒的表观密度与气体密度之差不同,本文所用颗粒直径为855 m可能为B类(鼓泡颗粒),也有可能为D类(喷动用颗粒).其中,D类颗粒流化时极易产生大气泡或节涌,使实验难以操作,然而数值模拟可以克服这一困难,而且D类颗粒粒度在1.5mm以下时,是完全可以流化的[3].文献[4]用粒径为3mm的颗粒进行了模拟与实验,研究了气体进口速度和温度对床内含湿量、颗粒温度等的影响,得出模拟与实验的结果大体是一致的.文献[5]研究了表观气速、床内有无管道及布风方式对大颗粒流动的影响.模拟和试验的结果都表明,布风方式对颗粒体积分率及速度径向分布有着很大的影响,而且不论有无管道,某些布风方式都有助于气固形成环核流动结构.文献[6]通过改变颗粒粒径(从0.25mm到1mm)、密度、进口气速等参数后进行了模拟,结果表明:颗粒的粒径和进口气速对颗粒滑移速度的影响较大;合适的进口气速对减少能耗起着很重要的作用.本文借助CFD软件FLUENT对大颗粒气固流化床进行了模拟计算.对比并分析了不同密度颗粒、曳力模型及布风装置对流化床流动特性的影响.有些曳力模型采用UDF(用户自定义函数)实现.通过这些研究,从数值计算的角度揭示出了一些大颗粒的流化及流动特性.1 控制方程及曳力系数模型1.1 流体控制方程由于气固间没有质量交换,且升力、附加质量力等对流化床的影响很小,故气固两相流动所遵循的连续方程和动量方程可以简化成如下形式:连续方程t(g!g)+ (g!g v g)=0(1)t(p!p)+ (p!p v p)=0(2)式中,是体积分率,v是速度,!是密度,下标p表示固相,g表示气相.动量方程t(g!g v g)+ (g!g v g v g)=-g P g+∀g+g!g g+#gp(v p-v g)(3)t(p!p v p)+ (p!p v p v p)=-p P g+∀p+p!p g+#gp(v g-v p)(4)式中,g是重力加速度,P是压力,∀为应力张量,#gp 是气体 流体相间的曳力系数.1.2 曳力系数模型颗粒在流场中受到的作用力包括曳力、重力、浮力和其他作用力(如Basset力、Mag nus力和Saff man力等).若忽略其他力的作用,则可认为气固间作用主要为曳力作用[1].Syamlal O Brien[7]、Arastoopour[8]和Gidaspow[9]等人先后对气固曳力作了大量的研究,并给出了反应相间作用强弱程度的曳力系数,如下所示:Gidaspow模型曳力系数#gp=34(1-g)g!g v g-v pd pC D0-2.7g(g!0.8)(5) #gp=150(1-g)2g d2p+74(1-g)!g v g-v pd p(g<0.8)(6)式中,C D0是单颗粒曳力系数.Syam lal O Br ien模型曳力系数#gp=3(1-g)g!g v g-v p2334第4期晁东海,等:大颗粒气固流化床内两相流动的CFD模拟式中,V r=12[a-0.06Re+(0.06Re)2+0.12Re(2b-a)+a2],Re=v g d p!g/ g,d p为颗粒直径, g为气体黏度.a= 4.14g,g∀0.85时,b=0.8 1.28g;g>0.85时,b= 2.65g.Arastoopour模型曳力系数#gp=(1-g)-2.8g!g v g-v pd p17.3Re+0.336(8)1.3 床层膨胀比床层膨胀比及床层膨胀率可按照吴占松等[10]的经验公式求得Ra=1-mf1-f(9)其中,操作条件下流化床层的空隙率f=Ar-0.21(18Re+0.36Re2)0.21阿基米德数Ar=d p3!g(!p-!g)g2g下标m f表示临界流态化.2 流化床几何结构和模拟参数设定图1(a)、(b)分别给出了两种不同布风板流化床的几何结构.两个流化床的床高和床宽分别为1.2m和0.186m(其中,凹型布风板的圆心距离床底0.116m).由于计算需要在一组离散的网格点上进行,所以模拟采用了GAM BIT2.2.30(Fluent公司,美国)进行网格生成.与表1等人相比,网格不但比其他的要小得多,而且数目也多.为了提高数值计算的速度,本次模拟使用了基于Linux操作系统的联想并行机进行求解计算.求解器采用Fluent 6.3.26(Fluent公司,美国),并选取∃-湍流模型及速度 压力耦合的#SIMPLE∃算法.初始时固定床床层的高度H为0.372m,固相体积分率为0.622,最大时达到0.65.流化床床内无任何构件,气体的密度为1.225kg/m3,黏度为1.785%10-5,颗粒直径为855 m.流化床的边界条件为:两壁面均假定成无滑移,下部为气体匀速进口,速度大小为1m/s;上部为压力出口,表压为0.依据颗粒动力论,将固相压力及黏度分别选用Lun et al和Syam lal O Brien模型,并把颗粒的弹性恢复系数假定为0.90(1代表完全弹性碰撞,0则表示完全非弹性碰撞).图1 流化床的几何结构Fig.1 Geometric structure o f the fluidized bed s表1 网格对比Tab.1 Co mp arison o f grid s作者X方向网格大小/m m,%xY方向网格大小/mm,%yZ方向网格大小/mm,%zX方向网格数Y方向网格数Z方向网格数文献[11]1010902626120文献[12]55-30200-文献[13]7.6276.2- 1072- 文献[14] 6.33 6.3318.753030160本文直型布风板 1.8612-100100-本文凹型布风板 1.866-100200-3 结果与讨论3.1 布风方式的影响图2(见下页)描述了颗粒在直型布风板流化床内不同时刻的体积分率分布.可以看出,流化床内部可以分为底部床层和上部自由空间.0s时颗粒静止地堆积在床底,当匀速气体从进口流进时,床内逐渐形成了典型的中心区颗粒体积分率小,边壁区颗粒体积分率大的环 核流动结构[15],图2中0.3s~0.6 s时间内显示了其部分变化过程.从图2中0.5s~ 1.117s可看出,气泡主要在床中心与壁面之间的区域内产生与上升.随着气泡的上升其形状和大小都在发生变化.在1.117s气泡发生破裂,颗粒被抛入335上海理工大学学报2010年第32卷自由空间.图2 直型布风板流化床瞬时颗粒(!p =906kg/m 3)体积分率的分布Fig.2 Instantaneous so lid(!p =906kg/m 3)volume fraction in the fluidized bed s with flat distributor图3显示了颗粒在不同时刻速度矢量场.其中t 表示时间,颗粒密度为2000kg /m 3,最右端t =2.5s 为床层表面情况,其他的均为床层底部.通过图形可以看出,随着时间变化,直型布风板流化床中颗粒及团聚物逐渐有了横向运动.除了t =0.2s 主要是由于重力与床壁作用外,其余都是由气泡引起的.t =0.5s 气泡在床层底部的产生;t =2.5s 气泡在床层底部的扰动及表面的抛洒作用.由此得出,气泡对流化床内颗粒及团聚物的横向运动有着重要的影响.图3 颗粒速度矢量图Fig.3 V elocity vector o f particles图4是凹型布风板流化床颗粒体积分率分布.可以看出,颗粒平行上升了仅0.08s,分布板上再次聚集了一层颗粒.与直型布风板流化床相比,由于颗粒重力及气体进口速度水平分量的影响,在底层中间位置形成了相对比较厚的颗粒层.也正是此原因,使得中间气体穿越这厚颗粒层后速度大大降低,这就加剧了其上方颗粒的下降,加速了床层在上升的过程中形成气泡.结合图2和图4,还可清楚地发现,直型布风板流化床里开始出现气泡的时间大概是在0.5s,颗粒体积分率分布开始出现明显不对称时刻大约在1.67s,而凹型分布板分别是0.33s 和1.5s,比相应的直型布风板流化床的时间都短.显然,横向对流混合起到了明显的作用.图4 凹型布风板流化床颗粒(!p =906kg/m 3)体积分率瞬时的变化Fig.4 Instantaneo us solid (!p =906kg/m 3)vo lume fractio n in the fluidized b ed s with co ncave distrib uto r流化系统中的流动结构不仅呈现局部非均匀性,而且呈现整体非均匀性.局部非均匀性表现为稀相和密相在同一点交替出现;整体非均匀性表现为系统内部不同空间位置可以出现稀相或密相两种完全不同的结构[1].图5显示了流化床内颗粒体积分率瞬时分布,其中颗粒密度为906kg/m 3,h 为床层轴向位置,H 为固定床高度,r 为床层径向位置,R 为流化床的半径.图中(a)、(c)在初始流化床高度(h/H =1)的中心位置处(r /R =0),(b)、(d)则在靠近壁面处(r /R =0.75).由图5(a)、(b)曲线可知:在直型布风板流化床里,颗粒体积分率在靠近壁面处几乎都在0.25以上波动,在床层中心位置多个时间段内小于0.1,甚至有几段几乎为0,这些说明颗粒和颗粒微团主要集中在靠近床层壁面处,气泡则位于中心位置处.相对直型布风板流化床,颗粒在凹型布风板流化床床层底部有向中心运动的能力,但通过分析固定床高处颗粒体积分率瞬时分布却得知:和直型布风板流化床内颗粒体积分率分布一样,颗粒微团靠近壁面,气泡位于中心处,如图5(c)、(d)所示.3.2 颗粒密度的影响图6显示了在初始流化阶段,不同密度颗粒在凹型布风板流化床内的体积分率分布和流线图.由图可见,随着颗粒密度的增大,其体积分率分布及运动趋势都呈现出了很大的差异.在颗粒密度!p =906kg/m 3时,颗粒体积分率在流化床内主要分布在中心及壁面处.通过云图可知,出现这种现象的原336第4期晁东海,等:大颗粒气固流化床内两相流动的CFD 模拟因是由于流化床内存在两个大气泡。