激光蚀除生物组织时的传热过程分析与数值模拟研究.pdf

中国工程热物理学会 传热传质学 学术会议论文 编号：073272 激光蚀除生物组织时的传热过程分析与 数值模拟 张纪庄，张学学 (热科学与动力工程教育部重点实验室，清华大学热能工程系，北京100084) 电话：01 摘要：本文对连续和脉冲激光蚀除生物组织时的传热过程进行了理论分析和数值模拟。数值结果表明， 即使是穿透深度较小的c0：激光蚀除，表面组织碳化前组织也具有较厚的脱水层，应该将激光辐照当 作体积内热源而不是表面热流处理。数值计算还得到厂备蚀除阶段的能量分布、温度分布、各组织层 边界位置随时问的变化规律。 关键词：激光蚀除：多层结构模型：stefan问题；动态光热模型；生物传热 1引言 激光对生物组织的蚀除在临床上具有高精度、微侵袭、手术时间短等优点，已广泛 应用于对皮肤、牙齿、心血管等疾病的治疗，但人们对蚀除过程中的热作用规律认识还 不足。以往研究者所提出的模型往往过于简化H“，忽略了组织对激光的容积吸收特性、 组织物性参数动态变化等因素对蚀除的影响。张纪庄等p山j在№kenzie三层结构模型¨j 的基础上，初步提出了描述激光蚀除生物组织过程的多层结构模型，在以下几个方面进 行了尝试性改进：(1)考虑激光能量在组织内的容积吸收特性，组织中存在具有～定厚 度的汽化区域；(2)考虑了蚀除时组织水分汽化产生的水蒸气流过高温多孔的干化层和 碳化层组织时与这些固体组织的对流换热；(3)考虑了蚀除过程中光热物性参数的动态 变化对激光能量分布和温度分布的影响。本文将进一步完善对此模型的理论分析。此外， 由于涉及具有多个待定边界的Stefan问题和需要耦合求解分布不均且动态变化的物性参 数、激光能量分布和温度分布，该模型的数值求解比较困难，本文也将对其做重点讨论。 2激光蚀除生物组织的多层结构模型 文献f71将激光蚀除生物组织的过程分为初始阶段和准稳定蚀除阶段两个阶段。在初 始阶段，组织吸收激光能量而温度升高，．当温度高于某阂值时，表面组织开始汽化脱水， 脱水后剩下的干组织的温度继续升高，当达到碳化温度时，于组织发生碳化反应，析出 碳黑，当达到燃烧温度时，表面组织开始燃烧而被去除。与此同时，更深层组织的温度 也在透射的激光能量和表面组织传递的热量的共同作用卜．而升高，深度方向上依次重复 上述脱水、碳化、燃烧等过程，蚀除界面和各层组织间的相界面均向深层组织推进。当 蚀除界面和相界面移动的速度基本相等时，各层组织的厚度也基本不再变化，进入准稳 定蚀除阶段。可见，在蚀除过程中，随着温度升高，组织经历不同的热物理过程，具有 多层结构的特点，应该用多层结构模型研究其作用的过程和规律。 图1为多层结构模型的示意图。以初始组织表面为原点、垂直于组织表面指向深层 组织方向为工轴建立坐标系，翔、孙粕、K。、和＆。。分别为蚀除界面、碳化层与干化层 界面、二=f二化层与脱水层界面、脱水层与凝结层界而和凝结层与L}：常组织层界面。 98l 2．1连续激光蚀除生物组织 在中红外波段，组织对激光能量的吸收系数远 大于散射系数，可以认为激光强度的衰减仅由组织 吸收造成，激光能量在各层组织中的衰减规律可以 用Beer定律描述： or／，．、 三型：一口(x)，(x) (1) 苏 其中口仅)为x处组织的吸收系数，不同层组织具有 不同数值；戤)为x处的激光强度。组织吸收激光能 量后转化热量，作为组织的体积内热源： ’ Qr=a(x)，(x) 图1：红外激光蚀除生物组织的一维多 层结构模型示意图 2．1．1初始阶段 一 在激光蚀除的初始阶段，被生物组织吸收的激光能量转换成热量在组织内传输，可 用传热方程