【精选】旋转对称椭球颗粒沿轴脉冲光声波波形的几何计算.pdfVIP

第3O卷 第 3期 计 算 物 理 CHINESEJOURNALOFCOMPUTA~ONALPHYSICS Vo1．30．No．3 2013年 5月 May，2013 文章编号：1001．246X(2013)03—0409-06 旋转对称椭球颗粒沿轴脉冲光声波 波形的几何计算 李 勇， 方 晖 (南开大学现代光学研究所 光学信息技术科学教育部重点实验室，天津 300071) 摘 要：针对旋转对称椭球颗粒脉冲光声波的产生，推导沿其对称轴方向的波形的解析表达式．结果表明，无论是 长椭球还是扁椭球颗粒，沿轴波形相对于球形颗粒所产生的N形波都有一定的形变，并且在近场条件下显示出能 反映椭球颗粒形状与大小的明显特征．因此通过对脉冲光声波的测量和分析，能对旋转对称椭球颗粒进行形貌判 定． 关键词：脉冲光声波；旋转对称椭球颗粒；光声成像 ；光声效应 中图分类号 ：0426．3 文献标识码 ：A 0 引言 近十年来，基于短脉冲激光激发脉冲光声波的光声成像技术在生物医学领域得到了迅猛发展…．特别 是在生物组织血管成像方面，其显示了其它生物医学影像技术所不具有的独特优势 ]．由于 自然的几何形 状，血管往往被考虑为光吸收圆柱体．这不仅体现于通常的模拟实验上 “，而且也反映在理论计算分析 中 [引 ． 但是光声成像近年来的一个重要发展端是提高成像的分辨率 ，如：光学分辨率光声成像 已经能够成 像单个的细胞核 、血红细胞和更小的亚细胞器 ；而光声流式细胞仪也 已能够检测出单个的肿瘤细胞、血 红细胞，以及纳米颗粒 ]．这就使得对脉冲光声波产生机理的研究不能停留在血管尺度 ，而要跟进到单个颗 粒水平． 值得注意的是如上所述的各种微颗粒并非是球形的，而是更接近于旋转对称椭球 (如血红细胞接近于 旋转对称扁椭球)，因此有必要发展相应的理论，至少能计算旋转对称椭球颗粒的脉冲光声波．据我们 了解 ， 这方面的理论还未见报道．目前研究比较充分的是液体 中的球形颗粒，已考虑的情况包括：与周 围液体声阻 抗相匹配的球形液滴 ¨ ，与周 围液体声阻抗不匹配的球形液滴 ¨¨，以及不可压缩的固体微球 ¨． 我们考虑相对简单的情况，即椭球颗粒具有与周围液体相同的声阻抗，集中研究椭球颗粒外部在旋转对 称轴上的脉冲光声波．通过推导相应的解析表达式并进行计算，结果表明：不论是对于长椭球或者是扁椭球 颗粒，其脉冲光声波的沿轴波形都不同于球形颗粒所产生的N形波，而是有所形变，我们特定义其为 “变形 N形波”．通过对这些计算结果的近一步分析，我们发现，相对于N形波的中心零点位置 ，长椭球颗粒的变形 N形波的中间零点位置在时间上略微超前 ，而扁椭球颗粒的在时间上略微延后 ，这些时间偏离与椭球的焦距 有关．因此在近场条件下，通过检测和分析沿轴脉冲光声波 ，可以获知椭球颗粒的具体几何形状及大小．所发 展的理论计算模型也适用于对聚焦脉冲光产生光声波的情况，因为在光焦点处的光能量分布接近于旋转对 称长椭球而并非球形 ¨． 收稿 日期 ：2012—09—06；修回 日期 ：2012—12—11 基金项 目：国家 自然科学基金和教育部新世纪优秀人才支持计划(NCET一1O一0502)资助项 目 作者简介：李勇 (1980一)，男，博士生，主要从事光声探测技术研究；E-mail：qq865704819@163．{2om 41O 计 算 物 理 第3O卷 1 计算方法 所考虑的条件是在透明液体中的光吸收椭球颗粒，并且颗粒的密度和声速均与周围液体的相同(这样 颗粒与液体的声阻抗相 同)．计算方法是在标准的格林函数求解法 基础上的一种扩展，简要描述如下． 假定光源脉冲为理想的时间6函数，颗粒受到均匀照明(这要求颗粒的光吸收系数远小于其尺寸的倒 数)，因此颗粒即时均匀受热，这时的热产生函数(单位时间单位体积内产生的热能)为