求解方向辐射强度的等效源项热流法.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：073252 求解方向辐射强度的等效源项热流法木 阮立明 齐宏 宋跃超王英立 (哈尔滨工业大学航空航天热物理研究所，黑龙江，哈尔滨150001) Tel：045l E．mail=qmong@hit．edu．cn 摘要：本文提出一种能够快速准确模拟任意方向辐射强度的等效源项热流法模型。该模型采用迎风格 式进行数值差分，将各向异性散射近似处理为与六个正交方向量相关的等效“源项”，建立六流模型， 计算该等效源项，然后根据源项不变的原则，建立任意方向的计算摸型。通过与反向蒙特卡罗法方法 (BMC)计算结果的比较，表明两者吻合较好，而计算效率远高于BMC方法。 关键词：等效源项热流法；任意方向辐射强度：反蒙特卡罗法 ’ 1 引言 对于辐射传输过程的数值模拟，国内外学者已提出了许多计算方法，如：区域法、 离散传递法、蒙特卡罗方法、反蒙特卡罗方法、离散坐标法、有限体积法等【lJ。但针对 任意方向的辐射强度计算而言，这些方法都具有局限性：反蒙特卡罗方法是公认的高精 度的可作为基准的方向辐射强度模型【l】，但需要大量的统计光线模拟，如果需要同时模 拟多个方向辐射强度：其计算速度相对较慢；离散坐标法(DOM)是一种具有潜力的 计算速度较快的方法，然而由于离散方向的固定，无法准确求解任意方向辐射强度而只 有通过捕值计算进行近似处理；华中科技大学的李宏顺等【z’3】提出的“带无限小加权离 散坐标格式”的DOM方法可求解任意方向的辐射强度，但如果需要计算很多不同方向 辐射强度，由于每个独立方向均需要进行重复的离散辐射传输方程求解，其计算量势必 将会很大。 本文提出一种基于辐射“源项”的新型六热流法，在保证计算精度的前提下，可快 速的求解任意方向的辐射强度，同时可以兼顾处理各向异性散射问题。其基本学术思想 是：将各向异性散射近似处理为与六个正交方向量相关的等效“源项”，假设各微小控 制体单元内的等效“源项”与所选计算方向无关【4J，即每个微元的源项是一不变的常标 量，通过求解六个正交坐标方向的辐射强度，计算出每一控制体单元的等效辐射“源项”， 然后根据“源项”不变假设，采用迎风差分格式对所求方向的辐射传递方程进行离散， 迭代求解沿该方向的辐射强度分布。该方法计算程序实现方便，可求解任意复杂几何形 状的辐射换热问题，计算结果与反向蒙特卡洛方法(BMC)模拟的结果基本吻合，计算 速度快，效率远高于BMC法，是一种适用于求解方向辐射的传输计算模型。 2理论模型 本文所提出的源项热流法模型的求解步骤分为两个部分：正交方向六流模型的求解 和任意方向辐射强度的求解。具体而言，首先建立沿六个正交坐标系轴方向的六流模型， 计算等效辐射‘：源项”，然后根据“源项”不变假设j将辐射“源项”代入所求任意方 向的辐射传输模型，通过数值迭代得到该方向的辐射强度。 2．1正交方向流模型 根据参与性介质内的辐射传递方程： ‘项目资助：哈尔滨工业大学“新世纪优秀人才支持计划”基金 708 (1) 以工轴为例，假设x正向分量的所有辐射强度与角度无关并且为一常数，x负向分量的 所有辐射强度与角度无关并且为一常数，。即在每个方向上的辐射强度为各向同性。则正 交坐标轴方向六流法模型如下： ‘ (2) ±笺}：一￡+国饿+∞6e+(1一国)，b+∞s(c+厶+，；+‘) ． dr朋 其中历上而上多，厂、6代表从正向(前向)、负向(后向)进入控制体的辐射散射份额， 用J表示从侧向进入控制体的散射份额，表达式为：厂=告f①(秒，2)嘲∥， 当历取