梯度折射率介质内辐射传递方程的P1近似研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093378 梯度折射率介质内辐射传递方程的P 近似 1 刘立君，刘林华 （哈尔滨工业大学能源科学与工程学院，黑龙江 哈尔滨 150001 ） Tel.: 0451，Email: lhliu@hit.edu.cn 摘要：考虑一个由不透明和漫射壁面包围的二维方形半透明介质内的辐射传递问题。利用P 1 近似对均 匀折射率及梯度折射率介质内的温度场进行了求解，求解结果分别与蒙特卡罗法进行了对比分析。研 究结果表明，介质的光学厚度对P 1 近似的求解结果影响较大，而且由于梯度折射率介质内光线的曲线 轨迹导致热流分布不均匀，致使P 1 近似的求解精度变差，但当介质为各向异性散射介质时，在特定情 况下散射特性与折射率变化所造成的计算误差可以部分抵消，从而在某些情况下可以提高P 1 近似的计 算精度。 关键词：P 近似；梯度折射率；蒙特卡罗法 1 1 前言 由于材料的结构特性以及可能出现的对温度的依赖性，介质的折射率可能是空间位 置的函数。在这种情况下，介质内光线不是沿直线传输，而是沿曲线传输[1,2] 。半透明梯 度折射率介质内的辐射传输在生物组织光学成像、等离子体诊断、玻璃加热及光学图层 等方面具有重要的研究意义，已经引起了研究人员的广泛关注。 由于很难对弯曲光线进行跟踪，因此梯度折射率介质内辐射传输方程的求解与均匀 折射率介质内的求解相比要困难得多。在辐射传热领域，蒙特卡罗法（MCM ）被普遍 认为是一种精确的求解方法，对于多维辐射传热问题，蒙特卡罗法得到的数值解经常作 为基准解。但正是基于这种复杂的曲线跟踪技术，蒙特卡罗法由于其非常耗时的计算过 程而很难得到普遍的应用。而P 1近似法计算速度快，而且在某些特殊的应用场合具有足 够的精度，因此，很多研究人员对梯度折射率介质内辐射传递方程的P 1近似解法投入了 大量的精力[3,4] 。Modest求解了均匀折射率介质的P 近似方程，并从物理学原理的角度分 1 析了该方程的物理意义[5] 。Bhuvaneswari和Wu利用不同的扩散近似方法求解了各向异性 散射介质内的瞬态辐射传递方程，并将所求得的时域反射率和时域透射率与离散坐标法 进行了比较[6] 。 本文推导了梯度折射率介质内辐射传递的 P 1 近似方程，并分别求解了二维梯度折 射率介质及均匀折射率介质内温度场的分布情况，求解结果与 Liu[7]提出的蒙特卡罗曲 线射线跟踪法进行了对比分析，并分析了介质的光学厚度、折射率的变化以及散射特性 对P 1 近似求解精度的影响。 资助项目：国家自然科学基金(No. -1- 2 数学公式 2.1 梯度折射率介质内的稳态辐射传递方程 [8] 梯度折射率介质内的稳态辐射传递方程可以写为 ： ( , , ) 1 ∂ ( , , ) ( ) 2  s ⋅∇I r s t +2n2 sin θ θ∂ {I r s t ξs −k ⋅∇n } 1 ∂ ( , , )  2  ( ) ( , , ) + ⋅∇ + +