梯度折射率介质中瞬态辐射问题的DRESOR法求解.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：123094 梯度折射率介质中瞬态辐射问题 的DRESOR法求解 1 1* 1,2 王志超 ，程强 ，周怀春 1华中科技大学煤燃烧国家重点实验室，430074 2清华大学热能动力仿真与控制研究所，100084 *Tel ：0278507，Email ：chengqiang@mail.hust.edu.cn 摘要：近年来，梯度折射率介质中瞬态辐射传递问题引起了人们越来越多的关注，而求解该问题的研 DRESOR ( ) 究相对较少。本文将 法被介质散射或壁面反射的能量分布方法 推广到求解梯度折射率各向 异性散射介质中的瞬态辐射传递问题。本文考察了一维无发射，吸收，各向异性散射梯度折射率介质 中瞬态辐射传递问题，对比结果显示DRESOR 法具有很好的求解准确性，同时讨论了光学厚度、散射 反照率、散射相函数以及梯度折射率对反射率和透射率分布的影响。 关键词：梯度折射率；瞬态辐射；DRESOR 法 0前言 由于激光脉冲技术的高速发展，瞬态脉冲已广泛应用于诸多领域。诸如雷达和光学 通信[1],海洋激光雷达[2],远距离遥感[3],热疗[4],光动力疗法[5],医学诊断[6],激 光加工的微结构材料[7]，等等。这些应用反过来引起并促进了科研工作者对瞬态辐射 传递过程的研究。瞬态脉冲可持续时间在皮秒甚至飞秒量级，辐射传递控制方程必须考 虑瞬态项，也即，辐射传递方程转化成瞬态辐射传递方程，这使得求解更加复杂。近年 来，许多科研工作者在此领域开展了很多科研工作并应用许多方法来求解瞬态辐射传递 问题，如，积分方程法[8-10]，蒙特卡罗法[11-13]，有限体积法[14-17]，有限元法[18, 19]， 逆蒙特卡洛法[20-22]，DRESOR 法[23, 24]以及离散坐标法[25, 26]等。 需要指出的是，上述方法均针对均匀折射率介质中的瞬态辐射传递问题。然而，Guo 等[13]发现介质的梯度折射率会对透射脉冲形状的预测产生强烈影响。Liu 和Hsu[27]注 意到，基于介质折射率均匀假设进行数值求解得到的瞬态结果与实验结果相比，在t=0 处会出现一个偏移，其中可能的原因是忽略了梯度折射率的影响。由于具有广泛的应用 背景，近些年来，梯度折射率介质中瞬态辐射传递问题引起了人们越来越多的关注。 Ferwerda[28]推导了变折射率散射介质的辐射传递方程，但并不满足能量守恒；接着 Tualle 和Tinet[29] ，Khan 和Jiang[30]基于Ferwerda 的方程推导了两种不同形式的弥散 近似形式；Shendeleva[31]用欧拉方程处理得到变折射率介质中的瞬态辐射传递方程； Premaratne 等[32]得到了各向同性散射变折射率介质中光子传递方程。然而，致力于求 解梯度折射率介质中瞬态辐射传递方程的研究相对较少。Wu[33]将离散坐标法推广应用 于求解短脉冲入射到变折射率介质的问题中；Liu 和Hsu[27]用间断有限元法分析了该问 题；Wu [34]接着用蒙特卡洛方法应用于该问题的研究；Bhuvaneswari 和Wu[35]对比了 求解梯度折射率介质中瞬态问题的P-1 ，P-1 抛物线近似，P1/3 和双热流法。Wang 和Wu[36] 用一阶离散坐标法和改进的离散坐标法来求解弥散基质中变折射率瞬态辐射传递问题。 介质的折射率变化时，辐射能量传递沿着由费马定理确定的弯曲轨迹传递且传递速 度随着介质折射率的变化而变化，使得该问题的求解非常复杂。本文中，Zhou 等[37] 提出的DRESOR 法(被介质散射或壁面反射的能量分布方法)被推广到求解梯度折射率 基金项目：国家自然科学基金（No.51021065 ），中央高校基本科研业务费专项 资金资