锅炉过热器进口集箱三通区域静压分布的数值模拟及三通的结构改进.docVIP

锅炉过热器进口集箱三通区域静压分布的数值模拟及三通的结构改进 　　摘要： 为解决过热器超温爆管问题，根据电站锅炉过热器带三通集箱的结构特点，采用数值模拟方法对采用径向引入方式的过热器进口集箱三通区域的静压分布和流动规律进行了研究，并在此基础上对三通区域结构进行了改进.结果表明：三通区域集箱中存在涡流区，使该区域的静压变低，导致布置在该区域的支管蒸汽流量减少；适当增大三通区域集箱内径可提高集箱涡流区的静压，但当三通区域集箱内径增大到一定程度后，继续增大三通区域集箱内径，集箱涡流区的静压已不再有明显升高. 　　关键词： 　　带三通集箱； 数值模拟； 静压分布； 涡流区； 改进 　　中图分类号： TK 22文献标志码： A 　　目前大容量电站锅炉过热器大多采用带三通结构的径向引入、引出集箱系统，当蒸汽引入集箱后会在三通区域附近产生涡流区，使得布置在该区域的管屏中蒸汽流量减小.若该部分管屏也正好处于锅炉烟温较高的区域，就很容易发生超温爆管事故.由于三通区域集箱内流动状况复杂，静压分布规律与轴向引入、引出方式的集箱差别很大，运用于轴向引入、引出集箱系统压力分布规律的水动力计算方法[1]已不再适用，美国CE公司的计算标准中关于径向引入、引出集箱系统的计算方法也被证实是错误的[2].因此，研究并掌握径向引入方式过热器集箱三通区域的压力分布和流动规律，对过热器的设计和锅炉的安全运行具有重要的意义. 　　Kreid等[3]研究发现，与充分发展的直管中流速呈抛物线分布的情况相比，三通区域流体速度分布曲线发生了严重的扭曲.陆方等[2]和罗永浩等[4]对进口集箱三通区域压力分布规律进行了冷态模化试验研究，找到了一些影响三通区域静压分布的因素，并得到了用于计算三通区域静压分布的经验公式.匡江红等[5]、刘进等[6]、张润卿等[7]和周云龙等[8]采用数值模拟方法对进口集箱三通区域进行了研究，得到了集箱三通区域的静压分布规律、流动特性及相对应的支管流量分配情况.其中，周云龙等[8]还对三通涡流区的支管入口形状进行了改进，改进后支管流量明显增大. 　　从以上分析可以看出，国内外学者通过试验和数值模拟方法，对过热器进口集箱三通区域的静压分布和流动特性进行了研究，也取得了一些成果，但并没有根据这些研究成果对三通区域结构的改进提出具体方案.本文将通过数值模拟方法研究三通区域静压分布规律，在此基础上对三通区域结构提出改进方案，以期为工程实际提供参考. 　　1数值模拟 　　1.1研究对象 　　本文研究对象为有支管等径三通模型，等径三通内径D=200 mm，1～15号支管从集箱正母线和右下40°方向引出，支管直径为Φ 63.5×12.5 mm，支管间距112 mm.由于涡流区主要集中在距离三通径向引入管中心线±（2～3）D范围内[2，6-7]，所以本文选择距离三通径向引入管中心线-4D～4D范围内布置支管，模型如图1所示，集箱支管位置如图2所示. 　　1.2湍流模型 　　由于过热器进口集箱内的流动为湍流流动，所 　　以本文采用标准k-ε方程进行数值计算.其湍流 　　动能k方程为 　　x（ρuk）+y（ρvk）+z（ρwk）=xμσkkx+ 　　yμσkky+zμσkkz+Gk-ρε 　　（1） 　　湍流?幽芎纳⒙师欧匠涛? 　　x（ρuε） +y（ρvε） +z（ρwε）= 　　xμσε εx+yμσεεy+ 　　zμσεεz + εk（C1Gk-C2ρε） 　　（2） 　　式中：u、v、w分别为x、y、z方向的速度；ρ为密度；μ为动力黏度；Gk为由平均速度梯度引起的湍动能的产生项；σk为对于k的湍流普朗特数；σε为对于ε的湍流普朗特数；C1、C2为常数. 　　方程的离散选择二阶迎风差分格式；压力与速度的耦合选择SIMPLE算法. 　　1.3边界条件 　　壁面边界条件： 　　集箱壁面采用无速度滑移和无质量渗透边界条件，即假定相对于固体壁面的气流切向分速度和法向分速度均为零. 　　入口边界为集箱径向引入管的入口横截面，入口处采用速度入口条件；密度设置为常数. 　　出口边界条件设置为自由出流边界条件. 　　1.4计算结果及分析 　　本文数值模拟工况：流动介质为空气，温度T=303 K，入口流速ω=62.5 m?s-1，密度ρ=1.165 kg?m-3，动力黏度μ=1.86×10-5 Pa?s.此时入口雷诺数Re=7.81×105，该工况已进入第二自模区. 　　为了达到通用性的目的，对数值模拟中集箱三通区域的静压进行无量纲化处理，将静压分布规律转变成欧拉数Eu的分布规律，即 　　Eu=pi-p0ρω2 　　（3） 　　式中：pi为数值模拟得到的集箱上任一点的静压，Pa；p0为三通径向引入管入口静压，Pa.