主要作者及发表信息
本研究由Wafaa Karaki、Ali Akyildiz、Suvranu De以及通讯作者Diana-Andra Borca-Tasciuc共同完成，作者均来自伦斯勒理工学院（Rensselaer Polytechnic Institute）的机械、航空航天与核工程系。该研究成果于2016年7月27日在线发表，最终刊载于2017年6月出版的《IEEE Transactions on Biomedical Engineering》第64卷第6期。
学术背景与研究目的
本研究属于生物医学工程与热物理的交叉科学领域，核心关注的是电外科手术（electrosurgery）过程中的能量耗散（energy dissipation）机制。电外科手术通过施加射频（radiofrequency，300 kHz至5 MHz）交流电来实现组织切割、消融（ablation）和干燥等效果。精确理解并预测手术过程中的热量传递及组织温度变化，对于模拟和评估组织损伤及坏死区域至关重要。目前，该领域的研究主要依赖于有限元建模（finite-element modeling）和文献中报告的组织特性来预测温度变化，但计算结果常与实验数据存在显著偏差。其中一个主要挑战在于，能量耗散和储存过程受温度、频率和含水量等多种因素影响，建模极为困难。更为关键的是，对于电外科手术所涉及的宽温度范围内组织的热物性参数，尤其是比热容（specific heat），缺乏准确的数据。现有研究大多假设热物性不随温度变化，且测量多局限于室温或体温附近，未能充分考虑高温下的相变（phase change）过程。由于软组织含水量超过70%，而水具有极高的汽化潜热（latent heat of vaporization），因此，高温下组织内水分的蒸发及其伴随的潜热损失，应是能量耗散的主导机制，但在文献中常被忽略。基于此，本研究旨在通过实验与数值模拟相结合的方法，精确测定猪离体肝脏在单极射频加热过程中，表观比热容（apparent specific heat）随温度的变化规律，以揭示汽化等过程对能量耗散的贡献。
详细研究流程
该研究构建了一套将红外热成像实验与逆有限元分析（inverse finite element analysis）相结合的完整工作流程，具体步骤如下。
首先，研究团队进行了离体组织实验。研究对象为四份采购自本地、屠宰后不超过48小时的新鲜猪肝的单叶样本。实验前，样本在运输及冷藏（4°C）过程中均采取隔热措施。实验时，将肝脏置于一个电外科手术用的分散回流垫上以确保良好电接触。核心加热装置为一个半径1.17毫米的弯曲球头电极（bent ball-tip electrode），连接至一台电外科发生器（electrosurgical unit），在400 kHz的“纯切”正弦波形下，于10、25、30和50瓦的不同功率设定下进行静态消融，每次激活持续5秒。为精确捕捉组织表面的瞬态温度场，研究团队使用了一台红外热像仪（infrared thermometry），其分辨率为480×640像素，采集频率为30 Hz，垂直置于组织上方进行记录。根据维恩定律估算，所捕获的红外辐射穿透深度仅为6至12微米，因此测量温度可视为表面温度。为防止组织表面干燥，样本在消融前被覆盖一层薄塑料膜，测试时才暴露电极周围至少30毫米半径的区域。电极由一个机械臂以特定角度夹持，以确保红外相机能无遮挡地俯拍。每次测量时，同时用示波器记录电外科发生器产生的电压波形及通过10欧姆串联电阻测量的电流，以精确计算实际输出功率，验证其与名义功率的偏差在5%以内。实验后，从红外图像中识别出电极位置，并对每一时间帧上距电极尖端相同径向距离的点的温度进行平均，从而计算出径向温度分布。
其次，研究团队构建了有限元模型（finite-element model）并实施了逆分析算法。他们使用商业软件Abaqus v6.13建立了一个30毫米×30毫米的轴对称模型，模拟组织内的耦合物理场。模型的控制方程包括描述能量守恒的Pennes类生物传热方程和描述电荷守恒的拉普拉斯方程。其中，热生成项主要由焦耳热主导，而能量存储项中的比热容被定义为“表观比热容（c_app）”，它并非真实的比热容，而是将汽化和凝固等导致的能量耗散效应包含在内。模型中使用了文献报道的温度依赖的电导率和热导率，并考虑了电极-组织界面的电接触电阻和热接触电阻。研究的核心是通过逆有限元法确定c_app随温度的变化。具体而言，他们将c_app建模为一个温度的连续分段线性函数，其形状在37°C、80°C（蛋白质凝固和汽化开始）、100°C（水的沸点）和120°C四个温度点发生变化，由四个参数（c1至c4）控制。拟合过程采用了一种顺序迭代策略以优化计算时间：先利用小于80°C的实验温度数据，通过最小化实验与模拟径向温度分布之间的差值，来拟合参数c1和c2；之后再利用80至120°C区间内的数据拟合c3和c4。此最小化过程通过Matlab中的信赖域反射算法（trust-region-reflective algorithm）实现。此外，研究还进行了灵敏度分析，通过定义一个误差度量，评估了溶液中温度分布对电极半径、热导率、电导率和表观比热容变化的敏感程度。
主要研究结果
灵敏度分析的结果表明，温度分布对热导率的变化不敏感，30%的热导率变化引起的平均温度变化小于1°C，这说明在快速加热过程中，热传导并非能量耗散的主导模式。相反，温度对电导率和比热容的变化非常敏感，30%的改变可导致超过7°C的平均温度变化，这凸显了获取准确的温度依赖性参数对于建模的重要性。
实验与反演拟合的核心结果是获得了表观比热容（c_app）随温度的显著变化规律。在37°C时，测量得到的平均c_app为3659 J/kg·K，标准差很小，且与文献中报道的3060至3870 J/kg·K的范围高度一致，验证了该方法的可靠性。然而，随着温度升高，c_app表现出戏剧性的增长。在80°C时，平均c_app增至约5330 J/kg·K，此值高于Choi等人通过密闭容器测量（排除了水分流失）所得的结果，这正是本方法计入汽化效应的必然体现。最引人注目的发现是在100°C，即水的沸点处，c_app达到了峰值，平均值高达19071 J/kg·K，是37°C时数值的五倍以上。这一巨大的峰值被明确归因于水的大量蒸发和伴随的潜热耗散。当温度继续升至120°C时，c_app又急剧下降至约7160 J/kg·K，这可以解释为组织内可供蒸发的自由水含量已显著减少。值得注意的是，这一实测峰值在形态上与Chen等人提出的数值模型一致，但在数值上要小得多，后者是基于不同实验方法（长时间加热后切片称重）获得的水分流失数据推导的。这种巨大差异突显了开发更精确描述电外科手术中相变能量损失模型的迫切需求。研究还观察到了功率依赖性，在10瓦和50瓦功率下未能获得120°C的c4参数，前者因功率密度过低无法达到沸点以上温度，后者则因功率密度过高导致电极-组织界面迅速碳化而中断了有效加热过程。
研究结论与价值
本研究成功开发并验证了一种结合红外测温和逆有限元分析的方法，能够有效测定离体组织在射频加热过程中宽温度范围内的表观比热容。研究结论明确指出，在电外科手术中，尤其是在温度接近水的沸点时，能量耗散过程由能量储存和水的汽化潜热损失所主导，而非热传导。研究提供的c_app随温度变化的定量数据，特别是其在沸点处的峰值及随后的下降趋势，为开发更精确的电外科手术计算模型提供了关键且目前缺失的参数。其科学价值在于深刻揭示了组织热损伤过程中相变的主导作用，挑战了以往忽略或简化此机制的模型；其应用价值则直接服务于新型电外科手术器械的设计、手术方案的优化以及手术安全性的提升，有助于更准确地预测组织坏死范围，减少并发症。
研究亮点
本研究的亮点主要体现在三个方面。第一，发现了一个重要且在被定量化的新现象：表观比热容在水的沸点（100°C）处出现数倍于正常值的巨大峰值，提供了汽化主导能量耗散的直接热力学证据。第二，在研究方法上具有创新性，它并非直接测量被隔绝水分的“真实”比热，而是通过瞬态加热实验和数值反演，提取了一个包含了相变潜热效应的“表观”参数，这种原位、动态的表征方式更能反映电外科手术的真实物理过程。第三，研究对象和问题的针对性极强，专注于文献中数据极为匮乏的高温区间（>80°C），直接回应了建模仿真领域对于准确高温热物性参数的迫切需求。