Next Article in Journal
Fundamental Mechanisms for Irradiation-Hardening and Embrittlement: A Review
Next Article in Special Issue
Correction: Ono, S., et al. A Heterothermic Kinetic Model of Hydrogen Absorption in Metals with Subsurface Transport. Metals, 2019, 9(10), 1131
Previous Article in Journal
On the Development Concept for a New 718-Type Superalloy with Improved Temperature Capability
Previous Article in Special Issue
Efficient Synthesis of Alkali Borohydrides from Mechanochemical Reduction of Borates Using Magnesium–Aluminum-Based Waste
 
 
Correction published on 9 January 2020, see Metals 2020, 10(1), 108.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Heterothermic Kinetic Model of Hydrogen Absorption in Metals with Subsurface Transport

1
Department of Chemical Engineering, Kyoto University, Nishikyo, Kyoto 615-8510, Japan
2
Department of Chemical Engineering, University of Manchester, Oxford Rd, Manchester M13 9PL, UK
3
Graduate School of Advanced Integrated Studies in Human Survivability, Kyoto University, Sakyo, Kyoto 606-8306, Japan
*
Author to whom correspondence should be addressed.
Metals 2019, 9(10), 1131; https://doi.org/10.3390/met9101131
Submission received: 30 September 2019 / Revised: 18 October 2019 / Accepted: 20 October 2019 / Published: 22 October 2019
(This article belongs to the Special Issue Metals in Hydrogen Technology)

Abstract

:
A versatile numerical model for hydrogen absorption into metals was developed. Our model addresses the kinetics of surface adsorption, subsurface transport (which plays an important role for metals with active surfaces), and bulk diffusion processes. This model can allow researchers to perform simulations for various conditions, such as different material species, dimensions, structures, and operating conditions. Furthermore, our calculation scheme reflects the relationship between the temperature changes in metals caused by the heat of adsorption and absorption and the temperature-dependent kinetic parameters for simulation precision purposes. We demonstrated the numerical fitting of the experimental data for various Pd temperatures and sizes, with a single set of kinetic parameters, to determine the unknown kinetic constants. Using the developed model and determined kinetic constants, the transitions of the rate-determining steps on the conditions of metal-hydrogen systems are systematically analyzed. Conventionally, the temperature change of metals during hydrogen adsorption and absorption has not been a favorable phenomenon because it can cause errors when numerically estimating the hydrogen absorption rates. However, by our calculation scheme, the experimental data obtained under temperature changing conditions can be positively used for parameter fitting to efficiently and accurately determine the kinetic constants of the absorption process, even from a small number of experimental runs. In addition, we defined an effectiveness factor as the ratio between the actual absorption rate and the virtually calculated non-bulk-diffusion-controlled rate, to evaluate the quantitative influence of each individual transport process on the overall absorption process. Our model and calculation scheme may be a useful tool for designing high-performance hydrogen storage systems.

1. Introduction

Understanding the kinetics of the hydrogen transport processes on the surface of and inside of metals is essential for manufacturing practical solid-state hydrogen storage devices based on hydrogen-absorbing metals and alloys [1,2,3]. Modeling the kinetics of the hydrogen absorption processes to numerically simulate the absorption rates for different conditions would be useful for efficiently designing and optimizing high-performance hydrogen storage systems.
Figure 1 illustrates the basic concept of hydrogen absorption by metals [4,5]. First, gaseous hydrogen molecules dissociate into atoms at the metal surface and then the atoms are adsorbed on the surface sites. The hydrogen atoms are transported from the surface sites to the subsurface sites, and then to the bulk region where they diffuse. A region known as subsurface has been identified and defined in the literature [4,5,6]. The subsurface is the region located at the transition from the surface area to the bulk region of a material. The subsurface typically consists of one or two atomic monolayers immediately below the surface and is a different environment from both the surface and the bulk region. The potential energy profile of metals for hydrogen absorption is schematically presented in Figure 2. A difference in the activation energies of the subsurface and bulk region can be noticed.
Research into the adsorption and desorption of hydrogen at metal surfaces has made significant progress, and the corresponding physical mechanisms for these processes have been established [7,8,9,10,11]. However, the behavior of hydrogen at the subsurface sites is still poorly understood [12]. Moreover, the atomic-level hydrogen absorption mechanism has not been established yet, and scientists have been debating between the penetration mechanisms, quantum tunneling, or a multi-atomic concerted process [12].
For metals featuring high surface activity for the dissociative adsorption of hydrogen molecules, such as Pd [12,13,14,15], the transport of hydrogen across the subsurface region often becomes the rate-determining process, in contrast to metals featuring lower surface activity, such as Mg [16,17,18]. A significant discrepancy has been frequently observed between the atomic coverage at the surface and that at the subsurface for high-surface-activity metals featuring low activation energy or potential barrier for the dissociative surface adsorption of hydrogen molecules [5,6]. Therefore, the hydrogen transport across the subsurface should be taken into account when developing numerical models that would be used for highly accurate simulations of the net absorption rates of high-surface-activity metals. In addition, temperature changes in metals during hydrogen absorption, owing to the heat of adsorption and absorption, could cause errors in the numerical simulations used to estimate the hydrogen absorption rates because many kinetic parameters are temperature-sensitive. Therefore, in this study, we developed a versatile kinetic model of hydrogen absorption into metals, which accounted for the subsurface transport process and temperature changes.

2. Theory and Calculation Methods

The basic concept of our calculation model is schematically depicted in Figure 1 and Figure 2. First, gaseous hydrogen molecules dissociate into hydrogen atoms at the metal surface, then the atoms are adsorbed on the surface sites. The hydrogen atoms migrate from the surface sites to the subsurface sites and then to the bulk region where they diffuse. In the model, we accounted for the hydrogen fluxes for the surface adsorption and desorption processes, the subsurface transport between the surface and the bulk region, and the diffusion in the bulk region. The mass balances for the fractional coverages (occupancy) at the surface and subsurface (θ and θss, respectively) and the hydrogen content X [mol-H/mol-Pd] could be expressed as follows:
N s d θ d t = J a d s J d e s J s b + J b s ,
N s d θ s s d t = J s b J b s J d i f ,
N b δ d X d t = J s b J b s ,
where Jabs and Jdes are the hydrogen fluxes for the surface adsorption and desorption processes, respectively; Jsb and Jbs are the hydrogen fluxes for the inward and outward subsurface transport, respectively; Jdif is the hydrogen flux for the diffusion in the bulk region; Ns (2.8 × 10−5 mol m−2 [19]) is the number of atomic sites per unit of surface area; Nb (1.1 × 105 mol m−3 [19]) is the number of atomic sites per unit of volume; and δ is the thickness of the bulk region. The Jabs, Jdes [19], Jsb, Jbs [20], and Jdif values could be calculated as follows [19,20]:
J a d s = 2 S ( θ ) P H 2 2 π M H 2 R T ,
J d e s = 2 k d 0 N s 2 θ H H exp ( 2 E d e s R T ) ,
J s b = k s b 0 N s θ ( 1 θ s s ) exp ( E s b R T ) ,
J b s = k b s 0 N s θ s s ( 1 θ ) exp ( E b s R T ) ,
J d i f = D H 0 C ( z , t ) z exp ( E d i f R T ) ,
where S is the capturing coefficient; PH2 and MH2 are the partial pressure and molar mass of hydrogen, respectively; R and T are the ideal gas constant and absolute temperature, respectively; kd0 (4.8 × 1017 m2 mol−1 s−1 [19]), ksb0, and kbs0 are the frequency factors for the surface desorption and the inward and outward subsurface transport of hydrogen atoms, respectively; Edes (4.2 × 104 J mol−1 [19]), Esb (5.6 × 104 J mol−1 [20]), Ebs (2.2 × 104 J mol−1 [20]), and Edif (2.2 × 104 J mol−1 [20]) are the activation energy values for the desorption, inward and outward subsurface transport, and diffusion of hydrogen in the bulk region, respectively; DH0 (2.9 × 10−7 m2 s−1 [20]) is the pre-factor of diffusion coefficient of hydrogen; and C is the local concentration of hydrogen in the bulk region. Two adjacent vacant surface atomic sites are required for the dissociative surface adsorption of hydrogen molecules. Similarly, two adjacent occupied surface sites are required for the recombinative surface desorption of hydrogen atoms. Therefore, we rigorously accounted for the locations of the adsorption and desorption sites on the surface, as follows [19,21]:
θ O O = 1 θ θ O H 2 ,
θ H H = θ θ O H 2 ,
where θOO and θHH are the probabilities that two adjacent sites are both vacant or occupied, respectively, and θOH is the probability that one of two adjacent sites is vacant. Assuming that θOO, θHH, and θOH reach chemical equilibrium, the following equation has been proposed [19,21]:
4 θ O O θ H H θ O H 2 = exp ( E e q R T ) ,
where Eeq (2.1 × 103 J mol−1 [19]) is the energy separation between the site-occupation states. From Equations (9)–(11), θOO and θHH can be expressed as:
θ O O = 1 θ 2 θ ( 1 θ ) 1 + 1 4 θ ( 1 θ ) { 1 exp ( E e q R T ) } ,
θ H H = θ 2 θ ( 1 θ ) 1 + 1 4 θ ( 1 θ ) { 1 exp ( E e q R T ) } .
The capturing coefficient, S(θ), which takes into account the hopping of the physisorbed hydrogen molecules at the surface, can be expressed as follows [19,21]:
S ( θ ) = S 0 1 + K ( 1 θ O O 1 ) ,
where S0 is the capturing coefficient for θ = 0, which we set to be 1.0 (dimensionless) [21] in this study, and K (5.0 × 10−2 (dimensionless) [19]) is a constant related to the surface adsorption and desorption. The value of S0 has little influence on the calculation results in this study because of the subsurface-transport-control or bulk-diffusion-control nature. When fitting our calculations to the experimental data, we assumed that the adsorption and desorption at the Pd surface were always in equilibrium [19,22]. Based on the above-described model, we calculated the time evolution of the hydrogen concentration profile in Pd. We considered ksb0 and kbs0 to be the fitting kinetic parameters for the transport of hydrogen atoms between the Pd surface and subsurface, and we determined them using numerical fitting.

3. Experimental Methods

3 cm × 3 cm flat Pd plates of various thicknesses (0.2, 0.5, 1, and 2 mm) were used in this study. First, we annealed the Pd samples at 1000 °C for 10 h in a furnace under a nitrogen atmosphere to thoroughly crystallize the material and to degrease the surface for good experimental reproducibility. Then each Pd sample was placed inside a stainless steel vacuum chamber with the capacity of 1.55 × 103 cm3. The chamber was evacuated (~10−8 Torr) and then hydrogen gas was promptly injected into the chamber to reach the pressure of 760 Torr. Afterward, the pressure in the chamber gradually decreased owing to the absorption of hydrogen into Pd. The chamber pressure was monitored in time to determine the hydrogen content in Pd. The temperature of Pd was also monitored using a thermocouple. For some experimental runs, the Pd sample was actively heated using a heater to investigate the temperature dependence of the hydrogen absorption rate.

4. Results and Discussion

4.1. Hydrogen Absorption Rate at Room Temperature

4.1.1. Experimental Data and Numerical Fitting

Figure 3 presents the time evolution of the hydrogen content in a 1 mm-thick Pd sample at the initial Pd temperature of 17 °C; Pd was not actively heated. The fitting result of our model calculation is also plotted in Figure 3. For the numerical fitting, we considered ksb0 and kbs0 to be the fitting parameters, and used literature-reported values for Esb and Ebs [20] and also for the other adsorption, desorption, and diffusion parameters [19,20]. Through numerical fitting, we determined the kinetic parameters ksb0 and kbs0 to be 5.7 × 1012 and 1.3 × 1012 s−1, respectively, for the transportation of atomic hydrogen between the Pd surface and subsurface. Because a single, monotonic curve had to be fit, multiple combinations of ksb0 and kbs0 values could be used to achieve favorable fits.

4.1.2. Case Study Simulation I: Effective Surface Area and Hydrogen Absorption Rate

After we obtained the kinetic parameters, we were able to perform simulations. Herein, we included some case studies to demonstrate the usefulness of our numerical model. First, we increased the hydrogen absorption rate by increasing the effective surface area of the metals using a roughening surface treatment, such as chemical wet etching. Figure 4 presents the simulation results for hydrogen absorption using 1 mm-thick Pd plates of various effective surface areas. As the effective surface area increased, the hydrogen absorption rate increased. However, the degree of increment decreased as the surface area continued to increase. This tendency could be attributed to the changes in the hydrogen flux. During the initial absorption stage, as the effective surface area increased, the hydrogen flux significantly increased. However, the flux gradually decreased as the bulk diffusion slowed owing to the decrease in the hydrogen-concentration slope particularly around the central region of Pd. From these simulation results, we concluded that to increase the hydrogen absorption rate, it would be effective to achieve an effective surface area several times larger than that of the smooth surface Pd. Moreover, the simulation results indicated that, to further increase the hydrogen absorption rate, it would be advantageous to decrease the influence of the bulk diffusion.

4.1.3. Case Study Simulation II: Plate Thickness and Hydrogen Absorption Rate

Next, we considered another strategy for increasing the hydrogen absorption rate: thinning the Pd plate to decrease the effect of bulk diffusion. Figure 5 depicts the numerical simulation results for the time evolutions of the hydrogen content in Pd plates of various thicknesses. In addition, Figure 6 presents the time evolution of the spatial hydrogen concentration profiles in Pd plates of various thicknesses. As the thickness of the Pd plate decreased, the hydrogen absorption rate increased. This could be attributed to the smaller influence of bulk diffusion for the thinner Pd plates owing to the sustainability of the hydrogen concentration gradient. The amount of hydrogen distributed in the bulk regions of the 500 μm and 1 mm Pd plates significantly varied from the surface to the center of the Pd plates owing to the influence of the bulk diffusion control on the overall hydrogen transport process relative to the influence of the surface adsorption and the subsurface transport. By contrast, hydrogen was uniformly distributed in the 100 μm Pd plate owing to the small influence of bulk diffusion. Therefore, we could shift the rate-determining process in the Pd-hydrogen system from bulk diffusion to surface-to-subsurface transport by thinning the Pd plates to 100 μm.

4.1.4. Case Study Simulation III: Particle Size and Hydrogen Absorption Rate

We considered spherical metal particles, which could be more common for practical hydrogen storage applications, in our model, and we investigated decreasing the particle size to increase the hydrogen absorption rate. For the calculations, we used the spherical form of the diffusion equation [23], instead of Equation (8). Figure 7 presents the results of the numerical simulation for the time evolution of the hydrogen content in spherical Pd particles of various diameters. Figure 8 presents the time evolution of the spatial hydrogen concentration profiles in Pd particles of various diameters. Similarly to the plate-shaped Pd, as the diameter of the Pd particles decreased, the hydrogen absorption rate increased. This tendency was again attributed to the smaller influence of the bulk diffusion control for smaller particles. The hydrogen absorption rate was higher for the Pd particles than for the Pd plates of the same dimensions. This was attributed to the total surface area of the particles being larger than that of the plates. Similar to the Pd plates, it was determined that decreasing the particle diameter to 100 μm would circumvent diffusion control for the kinetically efficient hydrogen absorption process, as observed in the spatially uniform hydrogen concentration profile in Figure 8.

4.2. Hydrogen Absorption Rate at Higher Temperature

4.2.1. Heterothermic Model Development

When the Pd temperature was higher, the initial hydrogen absorption rate was high, which caused the Pd temperature to significantly increase owing to exothermic adsorption and absorption processes, as described below. When adsorption and absorption occur at high temperature, temperature evolution in metals has to be accounted for in the numerical simulation model. To address the temperature change in metals, as we have previously conceptually proposed [24], we added a calculation routine to our kinetic model where the Pd temperature was monitored at each time step of the process, as reflected in the computation of the series of kinetic parameters illustrated in Figure 9. We updated the Pd and gas phase temperatures for each time step after using the mass balance to calculate status parameters, such as the surface coverage and hydrogen concentration for each spatial mesh. One might think that temperature evolution could be calculated using the adsorption and absorption enthalpies without utilizing the monitored temperature values. However, such calculations would require that the heat dissipation rates of the Pd samples for the heat diffusion, convection, and radiation components under non-steady-state conditions be determined. This would be a challenging approach and therefore we adopted the monitoring scheme in this study. In addition, for the experimental data reported hereafter, we have artificially modified the onset hydrogen content in such a way that the temporally saturated value equaled the literature-reported equilibrium value for Pd [25]; this should cover the initial fast uptake that was potentially missed while monitoring the pressure.

4.2.2. Calibration of Gas Phase Temperature

Once hydrogen absorption into Pd started, the gas phase (hydrogen, in this study) temperature would no longer be equal to or at equilibrium with the measured temperature of Pd. To precisely determine the hydrogen content in Pd from the measured chamber pressure, the spatially averaged hydrogen temperature was required. Therefore, as a preparatory experiment, the chamber pressure was measured when a 1 mm-thick Pd sample was heated. In advance of this preparatory experiment, hydrogen was absorbed by the Pd sample until equilibrium was reached. When the Pd temperature increased, a certain amount of hydrogen would be desorbed owing to the decrease in the equilibrium hydrogen content in Pd [19]. To account for this effect, we estimated the gas temperature, Tgas, using the initial pressure and temperature, P0 and T0, respectively, as:
T g a s ( t ) = n H , 0 n H , 0 + { X 0 X e q ( T P d ( t ) ) } n P d / 2 P ( t ) P 0 T 0 ,
where nH,0 is the initial number of moles of hydrogen in gas phase; X0 and Xeq are the initial and equilibrium [25] hydrogen contents in Pd, respectively; TPd is the temperature of Pd; nPd is the number of moles of Pd; and P is the temporally monitored pressure in the chamber. The time evolutions of the experimentally measured P and TPd, and estimated Tgas values are plotted in Figure 10. Then, the relationship between Tgas and TPd is presented in Figure 11. By using quadratic fitting, the following numerical relationship between Tgas and TPd was established to precisely determine the hydrogen content evolution and thus the hydrogen absorption rate:
T g a s [ ° C ] = 0.004366118 T P d 2 0.1198489 T P d + 25.48853 .

4.2.3. Experiments and Calculations at Various Temperatures

We performed hydrogen absorption experiments at various Pd temperatures to test the heterothermic kinetic model we developed. Figure 12 presents the time evolutions of the hydrogen content in 1 mm-thick Pd samples of various initial temperatures and also the fitting results from our model calculations. Using numerical fitting, we determined that ksb0 = 5.8 × 1012 s−1 and kbs0 = 7.6 × 1011 s−1. We employed this single set of values of the parameters ksb0 and kbs0 for all of the following calculations. The discrepancy for the initial stage of hydrogen absorption between the experimental and numerical results, particularly for the higher temperatures, can be partly attributed to the initial fast uptake that was potentially missed while monitoring the pressure.

4.2.4. Experiments and Calculations for Various Metal Sizes

We tested the validity of our model and the ksb0 and kbs0 fitting parameters reported in the previous section by comparing the experimental data with the calculation results for various metal sizes. Figure 13 presents the time evolutions of the hydrogen content in Pd samples of various plate thicknesses (Pd was not actively heated) and also the fitting results from our model calculations. The discrepancy for the initial stage of hydrogen absorption between the experimental and numerical results, particularly for the thinner Pd samples, can be, again, partly attributed to the initial fast uptake that was potentially missed while monitoring the pressure. Nevertheless, the time constants for the absorption process were well reproduced for all the thicknesses using a single set of modelling parameters.
In Figure 12 and Figure 13, we have demonstrated the numerical calculations fitted a series of hydrogen absorption experimental data for different temperatures and Pd sizes using a single set of material parameters. In another point of view, we have been able to finely extract the unknown kinetic parameters using the series of fittings. For the conventional Arrhenius plot scheme, multiple experimental runs at various temperatures are required. Moreover, metal temperature is not stable during the experimental runs owing to the heat of adsorption and absorption, which is not favorable for the Arrhenius plots. Our model could provide a precise method for determining unknown parameters using the numerical fitting of the experimental data for the time evolution of the hydrogen absorption rate even from a small number of experimental runs by positively utilizing the change in temperature of the hydrogen-absorbing media; this could subsequently enable highly accurate numerical simulations. Using our calculation scheme, one can determine the temperature-dependent rate constants, including the reaction activation energy, utilizing data fitting. This would be similar to using the Arrhenius plot method but it would only require a single experimental run where the temperature would spontaneously vary owing to the exothermic adsorption and absorption processes. In addition, while we used Pd as the hydrogen-absorbing medium in this study, our numerical model could be applied to various types of metal/alloy-hydrogen systems.

4.3. Analysis of Rate-Determining Steps

Using numerical investigations in conjunction with experimental data, we have determined the kinetic parameters of the hydrogen absorption process. In this section, we use the results to analyze the rate-determining steps of the overall hydrogen transport process. We define the effectiveness factor, η, as the ratio of the actual hydrogen absorption rate to the virtually calculated non-bulk-diffusion-controlled hydrogen absorption rate. The η factor would enable us to quantitatively evaluate the influence of each transportation process and determine the rate-determining step for different metal species and structures (e.g., sizes and shapes). Figure 14 depicts the concept of the hydrogen absorption rate vs. the hydrogen content plots for the actual absorption rate, which accounted for the bulk diffusion process, and the virtually calculated absorption rate assuming sufficiently fast bulk diffusion. The hydrogen absorption rate decreased linearly as the hydrogen content increased throughout the entire hydrogen content range for the non-bulk-diffusion-control case. This occurred because in such a surface-adsorption-control or subsurface-transport-control regime, the surface or subsurface coverage becomes approximately proportional to the hydrogen content, while the adsorption or subsurface-transport rate is approximately proportional to the surface or subsurface vacancy fraction, respectively. By contrast, the relationship between the actual absorption rate and hydrogen content becomes linear in the high-hydrogen-content range because in such a bulk-diffusion-control regime the hydrogen concentration gradient approximately linearly decreases as the hydrogen content increases owing to the quasi-steady state of diffusion. Therefore, the ratio between the actual and virtual rates becomes approximately constant in the high-hydrogen-content range. For practical applications, the total hydrogen absorption duration is mainly dominated by the duration of the absorption process in the high-hydrogen-content range, where the absorption rate becomes lower and therefore the process requires more time. Therefore, it would be important to analyze the rate-determining steps, particularly in the high-hydrogen-content range. Consequently, we defined η as the ratio of the actual rate to the virtual, non-bulk-diffusion-controlled rate for the range when the hydrogen content exceeds 0.4, where η becomes constant.
Figure 15 presents the calculated η values of the Pd plates and spherical Pd particles of various thicknesses and diameters, respectively, under the constant hydrogen pressure values of 1 and 1000 atm. As the size of the Pd plates and particles decreased, η approached unity, which indicated that the process was subsurface-transport controlled. In general, a high-η state corresponds to a surface-adsorption or subsurface-transport controlled process, but in this study it was always subsurface-transport controlled owing to the high surface activity of Pd. When the plate thickness or particle diameter was smaller than 100 μm, the influence of bulk diffusion could be neglected, which was consistent with the results in Figure 5, Figure 6, Figure 7 and Figure 8. Conversely, when the Pd size increased, η decreased and the hydrogen transport kinetics of the H-Pd system became bulk-diffusion controlled. The η values of the spherical Pd particles were larger than those of the Pd plates of the same sizes because of the larger specific surface area (surface area-to-volume ratio) of the particulates. In addition, the η values were higher at higher hydrogen pressure levels. Nevertheless, the differences in the η values were not significant despite the 1000-fold increase in pressure. This was attributed to the surface adsorption rate, which is highly influenced by the hydrogen pressure, being sufficiently high compared with the subsurface-transport and bulk-diffusion rates, even at 1 atm in the high-hydrogen-content region, and therefore the effect of the hydrogen pressure on η was not significant. In this way, using η, one can quantitatively evaluate the influence of the bulk diffusion on the entire hydrogen transport process in metals and also identify the rate-determining step of the system for different conditions.

5. Conclusions

In this study, we developed a versatile kinetic numerical model for calculating the time-dependent rates of hydrogen absorption into metals. Our hydrogen transportation model particularly addressed the kinetics of the surface dissociative adsorption, subsurface transport (which significantly affects metals with active surfaces), and bulk diffusion. Moreover, this model enables us to perform flexible numerical simulations that could be universally used for various conditions, including different hydrogen absorbing material species, absorbing materials of different dimensions and structures, and, also, different operating conditions. Furthermore, our calculation scheme reflected the relationship between the temperature changes in metals, owing to the heat of adsorption and absorption, and the temperature-dependent kinetic parameters for highly precise numerical simulations. We demonstrated numerical calculations fitted a series of experimental data on hydrogen absorption for various Pd temperature ranges and sizes by a single set of the material’s kinetic parameters. In another point of view, we have been able to finely extract the unknown kinetic parameters using a series of fittings. At the present stage, our numerical fitting was not necessarily excellent for the whole range of conditions. The fitting accuracy may be improved by increasing the precision of finely time-resolved pressure monitoring in the hydrogen absorption experiments. Nevertheless, the time constants for the absorption process were well reproduced for various Pd temperatures and sizes using a single set of modelling parameters. Using the developed model and the determined kinetic constants, the transitions of the rate-determining steps on the conditions of metal-hydrogen systems were systematically analyzed in a series of case studies. We showed the effectiveness of increasing the effective surface area, by several times relative to the smooth surface, to increase the hydrogen absorption rate for Pd through numerical simulations. Our simulation results also indicated that, to further increase the hydrogen absorption rate, it would be advantageous to decrease the influence of the bulk diffusion. Moreover, we found that decreasing the Pd size to 100 m would circumvent diffusion control for the kinetically efficient hydrogen absorption process. Conventionally, the changes in the temperature of metals during hydrogen adsorption and absorption have not been a favorable phenomenon because it could cause errors in the numerically estimated hydrogen absorption rates. By contrast, the model and numerical scheme we proposed in this study could represent highly accurate methods for determining unknown kinetic parameters, such as temperature-dependent rate constants, including the reaction activation energy, by numerically fitting the time-evolution experimental data of the hydrogen absorption rate. This would be similar to using the Arrhenius plot method; however, our model and scheme relies on a single experimental run as the temperature spontaneously varies owing to the exothermic adsorption and absorption processes. One can thus positively utilize the change in temperature of the hydrogen-absorbing media, to subsequently enable highly accurate numerical simulations. While we used Pd as the hydrogen-absorbing medium in this study, our numerical model could be applied to various types of metal/alloy-hydrogen systems. A potential disadvantage of our scheme may be the requirement of accurate monitoring of the time-transient temperature of the hydrogen-absorbing media, even for the temporal region with drastic temperature changes. In addition, we introduced an effectiveness factor: the ratio between the actual absorption rate and the virtually calculated one neglecting bulk diffusion, as a parameter that would allow us to evaluate the quantitative influence of each individual transport process and determine the rate-determining steps depending on the conditions of the metal-hydrogen systems, such as metal species and structures (e.g., sizes and shapes). The kinetic model and calculation scheme presented could allow for high-precision parameter determination and numerical simulations, and thus could be a useful tool for efficiently designing high-performance hydrogen storage systems.

Author Contributions

Y.H., Y.K., and K.T. conceived the idea of the study. E.Y. provided the experimental chamber and the Pd material. S.O., T.U., Y.H., Y.K., E.Y., and K.T. set up the experimental system. S.O., T.U., Y.H., Y.K., and G.Y. performed the experiments. S.O., T.U., Y.H., Y.K., and K.T. set up the theory. S.O., T.U., Y.H., and Y.K. carried out the calculations. S.O., T.U., Y.H., Y.K., G.Y., and K.T. analyzed the data. K.T. supervised the study. S.O., T.U., Y.H., and K.T. composed the manuscript.

Funding

This research was funded by the Japan Society for the Promotion of Science, grant number 19K22081, the Thermal & Electric Energy Technology Foundation, and the Research Foundation for Opto-Science and Technology.

Acknowledgments

The authors would like to thank Hiroshi Sugiura of Powdec K.K., Shinji Kai of Tanaka Kikinzoku Kogyo K.K., and Takenori Naito of Kyoto University for their technical support and advice.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Schlapbach, L.; Zuttel, A. Hydrogen-storage materials for mobile applications. Nature 2001, 414, 353–358. [Google Scholar] [CrossRef] [PubMed]
  2. Ley, M.B.; Jepsen, L.H.; Lee, Y.S.; Cho, Y.W.; von Colbe, J.M.B.; Dornheim, M.; Rokni, M.; Jensen, J.O.; Sloth, M.; Filinchuk, Y.; et al. Complex hydrides for hydrogen storage—New perspectives. Mater. Today 2014, 17, 122–128. [Google Scholar] [CrossRef]
  3. Mohtadi, R.; Orimo, S. The renaissance of hydrides as energy materials. Nat. Rev. Mater. 2016, 2, 16091. [Google Scholar] [CrossRef]
  4. Behm, R.J.; Penka, V.; Cattania, M.G.; Christmann, K.; Ertl, G. Evidence for “subsurface” hydrogen on Pd(100): An intermediate between chemisorbed and dissolved species. J. Chem. Phys. 1983, 78, 7486–7490. [Google Scholar] [CrossRef]
  5. Wilde, M.; Fukutani, K. Penetration mechanisms of surface-adsorbed hydrogen atoms into bulk metals: Experiment and model. Phys. Rev. B 2008, 78, 115411. [Google Scholar] [CrossRef]
  6. Greeley, J.; Mavrikakis, M. Surface and subsurface hydrogen: Adsorption properties on transition metals and near-surface alloys. J. Phys. Chem. B 2005, 109, 3460–3471. [Google Scholar] [CrossRef]
  7. Christmann, K. Interaction of hydrogen with solid surfaces. Surf. Sci. Rep. 1988, 9, 1–163. [Google Scholar] [CrossRef]
  8. Dino, W.A.; Kasai, H.; Okiji, A. Orientational effects in dissociative adsorption/associative desorption dynamics of H2 (D2) on Cu and Pd. Prog. Surf. Sci. 2000, 63, 63–134. [Google Scholar] [CrossRef]
  9. Pundt, A.; Kirchheim, R. Hydrogen in metals: Microstructural aspects. Annu. Rev. Mater. Res. 2006, 36, 555–608. [Google Scholar] [CrossRef]
  10. Du, A.J.; Smith, S.C.; Yao, X.D.; Lu, G.Q. Hydrogen spillover mechanism on a Pd-doped Mg surface as revealed by ab initio density functional calculation. J. Am. Chem. Soc. 2007, 129, 10201–10204. [Google Scholar] [CrossRef]
  11. Firmino, T.; Marquardt, R.; Gatti, F.; Dong, W. Diffusion rates for hydrogen on Pd(111) from molecular quantum dynamics calculations. J. Phys. Chem. Lett. 2014, 5, 4270–4274. [Google Scholar] [CrossRef] [PubMed]
  12. Fukutani, K.; Wilde, M.; Ogura, S. Nuclear dynamics and electronic effects of hydrogen on solid surfaces. Chem. Rec. 2017, 17, 233–249. [Google Scholar] [CrossRef] [PubMed]
  13. Lynch, J.F.; Flanagan, T.B. An investigation of the dynamic equilibrium between chemisorbed and absorbed hydrogen in the palladium/hydrogen system. J. Phys. Chem. 1973, 77, 2628–2634. [Google Scholar] [CrossRef]
  14. Conrad, H.; Ertl, G.; Latta, E.E. Adsorption of hydrogen on palladium single crystal surfaces. Surf. Sci. 1974, 41, 435–446. [Google Scholar] [CrossRef]
  15. Wagner, S.; Pundt, A. Quasi-thermodynamic model on hydride formation in palladium-hydrogen thin films: Impact of elastic and microstructural constraints. Int. J. Hydrogen Energy 2016, 41, 2727–2738. [Google Scholar] [CrossRef]
  16. Huot, J.; Liang, G.; Boily, S.; Van Neste, A.; Schulz, R. Structural study and hydrogen sorption kinetics of ball-milled magnesium hydride. J. Alloys Compd. 1999, 293, 495–500. [Google Scholar] [CrossRef]
  17. Pozzo, M.; Alfe, D. Hydrogen dissociation and diffusion on transition metal (= Ti, Zr, V, Fe, Ru, Co, Rh, Ni, Pd, Cu, Ag)-doped Mg(0001) surfaces. Int. J. Hydrogen Energy 2009, 34, 1922–1930. [Google Scholar] [CrossRef]
  18. Evard, E.; Gabis, I.; Yartys, V.A. Kinetics of hydrogen evolution from MgH2: Experimental studies, mechanism and modelling. Int. J. Hydrogen Energy 2010, 35, 9060–9069. [Google Scholar] [CrossRef]
  19. Ward, T.L.; Dao, T. Model of hydrogen permeation behavior in palladium membranes. J. Membr. Sci. 1999, 153, 211–231. [Google Scholar] [CrossRef]
  20. Bhargav, A.; Jackson, G.S.; Ciora, R.J.; Liu, P.T.K. Model development and validation of hydrogen transport through supported palladium membranes. J. Membr. Sci. 2010, 356, 123–132. [Google Scholar] [CrossRef]
  21. Behm, R.J.; Christmann, K.; Ertl, G. Adsorption of hydrogen on Pd(100). Surf. Sci. 1980, 99, 320–340. [Google Scholar] [CrossRef]
  22. Ali-Khan, I.; Dietz, K.J.; Waelbroeck, F.G.; Wienhold, P. The rate of hydrogen release out of clean metallic surfaces. J. Nucl. Mater. 1978, 76, 337–343. [Google Scholar] [CrossRef]
  23. Kitagawa, Y.; Tanabe, K. Development of a kinetic model of hydrogen absorption and desorption in magnesium and analysis of the rate-determining step. Chem. Phys. Lett. 2018, 699, 132–138. [Google Scholar] [CrossRef]
  24. Tanabe, K. Modeling of hydrogen/deuterium dynamics and heat generation on palladium nanoparticles for hydrogen storage and solid-state nuclear fusion. Heliyon 2016, 2, e00057. [Google Scholar] [CrossRef]
  25. Frieske, H.; Wicke, E. Magnetic susceptibility and equilibrium diagram of PdHn. Ber. Der Bunsen-Ges. Fur Phys. Chem. 1973, 77, 48–52. [Google Scholar] [CrossRef]
Figure 1. Conceptual schematic diagram of hydrogen absorption into metals. Jabs and Jdes are the hydrogen fluxes for the surface adsorption and desorption processes, respectively; Jsb and Jbs are the fluxes for the migration of hydrogen from the surface to the subsurface, and from the subsurface to the surface, respectively; and Jdif is the hydrogen flux for the diffusion in the bulk region.
Figure 1. Conceptual schematic diagram of hydrogen absorption into metals. Jabs and Jdes are the hydrogen fluxes for the surface adsorption and desorption processes, respectively; Jsb and Jbs are the fluxes for the migration of hydrogen from the surface to the subsurface, and from the subsurface to the surface, respectively; and Jdif is the hydrogen flux for the diffusion in the bulk region.
Metals 09 01131 g001
Figure 2. Schematic diagram of potential energy profile of hydrogen absorbing metals. Edes, Esb, Ebs, and Edif are the activation energy values for the desorption, migration of hydrogen from the surface to the subsurface and from the subsurface to the surface, and diffusion of hydrogen in the bulk region, respectively.
Figure 2. Schematic diagram of potential energy profile of hydrogen absorbing metals. Edes, Esb, Ebs, and Edif are the activation energy values for the desorption, migration of hydrogen from the surface to the subsurface and from the subsurface to the surface, and diffusion of hydrogen in the bulk region, respectively.
Metals 09 01131 g002
Figure 3. Experimental result and calculated fitted result of the time evolution of hydrogen content in Pd at 17 °C.
Figure 3. Experimental result and calculated fitted result of the time evolution of hydrogen content in Pd at 17 °C.
Metals 09 01131 g003
Figure 4. Simulated time evolution plots of hydrogen content absorbed in Pd plates of various effective surface areas, which were 1, 2, 5, 10, and 100 times larger than that of the smooth surface Pd. The inset depicts a magnified view for the first 30 min of the hydrogen absorption process. The accompanying schematic (top) is a conceptual drawing of the surface roughening process.
Figure 4. Simulated time evolution plots of hydrogen content absorbed in Pd plates of various effective surface areas, which were 1, 2, 5, 10, and 100 times larger than that of the smooth surface Pd. The inset depicts a magnified view for the first 30 min of the hydrogen absorption process. The accompanying schematic (top) is a conceptual drawing of the surface roughening process.
Metals 09 01131 g004
Figure 5. Simulated time evolutions of hydrogen content in Pd plates of various thicknesses.
Figure 5. Simulated time evolutions of hydrogen content in Pd plates of various thicknesses.
Metals 09 01131 g005
Figure 6. Spatial hydrogen concentration profiles in Pd plates at 5, 10, 30, 50, 100, and 300 min after the beginning of hydrogen absorption process for Pd plates of different thicknesses.
Figure 6. Spatial hydrogen concentration profiles in Pd plates at 5, 10, 30, 50, 100, and 300 min after the beginning of hydrogen absorption process for Pd plates of different thicknesses.
Metals 09 01131 g006
Figure 7. Simulated time evolution of hydrogen content in spherical Pd particles of various diameters.
Figure 7. Simulated time evolution of hydrogen content in spherical Pd particles of various diameters.
Metals 09 01131 g007
Figure 8. Spatial hydrogen concentration profiles in spherical Pd particles at 5, 10, 30, 50, and 100 min after the beginning of hydrogen absorption process for Pd particles of different diameters.
Figure 8. Spatial hydrogen concentration profiles in spherical Pd particles at 5, 10, 30, 50, and 100 min after the beginning of hydrogen absorption process for Pd particles of different diameters.
Metals 09 01131 g008
Figure 9. Summary of calculation routine accounting for temperature changes in metals.
Figure 9. Summary of calculation routine accounting for temperature changes in metals.
Metals 09 01131 g009
Figure 10. Time evolution of measured Pd temperature, gas phase pressure, and numerically determined gas phase temperature.
Figure 10. Time evolution of measured Pd temperature, gas phase pressure, and numerically determined gas phase temperature.
Metals 09 01131 g010
Figure 11. Relationship between Pd and gas phase temperatures.
Figure 11. Relationship between Pd and gas phase temperatures.
Metals 09 01131 g011
Figure 12. Experimental results (blank circles) and numerical fitting (solid lines) of time evolution of hydrogen content in Pd at different temperatures (a: low, b: middle, and c: high temperature ranges). The fitting parameters were ksb0 = 5.8 × 1012 s−1 and kbs0 = 7.6 × 1011 s−1 for all the calculations.
Figure 12. Experimental results (blank circles) and numerical fitting (solid lines) of time evolution of hydrogen content in Pd at different temperatures (a: low, b: middle, and c: high temperature ranges). The fitting parameters were ksb0 = 5.8 × 1012 s−1 and kbs0 = 7.6 × 1011 s−1 for all the calculations.
Metals 09 01131 g012
Figure 13. Experimental results (blank circles) and calculation fittings (solid lines) of time evolution of hydrogen content in Pd plates of various thicknesses. The single set of parameters ksb0 = 5.8 × 1012 s−1 and kbs0 = 7.6 × 1011 s−1 was used for all the calculations.
Figure 13. Experimental results (blank circles) and calculation fittings (solid lines) of time evolution of hydrogen content in Pd plates of various thicknesses. The single set of parameters ksb0 = 5.8 × 1012 s−1 and kbs0 = 7.6 × 1011 s−1 was used for all the calculations.
Metals 09 01131 g013
Figure 14. Conceptual schematic for definition of effectiveness factor η; η = robs/rideal, where robs and rideal are the actually observed hydrogen absorption rate that accounts for bulk diffusion and the virtual absorption rate that was calculated neglecting the bulk diffusion process. robs and rideal were plotted against the hydrogen content.
Figure 14. Conceptual schematic for definition of effectiveness factor η; η = robs/rideal, where robs and rideal are the actually observed hydrogen absorption rate that accounts for bulk diffusion and the virtual absorption rate that was calculated neglecting the bulk diffusion process. robs and rideal were plotted against the hydrogen content.
Metals 09 01131 g014
Figure 15. Calculated η values for Pd plates and spherical particles of various thicknesses and diameters, respectively, under constant hydrogen pressure of 1 and 1000 atm at 25 °C.
Figure 15. Calculated η values for Pd plates and spherical particles of various thicknesses and diameters, respectively, under constant hydrogen pressure of 1 and 1000 atm at 25 °C.
Metals 09 01131 g015

Share and Cite

MDPI and ACS Style

Ono, S.; Uchikoshi, T.; Hayashi, Y.; Kitagawa, Y.; Yeh, G.; Yamaguchi, E.; Tanabe, K. A Heterothermic Kinetic Model of Hydrogen Absorption in Metals with Subsurface Transport. Metals 2019, 9, 1131. https://doi.org/10.3390/met9101131

AMA Style

Ono S, Uchikoshi T, Hayashi Y, Kitagawa Y, Yeh G, Yamaguchi E, Tanabe K. A Heterothermic Kinetic Model of Hydrogen Absorption in Metals with Subsurface Transport. Metals. 2019; 9(10):1131. https://doi.org/10.3390/met9101131

Chicago/Turabian Style

Ono, Shunsuke, Takeru Uchikoshi, Yusuke Hayashi, Yuta Kitagawa, George Yeh, Eiichi Yamaguchi, and Katsuaki Tanabe. 2019. "A Heterothermic Kinetic Model of Hydrogen Absorption in Metals with Subsurface Transport" Metals 9, no. 10: 1131. https://doi.org/10.3390/met9101131

APA Style

Ono, S., Uchikoshi, T., Hayashi, Y., Kitagawa, Y., Yeh, G., Yamaguchi, E., & Tanabe, K. (2019). A Heterothermic Kinetic Model of Hydrogen Absorption in Metals with Subsurface Transport. Metals, 9(10), 1131. https://doi.org/10.3390/met9101131

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop