梯度磁场下磁流体三维结构的格子+Boltzmann方法模拟.pdfVIP

中困1-程热物理学会 传热传质学 学术会议论义 编0：063399 梯度磁场下磁流体三维结构的格子 ．Boltzmann方法模拟 宣益民 李强 张乐莹 (南京理T人学动力T程学院．南京。210094) nI：(025)843 模型．从介则角度)k研究外加梯度磁场时磁流体的三维结构。 关键词磁流体．格了-Boltzmann，三维．梯度磁场 1．引言 下，磁流体中的磁性粒子会沿磁场方向形成链；吠结构，磁流体的微观结构节现各向异性 特征，磁性粒子的分布结构对磁流体的物性参数(导热系数、枯度等)和流动传热特性 平面圆形，空间上作用的各种内力和外力也简化成了平面上的作用，造成了一定程度上 的结构，比较不同参数(时间、磁场梯度等)对梯度磁场作川。卜．磁流体分布形态的影响。 2．格子Boltzmann三维模型 基于两相的D3Q15格子一Boltzmann 模型【2i(见图1)，其粒子分布演化方程为： Z4(x+ce，At，，+At)-f,7(x,t)= a’ 一去防(x，，)一，a叫(x’f)] 式中，Z4(x，t)为在t时刻空间某点X处 速度为Ce，的第仃种粒子的密度分布函 数。，=m·■，m表示粒f质越，哆为 粒二f数密度分布函数。仃=1．…。S，这厦 幽I I)30I5述J：!}：模型 S=2。f仃为趋向局部sF衡分布单一松弛时间洲f-。 球岔项口：因家自然科学坫会重点顺口(No 1352 D3Q15模裂的离散速度为i格子速皮c：垒)： △f (o，0，o)c’ i=1 2 (2) ceJ (o，±l，o)c，(±l，o，o)c，(o，0，±1)c，i=2，…，7 5 (±1，±1，±1)c， i=8，…，1 密度平衡分布函数Z卯9(x，f)州3】： 俨～4[1+等+掣一蒡] (3) 叫·+竽+掣一虿3U2] 其中，M={ =2，…，7} 72／li =，…， l i≥8：175}I 韩M=触1／9 p(x，，)=∑f,Cx：，f)，pu(x,．，)=∑ce，Z(x，f)． (4) (5) 3．受力分析 3．1磁粒子在外梯度磁场的作用下的受力 (1)磁偶极．偶极作用力 磁偶极．偶极作用是指磁性粒子自身磁场间的相互磁作川。 磁性粒子受到周闱磁性 粒子的作．1=fj力表示如下： ¨ L，豢【m‘nl’一3(m‘ro)(m’‘ro)]ro (6) 2 其中，m及m’为磁性粒子(j)及周同磁性粒子(j)的磁矩，ro L=‘一L。 M，。为磁性粒子的饱和磁化强度．根据朗之万函