梯度折射率光纤热辐射特性的有限元分析研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083297 梯度折射率光纤热辐射特性的有限元分析 张琳，刘林华 （哈尔滨工业大学能源科学与工程学院，黑龙江 哈尔滨 150001 ） （Tel.: 0451Email: lhliu@hit.edu.cn ） 摘要：本文给出了柱坐标系下轴对称梯度折射率圆柱内的辐射传递方程。该方程易于处理具有 Fresnel 表面的圆柱形介质内的辐射传递。基于该方程，发展了一种分析半透明梯度折射率光纤热辐射特性的 有限元法。分析了等效发射率随着光学厚度的变化，并对折射率为均匀分布及线性分布的光纤热辐射 特性进行了比较。有限元法获得的结果与蒙特卡洛法得到的结果吻合，折射率分布的不均匀性使得介 质内部热流分布的非线性增强。 关键词：光纤 等效辐射特性 梯度折射率 有限元法 0 引言 光纤在其拉拔过程中一般会达到 1000-2000K 的高温，此时对光纤与系统换热的研 究中热辐射过程尤其重要。光纤的热辐射特性，如吸收率和发射率，是了解光纤辐射换 热特性的重要参数。同时对于光纤拉拔过程的数值模拟，也经常采用基于光纤等效热辐 射特性的简化辐射换热模型来进行光纤辐射换热的计算。由此可见，光纤热辐射特性对 于光纤相关传热分析具有重要意义。 通常光纤的直径为 125µm 左右，远远大于 1000K 或更高温度时热辐射能集中分布 的波长范围 0.5-30µm 左右，因而辐射传递方程可用于描述光纤内部的热辐射传递。刘 林华等[1]给出了一种射线踪迹法并研究了均匀折射率半透明圆柱的吸收特性。基于文献 [1]的射线踪迹法，Tian 等[2]研究了表面包覆薄膜的均匀折射率光纤的热辐射特性。对于 梯度折射率光纤，由于热射线沿曲线传播，其内部的辐射传递采用射线踪迹法来求解比 较困难。因而文献[2]中采用的方法只能用于均匀折射率光纤热辐射特性的分析。目前还 没有研究梯度折射率光纤热辐射特性的工作发表。 本文发展了一种分析梯度折射率光纤热辐射特性的有限元法。首先推导了圆柱坐标 下梯度折射率介质内的辐射传递方程，然后给出了基于该方程的最小二乘有限元法。最 后通过算例对本文的方法进行了验证并对折射率梯度分布对光纤热辐射特性的影响进 行了分析。 1 辐射传递模型 如引言所述，这里所研究的梯度折射率光纤，其内部的热辐射传递可以通过梯度折 射率介质内辐射传递方程来描述。梯度折射率介质内辐射传递方程的一般形式沿射线轨 迹可以写为[3, 4] ： 2 d I (r, Ω)  2 κs ′ ′ ′ n ( )I (r, Ω) n I (r) I (r, Ω) (Ω, Ω)d (1) + + =+ Φ Ω κ κ κ d s  n2  a s a b 4π ∫4π 资助项目：国家自然科学基金项目(No.50620120442 其中κ 、κ 分别为介质的吸收、散射系数，n 为介质折射率，Φ为散射相函数。式中d/ d s a s 称为流算子，为了采用有限元法进行求解，需要将其在对应的坐标系下的展开，以得到 对应坐标系下的辐射传递方程。这里考虑图 1 中以ρ,ψ, z 定义的圆柱坐标系，其中 e ,e ,e 分别为对应单位向量，光纤的轴向沿z 方向。在该圆柱坐标系中，辐射传递方 ρ ψ z 向向量Ω与局部坐标系位向量e 和e 的夹角分别为θ和ϕ 。