数值预测三维复杂纳米流体有效热导率的序参数方法.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093515 数值预测三维复杂纳米流体有效热导率 的序参数方法 1,2 1 何志祝 ，刘静* （1 中国科学院理化技术研究所, 北京 100190；2 中国科学院研究生院，北京 100039） Tel ：010，*E-mail: jliu@mail.ipc.ac.cn 摘要 本文发展一种用于预测复杂纳米流体有效热导率的新型数值计算方法。通过引入序参数使空 间不连续的热导率连续化，其物理意义在于刻画纳米颗粒和基体接触区域的热导率。本文采用对称半 隐式有限差分法求解三维纳米流体稳态温度场，继而计算出纳米流体的有效热导率。与传统的理论预 测模型方法比较，数值模拟预测方法具有不受限于纳米颗粒的体积份额配比大小，容易考察纳米颗粒 之间相互作用，且能够处理复杂形状的纳米颗粒等优点。本文工作为预测纳米流体热导率提供了一种 简单有效的新途径。 关键词：纳米流体，有效热导率，数值模拟 1．引言 实验研究表明，掺混了高热导率金属或非金属纳米颗粒的纳米流体热传导性能得到 了显著改善，例如体积含量为4%纳米氧化铜的氧化铜/ 乙二醇纳米流体的热导率提高了 22%[1] 。因此纳米流体在强化传热领域具有广阔的应用前景。近年来高热导率纳米颗粒 已被用于冷热疗医学中，特别是纳米冷冻手术概念的提出[2,3]，拓展了传统冷冻外科治 疗的极限。将一定量高热导率纳米颗粒溶液注射到目标组织部位后，可以显著提升其降 温速率及冷冻范围，从而大大扩展了常规冷冻探针的冻结能力。 目前已有大量理论模型用于预测纳米流体有效热导率，比如经典的Maxwell 模型， 考虑纳米颗粒之间相互作用的Jeffrey 模型等。近年来纳米颗粒的界面效应，聚集现象以 及微观运动等已被吸收到各种理论模型中[4]。然而大多理论预测模型只适用于纳米颗粒 体积份额较低的纳米流体，同时很难处理复杂形状的纳米颗粒，即使像椭球这样形状相 对简单的纳米颗粒都很难得出理论解。最近有学者采用数值模拟方法预测固液混合物的 有效热导率[5,6]，但其数值计算方案（有限元或有限体积）比较复杂，且很难考虑颗粒 的界面效应。本文通过引入序参数描述空间变化的热导率，使原本在纳米颗粒和基液接 触界面上不连续的热导率连续化。这样处理可使求解整个空间的温度场统一到一个方程 中，而不必对纳米颗粒和基液区域分别求解。据此可以采取有限差分法获得三维纳米流 体稳态温度场数值解，继而计算出纳米流体的有效热导率。结果表明新方法简单且能准 确预测纳米流体有效热导率。 2. 理论模型及计算方案 考察热导率为k 的纳米颗粒悬浮于热导率为k 的基液中，其温度场稳态方程可以 p m 写成 Ñ gkÑ T = 0 , （1） 资助项目：国家自然科学基资助 1 其中k = f k + 1- f k , 序参数f = 0.5(1- tanh 3d- 1r) （r 为空间任意一点到纳米 p ( ) m 颗粒与基液接触界面的距离）。f 沿界面法向呈双曲函数衰减（图 1 ），即f = 1对应于 纳米颗粒，f = 0 对应于基液，而在纳米颗粒与基液界面上存在一个长度为2d 狭窄区 域，f 从f = 1快速变化到f = 0 。此时k 是空间连续的。若不考虑界面效应，颗粒界 面温度应该满足边界条件 T + - T - = 0 ， ng k Ñ T