LiF-CaF2-Yb2O3熔盐及Ni-Yb合金表面张力研究

张永健， 蔡伯清， 王旭*， 王瑞祥， 石忠宁

（1. 江西理工大学材料冶金化学学部，江西 赣州 341000； 2. 东北大学冶金工程学院，沈阳 110006）

稀土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 表面张力测量装置示意Fig.1 Schematic diagram of the surface tension measuring device

式（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所列。

式（2）中：、分别为体相A、B 的摩尔分数；σid为理想纯组元的表面张力数值，单位N/m；σA、σB分别为纯组元A、B 的表面张力数值，单位N/m；VA、VB分别为纯组元A、B 的摩尔体积数值，单位m3/mol；(∂SA∂p)T、(∂SB∂p)T分别为恒温条件下纯组元A、B摩尔表面积与压力的偏导数；R为气体常数（8.314 J/K/mol）；T为绝对温度数值，单位K；KT,i为组元i 的等温压缩系数数值，单位Pa-1；N0为阿伏伽德罗常数（6.02×1023mol-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在表面相与体相含量关系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

2.1.2 LiF-CaF2-Yb2O3体系表面张力变化规律

廖春发等[16]研究表明在1 173～1 523 K 温度范围内，（LiF-CaF2）etu体系中Yb2O3溶解度低于2%（质量分数），因此，图3（a）中为一定Yb2O3含量下 （LiFCaF2）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的溶解度较小，（LiFCaF2）etu体系中主体结构不会发生根本性变化，但可以确定在Yb2O3的溶解量为0～2%（质量分数）时，体系内部会产生吸引作用力强的离子团簇，使其表面张力增加。随着Yb2O3含量的增加，体系表面层中离子结构发生变化，削弱了表面层离子或团簇的吸引作用，导致表面张力下降。拟合图3（b）中数据可以得到LiF-CaF2-Yb2O3熔盐体系表面张力的回归方程式（4）。

图3 温度（a）和Yb2O3含量（b）对LiF-CaF2熔盐体系表面张力的影响Fig.3 Effect of temperature （a） and Yb2O3 content （b） on the surface tension of the LiF-CaF2 molten salt system

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

图4 Ni-Yb合金与LiF-CaF2-Yb2O3体系表面张力曲面对比Fig.4 Comparison diagram of the surface tension curved surface of the LiF-CaF2-Yb2O3 system and Ni-Yb alloy

式（5）中：β=0.75，σ和σi为溶体和纯组分i 的表面张力数值，单位mN/m；R为气体常数（8.314 J/mol/K）；T为绝对温度数值，单位K；和分别为组分i 在表面相和体相中的摩尔分数；Si为组分i 纯物质的单层表面积数值，单位m2；和分别为表面相和体相中组分i 的偏摩尔吉布斯自由能数值，单位J/mol；ZS和ZB分别是表面相和体相的配位数；No为阿伏伽德罗常数（6.02×1023mol-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 范围内的LiFCaF2-Yb2O3体系的表面张力在220～260 mN/m 较小范围内波动，相对稳定，与Yb 摩尔分数高于10%的Ni-Yb合金表面张力差值较小，润湿性更好。