ISA 炉余热锅炉水平烟道流场、温度场优化模拟

陈飞; 宋波; 杨占兵; 李振祥

doi:10.13264/j.cnki.ysjskx.2016.04.001

ISA 炉余热锅炉水平烟道流场、温度场优化模拟

陈飞^{1, 2},
宋波^{1, 2, ,},
杨占兵^{1, 2},
李振祥³

1.
北京科技大学钢铁冶金新技术国家重点实验室，北京 100083
2.
北京科技大学冶金与生态工程学院，北京 100083
3.
赞比亚谦比希铜冶炼有限公司，基特韦 POB23558

基金项目:

国家自然科学基金资助项目 51234001

详细信息

通讯作者:
宋波（1963-　），男，教授，博导，主要从事夹杂物冶金、钢液净化与凝固细晶等方面的研究，E-mail: songbo@metall.ustb.edu.cn.

中图分类号: TK172;TF536
计量
- 文章访问数: 105
- HTML全文浏览量: 46
- PDF下载量: 6
出版历程
- 收稿日期: 2015-12-01
- 发布日期: 2016-08-30
- 刊出日期: 2016-06-30

Optimum simulation of velocity field and temperature field of horizontal passage of ISA waste heat boiler

CHEN Fei^{1, 2},
SONG Bo^{1, 2, ,},
YANG Zhanbing^{1, 2},
LI Zhenxiang³

1.
State Key Laboratory of Advanced Metallurgy, University of Science and Technology Beijing, Beijing 100083, China
2.
School of Metallurgical and Ecological Engineering, University of Science and Technology Beijing, Beijing 100083, China
3.
Chambishi Copper Smelt Co. Ltd., Kitwe POB23558, Zambia

摘要

摘要: 余热锅炉是工业上回收烟气余热的重要设备，其换热效率的高低与烟道内部温度场分布密切相关. 某公司的ISA炉余热锅炉，其结构呈三段式（上升烟道-下降烟道-水平烟道），烟气由下降烟道进入水平烟道后，受惯性作用，主要从对流换热区下部经过，难以充分发挥对流换热区的功能. 作者借助数值模拟的手段，在下降烟道与水平烟道连接处添加导流装置，强制引导烟气流向对流换热区.模拟结果表明，添加导流装置后，水平烟道竖直方向速度场分布趋于均匀：速度范围介于2.0~3.5 m/s的区域扩大，而大于3.5 m/s与小于2.0 m/s的区域相对减小，且速度最大值与最小值均有所收缩.单位时间内，进入对流换热区热量由5.27 MJ提升为7.70 MJ，而进入灰斗的热量由4.89 MJ减少为2.90 MJ. 同时，回流现象减弱，水平烟道上部回旋低温区范围缩小，有利于避免低温腐蚀的发生.
- 余热锅炉 /
- 数值模拟 /
- 流场 /
- 温度场 /
- 分流
Abstract: Waste heat boiler (WHB) is important equipment recovering waste heat of the flue gas in industry, the recovery efficiency of the WHB is closely related with its internal temperature field distribution. The ISA WHB has a three-passage structure（rising passage-descending passage-horizontal passage）, and exhaust gas, after entering from the descending passage, flows underneath the convection banks in the horizontal passage, under the effect of inertia, which has an adverse effect on heat exchange of the horizontal passage. With the help of the simulation method, a flow guide plate was installed around the corner to lead the gas to convection banks. The results show that with the help of the plate, the distribution of velocity in the horizontal passage tends to uniform: the range from 2.0 m/s to 3.5 m/s is expanded, and the maximum and minimum of the velocity are both contracted. The flow rate of enthalpy through convection banks rises from 5.27 MJ to 7.70 MJ, while the flow rate of enthalpy through underneath the banks declines from 4.89 MJ to 2.90 MJ. At the same time, the district of low temperature caused by backflow is narrowed, which is beneficial to low temperature corrosion.
- waste heat boiler /
- numerical simulation /
- flow field /
- temperature field /
- distributor

HTML全文

离子交换技术是一种液相组份分离技术，具有优异的分离选择性与很高的浓缩倍数，操作方便，效果突出。因此，离子交换技术作为一种优良的提取、浓缩、精制手段，在各种回收、富集与纯化等冶金、化工作业中得到了广泛的应用和迅速的发展^[1]。

目前，人们对离子交换过程的数模研究大都限于静态模型研究，而实际生产实践中广泛采用的柱操作过程动态数模及仿真研究迄今未见报道。笔者针对离子交换工艺柱操作过程的特点，基于离子交换动力学机理，运用菲克第一、第二定律，对柱内溶液流动微元层分析研究，从而建立了离子交换柱过程动态数学模型，并基于该模型以钨离子交换过程为例进行了仿真实验，通过淋洗曲线对比，结果表明该模型与实际离子交换过程吻合程度较好，用于指导离子交换生产实践、优化工艺条件具有一定的实用价值。

1 数学模型的建立

树脂在溶液中的离子交换过程可分为以下几个步骤：

(1) 溶液中的离子通过对流和扩散到达树脂表面的静止液膜；

(2) 离子通过静止液膜扩散到树脂表面；

(3) 离子在树脂内部进一步扩散；

(4) 扩散进入的离子与树脂上的相关离子发生交换；

(5) 交换下的离子在树脂内部扩散；

(6) 交换下的离子通过静止液膜扩散进入溶液；

(7) 交换下的离子在溶液中对流、扩散^[2]。离子交换过程的总速度取决于上述步骤中最慢的步骤，在一般情况下，步骤(1) (7)速度相对较快，可忽略不计。步骤(2) (6)一般称为外扩散，步骤(3) (5)为内扩散，步骤(4)为离子交换化学反应。笔者在综合考虑外扩散、内扩散、交换反应，以及柱内溶液流动情况的基础上，建立了离子交换工艺柱过程动态数学模型。

1.1 离子交换工艺柱过程外扩散数学模型

设离子交换柱内溶液流动状态为理想活塞流, 没有反混现象, 并且流速恒定。树脂床层高为H，在微元层dh中，溶液体积为v。考虑一体积为v的微元溶液层以线速度u自树脂床层上部H处向下流动，当流至h时，溶液中i离子的浓度从初始浓度Y_i变为y_i，设此时树脂层内树脂颗粒内部最外层溶液中i离子的浓度为z_i, 由菲克第一定律有：

(1)

由：，有：

代入(1)式，加上初始条件得：

(2)

(3)

式中：S—微元层中树脂颗粒总表面积，cm²；

δ—树脂颗粒表面静止液膜层厚度，cm；

u—柱内溶液流动线速度，cm/min；

v—微元层中树脂颗粒外部溶液体积，L；

D_i—溶液中i离子外扩散系数；

y_i—颗粒外部溶液中i离子浓度，mol/L;

z_i—颗粒表面层内部溶液中i离子浓度(待求)，mol/L。

1.2 离子交换工艺柱过程内扩散数学模型

设树脂颗粒为球形，半径为R，离子在树脂颗粒内部沿径向扩散，在球坐标下，暂不考虑离子交换反应，根据菲克第二定律，有：

(4)

式中：D'_i —i离子内扩散系数

g_i—颗粒内部半径为r处溶液中i离子浓度

将球形树脂颗粒沿径向划分为若干微元壳层，离散化(4)式得：

(5)

在最内层，即球心r=0处：

(6)

在最外层，即球表面r=R处：

(7)

1.3 离子交换反应数学模型

由(5)-(7)式知，扩散Δt时间后，颗粒微元层中i离子浓度的变化量为Δg_{i, r}由于在树脂颗粒中有离子交换反应存在，所以颗粒微元层中i离子浓度的实际变化量Δz_i，r(或者说i离子浓度的实际值)还与交换反应有关。

设有n种离子参于交换反应，根据树脂离子交换反应原理，其独立反应有n-1个，交换反应式通式可表示如下：

n_iL_i+m_jR_j=m_jL_j+n_iR_i反应平衡常数为Kij式中：Li，Lj—颗粒微元层溶液中的i，j离子；

Ri，Rj—交换到颗粒微元层树脂上的i，j离子；

j—在n种离子中的任选一种，设为n；

i—为1，2，3，…n-1。

由反应平衡，有n-1个方程

由物料平衡，有n个方程

式中：Li—颗粒微元层溶液中的i，j离子；

Ri—交换反应后，颗粒微元层树脂上的i离子；

Ri'—交换反应前，颗粒微元层树脂上的i离子。

由交换饱和容量，有1个方程

TR为树脂交换容量，ΔW为微元层树脂重量，由初始参数，联立解以上方程即可以求解出不同时刻，不同位置的离子浓度。

2 计算机仿真及数学模型验证

人们开发了离子交换法制取精制钨酸钠的工艺，其实质是：采用强碱性717阴离子交换树脂处理钨酸钠溶液以吸附钨，同时在净化除硅、磷、砷杂质的工艺基础上选择一定解吸剂解吸^[3]。通过以上研究我们已经建立了离子交换的数学模型，在此基础之上，以钨离子交换为例进行计算机仿真研究。该仿真软件采用Delphi编程工具，应用图型化界面，可以任意设置离子交换过程的不同参数，如交换液流速, 穿漏浓度，淋洗液浓度和淋洗速度等。该软件还具有动画显示功能，能形象地显示出交换柱内不同时刻不同离子的浓度分布情况。设置参数后，通过计算模拟，能够得出不同参数下的流出曲线，交换所需时间和不同时刻的离子浓度分布。

2.1 流程与算法

采用龙格-库塔(Runge-Kutta)解扩散微分方程，迭代法解交换反应方程，计算流程框图见图 1。

图 1 离子交换过程仿真计算框图

下载: 全尺寸图片幻灯片

2.2 验证

以下是在淋洗参数为:WO3料液浓度20g/L，树脂交换容量为3mmol/g，吸附期线速度5cm/min，不同淋洗浓度和淋洗速度，实验与仿真淋洗曲线对照图，见图 2~图 5。

图 2 不同浓度NH₄Cl实验淋洗曲线(3cm/min)

下载: 全尺寸图片幻灯片

图 3 不同淋洗速度实验淋洗曲线(4.5mol/LNH₄Cl)

下载: 全尺寸图片幻灯片

图 4 不同浓度NH₄Cl仿真淋洗曲线(3cm/min)

下载: 全尺寸图片幻灯片

图 5 不同淋洗速度仿真淋洗曲线(4.5mol/L NH₄Cl)

下载: 全尺寸图片幻灯片

从图 2和图 4可以看出在同一淋洗速度之下，随着淋洗浓度的增加, 淋洗峰值越大，淋洗区间越小，且不同浓度对比，淋洗峰值大致一致。从图 3和图 5可以看出在相同淋洗浓度下，淋洗速度越小淋洗峰值越大, 淋洗区间越小，不同淋洗速度峰值比较大致一致。从以上图可看出，仿真结果与实验结果基本相符。

3 结语

离散化离子交换柱过程，在已有的扩散动力学的基础之上，建立离子交换过程的数学模型，笔者在建立的离子交换数学模型之上，以钨离子交换为例, 运用Delphi语言编制的钨离子交换计算机仿真软件能灵活，方便，快速地再现出不同工艺参数对钨离子交换流出曲线的影响。通过钨离子交换拄过程的仿真研究，对深刻理解不同离子交换参数和流出曲线的关系有重要的意义，能够全面形象地反应钨离子交换过程不同时刻和不同树脂层位置各离子的浓度变化。

图 1 余热锅炉正视图

Fig 1. Front view of WHB

下载: 全尺寸图片幻灯片

图 2 对流管束结构示意

Fig 2. Sketch of convection bank structure

下载: 全尺寸图片幻灯片

图 3 余热锅炉几何模型

Fig 3. Geometric model of WHB

下载: 全尺寸图片幻灯片

图 4 截面z=0m流场云图

Fig 4. Filled contours of velocity at cross section z=0m

下载: 全尺寸图片幻灯片

图 5 速度矢量

Fig 5. Vectors of velocity

下载: 全尺寸图片幻灯片

图 6 截面z=0m温度云图

Fig 6. Filled contours of temperature at cross section z=0m

下载: 全尺寸图片幻灯片

图 7 水平烟道右视图

Fig 7. Right side view of horizontal passage

下载: 全尺寸图片幻灯片

图 8 截面x= 0 m,x= 4.476 m 上的速度矢量

Fig 8. Vectors of velocity at cross section x= 0 m, x= 4.476 m

下载: 全尺寸图片幻灯片

图 9 下降段与水平段交界面位置速度云图

Fig 9. Filled countours of velocity at interface between descenting and horizontal passages

下载: 全尺寸图片幻灯片

图 10 导流板形状/位置示意

Fig 10. Sketch of shape/positon of flow guide plate

下载: 全尺寸图片幻灯片

图 11 加装导流板前后速度矢量对比

Fig 11. Vectors of velocity with/without flow guide plate

下载: 全尺寸图片幻灯片

图 12 加装导流板前后袁截面x= 0 m 上的温度云图对比

Fig 12. Filled contours of temperature at cross section x= 0 m with/without flow guide plate

下载: 全尺寸图片幻灯片

图 13 加装导流板前后袁截面x= 4.476 m 上的温度云图对比

Fig 13. Filled contours of temperature at cross section x= 4.476 m with/without flow guide plate

下载: 全尺寸图片幻灯片

表 1 模型流体组成

Table 1 Component of flue gas

组分	SO₂	SO₃	CO₂	H₂O	O₂	N₂
体积分数/%	9.66	0.10	12.35	25.48	5.30	47.11

下载: 导出CSV

表 2 对流换热区多孔性参数

Table 2 Parameter of porous zone

对流管束	x 方向		y 方向		z 方向
对流管束	C₂	1/α	C₂	1/α	C₂	1/α
1	0.73167	21 652.32	0.731 67	21 652.32	73.167	2 165 232
2	9.16124	4 064.554	9.161 24	4 064.554	916.124	406 455.4
3	13.76	7 470.182	13.76	7 470.182	1 376	747 018.2
4	13.76	7 470.182	13.76	7 470.182	1 376	747 018.2
5	13.76	7 470.182	13.76	7 470.182	1 376	747 018.2

下载: 导出CSV

表 3 加装导流板前后不同截面上通过的烟气流率

Table 3 Flow rate at different cross sections with/without flow guide plate

质量流率 /(kg.s^-1)	截面x=0 m		截面x=4.476 m
质量流率 /(kg.s^-1)	灰斗上部	灰斗区域	灰斗上部	灰斗区域
优化前	12.623	2.813	7.709	7.766
优化后	12.531	2.903	10.580	4.904

下载: 导出CSV

参考文献(16)

[1]	史跃岗, 陆成. 内补燃余热锅炉的燃气-蒸汽联合循环效率计算[J]. 余热锅炉, 2013 (4): 11-14. http://www.baidu.com/link?url=vBih0ZYYAqUjxID9p357xVshLJGjhVHazeO83akoH3jGPJnF0Ct1dvHabbubCxPdNe6hyL0fZyKBreyVFRQmyVmUqu3GkKoMuUYS1T8BHNFPYR4_-QXCs_Rl7F76ApZuwbzW2wATT_UdRha0cvAKxbXca0bPLCowGfud3J3rup9flphH4uCpYeaiZC2GPO8-Y7zjVlp6Peh8WGxpNDbQLOZE2c-x4BFqGOnd1NBRNOMq8UJzA1M1W6RsdRboH0w4nWz5csvHeXBkR5-zHLPqd5ib641mq5dvkm0WJe2IF2XtQaS8h2Xm74QSCxgOHQVmf90HtpgJfQy0c0nJHoFgO2pPqX-bcsECJyScwE6da30hsooYwcVkH9AyM5f0aj1Xw8oEPh8v3yyH-wo7oXzW0mJmKNF-S7SN0MxcW5lGq8x9y_puN5o6pXIA-DuvdSU1AYG4IuK1T34qRfmLG9oLdyASdRiZi6cnEtNhDUHuLX_&wd=&eqid=90fe42ff0006ff3c00000005580ee6b5
[2]	王奇. 艾萨-余热锅炉的设计特点分析[J]. 余热锅炉, 2008 (3): 1-4. http://www.baidu.com/link?url=f8E3JkI4_a0EcDZhYwfAB4cqNNMps-ONHDZpNLYxMvsj7lIRuFakq8811me4LQpOnuUXFCqbPA6EqnaRNaFQyKRAVPyktmvDkSagDbl14x0TKf4Q1dOY_SS0_6e4PU5V5gXGO9vId0fTDUHUgGUytOSVEI_sJIbsevpbi086Re0QoZ93V8RBLjZ8bqYmMx-Ix8ZL0Z37u-jgSrAfLB_uTujXEAxDHGoxkDTKjtSe-ucCNJRUDDhY0DNBlvRBMKNwsh0pT2vFKHc8UV6sPjwnOHAL0WXnlb1vfGSKRN06i1JWprDCt-5-3rE5msdfSs1KCxe-u1RwXM6fInbaDyVY2E0CVRN4UEOQm6vNAPwA4FS7bcnH6Xaq4adwQ9KynT1BQ_k2tyhJfi55kzFlOJo6P_&wd=&eqid=9d03cd1900063e6a00000005580ee6c4
[3]	陈宇昂. 余热锅炉设计的若干研究[J]. 余热锅炉, 2015 (3): 30-33. http://www.baidu.com/link?url=eQNLJO6uvC8TYznunjSIXltwMdaV2hUzHRM8SXGH0MAGmvwXJ08faaSrZ2spkrUGiZ8BzB2Dj7zydol9r1IHTDfli9S473HXHipO9vrzUEW6Y76ogNlk74c_64Dxk-14IHMhW7KmJTf-rNbeyfyVulOVMchmSQr_194pS--pJj6TJltc68XpMj8DRAfL6KqUkCl9BL19SmESD-6utT0n4q91h40byfKNPHDv-tfCdfwS8jfj_fgeHHaGmn5lnbOamQKa0Te_mCb_Xcom3cqcFvOcxVztNiuj_iEFHHYZZemV_Xe_B0-X_GMapn6jFMv-VUGiv600mqC4aWytsH4Uc7p1jmLp88JOfB6RQdMkAl_kEtXUJ5gTU7gRqDx76BCp&wd=&eqid=94c34d8f00066a2f00000005580ee6d6
[4]	宋冬根，杨平，鞠霞. 闪速炼铜余热锅炉炉型改造与优化[J]. 有色冶金设计与研究，2013, 34(1): 17-20. http://www.cnki.com.cn/Article/CJFDTOTAL-YSYJ201301009.htm
[5]	赵钦新. 余热锅炉研究与设计[M]. 北京：中国标准出版社, 2010.
[6]	王晓瑜. 余热锅炉结构形式及运行特性研究综述[J]. 中国有色冶金，2012(1): 40-45. http://www.cnki.com.cn/Article/CJFDTOTAL-YSYL201201019.htm
[7]	钟崴, 丁晟, 宋冬根, 等. 闪速炼铜余热锅炉辐射室流场温度场数值模拟[J]. 浙江大学学报(工学版), 2012, 46(2): 321-326. http://www.baidu.com/link?url=GDnKFppeX_StcK0mNc5JgF13Qc27D9fNNoMH0thNYNbB3APgCRSEXQwQxjuvcUkIKi7-Ivr7QbRP1ohxjl5OHMYD71J6mlVZVrUz08F4fVHvhDjGgnaaf6r3jGCn7DvRURc_OErGncPccTL9E7cc5n0mud10MrIe_VrZzA9cF0Qg-vqN1_N3ItUuXNvhe0-bo_sgPBDipeYMJ4QlUh2y_ISGth1rSuBl11C1DtAFR47xHlhd3kjFST6ttB-doBi0rML-Z5gUvhYp2q8Apt8MmFZAl659gzsa7EXKGzS5El5OYQoz1C2TFa4rXr-6UCemfk_4Kkp97WV1Sw8E4XUifsPhZLb3iCRY9Bq8L1bA2tvFGaJ5Pj89Y4MxXnbAJz9txO8z9PfsZgGhODim81Ea_EaaievCcEXwuzWWhS3FJ6xqsgSvq9198DSNQqFS0FgZtS7U_alww84Y7gRWQqU9FdQtPYeFobIKtYbwHotd0se&wd=&eqid=9c33c5810006e6a100000005580ee6f1
[8]	王晓瑜, 袁艳平, 张冬洁, 等. 顶吹熔炼炉余热锅炉数值模拟[J]. 中国有色冶金, 2013, 42(3): 50-53. http://www.cnki.com.cn/Article/CJFDTOTAL-YSYL201303020.htm
[9]	LAUNDER B E, SPALDING D B. Lectures in mathematical models of turbulence[J]. London: Academic Press, 1972. http://cn.bing.com/academic/profile?id=180824322&encoded=0&v=paper_preview&mkt=zh-cn
[10]	LAUNDER B E, SPALDING D B. The numerical computation of turbulent flows[J]. Computer Methods in Applied Mechanics and Engineering, 1974, 3(2): 269-289. doi: 10.1016/0045-7825(74)90029-2
[11]	HENKES R, VAN DER VLUGT F F, HOOGENDOORN C J. Natural-convection flow in a square cavity calculated with low-Reynoldsnumber turbulence models[J]. International Journal of Heat and Mass Transfer, 1991, 34(2): 377-388. doi: 10.1016/0017-9310(91)90258-G
[12]	王学猛，张新庄. 炼铅炉渣SKS炉底烟化提锌的试验探索[J]. 中国有色冶金， 2015(3): 12-13. http://www.cnki.com.cn/Article/CJFDTOTAL-YSYL201503005.htm
[13]	CHUI E H, RAITHBY G D. Computation of radiant heat transfer on a nonorthogonal mesh using the finite-volume method[J]. Numerical Heat Transfer, 1993, 23(3): 269-288. doi: 10.1080/10407799308914901
[14]	RAITHBY G D, CHUI E H. A finite-volume method for predicting a radiant heat transfer in enclosures with participating media[J]. Journal of Heat Transfer, 1990, 112(2): 415-423. doi: 10.1115/1.2910394
[15]	赵钦新. 余热锅炉研究与设计[M]. 北京：中国标准出版社, 2010.
[16]	聂宇宏, 聂德云, 梁融, 等. 余热锅炉灰斗内两相流动的三维数值模拟及结构优化[J]. 江苏科技大学学报(自然科学版), 2013, 27(5): 448-492. http://www.cnki.com.cn/Article/CJFDTOTAL-HDCB201305017.htm

施引文献

资源附件(0)

图(13) / 表(3)

计量

文章访问数: 105
HTML全文浏览量: 46
PDF下载量: 6
被引次数: 0

1 数学模型的建立
1.1 离子交换工艺柱过程外扩散数学模型
1.2 离子交换工艺柱过程内扩散数学模型
1.3 离子交换反应数学模型
2 计算机仿真及数学模型验证
2.1 流程与算法
2.2 验证
3 结语

ISA 炉余热锅炉水平烟道流场、温度场优化模拟

通讯作者: 宋波（1963- ），男，教授，博导，主要从事夹杂物冶金、钢液净化与凝固细晶等方面的研究，E-mail: songbo@metall.ustb.edu.cn.

计量

出版历程