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

    The regulation of memory effect and its influence on discharge properties of a dielectric barrier discharge driven by bipolar pulse at atmospheric-pressure nitrogen

    2021-10-31 08:15:26RuiFAN樊瑞YaogongWANG王耀功XiaoningZHANG張小寧ZhentaoTU屠震濤andJunZHANG張軍
    Plasma Science and Technology 2021年10期
    關(guān)鍵詞:張軍

    Rui FAN(樊瑞),Yaogong WANG(王耀功),?,Xiaoning ZHANG(張小寧),?,Zhentao TU (屠震濤) and Jun ZHANG (張軍)

    1 Key Laboratory of Physical Electronics and Devices of the Ministry of Education, Xi’an Jiaotong University, Xi’an 710049, People’s Republic of China

    2 School of Electronic Science and Engineering, Xi’an Jiaotong University, Xi’an 710049, People’s Republic of China

    Abstract The regulation of memory effect that the residual charges generated during and after discharge act on the initiation and development of subsequent discharge is explored by adjusting the pulse parameters, which have an influence on the discharge characteristics.The memory effect is quantified by the measurement of ‘wall voltage’ through a series of reference capacitors.The influences of memory effect on the discharge properties corresponding to rising/falling time 50–500 ns,pulse width 0.5–1.5 μs,and frequency 200–600 Hz are analyzed.It is found that the‘wall voltage’ increases from 1.4 kV to 2.4 kV with rising/falling time from 50 ns to 500 ns, it varies in the range of 0.18 kV with frequency of 200–600 Hz, and 0.17 kV with pulse width of 0.5–1.5 μs.The propagation velocity of wavelike ionization under the negative pulse slows down from 2184 km s?1 to 1026 km s?1 as the rising/falling time increases from 50 ns to 500 ns due to the weakening of the electric field by the surface memory effect.More intense and uniform emission can be achieved through faster rising/falling time and higher frequency based on the volume memory effect, while pulse width has less influence on the emission uniformity.Furthermore, similar laws are obtained for spectral and discharge intensity.Therefore, the memory effect is most effectively regulated by rising/falling time, and the discharge properties are affected by the surface and volume memory effect.

    Keywords: memory effect, discharge properties, wall voltage, pulse parameters

    1.Introduction

    Low-temperature plasma has been proven to be useful for pesticide elimination, surface treatment, gas purification, and laser technology [1–9].Dielectric barrier discharges (DBDs)with simple scalability [10] excited by a nanosecond pulse can generate low-temperature plasma with high electron density, high energy efficiency, and better uniformity [11].The collective behavior of filamentary and patterned DBDs can be influenced by residual species [12].The residual species and neutral active species (vibrationally excited molecules,metastable molecules,positive ions,etc)generated during the previous discharge and after discharge cease have an influence on subsequent discharge that is called the memory effect[13–16].The investigation of regulation of the memory effect and its influence on discharge characteristics is valuable for both an improvement of discharge efficiency and a reduction of sustain voltage [17].

    Memory effect agents include residual charges, neutral active particles, the detachment of negative ions, and the photon-assisted detachment of weakly bonded electrons on the dielectric surface[18,19].The memory effect universally facilitated the subsequent discharge initiation.Specifically,the reproducibility of the micro-discharge channel position is governed by the volume memory effect, and the surface memory effect is responsible for micro-discharge jittering in time[20].The electric field offset before the discharge pulse,streamer initiation, and propagation are also attributed to the memory effect [21–24].In addition, these discharge properties related to memory effects, such as electron avalanche development and pre-ionization distribution, can be modulated by voltage waveform parameters and dielectric material[25–29].High-energy electrons, discharge current density,electron density,and streamer propagation velocities,with the applied voltage,pulse repetition rate,pulse width,and rising/falling time were investigated[30–35].The propagation speed of ionization waves decreases as the rising time slows down in a micro-DBD at low gas pressure [36].In addition, the pulse polarity also has an impact on the properties of discharge in asymmetrical arrangements, while the polarity dependence of volume DBDs in symmetrical arrangement is not known [37–40].

    In summary, dominant memory effect agents are significantly affected by many factors, such as the voltage pulse waveform, gas pressure, gas composition, discharge history/mode, electric field inhomogeneity, and ‘pulse off’ period[13].Major memory effect agents, their influential mechanisms, and decay pathways have been summarized [41].Moreover, the discharge peculiarities associated with the memory effect at the effects of waveform parameters were studied [42–44], while the modulation of memory effect by pulse parameters needs to be explored.

    In this work, the memory effect is quantified by measuring the‘wall voltage’—a term introduced from the plasma display panel, which has the same structure as DBD—by using a reference capacitor.The ‘wall voltage’ with pulse parameters (rising/falling time, pulse width, and frequency)between the positive and negative half cycle of a bipolar pulse is investigated.Moreover, the influences of memory effect(surface and volume) on the discharge properties, such as ionization propagation, discharge intensity, and uniformity emission, with pulse parameters in a DBD device at atmospheric-pressure N2are preliminarily explored.

    The rest of the paper is organized as follows.Section 2 describes the experimental setup and the principle of the memory effect.Section 3 elaborates on the influence of pulse parameters on wall voltage,presents the experimental results,and discussions on the memory effect acting on discharge characteristics.Section 4 gives the conclusions.

    2.Experimental setup

    A schematic of the memory effect in the DBD device used in our experiment is indicated in figure 1.The memory effect discussed in this work mainly refers to the effect of the residual charge on the negative pulse after the discharge of positive pulse when the DBD is driven by a bipolar pulse.The bipolar high-voltage pulse is input to the left electrode of the DBD and the right electrode is grounded.After a positive pulseVswith amplitude higher than the ignition voltage acts on DBD,micro-discharge channels are formed in the gas gap.Positive charges and electrons are accumulated on the surface of the dielectric layer close to the cathode and anode,respectively.The formed voltage of accumulated charge is called the wall voltageVw.Then, in stage 1 (the interval between positive and negative pulses), the residual charge mainly exists in the former discharge channel and the surface of the dielectric layer.In stage 2, the negative pulse ?Vsis applied, and the remaining wall voltageVwis superimposed on the negative pulseVs, which affects the discharge in stage 2.

    In order to measure the wall voltage and investigate the influence of the memory effect on discharge properties with pulse parameters,the experimental system is set up as shown in figure 2.The plasma reaction chamber is pumped to 10?1kPa by a mechanical pump, and then the reaction chamber is filled with nitrogen to atmospheric pressure in order to analyze the discharge properties in a static nitrogen environment.The closed chamber is made of plexiglass.Two pieces of 1 mm thick indium tin oxide(ITO)glass(glass area 10×10 cm2, ITO area 5×5 cm2) are fixed together,forming a DBD device with a 3 mm gap between the two dielectric layers.The DBD device is excited by a parameterized nanosecond pulse power supply that can be programmed to output a bipolar pulse with adjustable pulse waveform[40,45].The reference capacitorCcis connected in series with the DBD to measure the wall voltage.The pulse voltages and discharge currents are measured with an oscilloscope (DSO90254A), a 1000:1 high-voltage probe (Tektronix P6015A) and a current probe (Magnelab FCT).In addition, optical properties are measured by a spectrometer(Ocean optics 2000), discharge photos are taken by a camera(Canon EOS 550D),and the morphologies of DBD discharge propagation are captured by an intensified charge-coupled device (ICCD) (Andor DH334).The electric fields of the discharge gap parallel to the electrode surface are denotedExandEy, and that perpendicular to the electrode surface isEz.

    The applied voltage, current waveform, and voltage waveform on the reference capacitorCcfor typical parameters(rising/falling time 50 ns,pulse width 1 μs,frequency 500 Hz)of the bipolar pulse used in the experiment are shown in figure 3.Without specific declaration,the default time interval between the positive and negative pulse is 8 μs, and the amplitude is 14 kV, while the rising/falling time can vary in the range of 50–500 ns,pulse width in the range of 0.5–1.5 μs,and frequency in the range of 200–600 Hz.The voltage waveform on the reference capacitorCcis measured as shown in figure 3(b) when the high-voltage pulse in figure 3(a) is applied to the DBD.According to the measured voltageVcwonCcduring the time interval between the positive and negative pulse, the wall voltageVwcan be calculated by using equation (1).The specific principle has been described in our previous work [45].

    whereCm=ε0εrS/dis assumed as an equivalent capacitor of the DBD,ε0is the absolute permittivity of a vacuum,and εris the relative permittivity of glass.

    3.Results and discussion

    3.1.Influence of pulse parameters on wall voltage

    The memory effect agents may mainly include the positive particles(N2+and N4+)and metastable particlesandin atmospheric-pressure N2;some typical reactions[14] are listed:

    By applying the measured equivalent capacitorCmof the DBD and reference capacitorCcwhich are 15.38 pF and 633 pF, respectively, the wall voltage can be calculated.The calculated wall voltages and the measured minimum sustain voltages change with the rising/falling time,pulse width,and frequency as shown in figure 4.According to figure 4(a), the wall voltage increases from 1.4 kV to 2.4 kV when the rising and falling time increases from 50 ns to 500 ns, as indicated by the purple circle curve.The ionization rates will cause changes in electron density and ion density, which leads to different time evolutions of the gap voltage and electric field for different rising and falling times.At the falling time of the positive pulse, the space charges will be recombined under the action of the electric field.The faster the falling time, the faster the electric field changes, resulting in a greater probability of recombination and lower residual wall charges,and vice versa.Therefore,as the falling edge slows down,there is more residual charge after the positive pulse, thus the wall voltage is greater.Furthermore, the minimum sustain voltage to maintain discharge varying with the rising and falling time is reflected by the gray square curve.It can be seen that the sustain voltage decreases from 11.7 kV to 9.1 kV when the rising and falling time increases from 50 ns to 500 ns.It also verifies that the wall voltage gradually increases with the decrease of the rising and falling time.

    Table 1.Maximum transferred charge qmax,residual charge qres,and energy with rising and falling times of 50 ns, 100 ns, 200 ns, and 500 ns.

    In addition,Q–VLissajous curves at rising/falling times of 50 ns, 100 ns, 200 ns, and 500 ns are indicated in figure 5.The maximum transferred chargeqmaxafter the positive pulse is calculated by theQ–VLissajous curves [46] and residual chargeqres, corresponding to rising/falling times of 50 ns,100 ns,200 ns,and 500 ns,as listed in table 1;moreover, the energyEof each pulse can be calculated by voltage–current–time integration[47],as shown in equation(2).It can be seen that theqmaxincreases from 842 nC to 922 nC,qresincreases from 24 nC to 255 nC when the rising and falling time increases from 50 ns to 500 ns, whileEdecreases.

    whereMis the record length of data, Δtis the sample interval, andVn,Inare the values of measured voltage and current at the recording pointn.

    In figure 4(b), the wall voltage remains nearly unchanged with pulse width within a certain range.When the pulse width changed from 0.5 μs to 1.5 μs, the wall voltage varied from 1.26 kV to 1.43 kV, and the amplitude of the sustain voltage changed from 10.8 kV to 10.4 kV.The pulse width refers to the pulse flat top period for the convenience of independent regulation of pulse parameters.Since the process of the accelerated electron collision, excitation ionization, and recombination mainly occurs during the rising and falling time, and the metastable particles e.g.during the flat top have little effect due to the small range of pulse width[48].In addition,the wall voltages and sustain voltages are almost unaffected by frequency of 200 Hz to 600 Hz; the wall voltages and sustain voltages vary in the range of 0.18–0.2 kV, respectively, as indicated in figure 4(c).The frequency has an enhanced effect on the discharge times, i.e.the cumulative discharge intensity;however,it has little effect on the wall charge accumulation in a single pulse period.Therefore, the memory effect of DBD discharge,namely the wall voltage,can be effectively regulated by the rising and falling time of the pulse in a certain range.

    3.2.Discharge propagation characteristics

    Due to the recombination of positive ions, the larger the interval between adjacent discharges, the more the influence of the memory effect may be reduced.In this work, the influence of residual charges and metastable particles generated by previous discharge on the rising edge of negative pulse is mainly studied.In nitrogen, the time-resolved, timeintegrated ICCD images of the DBD from the front and side view subjected to a bipolar pulse during the rising time of the negative pulse when the rising/falling time is 50 ns, 100 ns,200 ns,and 500 ns are studied in figures 6(a)–(d).The left part of figure 6 shows the voltage and current waveforms of the rising edge of the negative pulse.The middle part is the timeresolved images from the front view.The gate width of the ICCD detector is fixed at 3 ns, and a series of false color intensity images are captured in a continuous time range.The timet=0 ns is defined as the zero-crossing preceding the positive pulse of the bipolar pulse ofV(t).The four ICCD images correspond to the four moments of the rising edge of the negative pulse.For a fixed value oft, the obtained timeresolved images are the average of the corresponding image of 100 cycles of the pulse waveforms.The right part of figure 6 is the time-integral image from the front and side view for 100 cycles, i.e.0.2 s.

    A wavelike ionization wave is observed in the timeresolved images at the negative pulse phase of the bipolar pulse.The discharge starts from the upper left corner and propagates along the diagonal to the bottom right corner because of the lead position of the electrode, as shown in figure 2, and the discharge propagation presents a smooth periphery of the packet.As the rising and falling time slows down,the wavelike ionization waves become more dispersed until 500 ns.This mechanism of wavelike ionization wave propagation can be mainly referred to the combined action of electric field distortion, surface memory effect, and photoemission.The discharge is weakened which is caused by the reduction of the electric field distortion, as the rising and falling time slows down.On the other hand, the surface memory effect is enhanced with the increase of residual charge.The residual charge particles, combined with the applied electric field, produce the horizontal electric fieldEy(the electric field parallel to the electrode), as shown in figure 2, to move the ionization wave forward.TheEydecreases with the increase of the residual charges distributed on the surface of the dielectric layer, resulting in weaker propagation characteristics [49].

    The false color image (in 3D view) of DBD discharge propagation in nitrogen during the rising time of the negative pulse is indicated in figure 7(a).The propagation velocities of the first packet changing with the rising and falling time of the negative pulse along theYaxis marked in figure 7(a) are shown in figure 7(b).The propagation velocities along theYdirection refer to the propagation velocities of the peak discharge intensity in theYdirection calculated by MATLAB.The general trend is that the propagation speed decreases with the increase of the rising and falling time during the negative pulse, since in the case of fast rising time, the electric field distortion is more significant and ionization velocity is faster,while in the case of slow rising time, the electric field distortion alleviates and the ionization velocity is relatively lower.In addition, as the rising/falling time slows down, the surface memory effect is enhanced and the electric fieldEyalong the electrode plane weakens, which leads to the decrease of propagation velocity [49].Moreover, the propagation speed at the falling time of the negative pulse is lower than that at the rising time of the negative pulse, because the current of the falling time of the negative pulse is mainly affected by the residual charge particles generated by the applied electric field during the rising time.

    3.3.Emission uniformity and discharge intensity

    DBD discharge photos are taken by a camera (Canon EOS 550D)and the emission uniformity of DBD discharge photos is evaluated by using the nine-point display uniformity method, which is usually used to evaluate the uniformity of plasma display panel displays [45], and the rules of emission uniformity, discharge intensity with rising and falling time,pulse width, and frequency in nitrogen are studied, as represented in figure 8.The emission uniformity is reflected by the average deviation of the brightness of each test area and the average brightness of the nine areas; therefore, the larger the value of nonuniformity, the worse the uniformity.With the increase of rising/falling time, the emission uniformities and discharge intensities decrease gradually as shown in figure 8(a).The discharge intensity and emission uniformity are both affected by the electric field distortion and volume memory effect.They are calculated from the discharge photos, which are the cumulative process of the discharge.Since the micro-discharges have a spatial ‘memory,’ every subsequent micro-discharge appears at exactly the same location occupied by the micro-discharge at the preceding half-period, which is caused by the volume memory effect[20].The higher the wall voltage, the stronger the discharge intensity at the position of the previous discharge due to the volume memory effect, which may result in increased discharge intensity differences between regions, and worse uniformity.In addition, as the discharge intensity is the average value of the whole discharge area, it is mainly affected by the electric field.

    The influence of pulse width on the emission uniformity and discharge intensity in nitrogen is indicated in figure 8(b).With the increase of pulse width, the emission uniformity does not vary obviously when the pulse width changes from 0.5 μs to 1.5 μs.According to the above results, the wall voltage hardly changes with the pulse width in a certain range, as does the electric field.However, the discharge intensity decreases slightly with the increase of pulse width.The discharge at the rising time of the pulse is excited by the impact of an applied electric field,while at the falling time the discharge is the recombination of residual charged particles.Since the discharge intensity is a cumulative process, the increase of the pulse width will cause the decrease of residual particles during the falling time due to the recombination of positive ions, resulting in a small drop in the overall cumulative discharge intensity.In addition, figure 8(c) shows the influence of the pulse frequency on emission uniformity and intensity in nitrogen.With the increase of pulse frequency,the emission uniformity and discharge intensity gradually increase.The discharge photo is a process of integration.The higher the discharge frequency, the more discharge times per unit time, and the greater the discharge intensity.Furthermore, more discharge times will make the discharge in each area more uniform per unit time.Although the wall voltage does not change much with the frequency of 200–800 Hz, the micro-discharge will slightly shift between each discharge, so that the emission uniformity per unit time is improved.

    For further analysis, the discharge intensity gradually decreases with the increase of the rising and falling time because of weaker electric field, although the wall voltage gradually increases.Figure 9 shows the schematic diagram of electric field distribution considering the memory effect.When the subsequent pulse is applied to the DBD, the actual electric field of the DBD is the superposition of the applied electric field and the electric field generated by the wall charge; therefore, the intensity of the subsequent discharge is determined by the applied electric field and the wall voltage.Since the amplitude of the wall voltage is 14.3% of the applied voltage, the discharge intensity is mainly affected by the applied voltage.The shorter rising time means the plasma is driven under significant overvoltage breakdown, plasma should be more reactive and have higher electron temperature,and more atoms should be excited to higher energy states [33].

    Figure 1.Schematic of the DBD device and memory effect.

    Figure 2.Schematic of the experimental setup.

    Figure 3.(a)Applied voltage and current waveform of rising/falling time 50 ns and (b)voltage waveform on the reference capacitor under the action of the waveform in (a).

    Figure 4.Wall voltage and sustain voltage change with (a) rising and falling time, (b) pulse width, and (c) frequency.

    Figure 5. Q–V Lissajous curves with rising/falling times of 50 ns,100 ns, 200 ns, and 500 ns.

    Figure 6.Voltage and current waveform, time-resolved and time-integrated ICCD images during the rising edge of the negative pulse half cycle at slopes of (a) 50 ns, (b) 100 ns, (c) 200 ns, and (d) 500 ns.

    Figure 7.(a)3D process of DBD discharge and(b)the propagation velocity of the first packet change with the rising/falling time during the negative pulse along the Y direction.

    Figure 8.Uniformity and discharge intensity changes with (a) rising/falling time, (b) pulse width, and (c) frequency in nitrogen.

    Figure 9.Schematic of electric field distribution considering the memory effect.

    Figure 10.Spectra change with (a) rising/falling time, (b) pulse width, and (c) frequency in N2.

    The variations of spectrum (N2) with rising and falling time, pulse width, and frequency are shown in figures 10(a)–(c).The wavelength corresponding to the peak of the spectrum almost does not change with the rising and falling time,pulse width, and frequency, while the peak intensity of the spectrum gradually decreases with the increase of the rising and falling time.A faster rising and falling time means higher dv/dt, higher electron energy, more violent collisions between the particles, and more intense discharge.Although the wall voltage increases as the rising and falling time slows down, the wall voltage 1.4–2.4 kV is expected to be smaller than the applied voltage 14 kV, and the applied electric field plays the main role.As shown in figure 10(b), the peak intensity of the spectrum basically does not change with the increase of the pulse width.Since the accelerated collision and excitation processes of particles mainly occur in the rising and falling time, the pulse width has little effect on spectral intensity in a certain range.In addition, with the increase of frequency,the spectral peak intensity increases gradually.The measurement of the emission spectrum is also a cumulative process.With the increase of pulse frequency, the number of active particles excited per unit time will increase; therefore,the discharge intensity and spectral intensity will increase.At present, the wavelength corresponding to N2+(B2∏u+) greater than 350 nm is mainly measured.Although the intensity ratio of N2+(B2∏u+)/N2(C2∏u)varies with the pulse parameters,the effect of the pulse parameters on the spectral intensity of N2(C2∏u) qualitatively should be similar to that of N2+(B2∏u+); specifically, the shortening of the rising edge not only increases the spectral intensity of N2+(B2∏u+), but also the spectral intensity of N2(C2∏u).The discharge device will be improved to measure the influence of parameters on the spectral intensity during the ultraviolet band further.

    4.Conclusions

    By adding the reference capacitance,the quantification of the partial memory effect (i.e.wall voltage) between the positive and negative half cycle of a bipolar pulse has been indirectly measured, and the influence of pulse parameters including rising/falling time, pulse width, and frequency on memory effect has been analyzed.It has been shown that the rising/falling time plays a dominant role in wall voltage.The steeper the falling time of the positive half cycle, the more likely leads to the recombination and deexcitation of relatively more positive ions and neutral active particles,and the wall voltage is reduced, while the memory effect is almost unaffected by pulse width and frequency within a certain range(0.5–1.5 μs,and 200–600 Hz).Furthermore, the influence of memory effect on propagation properties, discharge intensity, and emission uniformity has been investigated at atmosphericpressure N2.The discharge properties are affected by the complex interaction of the applied electric field and memory effect.The mechanism of wavelike ionization wave in the rising edge of the negative half cycle is not affected by rising/falling time,while the propagation speed decreases due to the reduction of the electric field distortion and the effect of surface residual charge on the electric field in the direction of parallel electrodes.When the rising/falling time ranges from 50 ns to 500 ns, the applied electric field plays a dominant role in the discharge intensity; that is, the steeper the rising/falling time, the stronger the discharge intensity.Moreover,the discharge intensity decreases with the increase of pulse width and the decrease of frequency because of the memory effect aging.Due to the enhancement of the memory effect and attenuation of frequency, the difference of micro-discharge intensity in different positions increases, which leads to the decrease of the emission uniformity.In addition, the spectral intensity of N2+(B2∏u+)decreases with the increase of rising/falling time and the decrease of frequency, while the pulse width hardly affects the spectral intensity.Further, we will investigate the regulation of a wider range of pulse parameters on the memory effect and mechanism of action.

    Acknowledgments

    The authors would like to acknowledge the financial support provided by National Natural Science Foundation of China(Nos.51807156 and 61771382), Projects of International Cooperation and Exchanges Shaanxi Province(No.2018KW-034), China Postdoctoral Science Foundation (No.2017M623174) and Central University Basic Scientific Research Operating Expenses (No.xpt012019041).

    猜你喜歡
    張軍
    Enhanced measurement precision with continuous interrogation during dynamical decoupling
    篆刻:張軍
    萬(wàn)松浦(2023年6期)2023-04-08 11:57:16
    Generation of elliptical isolated attosecond pulse from oriented H in a linearly polarized laser field
    張軍從脾論治痛風(fēng)性關(guān)節(jié)炎緩解期的經(jīng)驗(yàn)總結(jié)
    世上沒(méi)有卑微的工作
    37°女人(2020年7期)2020-07-14 03:15:15
    張軍的牙
    故事大王(2020年4期)2020-05-03 13:59:17
    男子吃“霸王餐”逃跑摔傷索賠未果
    華聲文萃(2019年11期)2019-09-10 07:22:44
    A Quantum Secret Sharing Scheme Using Orbital Angular Momentum onto Multiple Spin States Based on Fibonacci Compression Encoding?
    誰(shuí)動(dòng)了我的壽限
    貪小便宜毀人生
    国产淫片久久久久久久久 | 国产成人啪精品午夜网站| 国产精品伦人一区二区| 99热这里只有是精品在线观看 | 婷婷精品国产亚洲av在线| 久久精品夜夜夜夜夜久久蜜豆| 一二三四社区在线视频社区8| 极品教师在线免费播放| 嫁个100分男人电影在线观看| 男人狂女人下面高潮的视频| 91九色精品人成在线观看| av福利片在线观看| 国产野战对白在线观看| 亚洲av电影在线进入| 美女cb高潮喷水在线观看| 亚洲欧美精品综合久久99| 国产在线男女| 欧美一区二区精品小视频在线| 一进一出好大好爽视频| 亚洲aⅴ乱码一区二区在线播放| 三级毛片av免费| 九色国产91popny在线| 国产探花极品一区二区| 深夜a级毛片| 亚洲av美国av| 国产精品久久视频播放| 亚洲国产精品久久男人天堂| 好看av亚洲va欧美ⅴa在| 免费在线观看影片大全网站| 亚洲熟妇熟女久久| 男人和女人高潮做爰伦理| 精品99又大又爽又粗少妇毛片 | 少妇熟女aⅴ在线视频| 免费看a级黄色片| 90打野战视频偷拍视频| 国产单亲对白刺激| 白带黄色成豆腐渣| 国产精品久久久久久精品电影| 亚洲aⅴ乱码一区二区在线播放| 中出人妻视频一区二区| 亚洲va日本ⅴa欧美va伊人久久| 人人妻人人澡欧美一区二区| 亚洲七黄色美女视频| 欧美中文日本在线观看视频| 激情在线观看视频在线高清| 一a级毛片在线观看| 狠狠狠狠99中文字幕| 亚洲精品456在线播放app | 成人三级黄色视频| 夜夜爽天天搞| 久久国产精品人妻蜜桃| 成人av一区二区三区在线看| 少妇的逼好多水| 1024手机看黄色片| 丁香六月欧美| 日本与韩国留学比较| 舔av片在线| 免费在线观看成人毛片| 亚洲va日本ⅴa欧美va伊人久久| 一级a爱片免费观看的视频| 观看美女的网站| 又爽又黄无遮挡网站| 欧美+日韩+精品| 国产午夜精品久久久久久一区二区三区 | 欧美极品一区二区三区四区| 亚洲专区国产一区二区| 亚洲av免费在线观看| 国产亚洲精品久久久com| 久久人妻av系列| 九九在线视频观看精品| 免费电影在线观看免费观看| 能在线免费观看的黄片| 在线观看舔阴道视频| 久久久久亚洲av毛片大全| 中文字幕熟女人妻在线| 精品乱码久久久久久99久播| 一区二区三区免费毛片| 中文字幕久久专区| 亚洲av成人不卡在线观看播放网| 乱人视频在线观看| 中文字幕免费在线视频6| 国产精品国产高清国产av| 在线观看一区二区三区| 在线观看午夜福利视频| 亚洲在线自拍视频| 1000部很黄的大片| or卡值多少钱| 欧洲精品卡2卡3卡4卡5卡区| 身体一侧抽搐| 超碰av人人做人人爽久久| 国产精品伦人一区二区| 亚洲人成网站在线播放欧美日韩| 国产成年人精品一区二区| 日韩亚洲欧美综合| 国产午夜福利久久久久久| 97超视频在线观看视频| 精品人妻熟女av久视频| 国产探花极品一区二区| 丰满人妻一区二区三区视频av| 久久久久九九精品影院| 脱女人内裤的视频| 久久久久亚洲av毛片大全| 在线观看美女被高潮喷水网站 | 三级国产精品欧美在线观看| 变态另类丝袜制服| 女人被狂操c到高潮| 国内精品一区二区在线观看| 美女黄网站色视频| 中文字幕精品亚洲无线码一区| 亚洲午夜理论影院| avwww免费| 久久久久国内视频| 麻豆一二三区av精品| 在线观看66精品国产| 色吧在线观看| 国产成人av教育| 特级一级黄色大片| 在线免费观看的www视频| 最近在线观看免费完整版| 真人做人爱边吃奶动态| 97人妻精品一区二区三区麻豆| 欧美黑人巨大hd| 级片在线观看| 欧美又色又爽又黄视频| 欧美日韩乱码在线| 国产色爽女视频免费观看| 麻豆久久精品国产亚洲av| 欧美性猛交黑人性爽| 性欧美人与动物交配| 老熟妇乱子伦视频在线观看| 在线免费观看的www视频| 别揉我奶头~嗯~啊~动态视频| 免费观看精品视频网站| 波野结衣二区三区在线| 人人妻人人澡欧美一区二区| 久久精品影院6| 一个人观看的视频www高清免费观看| 国产亚洲欧美98| 搡老熟女国产l中国老女人| 少妇被粗大猛烈的视频| 日本一二三区视频观看| 亚洲av中文字字幕乱码综合| 在线观看一区二区三区| 日韩 亚洲 欧美在线| 亚洲,欧美,日韩| 精品国产三级普通话版| 欧美xxxx性猛交bbbb| 欧美绝顶高潮抽搐喷水| 88av欧美| 成人鲁丝片一二三区免费| 国产精品98久久久久久宅男小说| 久久久久久大精品| 精品99又大又爽又粗少妇毛片 | 亚洲av成人av| 婷婷丁香在线五月| 亚洲成人久久性| 99国产综合亚洲精品| 欧美极品一区二区三区四区| 麻豆久久精品国产亚洲av| 欧美性感艳星| 亚洲欧美日韩高清在线视频| 91av网一区二区| 中亚洲国语对白在线视频| 精品国内亚洲2022精品成人| 精品国产三级普通话版| 国内精品久久久久久久电影| 99视频精品全部免费 在线| 国产三级黄色录像| 亚洲人成网站在线播| 久久精品91蜜桃| 九九热线精品视视频播放| 精品日产1卡2卡| 欧美最新免费一区二区三区 | 99riav亚洲国产免费| 久久久久精品国产欧美久久久| 91久久精品电影网| 日本 欧美在线| 中文字幕久久专区| 国产主播在线观看一区二区| 99精品久久久久人妻精品| 国产免费男女视频| 国产三级黄色录像| 国产av一区在线观看免费| 又粗又爽又猛毛片免费看| 亚洲aⅴ乱码一区二区在线播放| 午夜精品一区二区三区免费看| av天堂在线播放| 最近最新免费中文字幕在线| 欧美日本视频| 五月玫瑰六月丁香| 蜜桃久久精品国产亚洲av| 色哟哟哟哟哟哟| 欧美日韩乱码在线| 免费一级毛片在线播放高清视频| 黄色女人牲交| 别揉我奶头~嗯~啊~动态视频| 99国产精品一区二区三区| 国产私拍福利视频在线观看| 如何舔出高潮| 欧美+亚洲+日韩+国产| 免费搜索国产男女视频| 免费人成在线观看视频色| 亚洲国产精品sss在线观看| 高清毛片免费观看视频网站| 人妻制服诱惑在线中文字幕| 国产一区二区激情短视频| 精品午夜福利视频在线观看一区| 午夜福利在线观看免费完整高清在 | 乱码一卡2卡4卡精品| av天堂在线播放| 亚洲精品久久国产高清桃花| 老司机深夜福利视频在线观看| 丰满人妻一区二区三区视频av| 桃色一区二区三区在线观看| 国产免费一级a男人的天堂| 91在线精品国自产拍蜜月| 国产高清三级在线| av欧美777| 可以在线观看的亚洲视频| 激情在线观看视频在线高清| 国产精品女同一区二区软件 | 亚洲中文字幕一区二区三区有码在线看| 少妇裸体淫交视频免费看高清| 观看免费一级毛片| 国产主播在线观看一区二区| 亚洲精品456在线播放app | 国产一区二区在线av高清观看| 老司机福利观看| 国产爱豆传媒在线观看| 最后的刺客免费高清国语| 又黄又爽又免费观看的视频| 亚洲狠狠婷婷综合久久图片| 国产探花在线观看一区二区| 91麻豆精品激情在线观看国产| 亚洲欧美激情综合另类| 成人毛片a级毛片在线播放| 久久6这里有精品| 午夜视频国产福利| 亚洲精品色激情综合| 又爽又黄无遮挡网站| 我的老师免费观看完整版| 美女cb高潮喷水在线观看| 欧美绝顶高潮抽搐喷水| 淫秽高清视频在线观看| 脱女人内裤的视频| 久久久久久久久久黄片| 久久久久精品国产欧美久久久| 国产野战对白在线观看| 美女 人体艺术 gogo| 又紧又爽又黄一区二区| 免费大片18禁| 免费一级毛片在线播放高清视频| 日本精品一区二区三区蜜桃| 一本精品99久久精品77| 久久精品久久久久久噜噜老黄 | 欧美日本视频| 久久久久久久午夜电影| 夜夜躁狠狠躁天天躁| 精品久久久久久久久亚洲 | 中文字幕av在线有码专区| 精品国产三级普通话版| 成年女人毛片免费观看观看9| 亚洲成人久久性| 久久午夜福利片| 成人永久免费在线观看视频| 噜噜噜噜噜久久久久久91| 国产色爽女视频免费观看| 国产一区二区在线观看日韩| 国产精品伦人一区二区| 亚洲乱码一区二区免费版| 99在线视频只有这里精品首页| 国产乱人伦免费视频| 嫩草影院精品99| 怎么达到女性高潮| 免费电影在线观看免费观看| www.熟女人妻精品国产| 久久亚洲精品不卡| 一卡2卡三卡四卡精品乱码亚洲| 好男人电影高清在线观看| 最后的刺客免费高清国语| 国产精品乱码一区二三区的特点| 日韩欧美免费精品| 国产 一区 欧美 日韩| 大型黄色视频在线免费观看| 制服丝袜大香蕉在线| 中文字幕人妻熟人妻熟丝袜美| 特大巨黑吊av在线直播| 成人欧美大片| 国产一区二区在线观看日韩| 久久性视频一级片| 国产精品一区二区三区四区免费观看 | 午夜激情福利司机影院| 一个人看视频在线观看www免费| 久久精品国产自在天天线| 欧美成狂野欧美在线观看| 欧美日本视频| 日韩欧美三级三区| 久久久久性生活片| 亚洲一区二区三区不卡视频| 亚洲欧美日韩东京热| 国产伦精品一区二区三区视频9| 别揉我奶头~嗯~啊~动态视频| 激情在线观看视频在线高清| 久久99热这里只有精品18| 欧美日韩福利视频一区二区| 欧美潮喷喷水| 丁香欧美五月| 中亚洲国语对白在线视频| 脱女人内裤的视频| 97碰自拍视频| 免费观看人在逋| 欧美色欧美亚洲另类二区| 在线十欧美十亚洲十日本专区| 丁香欧美五月| 俺也久久电影网| 窝窝影院91人妻| 亚洲精品在线观看二区| 久久亚洲真实| 日韩大尺度精品在线看网址| 久久精品人妻少妇| 亚洲av一区综合| 国产av一区在线观看免费| 国产乱人视频| 天美传媒精品一区二区| 蜜桃久久精品国产亚洲av| 麻豆久久精品国产亚洲av| 一级毛片久久久久久久久女| 麻豆av噜噜一区二区三区| 国语自产精品视频在线第100页| 亚洲国产欧美人成| 看免费av毛片| 久久久久国产精品人妻aⅴ院| 在线播放无遮挡| 欧美在线黄色| 欧美成人免费av一区二区三区| 99国产极品粉嫩在线观看| 熟女人妻精品中文字幕| 中文字幕av成人在线电影| 久久午夜亚洲精品久久| 两个人的视频大全免费| 久久久国产成人免费| 永久网站在线| 在现免费观看毛片| 两个人的视频大全免费| 国产成人aa在线观看| 12—13女人毛片做爰片一| 久久性视频一级片| 麻豆久久精品国产亚洲av| 国产成人影院久久av| 亚洲最大成人中文| 亚洲黑人精品在线| 国内精品一区二区在线观看| 成人午夜高清在线视频| 国语自产精品视频在线第100页| 久久6这里有精品| 伊人久久精品亚洲午夜| 91在线精品国自产拍蜜月| 内射极品少妇av片p| 18禁在线播放成人免费| 国产成+人综合+亚洲专区| 国产精品精品国产色婷婷| 99久久99久久久精品蜜桃| 久久久久久大精品| 91狼人影院| 中文字幕熟女人妻在线| 天堂网av新在线| 亚洲激情在线av| 久久草成人影院| 97人妻精品一区二区三区麻豆| 757午夜福利合集在线观看| 亚洲激情在线av| 尤物成人国产欧美一区二区三区| 久久这里只有精品中国| 国产大屁股一区二区在线视频| 久久久成人免费电影| 欧美高清性xxxxhd video| 国产高清激情床上av| 国产精品久久电影中文字幕| www.色视频.com| 国产久久久一区二区三区| 午夜a级毛片| 免费黄网站久久成人精品 | 日本撒尿小便嘘嘘汇集6| 动漫黄色视频在线观看| 成人av在线播放网站| 亚洲第一区二区三区不卡| 精品久久久久久久末码| 九九在线视频观看精品| 亚洲欧美清纯卡通| 人妻夜夜爽99麻豆av| 精品欧美国产一区二区三| 精品国产亚洲在线| 亚洲国产色片| 欧美又色又爽又黄视频| 日本黄色片子视频| 一级a爱片免费观看的视频| 欧美成人免费av一区二区三区| 国产一区二区在线观看日韩| 高清在线国产一区| 五月玫瑰六月丁香| 美女被艹到高潮喷水动态| 欧美一区二区国产精品久久精品| 在线国产一区二区在线| 欧美性感艳星| 欧美绝顶高潮抽搐喷水| 毛片一级片免费看久久久久 | 国内精品美女久久久久久| 1024手机看黄色片| 波多野结衣高清无吗| x7x7x7水蜜桃| 美女免费视频网站| 日本在线视频免费播放| www.色视频.com| 观看免费一级毛片| 国产三级黄色录像| 在线观看舔阴道视频| 国产伦人伦偷精品视频| 国语自产精品视频在线第100页| 成人无遮挡网站| 成年女人看的毛片在线观看| 国产淫片久久久久久久久 | 青草久久国产| 老熟妇仑乱视频hdxx| bbb黄色大片| 亚洲真实伦在线观看| 性欧美人与动物交配| 国产一区二区在线观看日韩| 一本精品99久久精品77| 国产成人福利小说| 91在线观看av| 69人妻影院| 婷婷精品国产亚洲av| 高清毛片免费观看视频网站| 久久久久久九九精品二区国产| 亚洲人成伊人成综合网2020| 五月伊人婷婷丁香| 尤物成人国产欧美一区二区三区| 国产国拍精品亚洲av在线观看| 国产av麻豆久久久久久久| 亚洲av电影在线进入| 一级av片app| 91av网一区二区| 日本三级黄在线观看| 国产免费av片在线观看野外av| 久久久久久久亚洲中文字幕 | 日韩精品中文字幕看吧| 嫁个100分男人电影在线观看| 成人av一区二区三区在线看| 日本黄色片子视频| 久久99热这里只有精品18| 国产成人a区在线观看| 少妇的逼好多水| 国产熟女xx| 欧美在线一区亚洲| 永久网站在线| 欧美色欧美亚洲另类二区| eeuss影院久久| 尤物成人国产欧美一区二区三区| 国产黄a三级三级三级人| 午夜日韩欧美国产| 狠狠狠狠99中文字幕| 在线观看免费视频日本深夜| 亚洲国产色片| 亚洲欧美日韩无卡精品| 日本在线视频免费播放| 亚洲精品粉嫩美女一区| 国产探花极品一区二区| 一本久久中文字幕| 波多野结衣高清无吗| 精品午夜福利视频在线观看一区| 国产精品,欧美在线| 99久久精品一区二区三区| 深爱激情五月婷婷| 三级男女做爰猛烈吃奶摸视频| 一本一本综合久久| 久久人人精品亚洲av| 精品久久久久久成人av| 久久精品影院6| 色5月婷婷丁香| 偷拍熟女少妇极品色| 真人做人爱边吃奶动态| 久久久久久久午夜电影| 日本熟妇午夜| 久久久久久久精品吃奶| 日韩精品中文字幕看吧| 99久久无色码亚洲精品果冻| 熟女人妻精品中文字幕| 18禁黄网站禁片免费观看直播| 久久久久精品国产欧美久久久| 九色成人免费人妻av| 久久久久久久久久成人| 男人的好看免费观看在线视频| 亚洲 国产 在线| 成人亚洲精品av一区二区| 日本精品一区二区三区蜜桃| 国产精品亚洲美女久久久| 亚洲精华国产精华精| ponron亚洲| 精品日产1卡2卡| 精品午夜福利在线看| 国产精品不卡视频一区二区 | 人妻夜夜爽99麻豆av| 美女高潮的动态| 免费av毛片视频| 少妇的逼水好多| av视频在线观看入口| 国产一区二区三区视频了| 国产69精品久久久久777片| 久久久久久久久大av| 免费一级毛片在线播放高清视频| 69av精品久久久久久| 日本黄大片高清| 日本黄色片子视频| 精品乱码久久久久久99久播| 国产午夜精品论理片| 99riav亚洲国产免费| 久久欧美精品欧美久久欧美| 99热这里只有精品一区| 国产一区二区激情短视频| 成人国产综合亚洲| 午夜精品久久久久久毛片777| 五月伊人婷婷丁香| 99热精品在线国产| 亚洲aⅴ乱码一区二区在线播放| 人妻夜夜爽99麻豆av| 如何舔出高潮| 国产成人福利小说| 免费无遮挡裸体视频| 国产精品美女特级片免费视频播放器| 国产不卡一卡二| 免费观看精品视频网站| 欧美zozozo另类| 一进一出抽搐gif免费好疼| 亚洲欧美精品综合久久99| 亚洲av成人av| 成人国产综合亚洲| 一进一出抽搐动态| 草草在线视频免费看| 丁香六月欧美| 夜夜爽天天搞| 神马国产精品三级电影在线观看| 日本在线视频免费播放| 中文字幕熟女人妻在线| 熟女电影av网| 又爽又黄无遮挡网站| 天堂动漫精品| 黄色一级大片看看| 亚洲人成电影免费在线| 久久精品夜夜夜夜夜久久蜜豆| 美女被艹到高潮喷水动态| 欧美在线一区亚洲| 嫩草影视91久久| 丰满人妻熟妇乱又伦精品不卡| 九九热线精品视视频播放| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美激情综合另类| 日韩欧美三级三区| 日本熟妇午夜| 国产精品美女特级片免费视频播放器| 久久精品国产亚洲av香蕉五月| 18+在线观看网站| 俺也久久电影网| 精品人妻视频免费看| 99热这里只有是精品50| av视频在线观看入口| 日本免费一区二区三区高清不卡| 国内久久婷婷六月综合欲色啪| 中文字幕熟女人妻在线| 日本黄色视频三级网站网址| 成年女人毛片免费观看观看9| 免费高清视频大片| 又黄又爽又免费观看的视频| 狂野欧美白嫩少妇大欣赏| 香蕉av资源在线| 亚洲一区二区三区不卡视频| 欧美色视频一区免费| 午夜福利欧美成人| 亚洲中文日韩欧美视频| 午夜免费激情av| 亚洲成人中文字幕在线播放| 精品国内亚洲2022精品成人| 亚洲精华国产精华精| 日本与韩国留学比较| 人妻丰满熟妇av一区二区三区| 两个人视频免费观看高清| 午夜福利高清视频| 久久精品国产亚洲av天美| av在线天堂中文字幕| 三级国产精品欧美在线观看| 精品99又大又爽又粗少妇毛片 | 美女高潮的动态| 久久人人爽人人爽人人片va | 午夜亚洲福利在线播放| 欧美一区二区精品小视频在线| 国产av一区在线观看免费| 日本免费a在线| 国产精品综合久久久久久久免费| 日本三级黄在线观看| 亚洲av一区综合| 老熟妇仑乱视频hdxx| 国产在视频线在精品| 亚洲七黄色美女视频| 欧美+亚洲+日韩+国产| 国产 一区 欧美 日韩| av视频在线观看入口| 亚洲精品在线观看二区| 国产色婷婷99| 午夜福利在线观看吧| 亚洲欧美日韩无卡精品| 天堂av国产一区二区熟女人妻| 97热精品久久久久久| av国产免费在线观看| 波多野结衣巨乳人妻| 给我免费播放毛片高清在线观看| 精品国产三级普通话版|