1 问题背景

    肾动脉消融（Renal denervation，简称RDN）是一种治疗高血压的办法，其基本原理为利用插入肾动脉的电极消融导管进行射频消融，使肾动脉血管壁附近的交感神经因高温而损伤失活，减少神经系统过度活跃的信号刺激。

RDN原理示意图（图源：www.cirse.org）

    RDN治疗过程对治疗区域的温度控制有很高要求，温度太高会导致肾动脉的血管受损；温度太低治疗无效。因此利用ANSYS Fluent仿真，了解血管和周围软组织的温度分布情况，确保温度处于合适范围内。

2 仿真思路概述

    计算域如图所示，包括血液、血管壁、周围组织。消融设备仅保留和血液接触的外表面，不考虑内部的结构。

    此问题同时涉及三个物理场：电场、流场、温度场。

2.1 电场处理

    电场激励源为射频信号，其频率为MHz量级，变化周期相对于秒级的加热过程太短，可简化为设备外表面的固定电压，其电压值为射频信号的有效值。计算域的外边界为0电压。

    简化后，整个电场为静电场，满足电势的扩散方程，且为定常。电场的传播速度为光速，可认为在通电的一瞬间立刻完成，不考虑电场的传播分布。

2.2 流场和温度场处理

    根据解剖学，人体肾动脉的入口直径范围为5 - 7mm，分叉处直径范围3 - 3.5mm，分叉角度约为50°。其入口速度大约为50cm/s左右。可推算，血流雷诺数在1000左右，为层流状态。

    组织温度分布的影响因素包括：电场的焦耳热、血流和周围组织中的热对流和热传导。

3 网格划分

    几何模型采用ANSYS Discovery软件进行建模处理。

计算域总览

    采用Fluent Meshing划分多面体网格，其尺寸范围0.04 - 1 mm，总数约130万

网格总览

消融设备区域网格细节

4 模型设置

    模型分为两步处理：

1 稳态仿真，计算电场和血流分布

2 瞬态仿真，计算温升过程 

    稳态仿真中，电场模型必须开启焦耳热选项。计算域温度统一为37℃。

   瞬态仿真中，将冻结第一步计算的电场和流场，成为一个单纯的热传导过程，以节约计算量。

稳态仿真中关闭能量方程

 瞬态仿真中仅开启能量方程

5 仿真结果

 血液速度分布

 计算域电势分布

计算域焦耳热功率分布

 血管最高温度变化

 周围组织最高温度变化

Media1

计算域温度分布变化动画