低温外科手术中多维相变问题的双倒易边界元模拟研究.pdfVIP

中国工.翻自口日啥 .号: 低沮外科手术中多雄相交闷扭的双侧.边界元棋拟， 邓中山，刘.*. 中国科学院低沮技术实验中心，北京2711信箱 100080 (申.联系作者:jliul.cryo.ac.cn电话:010 翻共一 本文提出了一个能用于具有复杂形状的区域的网格自动生成方法，从而将双 侧易边界元方法推广应用于涉及相变运动界面的生物传热问瓜，并成功地对低温外科手术 中正常或内含肿泪的生物组织的冷冻过程进行了数值棋拟，分析了其中的一些传热机制。 . ‘ 关悦月 低沮外科，相变，双侧易边界元方法刁生物传热 1.引， 权 . ‘ 解相变问瓜时，区城内变盈有间断及间断边界形状以及固、液界面位1未知是造 成非线性的主要因景。从物理意义上讲，这种间断是物质发生相变时所释放和吸收的 潜热造成的，也即有源 (汇)存在。目前见诸报导的关于相变问题精确解的文献不 多，且都只限于一些半无限大或无限大区域且具有简单边界条件与初始条件的理想化 情形。其中大部分还是对非生命材料相变问瓜的求解。低沮外科手术中生物组织相变 问反的边界条件必须通过运动边界上设定的条件给出，属反问橄inverseproblem)。其 求解难度更大。对于多维相变问硬，理论求解不太现实，一般采用数值方法。有限差 分法及有限元法在求解这类问反时，需要反复迭代计算相变界面的贴体坐标，计算效 率较差，边界元方法直接在边界上离徽求解，具有降维的优点，是求解相变运动边界 闷题较好的方澎”。本文将Liu和Xu发展的双倒易边界元方}pl推广到求解涉及双区 域相交的生物传热问瓜，并分析低沮外科手术中生物组织的多维冷冻过程。 2.扭攀橄祖与算法 2.1挂翻方粗 Pennes生物传热方程31的一般表达式为: P C 一 刃 汾 =V-kVT+wb(T)pbcJT一T)十4.十4, (1) 其中P，C，k分别为组织的密度、比热及热导率;几。c。分别为血液的密度及比 1 : ‘ 仁 . . 热;0)。为血液班注率;兀是动脉血沮度;4.为组织代谢率，4，为空间热源项. r U初婚妇tow 闭琳 不失一般性，以矩形区域内的二维问庄为例 (图1)，降温前整个区域的沮度呈翻 态分布。即T(x,y;to)二TO(x.Y);体核沮度为370C，即T(0.03,y;t)=370C;与 冷刀接触的体表处取等热流边界条件Q(0Y;t)=-k(aT/cc未一Cnnst s,5y5s2;体表其它部位给定为对流边界条件匆=low/.，-0C,兀=250C， 即竹(T一乙)=一(aTlan).当y=0或y二0.08或x二O,ys,或x=O,ys2; 在相 变移 动界面上 必须满足沮度连续条件和能2守恒条件 T(x,y;t}=T,(x,y;t)”T ，Y=s(x,t}以及气掀Ian)一两执Ian)=pLU， 为相变潜热.p为密度，n为界面法向，U为移动界面的法向速度。 2蛇