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 ℃等等
余热锅炉是生产蒸汽的设备。 但是由于余热资源条件的多样性与复杂性, 使余热锅炉比工业锅炉具有不同的特点。 人们很早就研究余热锅炉, 一般主要考虑针对不同的余热资源设计余热锅炉, 从余热锅炉的积灰、 腐蚀、 磨损机理研究余热锅炉, 并取得了很好的效果
本文作者采用计算机仿真模拟的方法, 用k-ε方程模拟炉内气体流动, 并对某铜业公司的铜闪速炉烟气余热锅炉进行了工业规模的数值仿真计算, 以研究余热锅炉内气体流动与传热特性, 得出了余热锅炉内流场和温度场的分布。
1 基本方程和边界条件
1.1 基本方程
求解流场, 采用k-ε模型
连续和动量方程:
式中 ρ为流体平均密度, U= (U, V, W) 为流体平均速度, p为压力, t为时间, B为体积力, μeff为有效粘度。
紊流动能k和紊流扩散率ε的传输方程:
式中 P为剪切应力产生项; G为体积力产生项; C1, C2和C3均为经验常数
能量方程:
式中 H是总热焓, 由静态焓h的公式给出:
H=h+
式中 ρ是流体的平均密度, U= (U, V, W) 是流体的速度, p是压力, T为温度, t为时间, λ是导热系数。
传递方程必须加上本构方程才能封闭, 也就是说, 要加上静态焓与温度和压力之间的关系, 即h=h (T, p) 。 假设流体为温度特征值 (thermally perfect) , 即静态焓只是温度的函数, 则本构方程可由比定压热容和基准温度Tref给出。 比定压热容是温度的函数:
Tref是静态焓为0时的温度。 那么静态焓可由上式积分得到:
h=∫
在程序内部, 本构方程被假设成具有下列解析形式:
式中 Tref取环境温度。
余热锅炉内混合气体的物性参数仅考虑与温度有关, 可写成温度的多项式:
ρ=0.005 388 p/T
cp (T) =793.13+0.208 T
ν (T) =2.0×10-6+3.0×10-9+6.0×10-11T2
λ (T) =2.277 45×10-3+5.157 1×10-5T (10)
式中 ρ为混合气体平均密度, cp为混合气体比定压热容, ν为动力学粘度, λ为导热系数。
传热计算的目的是获得计算区域内的温度场。 为了耦合温度场内的辐射, 将辐射项处理为方程的源项。 在控制体积i, j, k中的辐射散热为
Cijk=4σT
式中 σ是斯蒂芬常数, T为温度, kijk是吸收率, Vijk是单元体积。
在每个辐射区域辐射属性假设为常数, 因此对一个辐射区域, 我们有总散热为
因此
式中的求和是针对每个属于辐射区域的体积i, j, k而进行的。
为使能量守衡区域Z的辐射温度由总辐射散热计算, 则
CZ=4σkZT
这样,
使用这个温度场, 计算辐射加热量HZ。 在辐射体积中净总辐射加热量—焓源项, 由下式给出:
H
式中 aZ是保证能量守衡的体积标准化因子
焓源项Hr由因子ur欠松弛, 这样在对计算的第n次调入中, 有
1.2 边界条件
进口条件: 进口采用Dirichlet条件, 直接设定进口速度 (或流量) 。 进口湍动能k和动能耗散系数ε的值对于进口均匀截面可以作如下假设
kint=cp1u
式中 uint为进口平均速度; cp1和cp2为经验数据, 分别为0.05和0.3
出口条件: 出口采用压力边界。
壁面条件: 紊流传输方程只能在全紊流区适用, 对于近壁面, 我们选择对数律来处理底层区域。 对于水冷壁, 管内介质为水, 采用定温条件, 根据不同情况, 取为250~550 ℃, 辐射屏温度也采用定温条件, 取为700 ℃
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%, H2O含量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) 计算机模拟为余热锅炉的设计、 节能改造提供了详细的数据和有效的手段。 它能替代不少实验, 并获得许多实验难以提供的数据, 可降低实验费用, 缩短设计周期。
参考文献
[7] SchwarzMP .Lowsimulationinmineralengineering[J] , MineralsEngineering , 1 991 , 7 (4) :71 7.