内壁多凸起炉膛对增强辐射换热强度的Monte—Carlo方法计算.doc

内壁多凸起炉膛对增强辐射换热强度的Monte—Carlo方法计算 第32卷第6期 1997年6月 钢铁 lRONANDSTEEL V0I_32mN6 June1997 f一- 内. 卜 壁多凸起炉膛对增强辐射换热强度的 Monte—Carlo方法计算 任雁秩李义科 ——一 — —一 (包头钢铁学院) 仁.协 l\ 摘要用Monte—Carlo方法对内壁多凸起炉膛增强辐射换热强度的作用进行了定量计算,并讨 论了炉膛内壁黑度,内壁凸起物的排列方式,凸起物高度等因素对辐射换热强度增强程度的影响. ,b^一 关键{司内壁多凸起炉膛辐射换热强度Monte-Carlo椿虬——一一 ～. 一 \ cALcuLATIONoFHEATExcHANGEINTHEFURNAC]~CHAMBER WITHCONVEXINNER—WALLBYMONTE—CARLOMETHOD RENYanqiuLIYike (BaotouUniverskyofIronamp;SteelTechnology) ABSTRACTInthispaper,thefluxdensityofradiativeheatexchangeinthefurnacechamberwith convexinner-walliscalculatedbyMonte—Carlomethod.Theeffectoftheblacknessoflininginner— surGce,theheightandarrangementpatternofconvexpartsoninner—walluponheatexchangeis discussed. KEYWORDSfurnacechamberwithconvexinner-wall,fluxdensityofradiativeheatexchange, Monte—Carlomethod 1前言 内壁多凸起炉膛是在炉膛内壁表面上布 设许多凸起物以加大炉膛内表面积,即提高 炉围展开度以增强炉膛热交换,这样可以提 高炉膛的加热速度,降低能耗,提高加热质量 和加热设备的产量脚. 内壁多凸起炉膛结构增强炉膛热交换的 理论分析和实践应用均缺乏其对增强辐射换 热强度的定量结果,为此,本文用所建立的内 壁多凸起炉膛辐射换热数学模型进行了分析 计算. 2数学模型 封闭体系内任意两个区域间的净辐射换 热能流Q表示为区域黑体辐射力之差和其 间的辐射总交换面积(,S,)的乘积[: 0一(.S)(Eb—Eb)(1) 由式(1)可见,,两区域所发射的辐射 能流通过各传播途径的相互交换均在辐射总 交换面积(.S)中得以体现.内壁多凸起炉 膛对被加热物体辐射的增强程度在式(1)中 表现为使(.S)增加,故(S)的增加程度反 映了内壁多凸起炉瞳辐射换热的增强程度. 2.1辐射总交换面积的计算模型 采用Monte—Carlo方法直接计算辐射总 交换面积.计算模型为1m×1m×1m的室 状炉,其底面为被加热表面,在其炉顶内置有 联系凡任雁秋.副教授,包头(014010)包头钢铁学院环境工程系 1997年第6期 2cm×2cm×2cm正方形凸起物计算炉 顶,炉墙对炉底的辐射总交换面积.如图1所 示. 3 2l5 6 图1计算模型示意圈 Fig.1Diagramofcalculatingmodel 2,2Monte-Carlo方法计算原理与模型 2.2.1计算原理 由于辐射总交换面积和封闭体系的温度 水平无关,为简化起见,在一等温体系内进行 计算 根据辐射总交换面积的定义,某一表面 i对面的辐射总交换面积,是i面在其有效 发射率下所能发出的所有光子经各种途径 (直射,反射)到达J面且被其吸收的部分(百 分比)乘以i面的面积.因此,对i面发出的 光子要严格跟踪,直至其被任意面面所吸 收. 2,2.2光子的发射模型 在发射面Al上均匀分布Ⅳ个光子,每 个光子所携能量为 W—E?A./N(2) 设定发射表面的发射率,抽取随机发 射率E(￡一R…R——0～1均匀分布随机 数). 若E≤,光子从发射表面发出,跟踪发 射出的光子;若egt;,光子不能发出. 2.2.3光子的跟踪模 图2计算程序框图 Fig.2Blockdiagramofcalculatingprogram 橱铁 随机发射角口(经角),纬角)按下式求 取: 2.2.5辐射总交换面积的计算 发射面对接受面J的辐射总交换面积 口=arccos尺;(3)为 :2积, (4)(s .s):{(,』)+月r(,J)]/Ⅳ}A(5) 式中,尺厂0～1均匀分布随机数. 2.2.4光子吸收判断模型 设定接受表面的吸收率,随机吸收率 a=R(尼——0～1均匀分布随机数) 若≤,i表面发射光子被接受表面 吸收且计数,n(i,j)=n(f,j)+1I若gt