多喷嘴喷动流化床内气固两相...传热过程CFD-DEM模拟_刘雪峰.pdf
《多喷嘴喷动流化床内气固两相...传热过程CFD-DEM模拟_刘雪峰.pdf》由会员分享,可在线阅读,更多相关《多喷嘴喷动流化床内气固两相...传热过程CFD-DEM模拟_刘雪峰.pdf(9页珍藏版)》请在咨信网上搜索。
1、第 37 卷第 3 期 高 校 化 学 工 程 学 报 No.3 Vol.37 2023 年 6 月 Journal of Chemical Engineering of Chinese Universities June 2023 文章编号:1003-9015(2023)03-0383-09 多喷嘴喷动流化床内气固两相流动与传热过程 CFD-DEM 模拟 刘雪峰1,白谨豪1,吴 峰1,马晓迅1,2,危仁波1(1.西北大学 化工学院 陕北能源先进化工利用技术教育部工程研究中心,陕西 西安 710069;2.陕西省洁净煤转化工程技术中心,陕西 西安 710069)摘 要:为研究喷动床内颗粒的流体
2、力学与传热特性,采用计算流体力学-离散单元(CFD-DEM)耦合方法对多喷嘴喷动流化床(IMJSFB)内的颗粒进行数值模拟,将模拟结果与常规喷动床(CSB)对比。结果表明:IMJSFB 内颗粒混合与流动特性明显优于 CSB,表明侧开孔的存在改善了颗粒流动死区现象。IMJSFB 内的颗粒总势能和颗粒总平移动能均低于 CSB,侧开孔对主喷嘴气体实现了分流,削弱了喷射区的颗粒运动。当颗粒的传热系数低、热容小、进气温度较低时 IMJSFB 的优势更显著。关键词:多喷嘴喷动流化床;计算流体力学;颗粒混合特性;流动特性 中图分类号:TQ021.3 文献标志码:A DOI:10.3969/j.issn.10
3、03-9015.2023.03.006 CFD-DEM simulation of gas-solid two-phase flow and heat transfer in multi-jet spout-fluidized bed LIU Xuefeng1,BAI Jinhao1,WU Feng1,MA Xiaoxun1,2,WEI Renbo1(1.Northern Shaanxi Energy Advanced Technology with Chemical Engineering Research Center of the Ministry of Education,School
4、 of Chemical Engineering,Northwest University,Xian 710069,China;2.Shaanxi Research Center of Engineering Technology for Clean Coal Conversion,Xian 710069,China)Abstract:The characteristics of particle flow and heat transfer in multi-jet spout-fluidized bed(IMJSFB)were studied by the computational fl
5、uid dynamics-discrete element method(CFD-DEM),and the simulation results were compared with conventional spouted bed(CSB).The results showed that the mixing and flow characteristics of particles in IMJSFB were better than those in CSB,indicating that the presence of side holes improves the particle
6、flow in dead zone.The total potential energy and total translational kinetic energy of particles in the IMJSFB are lower than those in CSB,and the existence of side holes diverts the gas from the main jet and weakens the movement of particles in spout region.The advantages of IMJSFB are obvious when
7、 the heat transfer coefficient of particles is low,the material heat capacity is small and the inlet temperature is low.Key words:multi-jet spout-fluidized bed;computational fluid dynamics(CFD);particle mixing characteristics;flow characteristics 1 引 言 与常规喷动床相比,多喷嘴喷动流化床具有颗粒相与气相混合更均匀、流动更快、颗粒死区更小、壁面黏结
8、更少等优点1,具有良好的研发优势。计算流体力学-离散单元(CFD-DEM)耦合模型对连续相求解 N-S 方程,对颗粒相求解牛顿第二定律,在喷动床研究中广泛应用2。Yue 等3用 CFD-DEM 模型研究喷动床内颗粒密度对喷动偏转的影响,用喷动偏转角度量化喷动偏转行为,结果表明交替喷动偏转的流体力学特性与初始状态和颗粒密度关系不大。Liu 等4用 CFD-DEM 耦合方法研究了不同密度颗粒的喷动特性,得到了不同密度和流速下的流型图,结果表明,最小喷动气速、床层压降和稳态喷动气速范围随颗粒密度的增加而增加。目前,采用 CFD-DEM 方法对多喷嘴喷动流化床内颗粒流体力学和传热特性的 收稿日期:20
9、22-03-22;修订日期:2022-06-28。基金项目:国家自然科学基金(22178286)。作者简介:刘雪峰(1997-),女,山西吕梁人,西北大学硕士生。通信联系人:吴峰,E-mail: 引用本文:刘雪峰,白谨豪,吴峰,马晓迅,危仁波.多喷嘴喷动流化床内气固两相流动与传热过程 CFD-DEM 模拟 J.高校化学工程学报,2023,37(3):383-391.Citation:LIU Xuefeng,BAI Jinhao,WU Feng,MA Xiaoxun,WEI Renbo.CFD-DEM simulation of gas-solid two-phase flow and heat
10、 transfer in multi-jet spout-fluidized bed J.Journal of Chemical Engineering of Chinese Universities,2023,37(3):383-391.384 高 校 化 学 工 程 学 报 2023年6月 研究未见报道。本研究基于欧拉-拉格朗日框架下的 CFD-DEM 方法,采用 Fluent15.0 分析多喷嘴喷动流化床内颗粒的流动及混合特性,并引入 Li 和 Mason5传热模型研究不同操作参数对传热的影响,将模拟结果与常规喷动床比较,为多喷嘴喷动流化床的推广应用提供理论和应用参考。2 模型方程 2.
11、1 固相 CFD-DEM 模拟方法对固相采用离散单元法(DEM)求解。通过牛顿第二定律和角动量守恒对单个颗粒进行受力分析,计算公式如下6:ntdragpddiiimmt=+vgFFFF (1)d diiiIt=M (2)式中:im为颗粒 i 的质量,kg;iv为颗粒 i 的泊松比,t为时间,s;g为重力加速度,ms?2;iI为颗粒i 的旋转惯量,kgm?2;i为颗粒 i 的角速度,rads?1,Fn、Ft、Fdrag、Fp分别为法向力、切向力、阻力、流体压力梯度;Mi为接触面扭矩。nF是法向重叠量n的函数,定义为7*3/2*relnnnn45236ERS m=+Fnv (3)式中:n 为法向单
12、位矢量;relnv为相对速度的法向分量,ms?1。等效杨氏模量*E(Pa)、颗粒等效半径*R(m)、颗粒等效质量*m(kg)、阻尼系数(Pasm?1)、法向刚度nS(Nm?1)分别定义为()()22*111ijijvvEEE?=+(4)*111ijRRR=+(5)*111ijmmm=+(6)22lnlnee=+(7)*nn2SER=(8)式中:jv为颗粒j的泊松比;jm为颗粒j的质量,kg;e为颗粒恢复系数;iE、jE分别为颗粒i、j的杨氏模量,Pa;iR、jR分别为颗粒i、j的接触球体的半径。tF8、dragF8、接触面扭矩iM9和流体压力梯度力pF10的计算公式如下:*reltnttt58
13、26GRS m=?+Ftv (9)()3draggpg16id=?Fuv (10)rniiiR=?MF (11)pp,giVp=?F (12)式中:*G为颗粒等效剪切模量,Pa;t为切向重叠量,m;t为切向单位矢量,tS为切向刚度,Nm?1;g为气体体积分数,gu为气体速度,ms?1;pd为颗粒直径,m;r为滚动摩擦系数,m;p,iV为颗粒体积,m3;gp为气相压力梯度,Pam?1;reltv为相对速度的切向分量,ms?1。2.2 气相 气相被当作连续相,通过N-S方程获得每相的运动信息。质量和动量守恒方程如下11:()()ggggg0t+=u (13)第 37 卷第 3 期 刘雪峰等:多喷嘴
14、喷动流化床内气固两相流动与传热过程 CFD-DEM 模拟 385 ()()()ggggggggggggggdpt+=?+?uu ugF (14)()()()gg gggggggggtgg,1vkiic Tc TkkTQt =+=+u (15)式中:g为气相密度,kgm?3;gc为气体的比热容,Jkg?1K?1;gT为气体温度,K;gk为气相热导率,Wm?1K?1;tk为气相的湍流热导率,Wm?1K?1;湍流普朗特数tPr设置为0.85,g,iQ是计算网格中气体和颗粒i之间的热流,J;g为气相应力张量;g为气相雷诺应力张量;Fd为颗粒流体阻力,N。g 12、g 12、Fd 13分别为()()()
15、1gggggg23?=+?|uuuI (16)()()()1gggg2233tt?=+?|uuuI (17)()dgs=?Fuu (18)式中:g为气体黏度,Pas;I 为颗粒应力张量;t为湍流黏度,kgm?1s?1;us为固相速度,ms?1。同时,引入Gidaspow曳力模型进行计算,计算公式如下14-15:()()()ggggs2.7D0ggp2gggsgg2gpp13 0.84111501.75 0.8Cddd?|=?|+|uuuu (19)()0.687pppD0p2410.15 1 0000.44 1 000ReReReCRe+|=|(20)式中:为动量交换系数,kgm?1s?1;D
16、0C为单颗粒的曳力系数;pRe为颗粒的雷诺数。气相湍流用k?湍流模型描述,湍流动能k和耗散率可用如下方程表示16:()()ggggggggkggktkkkGt+=+?|u (21)()()2gggggggg1k2gtCGCtkk+=+?|u (22)式中:1C、2C、k、的值分别为1.44、1.92、1.0、1.3,并且,kG表示平均速度梯度引起的湍流动能,J。2.3 传热模型 气体的能量方程的求解形式为17()()()()effeffhjjjEEpkThSt+=?+|vJv (23)式中:eff kT、jjjhJ、()effv 分别表示由热传导、物质扩散和黏性耗散导致的能量转移,v表示速度,
17、ms?1;effk表示有效热导率,Wm?1K?1;jJ表示物质j的扩散通量,kgm?2s?1;jh表示物质j的焓,J;eff 表示有效应力张量,Pa;hS表示热源相。Chaudhuri18提出的传热模型对DEM模拟中颗粒间的传热表征是通用的。Wei等19通过该模型研究了不同形状颗粒的传热特性,结果表明模拟结果和实验数据具有一致性。颗粒之间的热通量由下式计算18:().,i ji jjiQhTT=?(24)式中:,i jh为颗粒间传热系数,Wm?2K?1;jT和iT分别表示颗粒j和颗粒i的温度,K。通过Li和Mason5提出的传热模型对对流换热系数pg,ih进行计算:p,gpg,p,iiiNu
18、khd=(25)386 高 校 化 学 工 程 学 报 2023年6月 表 1 模拟使用颗粒的详细属性和模拟参数Table 1 Detailed properties and particle parameters used in simulation Property Symbol ValueGas density g 1.225 kgm?3Gas viscosity g 1.789 410?5PasApparent gas velocity ug 1.5,2.5 ms?1Time step CFD tCED 110?5sDEM TDEM 110?6sParticle Number N 1
19、200,4 800Diameter dp 2.8210?3mDensity p 2 503 kgm?3Shear modulus Gp 5.8109PaPoisson ratio vp 0.25WallShear modulus Gw 71010 PaPoisson ratio vw 0.3 Recovery factor Particle-particle ep-p 0.75Particle-wall ep-w 0.9 Static friction coefficient Particle-particle s,p-p 0.5 Particle-wall s,p-w 0.24Rolling
20、 friction coefficient Particle-particle r,p-p 0.06Particle-wall r,p-w 0.05表2 初始和边界条件的设置 Table 2 Setting of initial and boundary conditions Calculation boundaries Parameter setting Gas inlet Velocity direction perpendicular to inlet boundary;turbulence intensity of 3%,viscosity rate of 0.001 9;Gas ou
21、tletPressure outlet boundary Wall Gas phase:no slip Solid phase:no slip Fluid The initial temperature is 298 K,the filling height of large particles is 0.107 m,the filling volume fraction is 0.55,the axial pressure of the particles is corrected图 1 多喷嘴喷动流化床的示意图(unit:mm)Fig.1 Schematic diagram of mult
22、i-jet spout-fluidized bed(unit:mm)15219 5060750 multi-jet 3(a)IMJSFB(b)structure of multi-jet(c)grid3.51/21/3gp,p,3.51/21/33.50.81/3p,gp,gp,p,3.51.8gp,20.6 20020.50.02 2001 50020.000 045 iiiiiiiRePrReNuRePrRePrReRe+=+|(26)g,ggpCPrk=(27)式中:p,iNu表示颗粒-流体相间传热的努塞尔数;p,id表示粒径,m;p,iRe表示颗粒-流体相间雷诺数;Pr表示普朗特数;,
23、gpC表示气体的比定压热容,Jkg?1K?1。单个颗粒的温度随时间的变化19与气体和颗粒之间的对流传热17公式表示为 ppheatddTm CQt=(28)(),gpg,p,gp,iiiiQhATT=?(29)式中:pm、pC、heatQ分别表示颗粒的质量、比热容以及对流和传导热通量的总和。p,iA表示颗粒-流体的接触面积,m2;p,iT表示颗粒i与流体相混合的温度,K。颗粒-流体-颗粒间的热传导的计算公式为18:()()outinpgp,gp,p,22222dRijjiRijijxQkTTxlrxrx=?+?(30)式中:ijl表示颗粒j与颗粒i之间的距离,m;inR和outR表示颗粒-流体
24、-颗粒热传导的积分边界,ir、jr分别表示颗粒i和颗粒j的半径,m;Tp,j表示颗粒j与流体相混合的温度,K;x为颗粒自身的半径。3 数值模拟及边界条件 本研究的多喷嘴喷动流化床为侧壁开孔,喷嘴数量为1对,喷嘴的直径为3 mm,如图1所示,颗粒属性以及模拟参数见表1。恢复系数及摩擦系数与文献20一致,通过FFT功率谱分析实验和模拟的压力波动情况,结果表明实验和模拟结果吻合,参数选取合理。在颗粒混合过程中,2个喷动床的总气体体积流量一致,进气速度为1.875 ms?1。为了更好地比较颗粒的流动特性,颗粒特性与杨春玲21的研究保持一致。初始与边界条件设置如表2所示。采用k?湍流模型分析湍流的影响,
25、对体积分数项采用一阶迎风离散方案,对与动量、湍流动能和湍流耗散率有关的方程采用二阶迎风离散方案。多喷嘴喷动流化床的侧开孔流量由主喷嘴喷流实现,侧开孔速度由喷动床的进口速度决定。第 37 卷第 3 期 刘雪峰等:多喷嘴喷动流化床内气固两相流动与传热过程 CFD-DEM 模拟 387 表 3 模拟工况汇总 Table 3 Summary of simulated working conditions Working conditionsType of spouted bedug/(ms?1)Number of particles Logarithmic opening Case A1 CSB 1.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 喷嘴 流化床 内气固 两相 传热 过程 CFD DEM 模拟 雪峰
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。