参与性介质方向辐射的广义多流法研究.doc

哈尔滨工业大学工学硕士学位论文 场角也很小，只能探测到很小角度内投射到探测头上的能量。MC方法可以 解决这个问题，但计算效率很低。 近年来，已经发展了很多求解参与性介质内辐射传输问题的数值方 法，包括FVM、DOM、MC、DTM、ZM等，这些方法在求解辐射换热问题 上均有优点也有缺点。当只求解热流和能量问题时，各种方法之间的差别不 会很大，但对于某一确定方向辐射强度的计算问题，很多方法的求解出现困 难，每种方法的求解结果可能大不相同。例如，离散坐标法中，其离散方向 大都是固定的，只能精确求解某些满足固定离散方向的辐射强度，方向不能 任意调节；BMC方法理论上可以求解任意方向的辐射强度，但其缺点是计 算速度较慢（虽然相对于MC方法已有很大改善），迫切需要一种兼顾计算 精度和计算效率的数值计算方法，以满足目标探测、大气遥感、燃烧诊断等 领域的实际需要。 考虑到实际中的定向辐射探测问题，本文的目的是提出一种能够满足 需要的简单快捷的、并且具有普适性的计算模型。需要寻找一种能够快速准 确的求解方向辐射强度的方法，因此本文提出广义多流法模型，可以较为精 确的计算方向辐射强度和表观定向发射率，且计算效率较高。 辐射传递求解方法的概述 到目前为止，针对辐射传输问题的求解已提出了许多方法，传统的有： 区域法、蒙特卡洛法、离散传递法、离散坐标法、球形谐波法、有限体积法 等。各种方法都各有其优缺点，通过这些方法所能获得的解的精度及提供信 息的详细程度也不相同。从基本原理来看，求解辐射传输的数值方法大致可 以分为两类：基于能量追踪的求解方法和基于辐射传输方程的求解方法。蒙 特卡洛法和射线踪迹法是第一类方法的典型代表。而区域法、离散传递法、 离散坐标法、球形谐波法、有限体积法等属于第二类方法。其中，第二类方 法又可以分为：基于辐射微分方程的求解方法和基于辐射积分方程的求解方 法。除区域法和离散传递法外，大部分第二类方法都是基于辐射微分方程的 求解方法。 - 2 - 哈尔滨工业大学工学硕士学位论文 区域法 区域法(Zone Method)是由Hottel和Cohen[6]最先提出的，该方法实质上 是计算表面间辐射换热的净辐射法的一种扩展[7-8]。它将封闭空腔及其表面 分成若干体元和面元，并假定每一区域的温度和物性均匀一致。进而，对每 个区域可以写出辐射能量平衡方程式，计算每两个区域之间的直接辐射交 换，最后得到每个区域的净辐射热流。由于需要计算并储存大量的交换面积 参数，所以这种方法比较复杂费时。因而区域法常用来计算封闭空腔几何形 状不是很复杂且体积区域和表面区域可以划分得较少的辐射换热问题。区域 法对无散射的辐射问题有较好的计算精度，但还很难处理非均匀、各向异性 散射介质内的辐射换热问题。 近十余年来，区域法作为一种基本的热辐射传递计算方法得到了广泛的 应用[9-10]，其方法本身也得到了改进和完善。例如：Saltiel等[11]用非结构化 网格对区域法进行改进；针对耗时量大，提出了并行处理的方法[12]。 离散传递法 离散传递法最先由Lookwood和Shah等[13]人提出，该方法同时兼有区域 法和蒙特卡洛法的特点。其主要思想是使用边界网格面作为辐射的吸收和发 射源，将边界网格面上向半球空间发射的辐射能离散成有限束能束并形成相 应的特征射线，沿特征射线积分辐射传递方程直至到达另外的边界面。 若设由图 1-1 边界面 ?Ai 发射的一根特征射线 l 在进入某控制体 P 时的 谱带辐射强度为 Ik, P ，离开该控制体的谱带辐射强度则为 Ik, P+1 。沿射线在 控制体 P 上积分辐射传递方程直至最终到达边界面，忽略散射，则射线穿过 内部控制体时谱带辐射强度的变化可由下式递推求得： n2 B σ T 4 ( ) ( ) k k,T ? ? P Ik, P+1 Ik, P exp a, k s 1 exp a, k s = ?κ δ + ? ? ?κ δ ? (1-1) P π - 3 - 哈尔滨工业大学工学硕士学位论文 网格界 线 边界网 格面元 图1-1 离散传递法的计算模型 式中：δ s 为特征射线 l 在控制体P 中的行程；nk 为?λk 谱带内介质的 ∞ 折射率；TP 为控制体 P 的温度； = ∫ ( ) ∫ ( ) 是 Bk, T Ebλ T dλ Ebλ T dλ P P P 0 ?λ k T 温度下?λk 谱带内黑体辐射能占总辐射能的份额。 P 离散坐标法 离散坐标法(Discrete Ordinate Methods,DOM)又名 S 方法，与球形谐波 N 法类似，也是将传输方程（对灰体或基于光谱）转化为一系列偏微分方程组, 其理论上可应用于任意阶数和精度。该方法基于对辐射强度的方向变化进行 离散，通过求解覆盖整个4π