• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Effect of ternary transition metal sulfide FeNi2S4 on hydrogen storage performance of MgH2

    2023-11-18 01:12:36YaokunFuLuZhangaYuanLiaSanyangGuoHanYuWenfengWangKailiangRenWeiZhangShuminHana
    Journal of Magnesium and Alloys 2023年8期

    Yaokun Fu, Lu Zhanga,,?, Yuan Lia,, Sanyang Guo, Han Yu, Wenfeng Wang,Kailiang Ren, Wei Zhang, Shumin Hana,,?

    a State Key Laboratory of Metastable Materials Science and Technology, Yanshan University, Qinhuangdao 066004, China

    b Hebei Key Laboratory of Applied Chemistry, School of Environmental and Chemical Engineering, Yanshan University, Qinhuangdao 066004, China

    Received 6 July 2021; received in revised form 5 November 2021; accepted 29 November 2021

    Available online 27 January 2022

    Abstract Hydrogen storage is a key link in hydrogen economy, where solid-state hydrogen storage is considered as the most promising approach because it can meet the requirement of high density and safety. Thereinto, magnesium-based materials (MgH2) are currently deemed as an attractive candidate due to the potentially high hydrogen storage density (7.6 wt%), however, the stable thermodynamics and slow kinetics limit the practical application. In this study, we design a ternary transition metal sulfide FeNi2S4 with a hollow balloon structure as a catalyst of MgH2 to address the above issues by constructing a MgH2/Mg2NiH4-MgS/Fe system. Notably, the dehydrogenation/hydrogenation of MgH2 has been significantly improved due to the synergistic catalysis of active species of Mg2Ni/Mg2NiH4, MgS and Fe originated from the MgH2-FeNi2S4 composite. The hydrogen absorption capacity of the MgH2-FeNi2S4 composite reaches to 4.02 wt% at 373 K for 1 h, a sharp contrast to the milled-MgH2 (0.67 wt%). In terms of dehydrogenation process, the initial dehydrogenation temperature of the composite is 80 K lower than that of the milled-MgH2, and the dehydrogenation activation energy decreases by 95.7 kJ·mol–1 compared with the milled-MgH2 (161.2 kJ·mol–1). This method provides a new strategy for improving the dehydrogenation/hydrogenation performance of the MgH2 material.

    Keywords: Hydrogen storage materials; Magnesium hydride; FeNi2S4; Hydrogen storage kinetics; Catalysts.

    1. Introduction

    With the rapid development of society,energy consumption has been intensifying and is expected to reach a peak in 2035[1–3]. For the sustainable use of energy in the future, hydrogen energy has been regarded as one of the most promising alternative energy sources of traditional fossil fuels because of its high energy density (142 MJ·kg–1). However, achieving safe, efficient and economical storage of hydrogen has been a challenge for hydrogen application [4–7]. Magnesium based hydrogen storage materials (MgH2) have been considered to be one of the most promising solid hydrogen storage materials due to the high hydrogen storage capacity (~7.6 wt%),good reversibility and low cost, but the high thermal stability of Mg-H bond leads to the high dehydrogenation temperature (above 623 K) and sluggish reaction kinetics of MgH2,which hinders its practical application as a hydrogen storage medium [8–11].

    In recent years, researchers have made in-depth studies to improve the thermodynamic and kinetic properties of MgH2.The specific methods include alloying, nano-confining and catalysts addition,etc. [12–15]. Among them, the addition of highly active catalysts has been proved an effectively approach to reduce the activation energy of hydrogen desorption reaction,thus improving the dehydrogenation kinetics[16,17].The explored application of additives with forward catalytic activity include modified carbon materials [18–20], transition metals [21–23] and their oxides [24–26], sulfides [27–29],halides [30–32], and alloy compounds [33–36]. Thereinto,transition metals can effectively weaken the Mg-H bond because hydrogen atoms tend to form covalent bonds with transition metals rather than ionic bonds with magnesium atoms[37,38], where the covalent bond between 3d transition metal elements and H atoms is relatively weaker [39]. Therefore,transition metals have been deemed as the main catalytic active substances and widely used to improve the hydrogen absorption and desorption performance of Mg/MgH2hydrogen storage systems [40,41].

    One of the effective transition metal catalyst for Mg/MgH2system is nickel (Ni). The catalytic effect of Ni on hydrogen storage performance of MgH2has been extensively studied.For instance, Chen et al. reported that Ni nanofibers exhibited enhanced catalytic effect on the dehydrogenation of MgH2,which rapidly released 7.02 wt% H2within 11 mins at 598 K[42].Lan’s group showed that through Ni@rGO nanocomposites modification, the apparent activation energy of dehydrogenation of MgH2was reduced by 60.4 kJ·mol-1[43].Except for metal Ni, its compounds, especially nickel sulfides, are also beneficial to the hydrogen absorption and desorption performance of Mg/MgH2system. Wang et al. reported that the hydrogen absorption kinetics of Mg/MgH2system had been significantly improvedviathe addition of NiS2owing to the in-situ formed catalysts of MgS and Mg2NiH4phases during heating process [15]. In Xie’s work, they found that NiS reduced the activation energy of hydrogen desorption of MgH2by 64.71 kJ·mol-1[44]. In addition, transition metal iron (Fe)is also emerged as a highly effective catalyst for hydrogen absorption and desorption of Mg/MgH2system [45,46]. As reported by the deepgoing study of Zhang’s group, both Fe3S4and FeS2could remarkably increase the re/dehydrogenation kinetics of Mg/MgH2system and lower its initial dehydrogenation temperature [47,48]. In details, the initial dehydrogenation temperature of MgH2was reduced by 90 K and 102 K, respectively, after Fe3S4and FeS2were added, and Fe3S4and FeS2changed the rate control steps and promoted H diffusion of MgH2during hydrogen absorption. Despite the progress made to these transition metal catalysts, there still lack suitable catalysts that facilitate the fast hydrogen absorption and desorption for Mg/MgH2system, particularly for the operation at low temperatures.

    Recently, studies have pointed out that ternary transition metal sulfides,e.g.FeNi2S4, NiCo2S4and FeCo2S4, are with excellent catalytic activity owing to the synergistic effect between different 3d transition metals, which have been confirmed in the fields of supercapacitors, oxygen and hydrogen evolution reactions [49–52]. Along this line, it is valuable to reveal whether the ternary transition metal sulfides promote hydrogen absorption and desorption properties of Mg/MgH2system. Therefore, herein, we reported the effect of ternary transition metal sulfides of FeNi2S4on hydrogen storage performance of MgH2for the first time, where the hydrogen absorption capacity of MgH2-FeNi2S4composite at 373 K ups to 4.02 wt%, about six times than that of pure MgH2(0.67 wt%), and its initial dehydrogenation temperature is 80 K lower than that of the MgH2. Furthermore, the phase transformation process of the in-situ formation of MgS, Mg2NiH4and Fe species by the reaction of FeNi2S4with Mg/MgH2and the catalytic mechanism of these catalytic active substances are discussed in detail.

    2. Experimental section

    2.1. Synthesis of FeNi2S4

    All materials and reagents in the experiment were used as received without any further purification. FeNi2S4was synthesized by a facile one-step solvothermal method.First, 2 mmol FeSO4·7H2O (Aladdin, ≥99.0%), 4 mmol Ni(CH3COO)2·4H2O (Aladdin, 99.0%) and 16 mmol Lcysteine (Aladdin, 99.0%) were dissolved in a 100 mL mixed solution of deionized water and ethylene glycol (Tianjin Kaitong Chemical Reagent Co. LTD, 99.0%) with a volume ratio of 1:1 and magnetically stirred rapidly for 30 min. The solution was then transferred to a 200 mL Teflon-lined stainless steel autoclave and reacted at 433 K for 8 h. Finally,the black precipitate FeNi2S4was washed several times with deionized water followed by ethyl alcohol and dried at 333 K for 12 h in vacuum.

    2.2. Preparation of MgH2-FeNi2S4 composite

    The synthesized FeNi2S4was directly mixed with commercial Mg powder (Aladdin, 99.0%) at a mass ratio of 1:10 and hydrogenated at 673 K at 4 MPa for 20 h to form a MgH2-FeNi2S4composite as in our previous study[53].After partial Mg was hydrogenated, high energy ball milling was carried out with a ball-to-powder ratio of 40:1 at 400 rpm for 5 h,during which intermittent method is used (milling for 15 min and then pause for 15 min) to avoid friction overheating. The hydrogenation and ball milling processes were alternately repeated three times to ensure the complete conversion of Mg to MgH2. For comparison, the Mg powder without FeNi2S4was treated under the same condition and denoted as-milled MgH2. All samples’ transfer operations were carried out in a glove box filled with high purity argon (H2O and O2content less than 0.1 ppm). Schematic diagram of synthesis process model of the MgH2-FeNi2S4composite is as shown in Fig.S1.

    2.3. Materials characterization

    Fig. 1. Microstructure and composition of the FeNi2S4. (a) XRD pattern, (b) SEM images and corresponding elemental mappings of (c) Ni, (d) Fe, and (e)S, (f) TEM image, XPS spectra of (g) Fe 2p, (h) Ni 2p, and (i) S 2p.

    X-ray diffraction (XRD) was performed on a Rigaku Smaratlab X-ray diffraction powder diffractometer using Cu Kαradiation (λ= 1.5418 ?A) in a range of 10°–90° (2θ) at a scanning rate of 5° min–1. X-ray photoelectron spectroscopy(XPS) was conducted on a Thermo ScientificTMK-AlphaTM+spectrometer equipped with a monochromatic Al KαX-ray source (1486.6 eV) operating at 100 W. All peaks would be calibrated with C1s peak binding energy at 284.8 eV for adventitious carbon. The microscopic morphology and element distribution of the samples were studied through scanning electron microscope (SEM, SUPRA 55) with energy spectrometer (EDS) and transmission electron microscopy (TEM,JEM-2010). The particle size distribution was measured from the corresponding SEM images by the Nano Measurer 1.2 software. Inductively coupled plasma-atomic emission spectrometer (ICP) was performed on an iCAP-6300. Differential scanning calorimetry (DSC, STA449C) was used to measure the thermal decomposition of the samples at different heating rates of 5, 10, 15 and 20 K·min–1from room temperature to 823 K under high purity argon condition.

    The hydrogen absorption and desorption properties of the MgH2-FeNi2S4and as-milled MgH2samples were measured by Sievert-type pressure-composition-temperature (PCT, Beijing General Research Institute of Nonferrous Metals) volumetric apparatus. The operating pressure of hydrogen absorption and desorption should be higher than the equilibrium pressure of hydrogen absorption of Mg and lower than the equilibrium pressure of hydrogen desorption of MgH2, respectively. Therefore, the hydrogen absorption and desorption pressures were set as 3 MPa and 0.01 MPa. The temperatureprogrammed desorption (TPD) test was conducted from room temperature to 823 K at a heating rate of 5 K·min–1.

    3. Results and discussion

    3.1. Characteristic of the synthesized FeNi2S4

    Fig. 1(a) shows that the XRD pattern of the synthesized FeNi2S4is a cubic structure with a space group Fd-3 m(JCPDS No. 47–1740) in pure phase. The FeNi2S4is corroborated as a hollow balloon structure by the SEM studies(Fig. 1(b)), where damaged mouths of the hollow balloons are observed as marked in red circles. EDS results further demonstrate the evenly distributed Fe, Ni and S elements in the FeNi2S4(Fig.1(c-e)),and the ICP result reveals an atomic ratio of Fe and Ni of approximately 1:2 (Table S1), verifying that material’s formula of FeNi2S4. The same EDS elemental mapping was done on the other parts of the FeNi2S4, and all the elements were overlaid with their corresponding colors,where uniform distribution of the elements could be clearly odserved (Fig. S2). TEM image in Fig. 1(f) also confirms that FeNi2S4has a hollow balloon, of which the diameter is~2 μm, as attested in the particle size distribution test (Fig.S3).

    XPS was further performed to identify the elements of FeNi2S4. The survey spectrum of FeNi2S4proves the presence of Fe, Ni and S, and the appearance of oxygen may be related to the air adsorbed on surface or impurities in the compound containing oxygen groups (Fig. S4). In the Fe 2p spectrum (Fig. 1(g)), there are two main peaks at 711.3 eV and 723.6 eV, which can be attributed to Fe 2p3/2of Fe3+and Fe 2p1/2of Fe2+, respectively. For the Ni 2p spectrum(Fig. 1(h)), the two main spin-orbit peaks 2p3/2(852.0 eV,869.4 eV) and 2p1/2(853.5 eV, 869.4 eV) indicate the presence of nickel oxidation states of Ni2+and Ni3+, respectively.The peak centered at 858.5 eV and 877.0 eV correspond to two shakeup satellites peaks of Ni 2p [50]. The exchange of electrons between different valence bands may provide an environment for electron movement on the phase surface, leading to higher catalytic activity [29, 54]. Fig. 1(i) shows the S 2p spectrum, where the binding energies of 161.1 eV and 162.5 eV correspond to the S 2p3/2and S 2p1/2, and the satellite peak of 163.3 eV corresponds to the metal-sulfur bond in the FeNi2S4; besides, the weaker peak at 167.0 eV is due to the O-S bond. The XPS results further confirm the chemical formula of the ternary transition metal sulfide.

    Fig. 2. Microstructure and composition of the MgH2-FeNi2S4 composite after ball-milling, dehydrogenation and rehydrogenation. (a) Rietveld refinement patterns of the XRD patterns of MgH2-FeNi2S4 composite after ball-milling and hydrogenation combustion process (BM-HCS), rehydrogenation and dehydrogenation, where the experimental data, calculated results, Bragg diffraction positions, and the difference profile are marked with black dots, a red curve,blue vertical bars, and a green curve, respectively, and the blue vertical bars from top to bottom belong to the phases listed in the corresponding figures in a same order; SEM images of the (b) ball-milling, (c) dehydrogenation and (d) rehydrogenation samples; Particle size distributions of the (e) ball-milling, (f)dehydrogenation and (g) rehydrogenation samples.

    3.2. Phase composition and morphology characterization of MgH2-FeNi2S4 composites

    The FeNi2S4compound was composited with Mg powder by 10 wt%to synthesis a MgH2-10wt%FeNi2S4composite by ball-milling and hydrogenation combustion processes (BMHCS). The schematic diagram of the preparation process of the MgH2-FeNi2S4composite is shown in Fig. S1. During the hydrogenation and ball milling processes, the chemical bond of the FeNi2S4will break and react with Mg and H2,and then in-situ generates new phases of MgH2, Mg2NiH4,MgS and Fe, thereby establishing a hydrogen storage system of MgH2/Mg2NiH4-MgS/Fe. Fig. 2(a) compares Rietveld refinement of the XRD patterns of the MgH2-10wt%FeNi2S4composite obtained by the BM-HCS, and after dehydrogenation and rehydrogenation processes. Results show that the obtained MgH2-10wt%FeNi2S4composite after BM-HCS composes MgH2, Mg2NiH4, MgS, MgO, Mg and Fe phases in contents of 93.30 wt%, 1.22 wt%, 0.62 wt%, 2.06 wt%,1.84 wt% and 0.92 wt%, respectively. The emergence of the Mg2NiH4,MgS and Fe phases proves the formation of heterogeneous catalytic system. Moreover, no new compound forms after rehydrogenation compared with the original composite,however the contents of Mg2NiH4(89.89 wt%), Mg (4.17 wt%) and Fe (2.53 wt%) increased slightly while the contents of MgH2(1.54 wt%), MgS (0.31 wt%) and MgO (1.57 wt%) decreased after multiple hydrogenation. Noteably, Mg is found in the XRD patterns for both samples, with a content of 1.84 wt% and 4.17 wt% respectively, which is caused by the inherent nature of Mg powder and the short time of experimental data collection. The appearance of MgO is inevitable because the sample is briefly exposed to air during testing,while Mg is very easily corroded by oxygen and water vapor in the air.

    Based on the XRD characterization of hydrogen absorption, the possible phase transformation reaction of FeNi2S4with Mg during hydrogenation can be described as follows:

    Fig. 3. Hydrogen absorption performance of the MgH2-FeNi2S4 composite and ball-milled MgH2. (a) Isothermal absorption kinetics at different temperatures(3.0 MPa); (b) Hydrogen absorption capacity at different temperatures (3.0 MPa); (c) Isothermal absorption kinetics at 473 K (3.0 MPa); (d) Isothermal absorption kinetics at 373 K (3.0 MPa).

    For dehydrogenation (Fig. 2(a)), MgH2and Mg2NiH4transform into Mg and Mg2Ni respectively. The diffraction of MgS and Fe still exist, and their contents increase to 1.84 wt% and 2.09 wt%, respectively, which may be that the reaction in Eq. (1) has not stabilized yet. Fig. 2(b-d) shows the SEM images of the above samples, where we find that the particle size of the sample noticeably increases after dehydrogenation and rehydrogenation, which are 1.62 μm and 3.29 μm, respectively, but just 0.84 μm for the original composite (Fig. 2(e-g)). One reason is that in the process of hydrogen absorption, the entry of hydrogen leads to the lattice expansion of the in-situ generated Mg and Mg2Ni. Table S2 shows that the cell volumes of Mg and Fe increase significantly after rehydrogenation, and the expansion rates were 0.252% and 0.516%, respectively. With the increase hydrogen absorption cycles, the lattice continues to expand, and the grain spacing gradually increases, eventually leads to the increase of particle size in a certain degree.However,the most important reason is that the agglomeration effect of small particles causes the increase of particle size, which is inevitable in the reaction process.

    3.3. Hydrogen absorption and desorption kinetics

    Fig. 3 displays the isothermal absorption kinetics curves of the MgH2-FeNi2S4composite and the as-milled MgH2at different temperatures. As can be seen from Fig. 3(a), the hydrogen absorption capacity of the ball-milled MgH2is higher than that of the MgH2-FeNi2S4composite at high temperatures (573 K–623 K), because the theoretical capacity of the composite is reduced after the addition of FeNi2S4and a small amount of Mg has been converted to MgS during hydrogenation. However, with the decrease of operating temperature,the hydrogen absorption of the ball-milled MgH2drops off in a cliff-like manner, while the MgH2-FeNi2S4shows excellent low temperature hydrogen absorption performance(Fig.3(b)).The MgH2-FeNi2S4composite exhibits faster hydrogen absorption rate and higher hydrogen absorption capacity when the temperature is below 523 K. At 473 K, the hydrogen absorption capacity of the MgH2-FeNi2S4composite in 1 h reaches 5.8 wt%, about 3.1 times than that of the pure Mg(1.89 wt%), and it can rapidly absorb 4.7 wt% H2just within 60 s (Fig. 3(c)). Further, as the temperature reduces to 373 K,the composite still absorbs 4.02 wt%H2in 1 h,while the pure Mg can only absorb 0.67 wt% H2(Fig. 3(d)).

    The improvement on the hydrogen desorption performance of MgH2viacompositing with FeNi2S4is also reflected in the rate and capacity. The dehydrogenation rate of the composite is much faster than that of MgH2at 623 K (Fig. 4(a)),moreover, when the temperature is lowered to 573 K, the ball-milled MgH2almost no longer desorbs hydrogen, while the MgH2-FeNi2S4composite still releases 1.92 wt% H2(Fig. 4(b)). The initial hydrogen desorption temperature is another important indicator to measure the hydrogen desorption property of hydrogen storage materials. Fig. 4(c) shows that the initial hydrogen release temperature of the MgH2-FeNi2S4composite is ~540 K, which is ~80 K lower than that of the ground MgH2(~620 K). In addition, the desorption peak temperature of the composite is ~622 K (Fig. 4(d)), ~50 K lower than that of the ball-milled MgH2(672 K). The above results indicate that the FeNi2S4has outstanding catalytic activity on dehydrogenation of MgH2, because Mg2NiH4can dehydrogenate in advance of MgH2, and Fe and MgS catalysis can work synergistically on the dehydrogenation of MgH2[55,56].

    Fig. 4. Hydrogen desorption performance of the MgH2-FeNi2S4 composite and ball-milled MgH2. Isothermal desorption kinetics of the MgH2-FeNi2S4 composite and ball-milled MgH2 at (a) 623 K (0.01 MPa) and (b) 573 K (0.01 MPa); (c) TPD and (d) differential curves of the MgH2-FeNi2S4 composite and ball-milled MgH2.

    Table 1.Theoretical kinetic models and goodness of fitting for the sorption kinetic data of MgH2-FeNi2S4. α= reacted fraction, t = time, k = reaction rate constant[59].

    To clarify the catalytic mechanism of FeNi2S4on hydrogen absorption and desorption reactions of MgH2, the kinetic mechanism of hydrogenation and dehydrogenation of the composite are analyzed. Theoretical mechanical equations of contracting volume (CV) model and Johnson-Mehl-Avrami(JMA) model are used to fit the isothermal hydrogen absorption and desorption data (Table 1), where the Johnson-Mehl-Avrami (JMA) model is suitable for the case that the nucleation and growth of the new phase randomly start from the bulk phase, but nucleation begins at the surface and continues to grow from the surface to bulk in the contracting volume(CV) model [57].

    Fig.5. JMA model and CV model fitting curves for isothermal hydrogen absorption and desorption data of the MgH2-FeNi2S4 composite.Hydrogen absorption at (a) 598 K, (b) 473 K and (c) 373 K; (d) Hydrogen desorption at 598 K.

    The absorption and desorption kinetic curves of the MgH2-FeNi2S4composite is calculated based on the data within the capacity of 0–80 wt%[31].Notably,the rate-limiting step differs at different temperatures during hydrogenation (Fig. 5(ac)). At 598 K, the rate-limiting step is a three-dimensional growth and a diffusion-controlled process by the decreasing interface velocity (CV 3D div) (Fig. 5(a)). When the temperature drops, the corresponding theoretical kinetic model changes. At 473 K, the rate-limiting step can be described by a three-dimensional growth model with constant interface velocity (CV 3D civ) (Fig. 5(b)), while at 373 K, the best equation describing the hydriding rate of the composite becomes uncertain (Fig. 5(c)). As can be seen in Fig. 5(c), the fit results calculated using the CV 2D formula and the corresponding formula for CV 3D civ are extremely close, with R2of 0.99678 and 0.99676, respectively. Due to the stubborn thermodynamic barrier of MgH2, the H diffusion rate and nucleation rate of the hydrogen storage system are affected as the temperature decreases.The above results suggest that temperature is an important factor limiting the hydrogen absorption rate except for the catalyst,which is due to the inherently stable thermodynamic properties of MgH2. For dehydrogenation of the composite, the rate-limiting step at 598 K is consistent with the three-dimensional growth model of contracting volume at constant interfacial velocity (Fig. 5(d)), which demonstrates that the dehydrogenation reaction starts from the surface, and the reaction is controlled by the interfacial speed of the grain boundary transformation phase. In other words,due to the large particle size of the hydride, hydrogen ion diffusion is the rate-limiting step during dehydrogenation reaction. It is worth noting that the dehydrogenation reaction is initiated by the H-anion binding, namely the decomposition of the Mg-H bond, but the strong Mg-H bond has a high energy barrier, therefore it is necessary to study the weakening effect of FeNi2S4on the Mg-H bond [57,58].

    3.4. Dehydrogenation kinetics and thermodynamics

    Dehydrogenation activation energy (Ea) is an important index to evaluate the kinetic performance of hydrogen storage materials, which can be calculated according to parameters obtained from the DSC test based on Kissinger method [63].Fig.6(a)shows the DSC curves during dehydrogenation of the MgH2-FeNi2S4composite. Apparently, the endothermic peak temperature of the composite is much lower than that of the ball-milled MgH2(Fig. S5), which is consistent with the TPD results.Moreover,the DSC curves of the MgH2-FeNi2S4composite present two endothermic peaks, which correspond the hydrogen release processes of Mg2NiH4and MgH2, respectively. TheEavalue is calculated by the following Kissinger’s equation as Eq. (2):

    Fig. 6. Dehydrogenation activation energy and cycling performance of the MgH2-FeNi2S4 composite. (a) DSC dehydrogenation temperature profiles; (b)Kissinger plots compared with MgH2; (c) Kinetic curves of isothermal hydrogen absorption and desorption at 598 K for within 10 cycles; (d) XRD patterns after cycling 10 times.

    whereαis the heating rate,Tmis the peak temperature, and R is the gas constant. TheEavalue of the ball-milled MgH2is 161.2 kJ·mol–1(Fig.6(b)),but it decreases to 90.9 kJ·mol–1and 65.5 kJ·mol–1for the MgH2-FeNi2S4composite, respectively, for the two stages of hydrogen desorption as circled by green ellipses and marked as 1 and 2 in Fig.6(a).TPD curves should also show two platforms corresponding the two hydrogen desorption stages according to the reports in the literature[29],but it is not obvious in this study.The appearance of two endothermic peaks may be due to that the hydrogen desorption of Mg2NiH4and MgH2is simultaneous in the first stage,but Mg2NiH4is dominant. In the second stage, as the temperature rises, most of MgH2absorbs heat and releases hydrogen. Meanwhile, in order to further explore the improvement of thermodynamic properties of MgH2by adding FeNi2S4,the dehydrogenation enthalpy change of MgH2-FeNi2S4composite was measured by integrating DSC peaks. The dehydrogenation enthalpy change of MgH2-FeNi2S4composite is determined to be 1552 J·g–1H2, which is lower than the ballmilled MgH2, indicating that FeNi2S4has a positive effect on the improvement of the thermodynamic properties of MgH2.In addition, it is worth noting that compared with the hydrogen storage system doped with other transition metal sulfides in related literatures as shown in Table 2, theEavalue of the MgH2-FeNi2S4composite is at an upstream level. In conclusion, the addition of transition metal sulfide FeNi2S4reduces the energy barrier for the binding of hydrogen anions in thesystem, thus reducing the hydrogen desorption activation energy of MgH2to a large extent.

    Table 2.The apparent activation energy(Ea)values of several MgH2-catalyst systems.

    3.5. Cycling performance

    Fig. 7. Structure and Morphology of the MgH2-FeNi2S4 composite after hydrogenation and dehydrogenation. (a) TEM image, (b) SAED pattern, and (c)HRTEM image after hydrogenation; (d) TEM image, (e) SAED pattern, and (f) HRTEM image after dehydrogenation.

    In order to study the cycle stability of the MgH2-FeNi2S4composite,the isothermal hydrogenation and dehydrogenation within 10 cycles were carried out at 325 K and 3.0 MPa. As shown in Fig. 6(c), the capacity and rate of hydrogen absorption and desorption of the MgH2-FeNi2S4composite hardly decay within 10 cycles, indicating its good cycle stability.However, the hydrogen desorption capacity of the composite is 1.5 wt% lower than the hydrogen absorption capacity on average(Fig.S6).This is due to the inherent stable thermodynamic performance of MgH2. The XRD patterns of the composite after hydrogenation and dehydrogenation at the 10thcycle are shown in Fig. 6(d), where no significant change occurrs compared with the previous XRD. Since the temperature of the hydrogen absorption and desorption cycle is 598 K and the period is only 30 min, the hydrogen absorption and desorption processes are not fully completed, so Mg and MgH2are inevitably found in the diffractions after hydrogenation and dehydrogenation, respectively. SEM shows the morphology of the MgH2-FeNi2S4compositeviahydrogenation after cycling 10 times (Fig. S6). Notably, the particles slightly agglomerate due to Van der Waals attraction after cycling [65,66], which is consistent with the previous description (Fig. 2(c)). In terms of element distribution, MgS, Fe and Mg2NiH4are distributed uniformly in the MgH2system;thereinto, the continuous catalytic action of MgS and Fe is the important factor for maintaining the cyclic stability of the composite.

    3.6. Catalytic mechanism

    To further investigate the catalytic reaction mechanism of FeNi2S4in Mg/MgH2system, TEM study was carried out. Fig. 7(a) and (d) show the TEM images of the MgH2-FeNi2S4composite after hydrogenation and dehydrogenation.As shown in the corresponding SEAD patterns, the diffraction rings of MgS and Fe appear in both hydrogenation and dehydrogenation samples, and the MgH2and Mg2NiH4in the hydrogenation sample become Mg and Mg2Ni after the hydrogen is released (Fig. 7(b) and (e)). The lattice fringes corresponding to the (111) plane of MgH2and the (220)plane of Mg2NiH4with interplanar spacing of 0.2194 nm and 0.2289 nm, respectively, are also clearly observed from HRTEM (Fig. 7(c)). After dehydrogenation, the interplanar spacings of 0.2451 nm and 0.2247 nm of the sample corresponding to the Mg (111) plane and the Mg2Ni (200)plane are observed in Fig. 7(f). The crystal planes of MgS(200) and Fe (110) remain the same after hydrogen absorption and desorption, but there is a slight measurement deviation for the lattice spacing. These results as well as the corresponding SAED pattern are consistent with the previous XRD results, revealing the reaction process of FeNi2S4in Mg/MgH2system. More importantly, we can observe that the catalytic active of MgS and Fe species disperse among the main phase Mg/MgH2from the HRTEM images. Such dispersion state is conducive to improving the hydrogen absorption and desorption performance of the composite. At the same time, Mg2Ni and Mg2NiH4, as inducible phases, can absorb and desorb hydrogen prior to Mg/MgH2system because of their lower hydrogenation and dehydrogenation temperatures, so they can be considered as “hydrogen pump” in the composite. Besides, XPS spectra of Fe, Ni and S elements of the composite after multiple rehydrogenation and dehydrogenation processes also confirm the stable catalytic system, where the peak location and peak strength are almost the same, and no significant changes in valence states(Fig. S7).

    Fig. 8. Schematic diagram of the hydrogenation/dehydrogenation processes of the MgH2-FeNi2S4 composite.

    According to the above results, the hydrogen absorption and desorption schematic diagram of the MgH2-FeNi2S4composite is used to reveal the catalytic mechanism of FeNi2S4on the MgH2, which is as shown in Fig. 8. In fact, in the previous preparation process, the Mg/Mg2Ni-MgS/Fe hydrogenation system and the MgH2/Mg2NiH4-MgS/Fe dehydrogenation system has been successfully constructed. During hydrogenation, because of the ahead hydrogen absorption of Mg2Ni before Mg and the existence of MgS/Fe, the diffusion path of hydrogen atoms is changed, so that the system can quickly absorb hydrogen at a lower temperature.During dehydrogenation, on one hand, Mg2NiH4with low hydrogen desorption enthalpy induces MgH2to dehydrogenate at a lower temperature [29,67]. On the other hand, since the transition metals Fe and Ni have a strong affinity for H, in order to fill the 3d orbitals lacking electrons, the transition metals can easily absorb hydrogen atoms and release H2molecules from the hydrogen storage system [68–70]. In addition, the catalytic MgS exists stably in the hydrogen storage system, and provides more interfaces and channels for the rapid diffusion of H with other catalytically active substances. It is worth noting that MgS and Fe also play a synergistic catalytic role on the kinetics of de/hydrogenation performance of MgH2[47]. In summary, Mg2Ni/Mg2NiH4, MgS and Fe form a heterogeneous catalytic system, which has a multiphase interface and provides more active sites and diffusion channels for hydrogen atoms. Mg2Ni/Mg2NiH4, MgS and Fe serve as catalytic active substances, in which Mg2Ni/Mg2NiH4acts as the “hydrogen pump” driving Mg/MgH2hydrogen absorption and discharge, and MgS/Fe provides the catalytic active center to accelerate the transfer of H [71,72]. Therefore, the synergistic catalytic action of Mg2Ni/Mg2NiH4,MgS and Fe improves the hydrogen storage performance of Mg/MgH2.

    4. Conclusion

    In this study, a hydrogen storage material of the MgH2-FeNi2S4composite were prepared by hydrogenation combustion synthesis and ball milling method. The MgH2/Mg2NiH4-MgS/Fe heterogeneous catalytic system has been successfully constructed and shows excellent hydrogen storage performance and cycling stability. The MgH2-FeNi2S4composite can absorb 6.09 wt% H2and release 1.92 wt% H2in 1 h at 573 K. Noticeably, the composite is able to absorb 4.02 wt% H2at 373 K under 3.0 MPa H2, and the initial hydrogen desorption temperature is also reduced by about 80 K.In terms of desorption kinetics, the desorption activation energy of the composite with the addition of catalyst decreases from 161.2 kJ·mol–1of the pure MgH2to 65.5 kJ·mol–1. After 10 cycles of hydrogen absorption and desorption, the hydrogenation and dehydrogenation capacity of the composite still reach 98.9% and 99.7%, respectively. The in-situ generated Mg2NiH4/Mg2Ni during hydrogenation/desorptionation acts as “hydrogen pump” in the composite, while the stable MgS and Fe can be used to provide continuous catalysis for the system, therefore, the synergistic catalytic action of Mg2Ni/Mg2NiH4, MgS and Fe improves the hydrogen storage performance of Mg/MgH2.This work extends the horizon of ternary transition metal sulfides with special structures as catalysts to improve the hydrogen absorption and release performance of MgH2, and provides new ideas for the development of efficient and stable catalysts to enhance the hydrogen storage performance of Mg-based hydrogen storage materials.Based on FeNi2S4, the influence of other transition metal sulfides on MgH2will also be investigated in future studies, and the introduction of carbon materials will be considered to disperse the entire heterogeneous catalytic system and provide more catalytic active sites.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    This work was supported by the National Natural Science Foundation of China (grant numbers 52071281 and 51971197), the Natural Science Foundation of Hebei Province (grant numbers E2019203161, E2019203414 and E2020203081), and Science and Technology Major project of Inner Mongolia (2020ZD0012).

    Supplementary materials

    Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.sipas.2020.100009.

    国产日韩一区二区三区精品不卡| 国产日韩欧美视频二区| 免费少妇av软件| 久久精品国产亚洲av涩爱| 久久人妻熟女aⅴ| 国产免费福利视频在线观看| 免费久久久久久久精品成人欧美视频 | 夜夜骑夜夜射夜夜干| 好男人视频免费观看在线| 女性被躁到高潮视频| 一本大道久久a久久精品| 女的被弄到高潮叫床怎么办| 99re6热这里在线精品视频| 亚洲精品乱码久久久久久按摩| 久久久久视频综合| 又黄又粗又硬又大视频| 午夜福利视频在线观看免费| 女人久久www免费人成看片| 亚洲精品,欧美精品| 欧美少妇被猛烈插入视频| 男人爽女人下面视频在线观看| 午夜福利视频在线观看免费| 中文字幕另类日韩欧美亚洲嫩草| 精品一区二区三区四区五区乱码 | 人妻系列 视频| 日韩制服丝袜自拍偷拍| 最近中文字幕2019免费版| 日韩中字成人| 亚洲 欧美一区二区三区| 欧美3d第一页| 亚洲精品一二三| 国产永久视频网站| 美女xxoo啪啪120秒动态图| 欧美xxxx性猛交bbbb| 欧美精品高潮呻吟av久久| 亚洲性久久影院| 国产精品蜜桃在线观看| 视频区图区小说| 男女边吃奶边做爰视频| 国产亚洲午夜精品一区二区久久| 少妇的丰满在线观看| 亚洲av.av天堂| 大码成人一级视频| 在线观看国产h片| 狂野欧美激情性bbbbbb| 18+在线观看网站| 久久免费观看电影| 美女xxoo啪啪120秒动态图| 建设人人有责人人尽责人人享有的| 久久毛片免费看一区二区三区| 精品福利永久在线观看| 丰满迷人的少妇在线观看| 亚洲欧美日韩另类电影网站| 18禁国产床啪视频网站| 久久99蜜桃精品久久| 搡老乐熟女国产| 在线免费观看不下载黄p国产| 少妇被粗大的猛进出69影院 | 黑人欧美特级aaaaaa片| 久久精品久久久久久久性| 久久亚洲国产成人精品v| 国产精品国产三级专区第一集| 男女国产视频网站| 黄网站色视频无遮挡免费观看| 国产又色又爽无遮挡免| 一区在线观看完整版| 久久久久久久精品精品| 自线自在国产av| 亚洲成色77777| 精品一区二区三区四区五区乱码 | 九草在线视频观看| 亚洲国产av影院在线观看| 大片免费播放器 马上看| 亚洲伊人色综图| 国产免费现黄频在线看| 大片电影免费在线观看免费| 曰老女人黄片| 人人妻人人爽人人添夜夜欢视频| 91国产中文字幕| 久久综合国产亚洲精品| 女人被躁到高潮嗷嗷叫费观| 日韩熟女老妇一区二区性免费视频| 国产一区二区在线观看日韩| 久久午夜福利片| 日本91视频免费播放| 欧美xxⅹ黑人| 国产高清国产精品国产三级| 妹子高潮喷水视频| 99国产综合亚洲精品| 美女国产视频在线观看| 国产片特级美女逼逼视频| 在线观看免费视频网站a站| 午夜日本视频在线| 久久精品久久精品一区二区三区| 男人舔女人的私密视频| 满18在线观看网站| av天堂久久9| 激情视频va一区二区三区| 亚洲,欧美,日韩| 成人手机av| 国产av国产精品国产| 久久久久久人妻| 夫妻性生交免费视频一级片| 女人精品久久久久毛片| 亚洲四区av| 又黄又爽又刺激的免费视频.| 你懂的网址亚洲精品在线观看| 国产精品国产三级国产av玫瑰| 国产亚洲精品久久久com| 青春草视频在线免费观看| 伊人亚洲综合成人网| 激情五月婷婷亚洲| 中文字幕亚洲精品专区| 中文字幕另类日韩欧美亚洲嫩草| 免费av中文字幕在线| 青春草国产在线视频| 91精品三级在线观看| 免费在线观看黄色视频的| 9191精品国产免费久久| 久久毛片免费看一区二区三区| 两性夫妻黄色片 | 国产成人一区二区在线| 少妇人妻 视频| 精品久久久精品久久久| 国产又色又爽无遮挡免| 精品久久久久久电影网| 亚洲国产精品一区三区| 女人精品久久久久毛片| 久久狼人影院| 十八禁网站网址无遮挡| 人成视频在线观看免费观看| 国产一区二区在线观看av| www.熟女人妻精品国产 | tube8黄色片| 老司机影院毛片| 一区二区av电影网| 欧美性感艳星| 男男h啪啪无遮挡| 五月天丁香电影| 男人舔女人的私密视频| 午夜福利,免费看| 欧美xxⅹ黑人| 亚洲成国产人片在线观看| 中文字幕av电影在线播放| 国产乱来视频区| 一区二区日韩欧美中文字幕 | 99久久中文字幕三级久久日本| 久久97久久精品| 亚洲欧美成人综合另类久久久| 男女国产视频网站| av在线app专区| 高清黄色对白视频在线免费看| 黑人猛操日本美女一级片| 999精品在线视频| 亚洲人成77777在线视频| 大片电影免费在线观看免费| 欧美亚洲日本最大视频资源| 插逼视频在线观看| 国产成人av激情在线播放| 日日啪夜夜爽| 18在线观看网站| 免费大片18禁| 99热6这里只有精品| 精品亚洲成国产av| 亚洲内射少妇av| 亚洲欧美日韩卡通动漫| 成年人午夜在线观看视频| 最近2019中文字幕mv第一页| 久久久久久久大尺度免费视频| 久久久a久久爽久久v久久| 亚洲欧美成人精品一区二区| 亚洲一码二码三码区别大吗| 国产精品国产三级国产av玫瑰| 大片免费播放器 马上看| 飞空精品影院首页| 日韩在线高清观看一区二区三区| 成人毛片a级毛片在线播放| 七月丁香在线播放| 99久久综合免费| 晚上一个人看的免费电影| 人人妻人人爽人人添夜夜欢视频| 汤姆久久久久久久影院中文字幕| 高清av免费在线| 黄网站色视频无遮挡免费观看| 波野结衣二区三区在线| 97在线视频观看| 午夜精品国产一区二区电影| 制服丝袜香蕉在线| 在线观看免费高清a一片| 免费少妇av软件| 香蕉国产在线看| 中文字幕精品免费在线观看视频 | √禁漫天堂资源中文www| 丝袜在线中文字幕| 亚洲精品乱码久久久久久按摩| 午夜福利视频精品| 18禁国产床啪视频网站| 五月玫瑰六月丁香| 在线亚洲精品国产二区图片欧美| 纵有疾风起免费观看全集完整版| av国产精品久久久久影院| 国产精品麻豆人妻色哟哟久久| 蜜臀久久99精品久久宅男| 久久国产亚洲av麻豆专区| 成年女人在线观看亚洲视频| 哪个播放器可以免费观看大片| 欧美 亚洲 国产 日韩一| 色5月婷婷丁香| 成人国产av品久久久| 色视频在线一区二区三区| 97在线人人人人妻| 日本色播在线视频| 久久精品夜色国产| 全区人妻精品视频| 日日爽夜夜爽网站| 国产精品人妻久久久久久| 少妇人妻久久综合中文| 成人二区视频| 亚洲一区二区三区欧美精品| 99九九在线精品视频| 少妇的逼水好多| 菩萨蛮人人尽说江南好唐韦庄| 欧美xxxx性猛交bbbb| 黑丝袜美女国产一区| 在线观看国产h片| 国产色爽女视频免费观看| 大香蕉久久成人网| 九色亚洲精品在线播放| 国产男女超爽视频在线观看| 国产片内射在线| 精品亚洲成a人片在线观看| 国产精品成人在线| 国产 精品1| 国产成人精品在线电影| 综合色丁香网| 18+在线观看网站| 国产精品国产三级国产av玫瑰| 亚洲成人一二三区av| 中文字幕人妻丝袜制服| 99久久综合免费| 欧美激情 高清一区二区三区| 亚洲一级一片aⅴ在线观看| 超碰97精品在线观看| 男女下面插进去视频免费观看 | 午夜福利网站1000一区二区三区| 在线观看免费高清a一片| 国产精品人妻久久久久久| 一本—道久久a久久精品蜜桃钙片| 激情视频va一区二区三区| 在线看a的网站| 高清不卡的av网站| 97超碰精品成人国产| 亚洲一级一片aⅴ在线观看| 亚洲精品456在线播放app| 成人毛片60女人毛片免费| av卡一久久| 欧美激情 高清一区二区三区| 亚洲精品色激情综合| 中文字幕制服av| 交换朋友夫妻互换小说| 老司机亚洲免费影院| 精品少妇黑人巨大在线播放| 精品一区在线观看国产| 男女无遮挡免费网站观看| 亚洲国产av新网站| 亚洲国产成人一精品久久久| 日产精品乱码卡一卡2卡三| 国产精品久久久久成人av| 国产欧美亚洲国产| 91久久精品国产一区二区三区| 男女午夜视频在线观看 | 巨乳人妻的诱惑在线观看| 国产一区二区在线观看日韩| 久久精品久久久久久久性| 天天操日日干夜夜撸| 精品国产国语对白av| 深夜精品福利| 国产高清国产精品国产三级| 在线观看免费视频网站a站| 美女国产视频在线观看| 99久久人妻综合| 亚洲av国产av综合av卡| 国产综合精华液| 99久国产av精品国产电影| 久久精品久久久久久噜噜老黄| 在线观看三级黄色| 国产精品一区www在线观看| 午夜视频国产福利| 欧美xxⅹ黑人| 十分钟在线观看高清视频www| 最近手机中文字幕大全| 日韩视频在线欧美| 2022亚洲国产成人精品| 日本爱情动作片www.在线观看| 老熟女久久久| 欧美成人午夜精品| a级片在线免费高清观看视频| 国产精品一国产av| 亚洲婷婷狠狠爱综合网| 日韩视频在线欧美| 国产精品成人在线| 激情视频va一区二区三区| 精品人妻偷拍中文字幕| 亚洲av免费高清在线观看| 国产高清不卡午夜福利| 卡戴珊不雅视频在线播放| 97超碰精品成人国产| 有码 亚洲区| 欧美精品高潮呻吟av久久| 国产麻豆69| 中文字幕制服av| 国产精品无大码| 最近2019中文字幕mv第一页| 国产女主播在线喷水免费视频网站| 9191精品国产免费久久| 成人亚洲欧美一区二区av| 亚洲 欧美一区二区三区| 国产毛片在线视频| 丰满迷人的少妇在线观看| 午夜视频国产福利| 国产亚洲最大av| 欧美97在线视频| 精品一区二区三区四区五区乱码 | 一区二区日韩欧美中文字幕 | av免费在线看不卡| 久久午夜福利片| 国产精品国产三级国产专区5o| 成年动漫av网址| 午夜日本视频在线| 国产精品三级大全| 999精品在线视频| 国产熟女欧美一区二区| 男女国产视频网站| 天天躁夜夜躁狠狠躁躁| 成年美女黄网站色视频大全免费| 狂野欧美激情性xxxx在线观看| 九九爱精品视频在线观看| 纵有疾风起免费观看全集完整版| 观看美女的网站| 丝袜喷水一区| 亚洲精品,欧美精品| 成人无遮挡网站| 久久精品aⅴ一区二区三区四区 | 高清黄色对白视频在线免费看| 一本大道久久a久久精品| 国产av精品麻豆| 美女福利国产在线| videosex国产| 精品国产国语对白av| 国产精品一区二区在线观看99| 91精品国产国语对白视频| www.av在线官网国产| 午夜免费鲁丝| 日韩,欧美,国产一区二区三区| 免费大片18禁| 十八禁网站网址无遮挡| 日日摸夜夜添夜夜爱| 久久久久精品人妻al黑| 一区二区三区四区激情视频| 看免费成人av毛片| 少妇猛男粗大的猛烈进出视频| 九草在线视频观看| 91精品伊人久久大香线蕉| 看免费成人av毛片| 99久久综合免费| 夫妻性生交免费视频一级片| 亚洲综合精品二区| 国产免费一级a男人的天堂| 国内精品宾馆在线| 22中文网久久字幕| 欧美成人午夜精品| 亚洲色图 男人天堂 中文字幕 | 日韩一本色道免费dvd| 亚洲欧美一区二区三区国产| 亚洲精品乱码久久久久久按摩| 狂野欧美激情性bbbbbb| 亚洲精品,欧美精品| 性色av一级| 91国产中文字幕| 成人二区视频| 亚洲欧洲国产日韩| 极品少妇高潮喷水抽搐| 亚洲国产欧美在线一区| 精品一品国产午夜福利视频| 日本猛色少妇xxxxx猛交久久| 久久这里有精品视频免费| 久久久精品94久久精品| 色5月婷婷丁香| 午夜免费鲁丝| 久久女婷五月综合色啪小说| 最后的刺客免费高清国语| 99视频精品全部免费 在线| 99九九在线精品视频| 熟女电影av网| 丝袜脚勾引网站| a级毛片在线看网站| 秋霞在线观看毛片| 国产熟女午夜一区二区三区| 精品99又大又爽又粗少妇毛片| 国产欧美日韩综合在线一区二区| 丝瓜视频免费看黄片| 日日撸夜夜添| 成年女人在线观看亚洲视频| 亚洲欧洲精品一区二区精品久久久 | 久久热在线av| 久久久久精品人妻al黑| 精品卡一卡二卡四卡免费| 国产有黄有色有爽视频| 欧美日韩av久久| 亚洲av中文av极速乱| 天堂俺去俺来也www色官网| 国产一区二区三区av在线| 国产成人欧美| 久久这里只有精品19| 91成人精品电影| 自拍欧美九色日韩亚洲蝌蚪91| 久久婷婷青草| 免费大片18禁| 波野结衣二区三区在线| 老司机影院成人| 91午夜精品亚洲一区二区三区| 麻豆精品久久久久久蜜桃| 国产精品 国内视频| 老司机亚洲免费影院| 亚洲国产日韩一区二区| 色哟哟·www| 精品亚洲成a人片在线观看| 国产成人精品久久久久久| av卡一久久| 熟女av电影| 亚洲精品乱久久久久久| 女人被躁到高潮嗷嗷叫费观| 国产福利在线免费观看视频| 在线观看美女被高潮喷水网站| 乱码一卡2卡4卡精品| 中文字幕人妻熟女乱码| 99热全是精品| 人成视频在线观看免费观看| 亚洲综合精品二区| 日韩熟女老妇一区二区性免费视频| 久久热在线av| 日韩欧美一区视频在线观看| 桃花免费在线播放| 天天躁夜夜躁狠狠躁躁| 欧美亚洲 丝袜 人妻 在线| www日本在线高清视频| 色94色欧美一区二区| 校园人妻丝袜中文字幕| 赤兔流量卡办理| 国产精品三级大全| 在线精品无人区一区二区三| 热re99久久国产66热| 一级毛片黄色毛片免费观看视频| 国产精品久久久久久精品电影小说| 18在线观看网站| 18禁动态无遮挡网站| 桃花免费在线播放| 国产午夜精品一二区理论片| 亚洲色图综合在线观看| 一级毛片黄色毛片免费观看视频| 少妇人妻精品综合一区二区| 少妇被粗大猛烈的视频| 又黄又爽又刺激的免费视频.| 黑人猛操日本美女一级片| 国产乱人偷精品视频| 日韩伦理黄色片| 在线观看人妻少妇| 国产片内射在线| 久久久久国产网址| 婷婷色综合www| 伊人亚洲综合成人网| 国产无遮挡羞羞视频在线观看| 精品国产一区二区久久| 老女人水多毛片| 国国产精品蜜臀av免费| 捣出白浆h1v1| 一本久久精品| 2022亚洲国产成人精品| 久久精品国产a三级三级三级| 777米奇影视久久| 国产伦理片在线播放av一区| videossex国产| 精品人妻偷拍中文字幕| 亚洲国产成人一精品久久久| 在线观看三级黄色| 99re6热这里在线精品视频| 午夜91福利影院| 日韩一本色道免费dvd| 晚上一个人看的免费电影| 18禁裸乳无遮挡动漫免费视频| 另类亚洲欧美激情| 嫩草影院入口| 天天操日日干夜夜撸| 999精品在线视频| 一本—道久久a久久精品蜜桃钙片| 另类亚洲欧美激情| 美女内射精品一级片tv| 日韩一本色道免费dvd| 中文字幕亚洲精品专区| 一级片'在线观看视频| 99国产精品免费福利视频| 国产在视频线精品| 大香蕉久久成人网| 精品亚洲成国产av| 两个人看的免费小视频| 熟女电影av网| 天堂俺去俺来也www色官网| 日本vs欧美在线观看视频| 80岁老熟妇乱子伦牲交| av卡一久久| 久久久久久久久久久久大奶| 亚洲av综合色区一区| 国产精品国产三级国产av玫瑰| 亚洲欧美成人精品一区二区| 午夜91福利影院| 亚洲精品日韩在线中文字幕| 久久精品国产a三级三级三级| √禁漫天堂资源中文www| 精品酒店卫生间| 老司机亚洲免费影院| 制服诱惑二区| 在线观看一区二区三区激情| 国产视频首页在线观看| 大陆偷拍与自拍| 亚洲欧美日韩另类电影网站| 极品人妻少妇av视频| 国产午夜精品一二区理论片| 久久午夜福利片| 亚洲av国产av综合av卡| av黄色大香蕉| 成人影院久久| 久久久久久久久久人人人人人人| 一边摸一边做爽爽视频免费| 七月丁香在线播放| 国产成人一区二区在线| 99久国产av精品国产电影| 国产精品一国产av| 国产色婷婷99| 久久精品久久精品一区二区三区| 午夜福利网站1000一区二区三区| 久久久久久人人人人人| √禁漫天堂资源中文www| 看十八女毛片水多多多| 欧美xxxx性猛交bbbb| 亚洲国产av新网站| 宅男免费午夜| 99久久精品国产国产毛片| 亚洲色图 男人天堂 中文字幕 | 看非洲黑人一级黄片| 男女无遮挡免费网站观看| 高清不卡的av网站| 男女啪啪激烈高潮av片| 极品少妇高潮喷水抽搐| 久久人人爽av亚洲精品天堂| 大话2 男鬼变身卡| 如何舔出高潮| 国产亚洲午夜精品一区二区久久| 国产精品一区二区在线观看99| 秋霞在线观看毛片| 18禁裸乳无遮挡动漫免费视频| 欧美激情国产日韩精品一区| 国产片特级美女逼逼视频| 久久精品久久久久久噜噜老黄| 人妻 亚洲 视频| 王馨瑶露胸无遮挡在线观看| 欧美亚洲 丝袜 人妻 在线| av在线老鸭窝| 亚洲欧美日韩卡通动漫| 一级,二级,三级黄色视频| 97在线人人人人妻| 观看美女的网站| 亚洲精品aⅴ在线观看| 亚洲美女搞黄在线观看| 中国国产av一级| 日本-黄色视频高清免费观看| 飞空精品影院首页| 欧美国产精品一级二级三级| 国产av精品麻豆| 亚洲精品中文字幕在线视频| 国产男人的电影天堂91| 看十八女毛片水多多多| 国产精品国产三级国产专区5o| 精品久久蜜臀av无| 伦理电影大哥的女人| 韩国av在线不卡| av又黄又爽大尺度在线免费看| 国产一区二区三区综合在线观看 | 69精品国产乱码久久久| 日韩电影二区| 久久精品国产综合久久久 | av黄色大香蕉| 午夜免费鲁丝| a级毛片在线看网站| 在线观看人妻少妇| 亚洲,一卡二卡三卡| 国产黄色免费在线视频| 国产高清三级在线| 黄色毛片三级朝国网站| 欧美日韩成人在线一区二区| av在线老鸭窝| 91精品国产国语对白视频| 欧美日韩成人在线一区二区| 欧美97在线视频| 好男人视频免费观看在线| 制服诱惑二区| av视频免费观看在线观看| 夜夜骑夜夜射夜夜干| 精品国产乱码久久久久久小说| 97超碰精品成人国产| 国产老妇伦熟女老妇高清| 精品一区二区三区四区五区乱码 | 久久ye,这里只有精品| 狠狠婷婷综合久久久久久88av| 中文字幕制服av| 成人18禁高潮啪啪吃奶动态图|