乙烯裂解炉辐射段燃烧流场三维数值模拟.docVIP

乙烯裂解炉辐射段燃烧流场的三维数值模拟 （北京化工大学化学工程学院，北京 100029） 摘要：本文用FLUENT6.2计算流体力学软件对燕山石化公司E-BA-107乙烯裂解炉辐射段内的流动和燃烧做了三维数值模拟研究。为保证计算的准确性，建立了与裂解炉实际尺寸完全一致的几何模型，计算中采用了标准湍流模型、DO辐射模型以及inite-rate/Eddy-dissipation燃烧反应模型等计算模型。通过计算,获得了炉膛内的温度、速度分布及烟气组分, 其中烟气的出口温度和组分与设计手册一致，表明本文的数值模拟有较高的精度，具有重要的工程实际意义。 关键词：乙烯裂解炉 燃烧 计算流体力学三维数值模拟言 裂解炉是乙烯生产装置中的一个关键设备，其性能对乙烯收率有重大影响。在裂解炉辐射段，既存在侧壁燃烧器和底部燃烧器的剧烈的燃烧反应，又存在裂解炉管中原料发生裂解反应造成的大量吸热，这就使其中的流动和燃烧现象极为复杂，用实验和现场测量的方法很难捕捉具体的燃烧细节，随着计算机硬件和计算流体力学( CFD )技术的发展, 近年来国内外已开始用数值模拟的方法研究乙烯裂解炉辐射段的详细流体流动、传热和燃烧反应过程。Plehiers等[1]用蒙特卡洛辐射换热模型模拟了乙烷裂解制乙烯的过程，并将炉膛内的辐射换热和炉管内的裂解反应进行了耦合计算，但在进行炉膛辐射计算时，把各个燃烧器简化为辐射点源。Stefanidis等模拟了裂解炉辐射段的燃烧流场，得出涡流耗散模型（EDC-DRK）比涡流破碎模型（EBU-SRK）模拟结果更为合理，但建立的模型只有底部燃烧器，没有侧壁燃烧器蓝兴英等和吴德飞等对乙烯管式裂解炉进行了数值模拟，得到了裂解炉内的速度、温度及烟气组份的分布情况，但，将反应管简化为直管没有考虑横跨段，烟气直接从裂解炉上部排出且对侧壁燃烧器也作了简化这些均与裂解炉的实际差别较大，会影响到速度、压力分布及传热计算精度，使计算结果与实际工况差别较大，也就失去数值模拟的工程实际意义。 基于上述情况，本文以燕山石化公司E-BA-107乙烯裂解炉辐射段为研究对象，完全按照裂解炉的工程图纸建立了用于数值模拟的几何模型，模型的建立 .1 几何模型 考虑到炉体结构的对称性，取其一半用GAMBIT2.2建立如图1所示的几何模型，模型尺寸为mm×5500 mm×12000 mm，炉体前后两侧墙上分别安置了12个侧壁燃烧器，燃烧器以一定距离排开，上下共3排，在炉体底部安置有两排共12个底部燃烧器，裂解炉管为GK5分枝型炉管，共有两组，安置在炉膛中央部位，图2炉膛底部炉管连接及底部燃烧器的布置情况。裂解原料气入口在炉体顶部，出口位于一侧墙上，炉体另外两个面一个是对称面，另一个是炉体壁面，炉体底面中心处（0，0，2.8 m），出口方向为正Z方向，各侧壁燃烧器的轴线方向为X方向。图3为侧壁燃烧器外形图，内设导流筒，燃烧器开孔有两个部位，开孔位置1有φ5 mm圆孔共108个，分3排，每排36个等距离布置，其中，圆孔与燃烧器轴心线有80°的倾角，开孔位置2方孔2.5 mm×20 mm共90个，与轴心线垂直。图4为底部燃烧器外形图，有5个喷头，中间部位的喷头有5个喷孔，Z向位置较低，位置较高的两个喷头各有3个喷孔，余下的两个喷头各有两个喷孔。为保证计算的准确性，在网格时对炉管和燃烧器部位分别进行了加密。 图1 裂解炉几何模型示意图 图2 炉管连接和底部燃烧器布置Fig.1 Geometry model of the cracking furnace Fig.2 Distribution of reactor tubes and bottom burners 图3 侧壁燃烧器外形 图4 底部燃烧器外形Fig.3 Side wall burner diagram Fig.4 Bottom burner diagram 2.2 计算模型 .2.1 基本方程 燃烧是包含激烈化学反应的流动过程，它们都遵守物理和化学基本定律，即质量、组分、动量和能量的守恒定律。Fluent6.2在计算过程中，需要解连续性方程、动量方程、能量方程以及组分方程等，压力和速度的耦合采用SIMPLE算法，总的方程表达式如下所示： （1） 其中，、、的含义详见文献［］ .2.2 湍流模型 在计算中,湍流模型方程采用Launder和Spalding提出的标准双方程模型。 湍动能方程： （2） 湍动能耗散率方程：