Multi-scale modelling to relate beryllium surface temperature, deuterium concentration and erosion in fusion reactor environment
- Univ. of Helsinki (Finland)
- Technical Univ. of Madrid (Spain). Inst. of Nuclear Fusion
- Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Beryllium (Be) has been chosen as the plasma-facing material for the main wall of ITER, the next generation fusion reactor. Identifying the key parameters that determine Be erosion under reactor relevant conditions is vital to predict the ITER plasma-facing component lifetime and viability. To date, a certain prediction of Be erosion, focusing on the effect of two such parameters, surface temperature and D surface content, has not been achieved. In this paper, we develop the first multi-scale KMC-MD modeling approach for Be to provide a more accurate database for its erosion, as well as investigating parameters that affect erosion. First, we calculate the complex relationship between surface temperature and D concentration precisely by simulating the time evolution of the system using an object kinetic Monte Carlo (OKMC) technique. These simulations provide a D surface concentration profile for any surface temperature and incoming D energy. We then describe how this profile can be implemented as a starting configuration in molecular dynamics (MD) simulations. We finally use MD simulations to investigate the effect of temperature (300–800 K) and impact energy (10–200 eV) on the erosion of Be due to D plasma irradiations. The results reveal a strong dependency of the D surface content on temperature. Increasing the surface temperature leads to a lower D concentration at the surface, because of the tendency of D atoms to avoid being accommodated in a vacancy, and de-trapping from impurity sites diffuse fast toward bulk. At the next step, total and molecular Be erosion yields due to D irradiations are analyzed using MD simulations. The results show a strong dependency of erosion yields on surface temperature and incoming ion energy. The total Be erosion yield increases with temperature for impact energies up to 100 eV. However, increasing temperature and impact energy results in a lower fraction of Be atoms being sputtered as BeD molecules due to the lower D surface concentrations at higher temperatures. Finally, these findings correlate well with different experiments performed at JET and PISCES-B devices.
- Research Organization:
- Univ. of Helsinki (Finland); Technical Univ. of Madrid (Spain); Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
- Sponsoring Organization:
- USDOE; European Commission (EC)
- Grant/Contract Number:
- AC05-00OR22725; 633053
- OSTI ID:
- 1376422
- Journal Information:
- Journal of Physics. D, Applied Physics, Vol. 50, Issue 20; ISSN 0022-3727
- Publisher:
- IOP PublishingCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Analytical bond order potential for simulations of BeO 1D and 2D nanostructures and plasma-surface interactions
|
journal | March 2018 |
Molecular dynamics simulation of beryllium oxide irradiated by deuterium ions: sputtering and reflection
|
journal | March 2019 |
ERO modeling and sensitivity analysis of locally enhanced beryllium erosion by magnetically connected antennas
|
journal | December 2017 |
Similar Records
Kinetic Monte Carlo Simulation of Hydrogen Diffusion in Tungsten
Integrated model predictions on the impact of substrate damage on gas dynamics during ITER burning-plasma operations