Investigation of the surface tension of LiF-CaF2-Yb2O3 molten salt and Ni-Yb alloy
-
摘要: LiF-CaF2-Yb2O3熔盐体系的表面张力是优化电解制备Ni-Yb合金的关键物理化学性质之一。本文采用拉筒法测定在1 173~1 523 K温度范围内LiF-CaF2-Yb2O3体系的表面张力,并分析其变化规律,估算Ni-Yb合金的表面张力值。结果表明:在1 173~1 573 K范围内,随着温度的升高,LiF-CaF2体系的表面张力呈线性降低;LiF-CaF2-Yb2O3体系表面张力随着温度的升高而降低;而随着Yb2O3浓度的增加,在1%~4%(质量分数)范围内,LiF-CaF2-Yb2O3体系的表面张力先增大后减小,在Yb2O3约为1% 时达到最高值;不同配比的熔融Yb-Ni合金表面张力随温度变化较小;当Yb含量为0~10%(摩尔分数)时,Ni-Yb合金的表面张力值较高,液态Ni-Yb合金易与LiF-YbF3-Yb2O3熔盐分离。Abstract: The surface tension of the LiF-CaF2-Yb2O3 molten salt system is one of the key physical and chemical properties to optimize the preparation of Ni-Yb alloys by electrolysis. In this paper, the surface tension of the LiF-CaF2-Yb2O3 system was measured by the pulling cylinder method in the temperature range of 1 173 K to 1 523 K, its changing laws were analyzed, and the surface tension of the Ni-Yb alloy was estimated. The results showed that the surface tension of the LiF-CaF2 system decreased linearly with the increase of temperature in the range of 1 173 K to 1 573 K. The surface tension of the LiF-CaF2-Yb2O3 system also decreased with the increase of temperature in the range of 1 173 K to 1 573 K, while the surface tension of the LiF-CaF2-Yb2O3 system increased at first and then decreased with the increase of Yb2O3 concentration in the range of 1% to 4% (mass fraction). The surface tension of the LiF-CaF2-Yb2O3 system reached the highest value when the content of Yb2O3 was about 1% (mass fraction) in the LiF-CaF2 system. The surface tension of the liquid Yb-Ni alloy changed little with different components and temperatures. When the content of Yb is in the range of 0 to 0.10% (mole fraction), the surface tension of the Ni-Yb alloy is higher, so liquid Ni-Yb alloy easily separates from LiF-YbF3-Yb2O3 molten salt.
-
Keywords:
- molten salt /
- Ni-Yb alloy /
- electrolysis /
- surface tension /
- pulling cylinder method
-
稀土Yb不仅对铝、镁等合金材料的力学、耐腐蚀性能的改善作用非常显著,还能与其他金属形成功能合金。开发Yb的功能合金材料,拓宽其应用领域的前景非常广阔[1-3]。目前,Yb的单质金属生产主要是通过真空镧热还原Yb2O3获得,产能低、作业不连续、成本高。同时混熔法制备Yb合金,烧损率高且易偏析,而熔盐电解法直接制备稀土金属合金能够克服热还原和混熔法的不足,提高稀土合金制备的效率[4-5]。前期研究表明采用LiF-CaF2-Yb2O3体系,Ni作为自耗阴极能制备出Ni-Yb合金,具有良好的应用前景[6]。为进一步改善电解工艺参数、降低能耗、提高电解效率,需要全面系统地研究LiF-CaF2-Yb2O3体系的物理化学性质。其中,表面张力是熔盐的一种界面性质,也是熔盐重要的物理参数之一,对界面反应以及熔盐电解过程均有较大的影响,直接关系电解过程中阳极效应以及金属产物的氧化过程[7-8]。对熔盐表面张力进行深入探讨,可为熔盐中质点间的相互作用力、熔盐电解机理、熔盐离子结构的研究提供关键数据,从而为电解工艺参数的选择提供理论依据。为此,本文采用拉筒法[9]对Yb2O3溶解度范围内的LiF-CaF2-Yb2O3体系的黏度进行测量,并通过数学模型对Ni-Yb合金表面张力进行分析和计算,从而为优化电解工艺参数提供依据,同时,为深入研究和分析LiF-CaF2-Yb2O3体系的结构及其电解机理提供必要的基础数据。
1 实验部分
分析纯CaF2、LiF、Yb2O3在423 K下烘干48 h;光谱纯石墨坩埚经乙醇清洗后,在353 K下烘干24 h。表面张力测量系统如图1所示,主要由精密电子天平与钨测头通过钢丝连接组成,当垂直的圆筒状测头与液体接触时,液体的表面张力对测头产生向下的拉力,通过测量液体表面上的圆筒测头拉离液体表面时的最大拉力(Fmax)以及圆筒周长,然后通过式(1)计算熔体的表面张力。测量流程包括:①钨测头(Φ15 mm)通过钢丝悬挂在精准电子天平上,调整电子天平和测头的位置使测头在炉膛中心处且侧头平稳没有倾斜;②通过测量高纯水在293 K下的表面张力值和在不同温度下的熔融NaCl的表面张力来校准实验装置;③校正结束后,将干燥脱水且混匀后的待测样品放入坩埚,然后利用仪器自带的高温电阻炉对样品坩埚进行加热至实验温度,并恒温20 min后,确保测头位置与校正实验时的位置相同,运行表面张力测量程序对待测样品进行表面张力测量。
(1) 式(1)中:σ为表面张力数值,单位mN/m;r为长度数值,单位m。
2 结果与讨论
2.1 LiF-CaF2-Yb2O3熔盐体系的表面张力研究
2.1.1 LiF-CaF2体系表面张力模型及数据评估
LU 等[10]及LIAO等[11]研究表明Santos方程[12]能用于二元熔盐体系表面张力的预测,且能够描述二元体系的表面张力及其表面相与体相组成的关系。通过拉筒实验测量值对Santos方程预测LiF-CaF2表面张力的适用性进行评估。Santos预测二元混合熔盐体系表面张力的模型如式(2)所示,AQRA 等[13]、HARA等[14]及MARCUS[15]研究报道的纯组元LiF和CaF2的标准表面张力值如表1所列。
表 1 纯组元LiF与CaF2的标准表面张力值Table 1. Standard surface tension values of pure components LiF and CaF2物质 表面张力/(mN/m) KT,i/GPa-1 LiF γ(LiF)=373.2-0.109T 0.093 CaF2 γ(CaF2)=459.0-0.095 6T 0.064 (2) 式(2)中:
、 分别为体相A、B的摩尔分数;σid为理想纯组元的表面张力数值,单位N/m;σA、σB分别为纯组元A、B的表面张力数值,单位N/m;VA、VB分别为纯组元A、B的摩尔体积数值,单位m3/mol; 、 分别为恒温条件下纯组元A、B摩尔表面积与压力的偏导数;R为气体常数(8.314 J/K/mol);T为绝对温度数值,单位K;KT,i为组元i的等温压缩系数数值,单位Pa-1;N0为阿伏伽德罗常数(6.02×1023 mol-1);Mi为组元i的摩尔质量数值,单位kg/mol;ρi为组元i的密度数值,单位kg/m3。 图2(a)中实验测量获得的表面张力数据表明,随着温度的升高,表面张力均可近似为呈线性降低,其中,LiF-CaF2共晶体系表面张力与温度呈线性关系如式(3)所示。因为熔盐体内部及表面离子或分子间的动能增加,间距增大,离子或分子间相互吸引力减小,分子或离子的相互作用也弱化,导致表面张力值降低。根据图2(a)中数据对比可知:Santos方程计算值与不同组分LiF-CaF2体系表面张力测量值的最大估算误差为6.55%,表明Santos模型能够较好地预测LiF-CaF2表面张力。进一步利用Santos模型分析LiF-CaF2体系表面相与体相组成关系,可以看出,表面张力值随着CaF2在表面相组成的增加而增大,温度升高导致CaF2在表面相的组成减少,CaF2相对于LiF不易与在熔体表面层表面相聚集,推测与CaF2在熔体中的结构形式有关。
图 2 LiF-CaF2熔盐体系表面张力(a)实验值与Santos模型计算值比较和(b) CaF2在表面相与体相含量关系Figure 2. (a) Comparison between the experimental values and the related values calculated by the Santos model and (b) relationship between the CaF2 content in the surface phase and bulk phase of the surface tension of the LiF-CaF2 molten salt system(3) 2.1.2 LiF-CaF2-Yb2O3体系表面张力变化规律
廖春发等[16]研究表明在1 173~1 523 K温度范围内,(LiF-CaF2)etu体系中Yb2O3溶解度低于2%(质量分数),因此,图3(a)中为一定Yb2O3含量下 (LiF-CaF2)etu -Yb2O3体系表面张力随温度的变化规律,可以看出,在1 173~1 523 K范围内, 一定Yb2O3含量的(LiF-CaF2)etu -Yb2O3体系表面张力随着温度的升高而降低,原因是温度的升高致使熔体表面层离子或团簇具备较高的动能和较低的解离能,阴阳离子之间的相互吸引作用减弱。图3(b)为特定温度下(LiF-CaF2)etu-Yb2O3体系表面张力随Yb2O3含量的变化规律,可以看出,表面张力与Yb2O3的溶解度之间密切相关。当Yb2O3的质量分数为0~2%时,Yb2O3在LiF-CaF2体系内被充分溶解,体系的表面张力在Yb2O3的质量分数为1%时达到最高值。随着Yb2O3的质量分数进一步增加,达到过饱和状态,体系的表面张力出现减小,分析原因认为Yb2O3的加入会导致(LiF-CaF2)etu体系的离子构成形式发生变化,从而使表面层和体相离子团簇构成形式发生变化。由于(LiF-CaF2)etu体系中Yb2O3的溶解度较小,(LiF-CaF2)etu体系中主体结构不会发生根本性变化,但可以确定在Yb2O3的溶解量为0~2%(质量分数)时,体系内部会产生吸引作用力强的离子团簇,使其表面张力增加。随着Yb2O3含量的增加,体系表面层中离子结构发生变化,削弱了表面层离子或团簇的吸引作用,导致表面张力下降。拟合图3(b)中数据可以得到LiF-CaF2-Yb2O3熔盐体系表面张力的回归方程式(4)。
(4) 2.2 Ni-Yb合金的表面张力
综合LiF-CaF2-Yb2O3体系和Ni-Yb合金表面张力变化规律,表明随电解温度的升高,LiF-CaF2-Yb2O3熔盐体系的表面张力减小,但仅在220~260 mN/m较小范围内变化。Yb2O3在LiF-CaF2-Yb2O3体系中的含量引起的张力变化范围同样不显著,在1 173~1 523 K温度区间表面张力波动范围小于20 mN/m。因此,Ni-Yb合金的表面张力值随组分和温度的变化决定合金与LiF-CaF2-Yb2O3体系的分离。
Butler方程[17]如式(5)所示,该方程经过完善后被广泛应用于二元及三元合金内的表面张力的预测[18-19],经对多种合金体系的验证,预测值与实验值较吻合。利用Butler模型,结合Yb、Ni纯金属的标准物理参数[20],见表2所列,得到图4所示高于Ni-Yb合金熔点100~600 K和Yb摩尔百分含量0~100%范围内的表面张力曲面预测值。 同时,图4中给出了LiF-CaF2-Yb2O3体系在1 173~1 523 K温度区间以及Yb2O3质量分数在0~4%范围内的表面张力曲面预测值,当合金中Yb的摩尔分数低于10%时, 合金的表面张力变化梯度大, 在600~ 1 700 mN/m迅速增大。当Yb的摩尔分数高于10%时,合金的表面张力变化梯度小,在200~600 mN/m范围内缓慢减小。对比LiF-CaF2-Yb2O3体系的张力曲面,当合金中Yb的摩尔分数低于10%时,合金与熔盐的表面张力差值范围在340~1 480 mN/m间波动,远高于熔盐体系。合金与熔盐体系润湿性差,液态合金易于团聚。合金中Yb的摩尔分数高于10%时,合金与熔盐的表面张力差值仅在20~340 mN/m小范围内波动,合金与熔盐体系润湿性好,液态合金不易团聚。总体上, 在温度1 173~1 523 K范围内的LiF-CaF2-Yb2O3体系下,熔盐体系的黏度相对稳定,Yb摩尔分数高于10%的Ni-Yb合金产物与LiF-CaF2-Yb2O3熔盐表面张力差值较小,润湿性更好,有利于收集。
表 2 纯金属Ni、Yb的标准物理参数Table 2. Standard physical parameters for the pure metals Ni and Yb组元 表面张力(σ)/(mN/m) Ni 1 834-0.376×(T-1 455) Yb 320-0.102×(T-824) (5) 式(5)中:
=0.75,σ和σi 为溶体和纯组分i的表面张力数值,单位mN/m;R为气体常数(8.314 J/mol/K);T为绝对温度数值,单位K;XiS和XiB分别为组分i在表面相和体相中的摩尔分数;Si为组分i纯物质的单层表面积数值,单位m2; 和 分别为表面相和体相中组分i的偏摩尔吉布斯自由能数值,单位J/mol;ZS和ZB分别是表面相和体相的配位数; No为阿伏伽德罗常数(6.02×1023 mol-1);Mi为组分i的摩尔质量数值,单位kg;ρi为密度数值,单位kg/m3。 3 结 论
1)在1 173~1 523 K范围内,随着温度的升高,LiF-CaF2体系的表面张力呈线性降低。Santos方程能较好预测LiF-CaF2体系二元熔盐体系表面张力和描述表面相与体相组成的关系;在1 173~1 523 K范围内,LiF-CaF2-Yb2O3体系表面张力随温度的升高而降低;在Yb2O3的质量分数为1%~4%范围内,体系的表面张力先增后减,在质量分数为1%时达到最高。
2)当高于熔点200~600 K范围内,液态Ni-Yb合金中Yb的摩尔分数低于10%时,其表面张力在600~1 700 mN/m迅速增大。当Yb的摩尔分数高于10%时,合金的表面张力在200~600 mN/m范围内缓慢减小。温度为1 173~1 523 K范围内的LiF-CaF2-Yb2O3体系的表面张力在220~260 mN/m较小范围内波动,相对稳定,与Yb摩尔分数高于10%的Ni-Yb合金表面张力差值较小,润湿性更好。
-
图 2 LiF-CaF2熔盐体系表面张力(a)实验值与Santos模型计算值比较和(b) CaF2在表面相与体相含量关系
Fig 2. (a) Comparison between the experimental values and the related values calculated by the Santos model and (b) relationship between the CaF2 content in the surface phase and bulk phase of the surface tension of the LiF-CaF2 molten salt system
表 1 纯组元LiF与CaF2的标准表面张力值
Table 1 Standard surface tension values of pure components LiF and CaF2
物质 表面张力/(mN/m) KT,i/GPa-1 LiF γ(LiF)=373.2-0.109T 0.093 CaF2 γ(CaF2)=459.0-0.095 6T 0.064 表 2 纯金属Ni、Yb的标准物理参数
Table 2 Standard physical parameters for the pure metals Ni and Yb
组元 表面张力(σ)/(mN/m) Ni 1 834-0.376×(T-1 455) Yb 320-0.102×(T-824) -
[1] TUAN N Q, ALVES A C, TOPTAN F, et al. Effects of substituting ytterbium for scandium on corrosion behaviour of Al-Sc alloy[J]. Materials and Corrosion, 2015, 66(12):1504-1511.
[2] TANG C L, ZHOU D J. Precipitation hardening behavior of dilute binary Al-Yb alloy[J]. Transactions of Nonferrous Metals Society of China, 2014, 24: 2326-2330.
[3] HODGES J A, BONVILLE P, OCIO M. Magnetic properties of YbNi5 from 170Yb Mossbauer and magnetisation measurements[J]. The European Physical Journal B, 2007, 57: 365-370.
[4] 王玉香, 赖华生, 文小强, 等. 氟化物体系熔盐电解制备YNi合金[J]. 有色金属科学与工程, 2021,12(1):126-130. [5] 郭探, 王世栋, 叶秀深, 等. 熔盐电解法制备稀土合金研究进展[J]. 中国科学:化学,2012,42(9): 1328-1336. [6] 王旭, 廖春发, 王瑞祥, 等. 氟化物介质熔盐电解制备Ni-Yb合金及其表征[J]. 材料导报, 2019, 33(5):750-753 [7] DANEK V. Physico-Chemical Analysis of Molten Electrolytes[M]. Netherlands: Institute of Inorganic Chemistry, 2006.
[8] 康佳, 闫奇操, 于兵, 等. LaCl3-KCl熔盐体系物化性质的研究[J]. 有色金属科学与工程, 2022,13(3):145-151 [9] 范建峰, 袁章福, 柯家骏. 高温熔体表面张力测量方法的进展[J]. 化学通报, 2004(11):802-807. [10] LU G X, HE M G, KANG Z Y. Study on surface tension and surface phase of molten MCl-CeCl3 (M=Li, Na, K, Cs): A comparison of Butler's equation and a newly-proposed model[J]. Fluid Phase Equilib, 2016, 427:345-352.
[11] LIAO C F, CAI B Q, WANG X, et al. Study on surface tension of Dy-Cu intermediate alloy prepared using LiF-DyF3 molten salt system[J]. Journal of Molecular Liquids, 2020, 301:1-11.
[12] SANTOS M, S C S, REIS J C R. Activity coefficients in the surface phase of liquid mixtures[J]. Chemphyschem, 2015, 16(2):470-475.
[13] AQRA F. The cohesive energy density and the isothermal compressibility: Their relationships with the surface tension[J]. Physica B: Condensed Matter, 2014, 446:28-31.
[14] HARA S, SHIHIBAIKE H, OGINO K. The molar volumes and the surface tensions of melts in the systems CaF2-MF (M: Li, Na) and CaF2-MF2 (M: Mg, Sr, Ba)[J]. ISIJ International, 2007, 30(4):298-304.
[15] MARCUS Y. The compressibility and surface tension product of molten salts[J]. The Journal of Chemical Physical, 2013, 139(12):1-4.
[16] 廖春发, 孙强超, 王旭,等. Yb2O3在LiF-CaF2熔盐中的溶解度研究[J]. 中国稀土学报, 2019(1):99-104. [17] BUTLER J A V. The thermodynamics of the surfaces of solutions[J]. Proceedings of the royal society of London A, 1932, 135(827):348-375.
[18] YAN L J, QIAO Z Y. Surface tension modeling of molten rare earth containing halide mixtures[J]. Journal of Rare Earths, 2000, 18(4):308-311.
[19] MEHTA U, KOIRALA I, YADAV S K, et al. Prediction of thermodynamic and surface properties of ternary Ti-Si-Fe liquid alloy[J]. Modelling and Simulation Materials Science and Engineering, 2020, 28(6):1-16.
[20] ZHU D M, LI C R, GUO C P, et al. Thermodynamic assessment of Ni-Yb binary system[J]. Transactions of Nonferrous Metals Society of China, 2016, 26(4):1055-1062.