摘要
对碱金属氟化物LiF、NaF、KF熔盐黏度进行分子动力学计算,并利用实测值对计算结果进行评估,结果表明:在范德华势的基础上,引入核壳模型库仑势,并加入极化能量后作为势函数,依据能量耗散原理,通过Gromacs软件对LiF、NaF、KF熔盐在温度1 100~1 700 K范围内选点进行黏度计算,能够准确表达温度、碱金属离子半径与黏度的定性关系;计算的结果有系统性负偏差,主要原因及解决方案包括:(1)LiF、NaF、KF熔盐体系黏度描述模型采用严格的牛顿流体模型,而实际LiF、NaF、KF熔盐体系为似牛顿流体,需要对体系离子加速度方程添加修正项;(2)库仑常数(β)、极化率(α)以及核电荷量(qcore)是决定计算精度的关键参数,库仑常数的基础值以及合理的核电荷量和极化率的匹配关系对计算结果有显著影响。
- 单位