各向异性散射梯度折射率介质中圆柱坐标系下DRESOR法的求解.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：123096 各向异性散射梯度折射率介质中圆柱坐标 系下DRESOR 法的求解 张险，程强，周怀春 （华中科技大学煤燃烧国家重点实验室，武汉430074 ） (Tel: 0278507, Email: chengqiang@mail.hust.edu.cn) ： DRESOR Distribution of Ratio of 摘要 本文提出了各向异性散射梯度折射率介质中圆柱坐标系下 法（ Energy Scattered by the medium Or Reflected by the boundary surface ）的求解方法。以三维圆柱形非黑壁 面、吸收、发射、各向异性散射介质为例，采用圆柱坐标系下DRESOR 法计算其任意方向的出射辐射 强度，并且和文献中的计算结果进行了对比，进行了详细的分析。结果表明，圆柱坐标系下DRESOR 法，在处理三维梯度折射率介质中的辐射换热分析中是一种行之有效的方法。 关键词：三维圆柱坐标系；DRESOR 法；各向异性散射；梯度折射率 0 前言 热辐射是能量传递的一种重要方式，是能源动力以及化工冶金设备结构性能、工艺 性能、使用性能设计的关键[1] 。DRESOR 法应用蒙特卡洛法计算了颗粒介质的散射作用， 能够以较高的空间方向分辨率计算封闭腔内任一点任一方向的辐射强度。传统的 DRESOR 法都是基于直角坐标系的，划分矩形网格，通常应用于矩形介质中。然而很多 情况下，工业实际辐射换热研究分析中会面临一些圆柱形的介质。例如光纤的拉丝成型 过程，火箭等飞行机器的尾部喷焰等。对于这些圆柱形的介质，建立直角坐标系不易精 确的描述圆柱形几何边界，而采用圆柱坐标系可以很好的解决这一问题。通常在参与性 介质内辐射传递的研究中，总是假设介质的折射率均匀分布，而实际介质的折射率受其 组分、温度、浓度等的影响，通常是非均匀分布的，即通常所说的梯度折射率介质。所 以，研究梯度折射率介质内的辐射换热对于更加准确的模拟实际的辐射换热过程具有重 要的指导意义。本文在传统的直角坐标系下DRESOR 法的基础上开展了三维圆柱坐标系 下各向异性散射梯度折射率介质中DRESOR 法的研究。 1 梯度折射率介质中圆柱坐标系下能束的追踪和DRESOR 数的计算 1.1 体积微元能束的发射方向 z P (r ,ϕ, z ) z r y ϕ x 1 图 圆柱坐标系示意图 本文建立的圆柱坐标系如图1 所示，方位角 的范围为[ ] 。分别沿着R 、 、 ϕ 0, 2π ϕ Z 方向划分网格。体积微元向四周发射能束时，可以认为在各个方向上的机会是相等的， 基金项目：国家自然科学基金（No.51021065 ），中央高校基本科研业务费专项 资金资助（No.2011TS143 ） [2] 每个微元体以发