CN105468799B - 预报高温废气循环烧结工艺热状态参数的仿真方法 - Google Patents
预报高温废气循环烧结工艺热状态参数的仿真方法 Download PDFInfo
- Publication number
- CN105468799B CN105468799B CN201410446270.5A CN201410446270A CN105468799B CN 105468799 B CN105468799 B CN 105468799B CN 201410446270 A CN201410446270 A CN 201410446270A CN 105468799 B CN105468799 B CN 105468799B
- Authority
- CN
- China
- Prior art keywords
- temperature
- phase
- model
- gas
- sintering
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005245 sintering Methods 0.000 title claims abstract description 135
- 238000000034 method Methods 0.000 title claims abstract description 130
- 230000008569 process Effects 0.000 title claims abstract description 95
- 239000002912 waste gas Substances 0.000 title abstract description 5
- 239000000463 material Substances 0.000 claims abstract description 72
- 239000002245 particle Substances 0.000 claims abstract description 60
- 238000004088 simulation Methods 0.000 claims abstract description 49
- 238000013178 mathematical model Methods 0.000 claims abstract description 41
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 38
- 235000019738 Limestone Nutrition 0.000 claims abstract description 27
- 239000006028 limestone Substances 0.000 claims abstract description 27
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 26
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 24
- 239000000446 fuel Substances 0.000 claims abstract description 5
- 239000007789 gas Substances 0.000 claims description 92
- 239000012071 phase Substances 0.000 claims description 83
- 239000007790 solid phase Substances 0.000 claims description 37
- 238000012546 transfer Methods 0.000 claims description 32
- 238000002485 combustion reaction Methods 0.000 claims description 30
- 238000006243 chemical reaction Methods 0.000 claims description 29
- 238000001704 evaporation Methods 0.000 claims description 27
- 230000008020 evaporation Effects 0.000 claims description 27
- 239000007788 liquid Substances 0.000 claims description 15
- 230000005012 migration Effects 0.000 claims description 15
- 238000013508 migration Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000000694 effects Effects 0.000 claims description 14
- 238000002156 mixing Methods 0.000 claims description 13
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 claims description 12
- 239000000292 calcium oxide Substances 0.000 claims description 12
- ODINCKMPIJJUCX-UHFFFAOYSA-N calcium oxide Inorganic materials [Ca]=O ODINCKMPIJJUCX-UHFFFAOYSA-N 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 12
- 238000004519 manufacturing process Methods 0.000 claims description 12
- 238000009833 condensation Methods 0.000 claims description 11
- 230000005494 condensation Effects 0.000 claims description 11
- 239000007787 solid Substances 0.000 claims description 11
- 238000009826 distribution Methods 0.000 claims description 10
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 9
- 229910052760 oxygen Inorganic materials 0.000 claims description 9
- 239000001301 oxygen Substances 0.000 claims description 9
- 238000004134 energy conservation Methods 0.000 claims description 8
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims description 6
- 229910000019 calcium carbonate Inorganic materials 0.000 claims description 6
- BRPQOXSCLDDYGP-UHFFFAOYSA-N calcium oxide Chemical compound [O-2].[Ca+2] BRPQOXSCLDDYGP-UHFFFAOYSA-N 0.000 claims description 6
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 5
- 239000003546 flue gas Substances 0.000 claims description 5
- 239000007791 liquid phase Substances 0.000 claims description 5
- 238000010438 heat treatment Methods 0.000 claims description 4
- 239000001569 carbon dioxide Substances 0.000 claims description 3
- 229910002092 carbon dioxide Inorganic materials 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 229920006395 saturated elastomer Polymers 0.000 claims description 2
- 230000007246 mechanism Effects 0.000 abstract description 15
- 230000005945 translocation Effects 0.000 abstract 1
- 239000000203 mixture Substances 0.000 description 10
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 8
- 238000011161 development Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 238000001816 cooling Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000009792 diffusion process Methods 0.000 description 4
- 239000004744 fabric Substances 0.000 description 4
- 229910052742 iron Inorganic materials 0.000 description 4
- 239000002994 raw material Substances 0.000 description 4
- 238000005338 heat storage Methods 0.000 description 3
- 230000000977 initiatory effect Effects 0.000 description 3
- 239000000843 powder Substances 0.000 description 3
- 238000004886 process control Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000001035 drying Methods 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 208000031872 Body Remains Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000010924 continuous production Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000005469 granulation Methods 0.000 description 1
- 230000003179 granulation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000004482 other powder Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004449 solid propellant Substances 0.000 description 1
- 239000012798 spherical particle Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
- 238000009736 wetting Methods 0.000 description 1
Landscapes
- Manufacture And Refinement Of Metals (AREA)
Abstract
本发明涉及一种预报高温废气循环烧结工艺热状态参数的仿真方法,本发明仿真方法采用的技术方案是:建立微观上的单颗粒微元体机理模型,并基于此建立高温废气循环烧结工艺的宏观料层数学模型,将水分迁移、石灰石分解、焦粉燃烧等显著影响烧结过程的热效应通过微观机理模型进行求解,并将结算结果加权为宏观模型的质量源项与热量源项,从而在机理上更加贴近实际烧结过程,提高模拟精度。用以解决现有的模拟方式精度低以及用于模拟时间较长、效率低的问题。
Description
技术领域
本发明涉及烧结工艺,具体是一种用于预报高温废气循环烧结工艺热状态参数的数值仿真方法。
背景技术
烧结,是人造富铁矿的主要生产方法之一。目前国际上90%以上的人造富铁矿都是通过烧结得到的,其基本过程为:将准备好的混合料在烧结机上进行点火并通过抽风的作用,使混合料中的固定碳燃烧,从而产生高温,混合物局部软化或熔化,发生一系列的化学反应,生成一定数量的液相。随后由于温度降低,液相冷却而凝固成块,经破碎、冷却即得到烧结矿。
烧结矿质量的好坏直接影响高炉生产的质量。随着烧结设备大型化和高炉对烧结矿质量要求的提高,烧结过程控制技术的应用就显得更加重要,而利用数值计算方法对烧结过程进行数值模拟并配有实验验证,对准确掌握烧结过程的传热规律,提高烧结矿产量和质量有重大的现实意义。
带式烧结机的混合料层的在抽风烧结过程中,整个料层沿高度呈现性质不同的5个带,从上到下分别是烧结矿带、燃烧带、加热带、干燥带和过湿带,各种复杂的物理化学反应主要发生在燃烧带,温度最高。在燃烧带以上的烧结矿带内,炽热的烧结矿物料加热时从料层表面吸入空气,燃烧带以下的3带统称未烧结矿层,在这3个带中,来自燃烧带的高温烟气将热量传给物料,使物料得到干燥和预热升温。
烧结矿料层的传热传质模型始于上世纪60年代。1963年,英国Leeds大学的P.A.Young(烧结矿生产,冶金工业出版社,1981)把烧结矿料层当做固定床,简单计算了加热气流与烧结矿之间的传热过程;1967年,M.Cross(粉末烧结理论,冶金工业出版社,1998)建立了回转窑中气流与烧结矿颗粒之间的传热模型,但没有考虑各种反应对传热过程的影响;H.Toda等人(Transactions ISIJ,1984)建立了模拟焦粉粒度、混合料预热、气流分布等操作因素对烧结料层热曲线影响的数学模型,得到了高温带冷却时间与落下强度之间的关系;中南大学的邹志毅(中南工业大学博士论文,1992)建立了一个以湿料带传热、传质机理为基础的烧结水分迁移全过程模型,得到料层中的温度场与水分迁移过程有关的所有变量的动态响应;东北大学力杰(东北大学学报,2010)应用Fluent软件,采用多孔介质模型,分析了进口流速、料层高度等操作参数以及风箱结构对冷却过程的影响;中南大学龙红明(中南工业大学博士论文,2007)提出了将烧结过程热状态描述为料层各带厚度和迁移速度的思想,在烧结层上建立了三维模型,采用三对角矩阵算法和欠松弛迭代算法进行了数值模拟,并通过模拟解和实测解的对比分析,论证所建立的模型;重启大学的白晨光(重庆大学学报,2008)建立了烧结料层的蓄热模型,研究了燃烧层厚度、蓄热量在料层中的分配比对烧结蓄热的影响;2009年,北京科技大学的司俊龙(中南工业大学学位论文,2009)将烧结料层视为非均质多孔介质,忽略烧结机长度及宽度方向上的传输过程,得到一维非稳态传热过程,简要分析了料层厚度、台车速度、混料预热温度、分层布料等对烧结过程的影响。
上述烧结模拟模型主要用于离线操作分析,在工业应用上尚有限制,主要缺点如下:(1)模型中一些主要技术参数难以确定,因此其模拟的计算精度不能满足工业过程控制的要求;(2)综合数学模型的控制方程组相当复杂,而且又不能用线性方程式来表示,模型求解时间过长;(3)对系统的输出变量,即质量、产量指标和过程状态变量之间的定量关系尚未用数学模型描述出来;(4)适用范围比较窄,不能对烧结过程做出超过已知测定变量或实验变量范围以外的预报。
由上述分析可见,由于烧结过程复杂,影响因素繁多,不易建立合适的数学模型。但纵观烧结模型的发展历史,它经历着由简单到复杂,由低级到高级的发展过程。在模型的范围方面,由选球、制粒模型,烧结矿冷却模型向烧结模型发展,由窄范围向宽范围发展。在模型研制水平方面,由经验模型向综合模型发展,由静态模型向动态模型发展,由单元操作向系统模型发展。在模型运用方面,由用于过程分析、预测向用于模拟、控制的方向发展。对烧结数学模型的发展,国外始于20世纪60年代,我国始于80年代。由于烧结技术和计算机技术的发展,烧结数学模型的发展很快,应用范围不断扩大,它从过程模拟、参数优化不断向过程控制和新技术开发等方向发展。
发明内容
本发明的目的是提供一种预报高温废气循环烧结工艺热状态参数的仿真方法,本发明将水分迁移、石灰石分解、焦粉燃烧等显著影响烧结过程的热效应通过微观机理模型进行求解,并将结算结果加权为宏观模型的质量源项与热量源项,在机理上更加贴近实际烧结过程,提高模拟精度。用以解决现有的模拟方式精度低以及用于模拟时间较长、效率低的问题。
为实现上述目的,本发明的方案是:一种预报高温废气循环烧结工艺热状态参数的仿真方法,所述的仿真方法包括如下步骤:
(1)建立烧结机内烧结料层的烧结过程物理模型,所述的物理模型是一种一维非稳态模型,用以模拟烧结燃料层的物料与气体在烧结机内交叉换热时进行的流动、传热、燃烧和化学反应过程;
(2)建立单颗粒微元体微观数学模型;
(3)根据所述的单颗粒微元体微观数学模型,求解烧结过程中水分迁移、石灰石分解和焦粉燃烧的热效应,并将求解结果加权为烧结层宏观数学模型的质量源项和热量源项,建立烧结床层宏观数学模型;
(4)根据所述的微观数学模型和宏观数学模型,计算得到单颗粒微元体内部温度与浓度梯度,混料内部温度与气体温度,尾部烟气温度,模拟烧结过程不同工艺参数和控制参数下的混料温度变化过程;
(5)根据步骤(4)中的仿真结果,对实际烧结过程中的控制参数和工艺参数进行调整,优化烧结机生产工艺。
根据本发明所述的仿真方法,所述的步骤(3)中,求解水分蒸发的热效应时,建立
水分迁移子模型,所述的水分迁移子模型包括水分冷凝速度模型和水分蒸发速度模型,所
述的水分冷凝速度模型为:所述的水分蒸发速度模型为:
其中:Revap1为水分冷凝速率,kg/(m3·s);
Revap2为水分蒸发速率,kg/(m3·s);
为气相中水蒸气密度,kg/m3;
为饱和水蒸气密度,kg/m3;
△τ为模拟计算的时间步长,s;
Tevap为水分蒸发临界温度,K;
△Hevap为水分蒸发焓变,J/kg。
根据本发明所述的仿真方法,所述的步骤(3)中,求解石灰石分解的热效应时,建立石灰石分解子模型,即:
其中:Rdeco为石灰石分解速率,kg/(m3·s);
R为通用气体常数,取值8.314J/(mol·K);
和mCaO分别为固相中碳酸钙和氧化钙的质量,kg/m3;
和MCaO分别为碳酸钙和氧化钙的摩尔质量,kg/mol;
为气相中二氧化碳的分压,Pa;
Keq为反应平衡常数,通过求解。
根据本发明所述的仿真方法,所述的步骤(3)中,求解焦粉燃烧的热效应时,建立焦粉燃烧子模型,即:
式中:Rcomb为焦粉燃烧速度;
为气相中氧气分密度,kg/m3。
根据本发明所述的仿真方法,当单颗粒微元体的温度小于水分蒸发温度时,若气相中水蒸气密度超过饱和水蒸气密度,则根据所述的水分冷凝速度模型进行水分蒸发热的处理;当检测到单颗粒微元体的实际温度大于水分蒸发温度时,根据所述的水分蒸发速度模型进行水分蒸发热的处理.
根据本发明所述的仿真方法,当检测到单颗粒微元体的实际温度达到分解温度时,发生石灰石分解,根据所述的石灰石分解子模型进行石灰石分解热的处理;
根据本发明所述的仿真方法,当检测到单颗粒微元体的实际温度达到燃烧温度时,发生燃烧反应,根据所述的焦粉燃烧子模型进行焦粉燃烧热的处理。
根据本发明所述的仿真方法,所述的工艺参数和控制参数包括布料厚度、台车速度、混料预热温度、进口气体流量/温度/氧含量、混料水分/石灰石/焦粉含量。
根据本发明所述的仿真方法,所述单颗粒微元体微观数学模型的控制方程如下:
微观固相质量守恒方程:
微观气相质量守恒方程:
微观能量守恒方程:
定解条件如下:
其中,
ρCp=(1-φ)ρslCsl+φρgCg;
λ=(1-φ)λsl+φλg;
φ为多孔介质的孔隙率,无量纲;
ρ、ρsl、ρg分别为气液固三相、固液两相和气相的密度,kg/m3;
τ为时间,s;
Msl、Mg分别为反应带来的固液两相和气相质量的变化,kg/(m3·s);
r为微观单颗粒微元体的半径,m;
Deff为气相的有效质量扩散率,m2/s;
Cp、Csl、Cg分别为气液固三相、固液两相和气相的比热,J/(kg·K);
T为温度,K;
λ、λsl、λg分别为气液固三相、固液两相和气相的热导率,W/(m·K);
△qv为反应带来的热量变化,W/m3;
βm为微观单颗粒微元体表面的传质系数,m/s;
ρg,∞为微观单颗粒微元体表面来流的密度,kg/m3;
hsa为对流换热系数,W/(m2·K);
Asa为比表面积,m2/m3;
Tg,∞为微观单颗粒微元体表面来流的温度,K;
εm为黑度,无量纲;
σ为Stefan-Bolzmann常数,W/(m2·K4);
△Hk为化学反应焓变,J/kg。
根据本发明所述的仿真方法,所述的烧结床层宏观数学模型的控制方程如下:
宏观气相质量守恒方程为:
宏观固相能量守恒方程为:
宏观气相能量守恒方程为:
定解条件如下:
τ=0时:
y=0时:Tg=Tg,in,ug=ug,in,ρg=ρg,in,
且在烧结不同阶段,Tg,in、ug,in、ρg,in具有不同的取值;
式中:
ug,in、ug分别为宏观气相进入宏观料层上表面和流过料层内部的速度,m/s;
y为料层高度,m;
ρs为宏观固相的密度,kg/m3;
Cps、Cpg分别为宏观固相和宏观气相的比热,J/(kg·K);
λs为宏观固相的热导率,W/(m·K);
Ts、Tg分别为宏观固相和宏观气相的温度,K;
为化学反应产生热量被宏观固相吸收的比率,无量纲;
Tg,in为宏观气相进入宏观料层上表面的温度,K;
ρg,in为宏观气相进入宏观料层上表面的密度,kg/m3;
本发明达到的有益效果:本发明通过建立微观上的单颗粒微元体机理模型,并基于此建立高温废气循环烧结工艺的宏观料层数学模型,将水分迁移、石灰石分解、焦粉燃烧等显著影响烧结过程的热效应通过微观机理模型进行求解,并将结算结果加权为宏观模型的质量源项与热量源项,从而在机理上更加贴近实际烧结过程,提高模拟精度。
利用数值计算方法对烧结过程进行数值模拟并实验验证,对准确掌握烧结过程的传热规律,提高烧结矿的产量和质量有重大的支撑作用。本发明的数学模型是在综合全面地考虑了烧结机内流动/传热/燃烧和化学反应等诸多热现象的基础上建立的,利用该模型可以对不同布料厚度、台车速度、混料预热温度、进口气体流量/温度/氧含量、混料水分/石灰石/焦粉含量等影响烧结过程的主要因素进行数值模拟计算和分析,从而模拟不同生产工况下的烧结过程,通过模拟计算得到烧结终点的位置,提高烧结矿产量和质量,优化烧结实际生产的过程控制。
附图说明
图1是本发明的方法流程图;
图2从实际燃料床到一维非稳态模型的转换示意图;
图3本发明的小球微元体机理模型求解过程示意图;
图4烧结床层宏观模型求解过程示意图;
图5应用本发明方法模拟的不同料层高度下的固相温度变化曲线;
图6应用本发明方法模拟的不同台车速度下的固相温度变化曲线;
图7应用本发明方法模拟的不同混料预热温度下的固相温度变化曲线。
具体实施方式
下面结合附图和实施例对本发明作进一步详细的说明。
本发明首先建立烧结机内流动/传热/燃烧和化学反应过程物理模型,在对其进行必要的简化假设的基础上,建立“烧结机内流动/传热/燃烧和化学反应过程数学模型”,最后通过实际测量的烧结饼内各层温度来对所建立的数学模型进行验证与调整,全面的分析并优化实际烧结过程中的控制参数与工艺参数。
以往的烧结料层模型,均是建立宏观意义上的数学模型,对于模型中反应热量的处理过于粗糙,因而计算准度稍显不足。针对这个问题,本发明拟建立微观上的单颗粒微元体机理模型,并基于此建立高温废气循环烧结工艺的宏观料层数学模型,将水分迁移、石灰石分解、焦粉燃烧等显著影响烧结过程的热效应通过微观机理模型进行求解,并将结算结果加权为宏观模型的质量源项与热量源项。从而在机理上更加贴近实际烧结过程,提高模拟精度。
如图1所示,本发明的具体实施过程如下:
步骤1,建立烧结机内烧结料层的烧结过程物理模型:
在烧结过程中,烧结料层经历了一个从上到下逐层烧结的过程,同时由于烧结是连续生产过程,台车移动速度很慢,且在布料均匀的情况下,宽度方向上温度较均匀,可以不考虑台车运动方向和台车宽度方向上的传输过程,所以,烧结可以简化为料层厚度方向上的层状传输过程。
图2为从实际烧结燃料床到一维非稳态模型的转换示意图。在烧结过程中,气体从上向下流动,物料水平移动。气体与物料交叉流动换热,其物理模型如下:具有相同温度Ts,高度为H0的物料以ws的速度与具有进口温度Tg,in的气体进行热交换。
步骤2,建立单颗粒微元体微观数学模型:
烧结料层在宏观上是由铁矿粉、熔剂、固体燃料等粉末堆积而成的具有一定孔隙率和当量粒径的散料层,在烧结过程中随着高温带的下移分成烧结矿带、燃烧带、预热带、干燥带和过湿带五个区,各个区域的空隙率和粒径等结构参数具有明显差异。微观上,烧结混料在经过混合制粒过程中,各种粉末均匀混合并相互粘结后形成了具有一定粒径分布的微元体颗粒,在烧结过程中该微元体内部发生强烈的物理化学反应,并通过微元体颗粒之间的传热传质实现了热量的传递。
单颗粒微元体微观数学模型控制方程如下:
(1)质量守恒方程
固相(包括液相):
气相:
(2)能量守恒方程
其中,ρCp=(1-φ)ρslCsl+φρgCg;λ=(1-φ)λsl+φλg。
(3)定解条件
将宏观料层均匀的划分为n等分,每个等分就是一个控制体。假设料层宏观模型的控制体内的小球在每一个宏观计算步长内的环境条件(状态)保持不变,即气相温度、流速、密度、压力、热导率、动力粘度、普朗特数均保持不变,并考虑局部热平衡,得到定解条件如下:
式中:
φ为多孔介质的孔隙率,无量纲;
ρ、ρsl、ρg分别为气液固三相、固液两相和气相的密度,kg/m3;
τ为时间,s;
Msl、Mg分别为反应带来的固液两相和气相质量的变化,kg/(m3·s);
r为微观单颗粒微元体的半径,m;
Deff为气相的有效质量扩散率,m2/s;
Cp、Csl、Cg分别为气液固三相、固液两相和气相的比热,J/(kg·K);
T为温度,K;
λ、λsl、λg分别为气液固三相、固液两相和气相的热导率,W/(m·K);
△qv为反应带来的热量变化,W/m3;
βm为微观单颗粒微元体表面的传质系数,m/s;
ρg,∞为微观单颗粒微元体表面来流的密度,kg/m3;
hsa为对流换热系数,W/(m2·K);
Asa为比表面积,m2/m3;
Tg,∞为微观单颗粒微元体表面来流的温度,K;
εm为黑度,无量纲;
σ为Stefan-Bolzmann常数,W/(m2·K4);
△Hk为化学反应焓变,J/kg。
单颗粒微元体微观数学模型的求解过程如图3所示。输入的参数包括基本的设备参数、工艺参数、控制参数以及物料的物性参数、几何参数,生成网格就是将宏观料层均匀的划分为n等分。Pt是一个用来表示模拟迭代了多少次的量,迭代一次,Pt自加1。建立所述的单颗粒微元体微观数学模型的简化假设条件如下:
1)虚拟单颗粒微元体主要成分为细铁矿、石灰石、焦炭粉末,忽略其他成分,且各成份混合均匀,具有一定孔隙率;
2)微观上,虚拟单颗粒微元体在一个稳定的环境中受热升温、发生化学反应;
3)虚拟单颗粒微元体内部各化学反应是相对独立的,忽略其相互影响;
4)虚拟单颗粒微元体为标准球形颗粒,认为周向温度和成分均匀,建立径向一维模型;
5)虚拟单颗粒微元体内部热效应以热传导为主,气固间满足局部热平衡;
6)气体在虚拟单颗粒微元体内部传递以扩散为主;
7)在反应过程中虚拟单颗粒微元体的形状保持不变,颗粒密度逐渐减小。
步骤3,建立烧结床层宏观数学模型:
将水分迁移、石灰石分解、焦粉燃烧等显著影响烧结过程的热效应通过单颗粒微元体微观数学模型进行求解,并将求解结果加权为宏观模型的质量源项与热量源项,建立烧结床层宏观数学模型。烧结床层宏观数学模型控制方程如下:
(1)气相质量守恒方程
(2)能量守恒方程
固相:
气相:
(3)定解条件
τ=0时:
y=0时:Tg=Tg,in,ug=ug,in,ρg=ρg,in,
且在烧结不同阶段,Tg,in、ug,in、ρg,in具有不同的取值。
式中:
ug,in、ug分别为宏观气相进入宏观料层上表面和流过料层内部的速度,m/s;
y为料层高度,m;
ρs为宏观固相的密度,kg/m3;
Cps、Cpg分别为宏观固相和宏观气相的比热,J/(kg·K);
λs为宏观固相的热导率,W/(m·K);
Ts、Tg分别为宏观固相和宏观气相的温度,K;
为化学反应产生热量被宏观固相吸收的比率,无量纲;
Tg,in为宏观气相进入宏观料层上表面的温度,K;
ρg,in为宏观气相进入宏观料层上表面的密度,kg/m3。
烧结床层宏观数学模型的求解过程如图4所示。烧结过程在烧结不断阶段,要先后经过点火、保温、吸风烧结等,这些阶段料层上表面的风温、风速、氧浓度等边界条件是不一样的,即不同的时间,混料随烧结机台车到达不同的位置,就需要读取这个位置处的边界条件。t指时间,t=0表示初始时刻。Pt表示模拟迭代的次数。dt和dt-2都是代表模拟计算的时间步长,dt是宏观模型的,数量级是秒级,dt-2是微观模型的,数量级为毫秒级或微秒级。
结合微元体机理模型的求解,对料层物料换热的计算和分析过程进行了如下简化:
1)除了固相中的水分蒸发外,挥发份的生成、逸出和分解,固相的软化、熔化和再硬化等过程均包括在比热和体积密度等值随温度的变化中;
2)对流传热占绝对主导地位,气相与固相的热传导可以忽略;
3)认为固相粒子内部的热传导足够强,忽略颗粒内部的温差;
4)气体与物料料层的运动在断面上是均匀的,不存在偏流情况;
5)放置料层装置的墙是非渗透的和绝热的。
步骤4,根据所述的微观数学模型和宏观数学模型,计算得到单颗粒微元体内部温度与浓度梯度,混料内部温度与气体温度,尾部烟气温度,模拟烧结过程不同工艺参数和控制参数下的混料温度变化过程;
步骤5,根据实际测量的烧结过程混料温度,对步骤,4中的模拟温度进行调整,使得模拟温度与实际温度相吻合,从而根据模拟结果对实际烧结过程中的控制参数和工艺参数,如布料厚度、台车速度、混料预热温度、进口气体流量/温度/氧含量、混料水分/石灰石/焦粉含量进行调整,优化烧结机生产工艺。
本发明中,所述的水分迁移子模型包括水分冷凝速度模型和水分蒸发速度模型,所述的水分冷凝速度模型为:所述的水分蒸发速度模型为:
其中:Revap1为水分冷凝速率,kg/(m3·s);
Revap2为水分蒸发速率,kg/(m3·s);
为气相中水蒸气密度,kg/m3;
为饱和水蒸气密度,kg/m3;
△τ为模拟计算的时间步长,s;
Tevap为水分蒸发临界温度,K;
△Hevap为水分蒸发焓变,J/kg。
石灰石分解子模型为:
式中:
Rdeco为石灰石分解速率,kg/(m3·s);
R为通用气体常数,取值8.314J/(mol·K);
和mCaO分别为固相中碳酸钙和氧化钙的质量,kg/m3;
和MCaO分别为碳酸钙和氧化钙的摩尔质量,kg/mol;
为气相中二氧化碳的分压,Pa;
Keq为反应平衡常数,通过求解。
焦粉燃烧子模型为:
式中:为气相中氧气分密度,kg/m3。
本发明在建立单颗粒微元体微观数学模型时,需要对一些关键参数进行确定,即:
(1)孔隙率。在烧结前后,固相介质从烧结原料变成烧结矿,原料和矿体具有不同的空隙结构,其取值不同。孔隙率根据前苏联学者萨斯科夫与洛别斯卡娅提出的经验公式进行计算:
φ=0.0005·(15.5a1+11.3a2+9.1a3+7.6a4+6.7a5+6.3a6)
式中:a1~a6为筛级组成(%),相应为颗粒当量直径分别为>80、60-80、40-60、25-40、10-25、0-10mm。
(2)热物性参数。在烧结前后,固体介质从烧结原料变成烧结矿,组成烧结原料和烧结矿的化学物质含量各不相同,其热物性参数也不尽相同,本系统通过测定或计算烧结混料和烧结矿的化学组成(即百分比),再通过各化学组分所占质量百分比加权求和计算各自的物性参数。气相的热物性参数类似于固相。
(3)颗粒表面对流换热系数。对流换热系数通过求解努赛尔数计算。对于本项目中的微元体颗粒介质对流换热过程而言,相当于小球绕流换热,努赛尔数可用下式确定:
式中:Nu为努赛尔数;
Re为雷诺数;
Pr为普朗特数;
μ∞和μs分别为环境来流气体和颗粒表面气体的动力粘度,kg/(m·s)。
(4)颗粒内部气相有效扩散系数。气相组分在微元体内部的有效扩散系数通过下式确定:
式中:Deff,i为气相有效扩散系数,m2/s;
Di为气相组分的扩散系数,m2/s。
μ为多孔介质内部空隙曲折因子,取值
(5)颗粒表面气相传质系数。传质系数通过舍伍德数修正计算,舍伍德数可通过下式计算:
Sh=2+0.69Re0.5Sc0.333
式中:Sh为舍伍德数;
Sc为斯密特准数。
按照本发明的方法,在保持其它参数不变的情况下,模拟不同布料厚度、台车速度、混料预热温度下的固相温度变化曲线,模拟结果如图5~图7所示。图5~图7中五条曲线分别代表在高度方向上距离台车篦条0.03m、0.18m、0.28m、0.38m、0.48m处固相温度随台车位置的变化。图5a模拟的布料厚度为0.6m,图5b模拟的布料厚度为0.68m,图5c模拟的布料厚度为0.76m。图6a模拟的台车速度为1.0m/min,图6b模拟的台车速度为1.24m/min,图6c模拟的台车速度为1.5m/min。图7a模拟的混料初始温度为30℃,图7b模拟的混料初始温度为50℃,图7c模拟的混料初始温度为80℃。通过与实际烧结过程进行对比,该模拟结果与烧结过程实际生产符合程度高,达到当前在线控制的精度要求。
本发明提出的“基于层状传输理论的烧结热过程数学模型”是在综合全面地考虑了烧结机内流动/传热/燃烧和化学反应等诸多热现象的基础上建立的。该模型以计算混料内部温度场为基准,以混料底部温度达到温度峰值的时刻作为烧结终点,利用该模型可以模拟不同生产工况下的烧结过程,即对不同布料厚度、台车速度、混料预热温度、进口气体流量/温度/氧含量、混料水分/石灰石/焦粉含量等影响烧结过程的工艺参数、控制参数进行数值模拟计算和分析,从而更好实现烧结机生产过程实时控制。
Claims (10)
1.一种预报高温废气循环烧结工艺热状态参数的仿真方法,其特征在于:所述的仿真方法包括如下步骤:
(1)建立烧结机内烧结料层的烧结过程物理模型,所述的物理模型是一种一维非稳态模型,用以模拟烧结燃料层的物料与气体在烧结机内交叉换热时进行的流动、传热、燃烧和化学反应过程;
(2)建立单颗粒微元体微观数学模型;
(3)根据所述的单颗粒微元体微观数学模型,求解烧结过程中水分迁移、石灰石分解和焦粉燃烧的热效应,并将求解结果加权为烧结层宏观数学模型的质量源项和热量源项,建立烧结床层宏观数学模型;
(4)根据所述的微观数学模型和宏观数学模型,计算得到单颗粒微元体内部温度与浓度梯度,混料内部温度与气体温度,尾部烟气温度,模拟烧结过程不同工艺参数和控制参数下的混料温度变化过程;
(5)根据步骤(4)中的仿真结果,对实际烧结过程中的控制参数和工艺参数进行调整,优化烧结机生产工艺。
2.根据权利要求1所述的预报高温废气循环烧结工艺热状态参数的仿真方法,其特征在于所述的步骤(3)中,求解水分迁移的热效应时,建立水分迁移子模型,所述的水分迁移子模型包括水分冷凝速度模型和水分蒸发速度模型,所述的水分冷凝速度模型为:所述的水分蒸发速度模型为:
其中:Revap1为水分冷凝速率,kg/(m3·s);
Revap2为水分蒸发速率,kg/(m3·s);
为气相中水蒸气密度,kg/m3;
为饱和水蒸气密度,kg/m3;
Δτ为模拟计算的时间步长,s;
Tevap为水分蒸发临界温度,K;
ΔHevap为水分蒸发焓变,J/kg;
Cp为气液固三相的比热;
ρ为气液固三相的密度;
Tsl为固液两相的温度。
3.根据权利要求1所述的预报高温废气循环烧结工艺热状态参数的仿真方法,其特征在于所述的步骤(3)中,求解石灰石分解的热效应时,建立石灰石分解子模型,即:
其中:Rdeco为石灰石分解速率,kg/(m3·s);
R为通用气体常数,取值8.314J/(mol·K);
和mCaO分别为固相中碳酸钙和氧化钙的质量,kg/m3;
和MCaO分别为碳酸钙和氧化钙的摩尔质量,kg/mol;
为气相中二氧化碳的分压,Pa;
Keq为反应平衡常数,通过求解;
Tsl为固液两相的温度。
4.根据权利要求1所述的预报高温废气循环烧结工艺热状态参数的仿真方法,其特征在于所述的步骤(3)中,求解焦粉燃烧的热效应时,建立焦粉燃烧子模型,即:
式中:Rcomb为焦粉燃烧速度;
为气相中氧气分密度,kg/m3;
Tsl为固液两相的温度。
5.根据权利要求2所述的预报高温废气循环烧结工艺热状态参数的仿真方法,其特征在于当单颗粒微元体的温度小于水分蒸发温度时,若气相中水蒸气密度超过饱和水蒸气密度,则根据所述的水分冷凝速度模型进行水分蒸发热的处理;当检测到单颗粒微元体的实际温度大于水分蒸发温度时,根据所述的水分蒸发速度模型进行水分蒸发热的处理。
6.根据权利要求3所述的预报高温废气循环烧结工艺热状态参数的仿真方法,其特征在于当检测到单颗粒微元体的实际温度达到分解温度时,发生石灰石分解,根据所述的石灰石分解子模型进行石灰石分解热的处理。
7.根据权利要求4所述的预报高温废气循环烧结工艺热状态参数的仿真方法,其特征在于当检测到单颗粒微元体的实际温度达到燃烧温度时,发生燃烧反应,根据所述的焦粉燃烧子模型进行焦粉燃烧热的处理。
8.根据权利要求1所述的预报高温废气循环烧结工艺热状态参数的仿真方法,其特征在于所述的工艺参数和控制参数包括布料厚度、台车速度、混料预热温度、进口气体流量/温度/氧含量、混料水分/石灰石/焦粉含量。
9.根据权利要求1所述的预报高温废气循环烧结工艺热状态参数的仿真方法,其特征在于所述单颗粒微元体微观数学模型的控制方程如下:
微观固相质量守恒方程:
微观气相质量守恒方程:
微观能量守恒方程:
定解条件如下:
其中,
ρCp=(1-φ)ρslCsl+φρgCg;
λ=(1-φ)λsl+φλg;
φ为多孔介质的孔隙率,无量纲;
ρ、ρsl、ρg分别为气液固三相、固液两相和气相的密度,kg/m3;
τ为时间,s;
Msl、Mg分别为反应带来的固液两相和气相质量的变化,kg/(m3·s);
r为微观单颗粒微元体的半径,m;
Deff为气相的有效质量扩散率,m2/s;
Cp、Csl、Cg分别为气液固三相、固液两相和气相的比热,J/(kg·K);
T为温度,K;
λ、λsl、λg分别为气液固三相、固液两相和气相的热导率,W/(m·K);
Δqv为反应带来的热量变化,W/m3;
βm为微观单颗粒微元体表面的传质系数,m/s;
ρg,∞为微观单颗粒微元体表面来流的密度,kg/m3;
hsa为对流换热系数,W/(m2·K);
Asa为比表面积,m2/m3;
Tg,∞为微观单颗粒微元体表面来流的温度,K;
εm为黑度,无量纲;
σ为Stefan-Bolzmann常数,W/(m2·K4);
ΔHk为化学反应焓变,J/kg。
10.根据权利要求1所述的预报高温废气循环烧结工艺热状态参数的仿真方法,其特征在于所述的烧结床层宏观数学模型的控制方程如下:
宏观气相质量守恒方程为:
宏观固相能量守恒方程为:
宏观气相能量守恒方程为:
定解条件如下:
τ=0时:
y=0时:Tg=Tg,in,ug=ug,in,ρg=ρg,in,
;
且在烧结不同阶段,Tg,in、ug,in、ρg,in具有不同的取值;
式中:
ug,in、ug分别为宏观气相进入宏观料层上表面和流过料层内部的速度,m/s;
y为料层高度,m;
ρs为宏观固相的密度,kg/m3;
Cps、Cpg分别为宏观固相和宏观气相的比热,J/(kg·K);
λs为宏观固相的热导率,W/(m·K);
Ts、Tg分别为宏观固相和宏观气相的温度,K;
为化学反应产生热量被宏观固相吸收的比率,无量纲;
Tg,in为宏观气相进入宏观料层上表面的温度,K;
ρg,in为宏观气相进入宏观料层上表面的密度,kg/m3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410446270.5A CN105468799B (zh) | 2014-09-03 | 2014-09-03 | 预报高温废气循环烧结工艺热状态参数的仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410446270.5A CN105468799B (zh) | 2014-09-03 | 2014-09-03 | 预报高温废气循环烧结工艺热状态参数的仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105468799A CN105468799A (zh) | 2016-04-06 |
CN105468799B true CN105468799B (zh) | 2018-10-02 |
Family
ID=55606498
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410446270.5A Active CN105468799B (zh) | 2014-09-03 | 2014-09-03 | 预报高温废气循环烧结工艺热状态参数的仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105468799B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106548470B (zh) * | 2016-10-19 | 2019-04-12 | 大连理工大学 | 一种用ct计算多孔介质内气液间舍伍德数经验公式的方法 |
CN106834662B (zh) * | 2017-01-19 | 2018-09-18 | 中国地质大学(武汉) | 一种基于多工况烧结过程的co/co2比值预测方法 |
CN107341289B (zh) * | 2017-06-05 | 2020-05-29 | 中冶华天工程技术有限公司 | 描述铁矿石烧结床层空隙度变化的计算方法 |
CN108460508B (zh) * | 2017-11-06 | 2021-07-30 | 山西大学 | 基于循环流化床锅炉工作特性在线节能监测方法 |
CN109063320B (zh) * | 2018-07-27 | 2023-07-21 | 江苏大学 | 一种预测链篦机鼓风干燥段球团干燥过程的数值计算方法 |
CN112130454B (zh) * | 2020-08-17 | 2021-05-14 | 广西柳钢东信科技有限公司 | 一种烧结混合料水分智能调节方法及系统 |
CN115034141B (zh) * | 2022-06-28 | 2024-07-19 | 南京航空航天大学 | 一种磨削离散热源的自适应与最优化温度场建模方法 |
CN115938512A (zh) * | 2022-11-26 | 2023-04-07 | 宝钢工程技术集团有限公司 | 一种运用软件算法优化烘烤炉rto工艺参数状态的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1970357A1 (en) * | 2007-03-14 | 2008-09-17 | Colorobbia Italia S.p.a. | Materials for coating ceramic bodies, preparation and use thereof and ceramic bodies comprising them |
CN101893384A (zh) * | 2009-05-21 | 2010-11-24 | 宝山钢铁股份有限公司 | 烧结低温余热循环及排放废气减量的方法及其装置 |
-
2014
- 2014-09-03 CN CN201410446270.5A patent/CN105468799B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1970357A1 (en) * | 2007-03-14 | 2008-09-17 | Colorobbia Italia S.p.a. | Materials for coating ceramic bodies, preparation and use thereof and ceramic bodies comprising them |
CN101893384A (zh) * | 2009-05-21 | 2010-11-24 | 宝山钢铁股份有限公司 | 烧结低温余热循环及排放废气减量的方法及其装置 |
Non-Patent Citations (2)
Title |
---|
Micro-macro simulation of sintering process by coupling Monte Carlo and finite element methods;K.Mori等;《International Journal of Mechanical Sciences》;20040630;第46卷(第6期);第841-854页 * |
烧结冷却废气余热回收循环系统设计;王为术 等;《华北水利水电学院学报》;20091231;第30卷(第6期);第67-70页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105468799A (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105468799B (zh) | 预报高温废气循环烧结工艺热状态参数的仿真方法 | |
Yang et al. | Numerical modeling of reaction and flow characteristics in a blast furnace with consideration of layered burden | |
Bluhm-Drenhaus et al. | A coupled fluid dynamic-discrete element simulation of heat and mass transfer in a lime shaft kiln | |
Da Costa et al. | Modelling a new, low CO2 emissions, hydrogen steelmaking process | |
Shen et al. | Three-dimensional modeling of flow and thermochemical behavior in a blast furnace | |
Barati | Dynamic simulation of pellet induration process in straight-grate system | |
CN109063320B (zh) | 一种预测链篦机鼓风干燥段球团干燥过程的数值计算方法 | |
CN104480300B (zh) | 一种基于预测回转窑内球团矿抗压强度的球团生产方法 | |
Mujumdar et al. | CFD modeling of rotary cement kilns | |
Chang et al. | Modeling on dry centrifugal granulation process of molten blast furnace slag | |
JP3033466B2 (ja) | 高炉の操業方法 | |
Zhao et al. | Modelling and analysis of the combustion behaviour of granulated fuel particles in iron ore sintering | |
Fan et al. | Predictive models and operation guidance system for iron ore pellet induration in traveling grate–rotary kiln process | |
Jiao et al. | CFD modeling and analysis of particle size reduction and its effect on blast furnace ironmaking | |
Copertaro et al. | A discrete-continuous approach to describe CaCO3 decarbonation in non-steady thermal conditions | |
de Paula Vitoretti et al. | Study of the induration phenomena in single pellet to traveling grate furnace | |
Xu et al. | Numerical investigation of the reacting flows of three adjunct raceways in an industrial-scale blast furnace | |
Qu et al. | Kinetic characterization of flash reduction process of hematite ore fines under hydrogen atmosphere | |
Yi et al. | Mathematic simulation of heat transfer and operating optimization in alumina rotary kiln | |
Feng et al. | Drying and preheating processes of iron ore pellets in a traveling grate | |
Ohno et al. | Effect of coke combustion rate equation on numerical simulation of temperature distribution in iron ore sintering process | |
Dave et al. | A Mathematical Model for Straight-Grate Iron Oxide Pellet Induration Furnace: Formulation, Plant Scale Validation, Implementation and Control | |
Yang et al. | Mathematical model of thermal processes in an iron ore sintering bed | |
Mahiques et al. | Simulation of reacting, moving granular assemblies of thermally thick particles by Discrete Element Method/Computational Fluid Dynamics | |
Park et al. | Review of computational fluid dynamics modeling of iron sintering process |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |