[自然科学]蒙特卡洛方法.pdf

热辐射传输中的 蒙特卡洛方法 航空航天热物理所 2007 年9 月21 日 1 蒙特卡洛方法概述 1 蒙特卡洛方法概述 蒙特卡洛方法是一种随机模拟方法。将 其用于辐射传热计算时，其基本思想： 将辐射能量看成由大量独立的光束（光子） 组成 将复杂的辐射传递问题分解为发射、反射、 吸收、散射等独立过程。每一光束在系统 内的传递过程，由一系列随机数确定 跟踪一定量的光束，可得较为稳定的统计 结果 1 蒙特卡洛方法概述 1.1 辐射传递因子 辐射传递因子RDij 的定义为：在一个 换热系统中，由单元i 辐射出去的能量被 单元j 吸收的份额。比较角系数Xij 的定义： 在一个换热系统中，由单元i 辐射出去的 能量到达单元j 的份额。 1 蒙特卡洛方法概述 辐射传递因子用于能量方程 由N 个介质单元、M 个表面单元组成的热 辐射系统中，介质单元i 的能量方程为 N 4 4 4σκ V T ∑4σκ V T RD a i i a j j ji j 1 M 4 +∑F σε T RD k k k ki k 1 1 蒙特卡洛方法概述 辐射传递因子用于能量方程 壁面单元i 的能量方程为 N 4 4 F εσT ∑4σκ V T RD i i i a j j ji j 1 M 4 +∑F σε T RD k k k ki k 1 1 蒙特卡洛方法概述 如何求辐射传递因子 （蒙特卡洛M-C法） 首先确定光子的发射点。需要3 个参数 R , R , R x y z 再确定2个方