Cu-Co-Fe合金雾化合金液滴凝固过程研究*
赵雷1,2, 江鸿翔2, AHMAD Tauseef2, 赵九洲2
1 辽宁石油化工大学机械工程学院, 抚顺 113001
2 中国科学院金属研究所, 沈阳 110016

赵雷, 男, 1981年生, 讲师, 博士

摘要

对Cu-10%Co-10%Fe (质量分数)亚稳液态组元不混溶合金开展了气体雾化快速凝固实验, 制备了富Fe-Co相球形粒子均匀分布于基体Cu的复合粉末, 建立了Cu-Co-Fe合金雾化液滴冷却过程中的温度场、浓度场和液-液相变动力学控制方程, 研发了耦合合金热力学和相变动力学的模拟方法, 模拟分析了Cu-10%Co-10%Fe合金雾化液滴的凝固组织形成过程. 实验和模拟结果表明, 在气体雾化快速凝固条件下, 液-液相变过程中富Fe-Co相液滴Marangoni迁移和Ostwald熟化的影响很弱, 粉末中心绝大部分区域内富Fe-Co相粒子的空间分布均匀. 对于直径小于220 mm的Cu-10%Co-10%Fe合金粉末, 富Fe-Co相粒子的平均半径Ra和数量密度N与雾化粉末直径d之间符合指数关系.

关键词: Cu-Co-Fe合金; 液相分离; 快速凝固; 建模及模拟
中图分类号:TG111.4, TF83
STUDY OF SOLIDIFICATION FOR GAS-ATOMIZED DROPLET OF Cu-Co-Fe ALLOY
Lei ZHAO1,2, Hongxiang JIANG2, Tauseef AHMAD2, Jiuzhou ZHAO2
1 School of Mechanical Engineering, Liaoning Shihua University, Fushun 113001
2 Institute of Metal Research, Chinese Academy of Sciences, Shenyang 110016
Correspondent: ZHAO Jiuzhou, professor, Tel: (024)23971918, E-mail:jzzhao@imr.ac.cn
Abstract
Keyword: Cu-Co-Fe alloy; liquid phase decomposition; rapid solidification; modeling and simulation

Cu-Co-Fe合金具有高强度、高导电和巨磁电阻效应等特性, 在电接触材料、集成电路引线框架材料、电线电缆材料和铁磁性材料等领域具有重要的潜在应用前景[1,2]. 然而, 该合金存在亚稳液态组元不混溶温度区, 合金熔体冷却时发生液-液相变, 极易形成偏析严重的凝固组织[3,4]. 近年来, Cu-Co-Fe合金体系常被作为一种研究具有亚稳液态组元不混溶温度区间合金凝固行为的模型合金[5,6,7,8,9,10,11,12,13,14]. 以往研究工作主要集中在热力学上. Kim和Abbaschian[5]测定了Cu-Co-Fe合金的液相线温度和亚稳液相分离温度, 发现当Co和Fe的含量在一定范围内时, 合金具有亚稳液态组元不混溶温度区. Cao和G?rler[6]对Cu-Co-Fe合金进行了差热分析(DTA)实验, 测定了Cu组元的原子分数在10%~84%之间且Co和Fe原子比为1∶3, 1∶1和3∶1 时的液相线温度和亚稳液相分离温度. Curiotto等[7]采用差示扫描量热法(DSC)测定了Cu-Co-Fe合金液相线温度和亚稳液相分离温度. Munitz等[8]对Cu-Co-Fe合金体系的热力学参数进行优化, 并计算了亚稳液相分离温度. Turchanin等[9]利用量热法研究了Cu-Co-Fe合金的液相混合焓, 并利用计算相图方法评价了该体系的热力学参数. Bamberger等[10]、Wang等[11]和Palumbo等[12]也分别对该体系的热力学参数进行了优化和评价, 并计算了相图. 在动力学上, Munitz等[8]对Cu-Co-Fe合金开展了深过冷凝固实验, 研究了过冷度和冷却速率对合金凝固组织的影响, 结果表明, 当过冷度较高时合金发生了液-液相变, 富Fe-Co相形貌随过冷度增大的变化顺序为: 枝晶组织→枝晶+球状的混合组织→球状组织. Munitz和Abbaschian[13]利用电子束表面熔化方法研究了高冷却速率下Cu-Co-Fe合金的液-液相分离行为. 代富平等[14]利用落管研究了Cu-Co-Fe合金的液-液相分离过程. 研究表明, 快速凝固技术可以使Cu-Co-Fe合金快速通过液态组元不混溶区, 能有效地抑制液相分离和偏析组织的形成.

通常, 液-液相变过程中的组织演变是弥散相液滴形核、长大、粗化及空间迁移等共同作用的结果, 影响因素多, 且各影响因素间又相互作用, 十分复杂, 仅通过实验手段很难深入探索其凝固过程, 计算机模拟可以有效地揭示难混溶合金的凝固机理[15,16]. 本工作以Cu-10%Co-10%Fe合金为对象, 用气体雾化快速凝固技术和计算机模拟相结合的方法, 研究快速凝固条件下Cu-10%Co-10%Fe合金的组织形成过程.

1 实验方法

以质量分数大于或等于99.9%的Cu, Co和Fe为原料配制合金, 成分为Cu-10%Co-10%Fe (质量分数). 实验过程如下: 首先将原料在KGPS100/2.5真空中频感应炉中加热到1420 ℃, 保温15~20 min, 确保组元完全互溶, 形成均一熔体; 然后用N2气体将由安放于坩埚底部导流管流出的合金熔体雾化成液滴, 实现快速冷却凝固, 雾化气体压力为3 MPa; 最后对雾化粉末进行筛分、镶嵌和抛光处理, 利用Axiovert 200MAT光学显微镜(OM)和S-3400N扫描电子显微镜(SEM)观察凝固组织, 并进行能谱(EDS)分析. 利用SISC IAS V8.0定量金相分析软件对微观组织中第二相粒子尺寸和分布进行统计分析.

2 结果与分析
2.1 显微组织

图1为不同直径Cu-10%Co-10%Fe合金粉末的显微组织和数量密度分布. 由Cu-Co-Fe体系的垂直截面相图[5]和EDS分析(图2)可知, 深色第二相为富Fe-Co相, 浅色基体为富Cu相. 当雾化粉末直径小于220 μm时, 粉末的显微组织由富Fe-Co相球形粒子和富Cu基体相组成, 随着粉末尺寸的增加, 富Fe-Co相球形粒子尺寸逐渐增大. 由文献[5]可知, 当Cu-10%Co-10%Fe合金熔体的过冷度大于24 K时, 雾化液滴发生液-液相变, 富Fe-Co相液滴在富Cu液相中形核、长大. 由于Cu-Co-Fe合金体系中组元互溶温度线(binodal line)位于液相线(liquidus)下方, 因此, 在亚稳液-液相变过程中, 富Fe-Co相小液滴形成后即处于过冷状态, 倾向于凝固. 如果雾化液滴的冷却速率足够快, 以致在雾化液滴冷却通过组元不混溶区时富Fe-Co相小液滴尚未凝固, 则这些富Fe-Co相小液滴将在包晶反应温度以下凝固, 形成富Fe-Co相球形粒子, 雾化液滴凝固后将呈现图1所示的凝固组织. 此外, 富Fe-Co相球形粒子在粉末中基本呈现均匀分布, 只是在除雾化粉末表面很小的区域内富Fe-Co相球形粒子数量密度减少(图1b和d).

Fig.1 SEM images of Cu-10%Co-10%Fe powders with diameters of 88~100 mm (a), 125~150 μm (b), 200~220 mm (c) and distribution of number density (N) of Fe-Co-rich phase particles along radial direction (r) of the Cu-10%Co-10%Fe powder with diameter of 125~150 mm (d)图1 不同直径Cu-10%Co-10%Fe合金粉末的SEM像和数量密度分布

Fig.2 EDS analysis of light (a) and dark (b) phases in Cu-Co-Fe alloy图2 Cu-Co-Fe合金的EDS分析

2.2 雾化液滴凝固过程模型

2.2.1 模型建立 在凝固过程中, 雾化液滴表面热量通过与周围气体环境进行对流换热, 雾化液滴内部热量通过扩散和弥散相液滴运动进行传输, 雾化液滴温度场满足:

(1)??(k?T)+QS/L-Qd-4π3??0uMF(ρβCPβ-ρmCPm)TR3dR=ρCP?T?t

式中, T为热力学温度; ρ, CPk分别为合金密度、比热容和热导率; ρmρβ分别为基体和弥散相的密度; CPmCPβ分别是基体和弥散相的比热容; QS/L为固/液界面上凝固潜热释放速率; Qd= h(Tds-Tg)为液滴表面对流换热速率[17], 其中, Tds为雾化液滴的表面温度, Tg为雾化气体的温度[18], h=kg(2+0.6Re?Pr)/d是雾化液滴表面与雾化气体间的界面换热系数[17], d为雾化液滴的直径, kg为雾化气体的热导率, Pr为气体的Prandtl数, Re为Reynolds数; uM为由温度梯度导致的弥散相液滴的Marangoni迁移速度, R为弥散相液滴的半径, F为弥散相液滴的半径分布函数, F(R,r,t)dR给出 t时刻、位置 r处单位体积熔体内半径在 RR+dR范围内弥散相液滴的数量. 温度场方程左边第1项表示热传导的贡献, 第2项表示固/液界面上凝固潜热释放的影响, 第3项表示表面对流换热的影响, 第4项表示弥散相液滴迁移导致的热量传输; 温度场方程右边为熔体温度随时间的变化.

合金液-液相变过程中溶质通过扩散和弥散相液滴的空间运动进行传输. 根据溶质守恒关系, 组元i的浓度场满足下面偏微分方程:

(2)??[φlDi?Si]-4π3??0uMFCβi-CmiR3dR=?Cmixi?t

式中, Si=Cmi-Cmei为基体熔体中组元i的过饱和度, Cmi为基体熔体中溶质i的摩尔浓度, Cmei为溶质i在基体熔体中的平衡摩尔浓度, Cβi为弥散相液滴中溶质i的摩尔浓度, φl为基体液相的体积分数, Di为溶质i的扩散系数, Cmixi为合金(两液相混合物)中组元i的浓度. 浓度场方程左边第1项描述溶质原子在基体中扩散, 第2项描述弥散相液滴运动对浓度场的影响; 等式的右边描述合金浓度随时间的变化.

在弥散相液滴形核、扩散长大以及空间运动等因素共同作用下, 弥散相液滴半径分布函数满足如下连续方程:

(3)?F?t+??(uMF)+??R(vF)=?I?RR=R*

式中, vI分别为弥散相液滴的扩散长大速度和形核率[19,20], R*为弥散相液滴的临界形核半径. 方程左侧第1项描述弥散相液滴分布函数随时间的变化, 第2项描述弥散相液滴运动对其分布函数的影响, 第3项描述弥散相液滴长大对其分布函数的影响; 方程右侧为源项, 描述弥散相液滴形核的影响. 由于在雾化液滴快速冷却过程中弥散相液滴间的碰撞凝并作用较弱[21], 因此计算中忽略了弥散相液滴间碰撞凝并的影响.

采用有限体积法对温度场、浓度场和弥散相液滴半径分布函数连续性方程进行离散、数值求解, 并与Cu-Co-Fe热力学相图计算程序相耦合, 可模拟Cu-Co-Fe合金亚稳液-液相变过程.

2.2.2 计算中所用的主要物性参数 模拟采用文献[9,22~25]中给出的有关Cu-Co-Fe合金系的热力学数据计算相图, 计算中采用液-液界面能的通用表达式[26]:

(4)σ=σ01-TTC1.26

式中, σ为液-液界面能; TC≈(1620±10) K, 为Cu-10%Co-10%Fe合金的临界温度; 拟合计算表明, 常数 σ0=1.78 J/m2.

液态金属的动力学黏度随温度的变化关系满足下式[27]:

(5)η=η0expQRgT

式中, η0为常数, Q为流动活化能, Rg为气体常数. 富Cu液相的动力学黏度近似采用纯Cu的动力学黏度, 富Fe-Co液相的动力学黏度近似采用纯Fe和纯Co的动力学黏度的算术平均值.

溶质的扩散系数可以通过Stokes-Einstein方程与基体熔体的黏度 ηm相联系, Fe和Co溶质原子在液相Cu基体中的扩散系数 Di可如下计算[28]:

(6)Di=kBT6πηmri

式中, kB是Boltzmann常数; i=Fe, Co; ri为Fe和Co的离子半径, 2种组元的离子半径分别取为 rFe=0.055 nm和 rCo=0.53 nm. 雾化气体的初速度取为300 m/s[29].

2.2.3 模拟与讨论 根据实验条件模拟计算了雾化液滴的快速凝固过程, 图3为直径为138 μm雾化液滴中心处温度和冷却速率随时间的变化关系. 可见, 雾化液滴冷却速率的变化与雾化液滴表面对流换热系数的变化规律一致.

Fig.3 Variation of temperature (T) and cooling rate at the center of an atomized Cu-Co-Fe droplet of 138 μm in diameter with time (t) (The inset shows the variation of the convective heat-transfer coefficient (h) at the surface of an atomized Cu-Co-Fe droplet of 138 mm in diameter with time)图3 直径为138 mm的 Cu-Co-Fe雾化液滴中心处温度(T)和冷却速率随时间(t)的变化关系

图4给出了雾化液滴中心处液-液相转变的相变驱动力、富Fe-Co相液滴的形核率和数量密度随时间变化关系. 当驱动力达到一临界值时, 富Fe-Co相液滴开始形核, 发生亚稳液-液相变. 形核过程持续时间很短. 形核期间, 富Fe-Co相液滴数量密度迅速上升、平均半径缓慢增大(图5), 形核刚结束时基体熔体的过饱和度仍然很高, 富Fe-Co相液滴迅速长大, 此后, 随过饱和度的降低, 富Fe-Co相液滴长大速度逐渐减小.

Fig.4 Number density (N), nucleation rate (I) and driving force (DG) for nucleated Fe-Co-rich phase droplets at the center of an atomized drop of 138 μm in diameter vs time图4 直径为138 μm的雾化液滴中心处富Fe-Co相液滴的数量密度(N)、形核率(I)和形核驱动力(DG)随时间的变化关系

Fig.5 Average radius (Ra) of Fe-Co-rich phase droplets at the center of an atomized drop with different diameters vs time图5 不同直径雾化液滴中心处富Fe-Co相液滴的平均半径(Ra)随时间的变化关系

图6给出了直径为138 μm雾化粉末中富Fe-Co相粒子平均半径和数量密度沿粉末径向位置的变化曲线. 由于雾化液滴表面区域内温度梯度较大, 富Fe-Co相液滴Marangoni迁移作用较强, 且弥散相液滴迁移后得不到补充, 因此, 富Fe-Co相液滴持续减少, 最终导致粉末表面富Fe-Co相粒子数量密度减少(图1b); 在雾化液滴其它区域内, 在弥散相液滴形核结束后液滴的数量密度几乎不再变化, 这也表明快凝条件下富Fe-Co相液滴的Ostwald熟化作用很弱. 粉末中心绝大部分区域内富Fe-Co相粒子的空间分布非常均匀. 通过对比图5和6中计算结果和实验统计结果, 可见, 实验和模拟结果吻合得很好, 表明所建模型的合理性和有效性.

Fig.6 Variations of average radius and number density of Fe-Co-rich phase particles along the radial direction r in an atomized powder with diameter of 138 μm图6 直径为138 mm的雾化粉末中富Fe-Co相粒子的平均半径和数量密度随径向位置r变化的曲线

图7为富Fe-Co相液滴最大形核率 Imax随形核阶段冷却速率的变化关系. 雾化液滴的尺寸越小, 形核阶段雾化液滴的冷却速率越大, 熔体获得的过冷度越大, 富Fe-Co相液滴的最大形核率 Imax和液滴的数量密度越大, 液滴的平均半径越小. Imax与形核阶段冷却速率 T˙Nuc符合 Imax=BT˙Nuc2.18的指数关系, B为常数. 进一步分析表明, 当雾化粉末直径小于220 mm时, 富Fe-Co相粒子的平均半径Ra和数量密度 N与雾化粉末直径d之间分别满足 Ra=C1d0.71N=C2d-2.1的指数关系(图8), 式中, C1和C2为常数. 由这2个关系式可得弥散相体积分数 φNRa3常数, 进一步证明模型的有效性.

Fig.7 Relationship between the maximum nucleation rate (Imax) of Fe-Co-rich phase droplets and the cooling rate during the nucleation period (B—constant)图7 富Fe-Co相液滴的最大形核率Imax随形核阶段冷却速率的关系

Fig.8 Variation of average radius and number density of Fe-Co-rich particles vs powder diameter (d)图8 富Fe-Co相粒子的平均半径和数量密度随雾化粉末直径(d)变化关系

3 结论

(1) 对Cu-10%Co-10%Fe亚稳液态组元不混溶合金开展了气体雾化快速凝固实验, 建立了Cu-Co-Fe合金雾化液滴亚稳液-液相变过程的温度场、浓度场和液-液相变动力学控制方程, 并实现了合金热力学和相变动力学的耦合计算. 实验和模拟结果吻合得很好, 验证了所建模型的合理性.

(2) 在气体雾化快速凝固条件下, Cu-10%Co-10%Fe合金能够过冷至亚稳液态组元不混溶温度区而发生液-液相变. 当雾化粉末直径小于220 mm时, 富Fe-Co相液滴的凝固发生在包晶反应温度以下, 形成富Fe-Co相以球形粒子形式均匀分布于Cu基体的复合粉末组织. 雾化粉末的尺寸越小, 富Fe-Co相粒子的平均半径越小. 富Fe-Co相粒子的平均半径Ra和数量密度 N与雾化粉末直径d之间满足指数关系.

(3) 气体雾化快速凝固条件下Cu-10%Co-10%Fe合金雾化液滴除表面很小的区域外, 弥散相液滴形核后液滴的数量密度几乎不变, 富Fe-Co相液滴Marangoni迁移和Ostwald熟化作用很小, 粉末中心绝大部分区域内富Fe-Co相粒子的空间分布非常均匀.

参考文献
[1] Song J S, Hong S I. J Alloys Compd, 2000;311: 265 [本文引用:1] [JCR: 2.726]
[2] Berkowitz A E, Mitchell J R, Carey M J, Young A P, Zhang S, Spada F E, Parker F T, Hutten A, Thomas G. Phys Rev Lett, 1992;68: 3745 [本文引用:1] [JCR: 7.728]
[3] Dai F P, Cao C D, Wei B B. Chin Sci Bull, 2009;54: 402
(代富平, 曹崇德, 魏炳波. 科学通报, 2009;54: 402) [本文引用:1] [CJCR: 1.541]
[4] He J, Zhao J Z. Acta Metall Sin, 2005;41: 407
(何 杰, 赵九洲. 金属学报, 2005;41: 407) [本文引用:1] [JCR: 0.548] [CJCR: 0.835]
[5] Kim D I, Abbaschian R. J Phase Equilib, 2000;21: 25 [本文引用:3]
[6] Cao C D, G?rler G P. Chin Phys Lett, 2005;22: 482 [本文引用:1] [JCR: 0.924] [CJCR: 0.432]
[7] Curiotto S, Battezzati L, Johnson E, Palumbo M, Pryds N. J Mater Sci, 2008;43: 3253 [本文引用:1] [JCR: 2.305]
[8] Munitz A, Bamberger A M, Wannaparhun S, Abbaschian R. J Mater Sci, 2006;41: 2749 [本文引用:2] [JCR: 2.305]
[9] Turchanin M A, Dreval L A, Abdulov A R, Agraval P G. Powder Metall Met Ceram, 2011;50: 98 [本文引用:1] [JCR: 0.202]
[10] Bamberger M, Munitz A, Kaufman L, Abbaschian R. Calphad, 2002;26: 375 [本文引用:1] [JCR: 1.394]
[11] Wang C P, Liu X J, Ohnuma I, Kainuma R, Ishida K. J Phase Equilib, 2002;23: 236 [本文引用:1]
[12] Palumbo M, Curiotto S, Battezzati L. Calphad, 2006;30: 171 [本文引用:1] [JCR: 1.394]
[13] Munitz A, Abbaschian R. J Mater Sci, 1998;33: 3639 [本文引用:1] [JCR: 2.305]
[14] Dai F P, Cao C D, Wei B B. Sci Chin, 2007;37G: 376
(代富平, 曹崇德, 魏炳波. 中国科学, 2007;37G: 376) [本文引用:2]
[15] Zhou F M, Sun D K, Zhu M F. Acta Phys Sin, 2010;59: 3394
(周丰茂, 孙东科, 朱鸣芳. 物理学报, 2010;59: 3394) [本文引用:1] [JCR: 0.845] [CJCR: 1.101]
[16] Cui H B, Guo J J, Su Y Q, Wu S P, Li X Z, Fu H Z. Acta Metall Sin, 2007;43: 907
(崔红保, 郭景杰, 苏彦庆, 吴士平, 李新中, 傅恒志. 金属学报, 2007;43: 907) [本文引用:1] [JCR: 0.548] [CJCR: 0.835]
[17] Ranz W E, Marshall W R. Chem Eng Prog, 1952;48: 141 [本文引用:2] [JCR: 0.632]
[18] Grant P S, Cantor B, Katgerman L. Acta Metall Mater, 1993;41: 3097 [本文引用:1]
[19] Zhao L, Zhao J Z. J Mater Res, 2013;28: 1203 [本文引用:1] [JCR: 1.815]
[20] Granasy L, Ratke L. Scr Metall Mater, 1993;28: 1329 [本文引用:1]
[21] Zhao J Z. Mater Sci Eng, 2007;A454-455: 637 [本文引用:1] [JCR: 0.5]
[22] Dinsdale A T. Calphad, 1991;15: 317 [本文引用:1] [JCR: 1.394]
[23] Turchanin M A, Agraval P G. Powder Metall Met Ceram, 2007;46: 77 [本文引用:1] [JCR: 0.202]
[24] Turchanin M A, Agraval P G, Nikolaenko I V. J Phase Equilib, 2003;24: 307 [本文引用:1]
[25] Ohnuma I, Enoki H, Ikeda O, Kainuma R, Ohtani H, Sundman B, Ishida K. Acta Mater, 2002;50: 379 [本文引用:1] [JCR: 3.94]
[26] Moldover M R. Phys Rev, 1985;31A: 1022 [本文引用:1] [JCR: 9.478]
[27] Brandes E A, Brook G B. Smithells Metals Reference Book. 7th Ed. , Oxford: Butterworth-Heinemann Ltd, 1992: 14 [本文引用:1]
[28] Roy A K. Chhabra R P. Metall Trans, 1988;19A: 273 [本文引用:1]
[29] Moir S A, Jones H. Mater Sci Eng, 1993;A173: 161 [本文引用:1] [JCR: 0.5]