DOE PAGES 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

Abstract

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 surfacemore » 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.« less

Authors:
 [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)
Publication Date:
Research Org.:
Univ. of Helsinki (Finland); Technical Univ. of Madrid (Spain); Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE; European Commission (EC)
OSTI Identifier:
1376422
Grant/Contract Number:  
AC05-00OR22725; 633053
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Physics. D, Applied Physics
Additional Journal Information:
Journal Volume: 50; Journal Issue: 20; Journal ID: ISSN 0022-3727
Publisher:
IOP Publishing
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; 36 MATERIALS SCIENCE

Citation Formats

Safi, E., Valles, G., Lasa, A., and Nordlund, K. Multi-scale modelling to relate beryllium surface temperature, deuterium concentration and erosion in fusion reactor environment. United States: N. p., 2017. Web. doi:10.1088/1361-6463/aa6967.
Safi, E., Valles, G., Lasa, A., & Nordlund, K. Multi-scale modelling to relate beryllium surface temperature, deuterium concentration and erosion in fusion reactor environment. United States. https://doi.org/10.1088/1361-6463/aa6967
Safi, E., Valles, G., Lasa, A., and Nordlund, K. Mon . "Multi-scale modelling to relate beryllium surface temperature, deuterium concentration and erosion in fusion reactor environment". United States. https://doi.org/10.1088/1361-6463/aa6967. https://www.osti.gov/servlets/purl/1376422.
@article{osti_1376422,
title = {Multi-scale modelling to relate beryllium surface temperature, deuterium concentration and erosion in fusion reactor environment},
author = {Safi, E. and Valles, G. and Lasa, A. and Nordlund, K.},
abstractNote = {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.},
doi = {10.1088/1361-6463/aa6967},
journal = {Journal of Physics. D, Applied Physics},
number = 20,
volume = 50,
place = {United States},
year = {Mon Mar 27 00:00:00 EDT 2017},
month = {Mon Mar 27 00:00:00 EDT 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 14 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Chapter 1: Overview and summary
journal, December 1999

  • Editors, ITER Physics Basis; Co-Chairs, ITER Physics Expert Group Chairs an; Unit, ITER Joint Central Team and Physics
  • Nuclear Fusion, Vol. 39, Issue 12
  • DOI: 10.1088/0029-5515/39/12/301

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

  • Duxbury, G.; Stamp, M. F.; Summers, H. P.
  • Plasma Physics and Controlled Fusion, Vol. 40, Issue 3
  • DOI: 10.1088/0741-3335/40/3/002

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

  • Fichthorn, Kristen A.; Weinberg, W. H.
  • The Journal of Chemical Physics, Vol. 95, Issue 2
  • DOI: 10.1063/1.461138

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

  • Martin-Bragado, Ignacio; Rivera, Antonio; Valles, Gonzalo
  • Computer Physics Communications, Vol. 184, Issue 12
  • DOI: 10.1016/j.cpc.2013.07.011

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

  • Allouche, A.; Oberkofler, M.; Reinelt, M.
  • The Journal of Physical Chemistry C, Vol. 114, Issue 8
  • DOI: 10.1021/jp910806j

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

  • Henkelman, Graeme; Uberuaga, Blas P.; Jónsson, Hannes
  • The Journal of Chemical Physics, Vol. 113, Issue 22, p. 9901-9904
  • DOI: 10.1063/1.1329672

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

  • Henkelman, Graeme; Jónsson, Hannes
  • The Journal of Chemical Physics, Vol. 113, Issue 22
  • DOI: 10.1063/1.1323224

SRIM – The stopping and range of ions in matter (2010)
journal, June 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
  • DOI: 10.1016/j.nimb.2010.02.091

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

  • Vörtler, Katharina; Nordlund, Kai
  • The Journal of Physical Chemistry C, Vol. 114, Issue 12
  • DOI: 10.1021/jp9054473

Erosion yields of deposited beryllium layers
journal, June 2009


Works referencing / citing this record:

Analytical bond order potential for simulations of BeO 1D and 2D nanostructures and plasma-surface interactions
journal, March 2018

  • Byggmästar, J.; Hodille, E. A.; Ferro, Y.
  • Journal of Physics: Condensed Matter, Vol. 30, Issue 13
  • DOI: 10.1088/1361-648x/aaafb3

Molecular dynamics simulation of beryllium oxide irradiated by deuterium ions: sputtering and reflection
journal, March 2019

  • Hodille, E. A.; Byggmästar, J.; Safi, E.
  • Journal of Physics: Condensed Matter, Vol. 31, Issue 18
  • DOI: 10.1088/1361-648x/ab04d7

ERO modeling and sensitivity analysis of locally enhanced beryllium erosion by magnetically connected antennas
journal, December 2017