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

    Catalytic effect comparison of TiO2 and La2O3 on hydrogen storage thermodynamics and kinetics of the as-milled La-Sm-Mg-Ni-based alloy

    2021-02-24 13:16:18YnghunZhngXinWeiWeiZhngZemingYunJinlingGoHuipingRen
    Journal of Magnesium and Alloys 2021年6期

    Ynghun Zhng ,Xin Wei ,Wei Zhng ,Zeming Yun ,Jinling Go ,Huiping Ren

    a Key Laboratory of Integrated Exploitation of Baiyun Obo Multi-Metal Resources,Inner Mongolia University of Science and Technology,No.7 Aerding Road,Kun District,Baotou 014010,China

    b Department of Functional Material Research,Central Iron and Steel Research Institute,No.76 Xueyuannan Road,Haidian District,Beijing 100081,China

    c Weishan Cisri-Rare Earth Materials Co.,Ltd.,Jining 277600,China

    Abstract In this investigation,mechanical grinding was applied to fabricating the Mg-based alloys La7Sm3Mg80Ni10+5wt.% M(M=None,TiO2,La2O3)(named La7Sm3Mg80Ni10–5M(M=None,TiO2,La2O3)).The result reveals that the structures of as-milled alloys consist of amorphous and nanocrystalline.The particle sizes of the added M(M=TiO2,La2O3)alloys obviously diminish in comparison with the M=None specimen,suggesting that the catalysts TiO2 and La2O3 can enhance the grinding efficiency.What’s more,the additives TiO2 and La2O3 observably improve the activation performance and reaction kinetics of the composite.The time required by releasing 3wt.% hydrogen at 553,573 and 593K is 988,553 and 419s for the M=None sample,and 578,352 and 286s for the M=TiO2 composite,and 594,366,301s for the La2O3 containing alloy,respectively.The absolute value of hydrogenation enthalpy change |?H| of the M(M=None,TiO2,La2O3)alloys is 77.13,74.28 and 75.28 kJ/mol.Furthermore,the addition of catalysts reduces the hydrogen desorption activation energy().

    Keywords: Mg-based alloy;Ball milling;Catalysts;Hydrogen storage kinetics.

    1.Introduction

    Hydrogen is considered the most promising fuel applied in fuel-cell vehicle because of the wide applicability,high application efficiency,safety,environmentally friendly and inexhaustible reserves characteristics of it[1–3].An accepted view considers that the widespread utilization of clean renewable energy vehicles is an effective measure to deal with environmental pollution and excessive consumption of the limited fossil fuels because transportation consumes about a quarter of global energy while there is about 23% carbon dioxide emission originating from the combustion of fossil fuels in the world [4–6].However,the chemical energy of unit volume hydrogen is low,which limits the commercialization of it.Therefore,a lot of investigations pay attention to dealing with this problem and finding the safe,suitable hydrogen storing materials that can meet the needs of using hydrogen in the field of motor vehicle and electronic products[7–11].The ideal hydrogen storage methods should have the following characteristics,large capacity,low-temperature fast adsorption rate and high cycle survivability [12].In present,hydrogen storage methods include high-pressure storage,liquid storage,physical adsorption storage and hydride storage.Generally,the gas storage needs 70 MPa pressure for reaching only 4.8wt.% H2capacity [13].The condition of liquid storage is low-temperature(20 K)that brings about huge energy consumption(about 30% for filling)and short boiling storage duration [14].Solid hydrogen storage is reckoned as the best one among the various methods [15,16].

    Among the metal hydrides,MgH2with high hydrogen storage capacity(7.6 wt.%),low cost,non-toxicity,abundant reserves and superior reversibility characteristics is the most promising material for hydrogen storage [17,18].Nonetheless,the attempt of the widespread application of commercial MgH2is retarded by its some disadvantages,such as high thermodynamic stability(i.e.the high strength of the Mg-H bonds),the complicated activation procedure,the high dissociation temperature and the sluggish hydrogenation/dehydrogenation kinetics [19,20].To copy with the dilemma,many studies have been conducted to reduce the decomposition temperature,accelerate the adsorption kinetics and change the reaction thermodynamics by the method of grain refinement [11,21],mixing catalysts [22–24],ballmilling [25,26],alloying [27,28],surface modification [29,30]and the other [31].The mechanical grinding and meltspinning [32,33]are very efficient means to reduce the particle or grain size.Especially,the mechanical milling can realize the refinement of grains and particles simultaneously[38].In addition,mechanical milling brings on the formation of crystal defects and modifies the superficial characteristics of alloy particles [1].Hence,mechanical grinding is reckoned as the best way to improve the hydrogen storing performance of magnesium based materials [34].Besides,refining the particle to very small size will introduce the capillarity effect that can ameliorate the hydrogen absorption and desorption thermodynamics in theory.The calculation results indicate that a particle radius on the order of 5 nm will reduce the dehydrogenation enthalpy of pure Mg by about 10%[35].As well known,there is considerable energy used for the hydrogen dissociation on the magnesium surface and this segment is regarded as a rate controlling factor in the procedure of hydrogen absorption [36].Transition metals either in their pure form(e.g.:Ni,Ti,Nb,Fe,Co,Al etc.)[37,38]or as oxides(e.g.:Nb2O5,Fe2O3,TiO2etc.)[39–42],hydrides(e.g.:TiH2,ZrFe2Hx,etc.)[43,44],fluoride(e.g.:FeF3,TiF3,NiF2,NbF5etc.)[45,46]or intermetallics [47,48],can act as the catalysts because they can weaken this dissociation energy.As studied by Du et al.[49]and Pozzo and Alfe [50],the hydrogen dissociation energy on magnesium surface is 1.15 eV but it can be decreased to 0.03,0.06,0.56,0.39 eV by the addition of Co,Ni,Cu,Pd,severally.Daryani et al.[51]researched that adding 6 mol% TiO2could improve the hydrogen absorbing kinetics and reduce the decomposition temperature of as-milled Magnesium hydride by 100 K.According to the research of Shahi et al.[52],the composite MgH2?5 wt.%Ni absorbs 5.0 wt.%H2hydrogen at the temperature of 443 K in 15 min and at the temperature 613 K it starts to decompose.As investigated by Hou et al.[53],Mg2NiH4with the catalyst composed of MWCNTs and TiF3has the 503 K(516.6 K for pure Mg2NiH4)hydrogen releasing temperature(TD)and the activation energy(Ea)of it is 53.24 kJ/mol(90.13 kJ/mol for pure Mg2NiH4).In particular,adding appropriate rare earths or their oxides can obviously make the Mg-based hydrides unstable and accelerate the rate of dehydrogenation reaction [54,56].Lass [56]found that the Mg85Ni15-xMx(M=La,x=0 or 5)alloys possess a lower enthalpy change in the reaction of producing MgH2and Mg2NiH4.On the basis of the investigation of Luo et al.[55],the element Y is beneficial to improve the thermodynamics property of magnesium based materials and the composite Mg90In5Y5has a lower?H(about 62.9 kJ/(mol H2))in comparison with the Mg95In5binary alloy(about 67.9 kJ/(mol H2))and pure Mg(about 74.9 kJ/(mol H2)).Sadhasivam et al.[57]found that the original desorption temperature of the composite MgH2?5 wt.% Mm-oxide was reduced by 76 K from 654 to 578 K.Kalinichenka et al.[58]researched the improved reaction kinetics of Mg90Ni8RE2(RE=Y,Nd,Gd)and found that the activated Mg90Ni8RE2could reversible absorb and release 5.5 wt.% H2within 20 min.

    According to our investigation on REMg11Ni(RE=Sm,Y)+5 wt.% M(M=MoS2,CeO2)composites,the additives MoS2and CeO2play a catalytic role in improving hydrogen storing performance [59,60]and 5 wt.% addition of catalysts is optimal as studied in this reference [33].It must be very interesting to compare the effects of TiO2and La2O3additives with high hardness on the hydrogen storing performance of ball milling magnesium based materials.Thereby,the alloys La7Sm3Mg80Ni10–5 M(M=None,TiO2,La2O3)were fabricated by mechanical milling.The thermodynamics and dynamics of the experimental alloys were investigated.A comprehensive comparison of the impacts of different catalysts on the structure and hydrogen storing performances of the alloys is conducted.

    2.Experimental

    The La7Sm3Mg80Ni10material was fabricated by inductive melting La,Sm,Mg,Ni(purity ≥99.9%)under 0.04 MPa He(purity ≥99.999%)to inhibit the volatilization of magnesium.To compensate the melting losses,additional magnesium(8 wt.%)and RE(RE=La,Sm)(5 wt.%)are required.The above-mentioned materials were all provided by CISRI Corporation.A Varian Liberty 100 inductively-coupled plasma(ICP)was applied to determining the chemical composition of experimental alloys.Then the ingot was mechanically crushed and ground to the 200–400 meshes powders.The obtained power with 5 wt.% TiO2or La2O3(purity ≥99.9%)catalyst was mechanically ground by a mill crusher at the speed of 350 rpm(the weight ratio of specimen and balls is 1:40).The milling duration is set at 20 h.Thus,the chemical compositions of the as-milled powder were La7Sm3Mg80Ni10–5 M(M=None,TiO2,La2O3).In order to heat dissipation and reduce the cold welding of powder in the process of milling,the working mechanism of ball mill is to stop half an hour every 3 h and the powder adhered to the milling chamber walls and grinding balls needs to be scrapped in due course of time all of these operations were operated under the protective atmosphere of Ar.

    X-ray diffraction(XRD)(D/max/2400)determined the phase structures and compositions of the alloys.The experimental parameters were 40 kV,160 mA,and 2°/min with 2θchanging from 20° to 90°.The radiation was CuKα1filtered by graphite.The particles morphology observation was completed by a scanning electron microscope(SEM)(QUANTA 400).A high resolution transmission electron microscope(HRTEM)(JEM-2100F,operated at 200 kV)was utilized to the characterization of microstructure and crystalline state.

    Fig.1.XRD profiles of the as-milled(20 h)La7Sm3Mg80Ni10?5M(M=none,TiO2,La2O3)alloys:(a)before hydrogen absorption,(b)after hydrogen absorption,(c)after hydrogen desorption.

    Hydrogenation and dehydrogenation kinetics curves of the as-milled specimens were tested by automatic Sieverts apparatus.Prior to measuring,the sample need to be activated by six hydriding/dehydriding cycles(633 K and original hydrogen pressure of 3 MPa for hydrogen absorption,633 K and 1×10?4MPa original pressure for hydrogen desorption).The temperature of hydrogen absorption was set as 473,513,533,553,573,593,613 and 633 K,severally,while 553,573,593,613 and 633 K for hydrogen desorption.The setting of initial hydrogen pressure is the same as activation.The sample mass required for every determination was 300 mg.Nonisothermal hydrogen desorption property was researched by utilizing thermogravimetry(TGA)and differential scanning calorimetry(DSC)(SDTQ600)whose heating rates were 5,10,15 and 20 K/min.

    3.Results and discussion

    3.1.Microstructure characteristics

    Fig.1 gives the X-ray diffraction of the as-milled La7Sm3Mg80Ni10–5 M(M=none,TiO2,La2O3)composites before and after hydrogen absorption and desorption under the condition of 633 K,3 MPa and 633 K,1×10?4MPa,severally.ICDD(International Centre for Diffraction Data)identification of X-ray diffraction patterns shows that the as-milledM=none specimen consists of the major phase La2Mg17and secondary phases Mg2Ni and La2Ni3.The addition of catalysts TiO2and La2O3do not introduce any new phase,indicating that these additives are not involved in the reaction with alloy.Moreover,it is visible to observe the broadened diffraction peaks representing the typical nanocrystalline and amorphous structures of as-milled specimens in comparison with that of the as-cast one(XRD patterns are not show here).After hydrogen absorption,the diffraction peaks get narrow and sharp.Meanwhile there are four hydrides become visible and emerge in the specimens,including MgH2,Mg2NiH4,LaH3and Sm3H7.The reaction relationship between the elements is as follows:

    La2Mg17+H2→LaH3+MgH2

    Mg2Ni+H2→Mg2NiH4

    La2Ni3+La2Mg17+H2→LaH3+Mg2NiH4+MgH2

    Sm+H2→Sm3H7

    After dehydrogenated,the four phases Mg,Mg2Ni,LaH3and Sm3H7can be found.Evidently,the LaH3and Sm3H7phases are not decomposed because of the high thermal steadiness of them.Hence,the hydrogen desorption reactions are summarized as the following two equations:

    MgH2→Mg+H2

    Mg2NiH4→Mg2Ni+H2

    According to the above inference,we can see that in the process of hydrogenation and dehydrogenation,the reversible reactions of activated composites include

    Mg+H2?MgH2

    Mg2Ni+H2?Mg2NiH4

    Fig.2.SEM morphologies of the as-milled(20 h)La7Sm3Mg80Ni10?5M(M=none,TiO2,La2O3)alloys:(a) M=none,(b) M=TiO2,(c) M=La2O3.

    Through a careful observation,we find that the width of the XRD peak narrows down after dehydrogenation compared with that after hydrogen absorption,which was owing to the cell volume reduction and stress relief rendered by hydrogen desorption.As found by Montone et al.[61],the volume of a metallic Mg atom is about 33% smaller than that of Mg atom in MgH2.It has been reported in the literature [62,63]that lattice distortion along with expansion and contraction of cell volume are inevitable in hydrogen storage materials during hydrogen absorption and desorption,which will cause many lattice defects such as vacancy and dislocation.The formation of the defects will have a beneficial effect on the hydrogen absorption and desorption property of the alloy.Mechanical milling of Mg-based alloy with TiO2and La2O3catalysts creates the defects on the surface and inside the magnesium matrix,which generate reactive clean surfaces and shrink the particle size of Mg.The creation of defects facilitates nucleation,the production of the reactive clean surface enhances the superficial reactivity,and the diminution of particle decreases the diffusion distances of hydrogen atoms.These effects ameliorate the hydriding and dehydriding kinetics of magnesium based alloy significantly.

    The morphologies of the as-milled La7Sm3Mg80Ni10–5 M(M=none,TiO2,La2O3)alloy powder are observed by SEM and presented in Fig.2.As can be observed,the alloy particles have the typical morphology of ball-milling powder and the size of them is in a range of 1–10 μm.After careful observation,it is failed to find the catalysts TiO2and La2O3particles,which indicates that the TiO2and La2O3particles are not appear on the superficial part of alloy,but is wrapped in them.Evidently,the agglomeration tendency of the as-milled particles of theM=TiO2andM=La2O3alloys was decreased(Fig.2b and c)with smaller size than that ofM=none alloy.It means that adding a certain amount of TiO2and La2O3can significantly improve the efficiency of ball milling.After comparing the particles with different catalysts,we found no obvious difference in particle size,suggesting that two catalysts have similar effect on the efficiency of ball milling.As considered by Floriano et al.[64],some catalysts with high hardness,e.g.La2O3,CeO2,TiO2Nd2O5,etc.can act as lubricants,dispersants and/or cracking agents in the procedure of milling and are helpful to further reduce refine the particles of as-milled alloy.A very similar result also appears in the investigation of Daryani et al.[51]and Aguey-Zinsou et al.[65].

    The HRTEM micrographs and SAED(Selected Area Electron Diffraction)patterns of the as-milled La7Sm3Mg80Ni10?5 M(M=none,TiO2,La2O3)materials are presented in Fig.3.We have noticed the nanocrystalline and amorphous structures of the as-milled alloys and the emergence of crystal defects.The SAED patterns also prove the existence of La2Mg17,Mg2Ni and La2Ni3phases,and there is no any new phase caused by adding TiO2and La2O3.After hydrogen absorption,the structures of as-milled alloy still are amorphous and nanocrystalline(Fig.3b,e and h),but there is an observably decrease in amorphous phase,meaning that the dehydrogenation promotes the crystallization reaction.Four hydrides MgH2,Mg2NiH4LaH3and Sm3H7also can be identified after hydrogen desorption by SAED patterns.According to Fig.3(c),(f)and(i),we can see that the alloys after hydrogen desorption exhibit an entirely crystal structure,and the size of grain evidently increase,and Pukazhselvan et al.[66]also had the similar report.The SAED rings of dehydrogenated alloys reflect the existence of Mg,Mg2Ni,LaH3and Sm3H7.Apparently,it is consistent with the result of XRD,the LaH3and Sm3H7phases still exist after dehydrogenation.In addition,it is found from Fig.2 that the LaH3,Sm3H7,TiO2and La2O3nanoparticals distribute in Mg matrix dispersedly and uniformly,which is considered to be the preferred nucleation sites for hydride formation/decomposition.The phase interfaces of LaH3(or Sm3H7,TiO2and La2O3)/Mg(or MgH2)provide channels for the diffusion of hydrogen atoms.Therefore,the additives can be regarded as catalysts to improve the hydrogen storage performance of magnesium and Mg-based alloy [67].

    Fig.3.HRTEM micrographs and SAD patterns of the as-milled La7Sm3Mg80Ni10?5M(M=none,TiO2,La2O3)alloys:(a),(b),(c) M=none alloy before hydrogen absorption,after hydrogen absorption,after hydrogen desorption;(d),(e),(f) M=TiO2 alloy before hydrogen absorption,after hydrogen absorption,after hydrogen desorption;(g),(h),(i) M=La2O3 alloy before hydrogen absorption,after hydrogen absorption,after hydrogen desorption.

    3.2.Activation behaviors and hydrogen absorption kinetics

    In this investigation,it was found that the alloy powder prepared by traditional mechanical milling can hardly absorb hydrogen because the alloy powder exposed to air easily forms an oxide film on the particle surface that blocks the contact between H2molecules and alloy surface and prevents H2from dissociating into H atoms.As well known,this dissociation process is the basic step in the phase transformation from metallic Mg to MgH2and it is necessary for the incorporation of H atoms into the Mg lattice.Fortunately,it is found that when the alloy sample is kept under proper temperature and hydrogen pressure for a long time,the oxide film formed can be broken gradually,which results in exposing the fresh alloy surface and restoring the hydrogen absorption capability of the alloys.This process is called as activation.The activation performance of specimens is greater if it needs less cycle numbers.Fig.4 demonstrates the isothermal hydrogen absorption and desorption curves of the activated La7Sm3Mg80Ni10–5 M(M=none,TiO2,La2O3)materials.It is found from Fig.4(a),(b)and(c),the alloys are almost fully activated after the first cycle since the activation curves of the next five cycles are almost identical.It is noted that the time spent on the first hydrogen absorption to a saturated state is long.It takes 21716 s for theM=none alloy,14568 s for theM=TiO2alloy and 13340 s for theM=La2O3alloy to achieve the saturated capacity of 5.15 wt.%,5.052 wt.% and 4.916 wt.%,severally,suggesting that the activation property of the alloy is considerably improved by adding TiO2and La2O3.For a given hydrogen absorption capacity of 4 wt.%,by which the time required is 4820,4624 and 4758 s corresponding to the as-milledM=none,M=TiO2andM=La2O3alloys,respectively.It indicates that the hydrogenation rate is in the orderM=TiO2>M=La2O3>M=none.The improved activation ability by adding TiO2and La2O3was attributed to the modified particles surface state and the increased defect density of the crystals resulted from adding catalysts TiO2and La2O3.The catalyst nanoparticles distributing on the particles surface significantly increase the dissociation rate(limiting factor of hydrogen absorption rate)of hydrogen molecules.

    Generally speaking,the first activation reaction is a long procedure [68]and in this process,the H atoms penetrate the formed oxide layer to form mental hydrides.What’s more,the attendant mechanical stress and lattice distortion are unfavorable to the absorption of hydrogen [69].As well known,the nucleation of MgH2on the superficial sites of alloy is retarded by the thin oxide layers [70].Although the operation is conducted at inert gas atmospheres the oxide layers with 3–4 nm thickness still can easily form [69].The sluggish dissociation of H2on the alloy surface is another reason to explain the slow hydrogen absorption rate [71].The dissociation on pure Mg surface needs high energy [48].Besides,the diffusion of H atoms in the metal hydrides is difficult[56,57].The growth rate of MgH2is decided by the hydrogen pressure due to the fact that higher pressure provides the greater the thermodynamics driving force for the hydrogen absorption.Nevertheless,if the original hydrogenation process is fast enough,a superficial layer of magnesium hydride will form to retard the hydrogen permeation [71].Because hydrogen diffuses along the interfaces but not along the Mg hydride layer[72],the MgH2hydride grows up in the form of slow Mg/Mg hydride interface movement.When the thickness reaches a certain value(30–50 μm),the hydrogen absorption reaction stops [73],indicating that powdered magnesium used for hydrogenation changes into massive magnesium.So the hydrogen absorption rate is affected by the powder size [74].The hydrogenation kinetics is markedly enhanced after the first hydrogen absorption and desorption cycle.With the increase in the cycle number,the hydrogen absorption kinetics curves have little change,which means the great activation performance of experimental alloys.Noticeably,the hydrogen absorbing capacity of all the specimens after first cycle first is no more than 4.8 wt.%,which represents a visible decline in the capacity.The formation of stable hydrides LaH3and Sm3H7is most likely responsible for the 0.25 wt.% capacity loss.

    Fig.4.The isothermal hydrogenation/dehydrogenation curves of the as-milled La7Sm3Mg80Ni10–5M(M=none,TiO2,La2O3)alloys in the cycle activation:(a),(b),(c) M=none, M=TiO2 and M=La2O3 hydrogen absorption;(d),(e),(f) M=none, M=TiO2 and M=La2O3 hydrogen desorption.

    The dehydrogenation curves of the alloys are provided in Fig.4(d),(e)and(f).It is visible that the dehydrogenation rate is fast and the dehydrogenation kinetics of the alloy was markedly improved by adding TiO2and La2O3.In particular,the first hydrogen desorption took less time.For a given hydrogen desorption capacity of 3 wt.%,by which the time required is 193,162 and 175 s corresponding to the as-milledM=none,M=TiO2andM=La2O3alloys,respectively.Evidently,the dehydrogenation rate is in the orderM=TiO2>M=La2O3>M=none.The improved activation performance is deemed to be related to the decreased particle size,the surface modification and the weakening effect of Mg-H bond strength resulted from the additive TiO2or La2O3.The improvement of thermodynamics performance is directly associated with the weakening of Mg-H bond strength.The co doping of multi elements,especially the transition element(or their compounds)and rare element(or their compounds)is beneficial to the reduction of thermal stability of MgH2[75,76].The reduction of particle observably enhances the decomposition rate of H2on the particle surface and is beneficial to the H atoms diffusion thus enhance the activation performance [77].Particularly,because TiO2and La2O3are high hardness particles,they are likely to cut into the alloy particles and form a new interface under the action of high impact stress in the ball-milling process,which may become the nucleation sites of hydride,acting as rapid paths for atoms diffusion [78].So the addition of TiO2and La2O3not only enhance the efficiency of mechanical milling,make the particle size decrease but also modify the surface of alloy particles,make the nucleation of hydrides more easy.

    After the activation treatment,the hydrogen absorbing and desorbing properties of experimental composites were improved significantly.It is necessary to explore the change of structures in the process of activation.With the help of SEM,the morphological variations of the experimental composites before and after activation process are provided in Fig.5.Clearly,the particles show irregular morphologies with the very rough surface.After six hydriding and dehydriding cycles,the particle morphologies of the alloy have a dramatically change.It is very evident that many cracks appear on powder surface due to the lattice stress forming in the process of hydrogen absorption.When the lattice stress exceeds the fracture strength of the material,the pulverization of the alloy is inevitable and results in the improved properties.Through the above structural analysis,we believe that the activation is significant to the formation and decomposition of hydrides,the oxide film on the surface of alloy particles breaks,and along with the cracking of alloy particles,the specific surface area of the alloy is increased,thus improving the hydrogen absorbing and desorbing properties.

    To explore the influence of the different catalysts,the hydrogen absorption curves of the as-milled La7Sm3Mg80Ni10–5 M(M=none,TiO2,La2O3)composites were tested at different temperatures from 423 to 633 K and 3 MPa,as shown in Fig.6.It is noted that at the initial stage corresponding to the rapid formation of hydride layer near the surface,the rate is very fast and hydrogen absorption capacity can reach more than 85% of saturated capacity in less than 200 s,while in the following stage it takes long time to achieve the saturated state due to the hindrance of formed hydride layer acting on the hydrogen diffusion.Freidlmeier et al.[79]considered that when the thickness of hydride layer reached to a certain value(100 nm),the rate of hydrogen absorption tended to 0.Aiming at investigate the hydrogen absorption kinetics more deeply,the time required to absorb 4 wt.% hydrogen was calculated and compared.As obtained from Fig.6,the spent time is 108,67 and 56 s at 473,513 and 533 K for theM=none specimen,96,62 and 48 s for theM=TiO2composite,and 103,64 and 53 s for theM=La2O3material,respectively.Apparently,the hydrogenation rate of the composites is in order ofM=TiO2>M=La2O3>M=none,which suggests that the added TiO2and La2O3notably ameliorate the hydrogen absorption kinetics,but this favorable effect decreases rapidly with hydrogen absorption temperature rising.Noticeably,the alloys almost display the same hydrogen absorption kinetics when the temperature exceeds 513 K,indicating the predominant role of temperature among many factors acting on the rate of hydrogen absorption.As we all know,there are three steps happening in the hydrogenation procedure of magnesium [80],namely a)the dissociation of superficial H2molecules,b)the diffusing of H atoms through grain boundaries,c)the combination of H atoms and Mg atoms to form MgH2on the Mg/catalyst interfaces.Because the hydrogen dissociation needs quite high energy,it is reckoned as the rate-controlling step[81].As confirmed by Sakintuna et al.[35],the additives transition metals or their oxides in magnesium can act as the catalysts to decrease the dissociation energy.Liu et al.[82]considered that theoretically,the substitution atoms weakened the stability of Mg-H bond owing to the interaction between the valence electrons of H and the unsaturated d/f electron shell of the transition metals or oxides,improving the hydrogen absorption performance.Agarwal et al.[47]reported that it is difficult to refine the grains of Mg by mechanical milling due to the inevitable agglomeration of particles.The additive of brittle oxides or intermetallics provides convenience for reducing the particle size of Mg.The refined particles ameliorate the hydrogen absorbing and releasing properties due to the decreased diffusion length and the larger reactive surfaces of H2caused by particle refinement [83].Compared with the experimental alloy,TiO2and La2O3have higher hardness.Therefore,the existence of TiO2or La2O3nanoparticles increases the brittleness of alloy and eventually makes the equilibrium between fragmentation and agglomeration change to a reduced particle size,as stated by Rafi-ud-din et al.[84].The shorter diffusion channels for H atoms and larger specific surface for H2dissociation caused by particle refinement facilitate to enhance hydrogen absorption kinetics [85].

    Fig.5.SEM morphologies of the as-milled(20 h)La7Sm3Mg80Ni10–5M(M=none,TiO2,La2O3)alloys:(a)as-milled M=TiO2 alloy,(b)After six hydriding and dehydriding cycles M=none alloy,(c)After six hydriding and dehydriding cycles M=TiO2 alloy,(d)After six hydriding and dehydriding cycles M=La2O3 alloy.

    Fig.6.Isothermal hydrogen absorption kinetic curves of the as-milled(20 h)La7Sm3Mg80Ni10?5M(M=none,TiO2,La2O3)alloys at different temperatures:(a) M=none;(b) M=TiO2;(c) M=La2O3;(d)Rietveld refinement on the XRD patterns of the as-milled La7Sm3Mg80Ni10 alloy hydrogenated at 3MPa and 633K.

    To investigate the hydrogenation degree of the alloys and the phase structure changes during hydrogenation,The Rietveld refinements of the XRD patterns of as-milled La7Sm3Mg80Ni10alloy hydrogenated at 3 MPa and 593 K are provided,as illustrated in Fig.6(d).The milled alloys cannot be fitted with the Rietveld method because their XRD detections are amorphous.After hydrogenation,the amorphous phase is completely crystallized.Thus,the Rietveld method can be used analyzed the evolution of the phases of the asmilled alloy after hydrogenation.The result reveals that the as-milled hydrogenated alloy is composed of four hydrides,viz.MgH2,Mg2NiH4,LaH3and Sm3H7and the relative content of each phase is 56.1,25.9,11.5 and 6.5%,respectively.It suggested that the alloy is in saturated hydrogenation state.

    3.3.Hydrogen desorption kinetics and hydrogen desorption activation energy

    In order to research the relationships between the catalysts TiO2,La2O3and hydride stability,the temperature programmed desorption and DSC of the as-milled La7Sm3Mg80Ni10–5 M(M=none,TiO2,La2O3)after complete hydrogen absorption were measured at the heating rate of 5 K/min,as presented in Fig.7.It is observed that the adding catalyst renders an obvious effect on the hydride stability.The onset dehydrogenation temperature of La7Sm3Mg80Ni10–5 M(M=none,TiO2,La2O3)composites after hydrogenation is 547.4,537.2 and 540.1 K,severally.The temperatures of endothermic peaks in DSC curves of the alloys are 553.2,546.4 and 548.9 K,respectively.The change of initial hydrogen desorption temperature can reflect the hydrides stability.Evidently,the stability of the hydrogenated La7Sm3Mg80Ni10–5 M(M=none,TiO2,La2O3)composites is the following orderM=none>M=La2O3>M=TiO2.The decreased stability is attributed to the decline in Mg-H bond energy.The additives transition metals [86]or rare-earth elements can reduce the Mg-H bond energy and act as the catalysts owing to the electronic exchange reaction between these catalysts and magnesium hydride [57].Pighin et al.[87]investigated the function of various catalysts on Mg-H bond energy and found that the addition of transition metals or their oxides effectively decreased this bond energy,the magnesium hydride stability and the dehydrogenation temperature.According to the research of Abdellaoui et al.[12],the existence of new bonds weakens the bond strength between Mg and H,which can help us to understand the system instability and the decreased hydrogen desorption temperature mentioned above.

    Fig.7.Temperature programmed desorption curve of the as-milled(20 h)La7Sm3Mg80Ni10?5M(M=none,TiO2,La2O3)alloys after hydrogen absorption at a heating rate of 5K/min.

    Fig.8.Isothermal hydrogen desorption kinetic curves of the as-milled(20 h)La7Sm3Mg80Ni10?5M(M=none,TiO2,La2O3)alloys at different temperatures:(a) M=none;(b) M=TiO2;(c) M=La2O3,(d)Rietveld refinement on the XRD patterns of the as-milled La7Sm3Mg80Ni10 alloy dehydrogenated at 1×10–4 MPa and 633 K.

    Aiming at studying the hydrogen desorption kinetics of the composites with additives TiO2and La2O3,the plots of capacity versus time were tested at 553,573,593,613 and 633 K and presented in Fig.8.As can be observed,the reaction temperature greatly affects the hydrogen desorption kinetics of experimental alloys.Under the high-temperature conditions,all the alloys have a very fast reaction rate.In addition,it is noted that the adding catalysts TiO2and La2O3generates a favorable impact on the isothermal dehydrogenation kinetics.For further making sense of the influence of adding TiO2and La2O3on the kinetics,the time inquired by releasing 3 wt.%hydrogen is regarded as a reference standard.As shown in Fig.8,the time required by releasing 3 wt.% hydrogen at the temperatures of 553,573,593,613 and 633 K is 988,553,419,227 and 152 s for theM=none alloy,and 578,352,286,188,and 112 s for theM=TiO2alloy,and 594,366,301,197 and 132 s for theM=La2O3specimen,respectively.Obviously,the dehydrogenation rate is in the orderM=TiO2>M=La2O3>M=none.Based on the above data,the relationship between the time needed by desorbing 3 wt.% H2and temperature can be constructed,as displayed in Fig.8(d).It indicates that the additives TiO2and La2O3markedly ameliorate the hydrogen absorption kinetics,but the positive contribution decreases rapidly with the increase of hydrogen desorption temperature,which suggests that among all the factors affecting the dehydrogenation kinetics of alloys,temperature is predominate.It has come to light that the hydrogen desorption of MgH2is completed through three stages:(a)magnesium phase nucleates and grows,(b)hydrogen diffuses from the magnesium hydride matrix to the surfaces,and(c)two adjacent hydrogen atoms combine to form hydrogen molecules[88].The improved hydrogen desorption kinetics of magnesium based alloys is most likely attributed to the decline in hydride stability,the diminution of the particles size and the increase of the defect density on the particle surface of the alloys caused by adding catalysts TiO2and La2O3.As mentioned above,the rare-earth elements and transition metals or their oxides can reduce the bond energy of Mg-H,thus,weaken the magnesium hydrides stability and accelerate the hydrides decomposition [47].High hardness catalyst particles are likely to be embedded into the interior of alloy particles under the action of repeated impact stress in the process of ball grinding,so the alloy particles broken and particle size greatly reduced,as considered by Jain et al.[89].Meanwhile,the additive high-hardness catalyst induces the surface defects and brings on the particle refinement of alloy in the procedure of mechanical milling [84].The induced defects are favorable to the nucleation and increase the surface reactivity.Besides,the decreased particle size makes the length of hydrogen diffusion channels shorten.These factors definitely accelerate the hydrogenation and dehydration of the Mg-based materials[90].Nevertheless,the improved kinetics because of adding TiO2and La2O3particles is not only due to the decline in particles size.Other factors such as the nature of the added oxides and the local electronic structure and the reduction of the oxide during heating should also be considered.Partially reduced oxides are expected to have different valence states and may act on ameliorating the hydrogen desorption property.Generally speaking,the main catalysis of transition metal-based catalysts is engendered by the transition metal ions and their ability to form hydrogen bonds.In this way,transition metal-based catalysts provide a faster route for H atoms diffusion.It has been reported that oxygen vacancies(also known as anoxic surfaces)on oxide surfaces also have catalytic activity [84].Hence,we believe that the hydrogen desorption results from thermodynamically induced surface vacancies.In summary,the improved kinetics may be due to the uniform dispersion of these anoxic oxide particles,which shorten the diffusion path between reaction ions.The oxygen vacancies act as the sites for nucleation and growth of dehydrogenation products,and promote the dehydrogenation process.

    To investigate the dehydrogenation degree of the hydrides and the phase structure changes during dehydrogenation,The Rietveld refinements of the XRD patterns of as-milled saturated hydrides dehydrogenated at 1×10?4MPa and 633 K are given,as illustrated in Fig.8(d).It reveals that the phase Mg,Mg2Ni,LaH3and Sm3H7exist in the dehydrogenated alloy.It is very clear that rare earth hydrides LaH3and Sm3H7remain undecomposed at experimental temperatures and pressures.The relative content of each phase in the alloy is 55.7,26.1,11.6 and 6.6%,respectively.

    Fig.9.JMA plots and Arrhenius plots for the dehydrogenation of the as-milled(20 h)La7Sm3Mg80Ni10–5M(M=none,TiO2,La2O3)alloys at different temperatures:(a) M=none;(b) M=TiO2;(c) M=La2O3,(d)The hydrogen desorption activation energy(E)of the experimental alloys.

    Generally,the occurrence of gas-solid reaction needs to overcome a total energy barrier that can be reflected in terms of the apparent activation energy.Hence,when the activation energy reaches a certain requirement,the reaction can take place smoothly.The apparent activation energy of the hydrogenated La7Sm3Mg80Ni10–5 M(M=none,TiO2,La2O3)alloys in hydrogen desorption is evaluated by the Arrhenius and Kissinger methods.As well known,the nucleation and growth of dehydrogenation products are the crucial factors that control the hydrogen desorption reaction of magnesium based materials [91].In general,the Johnson-Mehl-Avrami-Kolmogorov(JMAK)model can simulate this solid-state reaction [92]:

    a—the reaction fraction;

    n—the avrami index;

    k—the dehydrogenation rate constant;

    t—the reaction time.

    According to Fig.8,the fitting curves ln [-ln(1-α)]vs.lntat 573,593,613 and 633 K can be ploted,as provided in Fig.9.As can be observed,the JMAK sketches are almost linear,suggesting the dehydrogenation of the composite is composed of two steps,including the first stage instantaneous nucleation and the second stage 3D growth controlled by interface [93].Theηandηlnkvalues were obtained according to the slope and intercept in fitting curves at the corresponding temperature.Thus the value ofkcan be acquired.The apparent activation energy()of dehydrogenation reaction was estimated gained by using Arrhenius formula [57]:

    A—a temperature independent coefficient;

    R—the gas constant(8.3145 J/mol/K);

    T—the absolute temperature of reaction;

    k—the dehydrogenation rate constant.

    The Arrhenius plots of lnkvs.1/Tof the alloys are sketched,as presented in Fig.9.The apparent activation energyof the as-milled alloys was acquired from the slopes of the Arrhenius plots.68.1,62.1 and 63.6 kJ/mol correspond to the La7Sm3Mg80Ni10–5 M(M=none,TiO2,La2O3)alloys,severally.Obviously,the apparent activation energy of the dehydrogenation of the alloys is in sequencedM=TiO2

    Fig.10.DSC curves and Kissinger plots of the as-milled(20 h)La7Sm3Mg80Ni10?5M(M=none,TiO2,La2O3)alloys at different temperatures:(a)M=none;(b) M=TiO2;(c) M=La2O3,(d)The hydrogen desorption activation energy(E)of the experimental alloys.

    In order to compare with the JMAK model,the Kissinger method is employed to estimate the activation energy,as following equation [94]:

    β—the heating rate;

    TP—the absolute temperature corresponding to the maximal desorption rate,

    R—the gas constant(8.3145 J/mol/K).

    The as-milled La7Sm3Mg80Ni10–5 M(M=none,TiO2,La2O3)alloys need to absorb H2to the saturated state before the measure of DSC.Fig.10 demonstrates the non-isothermal dehydrogenation curves tested at the heating rates of 5,10,15 and 20 K/min,severally.As we can see,an endothermic peak exists in each DSC curve,suggesting the same reaction procedure of each specimen.According to Fig.10,the graphs of ln(β/)vs.1/TPcan be sketched,as presented in Fig.10.It is noted that ln(β/)vs.1/TPplot is almost linear,so from the slopes of it,the activation energywas obtained.According to the calculation,apparent activation energies of the La7Sm3Mg80Ni10–5 M(M=none,TiO2,La2O3)alloys are 65.5,60.2 and 61.7 kJ/mol and in order ofM=TiO2

    3.4.P-C-T curves and hydrogen storage thermodynamics

    Reducing the thermal stability of magnesium hydride is the main goal to improve its hydrogen storage performance and realize its practical application.To inspect the effect of adding TiO2and La2O3on the thermodynamics,the P-C-T curves of the as-milled specimens were tested at the temperature of 593,613 and 633 K and given in Fig.11.Obviously,the pressure platform is quite flat and the hysteresis coefficient(Hf=ln(Pa/Pd))is small.The catalyst TiO2or La2O3has no evidently change in the platform characteristic reflected in the P-C-T curves of alloys.As we can observe,two pressure plateaus emerge in every P-C-T curve and the higher and lower platform pressures stand for the formation/dissociation of the Mg2NiH4and MgH2hydrides,severally [96,97].According to the plateau pressures(PaandPd)in Fig.11,the thermodynamics parameters enthalpy change?Hand entropy change?Sare evaluated by Van’t Hoff equation [98]:

    Fig.11.P-C-T curves and Van’t Hoff plots of the as-milled(20 h)La7Sm3Mg80Ni10?5M(M=none,TiO2,La2O3)alloys at different temperatures:(a)M=none;(b) M=TiO2;(c) M=La2O3,(d)The absolute value of hydriding and dehydriding enthalpy of the experimental alloys.

    PH2—the equilibrium hydrogen gas pressure corresponding to MgH2;

    P0—the standard atmospheric pressure;

    R—the gas constant(8.3145J/mol/K);

    T—the absolute temperature of reaction.

    The Van’t Hoff graphs of lnPH2/P0vs.1/Tfor the as-milled La7Sm3Mg80Ni10–5 M(M=none,TiO2,La2O3)composites can be sketched.Hence,the?Hand?Scan be calculated according to the slopes and intercepts in Van’t Hoff diagrams and listed in Table 1.It uncovers that the addition of TiO2and La2O3has not notably impact on the improvement of the experimental materials’thermodynamics and the reduction of corresponding hydrides stability.The addition of catalysts TiO2or La2O3decreases the stability magnesium hydride.A similar result also emerged in the investigations of Anik et al.[99]and Bououdina et al.[39].Obviously,the absolute values of dehydrogenation enthalpy change?Hdeof the alloys are in following orderM=none>M=La2O3>M=TiO2.Based on the above results,we can find that both isothermal and non-isothermal analyses reveal that the catalysts TiO2and La2O3weaken the magnesium hydride stability and improve the hydrogen absorption and desorption kinetics.The positive contribution to the hydrogen storing thermodynamic and dynamics of the specimens caused by two catalysts is in following orderM=TiO2>M=La2O3>M=none.

    Table 1.Enthalpy ?H and entropy ?S of the as-milled alloys.

    4.Conclusion

    (1)The addition of TiO2and La2O3has no change in the phase composition but reduces the agglomeration tendency of particles in the process of mechanical milling and make the particle size of the as-milled alloy markedly decreased.It is this modification of the microstructure that remarkably enhances the hydrogen absorption and desorption performances.

    (2)The addition of TiO2and La2O3have obviously positive contribution to the hydrogenation and dehydrogenation kinetics of the experimental alloys,which is in orderM=TiO2>M=La2O3>M=none.It is ascribed to the decline in the size of grains and particles,the generation of the fresh surface and the creation of the various crystal defects derived from ball milling and adding catalysts.

    (3)The addition of TiO2and La2O3catalysts has a slightly favorable influence on the improvement of the thermodynamics of alloy and the stability of the hydride,which is in sequenceM=TiO2>M=La2O3>M=none.

    Declaration of Competing Interest

    We declare that we have no financial and personal relationships with other people or organizations that can inappropriately influence our work,there is no professional or other personal interest of any nature or kind in any product,service and/or company that could be construed as influencing the position presented in,or the review of the manuscript entitled.

    Acknowledgments

    This study was financially supported by the National Natural Science Foundation of China(Nos.51901105,51871125,and 51761032),Natural Science Foundation of Inner Mongolia,China(2019BS05005),Inner Mongolia University of Science and Technology Innovation Fund(2019QDL-B11)and Major Science and Technology Innovation Projects in Shandong Province(2019JZZY010320).

    欧美 日韩 精品 国产| 国产又色又爽无遮挡免| 亚洲国产av影院在线观看| 男女无遮挡免费网站观看| tube8黄色片| 99热6这里只有精品| 久久国内精品自在自线图片| 免费观看无遮挡的男女| 国模一区二区三区四区视频| 国产高清有码在线观看视频| 免费播放大片免费观看视频在线观看| 一级,二级,三级黄色视频| av.在线天堂| 国产熟女欧美一区二区| 亚洲欧美精品自产自拍| 国产av精品麻豆| av黄色大香蕉| 亚洲五月色婷婷综合| 老熟女久久久| 亚洲美女搞黄在线观看| 免费人成在线观看视频色| 另类精品久久| 黄片播放在线免费| 国产 一区精品| 亚洲人成网站在线观看播放| 国产视频内射| 亚洲欧美日韩卡通动漫| 欧美bdsm另类| 精品熟女少妇av免费看| 免费观看在线日韩| 国产精品成人在线| 精品久久久久久电影网| 久久精品国产自在天天线| 日韩成人伦理影院| 亚洲精品国产av蜜桃| 夜夜骑夜夜射夜夜干| 黑人猛操日本美女一级片| 亚洲色图 男人天堂 中文字幕 | 国产成人a∨麻豆精品| 99久久中文字幕三级久久日本| 午夜福利网站1000一区二区三区| 亚洲欧美成人精品一区二区| 午夜精品国产一区二区电影| 国产精品久久久久成人av| 国产国拍精品亚洲av在线观看| 亚洲美女黄色视频免费看| 久久久久久久久久久免费av| 欧美97在线视频| 欧美日韩亚洲高清精品| 夜夜骑夜夜射夜夜干| 色哟哟·www| 九草在线视频观看| 夫妻性生交免费视频一级片| 欧美日韩亚洲高清精品| 三级国产精品欧美在线观看| 黄色一级大片看看| av有码第一页| 国产成人精品在线电影| 国产精品一区二区在线观看99| 人人妻人人澡人人爽人人夜夜| 99国产综合亚洲精品| 人妻夜夜爽99麻豆av| 日产精品乱码卡一卡2卡三| 免费播放大片免费观看视频在线观看| 亚州av有码| 国产高清国产精品国产三级| 一区二区三区四区激情视频| 亚洲精品国产av成人精品| 熟女人妻精品中文字幕| 中文字幕亚洲精品专区| 国产日韩欧美在线精品| 尾随美女入室| .国产精品久久| 亚洲在久久综合| 另类亚洲欧美激情| 91成人精品电影| 另类精品久久| 女的被弄到高潮叫床怎么办| 午夜福利影视在线免费观看| 一级毛片我不卡| 少妇猛男粗大的猛烈进出视频| 免费久久久久久久精品成人欧美视频 | 91精品三级在线观看| 国产毛片在线视频| 99久久精品一区二区三区| 少妇精品久久久久久久| 91精品一卡2卡3卡4卡| a级毛片在线看网站| 亚洲美女视频黄频| 色视频在线一区二区三区| 亚洲综合精品二区| 亚洲精品日本国产第一区| 精品一区二区三区视频在线| 丝袜在线中文字幕| 久久精品熟女亚洲av麻豆精品| 99国产精品免费福利视频| 全区人妻精品视频| av在线观看视频网站免费| 亚洲av男天堂| 青春草亚洲视频在线观看| 一级毛片我不卡| 日本黄色片子视频| 日本91视频免费播放| 蜜桃在线观看..| 全区人妻精品视频| 国产黄频视频在线观看| 精品久久国产蜜桃| 亚洲国产最新在线播放| 久久狼人影院| 久久99蜜桃精品久久| 国产黄片视频在线免费观看| 乱码一卡2卡4卡精品| 亚洲国产精品一区二区三区在线| 日韩强制内射视频| 99热这里只有是精品在线观看| 丝袜喷水一区| 婷婷成人精品国产| 亚洲国产色片| 精品亚洲成a人片在线观看| 免费av不卡在线播放| 一区在线观看完整版| 有码 亚洲区| 性色avwww在线观看| 18禁观看日本| 不卡视频在线观看欧美| 日本色播在线视频| 亚洲精品av麻豆狂野| 亚洲av综合色区一区| 久久久精品免费免费高清| 制服丝袜香蕉在线| 成人黄色视频免费在线看| 有码 亚洲区| 免费观看无遮挡的男女| 国产精品一区www在线观看| 新久久久久国产一级毛片| 麻豆成人av视频| 男的添女的下面高潮视频| 九九久久精品国产亚洲av麻豆| 中国国产av一级| 亚洲国产色片| 在线观看人妻少妇| 人妻制服诱惑在线中文字幕| 亚洲精品视频女| 国产男人的电影天堂91| 毛片一级片免费看久久久久| 国产乱来视频区| 亚洲在久久综合| 少妇人妻久久综合中文| 午夜激情福利司机影院| 亚洲av二区三区四区| 日韩中文字幕视频在线看片| 国产精品久久久久久精品古装| 成人亚洲精品一区在线观看| √禁漫天堂资源中文www| 亚洲成人一二三区av| 岛国毛片在线播放| 免费av中文字幕在线| 嘟嘟电影网在线观看| 美女脱内裤让男人舔精品视频| 最近最新中文字幕免费大全7| 免费大片18禁| 久久久久久久大尺度免费视频| 少妇的逼水好多| 91久久精品电影网| 亚洲欧美日韩卡通动漫| 成人午夜精彩视频在线观看| 日韩中字成人| 下体分泌物呈黄色| 国产精品.久久久| 蜜臀久久99精品久久宅男| 久久97久久精品| 久久国产精品男人的天堂亚洲 | 亚洲熟女精品中文字幕| 在线 av 中文字幕| 精品酒店卫生间| 精品久久蜜臀av无| 午夜精品国产一区二区电影| 亚洲精品视频女| 毛片一级片免费看久久久久| a级片在线免费高清观看视频| 亚洲精品日韩av片在线观看| 国产 精品1| 国产成人精品福利久久| 欧美少妇被猛烈插入视频| 下体分泌物呈黄色| 国产在线免费精品| 久久ye,这里只有精品| 黑人巨大精品欧美一区二区蜜桃 | 久久人人爽人人爽人人片va| 精品99又大又爽又粗少妇毛片| 在线看a的网站| 国产综合精华液| 国产日韩欧美亚洲二区| 男人操女人黄网站| 国产精品久久久久久精品电影小说| 国产精品99久久久久久久久| 久久人人爽人人爽人人片va| 亚洲国产精品国产精品| 国产精品无大码| 我要看黄色一级片免费的| a级片在线免费高清观看视频| 99热网站在线观看| 久热久热在线精品观看| 久久久亚洲精品成人影院| 久久韩国三级中文字幕| 国产成人一区二区在线| 免费久久久久久久精品成人欧美视频 | 色婷婷久久久亚洲欧美| 女性生殖器流出的白浆| 多毛熟女@视频| 九九爱精品视频在线观看| 国产免费一级a男人的天堂| 亚洲色图综合在线观看| kizo精华| 一级片'在线观看视频| 亚洲国产成人一精品久久久| 这个男人来自地球电影免费观看 | 男的添女的下面高潮视频| 91久久精品国产一区二区成人| 夜夜爽夜夜爽视频| 精品一区二区免费观看| 久久久久久久大尺度免费视频| 人成视频在线观看免费观看| 成人无遮挡网站| 欧美老熟妇乱子伦牲交| 高清av免费在线| 自线自在国产av| 亚洲av中文av极速乱| 最后的刺客免费高清国语| 18+在线观看网站| 久久久久久久久大av| 国产老妇伦熟女老妇高清| 中文字幕亚洲精品专区| 亚洲精品av麻豆狂野| 妹子高潮喷水视频| 黄色一级大片看看| 久久人妻熟女aⅴ| 大话2 男鬼变身卡| 国产亚洲最大av| 美女内射精品一级片tv| 黄色视频在线播放观看不卡| 中文字幕亚洲精品专区| 国产成人一区二区在线| 久久久亚洲精品成人影院| 久久精品久久精品一区二区三区| 婷婷成人精品国产| 国产成人精品福利久久| 国产一区二区在线观看av| 亚洲,一卡二卡三卡| 丁香六月天网| 日韩av在线免费看完整版不卡| 免费大片18禁| 成年av动漫网址| 亚洲av日韩在线播放| 欧美国产精品一级二级三级| 欧美亚洲日本最大视频资源| av.在线天堂| av播播在线观看一区| 中文字幕av电影在线播放| 欧美精品高潮呻吟av久久| 22中文网久久字幕| 一级a做视频免费观看| 国产在视频线精品| 精品亚洲乱码少妇综合久久| 天堂中文最新版在线下载| 九色成人免费人妻av| 亚洲精品456在线播放app| 人体艺术视频欧美日本| 少妇的逼好多水| 国产黄色视频一区二区在线观看| 欧美人与善性xxx| 菩萨蛮人人尽说江南好唐韦庄| av电影中文网址| 日韩视频在线欧美| 亚洲精品日韩在线中文字幕| 亚洲精品av麻豆狂野| 极品少妇高潮喷水抽搐| 男女免费视频国产| 亚洲激情五月婷婷啪啪| 久久国产亚洲av麻豆专区| videos熟女内射| 寂寞人妻少妇视频99o| 国产精品欧美亚洲77777| 三上悠亚av全集在线观看| 亚洲成色77777| 大码成人一级视频| 午夜激情福利司机影院| 80岁老熟妇乱子伦牲交| 3wmmmm亚洲av在线观看| 亚洲精品av麻豆狂野| videos熟女内射| 一级,二级,三级黄色视频| 日日爽夜夜爽网站| 涩涩av久久男人的天堂| 男女啪啪激烈高潮av片| 亚洲国产欧美在线一区| 婷婷色av中文字幕| 女人久久www免费人成看片| 啦啦啦在线观看免费高清www| 午夜激情福利司机影院| 亚洲国产日韩一区二区| 26uuu在线亚洲综合色| av不卡在线播放| 中国美白少妇内射xxxbb| 观看av在线不卡| 嫩草影院入口| 极品少妇高潮喷水抽搐| 美女脱内裤让男人舔精品视频| 日韩中文字幕视频在线看片| 国产乱来视频区| 九草在线视频观看| 少妇的逼好多水| 晚上一个人看的免费电影| 伊人亚洲综合成人网| 亚洲怡红院男人天堂| 免费av中文字幕在线| 丝袜脚勾引网站| 国产精品一区www在线观看| 亚洲精品久久成人aⅴ小说 | 国产精品偷伦视频观看了| 一二三四中文在线观看免费高清| 欧美成人午夜免费资源| 亚洲第一区二区三区不卡| 精品少妇久久久久久888优播| 91成人精品电影| 久久久亚洲精品成人影院| 日韩成人av中文字幕在线观看| 亚洲美女搞黄在线观看| 大话2 男鬼变身卡| 99九九线精品视频在线观看视频| 91成人精品电影| 一区二区三区精品91| 一区二区三区免费毛片| 18+在线观看网站| 国产精品秋霞免费鲁丝片| av女优亚洲男人天堂| 免费看光身美女| 精品久久蜜臀av无| 成人二区视频| .国产精品久久| 草草在线视频免费看| 老司机影院成人| 97超碰精品成人国产| 亚洲第一区二区三区不卡| 成人影院久久| 欧美成人精品欧美一级黄| 肉色欧美久久久久久久蜜桃| 亚洲综合色网址| 如何舔出高潮| 午夜福利在线观看免费完整高清在| 高清av免费在线| 精品99又大又爽又粗少妇毛片| 我的老师免费观看完整版| 欧美日韩视频高清一区二区三区二| 91成人精品电影| 免费人妻精品一区二区三区视频| 丝袜美足系列| 80岁老熟妇乱子伦牲交| 视频区图区小说| 国产黄色视频一区二区在线观看| 国产成人精品福利久久| 一级片'在线观看视频| 中文字幕免费在线视频6| 你懂的网址亚洲精品在线观看| 国产精品女同一区二区软件| 免费黄频网站在线观看国产| 亚洲av中文av极速乱| 一本大道久久a久久精品| 色网站视频免费| videossex国产| 色网站视频免费| 国模一区二区三区四区视频| 国产有黄有色有爽视频| 老司机亚洲免费影院| 亚洲国产精品一区二区三区在线| 丝袜脚勾引网站| 国产永久视频网站| 最近的中文字幕免费完整| 亚洲激情五月婷婷啪啪| 青春草国产在线视频| 在线观看免费日韩欧美大片 | 满18在线观看网站| 午夜激情av网站| 久久免费观看电影| 国产视频内射| 婷婷色综合大香蕉| 免费黄网站久久成人精品| 亚洲欧洲国产日韩| 久久久久久久久久久丰满| 菩萨蛮人人尽说江南好唐韦庄| 乱人伦中国视频| 国产伦精品一区二区三区视频9| 欧美成人午夜免费资源| 哪个播放器可以免费观看大片| 啦啦啦啦在线视频资源| 看非洲黑人一级黄片| 中文字幕亚洲精品专区| 精品久久久久久电影网| 高清在线视频一区二区三区| 大码成人一级视频| 久久精品熟女亚洲av麻豆精品| 国产男女内射视频| a级片在线免费高清观看视频| 丰满迷人的少妇在线观看| 亚洲无线观看免费| 久久精品久久精品一区二区三区| 国产精品一区二区在线观看99| 好男人视频免费观看在线| 成人免费观看视频高清| 丰满少妇做爰视频| 我的老师免费观看完整版| 欧美精品亚洲一区二区| 欧美亚洲 丝袜 人妻 在线| 欧美3d第一页| 日韩不卡一区二区三区视频在线| tube8黄色片| 97超视频在线观看视频| 两个人免费观看高清视频| 男女高潮啪啪啪动态图| 午夜福利,免费看| 观看美女的网站| 高清视频免费观看一区二区| 国产成人免费无遮挡视频| 久久久久久久久久人人人人人人| 我要看黄色一级片免费的| 中文天堂在线官网| 丰满饥渴人妻一区二区三| 国产精品无大码| 纯流量卡能插随身wifi吗| 精品一区二区免费观看| 久久国产精品男人的天堂亚洲 | 熟女人妻精品中文字幕| freevideosex欧美| 晚上一个人看的免费电影| 亚洲av成人精品一二三区| 乱码一卡2卡4卡精品| 又粗又硬又长又爽又黄的视频| 两个人免费观看高清视频| 9色porny在线观看| 女性生殖器流出的白浆| 大香蕉97超碰在线| 亚洲激情五月婷婷啪啪| 18禁动态无遮挡网站| 极品人妻少妇av视频| 国产日韩欧美在线精品| 亚洲欧洲精品一区二区精品久久久 | av视频免费观看在线观看| 各种免费的搞黄视频| 中文字幕免费在线视频6| 只有这里有精品99| 日韩 亚洲 欧美在线| 边亲边吃奶的免费视频| 午夜91福利影院| 国产精品久久久久久久电影| 日韩 亚洲 欧美在线| 欧美日韩成人在线一区二区| 自线自在国产av| av免费在线看不卡| 久久毛片免费看一区二区三区| 久久人人爽av亚洲精品天堂| 人人妻人人澡人人爽人人夜夜| 我要看黄色一级片免费的| 母亲3免费完整高清在线观看 | 黑人猛操日本美女一级片| 亚洲,欧美,日韩| 国产成人a∨麻豆精品| 日本-黄色视频高清免费观看| 黄色视频在线播放观看不卡| 久久久久久久久久成人| 国产精品人妻久久久影院| 如何舔出高潮| 日韩一区二区三区影片| 满18在线观看网站| 18禁在线无遮挡免费观看视频| 欧美激情国产日韩精品一区| 国产高清国产精品国产三级| 18+在线观看网站| 精品一区在线观看国产| 在线观看一区二区三区激情| 最新中文字幕久久久久| 美女视频免费永久观看网站| 亚洲精品自拍成人| 人妻少妇偷人精品九色| 日本黄大片高清| 高清欧美精品videossex| 国产高清不卡午夜福利| 嫩草影院入口| av有码第一页| 午夜免费鲁丝| 99re6热这里在线精品视频| 成人漫画全彩无遮挡| 欧美老熟妇乱子伦牲交| 狂野欧美激情性xxxx在线观看| 国产av一区二区精品久久| 亚洲欧洲日产国产| 国产淫语在线视频| 国产国拍精品亚洲av在线观看| 国产成人精品一,二区| 中国三级夫妇交换| 乱码一卡2卡4卡精品| 色吧在线观看| 国产日韩欧美亚洲二区| 亚洲伊人久久精品综合| 精品亚洲乱码少妇综合久久| 亚洲熟女精品中文字幕| 婷婷色av中文字幕| 亚洲精品,欧美精品| 国产精品99久久99久久久不卡 | 亚洲av日韩在线播放| 老司机亚洲免费影院| 午夜福利在线观看免费完整高清在| 婷婷色综合www| 肉色欧美久久久久久久蜜桃| 亚洲婷婷狠狠爱综合网| 成人亚洲欧美一区二区av| 一区二区三区免费毛片| 久久久久久久久久人人人人人人| 热re99久久国产66热| 日韩熟女老妇一区二区性免费视频| 亚洲高清免费不卡视频| 久久热精品热| 高清毛片免费看| 国产精品欧美亚洲77777| 九草在线视频观看| a级毛片黄视频| 亚洲色图 男人天堂 中文字幕 | 亚洲精品一二三| 亚洲婷婷狠狠爱综合网| 这个男人来自地球电影免费观看 | 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久久久久成人| 夜夜骑夜夜射夜夜干| av播播在线观看一区| 老女人水多毛片| 国产精品女同一区二区软件| 成年美女黄网站色视频大全免费 | 国产成人免费观看mmmm| 日韩在线高清观看一区二区三区| 一级毛片我不卡| 欧美亚洲 丝袜 人妻 在线| 精品亚洲成a人片在线观看| 亚洲av电影在线观看一区二区三区| 成人亚洲精品一区在线观看| 亚洲五月色婷婷综合| 极品人妻少妇av视频| 97超碰精品成人国产| 丝袜美足系列| 国产黄色免费在线视频| 丁香六月天网| 国产成人av激情在线播放 | 久久韩国三级中文字幕| 男女边摸边吃奶| 国产乱人偷精品视频| 搡老乐熟女国产| 国产精品久久久久久久电影| 久久久久国产精品人妻一区二区| 国产午夜精品久久久久久一区二区三区| 91精品国产九色| 国产免费又黄又爽又色| 亚洲综合色惰| 丝瓜视频免费看黄片| 亚洲美女视频黄频| 十八禁高潮呻吟视频| 国产欧美另类精品又又久久亚洲欧美| 在线观看人妻少妇| 免费不卡的大黄色大毛片视频在线观看| 男人添女人高潮全过程视频| 国产精品三级大全| 欧美激情国产日韩精品一区| 国产免费视频播放在线视频| 永久网站在线| 一本大道久久a久久精品| 日韩av不卡免费在线播放| 国产精品蜜桃在线观看| 国产熟女午夜一区二区三区 | 免费大片黄手机在线观看| 综合色丁香网| 免费av不卡在线播放| 曰老女人黄片| 在线观看免费日韩欧美大片 | av播播在线观看一区| 九九爱精品视频在线观看| 97超视频在线观看视频| 成人综合一区亚洲| 狂野欧美激情性xxxx在线观看| 91精品三级在线观看| 久久国产亚洲av麻豆专区| 免费观看av网站的网址| 国产精品一区www在线观看| 免费高清在线观看视频在线观看| 性色av一级| 国产熟女午夜一区二区三区 | 亚洲色图 男人天堂 中文字幕 | 欧美日韩av久久| freevideosex欧美| 美女视频免费永久观看网站| 国产一区二区在线观看av| 国产淫语在线视频| 成人影院久久| 亚洲中文av在线| 免费观看的影片在线观看| 18在线观看网站| av线在线观看网站| 又粗又硬又长又爽又黄的视频| 日日摸夜夜添夜夜添av毛片| 久久 成人 亚洲| 你懂的网址亚洲精品在线观看| 99久久精品国产国产毛片| 欧美精品人与动牲交sv欧美| 啦啦啦啦在线视频资源| 欧美日韩视频精品一区| 中文字幕制服av| 国产成人91sexporn| av在线app专区| 国产精品欧美亚洲77777|