ANSYS对磁流体流动的仿真分析

1、,_ 1-"r-1=1JL*ANSYS对磁流体流动的仿真分析郭仁宁，武春晓，李锐平辽宁工程技术大学机械工程学院能源与动力工程系，辽宁阜新（123000）E-mail: wuchunxiao2005摘 要：使用传统的计算流体力学方法可以有效地对磁流体进行数值模拟，但是该方法缺 乏直观性，为提高对磁流体流动状态分析的准确性和直观性。现利用ANSYS软件对其流动进行二维绕流流场分析，重点模拟绕流状态流场速度和能量的变化，得到较理想的流场图 像模拟结果，对结论分析起到了重要作用。关键词：ANSYS ;磁流体；电解质溶液；绕流 中图分类号：TK1. 引言磁流体作为一种新型的液体功能材料，由于它

2、特殊的优点，无论是在理论研究还是在 实际应用上，都被人们广泛地关注。目前磁流体的发展主要有两个方向：一是磁性流体的基础研究，磁流体各种现象与性能的微观机理；二是磁流体的应用研究，主要有液体流动控制、 研磨、润滑、密封、热传导、振动、医疗、传感器、印刷、光学仪器等。电解质溶液能较好的模拟电磁流体在电磁场中的流动现象。采用计算流体力学方法建 立磁流体流动模型，并进行数值模拟，但是传统的数值模拟方法对流场的直观性显示比较 困难。本文运用ANSYS软件对电解质溶液流动状态进行有限元仿真模拟，对磁流体整体流动进行分析和计算，并建立直观的流场图像模型。2. 电解质溶液在电磁场下运动的流动模型的基本理论建立

3、流动模型基本假设： 电流密度J和电磁通密度 B均匀分布且相互垂直； 流体为不可压缩体，在水平放置的通道内流动，其重力影响可以不计； 流体的密度 p电导率6动力粘度卩均为常数；流动模型基本方程：1）电流密度方程J = oE =（r（Ed + u XB）2）电磁力密度方程fe = J XB3）连续性方程?知=也+竺+竺=0?x ?y ?z4）流体运动方程的矢量表达式为du2p = fe - ? P + （J? u dt流体运动方程组：-1-亠.1=1 NJ-AJL*LU国觀血ILJCtT口 - S "W>1?J = cE =(r(Ed +u XB) ?fe = J XB?Ux?Uy

4、?uz? ?u= _+Z?x?y?z?P著 fe-?P+"? dt由于方程中电磁参数和流体参数紧密联系，再加上电磁和流体两者不同的边界和初始 条件，从而导致求解比较困难。可以根据实际条件，对方程进行必要的简化，再根据边界 条件和初始条件求解，就可以得出电磁场、流场的分布状况以及流体的运动规律。3. 圆柱绕流的仿真模拟和分析根据计算流体力学理论建立的数学模型，通过ANSYS分析软件对电解质溶液的流动状态进行仿真模拟，以探讨其流体绕圆柱流动的速度，能量及流场的分布情况23。假设模拟条件：溶液密度p=1.12 X 10 kg/m 3 ;卩=0.018 XPfeD- s ; v=2 m/s壁

5、面粗糙度可 忽略。3.1单一圆柱绕流过程的仿真模拟3TTEAW“ 4 £(M1 jtimnriau aisiJTiw mmAWiff Mfl"rm iiMInr7.Q3I I.I2QB«. ID*图1圆柱绕流速度矢量图Figi Vector chart of column circumferential flow velocity图2圆柱绕流速度云图Fig2 Cloud chart of columncircumferential flow speedwmu +1tarn imtNSWANi Mtibf 訓;41 7TM-I3iirfnwdfi*ii iM-i.r

6、nLiu 血2LII-2-#-图3圆柱绕流能量分布图Fig 3 Energy distribution map of column circumferential flow图4圆柱后侧速度曲线图Fig4 Speed diagram of curves ofcolumn rear end图1-4是对单一圆柱绕流流动状态的仿真模拟56。图1、2是圆柱绕流的速度矢量图和云图。从图中可以看出，在圆柱和固定壁面之间及其附近区域，流动速度不断增大，在圆1=1 JA1-A A*Luis觀血i上xtr 口 ' ! 1-柱的偏后方达到最大值；由于圆柱对流动的阻挡，流体沿着圆柱壁面流动，并且在圆柱的 正后

7、方产生涡流现象，除边界层外，此处是速度最小的区域。在圆柱后，经过一定距离的 流动，流体的流动状态恢复到正常。图4的速度曲线图同样说明，在圆柱和固定壁面之间及其附近区域，流速增大，速度曲线达到顶点最大值，而在圆柱的正后方，速度曲线出现 最低点。图3是圆柱绕流能量分布图，根据图像颜色的变化反映出湍流能量的变化，在速 度增大的区域中，湍流能量也达到最大，在圆柱壁面处，速度变化剧烈，能量的变化也随 之剧烈。3.2并列圆柱绕流流动的仿真模拟ANIff 1Miff I 卄1PP tJAiHrihll图5并列圆柱绕流矢量图Fig 5 Vector chart of compound column circu

8、mferential flowIlMIi!Ki.M7bi9U图6并列圆柱绕流速度云图Fig5 Cloud chart of compound columncircumferential flow speed图7并列圆柱绕流能量分布图Fig 7 Energy distribution map of compoundcolumn circumferential flow图8第一个圆柱后侧速度曲线图Fig 8 Diagram of curves of first columnrear end speed图5-8是对并列圆柱绕流流动状态的仿真模拟56。从图5, 6显示出，同样在两圆柱与固定壁面之间及其

9、附近区域，流动速度不断增大，并且在这区域之间有一小段速度减小的 区域。但是在图 6中可以反映出，第一个圆柱的速度增大区域比后者的略大些。同样在两 圆柱正后方都出现了涡流现象，但是在第一个圆柱后的涡流区域明显大于后者，紊乱的程 度比较强烈。在其后流动恢复正常的距离要长。在能量分布图中，同样反映出在第一个圆 柱的能量变化强度比后者强烈些，而速度曲线图去单一圆柱时相似。4. 结论（1 ）在圆柱和固定壁面之间及其附近区域，出现流速增大区域，并出现最大值；由于圆柱-3-1=9 HU-AA*UJ 国 SS It 3C对流动的阻挡，在圆柱的正后方产生涡流现象。(2) 在圆柱壁面及附近的固定壁面处出现紊流能量

10、强烈变化区域。(3) 由于并列圆柱对流动阻碍作用的增强，流动出现更强烈的紊乱，并且在第一个圆柱后 的涡流区域明显大于后者，紊乱的程度比较强烈。本次利用ANSYS软件对电解质溶液流动进行的二维绕流流场分析，重点模拟了速度和 能量的变化。为改进实验和提高结论分析的准确性提供了有力的帮助。参考文献1 刘顺隆.计算流体力学M，哈尔滨：哈尔滨工程大学出版社，1998.2 刘涛，杨凤鹏.精通ANSYSM，北京：清华大学出版社，2002.3 龚曙光.ANSYS工程应用实例解析M，北京：机械工业出版社，2003.4 Myung-man Kim, Andrew L Zydney. Effect of Elect

11、co-static, Hydrodynamic, and Brownian Forces on Particle Trajectories and Sieving in Normal Flow FiltrationJ, Journal of Colloid and Interface Science, 2004, 269(2): 425-4315 D E Armin, J A Ritter, H J Ploehn. Magnetic Hetero-Flocculation of Paramagnetic Colloidal ParticlesJ, Journal of Colloid and

12、Interface Science, 2000, 225(1): 39-46ANSYS Inc.ANSYS/FLOTRAN计算流体力学分析指南Z，美国ANSYS公司北京办事处ANSYS to simulation analysis of magnetic fluid flowingGuo Renning, Wu Chun xiao, Li Ruip ingEn ergy and power engin eeri ng departme nt of mecha ni cal engin eeri ng college, Lia oning Tech nicalUn iversity, Fuxi

13、n, Liao ning (123000)AbstractThe traditi on method of computati on hydromecha nics can effectively carry on the value simulati on to the magnetic fluid, but this method lacks intuitivism. To enhance the accuracy and intuitive of the magnetic fluid flow analysis. Now, using the ANSYS software to analysis the two-dimensional circumferential flow field. The key of simulation is the speed and energy change of circumferential flow field, to get the ideal image result of