计算热辐射学的进展.pdfVIP

讲 学 叛 2009年 第54卷 第18期：2627～2637 《中国科学 》杂志社 WWW．scichina．com csb．scichina．com SCIENCElNCHlNAPRESS 计算热辐射学的进展 谈和平 ，刘林华① ，易红亮①，赵军明①，齐宏①，谭建宇② ① 哈尔滨T业大学能源科学与工程学院，哈尔滨 150001 ② 哈尔滨工业大学汽车工程学院，威海 264209 联系人，E—mail：tanheping@hit．edu．cn；lhliu@hit．edu．cn 2009-07—13收稿，2009—08—26接受 国家自然科学基金资助项 目(批准号50620120442) 摘要 描述辐射能传输过程的基本方程为辐射传递方程．在三维半透明介质中，辐射强度是七 关键词 维变量的函数，多数情况下只能通过数值计算的途径进行近似求解．数值模拟方法已成为半透 热辐射 辐射传递方程 明介质 内热辐射传递理论研究和工程应用的重要手段．本文对国内学者在计算热辐射学方面 计算方法 取得的主要研究进展进行了综述，提出了今后计算热辐射学研究的重点． 在半透 明介质内热辐射传递的理论研究和工程 计算的途径进行近似求解 ．随着高速计算机的普及， 应用中，描述辐射能传输过程的基本方程为辐射传 目前数值模拟方法已成为半透明介质 内热辐射传递 递方程 ．均匀折射率灰介质内的辐射传递方程可写 理论研究和工程应用的重要手段． 为 目前 已发展起来的辐射传递方程数值求解方法 主要可以分为以下两类[31：1)基于射线跟踪的方法； + ．v，(，，，f)：一( -I- )( ，f) df 2)基于微分形式辐射传递方程全局离散的方法．第 产 一 类方法一般需要通过跟踪光束传播轨迹来进行求 + b(，，f)+÷LJ，(，，s，f)(，，)d，， (1) 解，如区域法、蒙特卡洛法和离散传递法．由于需要 47[ 4 式中r为空间位置矢量，s为辐射传递方向单位矢量， 进行繁重的光束传播轨迹计算，这类方法计算速度 为真空中的光速，n为折射率，，为辐射强度，，h为黑 极慢，不适合用来与热流体 问题 中的其他过程进行 体辐射强度， 为吸收系数， 为散射系数，鳓 散射 耦合计算．第二类方法的求解过程类似于一般偏微 相函数，t为时间．辐射传递方程描述 了辐射强度在 分方程的离散和求解，如离散坐标法、有限体积法、 时问、位置空间和角度空间上的变化．方程(1)中右边 有限元法、谱元法、无网格法等．这类方法计算速度 最后一项包含对整个 4兀立体角空间所有方 向上辐射 较快，适合与对流和导热过程耦合计算，是 目前计算 强度 的积分，因而在散射性介质 中所有方向上 的辐 热辐射学研究的重点． 射传递方程耦合在一起，需联立求解． 近10年来，在 国家 自然科学基金的资助下，国内 辐射传递方程为一积分一微分方程．若将 看作 许多高校对与计算热辐射学相关的课题进行了广泛 速度矢量，则s．VI(r，，t)相当于对流项，因而辐射传 的研究．本文对 国内学者在与辐射传递方程数值模 递方程可 以看作一个特殊的对流扩散方程．由于辐 拟相关的计算热辐射学研究中所取得的进展进行回 射强度是波长、时间、空间位置坐标、角度方向的函 顾，包括：1)辐射传递方程的数值解法，如射线踪迹 数，在三维半透明介质中，辐射强度是七维变量的函 一 节点分析法、有限元法、谱方法、谱元法、问断元 数，一般难以获得理论解，多数情况下只能