收稿日期:２０２０Ｇ１０Ｇ１０.作者简介:王坤,男,２０２１年将毕业于中国石油大学(华东)动力工程专业毕业.在读硕士,研究方向为动力工程,助教.E m a i l :１０８４１４９１３４＠q q ．c o m .基金项目:国家自然科学基金(５１８０６２４８);山东省自然科学基金(Z R ２０１８B B ０６８).湍动流化床颗粒流动特性的数值模拟王㊀坤,朱丽云,李安俊,王振波,王国涛(中国石油大学(华东)新能源学院,山东青岛２６６５８０)㊀㊀摘㊀要:采用欧拉Ｇ拉格朗日研究方法,通过数值模拟对湍动流化床流化特性开展研究,得出了床层内压力㊁固含率㊁颗粒速度及粒度分布规律.结果表明:床层中心处的压力沿轴向高度逐渐减小;固含率沿轴向高度逐渐减小,呈 ε型分布,分为底部密相区㊁中间过渡段及上部稀相区,在稀相区呈现出中间低㊁边壁高的 环核 结构,且该结构沿轴向逐渐减小,直至消失;颗粒轴向速度沿轴向起伏变大,中心区域的颗粒速度普遍大于边壁处的颗粒速度,且颗粒下落位置由边壁向中心存在移动;颗粒粒径分布沿轴向逐渐降低,稀相区分布相对比较均匀.关键词:湍动流化床㊀数值模拟㊀计算颗粒流体力学㊀流化特性d o i :１０．３９６９/j．i s s n ．１００６－８８０５．２０２１．０１．００２㊀㊀当前我国的重油催化裂化工艺(R F C C )经过不断改革创新,逐渐趋于科学化㊁高效化,已经成为石油加工的核心技术和重质油轻质化的主要手段.随着原油日益重质化㊁劣质化,对再生器等设备的考验越来越大.催化裂化反应使焦炭沉积在催化剂表面上,导致催化剂活性急剧下降ʌ１ɔ,而再生器则通过烧去结焦催化剂上的焦炭来恢复催化剂的活性.为了提高目标产品收率㊁降低焦炭产率,工艺上通过提高系统催化剂循环量㊁改变再生器内气固流化形式来提高催化剂再生效率,促进裂化反应,保证催化剂活性.系统催化剂循环量增加和再生器内流化形式的变化均需要更高的流化气速,随着再生器内操作气速的提高,床层内气体夹带能力大大提高,大量的催化剂被带入稀相空间,降低了再生效率;同时气固分离设备入口浓度增大,分离负荷增加,设备效率降低,也降低了催化剂的回收率.针对上述问题,在深入研究湍动流化床内颗粒流动特性的基础上,对流化床及旋风分离器的结进行优化㊁改善催化裂化条件是十分必要的.湍动流化床因其固含率高㊁气固接触良好㊁热质传递高效等独有的特点在目前F C C 再生器中的应用已十分广泛ʌ２Ｇ３ɔ,而由于湍动流化床内气固运动复杂㊁各种测量技术的局限及实验条件存在限制,目前针对湍动流化床中颗粒运动规律的研究仍较少.随着数值算法的发展以及各种商用软件的兴起,计算颗粒流体动力学(C P F D )ʌ４ɔ已经逐渐成为研究气固两相流的重要手段.计算颗粒流体动力学(C P F D )数值模拟方法采用欧拉Ｇ拉格朗日方法来求解气固的流动特性ʌ５ɔ,该方法基于M P ＧP I C (m u l t i ph a s e p a r t i c l e －i n Ｇc e l l )ʌ６ɔ,可以解决流化床等大型设备中颗粒与流体的耦合问题.欧拉Ｇ拉格朗日方法是将流体视作连续相,在欧拉坐标系下考察流体相的运动;颗粒视作离散相,在拉格朗日坐标系下通过对每个颗粒进行跟踪来研究颗粒的运动,即颗粒轨道模型,可对具有各种颗粒类型㊁尺寸㊁形状和速度的流动进行分析ʌ７ɔ.从现象上看,颗粒轨道模型是最适合研究颗粒㊁流动两相系统的宏观结构,它可以给定单个颗粒的物理特性,跟踪颗粒的蒸发㊁挥发和反应过程,并可以考虑气体和颗粒间的耦合作用,以及颗粒间的相互碰撞行为.同时模拟计算中的颗粒并非等同于实际意义的颗粒,而是将具有相同物理性质的颗粒打包成颗粒群(称之为 计算颗粒 ),将上亿的颗粒系统进行简化,从而大大提高计算效率,为准确高效模拟气固系统中颗粒流动特性提供了可静设备㊀㊀㊀㊀石油化工设备技术,２０２１,４２(１)５ P e t r o c h e m i c a l E q u i p m e n tT e c h n o l o g y行的条件.现代工业中的技术条件已变得越来越专业化㊁技术化,因此通过数值模拟了解流化床再生器内的颗粒浓度㊁速度及粒度分布等特性,为再生器及分离设备结构的设计与优化提供有意义的参考,对如今催化裂化工艺非常有帮助.１㊀数学模型以C P F D技术为核心ʌ８ɔ,采用欧拉Ｇ拉格朗日方法模拟颗粒Ｇ流体动力学,其中固体颗粒相采用离散拉格朗日数值颗粒模型,流体相采用欧拉计算网格模型,其技术方法路线如图１所示.图１㊀技术方法路线㊀㊀采用欧拉Ｇ拉格朗日两相模型,求解流体的基本守恒方程如表１所示.表１中的符号解释见文后.表１㊀数学控制方程２㊀几何模型及参数设定２．１㊀几何模型如图２所示,流化床结构主要是由气体分布器㊁床体㊁旋风分离器及料腿组成的.该装置模拟石油工业中催化裂化单元的再生系统装置,空气通过鼓风机压缩并进行缓冲,然后经过控制阀与转子流量计进入到气体分布器中,经气体分布器均匀分布并穿过整个流化床后,在流化床顶部进入二级串联旋风分离器完成初步气固分离;分离后的气相进入布袋除尘器和湿式除尘器除尘后排放,固相则沿旋风分离器料腿进入床层内进行循环利用.为了简化计算,本文只选择流化床主体作为研究对象,其详细结构尺寸如图３所示,其中z代表浓度测量开孔高度,在浓度测量孔上下距离０．２５m各开有１个压力测量孔,高度用h表示.流化床结构为下小上大变径式,其中小径段由２个筒节组成,大径段由８个筒节组成.流化床小径段底部设有气体分布器.分布器采用管式结构,上部开口数为４８个,孔径为０．００６m,如图２右下角所示.采用S o l i d w o r k s进行建模,通过笛卡尔网格划分功能对流化床进行网格正交划分,网格示意如图４所示.对流化床局部结构进行网格加密,该模型共设有７０２２４０个网格.２．２㊀模拟参数与边界条件实验中的颗粒相采用F C C催化剂,颗粒密度为２６００k g/m３,平均直径为６４．６２μm,粒度分布P S D如图５所示;气体介质为空气,密度为１．１７６ 石㊀油㊀化㊀工㊀设㊀备㊀技㊀术２０２１年㊀k g /m ３,气体粘度为１．８４ˑ１０－６k g /(m s ).定义入口为速度边界条件,出口为压力边界条件,数值模拟参数如表２所示,模拟初始条件云图见图６.图２㊀流化床结构图３㊀流化床几何模型３㊀结果分析与讨论３．１㊀模型验证为检验采用的数学模型及使用的模拟计算方法是否准确㊁合理,本文沿轴向取任意８组高度差,将对应高度差下数值计算获得的压降值与实验结果进行对比,各高度差代表的具体位置见表３.图７为压降的实验值与模拟值对比.如图７所示,任意高度差下压差的实验值与模拟值吻合较好,且整体变化趋势基本一致,由此可见该模拟方法能较为准确地反映出湍动流化床内气固两相的流动特性.图４㊀流化床网格示意图５㊀催化剂粒度分布表２㊀数值求解参数㊀㊀通过将采用不同曳力模型获得的密相床层压降值与实验值进行对比,获得与实验情况最接近的气固曳力模型.如图８所示,模拟情况下曳力模型H a i d e rL e v e n s p i e l 与R i c h a r d s o n 获得的密相床层压降最接近实验值,而R i c h a r d s o n 模型仅７ ㊀第４２卷第１期王㊀坤等．湍动流化床颗粒流动特性的数值模拟适用于单颗粒尺寸分布的模拟,故选择H a i d e rL e v e n s p i e l 曳力模型进行湍动流化床的数值模拟.图６㊀数值模拟初始条件云图图７㊀轴向压降实验值与模拟值对比表３㊀ΔH 代表的高度差图８㊀不同曳力模型密相区压降对比３．２㊀床层内颗粒流动特性分布图９为流化床中心区域压力沿轴向的时均分布.由图９可见,中心区域的压力随着轴向高度的增加而减小,在流化床的密相段颗粒浓度高,气泡运动剧烈,所以此处压力变化较大,而稀相区及上部自由空域颗粒浓度很小,气体介质流化比较平稳,压力变化较小.图９㊀时均轴向压力㊀㊀图１０表示湍动流化床沿轴向的时均固含率分布.由图１０可见,随着轴向高度的增加,颗粒浓度呈 ε形分布,总体上颗粒浓度随着轴向的增高而逐渐减小.根据轴向颗粒浓度,可以将流化床分为底部密相区㊁中间过渡段以及上部稀相区,其中稀相区中固含率沿轴向高度不变的区域为自由空域.颗粒浓度在密相区较高,稀相区较低,而中间的波动区域则是由过渡段内存在不规则的气泡运动所造成的.图１０㊀轴向固含率分布㊀㊀图１１表示湍动流化床内沿径向的时均固含率分布.由图１１可见:在轴向高度z １＝３．３２m 处,颗粒浓度在边壁处较高,在中心处较低.这主要是因为湍动流化床内的大部分气体都是从床层中心通过的,气体对颗粒的脉动比较剧烈,而在壁８ 石㊀油㊀化㊀工㊀设㊀备㊀技㊀术２０２１年㊀面处颗粒浓度较大便会发生团聚,所以在此处沿径向形成了 环核 结构;在z ２＝４．３２m 处边壁处的颗粒浓度与中心处相差不大,说明此处 环核结构已经消失,即沿轴向 环核 结构呈现递减的趋势,直至消失;随着轴向高度的增加,颗粒浓度沿径向的分布变得更加均匀,z ３＝５．３２m 处整个径向截面处的浓度已经十分接近.图１２为轴径向固含率云图分布,图中z ＝２．３２m 处也可明显看出边壁处的颗粒浓度较高,且总体分布规律与图１１的描述一致.图１１㊀径向固含率的分布图１２㊀轴向固含率的云图分布㊀㊀图１３表示湍动流化床内不同床层高度下颗粒的时均轴向速度分布.由图１３可见:流化床中心区域的颗粒速度普遍要大于边壁处的颗粒速度,且随着轴向高度的增加,颗粒轴向速度在径向分布的起伏变大,但总体分布趋势相同;在轴向高度z ＝５．３２m 处,此时颗粒处于稀相区,颗粒轴向速度从中心到边壁呈现出先减小后增大的趋势,且在中心区域附近也存在颗粒下落的情况,壁面处的颗粒不再下落,说明颗粒下落位置由壁面向中心存在移动.图１４表示湍动流化床内不同轴向高度颗粒的瞬时轴向速度分布云图.由图１４可以看出:由于气泡的剧烈运动,使得分布器区域的颗粒速度产生不规则的变化;在远离分布器区域的密相区,颗粒速度明显呈现出中间高㊁沿壁面速度逐渐降低的趋势,且在壁面处颗粒为向下运动,所以颗粒在密相区就会形成颗粒的循环流动;而稀相区的颗粒速度径向分布相对比较均匀,在边壁处也很少出现颗粒下落的情况,但颗粒速度总体仍呈现出中心高㊁边壁低的趋势.图１３㊀径向颗粒速度分布图１４㊀径向颗粒速度的云图分布㊀㊀图１５表示流化床内颗粒粒径沿轴向的分布.由图１５可见,颗粒粒径沿轴向高度的增加逐渐降低,在稀相区粒径分布相对比较均匀.在密相床层表面,由于气泡破裂时的夹带,使大量颗粒被抛向自由空域,而颗粒的扬析现象使颗粒在流化气体的作用下因具有不同的终端沉降速度而分级,在不同的轴向高度出现颗粒回落,较大的颗粒及颗粒团便会集中在密相区及过渡区,此处由于气泡的作用存在颗粒循环,因此颗粒在密相区及过渡区的分布相对不均匀,而较细的颗粒便存在于稀相区顶部.４㊀结语基于欧拉Ｇ拉格朗日双流体模型,采用颗粒流体动力学方法,通过数值模拟研究了湍动流化床压力㊁固含率㊁颗粒速度㊁粒度分布等流动特性,模９ ㊀第４２卷第１期王㊀坤等．湍动流化床颗粒流动特性的数值模拟拟结果表明:图１５㊀颗粒粒度的云图分布㊀㊀１)通过对采用不同曳力模型数值模拟获得的密相床层压降值与实验值进行对比,最终确定曳力模型H a i d e rL e v e n s p i e l 最接近实际情况.２)湍动流化床中心区域的压力随着轴向高度的增加而减小,在密相段压力变化较大,而稀相区压力变化较小.３)根据流化床内轴向颗粒浓度的分布,可以将床层分为底部密相区㊁中间过渡段以及上部稀相区,中间的浓度波动区域是过渡段内存在不规则的气泡运动所造成的;在密相区上部,颗粒浓度沿径向呈现中间高㊁边壁低的 环核 结构,且该结构沿轴向高度逐渐减小,直至消失.４)流化床中心区域的颗粒速度普遍要大于边壁处的颗粒速度,且随着轴向高度的增加,颗粒速度在径向分布的起伏变大,但总体分布趋势相同;在稀相区上部,颗粒下落位置由壁面向中心存在移动.５)颗粒粒径沿轴向高度的增加逐渐降低,在稀相区粒径分布相对比较均匀.符号说明㊀㊀t时间,;m 颗粒质量,k g;εg气相体积分数;p气相压力,P a ;ρg 气体密度,k g /m ３;F单位体积内气固两相动量交换律;u g 表观气速,m /s;g重力加速度,m /s ２;τg气相应力张量,P a ;D相间交换系数;D p曳力系数;f概率分布函数;u p 颗粒速度,m /s;εs单位网格内颗粒体积分数;ρp 颗粒密度,k g /m ３;X p颗粒位置向量;εp固相体积分数;τp颗粒法相应力,N ;d u pd t颗粒加速度,m /s ２;m p 颗粒质量,k g ;μf 流体粘度,kg /m s ;r p颗粒半径,u m ;R e雷诺数;V p 颗粒体积,m３;F D 气固相间曳力,N ;f h与雷诺数相关的参数.参考文献:[１]㊀L IAJ ,WA N GGT ,WA N GK ,e t a l ．M a c r o s c o pi c a n d M i c r o s c o p i c C h a r a c t e r i s t i c s o f P a r t i c l e si n a N o v e lG a s Ｇs o l i dC yc l o n eR e a c t o r [J ]．C h e m i c a lE n Ｇg i n e e r i n g R e s e a r c h a n dD e s i gn ,２０１９,１５２:２８８Ｇ２９９．[２]㊀B IH T ,E L L I S N ,A B B AIA ,e t a l ．A S t a t e Ｇo f Ｇt h e Ｇa r tR e v i e w o fG a s Ｇs o l i d T u r b u l e n tF l u i d i z a t i o n [J ]．C h e m i c a lE n g i n e e r i n g Sc i e n c e ,２０００,５５(２１):４７８９Ｇ４８２５．[３]㊀D U ,B ,WA R S I T O W ,F A N LS ．B e dN o n h o m o Ｇg e n e i t y in T u r b u l e n t G a s Ｇs o l i d F l u i d i z a t i o n [J ]．A I C h EJ o u r n a l ．,２００３,４９(５):１１０９Ｇ１１２６．[４]㊀S N I D E RD M．A n I n c o m pr e s s i b l eT h r e e Ｇd i m e n s i o n Ｇa lM u l t i ph a s eP a r t i c l e Ｇi n Ｇc e l lM o d e l f o rD e n s eP a r t i Ｇc l eF l o w s [J ]．J o u r n a lo fC o m p u t a t i o n a lP h y s i c s ,２００１,１７０(２):５２３Ｇ５４９．[５]㊀C H E N C ,W E R T H E R J ,H E I N R I C H S ．C P F DS i m u l a t i o no fC i r c u l a t i n g Fl u i d i z e db e d R i s e r s [J ]．P o w d e rT e c h n o l o g y ,２０１３,２３５:２３８Ｇ２４７．[６]㊀M．J ．A n d r e w s ,P ．J ．O ＇R o u r k e ．T h e m u l t i ph a s e pa r t i c l e －i n Ｇc e l l (M P ＧP I C )m e t h o d f o r d e n s e p a r t i c Ｇu l a t e f l o w s [J ]．I n t e r n a t i o n a l J o u r n a l o f M u l t i p h a s e F l o w ,１９９６,２２(２):３７９Ｇ４０２．[７]㊀P ．A．C u n d a l l ,O．D．L ．S t r a c k ．T h ed i s c r e t en u Ｇm e r i c a lm o d e l f o r g r a n u l a r a s s e m b l i e s [J ]．G e o t e c h Ｇn i qu e ,１９７９,(２９):４７Ｇ５．[８]㊀张向阳．工业级流体颗粒系统分析软件B a r r a c u d a [J ]．C A D /C A M 与制造业信息化,２０１１,(１１):２９Ｇ３２．[９]㊀HA I D E R A ,L E V E N S P I E L O．D r a g Co e f f i c i e n t a n dT e r m i n a lV e l o c i t y o fS p h e r i c a l a n d N o n s p h e r i Ｇc a l P a r t i c l e s [J ]．P o w d e rT e c h n o l o g y ,１９８９,５８(１):６３Ｇ７０．０１ 石㊀油㊀化㊀工㊀设㊀备㊀技㊀术２０２１年㊀K E Y P O I N T S O F D E S I G N A N D M A N U F A CＧT U R E O F S E P A R A T E D O U B L E T U B E S H E E T S E X C H A N G E R[１]X uW e i(S I N O P E CE n g i n e e r i n g I n c o r p o r a t i o n, B e i j i n g,１００１０１)A b s t r a c t:S h e l la n d t u b e h e a te x c h a n g e r s a r e m o s t l y o fs i n g l et u b e s h e e ts t r u c t u r e,a n db o t h s i d e s o f t u b e s h e e t c o n t a c t w i t h t h e m e d i u m r e s p e c t i v e l y．O n c et h ec o n n e c t i o n b e t w e e nt h e t u b e s h e e ta n dt h et u b e s f a i l s,t h e m e d i u m w i l l l e a k．A n i s o l a t i o n c a v i t y i s s e t b e t w e e n t h e i n n e r a n do u t e rt u b e s h e e t si nt h ed o u b l et u b e s h e e t s h e a t e x c h a n g e r,w h i c h c a n e f f e c t i v e l y s o l v e s u c h p r o b l e m s．T h i s p a p e rd e a l sw i t ht h ea n a l y s i so f m a n y k e y p o i n t s w i t h e n g i n e e r i n g d e s i g n e x a m p l e s,s u c h a s o c c a s i o n s o f d o u b l e t u b e s h e e t sa p p l i c a t i o n,c a l c u l a t i o no f i n n e ra n d o u t e r t u b e s h e e t s t h i c k n e s s,f a c t o r s f o r d e t e r m i n i n g t h el e n g t ho ft h ei s o l a t i o nc a v i t y,a n d t ub e s h e e t s m a n u f ac t u r i n g,c o n n e c t i o nb e t w e e nt u b e s h e e t sa n d h e a te xc h a n g et u b e s, p r e s s u r e t e s t,e t c．T h e p a p e r a l s o p r o v ide s o p t i m i z a t i o n a n d i m p r o v e m e n t m e a s u r e s t o b e t t e r t h e d e s i g na n dm a n uf a c t u r i n gq u a l i t y．K e y w o r d s:d o u b l et u b e s h e e t;h e a te x c h a ng e r; s e p a r a t i o nw a l l;w e l d i n g;s t r e n g t he x p a n s i o n N U M E R I C A L S I M U L A T I O N O F P A R T I C L E F L O W C H A R A C T E R I S T I C S I N T U R B U L E N T F L U I D I Z E DB E D[５]W a n g K u n,Z h uL i y u n,L iA n j u n,W a n g Z h e n b o, W a n g G u o t a o(C o l l e g eo f N e w E n e r g y,C h i n a U n i v e r s i t y o f P e t r o l e u m(E a s tC h i n a),Q i n g dＧa o,S h a n d o n g,２６６５８０)A b s t r a c t:T h e E u l e r i a nＧL a g r a n g i a n t w oＧf l u i d m o d e l w a s u s e d t o s t u d y t h e f l u i d i z a t i o n c h a r a c t e r i s t i c s o f t h e t u r b u l e n t b e db y n u m e r i c a l s i m u l a t i o n．T h e p r e s s u r e,s o l i dh o l d u p,p a r t i c l e v e l o c i t y a n d p a r t i c l es i z ed i s t r i b u t i o no f t h eb e d w e r e s t u d i e d．T h er e s u l t ss h o wt h a t t h ec e n t r a l p r e s s u r e g r a d u a l l y d e c r e a s e s a l o n g t h e a x i a l h e i g h t,w h i l e t h e s o l i d h o l d u p g r a d u a l l y d e c r e a s e sa l o n g t h ea x i a le l e v a t i o n,s h o w i n g a ε Ｇt y p ed i s t r i b u t i o n．A c c o r d i n g l y,t h e b e di s d i v i d e d i n t o a b o t t o md e n s e p h a s e z o n e,am i d d l e t r a n s i t i o n z o n e,a n dat o p d i l u t e p h a s ez o n e．I n t h e d i l u t e p h a s e r e g i o n,t h e c o r eＧa n n u l a rＧs t r u c t u r e g r a d u a l l y d e c r e a s e s a l o n g a x i a l d i r e c t i o nu n t i l i td i s a p p e a r s．B e s i d e s,t h ea x i a l v e l o c i t y o f t h e p a r t i c l e sf l u c t u a t e s i nt h er a d i a l d i r e c t i o n,w h i c h i nt h ec e n t r a l z o n e i s g e n e r a l l y g r e a t e r t h a nt h a ta tt h e w a l l,a n dt h e p o s i t i o n w h e r e t h e p a r t i c l e s f a l l i n g m o v e s f r o ms i d ew a l l t o t h e c e n t e r．T h e a v e r a g e p a r t i c l e s i z e g r a d u a l l y d e c r e a s e sa l o n g t h ea x i a la x i s,a n dt h e p a r t i c l e d i a m e t e r d i s t r i b u t i o no f t h ed i l u t e p h a s ez o n e i s r e l a t i v e l y u n i f o r m．K e y w o r d s:t u r b u l e n tf l u i d i z e d b e d;n u m e r i c a l s i m u l a t i o n;c o m p u t a t i o n a l p a r t i c l eh y d r o d y n a mＧi c s/C P F D;f l u i d i z a t i o n c h a r a c t e r i s t i c sR E S E A R C H A N D A N A L Y S I S O FS W X S T A T I C M I X E RB A S E DO NF L U E N T[１１]H u J i q u a n,L iZ h i f u,X i e D a n l i n(S I N O P E C J i n l i n g C o m p a n y,N a n j i n g,J i a n g s u,２１００３３) A b s t r a c t:T h i s p a p e r u s e d P r o/E s o f t w a r et o p e r f o r m３D m o d e l i n g o f t h es t r u c t u r ea n df l u i d m o d e lo ft h eS WX s t a t i c m i x e r．B a s e do nt h e F L U E N T m o d u l eo f A N S Y S W o r k b e n c h,t h e p a p e r m a d e r e s e a r c h a n d a n a l y s i s．N u m e r i c a l s i m u l a t i o n s o f f l o wf i e l d,v e l o c i t y f i e l d,p r e s s u r e f i e l d a n d h e a t t r a n s f e r c h a r a c t e r i s t i c s w e r e c a r r i e do u t t o s t u d y t h e f l u i d f l o wc h a r a c t e r i s t i c sA B S T R A C T S㊀P E T R O C H E M I C A LE Q U I P M E N T T E C H N O L O G YS t a r t e dP u b l i c a t i o n i n１９８０．B i m o n t h l y．J a n．２０２１V o l．４２N o．１ Ⅲ。
