skip to main content


This content will become publicly available on June 4, 2019

Title: Expanding the capability of reaction-diffusion codes using pseudo traps and temperature partitioning: Applied to hydrogen uptake and release from tungsten

Simulating the implantation and thermal desorption evolution in a reaction-diffusion model requires solving a set of coupled differential equations that describe the trapping and release of atomic species in Plasma Facing Materials (PFMs). These fundamental equations are well outlined by the Tritium Migration Analysis Program (TMAP) which can model systems with no more than three active traps per atomic species. To overcome this limitation, we have developed a Pseudo Trap and Temperature Partition (PTTP) scheme allowing us to lump multiple inactive traps into one pseudo trap, simplifying the system of equations to be solved. For all temperatures, we show the trapping of atoms from solute is exactly accounted for when using a pseudo trap. However, a single effective pseudo trap energy can not well replicate the release from multiple traps, each with its own detrapping energy. However, atoms held in a high energy trap will remain trapped at relatively low temperatures, and thus there is a temperature range in which release from high energy traps is effectively inactive. By partitioning the temperature range into segments, a pseudo trap can be defined for each segment to account for multiple high energy traps that are actively trapping but are effectively not releasingmore » atoms. With increasing temperature, as in controlled thermal desorption, the lowest energy trap is nearly emptied and can be removed from the set of coupled equations, while the next higher energy trap becomes an actively releasing trap. Each segment is thus calculated sequentially, with the last time step of a given segment solution being used as an initial input for the next segment as only the pseudo and actively releasing traps are modeled. This PTTP scheme is then applied to experimental thermal desorption data for tungsten (W) samples damaged with heavy ions, which display six distinct release peaks during thermal desorption. Without modifying the TMAP7 source code the PTTP scheme is shown to successfully model the D retention in all six traps. In conclusion, we demonstrate the full reconstruction from the plasma implantation phase through the controlled thermal desorption phase with detrapping energies near 0.9, 1.1, 1.4, 1.7, 1.9 and 2.1 eV for a W sample damaged at room temperature.« less
 [1] ;  [1] ;  [2] ;  [1] ;  [1] ;  [3]
  1. Univ. of California, San Diego, CA (United States). Center for Energy Research
  2. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
  3. Univ. of California, San Diego, CA (United States). Center for Energy Research; Univ. of California, San Diego, CA (United States). Dept. of Mechanical and Aerospace Engineering (MAE)
Publication Date:
Report Number(s):
Journal ID: ISSN 0022-3115
Grant/Contract Number:
AC52-06NA25396; FG02-07ER54912; SC0001999
Accepted Manuscript
Journal Name:
Journal of Nuclear Materials
Additional Journal Information:
Journal Volume: 508; Journal ID: ISSN 0022-3115
Research Org:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org:
Country of Publication:
United States
OSTI Identifier: