模拟辐射传热离散坐标法改进.pdfVIP

维普资讯 计 算 物 理 第 22卷 第 3期 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Vo1．22．No．3 2005年 5月 May．2005 [文章编号] 1001—246X(2005)03—0277—06 模拟辐射传热的离散坐标法的改进 杨 占春 ， 武文斐 ， 李义科 (内蒙古科技大学 能源与环境学院，内蒙古 包头 014010) [摘 要] 针对离散坐标法模拟求解炉内辐射换热的问题 ，在对辐射传递方程采用控制容积法离散时，按辐射 方向上各界面之间的投影关系划分区域 ，对辐射传递方程在各子控制体上沿辐射传输方 向进行积分 ，求解各子 控制体的辐射强度 ，再取容积平均值作为最终节点的辐射强度的新离散坐标方法 ，定义 为分段积分离散坐标法 ， 并对三维矩型燃气炉 内辐射换热进行了数值模拟 ，与传统离散坐标方法及 区域法进行了比较 ，比较结果表 明新 离散坐标法较好的解决 了假散射 问题 ，提高了离散坐标法的精度 ． [关键词] 辐射传热 ；离散坐标法 ；假散射 ；分段积分 [中图分类号] TK124 [文献标识码] A O 引言 随着计算机技术和计算数学的迅速发展，辐射换热数值计算作为传热学 的重要组成部分，近 2O多年得 到了空前发展 ．离散坐标法的主要优点是将辐射传递方程转化成微分表达式，便于同一般输运方程耦合求 解 ，计算工作量小，而且近几年出现了几种新的离散坐标取法，如 SRAP_l法和 法 ，使得离散坐标法的精 度大幅提高 ，可 以和区域法媲美 ，因此离散坐标法将是今后一段时期 内发展和应用较多的一种辐射换热模 型 ．然而，其本身仍存在难 以克服的缺点，如射线影响、假散射、对差分格式敏感等，笔者认为假散射的根源在 于对离散方程进行积分时 ，是按坐标方 向进行的，没有按传播方 向进行，从表 1中可 以看 出，在很多情况下， 面根本影响不到E面，而是对 面和Ⅳ面产生影响，传统做法忽略了这一点， 面根本和 E面没有关系， 却用 面来推算E面，相反，Js面和B面必然对E面产生影响，却没有考虑进来，为此 ，本文提出一种按各面 之间在传播方向上的投影关系，将微元控制体划分区域 ，分段对辐射传递方程进行积分的方法，定义为分段 积分离散坐标法 ． 1 辐射传递方程的坐标离散 离散坐标法是基于对辐射强度的方向变化进行离散，通过求解覆盖整个 4兀立体角空间的一套离散方 向 上的辐射传递方程而得到问题的解 ，在三维直角坐标系中，对于吸收、辐射 、各 向同性散射介质的辐射传递方 程可写为 + +7=一(+ ，+ko，+』他(n，n)dn， (1) 漫灰表面的边界条件为 ，=￡， + I ，In·sIdn． (2) 在离散坐标法 中，积分项用数值积分代替，写为 + m +等++ 一：一((a+))lm++¨』6++∑COm,………lm ((33) 式中的方向余弦 ， ， 和权值 根据离散坐标取法确定，这里不做阐述 ，参见相关文献 ．对于漫 一灰体 边界壁面 ，相应的边界条件为 [收稿 日期]2003—11—10；[修回日期]2004—04—08 【作者简介]杨 占春 (1977一)，男 ，河北滦县 ，硕士生，主要从事热能工程方面的研究，北京钢铁研究总院工艺所 100081 维普资讯 278