原子动力学蒙特卡洛程序OpenKMC在反应堆压力容器钢缺陷损伤研究中的优化与应用

计算机工程与科学 ›› 2020, Vol. 42 ›› Issue (12): 2151-2162.

原子动力学蒙特卡洛程序OpenKMC在反应堆压力容器钢缺陷损伤研究中的优化与应用

尚子豪1,2,商红慧2,王东杰3,张云泉2,贺新福3,陈泽华1,王栋2,张广婷2

（1.太原理工大学大数据学院，山西太原 030024；

2.中国科学院计算技术研究所计算机体系结构国家重点实验室，北京 100190；

3.中国原子能科学研究院，北京 102413）

收稿日期:2020-06-08 修回日期:2020-07-15 接受日期:2020-12-25 出版日期:2020-12-25 发布日期:2021-01-05
基金资助:
国家重点研发计划(2017YFB0202302)

Optimization and application of atomic Kinetic Monte Carlo program OpenKMC in defect damage of reactor pressure vessel steel

SHANG Zi-hao1,2 ,SHANG Hong-hui2,WANG Dong-jie3,ZHANG Yun-quan2,HE Xin-fu3,CHEN Ze-hua1,WANG Dong2,ZHANG Guang-ting2

(1.College of Big Data,Taiyuan University of Technology,Taiyuan 030024;

2.State Key Laboratory of Computer Architecture，Institute of Computing Technology，

Chinese Academy of Sciences，Beijing 100190;

3.China Institute of Atomic Energy，Beijing 102413，China)

Received:2020-06-08 Revised:2020-07-15 Accepted:2020-12-25 Online:2020-12-25 Published:2021-01-05

摘要/Abstract

摘要： 将具有体心立方结构的Fe-Cu二元合金作为RPV模型材料，分别采用基于Pair势和嵌入原子势(EAM)的动力学蒙特卡洛方法，以引入空位点缺陷的方式模拟了热时效下系统中富铜团簇的析出过程。同时对程序实现了计算优化，并验证了算法的正确性与有效性。利用高性能计算资源对优化后的程序进行了性能分析。数值实验结果表明，通过引入一定数量的空位点缺陷，体系内可以同时析出富铜团簇和铜-空位复合体团簇，并且这种复合体团簇更易成为体系内尺寸最大的团簇。增加体系内的空位数量可以加快沉淀的析出过程。另外，增加空位数量不会对体系总的团簇数量密度产生显著影响，但是能够促进团簇粗化，使其成长为尺寸更大的沉淀。

关键词: Fe-Cu合金, Pair势, 嵌入原子势, 团簇, 空位点缺陷, 动力学蒙特卡洛

Abstract: The Fe-Cu binary alloy with BCC structure is used as the RPV simulation steel. The kinetic Monte Carlo method based on Pair potential and EAM potential is adopted to simulate the precipitation process of Cu-rich clusters in the system under thermal aging by introducing vacancy point defects. At the same time, the program is optimized, and the correctness and validity of the algorithm are verified. The performance of the optimized program is analyzed by using high performance computing resources. The numerical results show that, by introducing a certain number of vacancy point defects, Cu-rich clusters and Cu-vacancy complex clusters can be precipitated simultaneously in the system, and the complex clusters are more likely to become the largest clusters in the system. The precipitation process can be accelerated by increasing the number of vacancies in the system. In addition, increasing the number of vacancies does not have a significant impact on the total number and density of clusters in the system, but it can promote the coarsening of clusters and make them grow into larger size precipitates.

Key words: Fe-Cu alloy, pair , potential, embedded-auto method, cluster, vacancy point defect, kinetic Monte Carlo

尚子豪, 商红慧, 王东杰, 张云泉, 贺新福, 陈泽华, 王栋, 张广婷. 原子动力学蒙特卡洛程序OpenKMC在反应堆压力容器钢缺陷损伤研究中的优化与应用[J]. 计算机工程与科学, 2020, 42(12): 2151-2162.

SHANG Zi-hao, , SHANG Hong-hui, WANG Dong-jie, ZHANG Yun-quan, HE Xin-fu, CHEN Ze-hua, WANG Dong, ZHANG Guang-ting. Optimization and application of atomic Kinetic Monte Carlo program OpenKMC in defect damage of reactor pressure vessel steel[J]. Computer Engineering & Science, 2020, 42(12): 2151-2162.