Abstract
Deuterium retention behavior in W-Fe-Ni alloy was investigated by the gas-phase driven permeation system and thermal desorption tests. The deuterium permeability, diffusion coefficient, solubility, and diffusion activation energy of deuterium in the W-Fe-Ni alloys were investigated. The thermal deuterium charging and thermal deuterium desorption experiments were conducted. Combined with the microstructure characteristics and numerical simulation, the deuterium retention behavior in W-Fe-Ni alloys was studied and the diffusion model of hydrogen isotope in W-Fe-Ni alloy was established to predict the deuterium retention in W-Fe-Ni alloys with different shapes. Compared with the results of thermal desorption tests, the amount of hydrogen isotope retention in W-Fe-Ni alloy can be accurately estimated by the multi-physics field numerical simulation.
Science Press
W-Fe-Ni alloy prepared by the liquid-phase sintering in hydrogen atmosphere is a two-phase alloy consisting of tung-sten and binder phases (nickel and iron
In this research, the transport of deuterium in W-Fe-Ni alloy was studied by the gas-phase driven permeation. The permeability, diffusion coefficient, and solubility of deuterium in the W-Fe-Ni alloy were obtained, and the hydrogen isotope diffusion model of W-Fe-Ni alloy was also established. The specimen was exposed to D2 with the pressure of 500 kPa at 773 K for 5 h to obtain the deuterium-saturated W-Fe-Ni alloy. The deuterium retention in W-Fe-Ni alloy was then studied by the thermal desorption spectroscopy. The thermal desorption characteristics of deuterium in the W-Fe-Ni alloy were obtained. The relationship between the deuterium release and desorption temperature was investigated. The retention mechanism of deuterium in W-Fe-Ni alloy was discussed. The experimental deuterium desorption amount was compared with the data obtained from the simulated deuterium charging process to evaluate the accuracy of the diffusion model.
The W-Fe-Ni alloy with the dimension of Φ12 mm×0.5 mm (Institute of Materials, China Academy of Engineering Physics) was mirror-polished by the 400#~2000# silicon carbide abrasive paper. The thickness of specimen after polishing was 0.456 and 0.442 mm. Then, the polished W-Fe-Ni alloy was ultrasonically rinsed in distilled water and dried by the argon gas for the gas permeation tests and thermal desorption tests. The specimens were placed in a quartz tube preheated to 1273 K at heating rate of 10 K/min with the vacuum degree >5.0×1
The permeation experiments were conducted by the gas- phase driven permeation system, as shown in

Fig.1 Schematic diagram of gas-phase permeation system
P1-molecular pump Number 1~9-high pressure valve (20 MPa)
;P2-oil pump BED-LaNi5 deuterium storage bed
;VG1-precision pressure gauge (1 MPa) VCR-penetration tooling
;VG2-precision pressure gauge (1 MPa) VG3-vacuum gauge (500 kPa)
The experiment process was as follows. (1) Air tightness test. After the installation was completed, the helium gas with the pressure>101.3 kPa was measured in the upstream VCR equipment. The helium gas was measured in the downstream vessel by QMS. If the helium ion signal is at the same level as that of the helium background, the VCR equipment has a good sealing performance. (2) Degassing. The heating furnace was used to heat the VCR equipment to the initial temperature (873 K) of the penetration experiment, and the heating belt was used to heat the downstream pipeline to reduce H2 and H2O as much as possible. The downstream vacuum was reduced to 1×1
The deuterium charging experiment was conducted in the high-pressure gas-solid reaction experiment system (Science and Technology on Surface Physics and Chemistry Laboratory), and the schematic diagram of the system is shown in

Fig.2 Schematic diagram of high-pressure gas-solid reaction experimental system
1#~10#-high pressure valve (20 MPa) D-thin film capacitance gauge (1333 Pa) A-mechanical dry pump E-metal reactor B-precision pressure gauge (0.1 MPa) F-precision pressure gauge (6 MPa) C-standard tank (300 mL) G-ZrCo deuterium storage bed
The process of hot charging experiment was as follows. (1) Air tightness test. Firstly, the annealed tungsten alloy was placed in the stainless steel reactor, and the corresponding valve was closed. Then the vacuum was pumped to <5 Pa, and the air tightness of the system was tested. (2) Pipeline gas cleaning. D2 was stored in the ZrCo storage bed, heated, and then entered into the standard vessel. The pipe system was cleaned three times by D2 of 1×1
After the deuterium charging experiment was completed, the specimen was placed for 5 d, and then the gas-phase thermal desorption tests were conducted. The system was equipped with a QMS to detect and quantitatively analyze the released gas components after the specimen was heated. The specimen was placed in a vacuum quartz tube with the vacuum pressure of more than 1×1
The data processing was conducted according to
(1) |
where C is the concentration of solid solution hydrogen, J is the diffusion flux (the number of substances diffused through the section of unit area per unit time), and D is the diffusion coefficient. Combined with Sievert's law, the steady-state permeation flux of hydrogen in the material can be obtained, as follows:
(2) |
where S is the solubility of hydrogen in the material; d is the thickness of the material; and are the pressures of hydrogen at the upstream and downstream ends of the system equipment, respectively. Under the experimental conditions of ph>>pl, the permeability (ϕ) formula of hydrogen in the material is ϕ=DS, and it can be transformed into
(3) |
Before the permeation experiment, the ion signal strength of QMS was calibrated by the standard deuterium leak and the relationship between the deuterium ion intensity (I) and leakage rate (L) can be obtained, as follows:
(4) |
According to the definition of leakage rate, L can be expressed as the variation of pressure in unit time (dt) in the chamber with a certain volume (V):
(5) |
where p is the pressure; n is the reaction order; R is the ideal gas constant; T is the permeability temperature. Combined with the definition of permeability flux J=1/A·dn/dt,
(6) |
where A is the permeability area of the specimen. Therefore, the calculation formula of permeability can be obtained, as follows:
(7) |
The diffusion coefficient of deuterium (D) in the material can be determined by the time-lag metho
(8) |
where τ corresponds to the time when 63% of the steady-state permeability flux is reached according to the permeability flux curv
The solubility S, diffusion coefficient D, and permeability ϕ of hydrogen isotopes in the structural materials in a certain temperature range can be transformed into the Arrhenius formula type, as follows:
(9) |
(10) |
where K represents S, D, and ϕ indicators; A0 is the pre-exponential factor of S, D, and ϕ indicators, and it can be divided into S0, D0, and ϕ0, respectively; E is activation energy of S, D, and ϕ indicators, and it can be divided into ES, ED, and Eϕ, respectively.
The permeability, diffusion coefficient, and solubility of deuterium in the W-Fe-Ni alloy were obtained through the penetration experiments. In order to measure the retention of hydrogen isotopes in the specimens, the diffusion model of hydrogen isotopes in tungsten-based alloy was established by the parameters obtained from the permeation experiment, and the desorption parameters were obtained from the deuterium-thermal desorption experiment with thermal filling to verify the rationality of the model. This numerical simulation model provides technical support to evaluate the hydrogen isotope retention in W-Fe-Ni alloy under different working conditions.
The W-Fe-Ni disc specimen with the dimension of Φ12 mm×0.5 mm possesses the axisymmetric properties. During the process of thermal charging and standing, the hydrogen isotope gas was adsorbed on the surface and desorbed from the surface, respectively. Without the consideration of the effect of material defects, the hydrogen isotope transport process in the specimen is mainly dominated by the diffusion. Therefore, the model can be established, as follows:
(11) |
where C is the concentration of hydrogen isotope atoms in the specimen, and DT is the atomic diffusion coefficient in the material.
The process of hydrogen isotope gas dissociation and atomic recombination mainly occurs on the material surface, which can be expressed by Sievert's law. In the process of thermal charging, the process of molecular dissociation and atomic recombination is considered to be in the dynamic equilibrium. Thus,
(12) |
where S is the atomic solubility in the material and P is the thermal charging pressure. During the experiment, the pressure fluctuation in the reactor was ignored due to the weak hydrogen absorption capacity of the specimen and the basically unchanged thermal charging pressure.
The partial pressure of hydrogen in the atmosphere is low during the specimen standing process, resulting in the neglection of the recombination on surface.
Parameter | Value | Ref |
---|---|---|
Diameter, Ф/mm | 12 | - |
Thickness, d/mm | 0.5 | - |
Room temperature, T1/K | 293 | - |
Experiment temperature, T2/K | 773 | - |
Diffusion, D/ |
3.7×1 | - |
Solubility, S/mol‧ |
2.0×1 | - |
Recombination coefficient, R/ |
3.2×1 |
[ |
Note: k is the Boltzmann constant.

Fig.3 XRD patterns of W-Fe-Ni alloy before annealing and referenced alloys

Fig.4 XRD patterns of W-Fe-Ni alloy after annealing and referenced alloys

Fig.5 SEM morphologies of W-Fe-Ni alloy before annealing

Fig.6 SEM morphologies of W-Fe-Ni alloy after annealing at 1273 K
As shown in

Fig.7 SEM backscattered electron morphology (a) and element distributions of W (b), Ni (c), and Fe (d) in W-Fe-Ni alloy after annealing at 1273 K
The element content in different regions of the alloy was determined by the energy dispersive spectroscopy (EDS), as shown in

Fig.8 SEM morphology (a) and EDS spectra of points marked in Fig.8a (b~e) of W-Fe-Ni alloy after annealing at 1273 K: (b) point 2; (c) point 6; (d) point 3; (e) point 1
The deuterium ion current signals detected by QMS in the deuterium permeation process of tungsten alloys at different temperatures (673~873 K) are shown in

Fig.9 Reltionship between ion flow intensity and time
The penetration parameters of permeability (ϕ), diffusion coefficient (D), and solubility (S) of W-Fe-Ni alloys are shown in
(13) |
(14) |
(15) |

Fig.10 Relationship between deuterium permeability ϕ and temperature T

Fig.11 Relationship between deuterium diffusion coefficient D and temperature T

Fig.12 Relationship between deuterium solubility S and temperature T
For comparison, the penetration parameters of permeability (ϕ), diffusion coefficient (D), and solubility (S) of pure tungsten are also shown in
(16) |
(17) |
(18) |
The related data are summarized in
The variation of deuterium in W-Fe-Ni alloy was simulated by the finite element method (FEM) with thermal charging at 773 K and 500 kPa for 5 h and releasing at room temperature for 5 d, as shown in

Fig.13 Variation of deuterium concentration and retention deuterium content of W-Fe-Ni alloys during thermal charging process (a) and releasing process (b)

Fig.14 Deuterium concentration distribution during releasing process
The deuterium thermal desorption tests were conducted on the deuterium specimens after the thermally charged deuterium experiment and standing for 5 d. A high desorption peak appears at 620.15 K, and the peak shape is consistent with the Gauss distribution, as shown in

Fig.15 Relationship between ion current and temperature during thermal desorption
The thermal desorption energy of deuterium can be obtained by
(19) |
where r is the thermal desorption rate, v is the desorption rate constant, n is the reaction order, θ is the coverage degree, and Ed is the thermal desorption activation energy. The thermal desorption activation energy Ed of deuterium in W-Fe-Ni specimen is calculated to be 0.35 eV.
In order to verify the accuracy of the model based on FEM, the actual deuterium desorption amount measured in the ther-mal desorption experiment was compared with the simulated results. Firstly, the signal intensity of deuterium particle flow under different pressures was measured by the standard leak holes, and the relationship between ion current and the deuterium leakage rate Q is obtained by
Q=0.547I-0.06 | (20) |
The deuterium amount can be deduced from the measured deuterium ion current and the conversion relationship between signal intensity and leakage rate.
The calculation results are shown in
Experiment | Simulation |
---|---|
1.572 | 8.009 |
1) W-Fe-Ni alloy prepared by the liquid-phase sintering in hydrogen atmosphere is a two-phase alloy consisting of W and binder phase (Fe-Ni). During the annealing at 1273 K and cooling in furnace, W extends to the binder phase and a needle-like phase, which mainly consists of W, is precipitated. The layered phase consisting of binder phase, W, C, and O appears around the W particles. The presence of carbon impurity has an important effect on the formation of these phases.
2) The solubility constant, diffusion activation energy, and activation energy for the permeation of deuterium in W-Fe-Ni alloy are lower than those of pure tungsten, because the Fe-Ni alloy in the binder phase has a higher deuterium permeability, which plays a leading role in the deuterium permeation process of the alloy.
3) The saturated retention of deuterium in W-Fe-Ni alloy is about 2.50 mol/
References
Şahin Y. Journal of Powder Technology[J], 2014, 2014: 764 306 [Baidu Scholar]
German R M, Suri P, Park S J. Journal of Materials Science[J], 2009, 44(1): 1 [Baidu Scholar]
Pang Qianlie. China Tungsten Industry[J], 2000, 15(4): 34 (in Chinese) [Baidu Scholar]
Fan Jinglian, Huang Baiyun, Wang Denglong et al. Rare Metal Materials and Engineering[J], 2001, 30(6): 401 (in Chinese) [Baidu Scholar]
Donten M, Cesiulis H, Stojek Z. Electrochimica Acta[J], 2000, 45(20): 3389 [Baidu Scholar]
German R M, Bose A. 120th Annu Meet Miner Met Mater Soc[C]. New Orleans: TMS, 1991: 53 [Baidu Scholar]
Jiao Z H, Kang R K, Dong Z G et al. International Journal of Refractory Metals and Hard Materials[J], 2019, 80: 114 [Baidu Scholar]
Deng S H, Yuan T C, Li R D et al. Powder Technology[J], 2017, 310: 264 [Baidu Scholar]
Wang Peng, Liu Kai. Heat Treatment of Metals[J], 2015, 40(5): 10 (in Chinese) [Baidu Scholar]
Symons D M. Metallurgical and Materials Transactions A[J], 1997, 28(3): 655 [Baidu Scholar]
Powell G L. Analytical Chemistry[J], 1972, 44(14): 2357 [Baidu Scholar]
Zhou Guoan, Lin Guobiao, Lai Heyi et al. Heat Treatment of Metals[J], 1995, 20(5): 13 (in Chinese) [Baidu Scholar]
Zhu Guiseng, Liu Mingcheng, Lin Zhongjie et al. Metallurgica Sinca[J], 1981, 17(1): 39 (in Chinese) [Baidu Scholar]
Wu Aihua. Thesis for Master[D]. Changsha: Central South University, 2004 (in Chinese) [Baidu Scholar]
Wang Bo, Liu Lingbo, Xiang Xin et al. Journal of Nuclear Materials[J], 2016, 470: 30 [Baidu Scholar]
Liu Feng. Thesis for Doctorate[D]. Beijing: University of Chinese Academy of Sciences, 2015 (in Chinese) [Baidu Scholar]
Frauenfelder R. The Journal of Chemical Physics[J], 1968, [Baidu Scholar]
48(9): 3955 [Baidu Scholar]
Hu A, Hassanein A. Journal of Nuclear Materials[J], 2014, [Baidu Scholar]
446(1-3): 56 [Baidu Scholar]
Lei Bingqiang, Zhu Guisen. Acta Metallurgica Sinca[J], 1987, 23(3): 182 (in Chinese) [Baidu Scholar]
Ebihara K I, Kaburaki H. ISIJ International[J], 2012, 52(2): 181 [Baidu Scholar]