基于模拟方法生成血管树的人舌三维温度场数值模拟研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号： 083323 基于模拟方法生成血管树的人舌三维温度场数值模拟 1 1 1 1 2 王雅博 诸凯 张艳 解海卫 姜智浩 1．天津商业大学天津市制冷技术重点实验室， 天津 300134 ；2. 天津中医药大学，天津300072 Email:zhukai210@tju.edu.cn 摘要：采用分步模拟方法生成血管树并以此作为物理模型，进行舌体三维温度场的实验研究与数值计 算。计算分为舌体组织和血液温度场数值模拟两部分，相应数学模型也包括舌体组织控制方程、血液 控制方程及两者之间的耦合方程。考虑了血液与组织间的耦合换热，舌津液蒸发的影响，使用数值模 拟得到的舌面对流换热系数，使计算结果与红外热像验证极为近似。以模拟方法生成的血管树替代生 物体的血管铸型，在某种意义上说将生物体三维传热计算方法的研究推向了一个新的层面。 关键词：舌 血管树 血管铸型 三维温度场 生物传热 0 前言 以人舌血管铸型作为内热源，进行舌体三维温度场的数值计算，体现出了依据生物 体的真实形状和真实受热状态而进行的生物传热研究[1-2] 。但所带来的问题是：人舌血 管铸型的获取绝非易事。在此背景下本课题提出了用分步模拟的方法获取血管树[3]，即 先根据统计规律建立生物体动脉血管树的根血管和一级分支，然后根据生理组织结构特 征和优化原则生成血管树的其它分支。此方法的提出对于生物体三维温度场的计算研究 具有重要意义。本文将在用分步模拟方法生成人舌动脉血管树的基础上，利用数值模拟 [3] [4] 得到舌面对流换热系数 、并考虑到唾液蒸发的影响 ，对人舌三维温度场进行了数值 计算。 1物理模型 舌体模型包括舌体组织和血管树两部分。舌体可以分成左、右各半，其左右两侧具 有基本独立的血管树。本文只建立半舌的模型及其内部的血管树。按照文献[3]介绍的分 步模拟血管树方法，建立三维人舌动脉血管树模型。主要 参数为：舌深动脉半径为 0.9mm；一级分支个数在 10-15 之间，具体个数由随机函数确定。图 1 所示的舌体物理模 型是在舌体模型中挖去血管树部分，然后用绘制的所有血 管的轴线体现三维血管树。 图1人半舌及内部血管的物理模型 2数学模型的建立 2.1 舌体组织控制方程及血液控制方程 舌体组织的控制方程如下式： div λgradT +ρ c W T −T +Q 0 （1） ( t t ) b b b ( b t ) met 式中：ρt 组织密度；Qmet 生物代谢热；T 组织温度；c 组织定压比热；λ组织导热系数。 t t t 资助项目：国家自然科学基金(No. 根据血管树模拟过程可知，生成的血管树忽略血液的温度在血管横截面上的变化，即只 考虑血液温度沿管长的变化。用一维的血液温度场替代了三维的血液温度场。管树中每 个血管段的流量和半径是确定的，相应每个血管段的速度可以通过简单计算得到。因此 血液的控制方程为血液的一维传热方程： ∂ ⎛ ∂T ⎞ ∂T λ b −ρ c u b 0 （2 ） ⎜ b ⎟ b b ∂x ⎝ ∂x ⎠ ∂x 2.2 血液与舌体组织之间的耦合传热方程 血液与舌体组织之间的耦合传热方程是指血液与周围组织之间的换热关系，即血液 在血管壁面流