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

    Effects of O2 addition on the plasma uniformity and reactivity of Ar DBD excited by ns pulsed and AC power supplies

    2022-06-01 07:55:32FengLIU劉峰YueZHUANG莊越YuleiZHAO趙昱雷JieCHEN陳潔andZhiFANG方志
    Plasma Science and Technology 2022年5期
    關(guān)鍵詞:陳潔劉峰方志

    Feng LIU(劉峰),Yue ZHUANG(莊越),Yulei ZHAO(趙昱雷),Jie CHEN(陳潔)and Zhi FANG (方志)

    College of Electrical Engineering and Control Science,Nanjing Tech University,Nanjing 211816,People’s Republic of China

    Abstract Dielectric barrier discharges(DBDs)have been widely used in ozone synthesis,materials surface treatment, and plasma medicine for their advantages of uniform discharge and high plasmachemical reactivity.To improve the reactivity of DBDs, in this work, the O2 is added into Ar nanosecond (ns) pulsed and AC DBDs.The uniformity and discharge characteristics of Ar ns pulsed and AC DBDs with different O2 contents are investigated with optical and electrical diagnosis methods.The DBD uniformity is quantitatively analyzed by gray value standard deviation method.The electrical parameters are extracted from voltage and current waveforms separation to characterize the discharge processes and calculate electron density ne.The optical emission spectroscopy is measured to show the plasma reactivity and calculate the trend of electron temperature Te with the ratio of two emission lines.It is found that the ns pulsed DBD has a much better uniformity than AC DBD for the fast rising and falling time.With the addition of O2,the uniformity of ns pulsed DBD gets worse for the space electric field distortion by which promotes the filamentary formation.While, in AC DBD, the added O2 can reduce the intensity of filaments,which enhances the discharge uniformity.The ns pulsed DBD has a much higher instantaneous power and energy efficiency than AC DBD.The ratio of Ar emission intensities indicates that the Te drops quickly with the addition of O2 both ns pulsed and AC DBDs and the ns pulsed DBD has an obvious higher Te and ne than AC DBD.The results are helpful for the realization of the reactive and uniform low temperature plasma sources.

    Keywords: nanosecond pulse power supply, dielectric barrier discharge, discharge uniformity,oxygen addition

    1.Introduction

    Dielectric barrier discharge(DBD)can produce large area low temperature plasma at a large variety of conditions (gas pressure, electrode structure, and power supply parameter)[1-3], which makes it widely used in air pollution control[4, 5], material surface treatment [6-8], plasma biomedicine[9] and energy conversion.One of the advantages of DBD is that the dielectric layer between the electrodes can prevent the increase of the discharge current into arc and keep the discharge under uniformity in atmospheric pressure [3, 10].But at atmospheric pressure, especially with the addition of admixture to improve plasma-chemical reactivity, the uniformity of DBDs is still not satisfied for many applications because usually there are a lot of short-live conductive microdischarge filaments existing in DBDs [11].With most of the energy and active species restricted into the filaments, there are low energy efficiency, uneven treatment, and local overheating problems in the DBDs in the filamentary discharge mode.

    Traditionally,DBDs are excited by AC power supplies,and recently, with the development of pulse power technology, the application of nanosecond (ns) pulse power supply in the excitation in DBDs becomes one of important methods to improve the discharge uniformity.It is because the short discharge duration can further limit the formation of filaments and prevent the streamer-to-spark transition.Besides the improvement of the discharge uniformity,the energy efficiency of the pulsed DBDs can be improved compared with AC DBDs because more energy goes into energetic electrons production other than gas heating [12-15].The high reduced electric field E/n caused by the high voltage slew rate and the high density chemical species produced by energetic electrons would enhance the plasmachemical efficiency.Therefore, the investigations of the discharge characteristics and applications of the pulsed DBDs attract much attention to improve discharge uniformity, energy efficient,plasma activity,and treatment effects[16].Walsh et al compared the efficiency of pulse and high frequency excited Ar DBDs and concluded that the electron density and atomic oxygen content in pulsed plasma were 3.9 and 3.3 times higher than those excited by high frequency, respectively [17].Jiang and Yang et al respectively studied and compared the optical images of AC and pulsed DBDs under different conditions[18,19].Wang et al investigated the discharge characteristics of AC and ns pulsed coaxial DBDs and found that the outside wall temperature of the AC coaxial DBD reached 158 °C and was much higher than that of the ns pulsed DBD at 64.3 °C after 900 s operation [20].Zhang et al studied the regime transitions between streamer discharge, diffuse discharge, and surface discharge in ns pulsed DBD using high resolution temporal-spatial spectra and instantaneous exposure images[21].Yu et al studied the ns pulsed DBD temporally resolved evolution by ICCD images [22].

    Meanwhile, the amounts and kinds of the active species in plasmas are also critical important in many applications because those species participate physicochemical processes[23].Usually, to improve the uniformity of atmospheric DBDs, the discharges need operate in inert gases, such as Ar and He,and there is lack of active and functional radicals for plasma-chemical processing.One of the solutions is the addition of additives (such as O2, H2O, CF4) to produce reactive species(such as O,OH,H,NO,HO2,O3,F)[24-26].Among of those additives, O2plays a crucial role because it can produce oxidative radicals, such as O and O3, which are important in thin film deposition and etching, VOCs oxidation, and biomedical treatment applications.Fang et al used Ar DBD for PET film surface hydrophilic modification and found that appropriate amount of oxygen addition(0.3%)can enhance the treatment effect [27].Weng et al inactivated bacteria with radio-frequency Ar DBD and found the optimum of inactivation of E.coli treatment effect at the addition of 0.5% oxygen and 30 s exposure time [28].Hicks et al studied the process optimization of Ar/O2/CF4plasma treatment of diamond substrate smoothing and the effect of removing subsurface polishing damage [29].

    However, the added oxygen would affect the plasma properties more complicated.A proper amount of oxygen admixture could increase the discharge by Penning effect with metastable molecules of inert gases.However,with excessive oxygen addition, it will reduce the electron density by formingand quench the metastable molecules, which exacerbate the distortion of the electric field distribution and the non-uniformity of discharge.For the purpose of realizing uniform and reactive oxygen containing DBD plasma sources, the understanding of the discharge phenomena is very important, which reveals the effect of O2in the plasma processing.

    Considering the widely use of AC and pulsed power supplies in the excitation of DBDs,in this work,the effects of oxygen addition on the characteristics of ns pulsed DBD and AC DBD are compared and quantitative analyzed.The measured voltage and current waveforms are separated to reveal the ignition and extinction of the discharge processes, and to estimate ne.The tendency of Teunder different O2contents is derived from the two emission lines by optical emission spectroscopy (OES) measurement.

    2.Experimental setup

    Figure 1 shows the experimental setup, which includes a discharge reactor, electrical and optical measurement systems, gas supply system, and AC and ns pulse power supplies.The stainless steel reactor can be vacuumed with a rotary pump.A pair of 50 mm diameter parallel plates is used as electrodes, and an 80 mm diameter 2 mm thick quartz glass is used as dielectric barrier.The gas gap between electrode plates is fixed at 2 mm.The voltage range of the AC power supply is 0-30 kV and frequency is from 1 Hz to 15 kHz.The ns pulse power supply can provide 0-15 kV peak pulse high voltage and 1 Hz to 100 kHz pulse repetition frequency.The rising and falling times,pulse width of the pulse voltage are in the ranges of 50 ns to 500 ns, and 0-5000 ns, respectively.In the experiment, the frequencies of both the power supplies are fixed at 10 kHz,and the rising and falling times are fixed at 50 ns, which can produce discharges on the rising and falling edges.The pulse width of the pulse voltage is fixed at 800 ns for a better discharge uniformity and stability.

    Figure 1.The schematic of the DBD experimental setup.

    The electrical measurement system is composed of a digital oscilloscope (Tektronix 3054), a high voltage probe(North Star PVM-5)and a current coil(Pearson Electronics Inc.2877) for the voltage-current measurement from the power supplies,and a differential probe(LDP-6002)for the voltage measurement on a 2.2 nF measuring capacitor (Cr)for electrical parameters calculation.A Canon EOS 6D camera is used to record the discharge image with 1/6 s exposure time.The OES is measured by an Ocean Optic HR4000CG spectrometer (0.7 nm resolution).The purity Ar and O2are used as discharge gases and a bottle of mixture gas of Ar/O2with 1% O2is used when the O2content is below 1% for the purpose of precise ratio control.

    3.Results and discussion

    3.1.Discharge uniformity of argon DBD

    Figure 2 shows the optical images of the AC and ns pulsed Ar DBDs, respectively.It can be seen that the AC DBDs are in filamentary mode with the applied voltage from 2.4 to 3.6 kV.At a lower voltage,2.4 kV,the discharge cannot be full of the space between the electrodes and just with isolated filaments.With 3.2-3.6 kV applied voltage, the discharge occupies all of space but is still in filamentary mode.Similarly, the Ar DBD excited by the ns pulsed power supply shows isolated filaments in the discharge space at a lower voltage, 3.3 kV.With the increasing of the applied voltage,the plasma fills in the whole discharge space without obvious filaments at 3.6 kV.The Ar DBD excited by ns pulsed power supply is more uniform, which is due to the short rise time of pulse voltage, thus terminating the transition from Townsend glow discharge to serious filamentous phase.However, as the applied voltage further increases to 4.4 kV, filaments appear in the discharge due to high voltage.

    Figure 2.The discharge images of DBD excited by AC (a) and ns pulse (b) power supplies.

    The experimental results show that the applied voltage range of ns pulse Ar DBD is larger than that of AC,which is due to the over-voltage breakdown of pulse discharge,and the gas needs higher voltage to be ionized.From figure 2, when the voltage of AC is 2.4 kV,several clear and countable weak filaments are generated locally in the discharge space.With the increase of the applied voltage to 3.6 kV,the electric field intensity increases, the collision ionization increases, and the electron avalanche intensifies, resulting in the discharge filament filling the whole space,the discharge intensity increases,and the discharge is still uneven.When the applied voltage of ns pulse is 3.3 kV,the discharge does not fill the whole space as AC DBD when the voltage is 2.4 kV.While, except the filaments,more diffuse plasma fills the gap in ns pulsed DBD than AC DBD.With the applied voltage increasing, the filament gradually fills the whole space and fuses in the ns pulsed DBD.When the applied voltage reaches 3.6 kV, it is hard to observe the filament in the ns pulsed DBD,and the discharge is with good uniformity.With 4.4 kV applied voltage of ns pulsed DBD, a few bright streamers are generated in the space, and the discharge uniformity of ns pulsed DBD decreases.

    Figure 3.The discharge gap voltage separation waveforms of AC DBD(a)and ns pulsed DBD(c),and the current separation waveforms of AC DBD (b) and ns pulsed DBD (d).

    Although the discharge images can be qualitatively judged by eye,it is impossible to quantitatively get the value of the discharge uniformity.To quantitatively analyze the discharge uniformity, some researchers have developed different image processing methods based on fast Fourier transform, spatial correlation function, and gray level histogram[30-32].In our previous research,we have developed a method of gray value standard deviation (GVSD) to quantitatively analyze the discharge uniformity from the optical images [33].The pixel gray values of discharge area can be obtained with graying processing, and the gray values xnof the middle horizontal line of the discharge area are averaged to get the averaged gray value xm.Then,the normalized valuecan be obtained:

    where v is the variance of pixel value.Then, the gradient value ofcan be obtained with differentiation processing,which is the value of GVSD and can be used to indicate the uniformity of discharge.

    With the above GVSD method for discharge uniformity,the GVSD values of the AC DBD images at 2.4, 3.2, and 3.6 kV are 8.25, 4.34, and 4.68, and the GVSD values of the ns pulsed DBD images at 3.3,3.6,and 4.4 kV are 5.48,0.73,and 3.45.The lower GVSD value indicates the better discharge uniformity.The values of GVSD and the discharge images are consistent.For the GVSD below 1, the discharge seems very uniform and it can be defined as a uniform discharge, as shown in the ns pulsed DBD at 3.6 kV.

    3.2.Electrical characteristics of argon DBD

    Because the DBD uniformity is strongly associated with the electrical characteristics of discharge, the electrical parameters of the ns pulsed DBD and AC DBD are calculated.For DBD structure, the applied voltage Utis dropped on the gas gap and the dielectric barrier,as Ugand Ud,and the measured current Itwith current probe measurement is the sum of the displacement current, Id, and conduct current, Ig.As the method described in our previous studies[10,33,34],Udcan be estimated with an insert capacitor in the ground side and Ugis obtained with the subtraction of Udfrom Ut.Idcan be obtained without discharge under same applied voltage by pumping the reactor to vacuum or replace the inert gas with air.Igcan be obtained with the subtraction of Idfrom It.The separations of the applied voltage of AC and ns pulsed DBDs are shown in figures 3(a) and (c), and the separations of the measured currents are shown in figures 3(b) and (d).

    From the above figures that the peak values of Ugof AC and ns pulsed DBDs are much smaller than those of the applied voltages for part of the voltages dropped on the dielectric barriers.The peak value of Ugof AC DBD is about 1.2 kV and the peak value of Ugof ns pulsed DBD is about 1.6 kV.In ns pulsed DBD, the duration time of the Ugof the rising side of the pulse is very small, about 100 ns, which prevents the electron avalanches over development.At the falling side of the pulse, the accumulated charge causes the reverse stroke.In AC DBD, the duration time of the Ugis very long, about 25 μs, which results in many current pulses in a very long time scale.The electron avalanches are under lasting electric field and the filaments are easily formed.The peak value of Igseparated from ns pulsed DBD is about 2.4 A, which is much higher than that at mA level under AC DBD.The amount of transferred charge(Q)has an important influence on the reaction rate,efficiency and application effect of plasma.In general, the Q can be calculated by integrating Igover a period.However, due to the bipolar nature of Igduring AC and ns pulse DBD discharge, the Q at the rising edge and the falling edge will offset each other.Therefore,the method of calculating Q in the whole discharge process is to add the absolute value of Q in the first half cycle and the second half cycle, as shown in formula (2):

    According to the above formula, the values of Q of AC and ns pulsed DBDs are 142.5 nC and 151.1 nC,respectively.The results show that there are more active particles produced in ns pulsed DBD than in AC DBD.

    Discharge power is an important parameter to describe the intensity of discharge, including the instantaneous total power Pt, instantaneous gas gap power Pgand instantaneous dielectric layer power Pd, which can be obtained by equations (3)-(5):

    The Pt,Pgand Pdof Ar DBDs excited by AC and ns pulse are shown in figure 4.From figure 4 that the peak value of Ptunder AC DBD is around 30 W,and that of ns pulsed DBD is about 4.2 kW.The Pgof DBD excited by ns pulse is also larger than the value of AC DBD.Using the instantaneous power data,the corresponding average total powerPt,average air gap powerPg,average dielectric layer powerPdand energy efficiency η can be calculated by equations (6)-(9):

    Figure 4.The power separation waveforms of DBD excited by AC (a) and ns pulse (b).

    3.3.Optical emission spectra of argon DBD

    In DBD, a large number of electrons and ions are produced when the applied voltage exceeds the breakdown voltage and with the electric field in space, the electrons are accelerated and gain energy to become high-energy electrons, which eventually lead to the excitation, ionization, and dissociation of the background molecules into active species.Some excited species can be observed by the optical emission spectra.In Ar DBD, most of the active species are excited states of argon atoms and the excited states of argon atoms transit back to the lower level excited or ground states to release photon, which can be observed

    The intensities of the observed optical emission spectra represent the relative number densities of the excited species.Figure 5 shows the emission spectra of AC Ar DBD and ns pulsed Ar DBD.The applied voltage of AC power supply is fixed at 3.2 kV, and that of ns pulse power supply is fixed at 3.6 kV.From figure 5, it can be seen that the excited species of Ar DBD excited by AC and ns pulse discharge are the same.But the intensities of the emission lines are different.The optical emission intensities of excited Ar in ns pulsed DBD are higher than those in AC DBD.Taking the Ar line(696.5 nm) as an example, under the same conditions, the optical emission intensity in ns pulsed DBD is about 2.31 times that in AC DBD.By comparing the emission intensities of the excited Ar,it can be concluded that ns pulsed Ar DBD can produce plasma with more active particle and higher reactivity than AC DBD.The main reason is that the applied voltage of ns pulse power supply can rise to a value much higher than the breakdown voltage of gas in a short time,so it can provide a larger initial electric field strength, generate a greater instantaneous power, and significantly improve the reaction rate, so as to produce more active particles in the discharge space.

    Figure 5.The optical emission spectra of AC and ns pulse excited Ar DBD.

    3.4.Discharge uniformity with different oxygen contents

    With the addition of oxygen, it participates in the discharge process and affects the discharge uniformity and reactivity.Figure 6 shows the discharge images of Ar DBDs under AC and ns pulse excitations with different oxygen contents.The applied voltage of AC power supply is fixed at 3.2 kV, and that of ns pulse power supply is fixed at 3.6 kV.To quantitatively analyze the variation of Ar DBD uniformity with oxygen content under different power supply conditions, the GVSD of the discharge images of DBDs excited by different power supplies at various oxygen additions are calculated.Figure 7 is the GVSD value of figure 6.The quantitative analysis of the optical discharge images shows that the regularity of uniformity change reflected by the GVSD value is consistent with the luminescent images.

    Figure 6.The discharge images of DBD excited by AC (a) and ns pulse (b) power supplies at various oxygen additions.

    Figure 7.The gray value standard deviation of DBD excited by AC and ns pulse power supplies at various oxygen additions.

    It can be seen that under high frequency AC excitation,when the oxygen content is in the range of 0%-0.2%, the discharge is more intense and the discharge gap is full of bright filaments,and the uniformity is poor.With the increase of oxygen content, the electronegative gas O2will adsorb a large number of electrons and inhibit the development of electron avalanche, resulting in the decrease of discharge intensity,the decrease of bright discharge filaments.When the oxygen content increases to 1%, the discharge is relatively stable, the discharge filaments become dense, and the uniformity is improved.In this stage,the GVSD decreased from 5.34 to 2.22.After that,with the increase of oxygen content to 6%,due to serious space electric field distortion,the discharge intensity further weakens, the number of discharge filaments decreases, the uniformity decreases again, and the GVSD gradually increased to 4.58.When the ns pulse is used as the excitation power supply, under the condition of extremely low oxygen content (0%-0.03%), the discharge uniformity under different oxygen content has little change and the GVSD of the discharge image is less than 1.This is due to the Penning ionization between Ar metastable particles and O2,which produces a large number of electrons and leads to the cross fusion of electron avalanches, so that the uniformity of discharge can be maintained.With the increase of oxygen content, the uniformity of DBD discharge decreases gradually.When the oxygen content reaches 0.1%,bright filaments appear in the discharge gap, the discharge begins to be unstable, and the uniformity decreases rapidly.In this stage,the GVSD increased to 2.86.With the increase of oxygen content to 0.2%-0.6%, the discharge filaments become dark,the number of filaments increases, the distribution of filaments in the discharge space tends to be uniform, and the GVSD value is maintained at about 2.When the oxygen content increases to 6%, a large number of metastable particles of Ar were quenched and the number of space free electrons was further weakened.Therefore, the discharge intensity becomes weaker, the brightness of filament discharge decreases, and the GVSD value decreases.When O2was added further(0.6%-6%),O2adsorbed a large number of electrons to formions, which led to the electric field in space further distortion and the poor uniformity.The excess O2cannot fill the whole discharge space.In this stage, the GVSD further increased to 4.58.

    3.5.Electrical properties with different oxygen contents

    Figure 8 shows the variation trend of the transfer charge with oxygen addition in DBD under AC and ns pulse excitations.Because O2is a kind of electronegative gas, it can adsorb electrons, inhibit the development of electron avalanche, and lead to the decrease of discharge filaments.It can be seen that with the increase of oxygen addition, the transfer charge (Q)presents a continuous downward trend.When the oxygen content increases from 0%to 6%,the Q decreases from 142.5 nC to 37.9 nC of DBD under AC excitation, and from 151.1 to 70.6 nC of DBD under ns pulse excitation.Under the same oxygen content,the value of Q of DBD excited by ns pulse is larger than that excited by AC.The reason is that the DBD discharge time of ns pulse excitation is much shorter than that of DBD AC excitation in a single cycle, and ns pulsed DBD has higher Ugand Ig.Figure 9(a)shows the variation trend of average discharge power of DBD excited by AC and ns pulse with oxygen addition.As shown in figure 9,with the increase of O2content,andof AC and ns pulsed DBDs decreased.Among them, when the oxygen content increases from 0%to 6%,theof AC decreases from 8.6 to 6.8 W,andof AC DBD decreases from 1.5 to 0.7 W.Theof ns pulsed DBD decreases from 2.9 to 0.9 W,andof ns pulsed DBD decreases from 1.9 to 0.5 W.

    Figure 8.The variations of transported charges of DBD excited by AC and the ns pulse power supplies at various oxygen additions.

    Figure 9.The variations of averaged power (a) and energy efficiency (b) of DBD excited by AC and ns pulse power supplies at various oxygen additions.

    Figure 9(b) shows the curves of energy efficiency with oxygen addition of DBDs under the excitation of two kinds of power sources.It can be found from figure 9(b) that the η of DBDs of AC and ns pulse excitation decreases with the addition of O2.Among them, η of ns pulsed DBD is in the range of 56.68%-67.01%, which is much larger than that of DBD of AC excitation (11.03%-17.17%).Generally, the energy efficiency is mainly related to the charge utilization rate and dielectric heat loss.Due to the short time interval between the two discharges of ns pulse, the charge will not dissipate immediately after the discharge at the rising edge of the voltage, but will still be on the surface of the dielectric barrier or in the gap, which is conducive to the generation of discharge at the falling edge of the pulse voltage.Therefore,the charge utilization rate of ns pulse DBD is higher.Because the duration of ns pulse voltage is very short, the degree of energy concentration on the surface of dielectric is reduced,so the energy loss is less.On the other hand, because O2is a kind of electronegative gas with high electron affinity, the higher the oxygen content is, the more electrons it adsorbs and the faster the quenching speed of active particles is.Therefore, the energy efficiency decreases continuously with the increase of oxygen content.

    3.6.Optical emission spectra with different oxygen contents

    Figure 10 shows the emission spectra of DBD excited by AC and ns pulse with oxygen content of 0.4%.After adding oxygen into Ar, the composition of the main spectral lines observed from Ar/O2DBD excited by two kinds of power sources has no obvious change.However, a small amount of O I(777.4 nm)can be found(as shown in the enlarged part in figure 10 with the spectra of AC DBD in pure Ar as a reference), and the intensity of O line of AC DBD is smaller than that of ns pulsed DBD.The reason is that the energy for the production of the upper state of O I(777.4 nm),O(3p5P),is higher than the energy of the observed Ar emission lines at 4p states(around 13 eV).The energies for the ground state O atom production and O(3p5P)excitation from ground state O are 5.1 eV and 3.6 eV [35], according to reactions (11) and(12),and 10.7 eV, respectively.Then, the O(3p5P)is mainly produced by the direct electron impact dissociation excitation,which needs 15.8 eV or 14.3 eV, and ns pulsed DBD has a higher electron temperature than the AC DBD.Therefore,the O I (777.4 nm) intensity is much weaker than the Ar emissions and the intensity of O line of ns pulsed DBD is larger than that of AC DBD.The energy of metastable Ar is 11.72 eV [10], which can dissociate O2molecular but cannot ionize the O2molecular (12.06 eV ionization energy [36]).But the upper states of Ar lines at 4p states with around 13 eV can collide and ionize O2molecular to produce more electrons

    Figure 10.The optical emission spectra of(a)AC and ns pulse excited Ar DBD at 0.4%O2,and(b)the enlarged spectra of AC Ar DBD,AC and ns pulse excited Ar DBD at 0.4% O2.

    Figure 11.The optical emission intensities of the Ar lines of DBDs excited by AC and ns pulse power supplies at various oxygen contents.

    Figure 11 shows the variation trend of the peak intensity of typical particles with oxygen concentration of DBDs under the excitations of two kinds of power sources.The typical optical line intensities under AC and ns pulsed DBDs decreased with the oxygen content increasing,and the typical line intensities of AC DBD are smaller than those of ns pulsed DBD.Taking the Ar I (696.5 nm) as an example, the intensity of the line decreased by 95.3% under AC excitation and 96.8% under ns pulse excitation after 6%oxygen added.This shows that a large number of electrons are adsorbed by O2to formons,which reduces the electron number density in space.As a result, the discharge intensity is weakened and the optical emission intensity of Ar line decreases significantly.

    The emission intensity I of Ar excited by electron collision can be given by formula (14) [10, 33]:

    where Cλ, A and τ are the spectral response of the optical measurement system,Einstein coefficient,and effective lifetime,respectively.(ε) is the excitation cross-section andf(ε) is the electron energy distribution function (EEDF).Considering the EEDF can be seen as Maxwell’s function for the electron local equilibrium, EEDF determines the value of Te.is proportional to the excitation rate k and equation(14)can be present as equation(15).[Ar]is the ground state Ar density.

    The change of experimental conditions results in the changes of the electron number density neand electron temperature Te.The ratio of two Ar emission intensities at different wavelengths λ with different excitation energies can eliminate the effect of neand is only determined by the ratio of k, as shown in equation (16).The ratio of k reflects the electron energy distribution of energetic electrons for the two excitation processes and the equation(16)can be rewritten as equation (17).With the ratio of k and the Bolzmann solver BOLSIG + code [37], the electron excitation energy Tecan be obtained [10]:

    From equation(17),it can be seen that for the determination of Te, the spectral response of different emission lines should be measured and the effective lifetime of the upper states should be obtained.The effective lifetime τ can be expressed by equation (18).

    where τ0is the radiative lifetime, [Q] is the quencher q number density, andis the collisional quenching rate coefficient.Considering the small amount of O2addition and the uncertainty of collisional quenching rate coefficient,the accurate value of the ratio of k is hard to be obtained,which results in the uncertainty of the value of Te.However,it is still possible to get the tendency of Tewith variable parameters for the purpose of comparation[33].

    Because the Ar I (696.5 nm, excitation energy: 13.3 eV)and the Ar I (912.3 nm, excitation energy: 12.9 eV) lines are separated from other lines, and the excitation energy difference is about 0.4 eV.Therefore,it can use the ratio of Ar I(696.5 nm)and Ar I (912.3 nm) intensities to express the relative concentration of high excited state particles and low excited state particles, which reflects the trend of Tewith different O2contents.Figure 12 shows the ratio of Ar I (696.5 nm) and Ar I(912.3 nm)emission intensities of DBDs under AC and ns pulse excitations at different O2contents.With the increase of O2content,the ratios of Ar I(696.5 nm)to Ar I(912.3 nm)spectral intensities of DBDs under the excitations of the two kinds of power supplies decease continually and the ratio of spectral intensities of AC DBD is lower than that of ns pulsed DBD.It means ns pulsed DBD has a higher Tethan AC DBD,and it is more obvious under low oxygen content.Assuming DBD in uniform status, the electron density, ne, can also be calculated with the drift velocity of electron, discharge area, and the peak values of Ugas described in our former study [33].The result shows that neis 2.26×1017m?3for ns pulsed Ar DBD and 1.15×1017m?3for AC Ar DBD and with 6%oxygen content,nedrops to 1.52×1017m?3for ns pulsed Ar DBD and 3.67×1016m?3for AC Ar DBD.

    Figure 12.The ratio of the emission intensities of Ar I (696.5 nm)and Ar I (912.3 nm) in AC and ns pulsed DBDs at various oxygen contents.

    4.Conclusion

    In this work, the uniformity and reactivity of AC and ns pulsed Ar/O2DBD are studied with optical and electrical measurements.With GVSD analysis, the discharge uniformity is obtained quantitatively and it is shown that the ns pulsed DBD has a much better uniformity than AC DBD.Compared with the AC DBD, the ns pulse DBD has higher peak value of Ug, Ig, Q, and total instantaneous powers.The maximum energy efficiency of ns pulse excitation is 65.5%,while that of AC excitation is only 17.4%.With the 0.03%O2addition,the GVSD of the ns pulsed DBD still keeps below 1,which indicates the uniform discharge.With more O2addition, the uniformity of ns pulsed DBD gets worse for the space electric field distortion by, which promotes the filamentary formation.For AC DBD, the intensity of filaments is weakened by the attachment effect of the added O2,which enhances the discharge uniformity.The Teis indicated by the ratio of two Ar emission intensities.The result shows that the ns pulsed DBD has an obvious higher Tethan AC DBD.With the addition of O2, the Teof both ns pulsed and AC DBDs drops quickly.With the O2addition from 0% to 6%, the neof ns pulsed Ar DBD decreases from 2.26×1017m?3to 1.52×1017m?3, and the neof AC Ar DBD decreases from 1.15×1017m?3to 3.67×1016m?3.

    Acknowledgments

    This work is supported by National Natural Science Foundation of China (Nos.52037004 and 51777091).

    猜你喜歡
    陳潔劉峰方志
    Ultrafast magneto-optical dynamics in nickel(111)single crystal studied by the integration of ultrafast reflectivity and polarimetry probes
    《湖中漫行》樂譜
    The investigation of OH radicals produced in a DC glow discharge by laser-induced fluorescence spectrometry
    History of World Smile Day
    陌上
    黑龍江民國方志所刊名家墨跡選
    書法賞評(2019年2期)2019-07-02 12:10:50
    嘉絨藏族地區(qū)的舊方志編纂
    西藏研究(2017年1期)2017-06-05 09:26:11
    劉峰:從“發(fā)明狂人”到“創(chuàng)業(yè)達(dá)人”
    金色年華(2016年14期)2016-02-28 01:44:29
    喇叭花兒搖銅鈴
    Average Incremenral Correlarion Analysis Model and Irs Applicarion in Faulr Diagnosis
    十八禁国产超污无遮挡网站| 精品无人区乱码1区二区| 狂野欧美白嫩少妇大欣赏| 老司机午夜福利在线观看视频| 精品国内亚洲2022精品成人| 在线观看66精品国产| 老熟妇仑乱视频hdxx| 丰满的人妻完整版| 夜夜爽天天搞| 最近2019中文字幕mv第一页| 男女那种视频在线观看| h日本视频在线播放| av天堂在线播放| www.色视频.com| 九九热线精品视视频播放| 午夜福利视频1000在线观看| 国内久久婷婷六月综合欲色啪| 国产成人一区二区在线| 亚洲精品在线观看二区| 淫秽高清视频在线观看| 亚洲人与动物交配视频| 亚洲欧美日韩高清在线视频| 国产老妇女一区| 亚洲美女黄片视频| АⅤ资源中文在线天堂| 亚洲四区av| 亚洲av不卡在线观看| 真人做人爱边吃奶动态| 成人av一区二区三区在线看| 国产伦精品一区二区三区视频9| 成年女人永久免费观看视频| 亚洲最大成人手机在线| 国产av一区在线观看免费| 成人高潮视频无遮挡免费网站| 此物有八面人人有两片| 欧美性猛交黑人性爽| 午夜免费激情av| 色av中文字幕| 成人高潮视频无遮挡免费网站| 日产精品乱码卡一卡2卡三| 亚洲av成人av| 日韩亚洲欧美综合| 国内精品久久久久精免费| 亚洲精品一卡2卡三卡4卡5卡| 日韩精品青青久久久久久| 天堂动漫精品| 亚洲欧美日韩高清专用| 麻豆国产av国片精品| 晚上一个人看的免费电影| 熟妇人妻久久中文字幕3abv| 日日啪夜夜撸| 少妇裸体淫交视频免费看高清| 国产精品美女特级片免费视频播放器| 日韩一区二区视频免费看| 熟女电影av网| 国产男人的电影天堂91| 男插女下体视频免费在线播放| 一级av片app| 天堂av国产一区二区熟女人妻| 国产老妇女一区| 少妇被粗大猛烈的视频| 一级毛片电影观看 | 免费av不卡在线播放| 美女大奶头视频| 成人漫画全彩无遮挡| 久久久国产成人免费| 中文在线观看免费www的网站| 欧美不卡视频在线免费观看| 婷婷六月久久综合丁香| 日韩高清综合在线| a级毛色黄片| 亚洲人成网站在线观看播放| 日韩精品青青久久久久久| 床上黄色一级片| 99久国产av精品国产电影| 国产在视频线在精品| 老熟妇乱子伦视频在线观看| 婷婷六月久久综合丁香| 国产亚洲精品久久久com| 久久久a久久爽久久v久久| 人人妻人人澡人人爽人人夜夜 | 亚洲国产精品成人久久小说 | 三级国产精品欧美在线观看| 国产成人福利小说| 精品人妻熟女av久视频| 亚洲av中文字字幕乱码综合| 亚洲精品一区av在线观看| 欧美性猛交黑人性爽| 在线免费观看不下载黄p国产| 国产激情偷乱视频一区二区| av中文乱码字幕在线| 99九九线精品视频在线观看视频| 夜夜夜夜夜久久久久| 亚洲人成网站在线观看播放| 一级a爱片免费观看的视频| 成人漫画全彩无遮挡| 国产69精品久久久久777片| 尾随美女入室| 日韩中字成人| 97热精品久久久久久| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩综合久久久久久| 香蕉av资源在线| 18禁裸乳无遮挡免费网站照片| 国产aⅴ精品一区二区三区波| 在线天堂最新版资源| 亚洲av免费在线观看| 好男人在线观看高清免费视频| 国产精品电影一区二区三区| 欧美又色又爽又黄视频| 嫩草影院新地址| 性欧美人与动物交配| 国产精品99久久久久久久久| avwww免费| 男女视频在线观看网站免费| 国产在线精品亚洲第一网站| 在线观看av片永久免费下载| 偷拍熟女少妇极品色| 性欧美人与动物交配| 国产精品国产高清国产av| 色综合亚洲欧美另类图片| 成人无遮挡网站| 国产 一区 欧美 日韩| 色哟哟哟哟哟哟| 亚洲最大成人av| 欧美3d第一页| 一个人看的www免费观看视频| av免费在线看不卡| 一进一出抽搐动态| 日韩,欧美,国产一区二区三区 | 国产一区二区三区av在线 | 亚洲av电影不卡..在线观看| 免费观看在线日韩| 嫩草影院精品99| 久久草成人影院| 伊人久久精品亚洲午夜| 国产精品免费一区二区三区在线| 欧美xxxx性猛交bbbb| 麻豆一二三区av精品| 成年女人看的毛片在线观看| 欧美性猛交黑人性爽| 69av精品久久久久久| 久久精品国产亚洲av天美| 久久久久久久久久成人| 两性午夜刺激爽爽歪歪视频在线观看| 午夜a级毛片| 久久九九热精品免费| 欧美高清性xxxxhd video| 国产精品亚洲美女久久久| av福利片在线观看| 内射极品少妇av片p| 午夜福利在线在线| 99九九线精品视频在线观看视频| 最新在线观看一区二区三区| 亚洲专区国产一区二区| 精品国产三级普通话版| 亚洲高清免费不卡视频| 女的被弄到高潮叫床怎么办| 亚洲精华国产精华液的使用体验 | 午夜a级毛片| 一夜夜www| 麻豆乱淫一区二区| 午夜日韩欧美国产| 国产免费男女视频| 色哟哟哟哟哟哟| 97超视频在线观看视频| 性色avwww在线观看| 亚洲精品乱码久久久v下载方式| 一级黄色大片毛片| 国内精品一区二区在线观看| 国产一区二区在线av高清观看| 色视频www国产| 国产精品国产高清国产av| 一级a爱片免费观看的视频| 在线观看午夜福利视频| 亚洲天堂国产精品一区在线| 嫩草影视91久久| 日韩欧美精品v在线| 91久久精品国产一区二区成人| 色视频www国产| 国产老妇女一区| 日本 av在线| 国产精品嫩草影院av在线观看| 久久久久精品国产欧美久久久| 欧美激情久久久久久爽电影| 国内揄拍国产精品人妻在线| 嫩草影院新地址| 亚洲人成网站高清观看| 国产视频一区二区在线看| 日韩精品有码人妻一区| 六月丁香七月| 最近最新中文字幕大全电影3| 大型黄色视频在线免费观看| 亚洲av一区综合| 麻豆乱淫一区二区| 国产免费男女视频| av国产免费在线观看| 性色avwww在线观看| 亚洲美女搞黄在线观看 | 蜜臀久久99精品久久宅男| 国产精品人妻久久久影院| 高清毛片免费观看视频网站| 午夜精品一区二区三区免费看| 欧美日本亚洲视频在线播放| 日日摸夜夜添夜夜爱| 极品教师在线视频| 国产黄片美女视频| 久久热精品热| 久久久久国内视频| 精品人妻一区二区三区麻豆 | 久久久午夜欧美精品| 亚洲真实伦在线观看| 丰满乱子伦码专区| 级片在线观看| 免费在线观看影片大全网站| 国产精品三级大全| 一进一出抽搐动态| 国产毛片a区久久久久| 最新中文字幕久久久久| 亚洲精品影视一区二区三区av| 亚洲欧美精品综合久久99| 国产成人freesex在线 | 成年av动漫网址| 亚洲欧美日韩东京热| 色在线成人网| 欧美激情在线99| 毛片一级片免费看久久久久| 级片在线观看| 不卡一级毛片| 亚洲国产欧洲综合997久久,| 国产av一区在线观看免费| 国产视频一区二区在线看| 欧美一区二区亚洲| 精品久久久久久久久久久久久| 精品熟女少妇av免费看| 久久久久久久久久黄片| 狂野欧美白嫩少妇大欣赏| 中文字幕人妻熟人妻熟丝袜美| 中文亚洲av片在线观看爽| 成人漫画全彩无遮挡| 少妇熟女欧美另类| 亚洲性久久影院| 人人妻人人澡欧美一区二区| 国产一区二区亚洲精品在线观看| 最近的中文字幕免费完整| 看黄色毛片网站| 亚洲精品日韩在线中文字幕 | 欧美在线一区亚洲| 欧美另类亚洲清纯唯美| 夜夜看夜夜爽夜夜摸| 日韩三级伦理在线观看| 国产精品爽爽va在线观看网站| 又爽又黄无遮挡网站| 天天一区二区日本电影三级| 久久精品国产自在天天线| 两个人视频免费观看高清| 国产单亲对白刺激| 中国美白少妇内射xxxbb| 亚洲成人精品中文字幕电影| 日本免费a在线| 精品无人区乱码1区二区| 在线播放无遮挡| 人人妻人人澡欧美一区二区| 免费av毛片视频| 国产亚洲欧美98| 成人永久免费在线观看视频| 一级a爱片免费观看的视频| 两个人视频免费观看高清| 亚洲久久久久久中文字幕| 美女被艹到高潮喷水动态| 男人舔奶头视频| 狂野欧美激情性xxxx在线观看| 欧美色欧美亚洲另类二区| 99久久精品一区二区三区| 五月伊人婷婷丁香| 免费搜索国产男女视频| 色尼玛亚洲综合影院| 啦啦啦观看免费观看视频高清| 亚洲人成网站在线播放欧美日韩| 日韩欧美三级三区| 欧美成人a在线观看| 欧美中文日本在线观看视频| 婷婷亚洲欧美| 一进一出好大好爽视频| 久99久视频精品免费| 亚洲在线自拍视频| 99热只有精品国产| 99在线视频只有这里精品首页| 黄色欧美视频在线观看| 两个人的视频大全免费| 舔av片在线| 99久国产av精品国产电影| 婷婷精品国产亚洲av| av免费在线看不卡| 给我免费播放毛片高清在线观看| 久久久久精品国产欧美久久久| 日本免费一区二区三区高清不卡| av在线天堂中文字幕| 久久午夜福利片| 日本在线视频免费播放| 亚洲一区二区三区色噜噜| 久久午夜亚洲精品久久| 成年版毛片免费区| 99热只有精品国产| 美女 人体艺术 gogo| 亚洲精华国产精华液的使用体验 | 久久久成人免费电影| 亚洲四区av| 成人精品一区二区免费| 国产三级在线视频| 久久午夜福利片| 91在线观看av| 九九久久精品国产亚洲av麻豆| 99久国产av精品国产电影| 成人一区二区视频在线观看| 91久久精品电影网| 免费人成视频x8x8入口观看| 中文字幕久久专区| 免费看a级黄色片| 人妻少妇偷人精品九色| 又黄又爽又免费观看的视频| 国产精品免费一区二区三区在线| 亚洲图色成人| 12—13女人毛片做爰片一| 欧美一区二区国产精品久久精品| 色av中文字幕| 成人国产麻豆网| 久久九九热精品免费| 人人妻人人看人人澡| 午夜a级毛片| 一区二区三区四区激情视频 | 久久久久久久午夜电影| 日韩大尺度精品在线看网址| 精品人妻视频免费看| 毛片一级片免费看久久久久| 免费av观看视频| 国产 一区精品| 久久久久国内视频| 日韩欧美在线乱码| 成人三级黄色视频| 淫秽高清视频在线观看| 91在线观看av| 精品午夜福利在线看| 国产日本99.免费观看| 欧美一级a爱片免费观看看| 国产成人91sexporn| 中文字幕精品亚洲无线码一区| 精品少妇黑人巨大在线播放 | 中国国产av一级| 国产91av在线免费观看| 免费观看精品视频网站| 国产男靠女视频免费网站| 免费观看人在逋| 国产成人一区二区在线| 久久这里只有精品中国| 在现免费观看毛片| 免费av观看视频| 亚洲精品色激情综合| 热99re8久久精品国产| 亚洲最大成人中文| 国产精品国产三级国产av玫瑰| 色尼玛亚洲综合影院| 亚洲精品亚洲一区二区| 97热精品久久久久久| 18+在线观看网站| 3wmmmm亚洲av在线观看| 国产v大片淫在线免费观看| 精品一区二区三区人妻视频| 欧美激情在线99| 淫妇啪啪啪对白视频| 国产成年人精品一区二区| 日本黄色片子视频| 国产成人精品久久久久久| 成人无遮挡网站| 国产亚洲精品av在线| 听说在线观看完整版免费高清| 国产乱人偷精品视频| 菩萨蛮人人尽说江南好唐韦庄 | 久久久久国产精品人妻aⅴ院| 国产一区亚洲一区在线观看| 淫妇啪啪啪对白视频| 无遮挡黄片免费观看| 在线观看午夜福利视频| 你懂的网址亚洲精品在线观看 | 97碰自拍视频| 性插视频无遮挡在线免费观看| 久久人妻av系列| 亚洲av美国av| 免费看av在线观看网站| 看非洲黑人一级黄片| 狂野欧美白嫩少妇大欣赏| 97碰自拍视频| 最近手机中文字幕大全| 日本黄大片高清| 精品免费久久久久久久清纯| 老女人水多毛片| av免费在线看不卡| 国产精品久久久久久久久免| 最新中文字幕久久久久| 日本黄色片子视频| 97人妻精品一区二区三区麻豆| 一本一本综合久久| 成人高潮视频无遮挡免费网站| 丝袜美腿在线中文| 18+在线观看网站| 欧美一区二区亚洲| 亚洲av中文字字幕乱码综合| 乱人视频在线观看| 日韩一区二区视频免费看| 久久久成人免费电影| 免费黄网站久久成人精品| 日产精品乱码卡一卡2卡三| 欧美3d第一页| 天天躁夜夜躁狠狠久久av| 亚洲精品影视一区二区三区av| 午夜视频国产福利| 成人漫画全彩无遮挡| 联通29元200g的流量卡| 国产蜜桃级精品一区二区三区| 免费观看的影片在线观看| 亚洲中文字幕一区二区三区有码在线看| 精品免费久久久久久久清纯| 麻豆一二三区av精品| 欧美日韩在线观看h| 成年女人毛片免费观看观看9| 悠悠久久av| 嫩草影院入口| 看免费成人av毛片| 亚洲精品国产成人久久av| 亚洲图色成人| 黄片wwwwww| 嫩草影院入口| 欧美激情久久久久久爽电影| 18禁在线无遮挡免费观看视频 | 变态另类丝袜制服| 女的被弄到高潮叫床怎么办| 国产精品久久久久久精品电影| 成年免费大片在线观看| 国产三级中文精品| 免费在线观看成人毛片| 悠悠久久av| 18禁在线无遮挡免费观看视频 | 亚洲五月天丁香| 国产伦精品一区二区三区四那| 久久精品综合一区二区三区| 国产亚洲精品综合一区在线观看| 十八禁网站免费在线| 午夜爱爱视频在线播放| 国产一区二区三区av在线 | 麻豆乱淫一区二区| 精品久久久久久久久亚洲| 精品午夜福利视频在线观看一区| 免费看a级黄色片| 成人亚洲精品av一区二区| 人人妻人人澡欧美一区二区| 99热这里只有精品一区| 午夜影院日韩av| 国产精品久久电影中文字幕| 亚洲国产精品成人综合色| 国产色婷婷99| 国产精品av视频在线免费观看| 天堂影院成人在线观看| 少妇丰满av| 又爽又黄a免费视频| 欧美bdsm另类| 国产又黄又爽又无遮挡在线| 久久亚洲精品不卡| 女的被弄到高潮叫床怎么办| 99久久精品热视频| 亚洲图色成人| 亚洲av美国av| 国产精品永久免费网站| 亚洲av美国av| 久久久久精品国产欧美久久久| 国产av不卡久久| 中文字幕人妻熟人妻熟丝袜美| 国产成人福利小说| 日日摸夜夜添夜夜添小说| 深爱激情五月婷婷| 亚洲国产日韩欧美精品在线观看| 在线免费观看的www视频| 欧美国产日韩亚洲一区| eeuss影院久久| 男女那种视频在线观看| 亚洲一级一片aⅴ在线观看| 午夜视频国产福利| 亚洲精品粉嫩美女一区| 乱码一卡2卡4卡精品| 亚洲精品粉嫩美女一区| 免费人成在线观看视频色| 色吧在线观看| 午夜福利成人在线免费观看| 五月伊人婷婷丁香| 中出人妻视频一区二区| 搡女人真爽免费视频火全软件 | 人人妻人人看人人澡| 欧美精品国产亚洲| 久久久精品大字幕| 午夜福利18| 天堂影院成人在线观看| www.色视频.com| 如何舔出高潮| 美女黄网站色视频| 最近最新中文字幕大全电影3| 菩萨蛮人人尽说江南好唐韦庄 | 午夜免费男女啪啪视频观看 | 波多野结衣巨乳人妻| 国产午夜精品久久久久久一区二区三区 | 日本a在线网址| 99久久九九国产精品国产免费| 欧美zozozo另类| 久久久久国产网址| 中国美女看黄片| 亚洲色图av天堂| 老司机影院成人| 亚洲一区二区三区色噜噜| 成人鲁丝片一二三区免费| 午夜福利在线在线| 九九爱精品视频在线观看| avwww免费| 国产欧美日韩精品一区二区| 国产国拍精品亚洲av在线观看| 欧美日韩国产亚洲二区| 中国国产av一级| 日韩欧美在线乱码| 一区二区三区四区激情视频 | 成人精品一区二区免费| 国产精品,欧美在线| 日日摸夜夜添夜夜添av毛片| 露出奶头的视频| 精品国内亚洲2022精品成人| 久久6这里有精品| 婷婷亚洲欧美| 亚洲欧美日韩卡通动漫| 麻豆久久精品国产亚洲av| 国产淫片久久久久久久久| 不卡一级毛片| 老司机福利观看| 亚洲精品日韩av片在线观看| 性欧美人与动物交配| 日韩三级伦理在线观看| 日韩大尺度精品在线看网址| 国产成人福利小说| 网址你懂的国产日韩在线| 69av精品久久久久久| 三级经典国产精品| 狠狠狠狠99中文字幕| 精品福利观看| 一边摸一边抽搐一进一小说| 亚洲欧美日韩东京热| 麻豆成人午夜福利视频| 日韩欧美 国产精品| 麻豆国产av国片精品| 天天一区二区日本电影三级| 校园人妻丝袜中文字幕| av在线老鸭窝| 国产精品无大码| 成人鲁丝片一二三区免费| 美女免费视频网站| 国产色婷婷99| 黄片wwwwww| 麻豆一二三区av精品| 我要看日韩黄色一级片| 日日摸夜夜添夜夜添小说| 亚洲av第一区精品v没综合| 国产毛片a区久久久久| 狂野欧美白嫩少妇大欣赏| 亚洲激情五月婷婷啪啪| 色视频www国产| 99热只有精品国产| 精品一区二区三区av网在线观看| 亚洲av美国av| 欧美激情久久久久久爽电影| 色播亚洲综合网| 看十八女毛片水多多多| 国产日本99.免费观看| 亚洲美女黄片视频| 欧美中文日本在线观看视频| АⅤ资源中文在线天堂| 国产欧美日韩精品一区二区| 久久久久久久久大av| 久久久欧美国产精品| 免费看美女性在线毛片视频| 麻豆久久精品国产亚洲av| 国产成人aa在线观看| 日日撸夜夜添| 国产精品亚洲美女久久久| 国产精品爽爽va在线观看网站| 亚洲中文字幕日韩| 桃色一区二区三区在线观看| 少妇人妻一区二区三区视频| 青春草视频在线免费观看| 国产一区二区激情短视频| 人人妻人人澡人人爽人人夜夜 | 女的被弄到高潮叫床怎么办| 国产一区二区亚洲精品在线观看| 女同久久另类99精品国产91| 欧美3d第一页| 九色成人免费人妻av| 日韩精品中文字幕看吧| 夜夜爽天天搞| 亚洲乱码一区二区免费版| 99久国产av精品| 亚洲av成人av| 变态另类成人亚洲欧美熟女| 中国美女看黄片| 精品一区二区三区人妻视频| 亚洲人成网站在线播| 黄色视频,在线免费观看| 身体一侧抽搐| 免费黄网站久久成人精品| 偷拍熟女少妇极品色| 免费av不卡在线播放| 成年女人毛片免费观看观看9| www.色视频.com|