skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Multi-scale modelling to relate beryllium surface temperature, deuterium concentration and erosion in fusion reactor environment

Journal Article · · Journal of Physics. D, Applied Physics
 [1];  [2];  [3];  [1]
  1. Univ. of Helsinki (Finland)
  2. Technical Univ. of Madrid (Spain). Inst. of Nuclear Fusion
  3. 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
Citation Metrics:
Cited by: 14 works
Citation information provided by
Web of Science

References (33)

Chapter 1: Overview and summary journal December 1999
Engineering challenges and development of the ITER Blanket System and Divertor journal October 2015
Plasma wall interactions in ITER journal April 2006
Study of physical and chemical assisted physical sputtering of beryllium in the JET ITER-like wall journal September 2014
Observations and modelling of diatomic molecular spectra from JET journal March 1998
Response of beryllium to deuterium plasma bombardment journal September 1998
The role of beryllium deuteride in plasma-beryllium interactions journal June 2009
Erosion of beryllium under high-flux plasma impact journal July 2013
Erosion of Be as plasma-facing material journal November 1997
Simulation of the plasma-wall interaction in a tokamak with the Monte Carlo code ERO-TEXTOR journal May 2000
Variables affecting simulated Be sputtering yields journal August 2013
Molecules can be sputtered also from pure metals: sputtering of beryllium hydride by fusion plasma–wall interactions journal June 2013
Chemical sputtering of Be due to D bombardment journal December 2009
Modelling of W–Be mixed material sputtering under D irradiation journal April 2014
Atomistic simulations of the effect of reactor-relevant parameters on be sputtering journal August 2015
Theoretical foundations of dynamical Monte Carlo simulations journal July 1991
Simulation of radiation damage in Fe alloys: an object kinetic Monte Carlo approach journal October 2004
Comparative study of radiation damage accumulation in Cu and Fe journal January 2000
A new algorithm for Monte Carlo simulation of Ising spin systems journal January 1975
Introduction to the Kinetic Monte Carlo Method book January 2007
MMonCa: An Object Kinetic Monte Carlo simulator for damage irradiation evolution and defect diffusion journal December 2013
Defects and transport processes in beryllium journal October 2011
Ab initio study of small vacancy complexes in beryllium journal February 2007
Quantum Modeling of Hydrogen Retention in Beryllium Bulk and Vacancies journal February 2010
Hydrogen in beryllium: Solubility, transport, and trapping journal April 2009
A climbing image nudged elastic band method for finding saddle points and minimum energy paths journal December 2000
Ab initiomolecular dynamics for liquid metals journal January 1993
Hydrogen isotope retention and recycling in fusion reactor plasma-facing components journal February 2002
Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points journal December 2000
SRIM – The stopping and range of ions in matter (2010)
  • Ziegler, James F.; Ziegler, M. D.; Biersack, J. P.
  • Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, Vol. 268, Issue 11-12 https://doi.org/10.1016/j.nimb.2010.02.091
journal June 2010
Interatomic potentials for the Be–C–H system journal October 2009
Molecular Dynamics Simulations of Deuterium Trapping and Re-emission in Tungsten Carbide journal January 2010
Erosion yields of deposited beryllium layers journal June 2009

Cited By (3)

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

Atomistic simulations of deuterium irradiation on iron-based alloys in future fusion reactors
Journal Article · Fri Oct 14 00:00:00 EDT 2016 · Nuclear Materials and Energy · OSTI ID:1376422

Kinetic Monte Carlo Simulation of Hydrogen Diffusion in Tungsten
Journal Article · Wed Jun 15 00:00:00 EDT 2016 · Transactions of the American Nuclear Society · OSTI ID:1376422

Integrated model predictions on the impact of substrate damage on gas dynamics during ITER burning-plasma operations
Journal Article · Wed Oct 13 00:00:00 EDT 2021 · Nuclear Fusion · OSTI ID:1376422