病变组织切除过程中的动力学仿真.docx

病变组织切除过程中的动力学仿真对于生物体病变组织的切除问题研究，传统的模拟方法只考虑了病变组织的温度情况，而事实已经证明这 是远远不够的最近，科学界提出了一种新模拟方法，通过加入反馈回路来研究不同热量和组织损坏情况 像美国食品药品管理局（FDA）这样的管理机构必须对新的器械进行审核确保仪器运行的安全性和有效性 尽管在某些情况下，也存在少量老数据与研究人员可以比较新系统的数据不符的情况但是这些仪器却能 正常的工作，有时发明者都不知道这些仪器能够在人体内的正常工作的潜在机制食品药品管理局通过进 行基础研究来填补这些在机制理解方面的不足，让代理商更好的理解潜在的机制并帮助制定指导方针和研 究测试方法在这方面的研究中，COMSOL Multiphysics软件是非常有价值的工具射频消融切除肿瘤射频消融技术是利用最新的射频消融探针，医生将把它们插入肿瘤或者病变细胞处仪器运行频率范围是 460-550kHz，通过射频探针将交变电流传送到病变组织在病变组织处离子由于电磁场极化方向的快速变 化而剧烈运动，从而使得射频能量转化为热能产生的热能使得组织细胞坏死，随着温度升高接近100摄 氏度足以杀死病变细胞。

消融技术已经证实可以有效地根除病变组织，可以治疗心律失常，减轻帕金森症症状，治疗过渡性子宫出 血，以及与化疗结合治疗恶性肝癌这些仪器对于外科医生来说都非常重要，但是它们的工作机制直到现在才被发现因为在活体实验和临床 实验中的仪器参数不同，所以利用计算机对消融进行模拟是研究这些仪器最经济和最有效的方式然而， 计算模型只能研究已知的特征参数，如电导率、热容量、血流密度，以及一些外界参数，如激励电压、时 间、最大消融温度以前研究人员采用的是利用温度数据来预测组织坏死区域大小和形状的粗糙模型这些模型首先计算电磁 场，然后求解热场来预测温度这些模型中假设组织的电学和热学参数是常数，所以坏死组织体积的模拟 结果很难接近临床观测结果在组织消融温度，蛋白质发生凝固，从而导致局部血液流动停止，并引起组 织细胞的电学参数和温度发生不可逆改变因此陈旧模型不能在外科医生在预测治疗方案所形成的组织损 坏尺寸方面提供足够的帮助模拟组织变化的更优模型FDA放射卫生与仪器中心研究人员Isaac Chang开发了一个模拟消融器件的先进模型，其中考虑了组织行 为的变化在计算每一个空间位置温度的基础上，还计算了热辐射以及累计组织损坏。

利用这些结果来改 变组织固有的参数，并反馈到下一时刻温度计算过程中因为在消融过程中组织参数不断地改变，所以在利用时间相关模型来求解消融问题为了进行模拟，Isaac Chang需要一个能够同时处理电磁和热效应的工具在那段时间里，除了 COMSOL Multiphysics他没有找 到能进行多物理场求解的工具利用电磁模块中的静电学分析模式和化工模块中的热传导模式，并采用时 间相关的迭代求解器对模型进行求解模型包含8787个节点和442045个网格每一次迭代都要计算探针 产生的电场、电流密度、热流量以及组织温度的变化同时也计算组织电导率和SAR （比吸收率）模型 计算时间为30分钟，其中15分钟是器件通电状态，15分钟为电源关闭状态，求解的时间步长为2秒 COMSOL Multiphysics通过创建几何模型和非均匀网格可以高效的对物理场进行求解，但是组织细胞在模 型中每一点都是独立的，必须对每一个点进行跟踪追踪组织累计损坏需要借助均匀直线网格点有限元 网格通常和直线网格点不能容合，所以为了避免不协调，研究人员还利用有限差分方法求解然而，相对 有限元方法，使用有限差分求解耗时长，内存需求大，当模拟消融探针中许多曲面时会形成几何假象。

Chang Modef: Coupled EM-Thermal-Tissue Damage-HysteSourceVoltage*二傩)ElectromagneticsSolverV- a(T,Q)时=0Thermal SolverpC(dT/dt)= V*VT+g(T,Q) |V《+F■ FEMLABMATLABUpdate CouplingRelationshipo(T,Q)二 g(T,Q)]Augmented square GridSquareGrid图1 Chang的消融探针模型，他在模型中利用反馈来确定随时间变化的组织参数蓝方块表示在COMSOL Multiphysics中进行的计算步骤，橙色方块表示在Matlab中进行的计算步骤 当在两款软件间交互时，需要把有限元网格转化成线格点Chang发明了一项技术来矫正这些不协调（图1）他以有限元网格开始，并在特定时间步长计算温度之后利用COMSOL Multiphysics中MESHGRID后处理函数将网格转化为线格点利用Matlab With COMSOL Multiphysics将数据传递到Matlab中，然后利用Matlab算法计算热辐射和细胞累计损坏。

算法是将线格点 上每一个位置的细胞累积损坏与阈值进行比较；如果超过阈值，血液流动停止之后他利用细胞损坏比例 计算组织的电学参数一旦Matlab算法结束，组织灌流以及校正的电导率都以数组的形式保存到交换文件 中这些扩张的数组是必须的，用来保证利用Matlab interp3命令插值时，模型边界不会出现非数（NaN）随后，COMSOL Multiphysics调用Matlab interp3命令将节点位置上电导率的值以及血液灌流插入到有限元 网格中模型给出了数学原理仿真的结果证实了他的怀疑，传统模型不能模拟射频消融中的一些关键现象结果显示电导率非常不均匀， 在消融电极表面电导率变化会达到200% （图2）在没有组织灌流时，会导致严重低估SAR和温度差 组织灌流直接影响温度结果和间接改变电导率和SAR图3是温度分布和计算损伤尺寸对比图可以很容 易看出，细胞死亡比温度减少的快Temperature Dependent Conductivity (S/m) o o oP P — P 2 Pq 5 2 5 3 5■、Radfa- Posfbon Re-af-ve fprox-ma- Edge (cm)图2利用COMSOL Multiphysics计算的电导率分布，每一个点都有不同的温度，电导率是温度的函数。

从图中可以看出随着远离探针电导率逐渐下降模拟结果与真实试验的结果非常吻合，这让Chang感到非常高兴对于每一个关键参数，他都将模拟结果 与真实实验对比来确定模拟结果的正确性通过对组织机制的模拟，他能更好的去研究新技术的影响此 外，他能够预测核心技术如何向其他领域或者应用方向转移WITH 100% NORMAL TISSUE PERFUSIONo ooTemperature (deg c)图3电压源为17.5V，消融时间为15分钟时温度(左)和组织损坏(右)的模拟结果仿真中假设100% 的正常组织灌流COMSOL在中国，中仿科技公司(CnTechCo.,Ltd)凭借个性化的解决方案、成熟的CAE产品线、专业的市场 推广能力以及强有力的技术支持服务赢得了国内众多科研院所以及企业的一致认可，目前国内几乎所有知 名大学以及中国科学院下属各研究所都已选择使用COMSOL Multiphysics作为其科研分析的CAE主要工具 随着中仿科技公司(CnTechCo.,Ltd)在全国各地的分公司、CAE技术联合中心、CAE培训中心的成立，为广 大客户提供更专业、更周到的本地化技术服务，众多企业也纷纷选用COMSOL Multiphysics作为企业的分 析工具，应用全球最先进制造技术，最终增强企业的核心竞争力，保证了企业持续发展。

COMSOL Multiphysics是一款业界领先的科学仿真软件，主要是利用偏微分方程来对系统建模和仿真它 的特别之处在于它的多物理场耦合处理能力从事专业科学研究的科研人员也可以开发具有专业用户界面 和方程设置的附加模块；现在已经有的模块有化工、声学、地球科学、高频电磁场、准静态电磁场、传热、 微机电系统、结构力学等模块软件可以在多种操作系统上使用，包括Windows、Linux、Sola芯HP-UX 等系统其他可选软件包有CAD输入模块、以及COMSOL化学反应工程实验室等更详细的介绍可参看 中仿科技网站：关于COMSOLCOMSOL公司是1986年在瑞典斯德哥尔摩成立，现在已经在多个国家(比利时、荷兰、卢森堡、丹麦、 芬兰、法国、德国、挪威、瑞士、英国、美国)成立分公司及办事处详细信息请登陆 关于中仿科技中仿科技(CnTech)是中国内地、香港、澳门领先的仿真分析软件和项目咨询解决方案的供应商CnTech 是中仿集团旗下的旗舰公司，总部设在香港，目前在全国设有四个分公司，分别是上海、北京、武汉和深 圳除了强大的销售和技术支持网络之外，我们还设有专业的售后服务团队和培训中心，为了更好的服务 广大客户，公司将陆续在全国各大主要城市设置业务分支机构。

CnTech业务包括项目咨询和专业软件的集 成和实施等，凭借多年来广大客户的支持和信任以及CnTech员工们的奉献精神和责任心，已为国内外数 百家企业、高校及科研院所提供专业软件系统及项目咨询等服务，服务领域涉及教学科研、机械工业、土 木工程、生物医学、航空航天、材料科学、化学化工、冶金科学、汽车工业、电子电器、气象环保、采矿 和石油工程等行业CnTech将始终怀着"谦虚、诚实、敬业"的态度，秉承"关注客户需求、提升客户价值"的核心理念，始终遵 循"客户满意为止"的服务准则，真诚地为客户排忧解难、出谋划策、坚持不懈，及时为国内外客户提供全 球最前沿最顶端的科技服务，力争成为仿真技术行业的典范！更详细的信息请参考http://www.CnT。