超短脉冲激光在二维非均匀介质内的瞬态辐射传输研究.pdf

中国一J：程热物理学会 传热传质学 第十二届年会 编号：063203 超短脉冲激光在二维非均匀 介质内的瞬态辐射传输 阮立明‘ 安巍 谈和平 (口A尔滨下业大学航空航天热物理研究所，黑龙江， 哈尔滨150001) Tel：0451 E-mai1：ruanlm@hit．edu．cn 摘 要：本文采用最小二乘有限冗法模拟丫超短脉冲激光局部入射条件下，具钶高散射核的二维啦均 匀介质内的瞬态辐射传输，研究了彳i同边界位置上反射和透射信号随时间的变化情况。研究表明，对 于具自．高散射核的1F均匀介质，其反射信I=}存在明显的双峰现象，这种现象可以明显的标示ilj光学物 性发生突变的位置。瓶且，在边界处还存在“双峰合并点”，即反射信号出现最大值的点，该点也足双 峰现象出现与消失的临界点。而对于透射信弓，则没有明显的反映出光学物性的突变。此外．本文还 提出了“穿越距离”的概念，研究表明：在两倍于“穿越距离”之后，各点的反射和透射信号的哀减 速率几乎相同。凶此，在此时间后采样所能得到的介质信息较少。 关键词：超短脉冲激光，瞬态辐射传输，最小二乘有限元 1．引言 在通常的应用中，往往将辐射传输假定为稳态过程，不考虑其暂态效应。近期，随 着超快激光的飞速发展，及一些相关的新兴应用的出现，使这种假设不再成立。在这些 新兴的应用中，超快激光的脉冲时间总是在皮秒或飞秒量级，与系统中辐射传递的时间 量级相当。因此，不能忽略辐射的暂态效应。超快激光的应用从金属材料的加工剑粒子 探测，从人气红外遥感到医学光学成像，范围涉及许多前沿技术领域。由于诱人的应用 前景，医学光学成像成为这些领域中最引人注目的技术。然而，由于瞬态辐射的复杂性， 在以往光学成像的研究中，人们总是利用简化的扩散模型进行瞬态辐射的求解。Mitra 和Kuamr[1]研究了几种不同的扩散近似模型的计算精度，他们发现：不同的模型计算所 得到的结果中，散射的传播速度，透射辐射强度的火小，光学信号的时域形状都差异很 大。扩散近似方法对于非散射或者低散射区域内的辐射传输有着很大的误差。Elaloufi[2】 等研究者也表明扩散近似不能够准确模拟光学薄介质内的暂态甚至是稳态辐射传输。即 使住光学厚条什．卜．，扩散近似也不能准确模拟入射光源周围的辐射传输。进一步的研究 表明，生物体内的辐射传输并不总是遵循光学厚假设，一些特殊的生物组织或细胞聚集 体在近红外波段的光学物性并非光学厚。例如，在对人脑组织的光学成像中，人们发现 分布在脑壳周围的脑髓液对光学成像的分辨率有着很大的影响。研究表明：这层薄的液 体有很小的吸收和散射系数。这就要求在模拟中必须应用更准确的辐射传输模型进行精 确的模拟。冈此，开发高精度的，能够和原有的光学成像技术有效融合的瞬态辐射传输 有限元模型有着理论和应用上的重要意义。 在处理超短脉冲激光引起的辐射传输模拟中，稳态辐射传递方程不可能提供准确的 结果。这是由下描述瞬态辐射波传递的微分方程已经变为了双曲型的偏微分方程，辐射 强度同时是时间，空间，方向，波K四种参数的函数。冈此，瞬态辐射传递的模拟研究 成为近十年来光学和计算辐射学的研究热点。近十年中，已经有许多理论研究承l数值方 资助项日：困家自然科学埚企资助(NoNo 法被提出。通常求解瞬态辐射传递的方法包括：蒙特卡洛法(MCM)，积分方程法(IE)， 等将最小二乘有限元法(LSFEM)应．}}j丁瞬态扩散辐射的模拟【3】，其后又将这一模型扩 展到超短脉冲激光入射条件下的一维瞬态辐射传输问题的求解[4】。他们的研究表明，最 小二乘有限元法能够准确地求解瞬态辐射传输问题。Lu和Hsu[5]朋反蒙特卡洛法分析了 一维分层1F均匀介质内短脉冲激光的辐射传输，并提出了反射信号的“双峰现象”。Kumar P：A格式的离散坐标法分析了二二维非均匀介质内短脉冲激光的辐射传 和Mitra[6]应用P 输。然而，在他们的研究中，脉冲激光总是覆盖予介质边界的一侧。更具实际意义的局 部边界入射条件，仍没有被系统的研究。