Modeling hydrogen isotope behavior in fusion plasma-facing components

In this work, we focus on understanding hydrogen isotope retention in plasma-facing materials in fusion devices. Three common simulation methods are usually used to study this problem that includes Monte Carlo, molecular dynamics, and numerical/analytical methods. A system of partial differential equations describing deuterium behavior in tungsten under various conditions is solved numerically to explain recent data compared to other methods. The developed model of hydrogen retention in metals includes classic, intercrystalline and trapped-induced Gorsky effects. The bombardment and depth profile of 200 eV deuterium in single crystal tungsten are simulated and compared with recent work. The total deuterium retention at various temperatures and fluences are also calculated and compared with available data. The results are in reasonable agreement with data and therefore, this model can be used to estimate deuterium inventory and recovery in future fusion devices.