欢迎访问中国科学院大学学报,今天是
数学与物理学

铁在液态铅铋中团聚和扩散的分子动力学模拟

  • 韩雪峰 ,
  • 余新刚
展开
  • 中国科学院大学工程科学学院,北京 100049
E-mail:xgyu@ucas.ac.cn

收稿日期: 2025-03-25

  修回日期: 2025-04-27

  网络出版日期: 2025-05-26

基金资助

国家自然科学基金(12172356);国家自然科学基金(U23B20110)

Molecular dynamics simulation of the aggregation and diffusion of iron in liquid lead-bismuth

  • Xuefeng HAN ,
  • Xingang YU
Expand
  • School of Engineering Science,University of Chinese Academy of Sciences,Beijing 100049,China

Received date: 2025-03-25

  Revised date: 2025-04-27

  Online published: 2025-05-26

摘要

液态铅铋被视为第4代铅冷快堆的堆芯冷却剂,其与结构钢材的相容性问题是目前研究的重点之一。采用分子动力学方法,对铁原子在液态铅铋中的团聚过程和铁团簇的扩散行为进行模拟,分析体系温度、铁团簇尺寸和质量浓度对铁团簇之间聚合过程,以及铁团簇扩散行为的影响规律和微观机理。结果表明,铁原子在液态铅铋中会倾向于形成团簇,并且高温能够加快团聚过程。研究还发现在液态铅铋中铁团簇之间的聚合需要克服能垒,在473 K下半径为10 Å的2个铁团簇聚合需要克服的能垒约为2.14 eV,并且体系温度的升高和铁团簇尺寸的增大有利于铁团簇之间发生聚合。另外,通过计算可知,铁团簇的扩散系数随着其质量浓度的降低而增大,随着其尺寸的增大而减小,随着体系温度的升高而增大。铁团簇尺寸对扩散激活能的影响并不显著,半径为10、12.5和15 Å的铁团簇的扩散激活能分别为0.128、0.140和0.134 eV。

本文引用格式

韩雪峰 , 余新刚 . 铁在液态铅铋中团聚和扩散的分子动力学模拟[J]. 中国科学院大学学报, 2026 , 43(5) : 614 -622 . DOI: 10.7523/j.ucas.2025.028

Abstract

Liquid lead-bismuth is considered as the core coolant of the fourth-generation lead-cooled fast reactors, and its compatibility with structural steels is one of the focuses of current research. In this paper, the aggregation process of iron atoms in liquid lead-bismuth and the diffusion behavior of iron clusters were simulated by using the molecular dynamics method. The influence laws and microscopic mechanisms of the temperature of the system, the size and mass concentration of iron clusters on the coalescence process between iron clusters and the diffusion behavior of iron clusters were analyzed. The results show that iron atoms tend to form clusters in liquid lead-bismuth and that high temperatures can accelerate the aggregation process. In addition, it is found that the coalescence between iron clusters in liquid lead-bismuth needs to overcome energy barriers, the energy barrier to be overcome for the coalescence of two iron clusters of radii 10 Å at 473 K is about 2.14 eV, and the increase of system temperature and the increase of iron clusters size are favorable for the occurrence of coalescence between iron clusters. Furthermore, the diffusion coefficients of iron clusters were calculated, and found to increase with the decrease of their mass concentration, decrease with the increase of their size, and increase with the increase of the temperature of the system. The effect of iron clusters size on the diffusion activation energies was not significant, and the diffusion activation energies of iron clusters with radii of 10, 12.5, and 15 Å were 0.128, 0.140, and 0.134 eV, respectively.

参考文献

[1] Zhang J S, Li N. Review of the studies on fundamental issues in LBE corrosion[J]. Journal of Nuclear Materials2008373(1/2/3): 351-377. DOI: 10.1016/j.jnucmat.2007.06.019 .
[2] Zhang J S. A review of steel corrosion by liquid lead and lead-bismuth[J]. Corrosion Science200951(6): 1207-1227. DOI: 10.1016/j.corsci.2009.03.013 .
[3] Ballinger R G, Lim J. An overview of corrosion issues for the design and operation of high-temperature lead- and lead-bismuth-cooled reactor systems[J]. Nuclear Technology2004147(3): 418-435. DOI: 10.13182/NT04-A3540 .
[4] Wang W T, Yang C X, You Y H, et al. A review of corrosion behavior of structural steel in liquid lead-bismuth eutectic[J]. Crystals202313(6): 968. DOI: 10.3390/cryst13060968 .
[5] Zhang J S. Lead-bismuth eutectic (LBE): a coolant candidate for gen. IV advanced nuclear reactor concepts[J]. Advanced Engineering Materials201416(4): 349-356. DOI: 10.1002/adem.201300296 .
[6] Gorynin I V, Karzov G P, Markov V G, et al. Structural materials for atomic reactors with liquid metal heat-transfer agents in the form of lead or lead: Bismuth alloy[J]. Metal Science and Heat Treatment199941(9): 384-388. DOI: 10.1007/BF02469876 .
[7] Barbier F, Rusanov A. Corrosion behavior of steels in flowing lead-bismuth[J]. Journal of Nuclear Materials2001296(1/2/3): 231-236. DOI: 10.1016/S0022-3115(01)00521-9 .
[8] Hosemann P, Dickerson R, Dickerson P, et al. Transmission electron microscopy (TEM) on oxide layers formed on D9 stainless steel in lead bismuth eutectic (LBE)[J]. Corrosion Science201366: 196-202. DOI: 10.1016/j.corsci.2012.09.019 .
[9] Zhu Z G, Zhang Q, Tan J B, et al. Corrosion behavior of T91 steel in liquid lead-bismuth eutectic at 550 °C: effects of exposure time and dissolved oxygen concentration[J]. Corrosion Science2022204: 110405. DOI: 10.1016/j.corsci.2022.110405 .
[10] Song C, Li D D, Xu Y C, et al. Corrosion related properties of iron (100) surface in liquid lead and bismuth environments: A first-principles study[J]. Chinese Physics B201423(5): 056801. DOI: 10.1088/1674-1056/23/5/056801 .
[11] Xu Y C, Song C, Zhang Y G, et al. An energetic evaluation of dissolution corrosion capabilities of liquid metals on iron surface[J]. Physical Chemistry Chemical Physics201416(31): 16837-16845. DOI: 10.1039/c4cp01224k .
[12] Long X L, Shi J L, Zhu J, et al. First-principles calculation of the resistance to lead-bismuth eutectic corrosion on Fe (111) surface of austenitic stainless steel[J]. Surface Science2022725: 122132. DOI: 10.1016/j.susc.2022.122132 .
[13] Zhou T, Gao X, Ma Z W, et al. Atomistic simulation of α-Fe(100)-lead-bismuth eutectic (LBE) solid-liquid interface[J]. Journal of Nuclear Materials2021555: 153107. DOI: 10.1016/j.jnucmat.2021.153107 .
[14] Chen L M, Xu S, He X X, et al. Molecular dynamics study of corrosion behavior of iron with vacancies exposed to lead-bismuth eutectic[J]. Materials and Corrosion202374(5): 793-802. DOI: 10.1002/maco.202213375 .
[15] Gao Y, Takahashi M, Cavallotti C, et al. Molecular dynamics simulation of metallic impurity diffusion in liquid lead-bismuth eutectic (LBE)[J]. Journal of Nuclear Materials2018501: 253-260. DOI: 10.1016/j.jnucmat.2018.01.044 .
[16] Daw M S, Baskes M I. Embedded-atom method: derivation and application to impurities, surfaces, and other defects in metals[J]. Physical Review B198429(12): 6443-6453. DOI: 10.1103/physrevb.29.6443 .
[17] Belashchenko D K. Computer simulation of the properties of liquid metals: gallium, lead, and bismuth[J]. Russian Journal of Physical Chemistry A201286(5): 779-790. DOI: 10.1134/S0036024412050056 .
[18] Zhou X W, Johnson R A, Wadley H N G. Misfit-energy-increasing dislocations in vapor-deposited CoFe/NiFe multilayers[J]. Physical Review B200469(14): 144113. DOI: 10.1103/physrevb.69.144113 .
[19] Gao Y, Raos G, Cavallotti C, et al. Molecular dynamics simulation on physical properties of liquid lead, bismuth and lead-bismuth eutectic (LBE)[J]. Procedia Engineering2016157: 214-221. DOI: 10.1016/j.proeng.2016.08.359 .
[20] Nosé S. A unified formulation of the constant temperature molecular dynamics methods[J]. The Journal of Chemical Physics198481(1): 511-519. DOI: 10.1063/1.447334 .
[21] Hoover W G. Canonical dynamics: equilibrium phase-space distributions[J]. Physical Review. A, General Physics, 198531(3): 1695-1697. DOI: 10.1103/physreva.31.1695 .
[22] Parrinello M, Rahman A. Crystal structure and pair potentials: A molecular-dynamics study[J]. Physical Review Letters198045(14): 1196-1199. DOI: 10.1103/physrevlett.45.1196 .
[23] Müller-Plathe F. Reversing the perturbation in nonequilibrium molecular dynamics: an easy way to calculate the shear viscosity of fluids[J]. Physical Review. E, Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics, 199959(5 Pt A): 4894-4898. DOI: 10.1103/physreve.59.4894 .
[24] Xu B, Liu S, Yu X. Molecular simulation study of the effect of iron clusters on the viscosity of liquid lithium[J]. Journal of University of Chinese Academy of Sciences202340(2): 165-172. DOI: 10.7523/j.ucas.2021.0062 .
[25] Thompson A P, Aktulga H M, Berger R, et al. LAMMPS-a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales[J]. Computer Physics Communications2022271: 108171. DOI: 10.1016/j.cpc.2021.108171 .
[26] Stukowski A. Visualization and analysis of atomistic simulation data with OVITO-the open visualization tool[J]. Modelling and Simulation in Materials Science and Engineering201018(1): 015012. DOI: 10.1088/0965-0393/18/1/015012 .
[27] Weeks J R, Romano A J. Liquidus curves and corrosion of Fe, Ti, Zr, and Cu in liquid Bi-Pb alloys[J]. Corrosion196925(3): 131-136. DOI: 10.5006/0010-9312-25.3.131 .
[28] Xu C, Li Z Y, Chen L G, et al. Atomic-scale insights into the precipitation behaviors of copper atoms in liquid lithium[J]. Fusion Engineering and Design2023194: 113901. DOI: 10.1016/j.fusengdes.2023.113901 .
[29] Cui W Z, Shen Z J, Yang J G, et al. Effect of chaotic movements of nanoparticles for nanofluid heat transfer augmentation by molecular dynamics simulation[J]. Applied Thermal Engineering201576: 261-271. DOI: 10.1016/j.applthermaleng.2014.11.030 .
[30] Gan X L, Xiao S F, Deng H Q, et al. Clustering of Fe atoms in liquid Li and its effect on the viscosity of liquid Li[J]. Nuclear Fusion201656(4): 046004. DOI: 10.1088/0029-5515/56/4/046004 .
[31] He Y R, Jin Y, Chen H S, et al. Heat transfer and flow behaviour of aqueous suspensions of TiO2 nanoparticles (nanofluids) flowing upward through a vertical pipe[J]. International Journal of Heat and Mass Transfer200750(11/12): 2272-2281. DOI: 10.1016/j.ijheatmasstransfer.2006.10.024 .
[32] Zwanzig R. On the relation between self-diffusion and viscosity of liquids[J]. The Journal of Chemical Physics198379(9): 4507-4508. DOI: 10.1063/1.446338 .
文章导航

/