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

    Spatial characteristics of nanosecond pulsed micro-discharges in atmospheric pressure He/H2O mixture by optical emission spectroscopy

    2022-02-24 08:58:52ChuanjieChen陳傳杰ZhongqingFang方忠慶XiaofangYang楊曉芳YongshengFan樊永勝FengZhou周鋒andRugangWang王如剛
    Chinese Physics B 2022年2期
    關(guān)鍵詞:永勝

    Chuanjie Chen(陳傳杰) Zhongqing Fang(方忠慶) Xiaofang Yang(楊曉芳)Yongsheng Fan(樊永勝) Feng Zhou(周鋒) and Rugang Wang(王如剛)

    1School of Information Engineering,Yancheng Institute of Technology,Yancheng 224051,China2School of Automotive Engineering,Yancheng Institute of Technology,Yancheng 224051,China

    Atmospheric pressure micro-discharges in helium gas with a mixture of 0.5%water vapor between two pin electrodes are generated with nanosecond overvoltage pulses. The temporal and spatial characteristics of the discharges are investigated by means of time-resolved imaging and optical emission spectroscopy with respect to the discharge morphology,gas temperature, electron density, and excited species. The evolution of micro-discharges is captured by intensified CCD camera and electrical properties.The gas temperature is diagnosed by a two-temperature fit to the ro-vibrational OH(A2Σ+–X2Π,0–0)emission band and is found to remain low at 425 K during the discharge pulses. The profile of electron density performed by the Stark broadening of Hα 656.1-nm and He I 667.8-nm lines is uniform across the discharge gap at the initial of discharge and reaches as high as 1023 m?3. The excited species of He,OH,and H show different spatio-temporal behaviors from each other by the measurement of their emission intensities,which are discussed qualitatively in regard of their plasma kinetics.

    Keywords: atmospheric pressure micro-discharges,nanosecond repetitively pulsed discharge,helium,optical emission spectroscopy

    1. Introduction

    Atmospheric pressure micro-discharges exhibit distinctive merits in applications, such as plasma medicine, surface treatment, combustion enhancement, and environmental remediation, for their high chemical reactivity, easy operation,and low maintenance cost.[1–3]With increasing pressure the electron–neutral collision frequency increases up to 109s?1at atmospheric pressure, leading to a sufficient momentum transfer via elastic collisions, and thus the discharges have a tendency to approach the state of thermal equilibrium.[4]Such high gas temperature restricts their applications in the treatment of heat-sensitive materials,such as biological tissue and polymers,[5]and leads to discharge contraction. Over the years, several techniques have been proposed to prevent excessive gas heating.[4]One of the most efficient methods is to generate the discharge by applying a pulsed high voltage with short duration across the electrodes, as the typical timescale of gas heating is of the order of hundreds of nanoseconds or longer. Therefore,nanosecond repetitively pulsed(NRP)discharges are of great interest for applications due to their nonequilibrium state.[6,7]Numerous designs of electrode geometry have been realized and one of which is the pin-to-pin discharge for it has the simplest electrode configuration and well reproducible in time and space.[8–10]

    In the last two decades,many properties of the pin-to-pin NRP discharges have been reported. Paiet al.investigated the three regimes (corona, glow, and spark) and streamer-tospark transition of NRP discharge in air.[9,11]Rusterholtzet al.verified the fast gas heating mechanism in air discharge occurs at nanosecond time scales due to the dissociative quenching of oxygen by electronically excited states of nitrogen N2(B/C).[12]Shaoet al.found that the fast runaway electrons and x-rays produced from the streamer head are benefit to obtain a diffuse and homogeneous discharge.[13]Huanget al.studied the effect of the repetition rate on the discharge properties, which found that the residual charges from previous discharge pulse play a key role in the breakdown period.[14]Recently a thermal spark discharge in air was generated by applying a single high voltage pulse.[15]Previous works on the discharges in air, nitrogen or noble (He, Ar) gases have been studied extensively with regard to the time evolution of the densities of electrons, reactive species (such as oxygen atoms,excited nitrogen molecules N2(A,B,C)and N2(υ),nitrogen oxide NOx, OH radicals), the neutral and vibrational temperatures.[16–19]However, the properties of such microdischarges show steep gradients in both time and space due to their filamentary in nature, which have significantly impact on the plasma chemistry. In many applications, water vapor (H2O) is usually added as a reactive gas or already present in the operation environment as impurity to providea rich source of OH radial which plays a vital role in plasma chemistry.[5,20,21]The specific case of NRP generated in water containing helium gas is an interesting one due to its low gas temperature and breakdown voltage,and high electron excitation energy in comparison to other gases.[19,21]For He–H2O mixture, Verreyckenet al. studied the kinetics of OH radical in NRP discharges using spatially resolved laser-induced fluorescence (LIF) and the plasma parameters (such as electron density and gas temperature) by means of OES.[19,22]In these studies, the optical emission is mostly captured from the center of discharge but there is no description of the spatial distribution of the aforementioned plasma parameters and excited species. In this work, the temporal evolution of discharge properties along the inter-electrode axis is investigated in He–0.5%H2O mixture pin-to-pin discharges by applying short overvoltage pulses repetitively. This will provides quantitative and qualitative insight into the plasma parameters and discharge kinetics for better understanding of process of such nanosecond pulsed micro-discharges.

    2. Experimental setup

    Figure 1(a) presents a schematic diagram of the experimental setup of the nanosecond pulsed discharge source. The geometry of the electrodes consists of two tungsten pin electrodes separated by a gap distance of 4.2 mm. The discharge was produced between the upper and bottom electrodes with diameters of 1.0 mm and 3.2 mm respectively, which are tapered to sharped tips with a radius of approximately 0.15 mm.A high-voltage (HV) power supply (NPG-18/3500, Megaimpulse: Modern Semi. Technol.) with a fixed pulse amplitude of 16 kV was connected to the upper pin electrode via a matched 75-Ω transmission line and the HV pulses were triggered by a function generator (Tektronix AFG3022B). A resistor of 72.5 Ω was placed between the bottom electrode and ground to avoid a short cut.[23]The voltage applied on the upper pin was measured by a high voltage probe (Tektronix P6015A). A Rogowski coil (Pearson electronics 2878) was employed to record the current through the discharges,which was placed between bottom pin and resistor. The waveforms of the applied voltage and current were recorded by a digital oscilloscope(RIGOL DS 1204B).The time delay between the voltage and current probes was obtained by recording the voltage and current curves of a vacuum capacitor. And the voltage and current probe signals are synchronized by shifting the current curve of?4.4 ns. In this work,the discharges were generated in a chamber to avoid air diffusion and the chamber was pumped down and filled three times with He–0.5%H2O gas. The gas flow rate was controlled by a mass flow controller (MFC) with a constant value of 1.0 L·min?1. The gas outlet was opened to the ambient air through an overpressure valve in order to maintain the pressure of the chamber to be 1 atm(1 atm=1.01325×105Pa).

    Fig. 1. Schematics of (a) experimental setup of pin-to-pin nanosecond pulsed discharge and (b) the optical emission spectroscopic measurement scheme.

    The axial plane of the pin-to-pin discharge was imaged by an achromatic lens(UV-NIR Triple Lens,f=90 mm)on the entrance slit of monochromator(Andor Shamrock 303i)with a magnification of 1:1. The entrance slit width was set to be 20μm. The plasma emission is transmitted to the monochromator equipped with three gratings of 150, 1800, and 2400 grooves/mm, respectively. Then, the time-resolved emission spectra were recorded by an intensified charged coupled device(ICCD, Andor iStar DH340T-25U-3). The external triggering signal of the ICCD was synchronized to that of the power supply by a series of TTL signals from the same function generator. The instrumental broadening profile was determined by using a low-pressure mercury–argon lamp(Newport Model 6035). The wavelength calibration of the spectrometer was also performed by the lamp with several well-known emission lines given by the NIST database.[24]

    3. Results and discussion

    3.1. Electric properties and time-resolved imaging

    Figure 2 displays the applied voltage and total current waveforms of the discharge. Note that the total current is almost entirely conduction current, as the current monitor was located at the ground electrode shown in Fig. 1. Besides,the displacement current was measured to be negligible by applying low voltage pulses without discharge. The voltage across pin electrodes driven by the high power supply starts to increase att=?5 ns with a rise rate of 2.6 kV/ns and reaches the peak of 16 kV.Within a few nanoseconds,the discharge is ignited by the overvoltage pulse as the breakdownfield in helium requires just 4.0 kV/cm.[25,26]Then, the discharge current starts to increase and the voltage drops rapidly.We observed that the current reaches its peak up to~100 A during the second half of the voltage, which is characteristic of spark discharges.[9]From Fig. 2, a weaker secondary discharge pulse is generated att=35 ns after the primary pulse.This is attributed to the fact that the primary pulse reflects from discharge to the generator and then towards back to the electrodes.[26]Since the length of the coaxial cable is about 3 m,the propagation time through the cable filled with teflon as insulator is estimated to be~32 ns,which is consistent with the time interval between two pulses. The energy coupled to the discharge by the primary and secondary pulses is estimated to be approximately 1.30 mJ and 0.18 mJ respectively,which is typical for NRP discharges in spark regime.

    Fig. 2. Applied voltage and current waveforms for the nanosecond pulsed discharge.

    The time-resolved images of the pin-to-pin discharges in He–0.5%H2O gases are captured by the ICCD camera with an exposure time of 5 ns as presented in Fig.3. The first observable light at the initial breakdown (t ≤5 ns) shows that the spark discharge fills the gap with a strong emission intensity.The dynamics evolution of the streamer propagation followed by spark discharge across the gap is almost indistinguishable with a gate time of 5 ns,since the propagation velocity of ionization front is higher up to 105m/s.[27]The optical emission of discharge is only detected during the period of current pulse(t ≤15 ns), presenting a homogeneous filament. And due to the high electric field near the tips,ions and electrons could be accelerated toward the cathode and anode respectively,which produces a large amount of secondary electrons accompanying with a high current of 100 A.Afterward(t ≤15 ns),the emission intensity of discharge shrinks into two small bright balls.And a secondary discharge is observed at 35 ns

    Fig. 3. (a) Time-resolved images of nanosecond pulsed discharge in He–0.5%H2O gas. The optical emission pattern is accumulated over 200 discharges,the gate width is 5 ns and the gain is kept constant. The time indicated in each picture corresponds to the delay time of ICCD.The emission intensities in all images are normalized for the same intensity scale. (b)Emission intensity along the inter-electrode axis for different delay time.

    Fig. 4. Calibrated overview spectra of the He+0.5% H2O nanosecond pulsed discharges in the range of 250 nm to 750 nm recorded with a gate of 5 ns at 10 ns and 55 ns,respectively.

    In Fig.4,overview emission spectra between 250 nm and 750 nm are recorded with a grating of 150 g/mm from the bulk of discharge at 10 ns and 55 ns, and the spectral-sensitivity calibration is performed using a standard light source (Energetiq EQ-99). During the discharge phase (t= 10 ns), the spectrum is covered by the neutral atomic helium and hydrogen lines, and OH(A–X) band. Besides, strong continuum emission from the bremsstrahlung of electron is detected dueto the high electron density.[23]In the recombination phase(t=55 ns),emission from helium atomic and molecular bands become weak, but the intensity of OH(A–X) and hydrogen lines(Hαand Hβ)are still strong.

    3.2. Gas temperature

    The gas temperatureTg,which represents the kinetic energy of heavy particles(e.g.,helium atoms and ions),is one of the most fundamental parameters. The measurements of gas temperature are usually performed by fitting the ro-vibrational emission spectra of diatomic species to obtain the rotational temperatureTr.[28]In water containing helium discharges,the emission of OH(A2Σ+–X2Π, 0–0) band with its bandhead at 309 nm is the most often used candidate for estimatingTr.[19]However, in many cases, the excited state OH(A) is not directly excited by the electron impact excitation with the ground state OH(X), so the nascent rotational distribution of OH(A) is not a mapping of its ground state.[28]As for high pressure, the rotational levels of the excited state could be thermalized by collisions with the neutral gas, provided the rotational–translational energy transfer time is shorter than the effective lifetime of the excited state.

    In this section,the rotational population of OH(A)state is represented by a superposition of two Boltzmann distributions as

    whereJ′is the rotation quantum number of the upper state,Ris the fraction of the population of the OH(A) with high rotational temperature,FJ′is the rotational energy of levelJ′,, andare the rotational temperatures of the high and low temperature components of OH(A) which represent the high and low rotational excited OH(A)characterized by twotemperature Boltzmann distribution. This behavior is related to the formation and relaxation processes of OH(A).Then,the emission intensity of OH(A–X,0–0)can be given by[28]

    hereAJ′J′′is the Einstein emission coefficient,his the Planck constant,λis the wavelength of the rotational transition, andΦinsis the line profile function. As shown in Fig. 5, a fit of the observed OH(A2Σ+–X2Π, 0–0) spectra with a synthetic spectra is performed, assuming a two-temperature rotational population distribution. The fitting is carried out using the nonlinear curve-fitting function FIT in the CURVE FITTING TOOLBOX of MATLAB. The values of,, andRare obtained simultaneously by using the method of least squares.

    In Fig.5,it is clear that the rotational distribution deviates significantly from equilibrium as the population of high temperature component of OH(A) accounts for 77% of the total population,which can be attributed to two aspects. On the one hand, the quenching rate constant for OH(A) by H2O is four orders of magnitude larger than that by helium atoms.[28]In spite of very low concentration of water added, the effective lifetime of OH(A)reduces significantly from several hundred ns to ten nanoseconds. On the other, the rate coefficient of the rotational energy transfer(RET)of OH(A)decreases dramatically with the increase of the rotational quantum numberJ′,which means that the high rotational levels cannot achieve the state of translation-rotational equilibrium.[29]With the increasing ofJ′from 0 to 10, the ratio of the effective lifetime and the RET time constant (τeff/RET) decreases from 12 to 2, as listed in Table 1. It is verified in the literature that the low part of rotational temperatureis a good estimate of gas temperature due to the faster RET rate of low rotational levels,i.e.,Tg=.[28]

    Fig. 5. Best fit of the experimental measured and synthetic spectra of the OH(A–X)transition with minimum residual of the fit.

    Table 1. The time constants of effective lifetime and RET of OH(A)molecule in helium+0.5% water discharge for different rotational quantum numbers. The effective lifetime and the RET are calculated by τeff =1/(++τRad)and RET=1/(+).[28]

    Table 1. The time constants of effective lifetime and RET of OH(A)molecule in helium+0.5% water discharge for different rotational quantum numbers. The effective lifetime and the RET are calculated by τeff =1/(++τRad)and RET=1/(+).[28]

    Rotational quantum τeff (ns) RET(ns) τeff/RET He+0.5%H2O J′≤5 10.4 0.9 ~12 J′≥10 15.5 7.8 ≤2

    The time-resolved gas temperature of discharge at midgap determined fromof OH(A)is illustrated in Fig.6. The gas temperature during the discharge pulse remains 425 K with a maximum error of 50 K.A simplified model for the gas heating in nanosecond pulsed discharges proposed by Paiet al.is employed to estimate the upper limit of gas temperature.[9]Assuming that the energy of discharge is converted totally into gas heating and this process is infinitely fast and adiabatic,the temperature of neutral gas can be given as[26]

    with the discharge energyεtotal=1.48 mJ(in Subsection 3.1),the density of heliumρo=0.1249 kg·m?3, the specific heatcapacity of helium gasCV= 5193.2 J·kg?1·K?1, the room temperature≈300 K.The volume of dischargeVis calculated by assuming a cylindrical volumeV=πR2D,whereRis obtained from the HWHM of the plasma emission in Fig.3.Taking these parameters into Eq. (3), the maximum gas temperature cannot exceed~720 K.From Fig.6,there is no significant gas heating during the discharge pulse which is in contrast to the spark discharge in air. Similar results have also been reported by Verreyckenet al.[19]and Chenet al.[26]in nanosecond pulsed helium discharges driven by a low voltage.

    Fig.6.Time evolution of the gas temperature in the center of the nanosecond pulsed discharges.

    3.3. Electron density

    The Stark broadening of atomic lines is the most often used technique for determining the electron number density in atmospheric pressure discharges.[30]As shown in Fig. 4,the intensity of Hβline is too weak to be detected during the discharge pulse, thus another Balmer line Hαatλ=656.3 nm and helium line atλ=667.8 nm are employed to extract the Stark broadening width. The Voigt shape of atomic line is a convolution of the Gaussian and Lorentzian components, where the Gaussian part originates from the instrumental and Doppler broadenings,[29]being given by ΔλG=,and the Lorentzian part is mainly ascribed to the van der Waals broadening,resonance broadening and Stark broadening,ΔλL=Δλvan+ΔλR+ΔλS. The full width at half maximum (FWHM) of Doppler broadening is a temperature dependent parameter withand. The resonance FWHM of He I line is calculated by ΔλR(He)=26.26/Tg,while that of Hαis negligible, ΔλR(Hα)≈0.[31]For the van der Waals broadening, one obtains the FWHM relation aswhereCtakes the values of 2.42 and 1.79 for Hαand He lines, respectively.[31]Substituting the gas temperature into the above expressions,one obtains the calculated FWHMs of these broadening mechanisms,as listed in Table 2. It is noted that even in the case of discharges driven by a high voltage pulse,the broadening of Hαdue to the Stark effect of external electric field(<80 kV·cm?1)is negligible.[16]

    Then, the Stark FWHMs of Hαand He I lines are obtained by deconvolving the experimental line profiles using a Voigt function. Two examples of the fitting of Hαand He I lines are illustrated in Fig. 7. The Stark broadening formula of Hαline is derived from the tabulated values of FWHMs in Ref.[32]for the electron temperature of 1 eV and reduced massμ?=10,given by

    where ΔλSandneare in units of nm and m?3,respectively.[16]From Fig.7(a)the electron density is determined to be 4.96×1022m?3by substituting the Stark FWHM of HαΔλS=0.522 nm into Eq.(4).In Fig.7(b),the Stark broadening of He I at 667.8 nm is also applied to checkneobtained by Hα. For an electron density of 1022m?3–1023m?3,the Stark broadening formula for He I is given by Ref.[30]andneis estimated to be 4.66×1022m?3. The results indicate that the electron density obtained from these two lines correspond well with each other.

    Table 2.Calculated full widths at half maximum(FWHMs)of different broadenings for Hα and He lines with Tg=425 K and p=1 atm.

    Fig.7. Line profile of the Voigt fit of(a)Hα and(b)He I lines at 667.8 nm measured at t=5 ns.

    The temporal evolution of electron density along the inter-electrode axis by means of the HαStark broadeningtechnique is depicted in Fig. 8. During the initial discharge phase (t ≤10 ns), the axial profiles ofneare uniform (~1×1023m?3) across the gap. Then, theneprofile becomes non-uniform characterized by two maxima located in the region near the electrodes. Furthermore,we can observe that the electron density near the cathode is higher than that near the anode,which could be attributed to the secondary electrons by ion bombardment. This is in agreement with the results for air spark discharges measured by the Mach–Zehnder laser interferometer when the discharge is fully formed.[33]

    Fig. 8. Temporal and spatial evolutions of the electron density during the discharge pulses in He–0.5%H2O nanosecond pulsed discharges. The anode is at the top of the image.

    3.4. Emission spectra of the excited species and discharge kinetics

    Figure 9 shows the time and space resolved emission intensity of the excited levels of helium. During the discharge onset, the emission of helium is found to be strong and uniform along the inter-electrode axis which is similar to the profile of electron density(see Fig.8),but a peak in He emission intensity is clearly observed att=5 ns despite the fact that the electron density do not reach its maximum value. This suggests that the excited helium atoms He?are predominantly populated by the direct electron impact excitation(e+He→e + He?) when the electron temperature is high due to the high electric field. Meanwhile,the conductive channel across the gap has built up accompanying with the current peak as shown in Fig.2. Afterward, the emission intensity of helium decreases rapidly and shrinks to the tips of electrodes especially the cathode where the electron density is locally higher.Another possible reason is that a strong light scattering reflected from the surface of the thicker cathode electrode. In addition, the optical emission of He in this case is only observed during the discharge pulse and early afterglow,which is quite different from that observed in helium or in helium with 0.05% H2O mixture.[34]In these case, we saw an increase in He emission again when the discharge is switched off with a fluorescence decay time of several hundred nanoseconds. It is obvious that these different behaviors is related to the addition of water molecules as will be discussed later.

    The spatio-temporal evolution of the emission intensities of OH(A–X) bandhead at 309 nm and Hαline at 656 nm is plotted in Fig.10. Both the OH(A–X)and Hαemission show peaks during the recombination phase,which is in contrast to the emission of He shown in Fig. 9. As the applied voltage declines, the electrons cool down resulting in an increase in OH(A–X)and Hαemission by electron-ion recombination.

    Fig.9. Temporal and spatial evolutions of the emission intensities of neutral helium lines from different transitions at 706.5 nm and 728.1 nm along the inter-electrode axis during the discharge pulse. The anode is at the top in both images.

    Fig.10. Time resolved axial profiles of the emission intensity of OH(A–X)at 309 nm and Hα at 656.1 nm along the inter-electrode axis for the He–0.5%H2O nanosecond pulsed discharges.

    However,there are some differences between them as follows: (i) the OH emission could be observed during the discharge (t <10 ns), while the Hαemission is almost unde-tectable;(ii)the maximum in OH emission is produced in the early afterglow,which precedes the maximum in Hαemission;(iii)the OH emission originates mainly from the center of discharge channel both during and after the discharge,but the Hαemission is from the region near the electrodes in the recombination phase,especially the cathode.

    The spatially homogeneous emission of OH(A–X) and Hαduring the discharge is similar to the profile of electron density, when the electric field is high. This indicates that OH(A) and H(n=3) are mainly produced by the direct dissociative electron excitation of H2O:[28]

    and direct excitation from their ground states:

    Since the excitation energy of H(n=3)(12 eV)is much larger than that of OH(A)(4 eV),the production of H(n=3)is more difficult than OH(A),and as a result the Hαemission is very weak during the discharge.

    In helium discharges without water vapor, the He emission could last for several hundred nanoseconds in the afterglow (not shown here). Besides, the helium molecular bandranging from 462 nm to 469 nm is detected in the spectra. These He?andin the afterglow are mainly produced by the dissociative recombination reactions of helium molecular ions,e.g.,e+→He?+He and, since He+is readily converted into He+2andvia three body collisions (see R1 in Table 3).[35]However,when 0.5%H2O is added to helium gas,in our case,both the He and He2emission disappeared in the afterglow. It is inferred that even if He+andions produced are comparable to the electron density they will be quenched rapidly by water molecules.

    Table 3. Main reactions for the loss of helium ions in He–0.5%H2O nanosecond pulsed discharges.

    In the case of He–0.5%H2O discharge,the dominant ion produced initially is still He+, rather than H2O+, for the ionization rate of He is higher than that of H2O.[34]The loss of He+ions, apart from electron-ion recombination reactions,could be the charge transfer reactions R4 and R5. The rates for these reactions can be estimated by, andrespectively, whereis the total rate coefficient of R5.[36]In our experimental conditions, the value ofν1is approximately equal toν2, which indicates that He+is quenched by H2O as fast as it is converted into. Furthermore, theions will also be quenched rapidly by the charge transfer reaction R6, which is an order of magnitude faster than R5,to produce H+,OH+,and H2O+ions.Consequently,the addition of H2O provides an important loss channel of He+ions that should be responsible for the disappearance of He?andemissions during the afterglow.

    As a result of the fast charge transfer reactions, a large amount of ionized water fragments(OH+,H+)and water ions(H2O+) will be produced during the discharge. The rates for R5 suggest that the production of water fragments is more likely than water ions. However, it is verified in the work of Verreyckenet al.that the charge transfer processes of H+and OH+ions with water molecules followed by dissociative recombination seems to be important in the production of OH in the afterglow.[34]When the applied voltage collapses, the electron temperature decreases drastically to the temperature of neutral gas with a decay time within~10 ns as the frequency of electron energy transfer in atmospheric helium is about 108s?1.[28]Hence,in the early afterglow,both the(dissociative) electron excitation (R1) and electron–ion recombination reactions (R7, R8) contribute to the production of OH(A),resulting in the maximum of OH(A–X)emission during the second half of the voltage pulse.

    There are several processes that can produce hydrogen atoms in the recombination phase, as list in Table 4. The reactions R7 and R8 cana prioribe excluded as the available energy are 7.5 eV at most, which is not enough to excite the upper level of H(n=3). The significant Hαemission can thus only be created by the three-body electron–atomic ion recombination R9. Since the rate coefficient of R9 is proportional to,the H(n=3)production depends strongly on the electron temperature and so the Hαemission peak lags behind the OH(A–X)emission peak.

    Table 4. Main reactions for the electron–ion recombination in the afterglow.

    The spatial profiles of OH(A–X) and Hαemission infer that the plasma kinetics are strongly dependent on the axial position,as shown in Fig.11. Similar mechanism is responsible for the dip in the radial profile of OH density in the midgap as discussed in the work of Verreyckenet al.[34,37]In the vicinity of electrodes where the electric field and electron temperature are higher, the water molecules could be dissociated more highly than in the mid-gap. Correspondingly, the discharge products in the region near the electrodes are dominantby water fragments (H and O atoms) and their ions, while in the mid-gap the water ions and its clusters(H2O+,H3O+)are more abundant. Consequently, the production of OH(A) in the mid-gap is more likely than the production of H(n=3).As reaction R9 shows a quadratic dependence on the electron density, significant Hαemission is detected near the cathode where the electron density is much higher as shown in Fig.8.

    Fig. 11. Schematic drawing of the plasma kinetics of H2O that are proposed to explain the observed spatial profile of OH(A–X)and Hα emission in Fig.10.

    4. Conclusion and perspectives

    The nanosecond pulsed pin-to-pin micro-discharge is produced in helium gas with water containing of 0.5% at atmospheric pressure. The temporal and spatial evolutions of the discharge characteristics are investigated by means of the optical imaging and optical emission spectroscopy. A uniform emission profile in space is observed at the initial phase of discharge.The population of OH(A2Σ+,υ=0)deviates from the Boltzmann distribution,which is attributed to the fact that the water molecules could reduce the effective lifetime of OH(A)significantly and the RET time constant decreases with the increasing of rotational quantum number. The gas temperature is extracted by a two-temperature fitting of the ro-vibrational OH(A–X)band and remains to be low at 425 K.The electron density measured by the Stark broadening of Hαand He is found to be uniform across the gap at the initial of discharge forne~1×1023m?3, and afterward shrinks to the tips of cathode and anode due to the secondary electron emission induced by the bombardment of helium ions and electrons. The evolution of He emission is similar to that of electron density,which indicates that the excited helium atoms are produced dominantly by the electron impact excitation. The production of OH(A)is attributed to the charge transfer between H2O and helium ions and subsequent the dissociative electron recombination of water ions,resulting in its maximum intensity in the early afterglow and the disappearance of the emission from helium and its molecules. The maximum in Hαemission during the recombination phase indicates that the exited H atoms are produced by electron–ion recombination of H+. According to the spatial profiles of OH(A–X) and Hαemission, the water molecules near the electrodes are more easily to be dissociated into ionized water fragments (H+and O+) than in the centre of discharge due to the higher electric field, which means that the discharge in the region near the electrodes is atomic. Since the production rate of excited H atoms is proportional to the square of electron density,the intensity of Hαis much stronger near the electrodes,especially the cathode.

    Acknowledgements

    Project supported by the National Natural Science Foundation of China (Grant No. 51806186), the Natural Science Foundation of the Jiangsu Higher Education Institutions of China (Grant No. 20KJB140025), the Natural Science Foundation of Jiangsu Province, China (Grant No. BK20181050),and the Scientific Research Project for the Introduction Talent of Yancheng Institute of Technology(Grant No.XJR2020).

    猜你喜歡
    永勝
    韓永勝
    大江南北(2022年11期)2022-11-08 12:04:18
    唱一首祖國的贊歌
    北京華語視聽牛永勝 以客戶為主,20年老店的“生意經(jīng)”
    Dielectric breakdown properties of Al-air mixtures
    一種兩級(jí)雙吸管道輸油泵
    山東科技大學(xué) 胡永勝國畫作品選登
    謝永勝
    寶藏(2018年6期)2018-07-10 02:26:38
    永勝油茶的深情厚意
    中國西部(2017年4期)2017-04-26 03:49:43
    風(fēng)情永勝 邊屯明郡
    中國西部(2017年7期)2017-01-24 12:57:16
    活著
    母亲3免费完整高清在线观看| 久久久久免费精品人妻一区二区 | 欧美日韩亚洲国产一区二区在线观看| 国产v大片淫在线免费观看| 国产单亲对白刺激| 久久久久九九精品影院| 欧美在线一区亚洲| 亚洲欧美一区二区三区黑人| 国产精品av久久久久免费| or卡值多少钱| 国产黄a三级三级三级人| 亚洲精品国产精品久久久不卡| 在线观看一区二区三区| 久热这里只有精品99| 国产不卡一卡二| 久久婷婷人人爽人人干人人爱| 精品久久久久久,| av免费在线观看网站| 亚洲午夜精品一区,二区,三区| 亚洲国产精品合色在线| av电影中文网址| 久久国产亚洲av麻豆专区| 怎么达到女性高潮| 国产精品 欧美亚洲| 国产片内射在线| 亚洲精品国产一区二区精华液| 美女扒开内裤让男人捅视频| 黄色丝袜av网址大全| 看免费av毛片| 久久久久国产精品人妻aⅴ院| 美女国产高潮福利片在线看| 精品久久久久久久人妻蜜臀av| 日日摸夜夜添夜夜添小说| 亚洲天堂国产精品一区在线| 淫妇啪啪啪对白视频| 国产1区2区3区精品| 亚洲av中文字字幕乱码综合 | 伊人久久大香线蕉亚洲五| 人人妻,人人澡人人爽秒播| 亚洲 国产 在线| 在线视频色国产色| 久久国产精品男人的天堂亚洲| 国内揄拍国产精品人妻在线 | 精品久久久久久久末码| 国产精品久久视频播放| 18禁美女被吸乳视频| 免费女性裸体啪啪无遮挡网站| 亚洲国产欧美一区二区综合| 中国美女看黄片| 大香蕉久久成人网| 在线观看免费视频日本深夜| 亚洲午夜精品一区,二区,三区| 亚洲五月色婷婷综合| 亚洲人成网站在线播放欧美日韩| 妹子高潮喷水视频| 欧美又色又爽又黄视频| av欧美777| 亚洲精品粉嫩美女一区| 日本免费a在线| 亚洲五月天丁香| 两性夫妻黄色片| 非洲黑人性xxxx精品又粗又长| 老司机午夜福利在线观看视频| 成人18禁在线播放| 国产又黄又爽又无遮挡在线| 欧美大码av| 九色国产91popny在线| 黄色a级毛片大全视频| 欧美乱码精品一区二区三区| 亚洲一区二区三区不卡视频| 又紧又爽又黄一区二区| 悠悠久久av| 午夜影院日韩av| 一进一出好大好爽视频| 久久香蕉国产精品| 亚洲人成77777在线视频| 淫秽高清视频在线观看| 精品免费久久久久久久清纯| 中国美女看黄片| 日本一区二区免费在线视频| 狂野欧美激情性xxxx| 亚洲精品美女久久久久99蜜臀| 精品电影一区二区在线| 国产精品久久久av美女十八| 久久精品国产清高在天天线| 亚洲熟妇熟女久久| 最近最新中文字幕大全免费视频| 亚洲五月天丁香| 久久久精品国产亚洲av高清涩受| 精品欧美一区二区三区在线| 观看免费一级毛片| 亚洲av片天天在线观看| 不卡av一区二区三区| 女同久久另类99精品国产91| 国产精品综合久久久久久久免费| 亚洲欧美日韩无卡精品| 美女高潮喷水抽搐中文字幕| 亚洲精品国产精品久久久不卡| 亚洲欧美一区二区三区黑人| 亚洲国产欧美一区二区综合| 亚洲人成电影免费在线| 99热这里只有精品一区 | 国产精品国产高清国产av| 夜夜看夜夜爽夜夜摸| 男女午夜视频在线观看| 在线永久观看黄色视频| 男女床上黄色一级片免费看| 欧美日韩亚洲综合一区二区三区_| 老鸭窝网址在线观看| 成年人黄色毛片网站| 成年免费大片在线观看| 亚洲天堂国产精品一区在线| 99国产精品99久久久久| 国产精品国产高清国产av| ponron亚洲| 午夜福利18| 久久天堂一区二区三区四区| 国产精品爽爽va在线观看网站 | 午夜成年电影在线免费观看| 精品国产乱子伦一区二区三区| 欧美日韩瑟瑟在线播放| 国产99白浆流出| 久久精品夜夜夜夜夜久久蜜豆 | 国产高清激情床上av| 久久久久久国产a免费观看| 夜夜爽天天搞| 99热6这里只有精品| 久久热在线av| 欧美绝顶高潮抽搐喷水| 在线观看66精品国产| 神马国产精品三级电影在线观看 | 99在线人妻在线中文字幕| a级毛片在线看网站| 久久久国产成人免费| 一卡2卡三卡四卡精品乱码亚洲| 啦啦啦韩国在线观看视频| 欧美日韩精品网址| 国产亚洲av高清不卡| 午夜影院日韩av| 深夜精品福利| 久久中文字幕人妻熟女| 桃红色精品国产亚洲av| 老司机午夜福利在线观看视频| 精品国内亚洲2022精品成人| 无限看片的www在线观看| 中文亚洲av片在线观看爽| 午夜福利18| 中文字幕人妻丝袜一区二区| 免费搜索国产男女视频| 日韩一卡2卡3卡4卡2021年| 正在播放国产对白刺激| 国产亚洲精品av在线| 国产一级毛片七仙女欲春2 | 97超级碰碰碰精品色视频在线观看| 制服人妻中文乱码| 亚洲午夜理论影院| 在线十欧美十亚洲十日本专区| svipshipincom国产片| 日日爽夜夜爽网站| 久久久久国产精品人妻aⅴ院| 一进一出抽搐gif免费好疼| 美女午夜性视频免费| 亚洲精品美女久久久久99蜜臀| cao死你这个sao货| 日本一区二区免费在线视频| 级片在线观看| 亚洲va日本ⅴa欧美va伊人久久| 一边摸一边做爽爽视频免费| 亚洲成人精品中文字幕电影| 亚洲熟妇中文字幕五十中出| 色播亚洲综合网| 三级毛片av免费| 日韩大尺度精品在线看网址| 欧美绝顶高潮抽搐喷水| 国产一区二区三区视频了| 日韩精品免费视频一区二区三区| 国产精品精品国产色婷婷| 99精品久久久久人妻精品| 欧美另类亚洲清纯唯美| 看黄色毛片网站| 侵犯人妻中文字幕一二三四区| 国产成人系列免费观看| 啦啦啦 在线观看视频| 中国美女看黄片| 99久久综合精品五月天人人| 热re99久久国产66热| 午夜福利视频1000在线观看| 国产区一区二久久| 久久久久久久久免费视频了| 国产黄片美女视频| 自线自在国产av| 免费av毛片视频| 久久精品国产亚洲av高清一级| 亚洲成人久久爱视频| 法律面前人人平等表现在哪些方面| 久久久久免费精品人妻一区二区 | 亚洲成人国产一区在线观看| 欧美黑人巨大hd| 桃色一区二区三区在线观看| 可以在线观看毛片的网站| 在线观看免费视频日本深夜| 日本五十路高清| 国产精品九九99| 精品熟女少妇八av免费久了| 两个人看的免费小视频| 国产亚洲精品久久久久久毛片| 亚洲国产看品久久| 成在线人永久免费视频| 欧美人与性动交α欧美精品济南到| 亚洲欧洲精品一区二区精品久久久| 亚洲国产精品合色在线| 亚洲精品中文字幕在线视频| avwww免费| 日日干狠狠操夜夜爽| 99久久精品国产亚洲精品| xxx96com| 日韩精品中文字幕看吧| www日本在线高清视频| 亚洲欧美精品综合一区二区三区| 国产成人系列免费观看| 一二三四社区在线视频社区8| 日韩 欧美 亚洲 中文字幕| 深夜精品福利| 国产欧美日韩精品亚洲av| cao死你这个sao货| 免费看日本二区| 亚洲av成人一区二区三| 悠悠久久av| 亚洲精品国产精品久久久不卡| 成年人黄色毛片网站| 国产高清videossex| 日韩国内少妇激情av| 久久欧美精品欧美久久欧美| 亚洲成av片中文字幕在线观看| 在线观看午夜福利视频| 国产亚洲欧美精品永久| 亚洲一区二区三区不卡视频| 午夜视频精品福利| 国产精品久久久久久精品电影 | 亚洲 欧美一区二区三区| 国产精品1区2区在线观看.| 欧美在线黄色| 免费看美女性在线毛片视频| 无人区码免费观看不卡| 一本精品99久久精品77| 国语自产精品视频在线第100页| 亚洲成人精品中文字幕电影| 嫩草影院精品99| 三级毛片av免费| 国产成人精品久久二区二区免费| 在线观看免费视频日本深夜| 久久99热这里只有精品18| 免费观看人在逋| 久久久久久亚洲精品国产蜜桃av| 亚洲午夜理论影院| 精品电影一区二区在线| 在线天堂中文资源库| 成年人黄色毛片网站| 一卡2卡三卡四卡精品乱码亚洲| 亚洲va日本ⅴa欧美va伊人久久| 国产一区在线观看成人免费| 中文字幕另类日韩欧美亚洲嫩草| 在线av久久热| 国产成人精品久久二区二区免费| 91字幕亚洲| 亚洲精品粉嫩美女一区| 一本一本综合久久| 黄色a级毛片大全视频| 国产精品免费视频内射| 久久天堂一区二区三区四区| 又大又爽又粗| 中文在线观看免费www的网站 | 色婷婷久久久亚洲欧美| 天天躁狠狠躁夜夜躁狠狠躁| 日韩视频一区二区在线观看| 国产国语露脸激情在线看| 一级毛片女人18水好多| 999久久久国产精品视频| 久久久国产精品麻豆| 俄罗斯特黄特色一大片| 嫁个100分男人电影在线观看| 欧美国产日韩亚洲一区| 变态另类成人亚洲欧美熟女| 男女做爰动态图高潮gif福利片| 亚洲国产精品成人综合色| 午夜免费激情av| 90打野战视频偷拍视频| 97碰自拍视频| netflix在线观看网站| 嫩草影院精品99| 无遮挡黄片免费观看| 精品卡一卡二卡四卡免费| 亚洲精品在线美女| 久久久久久久久久黄片| tocl精华| 成人欧美大片| 亚洲精品在线观看二区| 一个人观看的视频www高清免费观看 | 午夜亚洲福利在线播放| 精品国产乱码久久久久久男人| e午夜精品久久久久久久| 在线十欧美十亚洲十日本专区| 久久99热这里只有精品18| 日韩高清综合在线| 最新美女视频免费是黄的| 每晚都被弄得嗷嗷叫到高潮| 国产一区二区三区在线臀色熟女| 久久久国产成人精品二区| 亚洲精品一卡2卡三卡4卡5卡| 国产三级黄色录像| e午夜精品久久久久久久| 国产精品一区二区精品视频观看| 午夜福利一区二区在线看| 两性夫妻黄色片| 一本综合久久免费| 婷婷六月久久综合丁香| 十八禁人妻一区二区| 精品国内亚洲2022精品成人| 欧美日本亚洲视频在线播放| 亚洲精品中文字幕一二三四区| 午夜成年电影在线免费观看| 动漫黄色视频在线观看| av有码第一页| 免费看日本二区| 亚洲精品国产区一区二| 少妇裸体淫交视频免费看高清 | av片东京热男人的天堂| 又黄又粗又硬又大视频| 欧美在线黄色| 满18在线观看网站| 最新美女视频免费是黄的| 午夜福利高清视频| 男女之事视频高清在线观看| √禁漫天堂资源中文www| 久久香蕉精品热| 亚洲av日韩精品久久久久久密| 久久草成人影院| 在线观看免费视频日本深夜| 亚洲av成人av| xxxwww97欧美| 亚洲午夜精品一区,二区,三区| 免费在线观看完整版高清| 亚洲成人久久爱视频| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久免费精品人妻一区二区 | 一夜夜www| 亚洲成国产人片在线观看| 在线观看66精品国产| 国产一区二区在线av高清观看| 久久久久久亚洲精品国产蜜桃av| 亚洲色图av天堂| 免费在线观看完整版高清| 18禁黄网站禁片午夜丰满| 久久久久国产精品人妻aⅴ院| 日本精品一区二区三区蜜桃| 免费在线观看亚洲国产| 成人永久免费在线观看视频| 亚洲片人在线观看| 两性夫妻黄色片| 搡老熟女国产l中国老女人| 日本 av在线| 后天国语完整版免费观看| 欧美绝顶高潮抽搐喷水| 亚洲激情在线av| 中文字幕精品亚洲无线码一区 | 两个人免费观看高清视频| 久9热在线精品视频| 成人永久免费在线观看视频| 精品乱码久久久久久99久播| 亚洲成人久久爱视频| 在线十欧美十亚洲十日本专区| 亚洲成人国产一区在线观看| 非洲黑人性xxxx精品又粗又长| 成年版毛片免费区| 亚洲精品粉嫩美女一区| 欧美另类亚洲清纯唯美| 国产高清视频在线播放一区| 国产成人精品久久二区二区91| 看免费av毛片| 国产又色又爽无遮挡免费看| 午夜精品在线福利| 国产熟女xx| 少妇被粗大的猛进出69影院| 麻豆成人午夜福利视频| 不卡av一区二区三区| 国产私拍福利视频在线观看| 嫩草影院精品99| 国产色视频综合| 国产亚洲精品一区二区www| 国产精品久久久人人做人人爽| 欧美乱妇无乱码| 亚洲最大成人中文| 国产成年人精品一区二区| 欧美中文日本在线观看视频| 欧美三级亚洲精品| 免费高清视频大片| 欧洲精品卡2卡3卡4卡5卡区| 一边摸一边做爽爽视频免费| 国产成人精品无人区| 中文字幕高清在线视频| 村上凉子中文字幕在线| 搡老妇女老女人老熟妇| 亚洲男人的天堂狠狠| 老司机靠b影院| 午夜a级毛片| x7x7x7水蜜桃| 不卡一级毛片| 性色av乱码一区二区三区2| 九色国产91popny在线| 无遮挡黄片免费观看| 欧美成狂野欧美在线观看| 老熟妇乱子伦视频在线观看| 欧美性猛交黑人性爽| 精品第一国产精品| 2021天堂中文幕一二区在线观 | netflix在线观看网站| 免费电影在线观看免费观看| 亚洲一区高清亚洲精品| 国产午夜福利久久久久久| 亚洲av熟女| 欧美乱码精品一区二区三区| 美女高潮喷水抽搐中文字幕| 欧美丝袜亚洲另类 | 精品国产美女av久久久久小说| 男人的好看免费观看在线视频 | 国产精品免费视频内射| 亚洲av电影在线进入| 久久热在线av| 国产真人三级小视频在线观看| 99久久国产精品久久久| 亚洲三区欧美一区| 黄频高清免费视频| 精品久久久久久久人妻蜜臀av| 天堂影院成人在线观看| 黑人操中国人逼视频| 香蕉av资源在线| 中文字幕精品亚洲无线码一区 | 草草在线视频免费看| 人人妻人人澡人人看| 亚洲中文字幕日韩| 国产精品精品国产色婷婷| 两个人看的免费小视频| www日本黄色视频网| 精品熟女少妇八av免费久了| 精品不卡国产一区二区三区| 香蕉丝袜av| 国产激情偷乱视频一区二区| 日韩 欧美 亚洲 中文字幕| 国产极品粉嫩免费观看在线| 在线观看免费视频日本深夜| 欧美日韩亚洲综合一区二区三区_| 精品国产乱子伦一区二区三区| 在线永久观看黄色视频| 18禁黄网站禁片午夜丰满| 777久久人妻少妇嫩草av网站| 俄罗斯特黄特色一大片| 99精品在免费线老司机午夜| 婷婷六月久久综合丁香| 久久久久久亚洲精品国产蜜桃av| 日韩欧美 国产精品| 免费看十八禁软件| 在线十欧美十亚洲十日本专区| 久久久久久久午夜电影| 亚洲天堂国产精品一区在线| 美女大奶头视频| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品亚洲精品国产色婷小说| 深夜精品福利| 亚洲avbb在线观看| 91字幕亚洲| 大香蕉久久成人网| 老司机深夜福利视频在线观看| 欧美激情极品国产一区二区三区| 久久久久久久久免费视频了| 91麻豆精品激情在线观看国产| 亚洲男人的天堂狠狠| 中文字幕人成人乱码亚洲影| 十分钟在线观看高清视频www| 国产精品久久久久久亚洲av鲁大| 欧美日韩一级在线毛片| 亚洲精品中文字幕在线视频| 一边摸一边抽搐一进一小说| 成年版毛片免费区| 少妇裸体淫交视频免费看高清 | 国产区一区二久久| 日韩精品青青久久久久久| 男人操女人黄网站| 精品少妇一区二区三区视频日本电影| 欧美日韩亚洲综合一区二区三区_| 久久久久久久久中文| 色综合站精品国产| 精品国内亚洲2022精品成人| 国产成年人精品一区二区| 成人18禁高潮啪啪吃奶动态图| 久久精品国产亚洲av香蕉五月| 亚洲第一欧美日韩一区二区三区| 女人高潮潮喷娇喘18禁视频| 别揉我奶头~嗯~啊~动态视频| 精品国产超薄肉色丝袜足j| www.www免费av| 国产精品98久久久久久宅男小说| 亚洲欧美一区二区三区黑人| 一本综合久久免费| 淫秽高清视频在线观看| 热re99久久国产66热| 女性生殖器流出的白浆| 可以在线观看毛片的网站| 国产欧美日韩精品亚洲av| 亚洲人成网站高清观看| 在线免费观看的www视频| 国产人伦9x9x在线观看| 一a级毛片在线观看| 国产成人欧美在线观看| 亚洲国产毛片av蜜桃av| e午夜精品久久久久久久| 神马国产精品三级电影在线观看 | 日本撒尿小便嘘嘘汇集6| www.自偷自拍.com| 非洲黑人性xxxx精品又粗又长| 亚洲av成人不卡在线观看播放网| 免费在线观看黄色视频的| 免费在线观看影片大全网站| 国产熟女午夜一区二区三区| 少妇粗大呻吟视频| 久9热在线精品视频| 国产精品综合久久久久久久免费| 真人一进一出gif抽搐免费| 身体一侧抽搐| 一本一本综合久久| 老汉色∧v一级毛片| 嫁个100分男人电影在线观看| 亚洲国产高清在线一区二区三 | 精品久久久久久久久久久久久 | 久久午夜亚洲精品久久| 亚洲中文字幕一区二区三区有码在线看 | 国内精品久久久久精免费| 国产av不卡久久| 中文字幕av电影在线播放| 成人亚洲精品一区在线观看| 波多野结衣巨乳人妻| 免费高清视频大片| 波多野结衣巨乳人妻| 极品教师在线免费播放| 黑人操中国人逼视频| 嫩草影院精品99| 一本大道久久a久久精品| 亚洲av第一区精品v没综合| 欧美绝顶高潮抽搐喷水| 国产高清有码在线观看视频 | 村上凉子中文字幕在线| 成人国语在线视频| 久久国产乱子伦精品免费另类| 亚洲第一欧美日韩一区二区三区| 18禁黄网站禁片免费观看直播| 人妻久久中文字幕网| 男女做爰动态图高潮gif福利片| 在线播放国产精品三级| 亚洲第一青青草原| 黄色a级毛片大全视频| 久久香蕉国产精品| 人人妻人人看人人澡| 国产精品精品国产色婷婷| 亚洲专区字幕在线| 久久久久久久精品吃奶| 亚洲欧美精品综合一区二区三区| 一级片免费观看大全| 精品少妇一区二区三区视频日本电影| 大香蕉久久成人网| 久久久久久人人人人人| 黄色片一级片一级黄色片| 婷婷精品国产亚洲av在线| 首页视频小说图片口味搜索| 国产又爽黄色视频| 法律面前人人平等表现在哪些方面| 香蕉av资源在线| 人妻丰满熟妇av一区二区三区| 美女高潮到喷水免费观看| 一区二区三区国产精品乱码| 91成年电影在线观看| 亚洲自偷自拍图片 自拍| 香蕉丝袜av| 757午夜福利合集在线观看| av超薄肉色丝袜交足视频| 国产免费男女视频| 51午夜福利影视在线观看| 亚洲精品一卡2卡三卡4卡5卡| 最近最新免费中文字幕在线| 国产精品亚洲一级av第二区| 69av精品久久久久久| 日本成人三级电影网站| 久久午夜亚洲精品久久| ponron亚洲| 在线av久久热| 午夜a级毛片| 最新美女视频免费是黄的| 免费高清视频大片| 9191精品国产免费久久| 国产黄色小视频在线观看| 久热爱精品视频在线9| 给我免费播放毛片高清在线观看| 欧美一区二区精品小视频在线| 很黄的视频免费| 国产精品香港三级国产av潘金莲| 色综合欧美亚洲国产小说| 1024香蕉在线观看| 国产99久久九九免费精品| 久久久精品欧美日韩精品| 亚洲最大成人中文| 日韩高清综合在线| 国内少妇人妻偷人精品xxx网站 | 亚洲成人精品中文字幕电影| 欧美另类亚洲清纯唯美| 免费电影在线观看免费观看| 亚洲av中文字字幕乱码综合 | 国产成+人综合+亚洲专区|