求解辐射传递的非结构混合有限体积有限元法研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编口-：063262 求解辐射传递的非结构混合有限体积／有’ ．。限元法 张琳，刘林华 (哈尔滨工业火学能源科学与工程学院，黑龙江哈尔滨 150001) Tel．：045lEmail：lhliu@hit．edu．cn 摘要：奉文给了一种适用于任意非结构网格的有限体积，有限元法的混合算法用于求解多维半透咧吸 收、发射、散射性灰矩形介质内的辐射传递。该方法使用有限元法进行角度离敬，有限体积法进行空 问离散。与基于辐射传递离散坐标方程的方法不同的足，该方法在迭代求解的过程中，针对每一个空 间体元，所自．角度方向的辐射强度同时耦合求Uj．通过两个算例验证了该解法的正确性． 关键词：自．限元法有限体积法辐射换热非结构网格 1引言 半透明介质内的辐射传热有广泛的工程应用。近几十年来，对于半透明介质内的辐 射换热，各国学者做了大量的研究，已发展了许多的辐射换热的计算方法，主要有蒙特 卡洛法、区域法、离赦坐标法、有限体积法等。目前发展的基于辐射传递方程的求鳃方 法主要包含两部分内容，一是角度离散，一是空间离散。有多种方法可以进行空间和角 度的离散，如离散坐标法角度项采用离散坐标离散，空间采用差分法离散；有限体积法 空间和角度均采用有限体积法离散。 在进行空间离散的过程中，由于其容易处理空间上不规贝wJJt,何形状、本身的守恒特 性的优点，有限体积法得到了广泛的应用。但对于角度离散，有限体积离散与离散坐标 离散常受到射线效应的影响。为了克服角度离散带来的射线效应，Coelhoill提出了一种 角度离散采用有限元法的方案，在空间结构化网格的基础上发展了一种有限体积，有限元 的混合方法用于辐射传递的求解。其初步结果表明，这种新的方法对于求解多维辐射传 递问题给出了很好的结果。但他给出的公式仅限于结构化网格。 本文在文献【l】的基础上，发展了一种非结构混合有限体积／有限元法米求解多维半 透明介质内的辐射传递，其中空间使用非结构化网格用有限体积离散，角度用有限元法 离散。最后通过两个算例验证该解法的正确性。 2数值方法 考虑充满吸收、发射、散射性介质封闭空腔内的辐射传递，辐射传递方程可以写为 12，31 (1) 对q-．不透明、漫射边界，给定如下边界条什： 563 (2) 。 ms)=气”)+孚dfllI,以U～s’№s，)耐万 ‘●‘’ 其中 角度项采用有限元近似，辐射强度，可以通过形函数表示为 Ⅳ 一 l(r，s)=Xl”(r)r胛(s) m=l 口 L 乎 ＼ L ＼黝