冷冻外科中一维相变问题的数值解.pdfVIP

第四届全国低温工程 低 温 工 程 1999年第4期 学术会议论文集 CRYOGENICS 总第110期 冷冻外科中一维相变问题的数值解 张洁 华泽钊 (上海理工大学制玲与低温研究所上海200093) 摘要 用变时间步长法求解了冷冻外科手术时模拟生物组织中的一维相变问 题。计算方法经与文献比较，简便可行。运用本文的计鼻方法计算并讨论了不同边 界条件下，相变区的温度场分布和无量纲相界面降温速率随Ste数的变化。计算分析 表明，相变区的温度基本呈线性分布，无量纲相界面降温速率随Ste数增大而降低。 所得结论可对冷刀设计和冷冻外科手术提供帮助。 主题词冷冻外科相变问题变时间步长法 I前言 冷冻外科是60年代发展起来的外科手段，将被低温介质冷却的手术器械(冷刀)插入病 变组织，通过特定的降一复温程序，达到冻伤、杀伤病变细胞和组织的目的“]。低温对生物细 胞可以产生两种截然不同的作用，一方面，一定的低温可以充分降低细胞的新陈代谢能力，达 到长期保存的目的；另一方面，在降低到某一温度和恢复到常温的过程中，由于细胞内外水 分的结晶和晶体的融化、再结晶，细胞膜中蛋白分子的变化；电解质的浓缩和热冲击等原因口]， 会造成组织和细胞的损伤或死亡。冷冻外科正是利用后一机理通过控制冻结区的温度场分布 探求最大限度杀伤病变组织的方法。因此，组织冻结过程的预测和分析对冷冻外科手术方案 的制定和准确实施起着十分重要的作用。 在于移动的相变界面。相变界面上，由于潜热的吸收或释放，存在热流的阶跃，因此这一问 题在数学上是非线性的。实践中的大量相变问题，只有少数简单情形可以获得解析解u5]，多 数问题则需借助于数值解法口]。数值求解Stefan问题的方法可分为固定步长法、焓法、坐标 变换法、变步长法等o]。固定步长法即空间和时间步长保持不变的求解方法，由于相界面不可 能恰好与网格节点重合，因而需要不断的插值以确定各个时刻的相界面位置；焓法对整个求 解区域建立统一的守恒方程，求出热焓分布后，再根据焓值来确定相界面位置，求得的相界 面并不是某一明确位置而是某一焓值对应的相变区域，此法适合多维问题的求解；坐标变换 法通过坐标变换使运动的相界面变成不动的边界，但这种方法使控制方程变得复杂；变步长 法包括变时间步长法和变空间步长法。变时间步长法用于一维相变闻题的求解，概念清晰简 张洁，女，z8岁，博士 240 低温工程 1999年第4期 便可行，虽在冶金、铸造等工程领域得到应用【6’7]，但未见用于冷冻外科的报道。本文尝试采 用变时间步长法分析冷冻外科中典型边界条件下的Stefan问题。 影响冷冻外科手术疗效的因素主要有冷刀壁面达到的最低温度、冷刀壁面温度变化规律、 手术持续时间和冻一融频率等啪。文献[5]曾报道冷刀壁面温度恒速下降时组织中相变过程 的摄动近似解和实验值，其他壁面降温规律对组织冷冻过程的影响则未见报道。本文用变时 间步长法数值求解上述问题并和文献E53的结果进行了比较，得到了较为一致的结果。在此 基础上，应用变时间步长法进一步计算和分析r其他壁面降温规律(边界条件)下，组织中 温度分布和相界面降温速率随Ste数的变化规律。 2数学模型 由于生物组织中水的含量为60％～90％，作为初步模拟，假设组织具有纯水的物理特性； 各物性参数不随温度变化；组织内物性均匀且各向同性；不考虑组织内的自然对流；组织的 初始温度为相变温度。 直角坐标系中的一维相变导热问题可描述为 肛警一kk鍪肛i一面 ㈩(1) 0≤z≤s(t)，t0 初始条件为 丁(z，0)一Tm (2) j-≥0 边界条件为