贴体坐标系下复杂几何形体内辐射传递研究的论文.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：023147 贴体坐标系下复杂几何形体内 辐射传递研究 郝金波 阮立明 谈和平 (哈尔滨工业大学能源科学与工程学院，150001) (0451—6412108，haojinbo@hotmai1．com) 摘要：本文研制了贴体坐标系下有限体积法计算程序，采用该程序研究了不规则形体、各向异性散射 介质内的辐射传递。经与国外文献计算结果比较表明，本文的程序计算结果准确。在此基础上，以椭 圆柱腔体内辐射传递为例，研究了衰减系数对温度场分布影响；计算精度的选取对计算结果、迭代次 数的影响，研究表明，合理的计算精度取值，可节省计算时间． 关麓词：辐射传递，贴体坐标系，有限体积法 1 引言 高温设备的几何形状一般都很复杂；且辐射换热通常与其它传热方式耦合出现，或 与其它物理现象(流体流动、燃烧)紧密相关。适体坐标系下辐射传递方程的生成与离 散，国外近几年来才有文献发表【l-3】。显然，辐射换热的数值模拟己成为制约高温复杂 几何形体内多场耦合准确求解的关键问题之一。Liu掣‘谰有限体积法和离散坐标法分 析了三维不规则体各向同性散射介质的辐射传递，结果发现，在某些工况下，有限体积 法的计算精度比离散坐标法的高。本文研制了贴体坐标系下有限体积法计算程序，研究 了不规则形体、三维椭圆柱腔体各向异性散射介质内的辐射传递，并与文献【3】进行了比 较；考察了计算精度选取对计算结果的影响。 2贴体坐标系下有限体积法 吸收、发射、散射性灰介质的辐射传递方程为 (i·V)，(r，i)= (1) 肌∞+口(，)驰搠+掣L。妣邶(j，i’脚7 立体角的离散方式与文献[31相同，天顶角与圆周角均匀划分，对式(1)积分 (2) k，L，协‘扁)幽扪7=b‘r(一filI+S7)d矿搠7 选用阶梯格式将控制体表面上的辐射强度与控制体中心的辐射强度关联起来 茹矗= ∑ 4巧+畦 (3) I=E，∥，N。S 离散方程组(3)采用迭代方法求解。令：，g为前一次迭代计算值，Er为计算精度，则 迭代收敛的判别式为 max峙昂l。怍1l。 (4) 直角坐标系(‘Y，z)转换为贴体坐标系(f，r／，f)的详细推导见文献【4】。 教育部高等学校博士学科专项基金资助 433 3结槊及分耩 3．1二维不惑掰瑟边形碹肉{#散鸯圣性余震豹辖菇传递 文献f3l{中篝了不攫刘甄边形腔体底强的笼因次热流圣，但来考虑教射。腔体内 参与性分成的漩度保持为咒，壁西为Ox黑体，吸收系数分别取O．1m一、Im…1和 lOm～。本文用贴体坐标系下有限体积法计算网一问题，空阐网捂为Ⅳ。×N。=10X10， 立体角网格^0x^乞=8x16，网格划分见闰1，无因次热流分布见图2．由图可见，本 文结聚与文献【31的结果吻合。 文【ll】 本文 =o．1 瞄l小娥则四