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

    Deflagration characteristics of freely propagating flames in magnesium hydride dust clouds

    2024-02-29 08:23:44QiweiZhngYngfnChengBeieiZhngDnyiLiZhowuShen
    Defence Technology 2024年1期

    Qiwei Zhng , Yngfn Cheng ,,*, Beiei Zhng , Dnyi Li ,c, Zhowu Shen

    a School of Chemical Engineering, Anhui University of Science and Technology, Huainan 232001, Anhui, PR China

    b School of Civil Engineering and Architecture, Anhui University of Science and Technology, Huainan 232001, Anhui, PR China

    c Anhui Engineering Laboratory of Explosive Materials and Technology, Anhui University of Science and Technology, Huainan 232001, Anhui, PR China

    d School of Engineering Science, University of Science and Technology of China, Hefei 230027, Anhui, PR China

    Keywords: Magnesium hydride dust Flame combustion mechanism Particle size Dust explosion Two-color pyrometer

    ABSTRACT The flame propagation processes of MgH2 dust clouds with four different particle sizes were recorded by a high-speed camera.The dynamic flame temperature distributions of MgH2 dust clouds were reconstructed by the two-color pyrometer technique, and the chemical composition of solid combustion residues were analyzed.The experimental results showed that the average flame propagation velocities of 23 μm, 40 μm, 60 μm and 103 μm MgH2 dust clouds in the stable propagation stage were 3.7 m/s,2.8 m/s, 2.1 m/s and 0.9 m/s, respectively.The dust clouds with smaller particle sizes had faster flame propagation velocity and stronger oscillation intensity, and their flame temperature distributions were more even and the temperature gradients were smaller.The flame structures of MgH2 dust clouds were significantly affected by the particle sinking velocity, and the combustion processes were accompanied by micro-explosion of particles.The falling velocities of 23 μm and 40 μm MgH2 particles were 2.24 cm/s and 6.71 cm/s,respectively.While the falling velocities of 60 μm and 103 μm MgH2 particles were as high as 15.07 cm/s and 44.42 cm/s,respectively,leading to a more rapid downward development and irregular shape of the flame.Furthermore, the dehydrogenation reaction had a significant effect on the combustion performance of MgH2 dust.The combustion of H2 enhanced the ignition and combustion characteristics of MgH2 dust, resulting in a much higher explosion power than the pure Mg dust.The micro-structure characteristics and combustion residues composition analysis of MgH2 dust indicated that the combustion control mechanism of MgH2 dust flame was mainly the heterogeneous reaction,which was affected by the dehydrogenation reaction.

    1.Introduction

    As a new type of weapon, thermobaric bombs have received extensive attention in recent years.The unique damage effect with pervasive high heat and shock wave produced by its explosion is incomparable to traditional ammunition, and it is especially suitable for killing the enemy in closed spaces such as caves, underground fortifications and buildings [1].Thermobaric bombs have the characteristics of both high explosives and fuel air explosives.Powders such as aluminum, titanium, magnesium, zirconium,silicon and boron are added to thermobaric bombs, which further enhances their thermal and pressure effects with a significant energy multiplication function [2-4].Thermobaric bombs are originated from early fuel air explosives, after more than 20 years of development, they have formed a series of products including not only the short-range thermobaric bombs for individual soldiers,but also the artillery-launched thermobaric grenades and thermobaric rockets, as well as the fighter-launched thermobaric bombs and thermobaric air-to-ground missiles.After continuous development and improvement, thermobaric bombs have been gradually expanding to the directions of miniaturization multi-carrier launching.In order to improve the energy density of thermobaric bombs, countries in the world have begun to look for novel combustion materials [5].Hydrogen storage alloys served as additives are widely used in energetic materials such as explosives,pyrotechnics and propellants because of their high energy density,combustion calorific value and energy release efficiency [6-9],among which MgH2is one of the most representative materials[10].Xue et al.[11]found that RDX with 10%MgH2powders had the optimal shock wave parameters,and compared with the reference sample, its peak overpressure, specific impulse and shock wave energy were increased by 5.7%, 7.0% and 8.4%, respectively.Cheng et al.[12,13] demonstrated that MgH2powders could significantly increase the shock wave overpressure,specific impulse,shock wave energy and bubble energy of emulsion explosive in the underwater explosion, and improve their anti-dynamic pressure ability in the delay blasting operation.Liu et al.[14]explored the effects of MgH2powders on the thermal decomposition, combustion rate and explosion heat of ammonium perchlorate-based composite solid propellant, and found that MgH2powders could reduce its peak temperature of thermal decomposition and increase its combustion rate.Chen et al.[15] studied the properties of novel hydrogen storage alloys dust clouds,and the result showed that MgH2in the dust mixture could effectively reduce its minimum ignition energy and increase its peak explosion pressure.Hydrogen storage alloy MgH2has the dual characteristics of metal fuel and hydrogen,so it also has a good application prospect in thermobaric bombs [16].

    Abbreviations mp Particle mass, kg dp Particle diameter, m g Gravitational acceleration, m/s2 t Time, s up u Particle velocity, m/s2 Fluid velocity, m/s2 ut Particle settling velocity, m/s v Fluid kinematic viscosity, m2/s ρp Particle density, kg/m3 or g/cm3 ρg Air density, kg/m3 or g/cm3 ζd Drag coefficient Re Reynolds number g gas p Particle

    Studying the dust explosion characteristics of hydrogen storage alloy MgH2is helpful to reveal the damage mechanism of thermobaric bombs containing MgH2powders.The flame structures and propagation characteristics of the dust clouds are extremely important parameters to characterize its detonation characteristics,and researchers have done a lot of works on them [17].Gao et al.[18] studied the effects of particle size distributions on the flame propagation characteristics of octadecanol dust clouds in an open space, and the experimental results showed that the flame front structures were changed with the varying particle size distributions, and the flame propagated in the dust cloud with a smaller particle size had a regular shape,which was similar to the premixed gas explosion.Chang et al.[19]observed the flame development of nano and micron aluminum dust clouds,and found that the flame propagation of nano dust clouds had higher temperature gradients in the preheating zone and lower minimum explosive concentrations as compared to micron aluminum dust clouds.Huang et al.[20]studied the combustion characteristics of aluminum dust with different particle size distributions in the laminar air flow, and found that with the particle size decreased from micron to nano scale,the flame propagating velocity increased and the combustion mechanism changed from a diffusion-controlled mode to a kinetically-controlled mode.Yu et al.[21] observed the flame propagation behaviors and micro-structures in micron and nano titanium dust explosions, and discovered that the flame of nanotitanium dust cloud propagated rapidly and had an obvious micro-explosion phenomenon.Jing et al.[22,23] studied the explosion mechanism of micron-sized flake aluminum powder, and found that the transient reaction rate of flake aluminum powder/air was mainly determined by oxygen diffusion rate and chemical kinetic rate.

    Despite the explosion parameters, flame temperature is also a very important parameter for evaluating the thermal damage effect of dust explosions.In recent years,with the continuous maturity of photoelectric technology, non-contact optical temperature measurement methods have been developed rapidly,among which the two-color pyrometer technique is the most popular method for it is no need to determine the optical properties(such as transmittance and emissivity)of the flame before the temperature measurement[24,25].In our previous works, we also carried out a lot of studies on the two-color pyrometer technique[26-28],which provided an experimental technical support for measuring the flame temperature distribution of MgH2dust cloud.

    In the study, flame dynamic evolution processes and temperature distribution in the micron scale MgH2dust cloud were studied using a modified Hartmann tube.The flame structures and propagation processes were recorded with a high-speed camera,and the transient two-dimensional temperature distribution of MgH2dust explosion could be reconstructed with the high-speed images and a Python code.Furthermore, the combustion residues of MgH2dust were analyzed by the scanning electron microscopy and x-ray diffraction method, and the combustion mechanism of MgH2dust was also discussed on basis of the deflagration characteristics and combustion residues analysis.

    2.Experiments

    2.1.Experimental materials

    In the experiments,MgH2(hydrogen storage capacity of 7.6 wt%,Alfa Aesar, USA) dust was sieved into four groups of samples with different particle size.The particle size distribution and microstructure of the four different samples were measured by a laser particle size analyzer (Mastersizer 2000, Malvern, UK) and a scanning electron microscope (SEM, VEGA3 SBH, TESCAN), as shown in Fig.1.The mean diameters (D50) of MgH2powders were 23 μm,40 μm,60 μm and 103 μm,respectively.Fig.1 shows that the particle surfaces of MgH2powders were smooth and spherical,and those withD50of 23 μm had a slight agglomeration phenomenon.

    2.2.Experimental apparatus

    A modified 1.2 L Hartmann tube was used to create an open space for studying the dust explosions, as shown in Fig.2.It consisted of a combustion tube, a gas supplying system, a dispersion system,an ignition system,a high-speed photography system and a time controller system,Before the ignition,the movable part in the center of Hartmann tube was fixed, and the inlet solenoid valve,ignition electrode, retractable electromagnet and camera system were synchronously controlled by a programmable logic controller.In the experiments, 0.5 MPa compressed air was used to disperse the MgH2particles.After 0.12 s delay of dust dispersion, the suspended particles were ignited with a 15 kV high voltage transformer.At the same time, a high-speed camera (Memrecam HX-3,NAC Japan)equipped with a macro lens(AF Micro-Nikkor 200 mm f/4D IF-ED, NAC Japan) was triggered to record the propagation process of the dust flame.At the moment of ignition, the movable tube was moved down to provide an open space for the combustible dust cloud, so that the flame could propagate in the open space without the restriction of the tube wall,facilitating the direct observation of flame structure and propagation behaviors.

    Fig.1.Particle size distribution and micro-structure of MgH2 powders: (a) 23 μm; (b) 40 μm; (c) 60 μm; (d) 103 μm.

    Fig.2.Schematic of the experimental apparatus.

    2.3.Calibration of temperature measurement system

    As a non-contact combustion diagnostic technique, two-color pyrometer technique has been applied in measuring particle temperature of metal dust [27,28].The two-color pyrometer system was calibrated before measuring the temperature field of MgH2.In the experiment, the tungsten filament lamp was selected for temperature calibration(tungsten is the metal with the highest melting point in the world,and its melting point is as high as 3695 K)[26].By calculating the corresponding relationship between theR/Gratios of the camera sensor and the temperature of the tungsten filament lamp, the variation trend of correction coefficient with temperature was acquired through the calibration (as shown in Fig.3) [29].

    Fig.3.Temperature calibration with tungsten filament lamp: (a) schematic diagram of experimental set-up; (b) temperature-(R/G) ratio fitting curve.

    3.Results and discussions

    3.1.Kinetics of MgH2 dehydrogenation and oxidation

    The thermal decomposition and dehydrogenation characteristics of MgH2samples were studied by the thermogravimetrydifferential scanning calorimetry under the argon atmosphere,and the thermogravimetry (TG) and differential scanning calorimetry(DSC)curves of MgH2with different particle sizes are shown in Fig.4.Since the thermal expansion coefficient of MgH2is much larger than that of MgO [17], MgH2expanded after being heated,which would break the MgO film layer outside, and the internal MgH2would be decomposed into H2and Mg.The TG curves and DSC curves had a descending peak and an endothermic peak at 455-475°C, respectively, indicating that the dehydrogenation reaction began in this temperature range.As the particle size decreased, the mass loss was greater, so the hydrogen content of the small sized MgH2particles was higher.By comparing the TG and DSC curves of MgH2with different particle sizes under the argon atmosphere,it could be found that the smaller particle size of MgH2had the lower initial decomposition temperature,maximum heat absorption and energy required for hydrogen release.

    Fig.5 shows the TG and DSC curves of MgH2powders in air atmosphere.MgH2powders with different particle sizes had two rapid weight gain stages in the range of 370-450°C and 600-650°C, respectively.The first rapid weight gain stage (dehydrogenation and oxidation stage of MgH2powders), in which the TG curve rose rapidly around 370°C with the increase of temperature, corresponding to the hydrogen release from the particle interior, the diffusion combustion of the particle surface, and the oxidation reaction of the decomposition product Mg particles,and comparing Figs.4(a)and 5(a),it could be seen that the weight gain of Mg particles oxidation was much greater than the weight loss of hydrogen release, that was why the weight of sample increased.The second rapid weight gain stage (high-temperature oxidation stage of Mg particles),as in Fig.4(a),there was no dehydrogenation reaction of MgH2powders at 600-650°C,so the weight gain of the sample was mainly caused by the complete oxidation of the remaining Mg particles(boiling combustion),as in Fig.5(a),the TG curve started to show a second rise at 600-650°C.In addition, as shown in Fig.5(b), the first exothermic peak of the DSC curve of MgH2powders in air atmosphere appeared at 370-450°C,but the endothermic peak of the DSC of MgH2powders in argon atmosphere did not appear at this stage in Fig.4(b),which was because the endothermic reaction of MgH2powders dehydrogenation overlapped with the oxidation reaction,and the heat release of the oxidation process was much greater than the heat absorption of the decomposition hydrogenation process, so the overall performance presented an exothermic reaction.The DSC curve of MgH2powders in air atmosphere had a second exothermic peak at 600-650°C,which was caused by the complete oxidation(boiling combustion)of the remaining Mg particles.

    Fig.4.Thermal analysis for pure MgH2 at the heating rate of 10°C/min in argon: (a) TG curves; (b) DSC curves.

    Fig.5.TG and DSC results for pure MgH2 at the heating rate of 10°C/min in air: (a) TG; (b) DSC.

    Fig.5 shows that in the first weight gain stage (370-450°C),compared with the larger particles with the particle sizes of 60 μm and 103 μm, the particles of 23 μm and 40 μm had lower initial oxidation temperatures,higher weight gain and lower heat release.The reasons for that may concluded as follows:smaller particle size resulted in lower apparent activation energy of MgH2,leading to a lower initial oxidation temperature, and smaller particles had larger contact area with air to be oxidized more completely and weight gained more obviously, furthermore, due to the higher hydrogen content of the smaller MgH2powders, as the dehydrogenation reaction needed to absorb heat, so the heat released by the first rapid weight gain stage of smaller MgH2powders were lower.In the second rapid weight gain stage (600-650°C), during which a high-temperature oxidation(boiling combustion)reaction occurred, and the Mg powders were completely reacted, so the smaller MgH2powders had less weight gain and heat release due to the relative lower Mg content.

    3.2.Flame propagation behaviors

    The flame propagation processes in the MgH2dust clouds(830 g/m3concentration) with four different particle size distributions are shown in Fig.6.The initial shape of the flames was approximately spherical,and the flame fronts were smooth.As the combustion progressed, the flame fronts gradually became irregular, and bright cluster flames appeared in the lower part of the combustion zones, which was attributed to the uneven dispersion of dust particles caused by the gravity effect and turbulence[28].It was notable that dust cloud with the smaller particle size had a brighter flame and a faster propagation velocity.According to the diffusion oxidation mechanism of Mg particles, particle size is the main factor affecting the oxidation rate of Mg particles [27,30].Smaller particle had a longer suspension time and a larger specific surface area, which led to the faster heat rate and more violent combustion reaction [31].Furthermore, for both hydrogen release rate and hydrogen storage capacity of MgH2powders with small particle size were higher than those with large particle size,and the combustible gas H2released around the particles would reduce the influence of the 'gap effect' between particles [32,33], thus accelerated the combustion reaction of particles.

    Fig.7 is the flame propagation process of 23 μm Mg dust cloud with concentration of 830 g/m3.Mg and MgH2dust clouds had different flame phenomena.Comparing Figs.6(a) and 7, it can be seen that the flame propagation velocity of Mg dust cloud was significantly lower than that of MgH2dust cloud,that was because MgH2dust cloud would release H2during its burning process,forming the gas-solid two-phase combustion of H2and dust mixture.The heat released by H2combustion would accelerate the combustion of MgH2powders(a'positive feedback'effect[32]),and the combustion product water vapor would continue to react with Mg to produce new H2, this reciprocating cycle of feedback effect would accelerate the total combustion reaction [34].In addition,unlike the combustion of Mg dust cloud,a large number of floating flocculus could be observed in the MgH2dust cloud flame.In the combustion process of MgH2dust cloud, the combustion product MgO powders would diffuse outward under the action of H2and water vapor, and they gradually cool ed and eventually formed the flocculent structure with water in the propagation process.

    Fig.8(a)and 8(b)show the flame front propagation distance and velocity as a function of time in the MgH2dust clouds with four different particle sizes, respectively.In order to reduce the influence of gravity and buoyancy, the horizontal displacement of the flame was used to calculate the flame propagation velocity.As shown in Fig.8,in the 23 μm,40 μm,60 μm MgH2dust clouds with the concentration of 830 g/m3,there was a significant acceleration of flame propagation in the initial stage, while that in the 103 μm MgH2dust cloud with the same concentration was relatively stable.The average flame propagation velocities of 23 μm, 40 μm, 60 μm and 103 μm MgH2dust clouds in the stable propagation stage were 3.7 m/s, 2.8 m/s, 2.1 m/s and 0.9 m/s, respectively.The dust cloud with a smaller particle size had the faster flame propagation velocity.This was because smaller particle had larger specific surface area, shorter preheating time, and faster particle decomposition/gasification rate, so it was much easier to release H2.The combustion of H2had a large improvement in the flame velocity and initial acceleration[35].The flame propagation velocities were not constant but fluctuated.The flame velocity fluctuation phenomenon could be attributed to the velocity slip between gases and particles.During the flame propagation, MgH2powders would release H2and Mg vapors when heated.Initially,the velocity of gas was higher than that of the particle because of strong expansion;thus the particle number density ahead of the flame front would decrease due to the higher inertia of particles.The reduced particle number density would decrease the rate of surface reactions and consequently the heat release by both gas-phase and particlephase reactions.This decrease in reaction intensity would lower both gas and the particle velocities.However, gas responses were much faster than those by the particles.The consequence was that the particles were moving much faster than gas,and this essentially increased the particle number density ahead of the flame front.As a result,the reactions become intensified because of a higher particle number density.This led to more heat release and enhanced the gas and particle velocities, but the gas velocity increased more than particle velocity.This alternating change of gas and particle velocities resulted in a periodic change of local particle number density (that was, the local fuel equivalence ratio), and consequently it caused flame-speed oscillation [36].Furthermore, the oscillation frequency of the flame was irregular due to the effect of turbulence and buoyancy [37], and the oscillation phenomenon became stronger with the decreasing particle size.The thermal decomposition rate of smaller particles was relatively faster, and the change of its gas phase velocity was greater than that of large particles,so the oscillation amplitude enhanced for the velocity slip between the smaller particles and gas was larger.

    Fig.6.Flame propagation behaviors of MgH2 dust clouds with different particle size: (a) 23 μm MgH2 dust cloud; (b) 40 μm MgH2 dust cloud; (c) 60 μm MgH2 dust cloud; (d)103 μm MgH2 dust cloud.

    Fig.7.Flame propagation of 23 μm Mg dust cloud with concentration of 830 g/m3.

    3.3.Flame micro-structures

    The local magnification view of the combustion processes of 103 μm MgH2dust cloud are shown in Fig.9.After the ignition of MgH2powders, the smaller particles in the flame front were decomposed or gasified to form a local premixing flame.The surrounding large particles absorbed heat through thermal radiation and thermal convection, then participated in the combustion reaction to gradually form a propagating cluster flame.The light blue flame appeared during the combustion of the MgH2dust cloud,which was caused by the combustion of H2released from MgH2powders.As the combustion intensified,the decomposition rate of MgH2powders accelerated to release more H2, and the light blue flames became more visible.Different from the Mg dust cloud,there were many micro-diffuse flames around the burning particles in the combustion processes of MgH2dust clouds.At the initial stage of particle combustion, the particles had a bright spherical structure and continued to form a micro-diffusion flame(the effect of H2combustion), which rapidly expanded and ignited the unburned dust particles in the preheating zone [38].

    As illustrated in Fig.10,the phenomenons of particle explosions and ejections of smaller luminescent substances often occurred during the combustion of MgH2powders,which was named microexplosion [28].The direct cause of the micro-explosion was the rapid increase of internal pressure due to the gas expansion inside of the particles, and the particles would be broke into small fragments as the internal stress exceeded the maximum bearing capacity of the liquid particle surface.The probability and intensity of micro-explosion depended largely on the speeds of generation and expansion of gas bubbles within the particles [39].There were three potential sources of gas inside MgH2particles: first, a large amount of H2was stored inside of the particles;second,due to the high temperature inside of the particles, the vapor of Mg or its oxides might be produced;third,the surrounding gas would diffuse and dissolve into the molten MgH2powders.The occurrence of micro-explosions inhibited particle agglomeration to a certain extent, which was conducive to improving the rate of hydrogen release and the combustion of MgH2powders [40,41].

    Fig.9.Flame micro-structure of 103 μm MgH2 dust cloud.

    Fig.10.Micro-explosion phenomenon occurred in MgH2 dust flame propagation.

    The micro-explosion process of a single MgH2particle is shown in Fig.11.The MgH2particle surface burned and formed a bright tail-like light.As the temperature increased,the thermal expansion of the liquid or vapor Mg inside the particles increased the internal pressure of the particles.When the stress reached a certain threshold, the particles began to expand and rupture, ejecting Mg vapor to the surroundings and eventually leading to microexplosion phenomenon [42].

    3.4.Flame temperature distribution

    In the experiments, the combustion processes of MgH2dust clouds were recorded by a high speed camera, and the corresponding temperature distributions were calculated and rebuilt by the two-color pyrometer technique, as shown in Fig.12.The temperature of each propagation direction in the early stage of the flame was almost the same,and the temperature of the internal and external flames were also very close.With the development of the flame propagation,the temperature of each stage increased,but the temperature of the flame front was relatively lower than that of inside, which was closely related to the dehydrogenation reaction of MgH2powders and the gasification of Mg.Fig.4(a)demonstrates that the dehydrogenation temperature of MgH2powder was low and its dehydrogenation rate was fast, and a large amount of H2would be quickly released after heating.H2diffused and burned rapidly at the flame front, and the decomposed product Mg particles inside the flame vaporized and underwent a combustion reaction at the same time.For the combustion temperature of Mg powders was higher than that of H2[27],and the decomposition of MgH2powders at the flame front would absorb part of the heat generated by H2combustion, so that the temperature at the dust flame front was lower than that inside.Furthermore,by comparing the temperature distributions of MgH2dust clouds with different particle sizes, it could be seen that local high temperature areas would appear at the bottom part of the dust cloud,and the overall flame temperature distributions in the later stage of the dust clouds with small particle sizes were more uniform and the temperature gradient were smaller than those in the dust clouds with large particle size.That was because the gravity effect could not be neglected with the increasing particle size [28].In this regard, we performed analysis from the gravitational sedimentation perspective on the spatial distribution characteristics of dust flame temperature.

    Firstly, assuming that these particles were spheres with a diameter ofdp, they were sedimenting under the forces including gravity, buoyancy, and airflow drag force [43].According to the Newton's second law,

    wherempis the particle mass;upis the particle velocity;ugis the flow velocity of air,ρpis the density of particles;ρgis the density of air;dpis the diameter of particles;ζdis the drag coefficient;gis the gravitational acceleration.

    The drag force increases as the air velocity increases.When a certain velocityutis reached, the sum of gravity, buoyancy and airflow drag force acting on the particles is zero, and the particles fall at a uniform velocity ofut,which is the final settling velocity of the particles.Based on Eq.(1), when dup/dt= 0, the final falling velocity of the particles can be calculated.

    Fig.11.Micro-explosion of a single MgH2 particle.

    Fig.12.Flame temperature distribution in micron MgH2 dust clouds:(a)23 μm MgH2 dust cloud;(b) 40 μm MgH2 dust cloud;(c)60 μm MgH2 dust cloud;(d)103 μm MgH2 dust cloud.

    Fig.12.(continued).

    The drag coefficient ζdchanges with Reynolds numberRe, and, when Reynolds numberRe≤1, the drag coefficient ζd=24/Re,so the particle sedimentation formula can be expressed as

    where ν is the kinematic viscosity of air, ν = 15.7 × 10-6m2/s;ρg= 1.20 kg/m3.

    The falling velocity of the particles were approximated calculated according to the density of MgH2(1450 kg/m3),and the results are listed in Table 1.Eq.(3)shows that the falling velocity of the particles was proportional to the square of the particle size.Smaller particle had lower falling velocity of particles in the air,which was more conducive to the suspension and diffusion movement of particles as well as the heat transfer between particles.As shown in Table 1,the falling velocities of 23 μm and 40 μm MgH2particles were 2.24 cm/s and 6.71 cm/s, respectively, which was much smaller than those of 60 μm and 103 μm MgH2particles.Therefore, the flame shape were more spherical and the temperature distribution were more uniform in the 23 μm and 40 μm MgH2dust clouds.While the falling velocities of 60 μm and 103 μm MgH2particles were as high as 15.07 cm/s and 44.42 cm/s, respectively,leading to a more rapid downward development and irregular shape of the flame, which also explained why the cluster flame usually appeared at the bottom of the MgH2dust cloud with large particle sizes.

    Table 1Falling velocity of MgH2 particles.

    3.5.Analysis of MgH2 combustion residues

    Chemical composition analysis of solid combustion residues could invert the development process of MgH2dust explosion[44].Fig.13 shows the micro-structures of solid residues with a scanning electron microscopy(SEM).As a reference,Fig.13(a)demonstrates the SEM picture of an unreacted MgH2particle.Fig.13(b) shows that the solid combustion residues of MgH2particle became loose,porous and larger than the unburned ones.Floating flocculus formed during the MgH2combustion covered the surface of the solid combustion product,as shown in Fig.13(c).With the increase of temperature and the deepening of dehydrogenation degree,the pores on the surface of the combustion residues became larger and penetrated deep into the particles or even produced cracks under the action of micro-explosion, as shown in Fig.13(d).To further understand the combustion reaction of the MgH2dust cloud,X-ray diffraction experiments were performed on its combustion residues, and the results are shown in Fig.14.Fig.14 shows that only MgO and Mg were contained in the solid residues after MgH2dust cloud combustion,indicating that MgH2particles were completely decomposed in the combustion reaction, and MgO was the oxidation product of the combustion reaction of MgH2particles, while Mg was the decomposition product of MgH2that did not participate in the combustion.

    Fig.13.Micro-structures of MgH2 combustion residues:(a)Pure MgH2 particle surface;(b)Holes on residue particle surface;(c)MgO floccules on residue particle surface;(d)Crack on residue particle surface.

    Fig.14.XRD patterns of MgH2 combustion residues.

    Combining the experimental results of the flame propagation behaviors, flame micro-structure characteristics and combustion residues composition analysis, the dehydrogenation and combustion processes of MgH2powders in air can be obtained,as shown in Fig.15.After the MgH2particles were heated, heterogeneous surface reaction occurred between particle surface and oxygen to form a protective MgO layer.Since H2would be released from MgH2particle, and the thermal expansion coefficient of MgH2particle was significantly larger than that of MgO [17], so holes and cracks would be formed on the surface of the MgH2particles.Meanwhile,the transformation of MgH2to Mg was accompanied by volume shrinkage, which led to further expansion of pores and cracks formed on the surfaces of MgH2particles [45].H2has the characteristics of low ignition point and wide explosion limit,and formed a gas-solid two-phase combustion reaction with MgH2particles,which enhanced the explosion intensity.In addition, the water vapor of H2combustion product could promote the hydrolysis reaction of Mg to released H2and large amounts of heat, thereby accelerating the propagation of MgH2dust flame.At the same time,oxygen could enter the interior of MgH2particles through the cracks, and the MgO layer gradually expanded to the inside of the particles, while some Mg vapor and H2could escape through the cracks of the MgO layer, so that the degree of MgH2dehydrogenation deepened gradually from the surface to the interior.MgH2particles were not completely melted in the initial heating stage,at this moment the interaction of MgH2particle with oxygen mainly occurred on its surface.After the ignition of MgH2dust cloud, the surface of the dust particle was completely melted and its vapor was generated, the released H2and metal vapor were mixed with oxygen and then vapor phase diffusion combustion occurred.The flocculent MgO generated by vapor phase combustion of particles diffused into the surrounding space, and part of them diffused inward and deposited on the surface of the solid residues.Therefore,the combustion control mechanism of MgH2particles was mainly heterogeneous surface reaction.Furthermore, the specific surface area,hydrogen release rate and hydrogen storage capacity of MgH2powders with small particle size were higher than those with large particle size, thus the dust clouds with smaller particle sizes had faster flame propagation velocity and more violent combustion reaction, and their flame temperature distributions were more uniform and the temperature gradients were smaller.

    Fig.15.Schematic of the dehydrogenation and passivation process of MgH2.

    4.Conclusions

    The experimental results of thermal analysis experiments showed that, under argon atmosphere, the dehydrogenation of MgH2particles occurred at 455-475°C.MgH2particles with smaller particle size had higher H2content, but the initial decomposition temperature, maximum heat absorption and the energy required for H2release were lower.Under air atmosphere, MgH2particles had two rapid weight gain stages in the range of 370-450°C and 600-650°C, respectively.In the first rapid weight gain stage,the initial oxidation temperature and heat release of the smaller particles were both lower but the weight gain was greater.In the second rapid weight gain stage,MgH2powders with smaller particle sizes had less weight gain and heat release due to less Mg content.

    Flame propagated in the MgH2dust clouds all presented an approximately spherical shape at the early stage, the average propagation velocity,initial acceleration,and oscillation amplitude of the dust flame were related to the particle size distributions of MgH2powders.The smaller particle sizes resulted in the greater flame propagation speed and initial acceleration.In addition, the velocity oscillation became more obvious as the particle sizes decreased,owing to that the thermal decomposition rate of smaller particles was relatively faster, and their gas phase velocities changed more significantly than larger particles, which made the velocity slip between particles and gas larger, and the oscillation amplitude increased.Furthermore, different from Mg dust clouds,MgH2dust clouds had many combustion particles with microdiffusion flame structures during the combustion processes.Micro-explosion occurred during the combustion of MgH2dust clouds to inhibit particle agglomeration in a certain extent, which was conducive to improving the combustion of MgH2powders.

    The flame temperature of MgH2dust cloud in the early stage was almost the same in all directions.With the development of the dust flame, the temperature of each stage increased, but the temperature of the flame front was lower than that inside.Local high temperature appeared at the bottom of the MgH2dust cloud due to the dust falling,and in the later stage of the dust clouds with small particle size, the overall flame temperature distribution was more uniform and the temperature gradient was smaller than those with large particle size.In addition, combining with the experimental results of flame propagation behaviors, flame micro-structure characteristics and combustion residues composition analysis, it could be concluded that the combustion control mechanism of MgH2particles was mainly heterogeneous surface reaction.

    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

    The work was supported by the National Natural Science Foundation of China (Grant Nos.12272001, 11972046), the Outstanding Youth Project of Natural Science Foundation of Anhui Province (Grant No.2108085Y02), and the Major Project of Anhui University Natural Science Foundation (Grant No.KJ2020ZD30),and the authors would like to thank these foundations for the financial supports.

    2021天堂中文幕一二区在线观 | 亚洲专区中文字幕在线| 在线观看www视频免费| 国产爱豆传媒在线观看 | 中文字幕另类日韩欧美亚洲嫩草| 十八禁人妻一区二区| 一区二区三区精品91| 99热6这里只有精品| 欧美一级毛片孕妇| 欧美不卡视频在线免费观看 | 男女之事视频高清在线观看| 俺也久久电影网| 一卡2卡三卡四卡精品乱码亚洲| 黄色丝袜av网址大全| 岛国在线观看网站| 国内揄拍国产精品人妻在线 | 看黄色毛片网站| 亚洲成人久久性| 国产精品,欧美在线| 男女视频在线观看网站免费 | 色老头精品视频在线观看| 国产野战对白在线观看| 亚洲精品中文字幕一二三四区| 午夜精品在线福利| 国产激情欧美一区二区| 欧美亚洲日本最大视频资源| 88av欧美| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久久久久免费视频| 每晚都被弄得嗷嗷叫到高潮| 欧美性猛交黑人性爽| 麻豆成人午夜福利视频| 两个人看的免费小视频| 最新美女视频免费是黄的| 日日干狠狠操夜夜爽| 大型黄色视频在线免费观看| 欧美另类亚洲清纯唯美| 欧美色视频一区免费| 国产高清有码在线观看视频 | 99国产精品一区二区蜜桃av| 久久久久国产精品人妻aⅴ院| 中文字幕人妻熟女乱码| 日韩精品免费视频一区二区三区| 久久青草综合色| 国产精华一区二区三区| 亚洲av美国av| 天堂动漫精品| 久久久水蜜桃国产精品网| av视频在线观看入口| av片东京热男人的天堂| 欧美日本视频| 巨乳人妻的诱惑在线观看| 国产伦人伦偷精品视频| 国产精品 欧美亚洲| 亚洲国产精品合色在线| 色综合亚洲欧美另类图片| av超薄肉色丝袜交足视频| 亚洲avbb在线观看| 国产黄片美女视频| 久久中文看片网| 大香蕉久久成人网| 久久国产亚洲av麻豆专区| 久久久久久九九精品二区国产 | 久99久视频精品免费| 中文字幕久久专区| 十八禁人妻一区二区| 成人18禁高潮啪啪吃奶动态图| 国产精品免费视频内射| 久久精品国产清高在天天线| 午夜福利在线观看吧| 欧美国产日韩亚洲一区| 国产精品亚洲一级av第二区| 国产精品二区激情视频| a级毛片在线看网站| 在线观看免费午夜福利视频| 国产精华一区二区三区| 国产一卡二卡三卡精品| 欧美三级亚洲精品| 观看免费一级毛片| www.www免费av| 久久精品影院6| 999久久久精品免费观看国产| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品中文字幕一二三四区| 巨乳人妻的诱惑在线观看| 久久99热这里只有精品18| 深夜精品福利| 久久久精品欧美日韩精品| 18禁美女被吸乳视频| 亚洲精华国产精华精| 制服人妻中文乱码| 日本免费一区二区三区高清不卡| 精品国产乱码久久久久久男人| 亚洲成国产人片在线观看| 日韩国内少妇激情av| 亚洲欧美日韩无卡精品| 最近最新免费中文字幕在线| 日日爽夜夜爽网站| 美女大奶头视频| 最近最新中文字幕大全免费视频| 99热这里只有精品一区 | 成熟少妇高潮喷水视频| 国产精品亚洲一级av第二区| 欧美日韩中文字幕国产精品一区二区三区| 欧美色欧美亚洲另类二区| 久久久久免费精品人妻一区二区 | 别揉我奶头~嗯~啊~动态视频| 国产在线精品亚洲第一网站| 黄片播放在线免费| 日韩欧美一区视频在线观看| 午夜福利一区二区在线看| 欧美黄色片欧美黄色片| 久久中文看片网| 成人亚洲精品一区在线观看| 亚洲七黄色美女视频| a在线观看视频网站| 中文字幕精品亚洲无线码一区 | 99久久99久久久精品蜜桃| 亚洲最大成人中文| 亚洲熟女毛片儿| 91成人精品电影| 国产蜜桃级精品一区二区三区| 久久久久久大精品| 国产精品久久久久久人妻精品电影| 国产真实乱freesex| 欧美另类亚洲清纯唯美| 熟妇人妻久久中文字幕3abv| 国产男靠女视频免费网站| 日本a在线网址| 亚洲国产精品久久男人天堂| 一区二区三区高清视频在线| 听说在线观看完整版免费高清| 性欧美人与动物交配| 欧美精品亚洲一区二区| 十八禁网站免费在线| 久久这里只有精品19| 亚洲一区二区三区不卡视频| 国产欧美日韩一区二区三| 亚洲成av人片免费观看| 国产伦一二天堂av在线观看| 欧美激情 高清一区二区三区| 久久久国产成人免费| 午夜日韩欧美国产| 日韩精品免费视频一区二区三区| 精品熟女少妇八av免费久了| 久久国产亚洲av麻豆专区| 久久久国产欧美日韩av| a在线观看视频网站| 丁香六月欧美| 无遮挡黄片免费观看| 男女视频在线观看网站免费 | 免费在线观看影片大全网站| 国产高清有码在线观看视频 | netflix在线观看网站| 久久久久久人人人人人| 99精品在免费线老司机午夜| 国产日本99.免费观看| 黑人欧美特级aaaaaa片| 老司机福利观看| 精品国产乱子伦一区二区三区| 亚洲全国av大片| 国产伦人伦偷精品视频| 日日夜夜操网爽| 国产精品影院久久| 国产精华一区二区三区| 99国产精品一区二区三区| 丰满的人妻完整版| 成年免费大片在线观看| 久久精品91蜜桃| 91成年电影在线观看| 亚洲成a人片在线一区二区| 啦啦啦观看免费观看视频高清| 久久人妻av系列| 国产免费男女视频| 亚洲中文字幕一区二区三区有码在线看 | 久久久国产成人精品二区| 老熟妇乱子伦视频在线观看| 90打野战视频偷拍视频| 日韩欧美国产在线观看| 午夜福利18| 免费看十八禁软件| 亚洲专区字幕在线| 在线观看午夜福利视频| 国产久久久一区二区三区| 成人国产综合亚洲| 国产精品亚洲一级av第二区| 亚洲 欧美 日韩 在线 免费| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品美女久久av网站| 一级片免费观看大全| 最新在线观看一区二区三区| 国产成人av激情在线播放| 亚洲av五月六月丁香网| 亚洲全国av大片| 亚洲国产欧洲综合997久久, | 婷婷亚洲欧美| 国产av在哪里看| 国产黄a三级三级三级人| 大香蕉久久成人网| 两个人看的免费小视频| 久久久久久免费高清国产稀缺| 美女免费视频网站| 宅男免费午夜| 99国产综合亚洲精品| 女生性感内裤真人,穿戴方法视频| 精品国产乱码久久久久久男人| 18禁裸乳无遮挡免费网站照片 | 亚洲精品国产一区二区精华液| 美女扒开内裤让男人捅视频| 中文在线观看免费www的网站 | 精品欧美一区二区三区在线| 欧美黄色淫秽网站| 国产97色在线日韩免费| 亚洲成国产人片在线观看| 亚洲人成伊人成综合网2020| 日本一本二区三区精品| 国产免费男女视频| www日本黄色视频网| 亚洲av电影在线进入| 国产成人av教育| 欧美精品啪啪一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 一边摸一边抽搐一进一小说| 天天一区二区日本电影三级| 久久久久久大精品| 免费在线观看日本一区| 国产私拍福利视频在线观看| 色av中文字幕| 午夜影院日韩av| 国产色视频综合| 熟妇人妻久久中文字幕3abv| 午夜福利免费观看在线| 国产在线精品亚洲第一网站| 午夜福利一区二区在线看| 国产日本99.免费观看| 中文亚洲av片在线观看爽| 听说在线观看完整版免费高清| 久久国产精品人妻蜜桃| 一区福利在线观看| 首页视频小说图片口味搜索| 18禁观看日本| 好男人在线观看高清免费视频 | 久久九九热精品免费| 一本大道久久a久久精品| 天堂动漫精品| 国产又黄又爽又无遮挡在线| 国产不卡一卡二| 无人区码免费观看不卡| 欧美日韩黄片免| 免费人成视频x8x8入口观看| 欧美乱码精品一区二区三区| 久久久水蜜桃国产精品网| www.精华液| 韩国精品一区二区三区| 国产真实乱freesex| 一级a爱视频在线免费观看| 黑丝袜美女国产一区| 久久热在线av| 午夜免费观看网址| 亚洲精品国产一区二区精华液| 此物有八面人人有两片| 欧美日韩福利视频一区二区| 少妇 在线观看| 哪里可以看免费的av片| 亚洲人成电影免费在线| 脱女人内裤的视频| 亚洲欧美精品综合一区二区三区| 欧美日韩精品网址| 国产蜜桃级精品一区二区三区| 宅男免费午夜| 精品乱码久久久久久99久播| 亚洲精品色激情综合| 哪里可以看免费的av片| 听说在线观看完整版免费高清| 国产亚洲精品久久久久5区| 精品久久久久久久毛片微露脸| 国产成人精品久久二区二区91| 免费人成视频x8x8入口观看| 后天国语完整版免费观看| 伦理电影免费视频| 身体一侧抽搐| 国产一区二区三区在线臀色熟女| 啦啦啦韩国在线观看视频| 国产主播在线观看一区二区| 成熟少妇高潮喷水视频| 美女大奶头视频| 非洲黑人性xxxx精品又粗又长| 夜夜看夜夜爽夜夜摸| 亚洲aⅴ乱码一区二区在线播放 | 国产在线观看jvid| 久久久国产成人精品二区| 男人舔女人的私密视频| 欧美日韩一级在线毛片| 一区二区三区激情视频| 一边摸一边抽搐一进一小说| 变态另类成人亚洲欧美熟女| 亚洲五月天丁香| 日韩国内少妇激情av| 一个人免费在线观看的高清视频| 真人做人爱边吃奶动态| 亚洲熟妇熟女久久| 啦啦啦韩国在线观看视频| 国产91精品成人一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 最好的美女福利视频网| 欧美色欧美亚洲另类二区| 老熟妇仑乱视频hdxx| 在线观看免费日韩欧美大片| 黄网站色视频无遮挡免费观看| 99久久无色码亚洲精品果冻| 亚洲国产精品合色在线| 成人午夜高清在线视频 | 欧美成人一区二区免费高清观看 | 欧美乱码精品一区二区三区| 天天添夜夜摸| 国产熟女xx| 久久久久久大精品| 久久 成人 亚洲| 久久狼人影院| 18禁裸乳无遮挡免费网站照片 | 亚洲欧美精品综合久久99| 黄网站色视频无遮挡免费观看| 久久久久久九九精品二区国产 | 久久人人精品亚洲av| 精品久久久久久,| 日本 av在线| 久久香蕉国产精品| 欧美色视频一区免费| 成人亚洲精品av一区二区| 麻豆一二三区av精品| 在线观看免费视频日本深夜| 天天躁夜夜躁狠狠躁躁| 天天躁狠狠躁夜夜躁狠狠躁| 黑丝袜美女国产一区| 成年免费大片在线观看| 国产成人精品久久二区二区免费| 男女下面进入的视频免费午夜 | 男人舔女人的私密视频| 99久久国产精品久久久| 欧美又色又爽又黄视频| 在线看三级毛片| www.自偷自拍.com| 久久久精品国产亚洲av高清涩受| 亚洲 国产 在线| 亚洲色图 男人天堂 中文字幕| 国产亚洲精品久久久久久毛片| 久久精品91无色码中文字幕| 在线观看66精品国产| a级毛片在线看网站| 免费无遮挡裸体视频| а√天堂www在线а√下载| 搡老熟女国产l中国老女人| 天天添夜夜摸| 中文字幕精品免费在线观看视频| 久久精品国产亚洲av高清一级| 美国免费a级毛片| 欧美性猛交╳xxx乱大交人| 亚洲国产欧美网| av有码第一页| 成人精品一区二区免费| 听说在线观看完整版免费高清| 色播在线永久视频| 俄罗斯特黄特色一大片| 日本熟妇午夜| 国产精品久久久久久人妻精品电影| 岛国在线观看网站| 国产精品久久久久久人妻精品电影| 色老头精品视频在线观看| 欧美成人性av电影在线观看| 午夜激情av网站| 欧美黑人巨大hd| 777久久人妻少妇嫩草av网站| 大型av网站在线播放| 亚洲自偷自拍图片 自拍| 免费看美女性在线毛片视频| 日韩大尺度精品在线看网址| 黄色成人免费大全| 男女床上黄色一级片免费看| 青草久久国产| 欧美日韩一级在线毛片| 国产99久久九九免费精品| 亚洲欧美日韩无卡精品| 中亚洲国语对白在线视频| 亚洲成a人片在线一区二区| 欧美大码av| 亚洲专区国产一区二区| 国产高清激情床上av| 成人午夜高清在线视频 | 身体一侧抽搐| 色综合亚洲欧美另类图片| 中亚洲国语对白在线视频| √禁漫天堂资源中文www| 曰老女人黄片| 久热这里只有精品99| 亚洲成av人片免费观看| 国产精品国产高清国产av| 两个人视频免费观看高清| 在线十欧美十亚洲十日本专区| 国产亚洲精品久久久久5区| 日本熟妇午夜| 中亚洲国语对白在线视频| 在线观看免费视频日本深夜| 1024手机看黄色片| 男人舔女人的私密视频| 欧美黄色淫秽网站| 国产主播在线观看一区二区| 两个人看的免费小视频| av有码第一页| 在线观看舔阴道视频| 日本免费a在线| 在线十欧美十亚洲十日本专区| 免费人成视频x8x8入口观看| 亚洲熟女毛片儿| 一本一本综合久久| 日韩欧美国产一区二区入口| 成年女人毛片免费观看观看9| 99re在线观看精品视频| 午夜久久久久精精品| 日本精品一区二区三区蜜桃| 亚洲av中文字字幕乱码综合 | 美女国产高潮福利片在线看| 又大又爽又粗| 久久久久久九九精品二区国产 | 亚洲一区二区三区不卡视频| 大型av网站在线播放| 亚洲专区国产一区二区| 精品久久久久久久毛片微露脸| www.www免费av| 最新美女视频免费是黄的| 夜夜看夜夜爽夜夜摸| 91成年电影在线观看| 欧美丝袜亚洲另类 | 国产亚洲av高清不卡| 亚洲男人天堂网一区| 校园春色视频在线观看| 精品电影一区二区在线| 久久久久久亚洲精品国产蜜桃av| 国产在线观看jvid| 日韩欧美国产在线观看| 国产精品九九99| 好看av亚洲va欧美ⅴa在| 成人午夜高清在线视频 | 人人妻人人澡人人看| 亚洲美女黄片视频| 国产精品永久免费网站| 长腿黑丝高跟| 精品国产国语对白av| 国产高清有码在线观看视频 | 亚洲全国av大片| 最新在线观看一区二区三区| 国产av一区二区精品久久| 黄色a级毛片大全视频| 女生性感内裤真人,穿戴方法视频| 中文资源天堂在线| 亚洲黑人精品在线| 欧美三级亚洲精品| 亚洲欧美激情综合另类| 中文字幕人成人乱码亚洲影| 91大片在线观看| 欧美激情极品国产一区二区三区| 国产亚洲精品一区二区www| 欧美成人性av电影在线观看| 一区二区三区激情视频| 很黄的视频免费| 国产不卡一卡二| 九色国产91popny在线| 国产精品九九99| 又黄又粗又硬又大视频| 国产成人影院久久av| 亚洲成人久久性| 桃色一区二区三区在线观看| av欧美777| 国产精品一区二区免费欧美| 精品久久久久久久毛片微露脸| 757午夜福利合集在线观看| 日韩大码丰满熟妇| 精品电影一区二区在线| tocl精华| 日本五十路高清| 欧美日韩亚洲国产一区二区在线观看| 日本精品一区二区三区蜜桃| 亚洲一区二区三区不卡视频| 久久精品成人免费网站| 999久久久国产精品视频| 精品久久蜜臀av无| 精品乱码久久久久久99久播| 国产亚洲精品综合一区在线观看 | 国产亚洲欧美98| 国产男靠女视频免费网站| 欧美日韩福利视频一区二区| 国产亚洲av高清不卡| 69av精品久久久久久| 日日夜夜操网爽| 午夜福利18| 国产真实乱freesex| 少妇熟女aⅴ在线视频| 十八禁网站免费在线| 正在播放国产对白刺激| 国产一区在线观看成人免费| 日韩一卡2卡3卡4卡2021年| 成人国产一区最新在线观看| 18禁黄网站禁片午夜丰满| 国内少妇人妻偷人精品xxx网站 | 精品人妻1区二区| 亚洲人成电影免费在线| 久久久久久人人人人人| 午夜精品在线福利| 亚洲电影在线观看av| 在线观看免费日韩欧美大片| 亚洲三区欧美一区| 视频在线观看一区二区三区| 侵犯人妻中文字幕一二三四区| 中亚洲国语对白在线视频| 变态另类成人亚洲欧美熟女| 免费看十八禁软件| 亚洲五月色婷婷综合| 久久久久国内视频| 国产亚洲精品久久久久久毛片| 亚洲av第一区精品v没综合| 窝窝影院91人妻| 久久婷婷人人爽人人干人人爱| 老司机午夜十八禁免费视频| 欧美三级亚洲精品| 亚洲av熟女| 欧美日本视频| 日韩三级视频一区二区三区| 天堂影院成人在线观看| 久久人人精品亚洲av| 久9热在线精品视频| 一二三四在线观看免费中文在| 国产亚洲精品久久久久5区| 精品久久久久久久人妻蜜臀av| 亚洲,欧美精品.| 久久伊人香网站| 国产精品,欧美在线| 日韩欧美国产在线观看| 日韩精品免费视频一区二区三区| 国产激情久久老熟女| 美女午夜性视频免费| 首页视频小说图片口味搜索| 欧美av亚洲av综合av国产av| 性欧美人与动物交配| 久久国产乱子伦精品免费另类| 97超级碰碰碰精品色视频在线观看| 国产男靠女视频免费网站| 欧美丝袜亚洲另类 | 国产精品久久久久久人妻精品电影| 老熟妇乱子伦视频在线观看| 母亲3免费完整高清在线观看| 国产精品1区2区在线观看.| 男女视频在线观看网站免费 | 少妇粗大呻吟视频| 99精品久久久久人妻精品| www.自偷自拍.com| 亚洲免费av在线视频| 丁香欧美五月| 午夜成年电影在线免费观看| 久久精品国产清高在天天线| 久久午夜综合久久蜜桃| 久久中文字幕一级| 不卡一级毛片| 亚洲精品国产一区二区精华液| 国产精品日韩av在线免费观看| 国产精品一区二区精品视频观看| 无人区码免费观看不卡| 久久精品国产清高在天天线| 国产伦一二天堂av在线观看| 好男人在线观看高清免费视频 | 看免费av毛片| 两个人视频免费观看高清| 久久久久久九九精品二区国产 | 亚洲人成网站在线播放欧美日韩| 国产亚洲欧美98| 国产精品免费一区二区三区在线| 可以在线观看毛片的网站| 亚洲性夜色夜夜综合| 色婷婷久久久亚洲欧美| 白带黄色成豆腐渣| 精品久久久久久久毛片微露脸| 久久国产亚洲av麻豆专区| 国产av不卡久久| 真人做人爱边吃奶动态| a级毛片a级免费在线| 亚洲午夜理论影院| 性色av乱码一区二区三区2| 白带黄色成豆腐渣| 俺也久久电影网| 国产成人欧美| 1024视频免费在线观看| 岛国在线观看网站| a在线观看视频网站| 国产一级毛片七仙女欲春2 | 亚洲欧美精品综合一区二区三区| 久久午夜亚洲精品久久| 精品电影一区二区在线| 国产v大片淫在线免费观看| 他把我摸到了高潮在线观看| 欧美最黄视频在线播放免费| 欧美日韩亚洲综合一区二区三区_| 老熟妇乱子伦视频在线观看| 国产色视频综合| 国产精品久久久久久亚洲av鲁大| 精品日产1卡2卡| 免费在线观看完整版高清| 国产1区2区3区精品| 免费看美女性在线毛片视频| 白带黄色成豆腐渣| 欧美日韩一级在线毛片| 久久精品国产清高在天天线| 国产精品乱码一区二三区的特点| 三级毛片av免费| 在线观看一区二区三区| 夜夜躁狠狠躁天天躁| 久久中文字幕一级| 最新在线观看一区二区三区|