DOI:10.19476/j.ysxb.1004.0609.2000.02.029
余热锅炉流场温度场数值计算
杨卫宏 赖亚欣 萧泽强
中南工业大学物热系!长沙410083
摘 要：
用k ε方程模拟了余热锅炉内流场和温度场分布。几何拓扑结构用多块结构以减少计算机内存的浪费 , 用贴体坐标系统来描述复杂的几何边界 , 网格系统用交错网格。使用二阶中心差分法离散扩散项和源项 , 对于控制方程的对流项的离散采用混合差分 (HDS) 格式。求解方程使用SIMPLE算法。给出了余热锅炉内流场和温度场的分布 , 为研究余热锅炉炉内流动和传热、传质提供了理论依据 , 对余热锅炉的设计、节能改造提供了详细的数据和有效手段
关键词：
余热锅炉 ;流场 ;温度场 ;数值计算 ;
中图分类号： TK224
收稿日期： 1999-02-10
基金： 教育部博士点基金资助项目;
Numerical simulation of flow and temperature fields in waste boiler
Abstract：
The flow and temperature fields in the waste boiler have been simulated. The computing method is SIMPLE method. The multi block grid methodology involving unstructured set of block was used, and the Body Fitted grid was specified in the simulation. The method of discretization of diffusion terms and source terms in the control equations is central differencing (CDS) , and the method of discretization of advection terms in the control equations is Hybrid Differencing (HDS) . The computation results described by the distribution of the velocity and temperature are good in accordance with phenomena observed in prototypes. The simulation results can be used in design of waste boiler.
Keyword：
waste boiler; numerical calculation; flow field; temperature field;
Received： 1999-02-10
在冶金工艺过程中, 有着极其丰富的余热资源, 如铜、 镍、 锌闪速炉的排烟温度一般为1 200～1 400 ℃, 炼铜、 炼镍反射炉排烟温度为1 100～1 370 ℃等等
[1 ]
。 我国冶金行业每年消耗原煤一亿多吨, 其中可利用的余热资源就相当于3 300×104 t原煤
[1 ]
。 这些余热排放, 一方面造成能源浪费, 同时由于高温烟气中含尘量高、 腐蚀性大、 具有粘着性等也产生环境污染。 目前, 国内外为回收烟气中的大量余热和降尘一般采用余热锅炉。 如芬兰奥托坎普铜闪速炉的烟气就采用余热锅炉回收, 我国贵溪冶炼厂铜闪速炉的烟气也采用余热锅炉回收, 都收到了较好的效果。
余热锅炉是生产蒸汽的设备。 但是由于余热资源条件的多样性与复杂性, 使余热锅炉比工业锅炉具有不同的特点。 人们很早就研究余热锅炉, 一般主要考虑针对不同的余热资源设计余热锅炉, 从余热锅炉的积灰、 腐蚀、 磨损机理研究余热锅炉, 并取得了很好的效果
[2 ]
。 最近几年, 开始有学者进行余热锅炉内流动和传热、 传质现象的研究, 包括余热锅炉内管壁上灰渣沉积对传热的影响等
[3 ,4 ]
; 芬兰奥托坎普铜业公司则建立了余热锅炉模型, 并进行了实验研究。 国内这方面的研究, 特别是对采用火法提取有色金属的尾部烟气的余热锅炉的内流动传热研究尚甚少见。
本文作者采用计算机仿真模拟的方法, 用k -ε 方程模拟炉内气体流动, 并对某铜业公司的铜闪速炉烟气余热锅炉进行了工业规模的数值仿真计算, 以研究余热锅炉内气体流动与传热特性, 得出了余热锅炉内流场和温度场的分布。
1 基本方程和边界条件
1.1 基本方程
求解流场, 采用k -ε 模型
[5 ]
, 控制方程如下。
连续和动量方程:
?
ρ
?
t
+
?
?
(
ρ
U
)
=
0
?
?
?
(
1
)
?
ρ
U
?
t
+
?
(
ρ
U
?
U
)
-
?
?
(
μ
e
f
f
?
U
)
?
?
=
-
?
p
+
?
?
(
μ
e
f
f
(
?
U
)
Τ
)
+
B
?
?
?
(
2
)
式中 ρ 为流体平均密度, U U , V , W ) 为流体平均速度, p 为压力, t 为时间, B μ eff 为有效粘度。
紊流动能k 和紊流扩散率ε 的传输方程:
?
ρ
k
?
t
+
?
?
(
ρ
U
k
)
-
?
?
(
(
μ
+
μ
Τ
σ
k
)
?
k
)
=
?
?
Ρ
+
G
-
ρ
ε
?
?
?
(
3
)
?
ρ
ε
?
t
+
?
?
(
ρ
U
ε
)
-
?
?
(
(
μ
+
μ
Τ
σ
ε
)
?
ε
)
=
?
?
C
1
ε
k
(
Ρ
+
C
3
max
(
G
,
0
)
)
-
C
2
?
ρ
ε
2
k
?
?
?
(
4
)
式中 P 为剪切应力产生项; G 为体积力产生项; C 1 , C 2 和C 3 均为经验常数
[5 ]
。
能量方程:
?
ρ
Η
?
t
+
?
?
(
ρ
U
Η
)
-
?
?
(
λ
?
Τ
)
=
?
p
?
t
?
?
?
(
5
)
式中 H 是总热焓, 由静态焓h 的公式给出:
H =h +
1
2
U
2
?
?
?
(
6
)
式中 ρ 是流体的平均密度, U U , V , W ) 是流体的速度, p 是压力, T 为温度, t 为时间, λ 是导热系数。
传递方程必须加上本构方程才能封闭, 也就是说, 要加上静态焓与温度和压力之间的关系, 即h =h (T , p ) 。 假设流体为温度特征值 (thermally perfect) , 即静态焓只是温度的函数, 则本构方程可由比定压热容和基准温度T ref 给出。 比定压热容是温度的函数:
c
p
(
Τ
)
=
(
?
h
?
Τ
)
p
?
?
?
(
7
)
T ref 是静态焓为0时的温度。 那么静态焓可由上式积分得到:
h =∫
0
Τ
c p T ′) dT ′-∫
0
Τ
r
e
f
c p T ′) dT ′ (8)
在程序内部, 本构方程被假设成具有下列解析形式:
h
=
c
?
p
(
Τ
)
Τ
-
c
?
p
(
Τ
r
e
f
)
Τ
r
e
f
?
?
?
(
9
)
式中 T ref 取环境温度。
余热锅炉内混合气体的物性参数仅考虑与温度有关, 可写成温度的多项式:
ρ =0.005 388 p /T
c p T ) =793.13+0.208 T
ν (T ) =2.0×10-6 +3.0×10-9 +6.0×10-11 T 2
λ (T ) =2.277 45×10-3 +5.157 1×10-5 T (10)
式中 ρ 为混合气体平均密度, c p ν 为动力学粘度, λ 为导热系数。
传热计算的目的是获得计算区域内的温度场。 为了耦合温度场内的辐射, 将辐射项处理为方程的源项。 在控制体积i , j , k 中的辐射散热为
C ijk σT
i
j
k
4
k ijk V ijk
式中 σ 是斯蒂芬常数, T 为温度, k ijk V ijk
在每个辐射区域辐射属性假设为常数, 因此对一个辐射区域, 我们有总散热为
C
Ζ
=
∑
i
j
k
∈
Ζ
4
σ
Τ
i
j
k
4
k
i
j
k
V
i
j
k
?
?
?
(
1
2
)
因此
C
Ζ
=
4
σ
k
Ζ
∑
i
j
k
∈
Ζ
Τ
i
j
k
4
V
i
j
k
?
?
?
(
1
3
)
式中的求和是针对每个属于辐射区域的体积i , j , k 而进行的。
为使能量守衡区域Z 的辐射温度由总辐射散热计算, 则
C Z σk Z T
Ζ
4
V Z
这样,
Τ
Ζ
=
(
∑
i
j
k
∈
Ζ
Τ
i
j
k
4
V
i
j
k
V
Ζ
)
1
4
?
?
?
(
1
5
)
使用这个温度场, 计算辐射加热量H Z
H
i
j
k
r
V ijk a Z H Z V ijk C ijk
式中 a Z
a
Ζ
=
V
Ζ
∑
i
j
k
∈
Ζ
V
i
j
k
?
?
?
(
1
7
)
焓源项H r u r n 次调入中, 有
Η
i
j
k
r
(
n
)
=
u
r
(
a
Ζ
Η
Ζ
-
C
i
j
k
V
i
j
k
)
+
(
1
-
u
r
)
Η
i
j
k
r
(
n
-
1
)
?
?
?
(
1
8
)
1.2 边界条件
进口条件: 进口采用Dirichlet条件, 直接设定进口速度 (或流量) 。 进口湍动能k 和动能耗散系数ε 的值对于进口均匀截面可以作如下假设
[4 ]
:
k int =c p 1u
int
2
ε
int
=
k
int
1
.
5
c
p
2
d
i
n
式中 u int 为进口平均速度; c p 1c p 2
[4 ]
; d in 为水力直径, d in =4A /L , A 为进口截面面积, L 为进口截面周长。
出口条件: 出口采用压力边界。
壁面条件: 紊流传输方程只能在全紊流区适用, 对于近壁面, 我们选择对数律来处理底层区域。 对于水冷壁, 管内介质为水, 采用定温条件, 根据不同情况, 取为250～550 ℃, 辐射屏温度也采用定温条件, 取为700 ℃
[8 ,9 ]
。
2 数值模拟及分析
2.1 数值模拟
网格划分: 采用贴体坐标系统 (BFC) , 它不但能准确描述复杂的几何边界, 还可以利用很密的网格求解而不占用计算机过多的内存; 流道几何拓扑结构采用多块结构化, 以减少内存的使用, 网格数为10×40×40, 见图1。
图1 余热锅炉计算网格
Fig.1 Computation grid for waste-heat boiler
数值求解: 当余热锅炉运行稳定后, 各参数基本不变, 故可以假设流动传热为稳态, 这样也可减少计算工作量。 同时, 针对流动控制方程的特点, 使用二阶中心差分法离散扩散项和源项; 对于控制方程的对流项的离散, 为了克服中心差分当网格Pe数大于2后解的不稳定和迎风差分不考虑Pe影响的缺陷, 采用混合差分 (HDS) 格式, 即当网格|Pe|>2时采用迎风差分, 网格|Pe|<2时采用中心差分, 对时间采用二阶向后差分; 采用SIMPLE算法解离散方程, 考虑收敛, 对速度三个分量U , V , W 和k , ε 采用欠松弛因子方法处理和假时间步 (False Time Step Relaxation) 的方法, 取得了较好的收敛结果。
2.2 计算结果及分析
所模拟的余热锅炉具体结构如图2所示, 基本尺寸为23 m×5.4 m×11 m, 辐射炉膛内有8块辐射屏, 尾部烟道内有5组对流管束。 考虑计算机内存和速度的限制, 将辐射屏和对流管束简化为一组平盘。 烟气体积流量为19 000 m3 /h, 温度为1 300～1 350 ℃; 烟气成分为SO2 30%, CO2 含量10%, H2 O含量3%, O2 含量2%, N2 含量55%。 通过计算机模拟, 得出了炉内流场与温度场的分布图。
图2 余热锅炉基本结构
Fig.2 Schematic drawing of waste-heat boiler
图3为余热锅炉中心纵截面速度分布。 在辐射炉膛管束区域内, 流动速度较大且方向较水平, 速度比较均匀。 这是因为辐射屏部分流通面积减少造成的。 辐射屏区域对主流的影响不大, 但在辐射屏左下部, 能看到一个旋流, 而在下部和右下部也有一个回流。
图3 余热锅炉中心纵截面速度分布图
Fig.3 Velocity distribution of waste-heat boiler near symmetric center line
图4为余热锅炉横截面速度分布, 在辐射屏下部出现了两个方向相反位置对称的旋涡 (对涡) , 对流管束区域内流动的影响和辐射屏类似。
图5为炉内温度分布。 由于气流有强烈的向上流动 (因为进口气流向上) , 靠近辐射屏进口部位的顶墙的热流比较大, 这部分炉膛降温较快。 辐射区域下部的温度分布较均匀, 这是由于气体在此存在回流, 气体停留时间较长, 气体冷却比较充分造成的。
图4 余热锅炉A-A截面速度分布
Fig.4 Velocity distribution across radiationscreens of waste-heat boiler
图5 余热锅炉中心纵截面温度等值线图
Fig.5 Temperature contours near symmetric center-line of boiler
3 结论
1) 炉内气流在进口顶部热流较大, 炉内下部存在旋涡, 使整个流场分布欠佳, 并影响了温度分布。 可以考虑通过改变流场而使温度分布更合理 (具体的改进措施见作者另一篇文章) 。
2) 计算机模拟为余热锅炉的设计、 节能改造提供了详细的数据和有效的手段。 它能替代不少实验, 并获得许多实验难以提供的数据, 可降低实验费用, 缩短设计周期。
参考文献
[] [1 ] LIURong (刘 荣 ) , LIZhi pu (李支普 ) .PowerEquip mentofMetallurgy (冶金动力设备 ) [M ] .Shenyang :NorthEastUniversityPress , 1 993 .
[] [2 ] BackmanR , HupaHandMakinenJK .Formationandcorrosioneffectsofsulphurtrioxideincopper, smeltingprocesses [R] .MetallurgicalSocietyofAIME , TMSTechnicalPaper, PaperNo .A56 86 , 1 986 .
[3] LiliaandRajainmaki.Ontheexperimentalandcomputa tionalmodelingatOutokumpureasearch [A] .The 2ndcolloquiumon processsimulation [C] .Jokilakso , 6 8June 1 995 , Finland , 2 1 7～ 2 30 .
[4] KlausAHandSteveTC .ComputationalFluidDynam icsforEngineers [M] .VolumeⅠ /Ⅱ , APublicationofEngineeringEducationSystem , Kansas, USA , 1 993 .
[5] PatankarSV .Transl.ZHANGZheng (张 政 ) .Nu mericalHeatTransferandFluidFlow (传热与流体流动的数值计算 ) [M] .Beijing :ScienceTechnololgyPress , 1 984.
[6] MEIChi (梅 炽 ) andWANGQian pu (王前普 ) .有色冶金炉窑的仿真与优化 [J]。TheChineseJournalofNonferrousMetals (中国有色金属学报 ) , 1 996 , 6 (4) :1 9.
[7] SchwarzMP .Lowsimulationinmineralengineering[J] , MineralsEngineering , 1 991 , 7 (4) :71 7.
[8] YANGWei hong (杨卫宏 ) andXIAOZe qian (肖泽强 ) .双级旋风分离器的计算机预报 [J] .JournalofCentralSouthUniversityofTechnology (中南工业大学学报 ) , 1 999, 30 (5) :493 .
[9] YANGWei hong (杨卫宏 ) andXIAOZe qian (肖泽强 ) .余热锅炉辐射冷却室三维数值模拟 [J] .BoilerTech nology (锅炉技术 ) , 1 999, 30 (1 1 ) :8
