自适应边界元-分布残差法研究含局部相变的瞬态低温生物传热

《Case Studies in Thermal Engineering》：Transient cryogenic bioheat transfer with localized phase change: An adaptive boundary element–distributed residual case study

　　

局部低温冷却会产生陡峭的温度梯度、移动相变前沿、灌注与代谢源项以及潜热缓冲效应，这些对传统体积分求解器构成挑战。本案例研究采用自适应边界元-分布残差（BEM-DR）格式，考察二维圆形组织域中径向对称加载下瞬态冻结前沿的演化。边界积分处理平面导热，残差驱动的背景基函数与局部富集基函数重构分布源与集中探针冷却。表观热容（apparent heat capacity）模型表示有限糊状区间（mushy interval），后向欧拉（backward-Euler）积分推进瞬态解。光滑制造解（manufactured solution）验证给出底层平面BEM离散近二阶空间收敛；该测试不声称确立穿越相变的非线性收敛率。时间细化给出一阶行为。基函数细化使120 s冰球半径仅变化0.1%，最大全局能量残差为0.84%。在匹配几何、物性、边界条件与冷却历史下，与有限差分焓法（enthalpy solver）跨代码比较，冰球半径差异为3.2-4.2%，探针中心温度差异4.8%，平均边界热流差异4.6%。参数化算例显示导热系数、灌注、代谢与探针功率如何影响冻结前沿传播、近探针梯度与热流恢复。结果展示了用于重复低温热工程分析的网格缩减框架。本研究确立数值验证与跨方法一致性；临床预测前仍需实验验证。

低温冷冻消融（cryoablation）中组织温度场的准确预测是生物医学传热的核心问题，其疗效评估、周围组织保护及热损伤判读均依赖于空间与时间上非均匀的温度历史。经典Pennes模型建立了耦合导热、血液灌注与代谢产热的工程框架，此后边界元格式被扩展至广义生物传热问题。生物组织的热响应因温度相关的热物性与力学性质而进一步复杂化，导热系数、热容、密度在生理、热疗与冷冻状态下可发生显著变化。冷冻消融是尤为苛刻的热工程工况，能量被快速移出，施用器周围形成冻结区，细胞损伤源于热、渗透、血管与结构机制的相互作用。从数学角度看，低温生物传热是非线性的，因为相态与有效热存储随温度变化，同时灌注与内源项共同作用。焓基固定网格格式因无需显式追踪尖锐移动边界而成为处理冻结与融化的基础方法。边界积分法在传导主导问题中具有吸引力，Green函数表示将主要几何离散缩减至边界，同时保持精确的边界热流恢复。然而，双互易（dual reciprocity）与分布源格式的精度关键取决于域项的近似方式，固定全局基函数可能对集中低温汇欠分辨。

针对上述问题，研究人员开展了二维平面圆形组织域中径向对称加载的瞬态冷冻前沿演化案例研究，采用自适应边界元-分布残差（BEM-DR）格式。研究目标包括：(i) 推导可复现的BEM-DR格式用于准稳态与瞬态低温传热；(ii) 构建针对集中冷却与焓相关源项的自适应表示；(iii) 验证底层空间离散并评估时间细化、基函数独立性、条件数与能量闭合；(iv) 在匹配假设下与有限差分焓格式进行跨代码数值比较；(v) 量化冻结前沿传播、近探针温度、边界热流、参数敏感性及双探针相互作用。研究得出：制造解测试显示近二阶空间收敛，后向欧拉时间细化产生一阶时间收敛；基函数细化使120 s冰球半径变化降至0.1%，最大全局能量残差低于1%；与有限差分焓基准的差异为冰球半径3.2-4.2%、探针中心温度4.8%、平均边界热流4.6%。参数化算例阐明了导热系数、灌注、代谢产热与探针功率对冻结前沿运动、潜热缓冲、近探针温度梯度与边界热流恢复的影响。该方法最适用于需要在固定二维几何下重复模拟的传导主导低温研究。本研究意义在于建立了数值一致性与跨方法一致性证据，提出了在可重用边界积分传导框架内对源表示进行残差控制局部富集的方法学贡献；但光滑制造解验证的是底层平面离散而非非线性相变收敛，有限差分研究提供跨代码一致性但缺乏独立收敛档案，实验体模与离体测量、温度相关组织物性及三维扩展仍是患者特异性预测前的必要步骤。

研究人员开展研究用到的主要关键技术方法包括：自适应边界元-分布残差（BEM-DR）格式，以二维Laplace与修正Helmholtz基本解构建边界积分算子；表观热容模型表示有限糊状区间并嵌入潜热；后向欧拉隐式时间积分推进瞬态解；固定点（Picard）线性化处理非线性相变更新；Tikhonov正则化最小二乘求解分布残差基函数系数；残差指示子驱动的两级基函数（高斯径向基与紧支撑Wendland基）局部富集；制造解与时间/基函数细化进行代码与解验证；与独立实现的有限差分焓法进行跨代码比较。本研究为数值案例研究，未使用实验体模、离体或患者特异性数据集。

研究结果方面，代码验证通过制造解给出平面BEM-DR离散的空间误差随网格细化近二阶下降，L2与L∞误差曲线在对数坐标下平行于N?2参考线，但该测试不构成非线性相变前沿的收敛证明。解验证通过时间与基函数细化表明：后向欧拉时间细化给出一阶时间收敛，糊状区附近无振荡；分布残差基函数从16增至256时，120 s冰球半径从12.82 mm增至13.64 mm，最终细化变化仅0.1%，边界热流相应稳定，支持所报告指标的基函数独立性。跨方法数值基准在匹配几何、物性、边界条件与冷却历史下，BEM-DR与有限差分焓法在60 s、120 s、600 s冰球半径差异分别为3.2%、4.2%、3.5%，120 s中心温度差异4.8%，平均边界热流差异4.6%；由于次级FD代码的独立网格与时间收敛档案未提供，该差异仅解释为两种独立离散格式在匹配假设下的一致性。守恒、条件数与数值质量方面，最大报告步的全局能量残差为0.84%，边界热流闭合不平衡为1.1%，糊状区无可见过冲，正则化后条件数为3.2×10?，均满足验收目标。物理验证状态方面，本研究未使用实验体模、离体或患者特异性数据，结果仅确立实现正确性、数值稳定性及与独立数值格式在匹配假设下的一致性，不确立物理或临床有效性。热物理结果与案例解释方面，基线低温案例从均匀生理温度开始，在1 mm探针核心内施加强集中冷却负荷，响应呈双相：早期快速冷却阶段由局部汇主导，随后冻结前沿扩展因潜热与传导-灌注热流入而减慢；给出了径向温度剖面、?1 °C等温线（冰球半径）演化、探针中心冷却历史及基线单探针热响应指标表。敏感性分析表明：纳米颗粒使近探针导热系数增加30%时120 s冰球半径增加1.5 mm（+11.0%）；代谢产热加倍使半径减少0.8 mm（?5.9%）；灌注汇激活使半径减少1.5 mm（?11.0%）；探针功率降低20%使半径减少2.2 mm（?16.2%）。双探针演示表明8 mm间距在85 s形成连续冻结桥，12 mm间距在176 s形成中等热重叠，18 mm间距在600 s内冰球保持分离。计算效率方面，二维圆形基准中BEM-DR使用128个边界单元与128个DR基函数，有限差分比较使用160,000个网格单元，BEM-DR每案例运行时间18 s、内存0.18 GB，FD焓法为420 s、1.60 GB，瞬态有限元为760 s、3.10 GB，BEM-DR在固定几何重复参数扫描中具有计算优势。

讨论部分指出，作为案例研究，规范二维平面圆形配置在引入几何复杂性之前隔离了局部源重构与相变缓冲的热后果。主要工程发现是残差驱动富集将近似努力集中于探针附近，而边界算子在载荷与材料参数变化间保持可重用。本数值证据支持三项发现：BEM-DR格式保持边界元的主要优势，即无需密集域网格即可实现精确边界热流恢复与场重构；分布残差表示容纳低温生物传热所需的体积项，包括灌注、代谢、潜热缓冲与强局部汇；该方法产生热可解释量如冻结前沿半径、中心温度、边界热流与温度梯度，可直接与焓基数值解比较并未来与受控实验比较。纳米颗粒增强案例说明局部导热系数修改可在不改变边界离散的情况下表示。该方法并非所有体积求解器的替代品，解剖学强非均匀组织、轴对称或全三维几何、物性温度依赖性、界面运动、各向异性灌注网络与力学损伤模型可能需要多区域格式、域分解技术、BEM-FEM组合方法或已知的焓/水平集/相场算法。局限性方面，规范几何尚不代表应用特异性解剖，表观焓模型平滑冻结区间，弱非均匀性假设可能欠分辨层状组织，合成计算基准不确立物理有效性，未包含逆控制器，二维平面圆形基准的x-y平面径向对称并非真正三维轴对称模型，常数导热系数未解析冻结/未冻结的k(T)变化，有限7 °C糊状区间在更窄区间下可能使表观热容刚化并允许时间步跳过潜热，跨代码FD比较缺乏次级FD网格/时间收敛历史，求解器溯源元数据（迭代计数与β日志）未保留。结论部分翻译如下：本研究提出了二维平面圆形域中带局部相变的瞬态低温传热的自适应边界元-分布残差案例研究。Green函数边界算子与灌注、代谢、焓相关源及集中探针冷却的残差驱动重构相耦合。所声称的数值新颖性是在可重用边界积分传导框架内对源表示进行残差控制局部富集，而非表观热容概念或径向基近似本身。制造解测试展示了近二阶空间收敛，后向欧拉时间细化产生一阶时间收敛。基函数细化使120 s冰球半径变化降至0.1%，最大全局能量残差保持在1%以下。在匹配假设下，与有限差分焓基准的差异为冰球半径3.2–4.2%、探针中心温度4.8%、平均边界热流4.6%。这些结果确立数值一致性而非物理验证。参数化算例阐明了导热系数、灌注、代谢产热与探针功率对冻结前沿运动、潜热缓冲、近探针温度梯度与边界热流恢复的影响。因此该方法最适用于需要在固定二维几何下重复模拟的传导主导低温研究。光滑制造解验证底层平面离散而非非线性相变收敛；有限差分研究提供跨代码一致性但缺乏独立收敛档案。实验体模与离体测量、温度相关组织物性及三维扩展仍是患者特异性预测前的必要步骤。