脉冲袋式除尘器喷吹气流的数值模拟袋式除尘器脉冲喷吹的清灰效果是影响设备运行阻力和滤袋使用寿命的主要因素之一. 如图1 所示,含尘气体进入中箱体,经滤袋过滤后洁净气体由上箱体出口排出. 喷吹时,压缩气体由喷嘴高速喷出,诱导喷嘴周边的数倍气体在短时间内进入滤袋,滤袋急剧膨胀、振动,从而使附着在滤料表面的粉尘层剥落[122 ] . 目前有关脉冲喷吹气流及其性能参数优化研究的报道文献不多,设计中对喷吹装置喷吹效果的判断,大多是根据经验,或者通过实物试验进行的[324 ] . 为了更好的了解喷吹气流与滤袋间的作用状况,本文采用计算流体动力学(Comp utational Fluid Dynamics ,简称CFD) 方法,使用Fluent 软件对脉冲喷吹气流的压力分布进行了数值模拟,对于滤袋清灰装置的优化设计具有重要的意义.1 计算模型喷吹气流在滤袋内的流场非常复杂,为了便于建立数值计算模型,对滤袋内外流场进行如下简化:(1) 以单条圆筒滤袋作为研究对象,并将其内外气流假定为不可压缩流体;(2) 在喷嘴出口处,由于喷嘴直径远小于滤袋直径,可以假定喷吹气流速度沿喷嘴径向均匀分布;(3) 常温计算条件,不涉及温度对流场的影响;(4) 不考虑滤袋壁面的纵向位移;(5) 喷吹气流在滤袋内的流场是三维流动问题,喷吹气流沿滤袋长度方向的轴线是对称的,可简化为二维问题.图2 (a) 为单条滤袋喷吹清灰的几何模型,X 为喷吹距离(喷嘴距袋口距离) ,L 为滤袋长度,喷嘴直径d ,D 为计算区域宽度.将引射空间、滤袋内外空间作为计算区域.网格划分时,沿袋口方向网格间距为D/ 20 ,沿滤袋长度方向网格间距为L/ 10.而对于喷嘴出口部分,因其对网格比较敏感,在网格划分时进行了适当的加密,以减小计算误差,提高结算精度. 图2 (b) 为计算区域的网格图,以及喷嘴部分的网格细化图.图2 单条滤袋计算区域Fig. 2 Calculated section of single filter bag2 数值计算方法2. 1 控制微分方程脉冲喷吹气流的流动过程采用的控制方程如下:2. 2 边界条件与初始条件固体壁面包括上箱体壁面、喷管外壁、花板、袋底,其边界条件采用壁面函数法[5 ] . 中箱体壁面采用恒压边界条件,依据压力的大小确定出滤袋的外表面过滤气速. 净气出口采用压力出口边界条件,滤袋采用多孔跳跃边界条件.滤袋介质作为渗流壁,其内部沿半径方向的流动方程由非稳态的Darcy 公式确定[6 ] .式中: v 为气体通过袋壁的径向速度,m/ s ; K 为滤袋壁渗透系数,取决于滤袋和粉尘层的几何结构和化学性质,m2 ;μ为黏性系数,Pa ·s ; p ( t) 为不同时刻的滤袋壁面压力值,Pa ; r 为径向距离,m.该流场的初始化是从恒压面开始初始化,恒压面的压力采用中箱体的压力.脉冲喷吹是一个非稳态的过程,非稳态湍流流场的计算采用SIMPL E 算法[7 ] . 工程上压缩空气的喷吹时间一般设定在80～150 ms 之间,在模拟计算时取100 ms ,求解步长取4/ 1 000 s.3 数值模拟结果3. 1 Φ160 ×6 000 滤袋的模拟滤袋的几何尺寸为直径Φ= 160 mm ,长度L = 6 000 mm. 喷吹时间T = 100 ms ,喷吹距离X = 200mm ,喷嘴d = 26 mm. 实验气包压力为0. 4 MPa ,实验采用的Goyen 淹没式脉冲阀的阻力损失在0. 06～0. 07 MPa之间[829 ] ,模拟时采用的压力为0. 33 MPa (数值模拟是从喷嘴处开始计算) . 由于该实验袋口部分采用了文丘里管,故在模拟时也增加了文丘里管.数值模拟得出滤袋各点的峰值压力曲线,如图3所示. 图中同时给出了文献[10 ] 在相同条件下的实验检测数据.由图3 可以看出,两条曲线基本吻合.在袋口0～0. 3 m 处,由于增加了文丘里管,峰值压力较小,且第二个测试点(0. 3 m处)的压力值有所减少; 而后压力开始上升,在滤袋1m 处达到最大值;随着滤袋长度的增加,压力值不断衰减; 到滤袋中下部,压力有所回升. 在距袋口0～3m 段,两条曲线完全重合,在滤袋中下游部分,模拟值高于实测值.这一偏差主要是因为实验测定与数值模拟的有效喷吹时间和滤袋的渗透系数取值不一致造成的,可以通过调整参数取值来修正.对比结果表明模拟值与实测所得峰值压力曲线基本吻合,滤袋内喷吹气流压力分布可以利用数值模拟的方法来分析.图3 滤袋侧壁峰值压力对比Fig. 3 Comparison of t he peak pressure between t he experimentalvalue and t he numerical calculations3. 2 Φ130 ×6 000 滤袋的模拟圆筒滤袋通常的尺寸为Φ130 、Φ150 、Φ160 mm 等几种,长度一般在5 000～8 000 mm 之间. 其中Φ130 ×6 000 mm的滤袋在工程应用中较多,对于长滤袋不宜加文丘里管,本文选择Φ130 ×6 000 mm 滤袋进行喷吹气流压力场的模拟.图4 脉冲喷吹不同时刻滤袋侧壁压力曲线Fig. 4 Pressure inside a filter bag at different instant s图4 是常温条件下Φ130 ×6 000 mm 的滤袋,在T = 100 ms , X = 200 mm , d = 24 mm , P = 0.2 MPa时脉冲喷吹不同时刻滤袋壁面的模拟结果.从图4 中可以看出,滤袋壁面不同部位达到最大压力值的时间是有差别的. 距袋口越远,达到最大压力值所需的时间越长. 40 ms 左右时,距滤袋口0～1m 处压力值达到最大;喷吹到60 ms 时,距袋口1～3 m 处达到压力最大值;70 ms 左右,袋口压力减小,距袋口3～6 m处滤袋壁面压力达到最大值;而在80 ms ,100ms时滤袋壁面的峰值压力值不断下降.在喷吹过程中,滤袋壁面受到的最大静压值即为峰值压力. 峰值压力是喷吹过程中一个最直观的表现形式,可以通过实验测量获得. 据图4 所示各时刻的压力分布,经计算后得图5所示的滤袋壁面峰值压力曲线. 由图5可见,峰值压力不断上升,距袋口约0. 5 m处出现最大值,随后由于喷吹气流泄漏和运动压力损失,峰值压力开始下降,喷吹气流冲击到封闭的滤袋底部,产生返流,导致滤袋底部的压力值回升.压力上升速率是滤袋壁面的峰值压力与压力从零值上升到峰值所经历的时间的比值.压力上升速率是喷吹效果的一个重要评价指标. 图6 是根据图5 所示的模拟结果得出的滤袋壁面压力上升速率曲线图. 由图6可见,压力上升速率从袋口开始上升,距袋口约0. 5 m 处出现最大值,随袋长的增加,压力上升速率逐渐减小.图5 滤袋侧壁峰值压力Fig. 5 Peak pressure during cleaning at variouslocations inside t he filter bag图6 滤袋最大压力上升速率Fig. 6 Most increscent velocity of pressure at variouslocations inside t he bag4 结论(1) 利用CFD 方法对Φ160 ×6 000 mm滤袋的脉冲喷吹清灰过程进行数值模拟,结果与文献实验测定结果的基本趋势一致,可以采用数值模拟的方法分析滤袋内喷吹气流的压力分布.(2) 由Φ130 ×6 000 mm 滤袋的数值模拟结果可以看出,在设定为100 ms的喷吹时间里,随着喷吹气流向下运动,滤袋壁面依次达到最大压力值. 70 ms以后,尽管还在继续喷吹,但壁面的峰值压力下降,表明延长压缩空气喷吹时间对峰值压力的提高不起作用.(3) 喷吹过程中,滤袋上部受到的峰值压力最大,底部其次,中部最小. 压力上升速率从袋口开始上升,距袋口约0. 5 m处出现最大值,随袋长的增加,压力上升速率逐渐减小.参考文献References[1 ] L ; FFL ER F ,SIEVERT J . Cleaning mechanisms in pulse2jet fabricfilters[J ] . Filt ration & Separation ,1987 ,24 (3/4) :1102113.[2 ] MORRIS W J . Cleaning mechanisms in pulse jet fabric filters[J ] .Filt ration & Separation ,1984 ,21 (1/ 2) :50254.[3 ] 杨复沫. 脉冲袋式除尘器清灰能力的评价方法和手段的试验研究[D] . 武汉:冶金部安全环保研究院,1992.Y AN G Fu2mo. Research of the assessment method on the cleaning ability ofpulse2jet fabric filters [D ] . Wuhan :Metallurgy Department Safety and Environmental Protection ResearchInstitute ,1992.[4 ] 陈隆枢. 脉冲喷吹装置的清灰特性及相关技术[ C]/ / 全国袋式过滤技术研讨会论文集. 2005 ,厦门. ,200322005(12) :2522259.CHEN Long2shu. The characteristic of cleaning and the interrelatedtechnology of pulse2jet equipment [C]/ / Discus2sion Forum Nemoir of Bag Filt ration of China (2005 ,Xiamen) ,2003 - 2005(13) :2522259.[5 ] 陶文铨. 数值传热学[M] . 第2 版. 西安:西安交通大学出版社,1999.TAO Wen2quan. Numerical heat t ransfer [M] . Second Edition. Xi′an : Xi′anJiaotong University Press ,1999.[ 6 ] KUNZ R F ,LA KSHMINARA Y ANA B. Explicit navier2stokes computationof cascade flows using the k2turbulencemodel [J ] . AIAA. 1992 ,30 (1) :13221.[7 ] 佟桂芳,徐德龙,张强,等. 一种新改进的SIMPL E 算法[J ] . 西安建筑科技大学学报:自然科学版,2001 ,33 (3) :2212224.TON G Gui2fang ,XU De2long ,ZHANG Qiang ,et al. A new modification tosimple algorithm[J ] . J . Xi′an Univ. ofArch. & Tech. (Natural Science Edition) ,2001 ,33 (3) : 2212224.[8 ] 李兴华. 脉冲阀性能检测装置的研究[D] . 武汉:武汉安全环保研究院,2002.L I Xing2hua. Study on detection equipment s of diaphragm valveperformance[D] . Wuhan : Wuhan Safety and Envi2ronmental Protection Research Institute ,2002.[9 ] 陈志炜. 大型低压脉冲阀的研究[D] . 武汉: 武汉科技大学,2003.CHEN Zhi2wei. Study on the large2scale and low2pressure diaphragm valve[D] . Wuhan : Wuhan University of Sci2ence and Technology ,2003.[10 ] 田玮. 脉冲喷吹袋式除尘器清灰的研究[D] . 西安: 西安建筑科技大学,2005.TIAN Wei. Study on the cleaning process of pulse2jet fabric filters[D].Xi′an : Xi′an Univ. of Arch. & Tech . ,2005.3Biography :DANG Xiao2qing ,Associate Professor ,Ph. D. ,Xi′an 710055 , P.R. China , Tel :02729 ,Email :dangxq #163.将#换成@com基金项目:国家863 计划项目资助课题(2005AA642010)作者简介:党小庆(19642) ,男,副教授,博士研究生,主要研究方向为大气污染控制技术和设备。
