回转窑预热段温度优化及数值模拟.pdf

doi:10.16018550/n.202302006Journal of Yancheng Institute of Technology(Natural Science Edition)Jun.20232023年0 6 月Vol.36 No.2盐城工学院学（自然科学版）第36 卷第2 期回转窑预热段温度优化及数值模拟金欣玮1，薛长峰,邱海龙，张鸣鑫1.盐城工学院电气工程学院，江苏盐城224051;2.盐城工学院数理学院，江苏盐城224051摘要：为改善回转窑内部温度分布不均的问题，建立了一个基于FLUENT的回转窑预热段三维数值模型，通过正交试验对回转窑主要操作参数进行优化，给出了正常运行工况下回转窑预热段的最佳参数组合方案。结果表明：回转窑进气口高温气体流速3m/s，回转窑转速0.2 r/min，回转窑抬起角度2 时，仿真模拟结果为窑内最低温度10 9 2.2 0 6 K，与窑内最高温度的温差为30.794K，与初始方案相比，窑内最低温度提升了3.35%，温差降低了53.54%，改善了回转窑内温度分布不均的问题关键词：回转窑；数值模拟；正交实验中图分类号：TF355.1文献标志码：A文章编号：16 7 1-532 2(2 0 2 3)0 2-0 0 32-0 5近年来，随着水泥行业工艺设备的不断发展，人们开始广泛关注熟料熳烧流程中的高能耗和高排放问题。可以通过研究回转窑工况对窑内流场、温度场以及煤粉燃烧等因素的影响，为混凝土厂节能减排提供指导意见。由于回转窑内部热传递流程的复杂性，目前的检测技术手段尚不成熟，难以精确获取回转窑内部流场和温度场的分布数据。因此，需要研究和开发更为准确的检测技术，而数值模拟能够研究回转窑内部的流场和温度场分布情况，从而更有效地控制混凝土生产流程中的能耗和排放2-3,对于提高工程效率具有重要意义。国内外学者针对回转窑内部流场、温度场的数值模拟进行了许多研究。Elattar等4 建立了回转窑二维模型，研究了回转窑主要操作参数对不同燃料下火焰行为的影响；Manju等5 将回转窑内熟料看作假流体，对窑内气流和熟料进行三维建模，建立了CFD模型，模拟了固体与固体之间的传热分布情况；Sadighi等6 应用能量守恒和质量守恒定律，给出了一种新的数值模拟方法，通过测量过程变量和扫描壳体温度建立了水泥回转窑烧成带涂层厚度的预测模型；刘彬等7 利用FLUENT建立了回转窑的数值模拟模型，研究了回转窑倾斜角度和不同燃烧器参数对回转窑内部流场分布的影响。上述研究所建立的模型能够较好地模拟回转窑内物料燃烧的特性，但是未能综合考虑回转窑倾斜角度、高温气体进口速度、筒体转速等参数的影响，使得窑内整体温度场的数值模拟结果与实际情况存在较大差异。本文建立回转窑内部流流动、高温气体传热和气固两相的对流换热数学模型，利用有限元方法，应用正交试验，研究气体进口速度、回转窑抬起角度和回转窑转速等对窑内热流场与温度场分布的影响。1回转窑预热段模型针对回转窑内壁温度和整体温度分布不均性，设计回转窑预热段简化模型如图1所示。回转窑预热段简化模型有效长度6 0 0 mm，直径为50mm，内壁厚度为0.4mm，回转窑预热段模型相关物料参数如表1所示。高温气体从回转窑预热段左侧进入，经过内壁与从右侧落人的碳酸钙颗粒进行热交换，预热后的碳酸钙颗粒由左侧出收稿日期：2 0 2 3-0 2-2 2作者简介：金欣玮（19 9 8 一），男，江苏无锡人，硕士生，主要研究方向为回转窑数值模拟。通信作者：薛长峰（19 6 6 一），男，江苏盐城人，教授，博士，主要研究方向为偏微分方程数值解。33金欣玮，等：回转窑预热段温度优化及数值模拟第2 期口流出，达到热量传递的目的。根据实际生产经验将回转窑预热段高温气体进口温度设置为1123K（即窑内最高温度），根据窑内最低温度的高低可以判断回转窑预热段模型传热效果的优劣。最低温度值越高，与窑内最高温度的温差越小，窑内温度分布越均匀。碳酸钙进口Z高温气体进口图1回回转窑预热段简化模型结构示意图Fig.1Structure diagram of simplified model of preheating section of rotary kiln表1回转窑预热段模型相关物料参数Table1Related material parameters of preheatingsection model of rotarykiln密度/比热容/热导率/物料(kgm)(Jkg.K-l)(W mI.K-1)耐火砖320010201.15筒体780043430托轮775048015.1窑皮31008751.16环境空气1.161410070.262回转窑传热数值模拟2.1条件假设实际生产中回转窑内的运动及反应状况复杂，现对模拟过程做出如下假设：（1）模拟流体区域视为稳态流场，即将整个回转窑内部视为运转稳定的结构，运转时其温度、气体流速与物料进口速度等参数无关；（2）回转窑内部密封性完整，所有高温气体由右侧出口排放出回转窑预热段，进入燃烧准备阶段；（3）高温气体温度达不到产生热辐射所需的最低温度；（4）气体为不可压缩流体；（5）不考虑回转窑内的化学反应。2.2控制方程2.2.1质量守恒方程a(pu)a(pu)a(pw)0(1)otdy式中：p为碳酸钙密度，kg/m;t为碳酸钙物料在窑内的停留时间,s;x、y z 为三维坐标;uvuw是坐标所对应的高温气体流速沿X、Y、Z坐标轴方向的分量。2.2.2能量守恒方程a(pT+div(puT)=div(grad T)+Sr(2)at式中：T为窑内温度，K；入为窑内壁热导率，WI（mK）；S表示黏性耗散项。2.2.3对流换热模型高温气体、料层表面、未被物料覆盖的内壁之间的对流换热通过如下方程来表示：Q=hg(T,-T.)A,(3)名+esesQ=hgw(T,T.)A(4)名+ewgS-ewQhWc(T-T.)A(5)cW式中：Q.QQc+c分别表示窑内壁面与物料表面、高温气体与回转窑内壁、高温气体与料床表面的传输热量，J;hehwhcw-cs分别表示窑内壁面与物料表面、高温气体与回转窑内壁、高温气体与料床表面的对流换热系数，I（mK）;T、T、T 分别表示高温气体、窑内壁面和物料的温度，K;AgevAgwAew=c分别表示窑内壁面与物料表面、高温气体与回转窑内壁、高温气体与料床表面的换热面积，m。单位长度内，窑内壁面与物料表面、高温气体与回转窑内壁、高温气体与料床表面的换热面积，由回转窑内径和抬起角度所决定，计算公式分别为：A二Dsin(6)+es2DA(2T-)(7)ew2DA(8)二CWCS234.第36 卷盐城工学院字自然科学版）式中：D,为回转窑内径，m;为回转窑抬起角度，（）2.3初始条件和边界条件数据模拟时选用文献8 中提供的流体边界条件和FLUENT软件中提供的热力学边界条件，选用k-8端流模型9 ,并根据实际生产经验和表1设定相关参数。其中，换热方式为对流换热，传热系数为1.15W/（m.K），来流温度30 0 K，壁面厚度0.0 0 4m，热源功率为0。2.4网格划分为提高回转窑预热段模型计算的准确度和网格质量，利用FLUENT软件自动生成结构性网格，网格单元数为4450 2 7，网格节点数为432936。2.5计算方法利用ANSYSFLUENT2022R1数值软件进行模拟仿真求解，采用有限元方法进行计算，其中涉及的数学模型有控制方程公式（1）公式（8）、k-8模型和离散项模型，选择耦合算法对模型进行求解。模拟过程中，设定高温气体进口为速度进口，高温气体出口为压力出口。在抬起角度为1的情况下，设定回转窑运行时间为30 0 步长，步长为1s，并且每个步长迭代10 次。2.6仿真结果分析将相关参数代人FLUENT求解，可以得到回转窑预热段内部温度场分布如图2 所示。由图2可知，回转窑预热段内温度从左侧高温气体进口开始沿着窑长方向逐渐降低，在右侧碳酸钙进口处温度最低。这是由于从回转窑预热段左侧进入的高温气体（112 3K）与从右侧进入的室温下的碳酸钙物料（30 0 K)逆流接触，高温气体将热量传递给碳酸钙物料和窑内壁，高温气体温度下降，碳酸钙物料和窑内壁温度上升，回转窑预热段内部形成了温度分布。contour-1StaticTemperature1.12e+031.12e+031.12e+031.11e+031-11e+03111e+031.10+031.10e+031.10e+031.10e+031.0.9e+03K图2回转窑温度分布云图Fig.2Nephogram of temperature distribution of rotarykiln3正交试验分析3.1正交试验设计为进一步改善回转窑预热段内物料温度的均匀性，提高回转窑的生产效率和产品品质，综合考虑各因素的影响水平后，选取进气口高温气体的流速、回转窑的转速和回转窑的抬起角度等3个影响因素进行正交试验，以寻找最佳参数组合条件，优化仿真模拟结果。每个因素选取3个水平，高温气体流速（因素A)分别取1、2、3m/s；回转窑转速（因素B)分别取0.2、1.6、3.5r/min；回转窑抬起角度（因素C)分别取1、2、4,正交试验因素水平如表2 所示，其中，A1B1C1方案为生产厂家常用的初始方案。因素ABC水平高温气体流回转窑转速/回转窑抬起速/(ms)(rmin)角度）110.21221.62333.543.2正交试验分析正交试验采用L。（34)正交表进行设计，试验方案及仿真结果数据如表3所示。35金欣玮，等：回转窑预热段温温度优化及数值模拟第2 期表3正交试验方案及仿真结果Table3Orthogonal test scheme and simulation results序号因素A因素B因素C空列最低温度/K11(1)1(0.2)1(1)11056.717212(1.6)2(2)21072.013313(3.5)3(4)31055.58042(2)1231087.053522311059.434623121062.59673(3)1321 092.206832131 078.406933211 089.9171068.3101085.3251 072.5731075.3561 069.6941 069.9511082.9941075.6051086.8431 069.3641 069.0731073.680极差R10.3919.0507.2421.047采用极差法分析正交试验结果，极差越大该因素对窑内最低温度的影响越大。由表3可知，极差R从大到小依次为R(A）、R(B）、R(C),因此，各因素的影响由大到小分别是高温气体流速、回转窑转速和回转窑抬起角度；各因素水平最佳组合为A3B1C2，即高温气体流速3m/s、回转窑转速0.2r/min、回转窑抬起角度2 时，窑内最低温度最高，窑内温差最小。利用方差分析对上述结果进行验证，方差分析结果如表4所示。由表4可知，因素A对窑内最低温度的影响最为显著，因素C对窑内最低温度的影响最不显著，极差分析与方差分析结果一致。表4方差分析表Table 4Tableof varianceanalysis方差型显著方差均方F来源平方和性模型1 453.8976242.31673.663*A647.7832323.89298.462*B491.4682245.73474.702*C314.6462157.32347.826*注：*为 P0.01,*为 P0.05。采用生产厂家初始方案（A1B1C1)时，窑内最低温度为10 56.7 17 K,与窑内最高温度（112 3K）的温差为6 6.2 8 3K；采用正交试验得到的最优方案（A3B1C2)时，窑内最低温度为10 9 2.2 0 6 K，与窑内最高温度的温差为30.7 9 4K；优化后，窑内最低温度提高了3.35%，温差下降了53.54%，有效改善了回转窑预热段温度分布不均的问题。4结论本文采用正交试验进行仿真模拟，针对回转窑预热段传热问题进行研究，得出以下结论：(1)由正交试验结果可以看出，影响回转窑预热段窑内最低温度的最显著因素为回转窑进气口高温气体的流速，其次为回转窑转速，最后是回转窑抬起角度。(2)将优化后的方案与生产厂家初始方案进行比较可以发现，采用优化后的方案，回转窑预热段温度达到稳态时，窑内最低温度为1092.206K，较初始方案提升了3.35%，与最高温度的温差为30.7 9 4K，较初始方案下降了53.54%，改善了回转窑预热段窑内温度分布不均的问题。参考文献：1】李增慧，吴家欢，刘强强新型耐火材料的应用对水泥行业节能减排的影响J.能源研究与管理，2 0 2 2，14（4)：178-182.2】刘彬，周武洲，赵朋程，等水泥回转窑物料流动与传热特性数值模拟J中南大学学报（自然科学版），2 0 19，50责任编辑：熊璐璐36第36 卷盐城工学院学报自然科学版）(5)：12 35-12 43.3范月华，刘恩海，潘嘉信，等水泥熟料生产线余热回收利用系统与节能分析J。节能，2 0 2 0，39(1)：10 1-10 3.4 ELATTAR H F,STANEV R,SPECHT E,et al.CFD simulation of confined non-premixed jet flames in rotary kilns forgaseous fuelsJ.Computers&Fluids,2014,102:62-73.5 MANJU M S,SAVITHRI S.Three dimensional CFD simulation of pneumatic coal injection in a direct reduction rotary kilnJ.Fuel,2012,102:54-64.6 SADIGHI S,SHIRVANI M,AHMAD A.Rotary cement kiln coating estimator:integrated modelling of kiln with shell tem-perature measurementJ.The Canadian Journal of Chemical Engineering,2011,89(1):116-125.7刘彬，周武洲，赵朋程，等。水泥回转窑传热特性数值模拟J.燕山大学学报，2 0 19，43（3）：19 9-2 0 7.8】张志霄，池涌，李水清，等.回转窑传热模型与数值模拟J.化学工程，2 0 0 3,31(4)：2 7-31.9】雷林，王智祥，孙鹏，等。计算流体力学k-8二方程流模型应用研究J.船舶工程，2 0 10，32（3）：5-8.Temperature Optimization and Numerical Simulation of PreheatingSection of Rotary KilnJIN Xinweil,XUE Changfeng,QIU Hailong,ZHANG Mingxin(1.College of Electrical Engineering,Yancheng Institute of Technology,Yancheng Jiangsu224051,China;2.College of Mathematics and Physics,Yancheng Institute of Technology,Yancheng Jiangsu224051,ChinaAbstract:In order to improve the uneven temperature distribution in the rotary kiln,a three-dimensional numerical model of thepreheating section of the rotary kiln based on FLUENT was established.The main operating parameters of the rotary kiln were op-timized through orthogonal tests,and the optimal parameter combination scheme of the preheating section of the rotary kiln undernormal operating conditions was given.The results show that when the flow rate of high temperature gas at the inlet of the rotarykiln is 3 m/s,the rotating speed of the Rotary kiln is 0.2 r/min,and the lifting angle of the rotary kiln is 2,the simulation resultsareas follows:the lowest temperature in the kiln is 1 092.206 K,and the temperature difference from the highest temperature inthe kiln is 30.794 K.Compared with the initial scheme,the lowest temperature in the kiln is increased by 3.35%,and the tem-perature difference is reduced by 53.54%,improving the uneven temperature distribution in the rotary kiln.Keywords:rotarykiln;numerical simulation;orthogonal experiment