The reactor is an extreme environment of high temperature, high pressure, and high radiation dose. Damage accumulates in structural materials over service time. It will eventually lead to material failure, such as hardening, embrittlement, and swelling. The radiation damage accumulation is an inherent, complex, and multiscale process, which has been studied extensively with multiscale modeling and simulation methods. The rapid development of high-performance computing makes it possible to accurately operate multiscale simulation for the microstructure evolution of irradiated reactor materials. The European Union, the USA, and China have put great effort into this, and many related works have been carried out. This paper first outlines the basic application of multiscale modeling and simulation technology in understanding the effects of radiation on reactor structural materials. Then, some relevant projects carried out by the USA, the European Union, and China in recent years are summarized. Next, the paper focuses on three widely used simulation techniques at different scales: molecular dynamics, kinetic Monte Carlo, and cluster dynamics. For each method, some key developments in algorithms and computer implementations are reviewed. Finally, the comparison between them is discussed.
GuérinYWasGSZinkleSJ.Materials challenges for advanced nuclear energy systems. Mrs Bull2009; 34: 10–19.
2.
MarquesJ.Evolution of nuclear fission reactors: third generation and beyond. Energ Convers Manag2010; 51: 1774–1780.
3.
MascittiJMadariagaM.Method for the calculation of DPA in the reactor pressure vessel of Atucha II. Sci Technol Nucl Install2011; 2011: 1–6.
4.
DunnAY.Radiation damage accumulation and associated mechanical hardening in thin films and bulk materials. PhD Thesis, Georgia Institute of Technology, 2016.
5.
VladimirovPBouffardS.Displacement damage and transmutations in metals under neutron and proton irradiation. Comptes Rendus Physique2008; 9: 303–322.
6.
AzevedoCRdF. Selection of fuel cladding material for nuclear fission reactors. Eng Fail Anal2011; 18: 1943–1962.
7.
OdetteGWirthBBaconD, et al. Multiscale-multiphysics modeling of radiation-damaged materials: embrittlement of pressure-vessel steels. Mrs Bull2001; 26: 176–181.
8.
WirthBOdetteGMarianJ, et al. Multiscale modeling of radiation damage in Fe-based alloys in the fusion environment. J Nucl Mater2004; 329: 103–111.
9.
KnasterJMoeslangAMurogaT.Materials research for fusion. Nat Phys2016; 12: 424.
10.
DownarTSiegelAUnalC, et al. Science based nuclear energy systems enabled by advanced modeling and simulation at the extreme scale. Technical report, 2009.
11.
BishopAMessinaP, et al. Scientific grand challenges in national security: the role of computing at the extreme scale. Technical report, 2009.
WirthBDHuXKohnertA, et al. Modeling defect cluster evolution in irradiated structural materials: focus on comparing to high-resolution experimental characterization studies. J Mater Res2015; 30: 1440–1455.
14.
HafnerJ.Ab-initio simulations of materials using VASP: density-functional theory and beyond. J Computat Chem2008; 29: 2044–2078.
15.
MartínezECaturlaMJMarianJ.DFT-parameterized object kinetic Monte Carlo simulations of radiation damage. In: AndreoniWYipS (eds) Handbook of materials modeling: applications: current and emerging materials, 2018, pp.1–32. Berlin: Springer.
16.
HafnerJWolvertonCCederG.Toward computational materials design: the impact of density functional theory on materials research. Mrs Bull2006; 31: 659–668.
17.
KulikovDMalerbaLHouM.On the binding energies and configurations of vacancy and copper–vacancy clusters in bcc Fe–Cu: a computational study. Philos Mag2006; 86: 141–172.
18.
DjurabekovaFMalerbaLDomainC, et al. Stability and mobility of small vacancy and copper-vacancy clusters in bcc-Fe: an atomistic kinetic Monte Carlo study. Nucl Instrum Meth Phys Res B Beam Interact Mater Atom2007; 255: 47–51.
19.
WangLW.A survey of codes and algorithms used in NERSC material scienceallocations. Technical report, Ernest Orlando Lawrence Berkeley National Laboratory, Berkeley, CA (US), 2006.
20.
NordlundK.Historical review of computer simulation of radiation effects in materials. J Nucl Mater2019; 520: 273–295.
21.
PatelMKTallmanDJValdezJA, et al. Effect of helium irradiation on Ti3AlC2 at 500 °C. Scripta Materialia2014; 77: 1–4.
22.
TrinkausHSinghB.Helium accumulation in metals during irradiation–where do we stand?J Nucl Mater2003; 323: 229–242.
23.
YoshidaNIwakiriHTokunagaK, et al. Impact of low energy helium irradiation on plasma facing metals. J Nucl Mater2005; 337: 946–950.
24.
LasaAHenrikssonKNordlundK.MD simulations of onset of tungsten fuzz formation under helium irradiation. Nucl Instrum Meth Phys Res B Beam Interact Mater Atom2013; 303: 156–161.
25.
IwakiriHYasunagaKMorishitaK, et al. Microstructure evolution in tungsten during low-energy helium ion irradiation. J Nucl Mater2000; 283: 1134–1138.
Rapaport DC. The art of molecular dynamics simulation. Cambridge: Cambridge University Press, 2004.
28.
MalerbaL.Molecular dynamics simulation of displacement cascades in -Fe: a critical review. J Nucl Mater2006; 351: 28–38.
29.
SetyawanWNandipatiGRocheKJ, et al. Displacement cascades and defects annealing in tungsten, Part I: defect database from molecular dynamics simulations. J Nucl Mater2015; 462: 329–337.
30.
FuCCDalla TorreJWillaimeF, et al. Multiscale modelling of defect kinetics in irradiated iron. Nat Mater2005; 4: 68.
31.
SonedaNIshinoSTakahashiA, et al. Modeling the microstructural evolution in bcc-Fe during irradiation using kinetic Monte Carlo computer simulation. J Nucl Mater2003; 323: 169–180.
32.
BecquartCSBarbuABocquetJ, et al. Modeling the long-term evolution of the primary damage in ferritic alloys using coarse-grained methods. J Nucl Mater2010; 406: 39–54.
33.
JanssonVMalerbaL.Simulation of the nanostructure evolution under irradiation in Fe-C alloys. J Nucl Mater2013; 443: 274–285.
34.
VoterAF.Introduction to the kinetic Monte Carlo method. Berlin: Springer, 2007, pp.1–23.
35.
GhoniemNM.Clustering theory of atomic defects. Radiat Effect Defect Solid1999; 148: 269–318.
36.
DuparcAHMoingeonCSmetniansky-de GrandeN, et al. Microstructure modelling of ferritic alloys under high flux 1 MeV electron irradiations. J Nucl Mater2002; 302: 143–155.
37.
SemenovAWooC.Phase-field modeling of void formation and growth under irradiation. Acta Materialia2012; 60: 6112–6119.
38.
LiYHuSSunX, et al. A review: applications of the phase field method in predicting microstructure and property evolution of irradiated nuclear materials. npj Computat Mater2017; 3: 16.
39.
TonksMRGastonDMillettPC, et al. An object-oriented finite element framework for multiphysics phase field simulations. Computat Mater Sci2012; 51: 20–29.
40.
JuLZhangJDuQ.Fast and accurate algorithms for simulating coarsening dynamics of Cahn–Hilliard equations. Computat Mater Sci2015; 108: 272–282.
41.
LiYHuSZhangL, et al. Non-classical nuclei and growth kinetics of cr precipitates in FeCr alloys during ageing. Model Simulat Mater Sci Eng2014; 22: 025002.
42.
BinerSRaoWZhangY.The stability of precepitates and the role of lattice defects in Fe-1at% Cu-1at% Ni-1at% Mn alloy: a phase-field model study. J Nucl Mater2016; 468: 9–16.
43.
ShortMGastonDJinM, et al. Modeling injected interstitial effects on void swelling in self-ion irradiation experiments. J Nucl Mater2016; 471: 200–207.
44.
BelytschkoTGracieRVenturaG.A review of extended/generalized finite element methods for material modeling. Model Simulat Mater Sci Eng2009; 17: 043001.
45.
EtienneAHernández-MayoralMGenevoisC, et al. Dislocation loop evolution under ion irradiation in austenitic stainless steels. J Nucl Mater2010; 400: 56–63.
46.
OkaHWatanabeMKinoshitaH, et al. In situ observation of damage structure in ODS austenitic steel during electron irradiation. J Nucl Mater2011; 417: 279–282.
47.
JitsukawaSKimuraAKohyamaA, et al. Recent results of the reduced activation ferritic/martensitic steel development. J Nucl Mater2004; 329: 39–46.
48.
RothJTsitroneELoarteA, et al. Recent analysis of key plasma wall interactions issues for ITER. J Nucl Mater2009; 390: 1–9.
49.
HiraiTEscourbiacFCarpentier-ChouchanaS, et al. ITER tungsten divertor design development and qualification program. Fusion Eng Des2013; 88: 1798–1801.
50.
UedaYCoenenJDe TemmermanG, et al. Research status and issues of tungsten plasma facing materials for ITER and beyond. Fusion Eng Des2014; 89: 901–906.
TurnerJAClarnoKSiegerM, et al. The virtual environment for reactor applications (VERA): design and architecture. J Computat Phys2016; 326: 544–568.
53.
JumelSDomainCRusteJ, et al. Simulation of the irradiation effects in reactor materials: the REVE project. Journal de Physique IV2000; 10: Pr6–191.
54.
JumelSDomainCRusteJ, et al. Simulation of irradiation effects in reactor pressure vessel steels: the Reactor for Virtual Experiments (REVE) project. J Test Eval2002; 30: 37–46.
55.
JumelSVan-DuysenJC.RPV-1: a virtual test reactor to simulate irradiation effects in light water reactor pressure vessel steels. J Nucl Mater2005; 340: 125–148.
56.
JumelS.RPV-1: a first virtual reactor to simulate irradiation effects in light water reactor pressure vessel steels. PhD Thesis, Université des Sciences et Technologies de Lille I, 2005.
57.
AdjanorGBugatSDomainC, et al. Overview of the RPV-2 and INTERN-1 packages: from primary damage to microplasticity. J Nucl Mater2010; 406: 175–186.
58.
BecquartCDeckerKDomainC, et al. Massively parallel molecular dynamics simulations with EAM potentials. Radiat Effect Defect Solid1997; 142: 9–21.
59.
DomainCBecquartCMalerbaL.Simulation of radiation damage in Fe alloys: an object kinetic Monte Carlo approach. J Nucl Mater2004; 335: 121–145.
60.
OdetteGWirthB.A computational microscopy study of nanostructural evolution in irradiated pressure vessel steels. J Nucl Mater1997; 251: 157–171.
61.
LiuCOdetteGWirthB, et al. A lattice Monte Carlo simulation of nanophase compositions and structures in irradiated pressure vessel Fe-Cu-Ni-Mn-Si steels. Mater Sci Eng A1997; 238: 202–209.
62.
LeclercqSLidburyDVan DyckS, et al. PERFORM 60–prediction of the effects of radiation for reactor pressure vessel and in-core materials using multi-scale modelling–60 years foreseen plant lifetime. J Nucl Mater2010; 406: 193–203.
63.
Al MazouziAAlamoALidburyD, et al. PERFORM 60: prediction of the effects of radiation for reactor pressure vessel and in-core materials using multi-scale modelling–60 years foreseen plant lifetime. Nucl Eng Des2011; 241: 3403–3415.
64.
FazioCBricenoDGRiethM, et al. Innovative materials for Gen IV systems and transmutation facilities: the cross-cutting research project GETMAT. Nucl Eng Des2011; 241: 3514–3520.
65.
FazioCAlamoAAlmazouziA, et al. European cross-cutting research on structural materials for generation IV and transmutation systems. J Nucl Mater2009; 392: 316–323.
TchipevNSecklerSHeinenM, et al. Twetris: twenty trillion-atom simulation. Int J High Perform Comput Appl2018; 33: 838–854.
69.
EckhardtWHeineckeABaderR, et al. 591 TFLOPS multi-trillion particles simulation on SuperMUC. In: MartinJLudwigKTMeuerHW (eds) international supercomputing conference, 16–20 June 2013, pp. 1–13. Berlin: Springer.
VoterAF.Parallel replica method for dynamics of infrequent events. Phys Rev B1998; 57: R13985.
73.
SϕrensenMRVoterAF. Temperature-accelerated dynamics for simulation of infrequent events. J Chem Phys2000; 112: 9599–9606.
74.
FanYKushimaAYipS, et al. Mechanism of void nucleation and growth in bcc Fe: atomistic simulations at experimental time scales. Phys Rev Lett2011; 106: 125501.
75.
UberuagaBHoaglandRVoterA, et al. Direct transformation of vacancy voids to stacking fault tetrahedra. Phys Rev Lett2007; 99: 135501.
76.
BaiXMVoterAFHoaglandRG, et al. Efficient annealing of radiation damage near grain boundaries via interstitial emission. Science2010; 327: 1631–1634.
77.
HuCWangXLiJ, et al. Kernel optimization for short-range molecular dynamics. Comput Phys Commun2017; 211: 31–40.
78.
TarmyshovKBMüller-PlatheF.Parallelizing a molecular dynamics algorithm on a multiprocessor workstation using OpenMP. J Chem Inform Model2005; 45: 1943–1952.
79.
KunasethMRichardsDFGlosliJN, et al. Analysis of scalable data-privatization threading algorithms for hybrid MPI/OpenMP parallelization of molecular dynamics. J Supercomput2013; 66: 406–430.
80.
HockneyRGoelSEastwoodJ.Quiet high-resolution computer models of a plasma. J Computat Phys1974; 14: 148–158.
81.
KadauKGermannTCLomdahlPS.Molecular dynamics comes of age: 320 billion atom simulation on BlueGene/L. Int J Mod Phys C2006; 17: 1755–1761.
HuCBaiHHeX, et al. Crystal MD: The massively parallel molecular dynamics software for metal with BCC structure. Comput Phys Commun2017; 211: 73–78.
84.
BaiHHuCHeX, et al. Crystal MD: molecular dynamic simulation software for metal with BCC structure. In: ChenWYinGZhaoGHanQJingWSunGLuZ (eds) national conference on big data technology and applications, Harbin, 25–26 December 2015, pp.247–258. Berlin: Springer.
85.
StadlerJMikullaRTrebinHR.IMD: a software package for molecular dynamics studies on parallel computers. Int J Mod Phys C1997; 8: 1131–1140.
86.
BuchholzMBungartzHJVrabecJ.Software design for a highly parallel molecular dynamics simulation framework in chemical engineering. J Computat Sci2011; 2: 124–129.
RothJGahlerFTrebinH.A molecular dynamics run with 5 180 116 000 particles. Int J Mod Phys C2000; 11: 317–322.
89.
NiethammerCBeckerSBernreutherM, et al. ls1 mardyn: the massively parallel molecular dynamics code for large systems. J Chem Theor Computat2014; 10: 4455–4464.
90.
XuHOsetskyYNStollerRE.Self-evolving atomistic kinetic Monte Carlo: fundamentals and applications. J Phys Condens Matt2012; 24: 375402.
91.
WasGS. Fundamentals of radiation materials science metals and alloys. Berlin: Springer, 2017.
92.
LiSWuBZhangY, et al. Massively scaling the metal microscopic damage simulation on Sunway Taihulight supercomputer. In: MalonyADShendeSWeilandMVenkataMG (eds) proceedings of the 47th international conference on parallel processing, Eugene, OR, 13–16 August 2018, p.47. New York: ACM.
93.
DomainCBecquartCVan DuysenJ.Kinetic Monte Carlo simulations of cascades in Fe alloys. MRS Online Proc Libr Arch2000; 650: R3.25.
94.
YoungWElcockE.Monte Carlo studies of vacancy migration in binary ordered alloys: I. Proc Phys Soc1966; 89: 735.
95.
VincentEBecquartCDomainC.Atomic kinetic Monte Carlo model based on ab initio data: Simulation of microstructural evolution under irradiation of dilute Fe–CuNiMnSi alloys. Nucl Instrum Meth Phys Res B Beam Interact Mater Atom2007; 255: 78–84.
96.
GillespieDT.A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. J Computat Phys1976; 22: 403–434.
97.
Martin-BragadoIRiveraAVallesG, et al. MMonCa: an object Kinetic Monte Carlo simulator for damage irradiation evolution and defect diffusion. Comput Phys Commun2013; 184: 2703–2710.
98.
ShimYAmarJG.Semirigorous synchronous sublattice algorithm for parallel kinetic Monte Carlo simulations of thin film growth. Phys Rev B2005; 71: 125432.
99.
EickSGGreenbergAGLubachevskyBD, et al. Synchronous relaxation for parallel simulations with applications to circuit-switched networks. ACM Trans Model Comput Simulat1993; 3: 287–314.
100.
MerrickMFichthornKA.Synchronous relaxation algorithm for parallel kinetic Monte Carlo simulations of thin film growth. Phys Rev E2007; 75: 011606.
101.
NandipatiGShimYAmarJG, et al. Parallel kinetic Monte Carlo simulations of Ag (111) island coarsening using a large database. J Phys Condens Matt2009; 21: 084214.
102.
MartínezEMarianJKalosMH, et al. Synchronous parallel kinetic Monte Carlo for continuum diffusion-reaction systems. J Computat Phys2008; 227: 3804–3823.
103.
BortzABKalosMHLebowitzJL.A new algorithm for Monte Carlo simulation of Ising spin systems. J Computat Phys1975; 17: 10–18.
104.
Martin-BragadoIAbujasJGalindoP, et al. Synchronous parallel kinetic Monte Carlo: Implementation and results for object and lattice approaches. Nucl Instrum Meth Phys Res B Beam Interact Mater Atom2015; 352: 27–30.
WuBLiSZhangY, et al. Hybrid-optimization strategy for the communication of large-scale Kinetic Monte Carlo simulation. Comput Phys Commun2017; 211: 113–123.
107.
BattaileCC.The kinetic Monte Carlo method: foundation, implementation, and application. Comput Meth Appl Mech Eng2008; 197: 3386–3398.
108.
LeetmaaMSkorodumovaNV.KMCLib: a general framework for lattice kinetic Monte Carlo (KMC) simulations. Comput Phys Commun2014; 185: 2340–2349.
109.
StollerREGolubovSIDomainC, et al. Mean field rate theory and object kinetic Monte Carlo: a comparison of kinetic models. J Nucl Mater2008; 382: 77–90.
110.
ChiapettoMMalerbaLBecquartC.Effect of Cr content on the nanostructural evolution of irradiated ferritic/martensitic alloys: an object kinetic Monte Carlo model. J Nucl Mater2015; 465: 326–336.
111.
ŽeníšekJSvobodaJKozeschnikE, et al. A combined cluster dynamics/kinetic Monte Carlo model for precipitate nucleation in interstitial/substitutional alloys. Model Simulat Mater Sci Eng2008; 16: 075005.
112.
BrimbalDFournierLBarbuA.Cluster dynamics modeling of the effect of high dose irradiation and helium on the microstructure of austenitic stainless steels. J Nucl Mater2016; 468: 124–139.
113.
DunnAYCapolungoLMartinezE, et al. Spatially resolved stochastic cluster dynamics for radiation damage evolution in nanostructured metals. J Nucl Mater2013; 443: 128–139.
114.
GolubovSOvcharenkoABarashevA, et al. Grouping method for the approximate solution of a kinetic equation describing the evolution of point-defect clusters. Philos Mag A2001; 81: 643–658.
115.
WolferWMansurLSpragueJ.Theory of swelling and irradiation creep. Technical report, Madison (USA). Department of Nuclear Engineering; Oak Ridge National Lab., TN (USA); Naval Research Lab., Washington, DC (USA), 1977.
116.
TurkinABakaiA.Formation of steady state size distribution of precipitates in alloys under cascade-producing irradiation. J Nucl Mater2006; 358: 10–25.
117.
MarianJBulatovVV.Stochastic cluster dynamics method for simulations of multispecies irradiation damage accumulation. J Nucl Mater2011; 415: 84–95.
118.
KiritaniM.Analysis of the clustering process of supersaturated lattice vacancies. J Phys Soc Jpn1973; 35: 95–107.
119.
OvcharenkoAGolubovSWooC, et al. GMIC++: grouping method in C++: an efficient method to solve large number of master equations. Comput Phys Commun2003; 152: 208–226.
120.
GolubovSStollerRZinkleS, et al. Kinetics of coarsening of helium bubbles during implantation and post-implantation annealing. J Nucl Mater2007; 361: 149–159.
121.
JourdanTStoltzGLegollF, et al. An accurate scheme to solve cluster dynamics equations using a Fokker–Planck approach. Comput Phys Commun2016; 207: 170–178.
122.
JourdanTBencteuxGAdjanorG.Efficient simulation of kinetics of radiation induced defects: a cluster dynamics approach. J Nucl Mater2014; 444: 298–313.
123.
RempelJYBawendiMGJensenKF.Insights into the kinetics of semiconductor nanocrystal nucleation and growth. J Am Chem Soc2009; 131: 4479–4489.
124.
HoangTLMarianJBulatovVV, et al. Computationally-efficient stochastic cluster dynamics method for modeling damage accumulation in irradiated materials. J Computat Phys2015; 300: 254–268.
125.
XuDWirthBDLiM, et al. Combining in situ transmission electron microscopy irradiation experiments with cluster dynamics modeling to study nanoscale defect agglomeration in structural metals. Acta Materialia2012; 60: 4286–4302.
126.
XuDCertainALee VoigtHJ, et al. Ballistic effects on the copper precipitation and re-dissolution kinetics in an ion irradiated and thermally annealed Fe–Cu alloy. J Chem Phys2016; 145: 104704.
127.
XuDVanCoeveringGWirthBD.Defect microstructural equivalence in molybdenum under different irradiation conditions at low temperatures and low doses. Computat Mater Sci2016; 114: 47–53.
XuDHuXWirthBD.A phase-cut method for multi-species kinetics: Sample application to nanoscale defect cluster evolution in alpha iron following helium ion implantation. Appl Phys Lett2013; 102: 011904.
131.
DunnAAgudo-MeridaLMartin-BragadoI, et al. A novel method for computing effective diffusivity: application to helium implanted -Fe thin films. J Nucl Mater2014; 448: 195–205.
132.
DunnACapolungoL.Simulating radiation damage accumulation in -Fe: a spatially resolved stochastic cluster dynamics approach. Computat Mater Sci2015; 102: 314–326.
133.
TerrierPAthènesMJourdanT, et al. Cluster dynamics modelling of materials: a new hybrid deterministic/stochastic coupling approach. J Computat Phys2017; 350: 280–295.
