Tian Heqing, Zhang Wenguang, Guo Chaxiu
School of Mechanical and Power Engineering, Zhengzhou University, Zhengzhou 450001, China.
J Chem Phys. 2024 Oct 14;161(14). doi: 10.1063/5.0219401.
Thermal energy storage and utilization has been widely concerned due to the intermittency, renewability, and economy of renewable energy. In this paper, the potential energy function of binary Na2CO3-K2CO3 salt was first constructed using the Deep Potential GENerator (DPGEN) enhanced sampling method. Deep potential molecular dynamics simulations were performed to calculate the thermal properties and structural evolution of binary carbonates. The results show that as the temperature increases from 1073 to 1273 K, the viscosity and thermal conductivity decrease from 5.011 mPa s and 0.502 W/(m K) to 2.526 mPa s and 0.481 W/(m K), respectively. The decrease in viscosity is related to the distance and interaction between the molten salt ions. In addition, the diffusion coefficients, energy barriers, ionic radius, angular distribution function, and coordination number of molten salt were calculated and analyzed. The CO32- exhibits a stable planar triangular structure. The ionic radius of Na+ is smaller than that of K+, which makes Na+ suffer less spatial hindrance during motion and has a higher diffusion coefficient. The energy barriers that Na+ needs to overcome to escape the Coulomb force is greater than that of K+ ions, so molten salt containing Na+ may possess greater heat storage potential. We believe that the potential function constructed with DPGEN enhanced sampling strategy can provide more convincing results for predicting the thermal properties of molten salts. This paper aims to provide a technical route to develop the novel complex molten salt phase change material for thermal energy storage.
由于可再生能源的间歇性、可再生性和经济性,热能存储与利用受到了广泛关注。本文首先采用深度势能生成器(DPGEN)增强采样方法构建了二元Na2CO3-K2CO3盐的势能函数。进行了深度势能分子动力学模拟,以计算二元碳酸盐的热性质和结构演变。结果表明,随着温度从1073 K升高到1273 K,粘度和热导率分别从5.011 mPa·s和0.502 W/(m·K)降至2.526 mPa·s和0.481 W/(m·K)。粘度的降低与熔盐离子之间的距离和相互作用有关。此外,还计算并分析了熔盐的扩散系数、能垒、离子半径、角分布函数和配位数。CO32-呈现出稳定的平面三角形结构。Na+的离子半径小于K+,这使得Na+在运动过程中受到的空间位阻较小,扩散系数较高。Na+逃离库仑力所需克服的能垒大于K+离子,因此含Na+的熔盐可能具有更大的蓄热潜力。我们认为,采用DPGEN增强采样策略构建的势能函数能够为预测熔盐的热性质提供更具说服力的结果。本文旨在提供一条开发新型复合熔盐相变储热材料的技术路线。