《武汉工程大学学报》 2013年12期
63-67
出版日期：2013-12-31
ISSN:1674-2869
CN:42-1779/TQ
空气雾化硫磺喷枪的数值模拟
0引言 当代硫酸工业主要以硫铁矿和硫磺为原料制取硫酸.硫磺制酸比硫铁矿制酸更具有优势，主要体现在：减少工艺，消除污染源、能耗消耗下降、生产场地缩小、为企业提供了发展空间，所以硫磺制酸工艺应用较为广泛［1］.在硫磺制酸工艺中，对原料的利用及节能降耗有着重大影响的一个关键设备就是硫磺喷枪，对液硫的雾化有着直接影响.硫酸工业大多数采用的是直流单相喷枪雾化，其雾化效果不是很理想，对液硫的利用率不高.喷枪是一种特殊的喷嘴，相关学者对喷嘴的雾化机理探索多以实验以及经验理论为主，从直观的角度获取喷枪雾化特性参数，具有较大的局限性［2］.随着计算机技术和计算流体力学技术的发展，可以利用Fluent软件对喷嘴外部环境进行数值模拟，得出更宽阔更准确雾化特性参数. 喷枪按照硫磺的雾化方式，通常分为机械雾化喷枪和空气雾化喷枪.机械雾化是利用高液体压力经小孔高速射出而雾化，空气雾化是利用空气产生的直流、旋流或者交叉射流冲击液流束而雾化.本文以机械喷枪外加辅助空气为研究对象，以工艺参数为参考对象，利用计算流体力学软件Fluent对焚硫炉内喷枪的外部流场及雾化过程进行数值模拟分析，为生产实践提供参考.1模拟对象的几何结构焚硫炉及喷枪结构如图1所示（各部分结构未按比例绘制）.以某焚硫炉及配套机械喷枪（SG\|Ｉ）为例，对喷枪结构进行局部改进，在喷嘴四周通入压缩空气，结构参数如表1所示.图1焚硫炉Fig.1Burning sulfur furnace表1模拟对象基本尺寸表Table 1The mock object basic sizemm硫磺喷枪喷嘴直径d1辅助空气进口直径d2排气口直径d3焚硫炉直径d4焚硫炉有效长度L102003002 0006 000给定的工艺参数：该配套装置硫酸产量为800 t/d，计算出喷枪液硫流速为5 m/s.其物料的物性参数如表2所示.表2物性参数表Table 2Physical property parameter温度t/℃密度ρ/(kg·m-3)粘度μ/(Pa·s)喷枪出口处液态硫磺1301 797 8.66×10-3空气进口处1300.873 23.15×10-62仿真过程 （1）模型建立.利用Gambit创建模拟对象的几何模型，模型的具体尺寸以实际尺寸为准（见图1）.该焚硫炉有两个进口，一个出口.工作过程如下：液态硫磺经过喷枪高速射入焚硫炉内，形成机械雾化，喷嘴四周设有空气进口，辅助液硫雾化.该模型为对称结构，故采用2d模型求解即可，模型比较单一采用结构性网格（即四边形网格），由于进出口附近流体的流动非常复杂，流体流动的梯度较大，而且是数值模拟的最重要区域［3］，所以在这些区域采用了加密网格的方法，从而细化网格提高计算精度，模型最后的网格总数目为133 000个. （2）选择计算模型.该模拟对象只涉及到液体和气体的两相流动，故选用VOF模型，由喷嘴处液硫速度及空气流速计算出雷诺数，以及相关学者得出的这类射流流动的经验，选择k\|epsilon湍流模型［4］. （3）边界条件设置.焚硫炉内气液边界条件较为简单，进口处以液硫、空气的实际速度为基准，出口处采用压力出口，且空气、液流的进口温度与出口温度相同，焚硫炉的壁面选用标准无滑移壁面并假设其绝热. （4）求解器控制参数的设置.由相关定义计算后可确定该模拟对象为非定常流动，宜采用压力\|速度耦合方程求解，并选用PISO算法，其计算收敛性较好，结果较为准确［5］.第12期杨红，等：空气雾化硫磺喷枪的数值模拟武汉工程大学学报第35卷3计算结果及对比分析 以喷嘴处液硫作为离散相取速度5 m/s为定值的参考对象，液硫喷入焚硫炉内，其目的是液硫能充分燃烧，故需要一定量的空气.由相关化学反应方程式可计算出当空气的速度大于0.07 m/s时，即可满足燃烧要求，其中过量空气可循环利用.本文只进行喷枪的雾化数值模拟，对空气流速是一个试探性的取值，通过分析对比雾化的速度矢量图，判断空气流速取值的范围.所以空气进口的速度分别取0 m/s、2 m/s、4 m/s、6 m/s、8 m/s、10 m/s进行数值模拟求解.图2～7所示为采用k\|epsilon模型，在液硫速度为5 m/s的条件下，不同空气速度下计算出的喷枪外部流场的速度矢量图.图2空气流速v=0 m/s时的速度矢量图Fig.2Air velocity v=0 m/s图3空气流速v=2 m/s时的速度矢量图Fig.3Air velocity v=2 m/s图4空气流速v=4 m/s时的速度矢量图Fig.4Air velocity v=4 m/s图5空气流速v=6 m/s时的速度矢量图Fig.5Air velocity v=6 m/s图6空气流速v=8 m/s时的速度矢量图Fig.6Air velocity v=8 m/s图7空气流速v=10 m/s时的速度矢量图Fig.7Air velocity v=10 m/s（1）图2为没有改进前（按输入空气的速度为0）进行的数值模拟，由速度矢量图可以看出：该喷枪的雾化角小，雾化行程较长，而且速度分布不均匀，液滴主要分布在喷嘴附近，此状况表明该喷枪的雾化效果较差. （2）图3～7为改进后（该结构能从喷嘴附近输入一定速度的空气）进行的数值模拟，由速度矢量图可以看出：当输入一定速度的空气之后，喷枪的雾化角、雾化行程变大，液滴较为均匀的分布在整个焚硫炉内.这是因为当有一定速度的空气从喷枪周围输入时，在焚硫炉的喷枪附近形成一个涡流负压区，加速液硫向四周扩算，增大雾化角及雾化行程，这一现象与伴随射流理论一致［6］，从而使液硫能够快速雾化. （3）由图6～7的速度矢量图可以看出：当输入空气的速度过大时，导致喷枪附近的涡流强度较大，使液硫雾滴回流，雾滴甩向焚硫炉的壁面，造成硫磺挂壁现象.这种情况不利于雾化，在焚硫炉内要避免发生. （4）综合图2～7可以得出：当喷枪周围输入一定速度的空气后，有利于雾化.在喷枪出口处液硫速度为5 m/s，空气进口速度为4~6 m/s时，焚硫炉内喷枪的雾化效果较为理想.4结语 利用fluent软件对焚硫炉内喷枪的雾化情况进行了数值模拟，通过改进局部结构，在特定的工艺参数下，进行了多组数值模拟，得到多组速度场的分布图，分析雾化角及雾化行程的变化关系，得到了该工艺条件下液硫与空气速度相关的耦合参数.即在喷枪出口处液硫速度为5 m/s，空气进口速度为4~6 m/s时，该工艺中喷枪的雾化效果较好，对生产实践有一定的指导意义，同时也为以后焚硫炉内的燃烧情况的数值模拟研究奠定了坚实的基础.致谢：感谢武汉工程大学研究生创新基金提供的资助，以及武汉都市环保工程技术股份有限公司王呈祥在数值模拟过程中给予的指导.