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

    Plasma-assisted ammonia synthesis in a packed-bed dielectric barrier discharge reactor:roles of dielectric constant and thermal conductivity of packing materials

    2022-03-10 03:50:20JinLIU劉進(jìn)XinboZHU竺新波XueliHU胡學(xué)理andXinTU屠昕
    Plasma Science and Technology 2022年2期
    關(guān)鍵詞:學(xué)理

    Jin LIU(劉進(jìn)),Xinbo ZHU(竺新波),*,Xueli HU(胡學(xué)理) and Xin TU(屠昕)

    1 Faculty of Maritime and Transportation,Ningbo University,Ningbo 315211,People’s Republic of China

    2 Department of Electrical Engineering and Electronics,University of Liverpool,Liverpool L69 3GJ,United Kingdom

    Abstract In this article,plasma-assisted NH3 synthesis directly from N2 and H2 over packing materials with different dielectric constants(BaTiO3,TiO2 and SiO2)and thermal conductivities(BeO,AlN and Al2O3)at room temperature and atmospheric pressure is reported.The higher dielectric constant and thermal conductivity of packing material are found to be the key parameters in enhancing the NH3 synthesis performance.The NH3 concentration of 1344 ppm is achieved in the presence of BaTiO3,which is 106% higher than that of SiO2,at the specific input energy(SIE)of 5.4 kJ·l-1.The presence of materials with higher dielectric constant,i.e.BaTiO3 and TiO2 in this work,would contribute to the increase of electron energy and energy injected to plasma,which is conductive to the generation of chemically active species by electron-impact reactions.Therefore,the employment of packing materials with higher dielectric constant has proved to be beneficial for NH3 synthesis.Compared to that of Al2O3,the presence of BeO and AlN yields 31.0% and 16.9% improvement in NH3 concentration,respectively,at the SIE of 5.4 kJ·l-1.The results of IR imaging show that the addition of BeO decreases the surface temperature of the packed region by 20.5% to 70.3°C and results in an extension of entropy increment compared to that of Al2O3,at the SIE of 5.4 kJ·l-1.The results indicate that the presence of materials with higher thermal conductivity is beneficial for NH3 synthesis,which has been confirmed by the lower surface temperature and higher entropy increment of the packed region.In addition,when SIE is higher than the optimal value,further increasing SIE would lead to the decrease of energy efficiency,which would be related to the exacerbation in reverse reaction of NH3 formation reactions.

    Keywords:dielectric barrier discharge,ammonia synthesis,packing materials,thermal conductivity,dielectric constant

    1.Introduction

    Ammonia(NH3)is the world’s second largest chemical product and a vital precursor in the synthesis of many products.It is also an environmentally friendly fuel and excellent hydrogen storage carrier.Therefore,ammonia is an ideal alternative to fossil fuel[1-6].Conventional industrial-scale NH3synthesis via Haber-Bosch(HB)method consumes huge amount of energy(1%-2%of world’s primary energy supply)and requires large occupied area due to the severe working conditions,e.g.lasting high temperature(500°C-600°C)and high pressure(20-40 MPa),which limit the application HB method in medium and small scale NH3synthesis[7].In the past 100 years,researchers have devoted great efforts to develop more sustainable alternatives for HB method,including biochemical process,electro-chemical process and non-thermal plasma(NTP)assisted process[8,9].

    Among the NH3synthesis methods,NTP is particularly promising due to its characteristics of compact system,high energy efficiency(EE)and flexibility to combine with renewable energy sources.A variety of NTPs have been investigated for NH3synthesis,including dielectric barrier discharge(DBD)[10-14],radio frequency discharge[15,16]and glow discharge[17]plasmas.Recently,the DBD plasma has attracted a wide spread attention.The DBD reactor can be combined with packing materials,i.e.the packed-bed DBD reactor,to increase the reaction yield and EE.The combination of plasma and packing materials often leads to a synergistic effect,i.e.the enhancement in the reaction performance generated by the interactions between the plasma and the packing materials[18-20].To enhance the plasma-assisted NH3synthesis process,researchers have investigated the effect of different packing materials in the packed-bed DBD reactor.Akayet alemployed microporous silica-supported nickel catalysts packed with BaTiO3for NH3synthesis in a DBD reactor.When the ratio of catalyst increased from 0%to 100% in the catalyst-BaTiO3pellet mixture,the nitrogen conversion rate increased from 7.3% to 11.2%[21].Mizushimaet alevaluated the effects of Al2O3membranesupported transition metal catalysts(Ru,Ni,Pt,and Fe)on the synthesis of NH3using DBD plasma.Ru showed significant catalytic effects on the NH3formation in N2-H2plasma,although it had a lack contribution of the N2dissociation[22].Wanget aldemonstrated that the synergistic effect of plasma and a Ni/Al2O3catalyst appeared to be evident from 400°C,and the loaded Ni supported more chemically active sites[23].

    Although the synergistic effect of plasma-assisted NH3synthesis in a packed-bed DBD has long been investigated,there are only few studies on the influences of packing materials’ electrical properties on the discharge behavior and its contribution to NH3synthesis.In the presence of packing materials,the electrical characteristics of the DBD reactor would be affected together with the chemical processes in plasma region,e.g.the onset voltage,electron energy of plasma and uniformity of filamentary discharges,etc[7,8,24-29].The dielectric constant of packing materials played a dominant role in the plasma discharge process.Guet alreported that the amplitude of the current peak increased in presence of materials of higher dielectric constant in a coaxial DBD reactor.At the AC discharge frequency of 1 kHz and voltage amplitude of 2 kV,the maximum current amplitude increased 20 times when the relative dielectric constant of packing materials increased from 10 to 15[24].Duanet alinvestigated the effects of dielectric materials on the decomposition of CO2in a DBD reactor.CaO exhibited a superior promotional effect on the conversion of CO2compared a reactor with other packing materials,the conversion of CO2was highly improved in the CaO packed reactor,and reached the highest value of 32.8%,which was approximately twice the conversion of CO2in the SiO2packed reactor and the Al2O3packed reactor.The results indicated that the dielectric constant would affect the conversion of CO2[30]

    Another key factor is the thermal conductivity of the packing materials,which greatly affects the heat transfer process in the packed-bed DBD reactor.It has been reported that heat generates inside of packed region then transfers to environment and the flowing gas when the plasma is lighted in the packed-bed DBD,and heat production is proportional to the discharge power in a packed-bed DBD reactor[31].The temperature of packed region would vary depending on different thermal conductivity of the packing materials.As an exothermic reaction,low temperature is beneficial for the forward reaction for NH3synthesis(R1).Previously,Sunet alinvestigated the effects of ambient temperature and other operation parameters and reported high temperature would enhance the heat build-up effect,which is beneficial for NH3synthesis at the beginning of discharge[32].However,very few works have reported the influence of thermal conductivity of packing materials on NH3synthesis in a DBD reactor when the reaction has been in the steady state.

    In this work,several packing materials(BeO,AlN,SiO2,TiO2,BaTiO3and SiO2)with different relative dielectric constant and thermal conductivity are selected to study the effects of the dielectric constant and thermal conductivity of these packing materials on the process of plasma-assisted NH3synthesis.The effects of packing materials on the fundamental characteristics of the NH3synthesis in the DBD reactor are investigated through electrical and thermal diagnostics.We also illustrate that increasing specific input energy(SIE)or packing the materials of higher dielectric constant could enhance the reduced electric field(E/N)and change the electron energy distribution function(EEDF)by calculation through solving the Boltzmann equation using BOLSIG+.A complete understanding on the effects of dielectric constant and thermal conductivity of packing materials on NH3synthesis has been proposed.

    2.Experimental

    2.1.Experimental setup

    The DBD reactor is a coaxial cylindrical structure,as shown in figure 1.A 30 mm long stainless-steel mesh is wrapped around the quartz tube with an external diameter of 10 mm and an inner diameter of 8 mm.A stainless-steel rod with an outer diameter of 4 mm is placed in the center of the quartz tube and it acts as an inner electrode.The inner electrode is connected to the high voltage output of the power supply.An external capacitance(Cm=0.47 μF)is connected between the outer electrode and the ground.

    Zero grade(purity>99.999%)N2and H2are provided from gas cylinders and regulated by mass flow controllers(Sevenstars D07-B).Before being introduced to the reactor,the two gas streams are mixed in a mixing chamber.The DBD reactor is supplied by an AC high voltage power supply(Suman CTP-2000K).The gas feed ratio is the molar ratio of N2/H2,which is fixed at 1:1.The total gas flow rate is fixed at 100 ml·min-1.The frequency of AC discharge(f),which is kept in the range of 9.2-9.8 kHz in this work.The applied voltage(Ug)and the voltage across measuring resistance are measured by high voltage probe(Tektronix TPP0502)to obtain theU-Icurve.The voltage(Uc)on the external capacitor is measured to obtain the charge generated in the discharge.All of the electrical signals are sampled by a fourchannel digital oscilloscope(Tektronix TDS2024C).Various electrical parameters of plasma could be obtained from the typicalQ-ULissajous figure as illustrated by figure 2.

    Table 1.Relative dielectric constant and thermal conductivity of the packing materials.

    A series of packing materials with different dielectric constants and thermal conductivities include BaTiO3,TiO2,SiO2,Al2O3,AlN and BeO.These materials are nonporous ceramics with high purity(>99%)and purchased from Beilong Electronic Co.Ltd China.The relative dielectric constant and thermal conductivity of the packing materials are measured by the manufacturer and the value is given in table 1.Before using,these materials are crushed and sieved to 40-60 meshes,followed by ultrasonic cleaning in ethanol for 1 h.Then the materials were heated in an argon flow at 350°C for 2 h,and then cooled to room temperature in argon flow.After pretreatment,the packing materials are placed into the plasma region while the length of packed region is 12 mm.Quartz wool is evenly placed on both ends of the packing materials to hold them.

    The surface temperature of the reactor is measured by an IR camera(UNI-T UTi165A),while the average surface temperature of DBD reactor is calculated by UTi thermal analyzer software.The temperature of inlet and outlet gas mixtures could be measured by dual-channel digital thermometer(UNI-T UT320D).The outlet NH3concentration is analyzed by dualchannel gas chromatography(Fuli 9790Ⅱ),which is equipped with a thermal conductivity detector.The chromatography column was a packed column for ammonia separation produced by Wuhan Puli Technology.Co.Ltd China.The outlet NH3concentration is measured three times and the average value is given.The EEDF is calculated by BOLSIG+[33].

    2.2.Parameter calculation

    To better understand the effects of dielectric constant and thermal conductivity on NH3synthesis,a simplified model of the packed DBD reactor is established.The regions of packing materials and gas phase are equivalently considered as two parallel plate capacitors,as shown in figure 3.

    The relationship of relative dielectric constant and dielectric constant is as follow:

    where ε,εrand ε0are dielectric constant,relative dielectric constant and the dielectric constant of vacuum.

    The plasma discharge power(Pdis)can be calculated by:

    whereCmis external capacitance(0.47 μF),fis the frequency of AC discharge andAis the area of the Lissajous figure.

    The equivalent gap thicknesses in the gas(dg)and solid(dp)region are determined by equations(2)and(3):

    wheredgapis the gap width between internal face of quartz tube and the surface of inner electrode,lis the length of packed region andDis the outer diameter of quartz tube.

    The volume of the void of the packed region can be calculated by:

    whereVgasis the volume of the void in the packed region,Vpackingis the volume of the packed region.In this work,the value of measured void fraction coefficient(β)is measured at 0.43±0.01.

    Cdrepresents the capacitance of the quartz tube and the capacitance of the gap isCgap.The equivalent capacitance of the quartz tube is calculated according to the topology of a coaxial capacitor,as shown in equation(5)

    where ε0is the dielectric constant of vacuum(8.85×10-12F·m-1),ε is the dielectric constant of quartz tube andxis the wall thickness.

    The equivalent capacitance of the gas(Cg)and solid(Cp)regions is calculated in the form of a parallel plate capacitor:

    where εgis the relative dielectric constant of mixed gas(εg=1 when the gas feed ratio is 1:1);dganddpare the equivalent gaps of gas and packing materials respectively;Sis the surface area of outer electrode of packed region

    The relationship betweenCcellandCgapcan be calculated as follows:

    The gas voltage(Ug)can be calculated as follows:

    The average strength of electric field(E)is defined as:

    The average electric field strength of gas phase(Egas)can be calculated by:

    The breakdown voltage of DBD reactor could be calculated by:

    whereUmin,CcellandCeffcan be determined by Lissajous figure(as shown in figure 2).

    SIE is defined as the ratio of the plasma absorbed power to the gas flow rate:

    whereQis the total gas flow rate.

    The EE of the plasma-assisted NH3synthesis process is defined as follow:

    whereMis the molar mass of NH3,Coutdenotes the outlet NH3concentration andQafteris the outlet gas flow rate after the reaction and is measured by a soap film flow-meter.

    3.Results

    3.1.Effect of SIE on NH3 synthesis

    Figure 4 presents the effects of SIE on NH3synthesis over different packing materials.As presented in figure 4(a),when the SIE increases from 3 to 5.4 kJ·l-1,the NH3concentration increases from 112 to 1344 ppm in the BaTiO3packed reactor,and the increment of NH3concentration is only 283 ppm when the SIE increases in the range of 5.4-10.8 kJ·l-1.Meanwhile,the increment of NH3concentration in the TiO2packed reactor is 1019 ppm when the SIE increases from 3 to 5.4 kJ·l-1,which is 351 ppm when the SIE increases from 5.4 to 10.8 kJ·l-1.As shown in figure 4(b),for the case of packing BeO,the concentration increases from 478 to 1707 ppm with the SIE increasing in the range of 3-10.8 kJ·l-1.

    Figure 1.Experimental setup for the plasma-assisted NH3 synthesis.

    Figure 2.Schematic figure for the derivation of electrical parameters from a Q-U Lissajous.

    Figure 3.The simplified model of packed-bed DBD reactor.

    Figure 4.Effects of SIE on NH3 concentration over different packing materials.

    Figure 5.Effects of SIE on energy efficiency over different packing materials.

    In summary,higher NH3concentration is achieved over BaTiO3and BeO peaked reactors.The NH3concentration increases rapidly,but when the SIE is higher than a certain value(the optimal value),the trend of NH3concentration changes from rapid increase to slight increase with the increasing SIE,this trend is pronounced as shown in figure 4(a).

    Figure 5 shows the effect of SIE on EE in the presence of different packing materials.As shown in figure 5(a),the highest EEs of 0.68,0.56 and 0.33 g·kWh-1for BaTiO3,TiO2and SiO2packed reactors,respectively,are observed at the SIE of 5.4 kJ·l-1.Meanwhile,the EE is always higher in presence of BaTiO3,compared to that of TiO2and SiO2.As presented in figure 5(b),the EE follows the ascending order of Al2O3<AlN<BeO,indicating higher thermal conductivity of packing materials is beneficial for NH3synthesis.Therefore,it can be deduced that higher EE is achieved by employing packing materials with higher thermal conductivity or higher dielectric constant.

    Figure 6.Relationships between voltage and current for different packed materials reactors(SIE=5.4 kJ·l-1).

    For all of the cases,the EE increases initially,however,when the SIE is higher than the optimal value,the energy efficient would decrease.As shown in figure 5(a),the EE decreases dramatically when the SIE is in the range of 5.4 and 10.8 kJ·l-1,which decreases from 0.68 to 0.41 g·kWh-1in presence of BaTiO3and from 0.56 to 0.37 g·kWh-1in presence of TiO2.Similar phenomenon is also presented in figure 5(b).This observation may be related to electronimpact NH3decomposition and the inhibition of NH3synthesis by high temperature in the packed region.

    3.2.Discharge behaviors in different materials packed reactors

    Figure 6 shows the current profiles in different materials packed reactors.In packed-bed DBD reactors,the discharge would transfer from filamentary discharge to a combination of surface discharge on the surface of the packing materials and spatially limited microdischarges generated in the void space[34,35].At the SIE of 5.4 kJ·l-1,the intensity of time-averaged microdischarges current in BaTiO3packed reactor is 18.3 mA,which is about 2 times and 1.8 times higher than those of SiO2and TiO2,respectively.This phenomenon was also reported in various plasma-assisted processes[35-37].Obviously,the current intensity increases in presence of the materials with higher dielectric constant.Moreover,the breakdown voltage is increased in presence of packing materials with higher dielectric constant as determined from both Lissajous figures and current profiles.Compared to that of BaTiO3(εr=7200),the breakdown voltage decreases from 1.15 to 0.77 kV in presence of SiO2(εr=3.80).

    Figure 7.Q-U Lissajous figures for different materials packed reactors(SIE=5.4 kJ·l-1).

    The time-averaged microdischarge current intensity is approximately equal when Al2O3,AlN and BeO are employed as shown in figure 6(b).The current intensities in BeO and AlN packed reactors are 9.2 mA and 9.1 mA,respectively,which are close to the case of Al2O3packed reactor(8.8 mA).A slight increasing number of microdischarge is obtained in the presence of BeO and AlN,which may be related to the increasing electron density at lower temperature of discharge region.

    Figure 7 shows the Lissajous figures for different materials packed reactors at the SIE of 5.4 kJ·l-1.The shape of the Lissajous figure is varied when different packing materials are employed,which suggests the changes in the discharge characteristics includingVpk,Ceff,andQdischarge(quantity of the discharged charges per half-cycle).At the SIE of 5.4 kJ·l-1,the applied voltage of the DBD reactor increases from about 13.6 kVpk-pkfor SiO2packed reactor to about 15.8 kVpk-pkfor the case of BaTiO3.TheQdischargeincreases in the presence of packing materials with higher dielectric constant,which is 0.32 μC for BaTiO3packed reactor and 0.23 μC for SiO2packed reactor.Meanwhile,higherCeffof 88.39 pF is achieved in BaTiO3packed reactor,which is 11.8% higher than that of SiO2.The large effective capacitance suggests that the addition of higher dielectric constant materials benefits to achieve the fully-bridged discharge.The higherCeffproves that the employment of BaTiO3could bridge the gap between the electrodes and enhance the charge transfer between them[38].In contrast,the Lissajous figures for BeO,AlN and Al2O3packed reactors are approximately the same,indicating the employment of these materials is of no effect in varying electrical characteristics of plasma.

    Figure 8.(a)Surface temperature distribution of SiO2 packed reactor at different SIE(discharge time=1.5 min),(b)reactor surface temperature distribution at SIE of 5.4 kJ·l-1(discharge time=5 min,1-6 reactors are packed with BaTiO3,TiO2,SiO2,Al2O3,AlN and BeO,respectively).

    3.3.Surface temperature distribution of the reactors

    The surface temperature of plasma region in steady state is notably affected by the SIE as shown in figure 8(a).When the SIE increases from 3 to 10.8 kJ·l-1,the average reactor surface temperature of packed region(inside of the dotted box)increases from 19.1°C to 78.3°C,which increases 4.1 times when the SIE increases in the range of 3-10.8 kJ·l-1.Similar observations were also reported by Kimet aland Tirumalaet al[31,39].We find the average surface temperature of the packed region is significantly higher than other area covered by the outer electrode,indicating that the presence of packing materials contributed to generate intensified microdischarge.Figure 8(b)presents the IR imagings of the packed region of the plasma reactor in the presence of different packing materials at the SIE of 5.4 kJ·l-1.The average surface temperature of the packed region is obtained from figure 8(b)and summarized in figure 9.The reactors packed with BaTiO3,TiO2and SiO2are of similar average surface temperature of the packed region,which is 96.2 °C,91.2 °C and 90.6 °C,respectively.In contrast,average surface temperature of the reactor packed region decreases dramatically in presence of packing materials with higher thermal conductivity.Compared to the Al2O3packed reactor,the presence of BeO leads to a decrease of average surface temperature from 88.4 °C to 70.3 °C,this suggests that the heat transfer is significantly enhanced by employing packing materials with higher thermal conductivity.

    Figure 9.Relationships between NH3 concentration and average surface temperatures in steady state with different packing materials(discharge time=5 min,SIE=5.4 kJ·l-1).

    Figure 10.(a)Mean energy as a function of reduced electron field strength(E/N),(b)electron energy distribution function.

    4.Discussions

    4.1.Correlations of packing materials with higher dielectric constant to the ammonia synthesis

    The NH3concentration and EE in the presence of BaTiO3were 1344 ppm and 0.68 g·kWh-1,respectively,which were about 106% higher than those of SiO2,respectively,at the SIE of 5.4 kJ·l-1.Meanwhile,electrical characteristics of DBD reactors changed dramatically in presence of packing materials with different dielectric constants.As described in figure 7,theCeffof reactor,as well asQdischarge,was in an ascending order of SiO2<TiO2<BaTiO3,indicating that the presence of BaTiO3leads to an expansion of microdischarge in the gas phase.Thus,more energy could be injected into the plasma for chemical reactions[38,40].

    According to our calculation,Egaswould be improved by the higher dielectric constant of packing materials.To further understand the relationships between electron energy and reduced electric field strength,EEDF and mean electron energy are calculated by BOLSIG+.As shown in figure 10,with increasingE/N,the mean electron energy is improved as well.Meanwhile,the EEDF is found to shift with an increase in electron density in the high-energy tail of the distribution function,indicating the increase fraction of highly energetic electrons.In our case,at the SIE of 5.4 kJ·l-1,theE/Nis 1.62 times higher in the presence of BaTiO3,compared to that of SiO2.The higherE/Ndirectly results in the increase of mean electron energy,which is about 1.41 times while the proportion of energetic electrons would also increase as confirmed in figure 10(b).Similar phenomenon was also reported by Meiet alin CO2decomposition experiment,at the discharge power of 50 W,the reduced electric filed strength generated by packing BaTiO3was 1.37 times higher than which was generated by SiO2packed and 1.53 times to the non-packed reactor.Compared to that of the SiO2and nonpacked reactors,The EE of the BaTiO3packed reactor was 1.75 times and 1.35 times higher than those of the SiO2and non-packed reactors,respectively[35].It was also reported by Chenet al[29]and Duanet al[30]that higher electric field and consequently higher mean electron energy would be obtained in presence of packing materials with higher dielectric constant in a packed-bed DBD reactor.

    Since the quantity of energetic electrons and energy injected into plasma were improved in presence of materials with higher dielectric constant,electron-impact reactions would be enhanced and consequently contribute to the generation of chemically active species including excited species likeand(e.g.(R2)-(R5)).These chemically reactive species would recombine and form NHx(x=1,2)species,which were regarded as the most important intermediate products of plasma-assisted NH3synthesis process[23],and accelerate NH3formation(e.g.(R4)).

    Hence,the NH3concentration,as well as EE,was observed to follow the ascending order of dielectric constant of packing materials,i.e.SiO2<TiO2<BaTiO3.However,in plasma-assisted NH3synthesis process,the positive reactions towards NH3formation would coexist with the reverse reactions of NH3decomposition.Thus,NH3decomposition is inevitable during the whole process and unfavorable for NH3formation.These reverse reactions would cause the decomposition of NHx·(x=1,2 and 3)species and directly lead to the decrease of EE,which would be exacerbated in the presence of materials with higher dielectric constant.In particular,when the SIE is higher than the optimal value,e.g.5.4 kJ·l-1for BaTiO3,TiO2and SiO2packed reactor in this work,the EE decreases dramatically in presence of BaTiO3,compared to that of SiO2.This phenomenon is aligned with many works[7,41,42].

    Since the dielectric constant of Al2O3,BeO and AlN was similar,an increasing amount of microdischarge were generated when materials with higher thermal conductivity were packed at the same SIE.This phenomenon was related to the decrease of average temperature in discharge zone,which would contribute to improve the electron density,promote the formation of filamentary microdischarge and further contribute to the electron-impact reactions for NH3synthesis[37,43].

    4.2.Correlations of packing materials with higher thermal conductivity to the ammonia synthesis

    It has been reported the NH3concentration would increase with the increasing discharge time and eventually stabilize,then the maximum NH3concentration would be achieved[7,32].At the beginning of plasma discharge duration,heat storage inside the packed region is conducive to accelerate the synthesis reaction to reach the steady state[44].However,NH3synthesis reaction(R1)is an exothermic reaction.According to Le Chatelier’s principle,lower ambient temperature is conductive to the NH3synthesis process.To further understand the effect of thermal conductivity on NH3synthesis,our discussion is carried out under the condition that NH3synthesis reaction has been in the steady state.

    Since the NH3synthesis process is in steady state,the heat generation and dissipation in the whole reactor would reach thermodynamic equilibrium.Meanwhile,both inlet and outlet gas flow rates are fixed,and the pressure inside the reactor is also a constant(atmosphere pressure).The reactor is simplified as a thermodynamic steady flow system with a stable temperature field.

    According to the second low of thermodynamic,the entropy increment can be calculated as follows:

    whereSis entropy of the plasma region(J mol-1·K-1),Cvdenotes specific heat capacity at constant volume of the mixed gas(J kg-1·K-1)andTis the average temperature of the packed region(°C).

    The relationship between the NH3generation and the entropy increment is as follows:

    whereJis the relative reaction react rate of per molar NH3generation(mol ml·min-1);Ais the relative chemical entropy increase of per molar NH3generation(J mol·K-1).

    According to the simplified model,the heat transform quantity(Qtrans)could be calculated by:

    where λtis the total thermal conductivity of solid phase and gas phase(W m-1·K-1),Aisois the area of the isothermal surface of packed region(m2),tis the temperature of isothermal surface(°C),σ is the distance from isothermal surface to reactor outer surface(m),τ is the heat transfer time(s).

    In plasma-involved cases,the energy balance equation can be described as follows:

    wherePreactionis the power consumption of the reactions andPgasis the power absorbed by the gases.

    The power consumption of reactions(Preaction)can be calculated as:

    whereQmis the mass flow rate of the mixed gas,αiis the mass fraction of the one gas components,hiis the specific enthalpy of one of the gas components,imeans one of the gas components(e.g.N2,H2and NH3)andgis the total species number of the gas components.

    The discharge power absorbed by feed gas can be calculated by:

    whereTis the temperature of gas,Cp(i)is the heat capacity at constant pressure of the gas.

    In this work,the presence of BeO achieved both a higher NH3concentration of 2102 ppm and EE of 1.02 g·kWh-1,which were about 30%higher than that of Al2O3.Meanwhile,the heat transfer process would be enhanced in the presence of materials of higher thermal conductivity.At the SIE of 5.4 kJ·l-1,the heat transfer quantity(Qtrans)was increased by 32.5% in the BeO packed reactor compared to the reactor packed with Al2O3.As a result,the addition of BeO significantly decreases the average temperature to 70.3 °C.Therefore,with the enhanced heat transfer process,compared to that of Al2O3,a higher entropy increment would be achieved with BeO,according to equations(15)-(18).

    NH3synthesis reaction is a process of entropy reduction.According to the second law of thermodynamics,this process is spontaneous and needs an extra entropy increase process with the aim of increasing the system entropy to benefit the forward reaction of(R1).In plasma-assisted process,the processes of gas molecules electron-impact decomposition and the heat transfer process were considered as the entropyincreasing process[45,46].Considering that the temperature difference between the inlet and outlet gases of the reactor was small(<0.1 °C)and could be neglected in this work,furthermore,the ambient temperature was almost the same,it would be thermochemically feasible that the process which heat transfers to the environment was the dominate process of entropy increase of the plasma-assisted NH3synthesis process.Therefore,compared with that of Al2O3,both higher NH3concentration and EE were achieved in presence of BeO because of both the entropy increment and lower temperature of the packed region,which contributed to NH3synthesis reaction.

    5.Conclusions

    Plasma-assisted NH3synthesis has been carried out in different materials(BeO,AlN,SiO2,TiO2,BaTiO3and SiO2)packed DBD reactors.The effects of the dielectric constant and thermal conductivity of packing materials on NH3synthesis have been investigated.

    The presence of materials with higher dielectric constant materials could improve the electron energy and energy injected to plasma,which would be beneficial for the generation of chemically active species by electron-impact reactions and accelerate NH3formation.Therefore,both NH3concentration and EE were following the ascending trend of SiO2<TiO2<BaTiO3.The higher NH3concentration of 1344 ppm was achieved in presence of BaTiO3,followed by that of TiO2of 1106 ppm and SiO2of 654 ppm.Meanwhile,the NH3synthesis EE of 0.68 g·kWh-1was achieved with BaTiO3,which was 105%higher than that of SiO2at the SIE of 5.4 kJ·l-1.

    The presence of materials with higher thermal conductivity significantly enhanced the heat transfer process and decreased the average temperature of packed region.Compared to that of Al2O3,the quantity of heat transfer was increased by 32.5%in presence of BeO.The entropy of the packed region was also increased as result.Meanwhile,the lower temperature of the packed region(70.3°C)was also obtained in the BeO packed reactor,which was 20.5%lower than that of Al2O3,at the SIE of 5.4 kJ·l-1.Both the entropy increment and the lower temperature have proved to be conductive to forward reaction of NH3synthesis.Therefore,the presence of materials of higher thermal conductivity achieved higher NH3concentration and EE,which increased with different packing materials in the order of Al2O3<AlN<BeO.

    Acknowledgments

    The authors are grateful for the financial support from National Natural Science Foundation of China(No.51976093)and K C Wong Magna Fund in Ningbo University.

    ORCID iDs

    猜你喜歡
    學(xué)理
    論中國民族樂器“技法母語”的學(xué)理內(nèi)涵
    “歌劇學(xué)”學(xué)科建構(gòu)的學(xué)理闡釋
    冰雪難美性項群跨項選材的學(xué)理基礎(chǔ)、成功案例及實(shí)踐路徑
    政府信息依申請公開收費(fèi)的反思——對新修訂《政府信息公開條例》的學(xué)理支撐
    法大研究生(2019年1期)2019-11-16 00:36:58
    學(xué)理審思:真實(shí)情境寫作之中考命題
    司法義務(wù)視域下指導(dǎo)性案例的學(xué)理分析
    雙回單訓(xùn)釋教學(xué)新法的學(xué)理依據(jù)
    “新常態(tài)”概念的演進(jìn)與學(xué)理支撐
    張學(xué)理水墨漫畫選
    西漢道義觀的學(xué)理形成
    古代文明(2015年3期)2015-07-07 07:31:41
    纵有疾风起免费观看全集完整版| 国产淫片久久久久久久久| 久热久热在线精品观看| 国产探花极品一区二区| 国产午夜精品久久久久久一区二区三区| 麻豆国产97在线/欧美| 能在线免费看毛片的网站| 亚洲精品中文字幕在线视频 | 久久久精品免费免费高清| 高清午夜精品一区二区三区| 久久99热这里只有精品18| 97在线视频观看| 国产精品熟女久久久久浪| 亚洲国产精品国产精品| 亚洲成色77777| 精华霜和精华液先用哪个| 偷拍熟女少妇极品色| 人妻少妇偷人精品九色| 小蜜桃在线观看免费完整版高清| 亚洲av中文av极速乱| 99久久中文字幕三级久久日本| 91精品伊人久久大香线蕉| 妹子高潮喷水视频| 国产日韩欧美在线精品| 国产高潮美女av| 丰满人妻一区二区三区视频av| 国产精品久久久久久av不卡| 欧美极品一区二区三区四区| 亚洲av二区三区四区| 精品国产乱码久久久久久小说| 天堂8中文在线网| 秋霞在线观看毛片| 高清视频免费观看一区二区| 只有这里有精品99| 最近中文字幕高清免费大全6| 三级国产精品片| 国产在线男女| 国产永久视频网站| 国产精品爽爽va在线观看网站| 国产美女午夜福利| 91狼人影院| 一区二区三区精品91| 亚洲精品国产成人久久av| 国产精品一区www在线观看| 夜夜爽夜夜爽视频| 亚洲精品自拍成人| 亚洲av电影在线观看一区二区三区| av女优亚洲男人天堂| 久热久热在线精品观看| 久久精品国产亚洲av天美| 精品久久久噜噜| 18+在线观看网站| 国产男女超爽视频在线观看| 插逼视频在线观看| 交换朋友夫妻互换小说| 国产精品久久久久久久电影| 不卡视频在线观看欧美| 国产一区有黄有色的免费视频| 日本黄色片子视频| 涩涩av久久男人的天堂| 亚洲在久久综合| 亚洲第一区二区三区不卡| 国产真实伦视频高清在线观看| 欧美xxⅹ黑人| 久久久久国产网址| 精品视频人人做人人爽| 亚洲人成网站在线播| 国产成人精品一,二区| 亚洲精品乱码久久久v下载方式| 精品午夜福利在线看| 亚洲国产日韩一区二区| 2021少妇久久久久久久久久久| 自拍偷自拍亚洲精品老妇| 国产大屁股一区二区在线视频| av一本久久久久| 久久久久久久久久人人人人人人| 久久精品夜色国产| 久久久成人免费电影| 蜜桃在线观看..| 一本久久精品| 精品人妻视频免费看| 一级爰片在线观看| 精品少妇久久久久久888优播| 亚洲欧美一区二区三区黑人 | 亚洲人成网站高清观看| 观看美女的网站| 亚洲av日韩在线播放| 久久久久久久大尺度免费视频| 亚洲av福利一区| 黄色配什么色好看| 黄色视频在线播放观看不卡| 国内揄拍国产精品人妻在线| 秋霞伦理黄片| videossex国产| 色吧在线观看| 日韩成人伦理影院| 中文字幕精品免费在线观看视频 | 亚洲精华国产精华液的使用体验| 毛片一级片免费看久久久久| 91午夜精品亚洲一区二区三区| 伊人久久国产一区二区| 日韩成人av中文字幕在线观看| 永久免费av网站大全| 晚上一个人看的免费电影| 夫妻午夜视频| 国产熟女欧美一区二区| 欧美成人精品欧美一级黄| 多毛熟女@视频| 日日摸夜夜添夜夜爱| 国产精品欧美亚洲77777| 亚洲国产精品一区三区| 亚洲aⅴ乱码一区二区在线播放| 中文欧美无线码| 美女中出高潮动态图| 久久韩国三级中文字幕| 亚洲av日韩在线播放| 国产深夜福利视频在线观看| 观看av在线不卡| 国产伦精品一区二区三区四那| 搡老乐熟女国产| 国产免费福利视频在线观看| 中文乱码字字幕精品一区二区三区| 久久99精品国语久久久| 精品久久久噜噜| 美女内射精品一级片tv| 国产高潮美女av| 国产 一区 欧美 日韩| 丰满迷人的少妇在线观看| 制服丝袜香蕉在线| 国产乱人偷精品视频| 亚洲国产欧美人成| 国产av码专区亚洲av| 99热全是精品| 国产亚洲欧美精品永久| www.av在线官网国产| 三级国产精品片| 一级毛片电影观看| 日本-黄色视频高清免费观看| 久久久久久久大尺度免费视频| 一级二级三级毛片免费看| 在线观看av片永久免费下载| 国产色婷婷99| 日本欧美国产在线视频| 少妇人妻 视频| 精品酒店卫生间| 欧美成人精品欧美一级黄| 欧美日韩视频高清一区二区三区二| 久久热精品热| 精品熟女少妇av免费看| 一区二区av电影网| 久久久久久伊人网av| 日韩不卡一区二区三区视频在线| 不卡视频在线观看欧美| av在线老鸭窝| 国产乱人偷精品视频| 老司机影院成人| 一级片'在线观看视频| 亚洲欧美一区二区三区黑人 | 欧美区成人在线视频| 欧美高清性xxxxhd video| h日本视频在线播放| 色视频www国产| 精品99又大又爽又粗少妇毛片| 日韩欧美一区视频在线观看 | av黄色大香蕉| 国产深夜福利视频在线观看| 亚洲av在线观看美女高潮| 丝瓜视频免费看黄片| 99热这里只有精品一区| 国模一区二区三区四区视频| 日日撸夜夜添| 精品熟女少妇av免费看| 国产日韩欧美在线精品| 国产片特级美女逼逼视频| 自拍欧美九色日韩亚洲蝌蚪91 | 成人午夜精彩视频在线观看| 少妇人妻精品综合一区二区| 免费大片18禁| 国内揄拍国产精品人妻在线| 欧美日韩综合久久久久久| 亚洲人成网站高清观看| 丝袜脚勾引网站| 亚洲精品国产色婷婷电影| 久久99精品国语久久久| www.av在线官网国产| 涩涩av久久男人的天堂| 久热这里只有精品99| 高清不卡的av网站| 只有这里有精品99| 少妇高潮的动态图| 国产 一区 欧美 日韩| 五月天丁香电影| 麻豆国产97在线/欧美| 久久99精品国语久久久| 伦理电影大哥的女人| 国产高清三级在线| 在线免费十八禁| 久久综合国产亚洲精品| 成人国产麻豆网| 午夜福利在线观看免费完整高清在| 成人综合一区亚洲| 中文乱码字字幕精品一区二区三区| xxx大片免费视频| 亚洲精品自拍成人| 久久久久国产网址| 久久久a久久爽久久v久久| 国产精品久久久久久精品电影小说 | 成人黄色视频免费在线看| 久久久久久久久久久免费av| 中文字幕人妻熟人妻熟丝袜美| 男男h啪啪无遮挡| 大码成人一级视频| 国产亚洲5aaaaa淫片| 精品少妇黑人巨大在线播放| 韩国高清视频一区二区三区| 天美传媒精品一区二区| 99久久中文字幕三级久久日本| 22中文网久久字幕| 国产毛片在线视频| 在线观看人妻少妇| 欧美xxxx性猛交bbbb| 国产成人免费无遮挡视频| 一本—道久久a久久精品蜜桃钙片| 国产精品一区二区性色av| 久久精品熟女亚洲av麻豆精品| 秋霞在线观看毛片| 久久午夜福利片| 国产亚洲午夜精品一区二区久久| 最黄视频免费看| 永久免费av网站大全| 亚洲国产日韩一区二区| 欧美一级a爱片免费观看看| 妹子高潮喷水视频| 亚洲精品国产av成人精品| 噜噜噜噜噜久久久久久91| 草草在线视频免费看| 色网站视频免费| 一级毛片aaaaaa免费看小| 中文字幕亚洲精品专区| 在现免费观看毛片| 人人妻人人添人人爽欧美一区卜 | 成人特级av手机在线观看| 99热国产这里只有精品6| 精品人妻偷拍中文字幕| 三级经典国产精品| 欧美高清成人免费视频www| 精品亚洲乱码少妇综合久久| 久久精品夜色国产| 一级毛片我不卡| 国产综合精华液| av免费在线看不卡| 亚洲美女黄色视频免费看| 午夜福利在线观看免费完整高清在| 国产成人freesex在线| h视频一区二区三区| 亚洲国产精品国产精品| 国产视频内射| 在线 av 中文字幕| 国产成人a∨麻豆精品| 久久精品国产a三级三级三级| 99re6热这里在线精品视频| 成年美女黄网站色视频大全免费 | 久久久久国产精品人妻一区二区| 精品99又大又爽又粗少妇毛片| 久久人人爽人人片av| 日韩亚洲欧美综合| 18禁在线无遮挡免费观看视频| 视频中文字幕在线观看| 国产精品三级大全| 黑人高潮一二区| 日本wwww免费看| 国产真实伦视频高清在线观看| 午夜福利视频精品| 人人妻人人添人人爽欧美一区卜 | 久久6这里有精品| 少妇裸体淫交视频免费看高清| 大香蕉97超碰在线| 国产精品伦人一区二区| 成人综合一区亚洲| 久久久久久久久久久免费av| 亚洲精品乱码久久久久久按摩| 日本wwww免费看| 国产精品精品国产色婷婷| 国产精品一区二区在线不卡| 国国产精品蜜臀av免费| 大片电影免费在线观看免费| av免费在线看不卡| www.色视频.com| 在线观看一区二区三区| 亚洲精品日韩av片在线观看| 久久99蜜桃精品久久| 汤姆久久久久久久影院中文字幕| 男的添女的下面高潮视频| 国产精品一区二区性色av| 欧美三级亚洲精品| 国产精品久久久久久久电影| 老司机影院毛片| 三级经典国产精品| 女性生殖器流出的白浆| 亚洲最大成人中文| 在线精品无人区一区二区三 | 久久久成人免费电影| 3wmmmm亚洲av在线观看| 哪个播放器可以免费观看大片| 中文欧美无线码| 亚洲精品乱码久久久久久按摩| av网站免费在线观看视频| 1000部很黄的大片| 男女免费视频国产| 我的老师免费观看完整版| 天堂中文最新版在线下载| 欧美三级亚洲精品| 制服丝袜香蕉在线| 身体一侧抽搐| 成人毛片60女人毛片免费| 两个人的视频大全免费| 男女免费视频国产| 国产免费一级a男人的天堂| 一级毛片黄色毛片免费观看视频| 欧美区成人在线视频| 全区人妻精品视频| 亚洲国产高清在线一区二区三| 你懂的网址亚洲精品在线观看| 最近中文字幕2019免费版| 亚洲电影在线观看av| 国产免费福利视频在线观看| 午夜福利网站1000一区二区三区| 大话2 男鬼变身卡| 男的添女的下面高潮视频| 男女啪啪激烈高潮av片| 日日啪夜夜撸| 日本色播在线视频| 少妇丰满av| 久久99热这里只有精品18| 七月丁香在线播放| 亚洲自偷自拍三级| 精品久久国产蜜桃| 欧美性感艳星| 91久久精品国产一区二区成人| 国产成人a∨麻豆精品| 亚洲国产最新在线播放| av不卡在线播放| 人妻少妇偷人精品九色| 日韩一区二区三区影片| 国产精品久久久久久久久免| 国产欧美日韩一区二区三区在线 | 成人亚洲精品一区在线观看 | 欧美亚洲 丝袜 人妻 在线| 亚洲欧美日韩无卡精品| 久久精品国产自在天天线| 亚洲综合色惰| 中文精品一卡2卡3卡4更新| 夫妻午夜视频| 日本av手机在线免费观看| 91久久精品国产一区二区成人| 一级毛片我不卡| 人妻少妇偷人精品九色| 久久99热6这里只有精品| 亚洲人成网站在线播| 亚洲国产毛片av蜜桃av| 久久这里有精品视频免费| 99re6热这里在线精品视频| 亚洲经典国产精华液单| 男女免费视频国产| 免费人妻精品一区二区三区视频| 少妇被粗大猛烈的视频| 亚洲av男天堂| 多毛熟女@视频| 一区二区三区四区激情视频| 嫩草影院入口| 黄色视频在线播放观看不卡| 中文在线观看免费www的网站| 一个人看的www免费观看视频| 日本-黄色视频高清免费观看| 亚洲精品国产成人久久av| 一区二区三区四区激情视频| 我要看黄色一级片免费的| 少妇熟女欧美另类| 国产在线一区二区三区精| av专区在线播放| 久久久久国产网址| 身体一侧抽搐| 日本色播在线视频| xxx大片免费视频| 视频中文字幕在线观看| 一区二区三区免费毛片| 色婷婷久久久亚洲欧美| 少妇的逼水好多| 久久久久久伊人网av| 亚洲国产精品成人久久小说| 夜夜爽夜夜爽视频| 亚洲欧美精品自产自拍| 日韩制服骚丝袜av| 春色校园在线视频观看| 一级爰片在线观看| 深夜a级毛片| 久久人人爽人人爽人人片va| 午夜福利视频精品| 麻豆乱淫一区二区| 老师上课跳d突然被开到最大视频| 欧美日本视频| 亚洲第一区二区三区不卡| 日本午夜av视频| 亚洲欧洲日产国产| 丰满少妇做爰视频| 亚洲精品中文字幕在线视频 | 国产男女内射视频| 久久 成人 亚洲| 99热6这里只有精品| 最近最新中文字幕大全电影3| 国产永久视频网站| 卡戴珊不雅视频在线播放| 九九爱精品视频在线观看| 在线天堂最新版资源| 国产男女内射视频| 国产乱来视频区| 九草在线视频观看| 我的女老师完整版在线观看| 成年人午夜在线观看视频| 日韩,欧美,国产一区二区三区| 2022亚洲国产成人精品| av又黄又爽大尺度在线免费看| 一二三四中文在线观看免费高清| 99热全是精品| 成年女人在线观看亚洲视频| av福利片在线观看| 哪个播放器可以免费观看大片| 丰满乱子伦码专区| 一本一本综合久久| 亚洲国产精品国产精品| 久久精品国产亚洲av涩爱| 这个男人来自地球电影免费观看 | 女性生殖器流出的白浆| 一区二区三区精品91| 免费av中文字幕在线| 亚洲av免费高清在线观看| 精品99又大又爽又粗少妇毛片| 欧美老熟妇乱子伦牲交| 国产高潮美女av| 亚洲一级一片aⅴ在线观看| 久久97久久精品| 婷婷色av中文字幕| 小蜜桃在线观看免费完整版高清| 男女免费视频国产| 日日撸夜夜添| 久久 成人 亚洲| 如何舔出高潮| 熟女av电影| 亚洲精品aⅴ在线观看| 我的女老师完整版在线观看| 久久6这里有精品| 嘟嘟电影网在线观看| 欧美日本视频| freevideosex欧美| 亚洲欧洲日产国产| 人妻少妇偷人精品九色| 联通29元200g的流量卡| 一级毛片黄色毛片免费观看视频| 久久99热这里只有精品18| 久久 成人 亚洲| 欧美极品一区二区三区四区| 精品少妇黑人巨大在线播放| 街头女战士在线观看网站| 日韩强制内射视频| 国产成人aa在线观看| 成年女人在线观看亚洲视频| 寂寞人妻少妇视频99o| 2018国产大陆天天弄谢| 精品人妻熟女av久视频| 成人综合一区亚洲| 国产免费福利视频在线观看| 黑人高潮一二区| 99久国产av精品国产电影| av在线老鸭窝| 亚洲欧美一区二区三区黑人 | 久久久久久久精品精品| 妹子高潮喷水视频| 91狼人影院| 一区二区三区精品91| 久久影院123| 国产精品伦人一区二区| 久久久久国产网址| 日韩精品有码人妻一区| 在线观看免费高清a一片| 国产深夜福利视频在线观看| 亚洲成人av在线免费| 在线看a的网站| av又黄又爽大尺度在线免费看| 国产av精品麻豆| 国产精品福利在线免费观看| 国产成人精品久久久久久| 久久人人爽人人爽人人片va| 少妇熟女欧美另类| 寂寞人妻少妇视频99o| 99九九线精品视频在线观看视频| 中文在线观看免费www的网站| 精品人妻视频免费看| 亚洲婷婷狠狠爱综合网| 欧美97在线视频| 久久久午夜欧美精品| 一区二区三区乱码不卡18| 一个人看视频在线观看www免费| 日韩三级伦理在线观看| 极品教师在线视频| 在线观看av片永久免费下载| 国产av精品麻豆| 精品久久久久久电影网| 婷婷色av中文字幕| 男女免费视频国产| 毛片一级片免费看久久久久| 亚洲,欧美,日韩| 多毛熟女@视频| 国产爽快片一区二区三区| 国产女主播在线喷水免费视频网站| 亚洲欧美日韩卡通动漫| 国产伦精品一区二区三区四那| 一级毛片 在线播放| 三级国产精品片| 日本-黄色视频高清免费观看| 欧美zozozo另类| 亚洲人成网站在线观看播放| 日韩中文字幕视频在线看片 | av专区在线播放| 一本—道久久a久久精品蜜桃钙片| 夜夜看夜夜爽夜夜摸| 在线播放无遮挡| av天堂中文字幕网| 国产亚洲91精品色在线| 久久精品国产自在天天线| 日韩欧美精品免费久久| 国产高清三级在线| 91久久精品电影网| 中文字幕精品免费在线观看视频 | 欧美老熟妇乱子伦牲交| 日韩中文字幕视频在线看片 | 午夜福利在线观看免费完整高清在| 久久久久久久精品精品| 新久久久久国产一级毛片| 亚洲精品日韩在线中文字幕| 亚洲精品日本国产第一区| 国产一区二区三区综合在线观看 | 亚洲精品乱久久久久久| 晚上一个人看的免费电影| 中文天堂在线官网| 欧美日韩国产mv在线观看视频 | 黑人猛操日本美女一级片| 又爽又黄a免费视频| 日产精品乱码卡一卡2卡三| 国产亚洲最大av| 最近最新中文字幕大全电影3| av在线观看视频网站免费| 2021少妇久久久久久久久久久| 亚洲欧美日韩东京热| 狂野欧美激情性xxxx在线观看| 极品少妇高潮喷水抽搐| 亚洲精品乱码久久久久久按摩| 国产成人a∨麻豆精品| 一个人看视频在线观看www免费| 97在线视频观看| 成人一区二区视频在线观看| 美女主播在线视频| 人妻少妇偷人精品九色| 老司机影院毛片| 啦啦啦视频在线资源免费观看| 人妻一区二区av| 啦啦啦在线观看免费高清www| 狠狠精品人妻久久久久久综合| 国产av码专区亚洲av| 日本av手机在线免费观看| 免费人成在线观看视频色| 国产成人91sexporn| a级毛色黄片| 午夜福利在线观看免费完整高清在| 精品久久久久久久末码| 国产毛片在线视频| 成人影院久久| 久久婷婷青草| 99久久中文字幕三级久久日本| 一级av片app| 国产极品天堂在线| 下体分泌物呈黄色| 永久网站在线| 一二三四中文在线观看免费高清| 亚洲在久久综合| 舔av片在线| 黄色配什么色好看| 嫩草影院入口| 国产av码专区亚洲av| 亚洲国产精品成人久久小说| 男女啪啪激烈高潮av片| 丝袜脚勾引网站| av天堂中文字幕网| 97热精品久久久久久| a级毛色黄片| av网站免费在线观看视频| 久久久色成人| 午夜福利在线观看免费完整高清在| 春色校园在线视频观看| 老司机影院毛片| 日韩欧美精品免费久久| av天堂中文字幕网| 亚洲国产最新在线播放| www.av在线官网国产| 免费大片18禁| 熟女av电影| 在线观看国产h片| 亚洲欧美成人精品一区二区| 有码 亚洲区| 日本欧美视频一区| 日本wwww免费看| 青春草国产在线视频| 久久毛片免费看一区二区三区| 国产黄频视频在线观看| 国精品久久久久久国模美| 国产在线一区二区三区精|