小型增压锅炉过热器动态性能仿真研究 谢海涛.docVIP

中国工程热物理年会 传热传质学术会议论文 编号： 小型增压锅炉过热器动态性能仿真研究 谢海涛，李彦军，张国磊，史智俊，郭家敏 (哈尔滨工程大学 动力与能源工程学院 黑龙江 哈尔滨 150001) （TelEmail：xht037030@163.com） 摘要：本文根据热力学和传热学基本原理，结合增压锅炉过热器实际工作过程，建立增压锅炉过热器动态三维数学模型。采用四阶龙格—库塔法对数学方程求解，并编写仿真程序，对多变量扰动下的小型增压锅炉过热器动态过程进行仿真实验，结果表明过热器出口蒸汽温度响应存在延缓，其响应时间远超过蒸汽压力、烟气温度的响应时间。 关键词：增压锅炉；过热器；动态仿真 前言 增压锅炉在舰用蒸汽动力装置中得到广泛应用。与陆用电站锅炉相比，舰用增压锅炉存在频繁、大范围的变负荷操作。同时由于空间限制，船用增压锅炉不设置喷水减温装置。过热器的动态特性决定了增压锅炉变负荷运行性能。因此，对过热器的动态响应性能进行研究，为过热器优化设计、增压锅炉操作运行提供参考，具有重要的实际指导意义。 目前仿真研究在过热器的研究中得到广泛应用。文献[1]、[2]中建立过热器集总参数动态数学模型，研究过热器压力和温度在喷水流量和主汽调门开度扰动下的动态响应。MadejskiPawel、TalerDawid[3]采用数值模拟软件对减温器和吹灰器运行时的过热蒸汽温度以及管壁应力的动态特性进行研究。Ray, Tapan K[4]等在Matlab/Simulink中研究减温器喷水流量扰动对过热系统火用损的影响。上述研究中，只针对单一因素扰动条件下过热器动态响应特性进行研究。实际的过热系统是一个多参数耦合的系统，在船用增压锅炉变负荷运行中，通常是多个参数同时协调变化，这使得上述研究成果无法精确的反映出船用增压锅炉过热器的动态特性。 为了得到增压锅炉过热器在多参数扰动下的动态响应特性，本文建立过热器三维数学模型。采用四阶龙格—库塔法对模型求解，并以小型增压锅炉过热器为例，进行变负荷动态特性仿真研究。 过热器数学型 .1 过热器数学模型建立 根据过热器实际工作特点，可以做以下假设：烟气不发生横向混合，烟气之间没有传热；不考虑并联管的热力和水力不均匀性；忽略辐射换热和蒸汽、管壁的轴向导热；蒸汽在流动方向上参数分布均匀；烟气入口参数均匀。 图1 过热器流程示意图 图1为某小型增压锅炉过热器流程示意图。过热蒸汽被上、下集箱中的档板分成三个流程，管外烟气依次横向冲刷过热器第三、第二、第一流程。 （a）主视图 （b）俯视图 图2过热器离散微元体 过热器管束沿烟气流动方向分成12排，用=1,2,3......,12表示管排序号；沿烟道宽度方向分成3列，用=1,2,3表示管列序号；沿蒸汽流动方向均匀划分成 段，用=1,2,3......,表示。至此，过热器在空间上已完全离散。过热器内部参数可表示为。图2为过热器一个离散的微元体，对该微元体分析可以得到以下方程。 1、蒸汽侧热平衡方程： （1） 2、管壁热平衡方程： （2） 3、烟气侧热平衡方程： （3） 式中：下标、、分别代表蒸汽、烟气和管壁；，温度，；，表面传热系数，W/(m2.K)；，质量流量，kg/s;，沿蒸汽流动方向单位长度质量，kg/m； 、，传热管内、外径，m；，定压比热容，J/(kg.K)。联立方程（1）、（2）、（3）即可得到描述过热器动态过程的微分方程组 微元体内蒸汽压降为： （） 式中：为流动阻力系数。 .2 过热器数学模型求解 为了获得较高的求解精度，本文采用四阶龙格—库塔公式对方程组进行求解，即有： （） 其中， （6） （7） （8） 为了得到更加精确的解，本文在计算过程中根据温度、压力实时调用工质物性计算函数。烟气物性采用自编程序，蒸汽物性根据IAPWS-IF97标准计算。 边界条件及仿真模型验证 仿真边界条件包括：过热蒸汽产量、过热蒸汽入口温度和压力、燃料消耗量、烟气入口温度和压力。其值由小型增压锅炉试验平台运行数据给出，如表1所示。 表1 小型增压锅炉试验平台运行数据 序号 项目 单位 54%负荷 1 油量 kg/h 64.0 81.7 2 蒸汽进口温度 K 447.1 468.8 3 蒸汽进口压力（表压） MPa 0.7012 1.2616 4 蒸汽流量 t/h 0.6627 0.8059 5 烟气进口温度