磁性流体导热系数的数值模拟研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编譬：063320 磁性流体导热系数的数值模拟 李强黄璐宣益民 (南京理T大学动力工程学院．南京。210094) 摘要基于磁流体的三三维微观结构。计算模拟丫自．、尤外加磁场作用条件下磁流体的导热系数， 分析了粒子聚集程度、粒子体积份额、外加磁场作用势对磁流体导热系数的影晌。结果表明：不 论有无外加磁场作用，磁流体的导热系数均大子基液且随磁性粒了体积份额的增加而增大。在无 外加磁场作用时，磁流体的导热系数随着粒了聚集程度的增人而减小。=fI=肯外加磁场作用时，外 加磁场作用势强度越大，磁流体的导热系数越人；沿外加磁场方向的磁流体导热系数将大于je它 方向的磁流体的导热系数，磁流体导热系数呈现各向异性特征。 关键词、磁流体导热系数微观结构各向异性 1．引言 在外加磁场的作用下，磁流体中的磁性粒子会沿磁场方向形成链状结构，磁性粒子 的链状分布结构对磁流体的导热系数起着至关重要的影响作用，磁流体的物性参数会随 着外加磁场的强度、方向以及磁场性质的变化而改变，呈现各向异性特征。实验研究表 明，磁流体的导热系数人予基液导热系数【I’2】，当外加与热流传播方向平行的磁场，磁 流体的导热系数大于相同磁场强度的与热流传播方向垂直的磁场作用．卜．的导热系数121。 本文基于分子动力学(MD)方法所模拟的磁流体在有、无外加磁场作用下的三维 微观结构f3】，计算模拟有、无外加磁场作用条件下磁流体的导热系数，分析粒子聚集程 度、粒子体积份额、外加磁场作用势对磁流体导热系数的影响，探讨磁流体导热系数与 磁流体微观结构、外加磁场间的本构关系和磁流体导热系数的控制方法。 2．数值模拟 边长为￡的立方体，其中含有Ⅳ个直径为以的球 热漉方向 形磁性粒子，悬浮在粘度为刁的基液中，粒子体积份 额为≯。如图1所示，假定这个三维区域左边界和以 Z 边界的温度已知，上下前后边界为绝热肇面，无内热 Y 源。考虑整个系统的能量守恒。导热方程组为【4】： X o (1)‘l’ 幽I磁流体热传导，J÷意幽 宴+宴+OZT；：o 丽+萨+万5 堆会项目：图家自然科学荣合重点项日(No． 1242 一七鼍L。=0；“豢l二。=。：．七鼍l：。。=。；一后鼍l：；。’_，o (2) rb=L孔：。=瓦 (3) 针对文献[3]所模拟的磁流体三维微观结构，进行网格划分，根据能量平衡法得 到各点温度的离散方程。控制容积分界面上的当餐导热系数用调和平均公式计算14l，绝 热壁面上的边界：肖点采用附加源法【41。在模拟区域达到稳态后，任取沿热流传播方向的 某一截面，将这一截面上各网格的导热量由Fourier定律算出后叠加，得出总的导热量， 从而计算出磁流体的导热系数。 3．结果分析 3．1 粒子聚集对磁流体导热系数的影响 图2给出了l％体积份额的磁流体在有无外加磁场条件下、不同时间步长时的导热系 数。从计算结果可以看出，无论有无外加磁场’，磁流体的导热系数始终人于基液的导热 系数。在无外加磁场条什下，随着时间步数的增加，粒子聚集程度增人，导热系数比* 明显的下降趋势，不同方向的磁流体导热系数无明显差距，磁流体的