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

    Effects of the transverse electric field on nanosecond pulsed dielectric barrier discharge in atmospheric airflow

    2020-06-14 08:45:30YongfengXU徐永鋒HongfeiGUO郭宏飛YuyingWANG王玉英ZhihuiFAN樊智慧andChunshengREN任春生
    Plasma Science and Technology 2020年5期
    關鍵詞:春生智慧

    Yongfeng XU (徐永鋒),Hongfei GUO (郭宏飛),Yuying WANG (王玉英),Zhihui FAN (樊智慧) and Chunsheng REN (任春生),?

    1 Key Laboratory of Materials Modification by Laser,Ion,and Electron Beams,Dalian University of Technology,Dalian 116024,People’s Republic of China

    2 School of Science,Hebei University of Engineering,Handan 056038,People’s Republic of China

    ? Deceased.

    Abstract

    Keywords:dielectric barrier discharge,high-voltage,atmospheric pressure air,airflow,lowtemperature plasma,nanosecond pulsed discharge

    1.Introduction

    Dielectric barrier discharge (DBD) is a typical method for generating non-equilibrium plasma at atmospheric pressure.Due to its moderate electron temperature and plasma density,it has been widely used for a variety of industrial applications,such as ozone generation,environmental treatment,surface modification,and plasma display panels[1–5].DBD is usually driven by a pulse power supply or alternating current (AC)power source with a frequency from tens of hertz to several tens of kilohertz.Compared with the DBD excited by an AC power source,the nanosecond pulsed DBD is more diffuse[6],which is mainly ascribed to the over-voltage breakdown based on the fast rising time in tens of nanosecond scale [7,8].This results in a larger reduced electric field E/N (where E is the electric field,N is the density of neutral particles),which can generate higher energy electrons and numerous electron avalanches simultaneously [7,9].In addition,nanosecond pulses can improve the discharge stability when the narrow duration of pulses is shorter than the characteristic time needed for the development of discharge instability [10].Therefore,a nanosecond pulse power supply is used to excite the discharge in our experiment due to the above advantages.

    Although a nanosecond pulse power supply makes it possible to realize uniform discharge under strict conditions,such as a narrow gap[8],DBD usually presents a filamentary discharge.The reason is that the Pd value (product of gas pressure P and gas gap d) is so large at atmospheric pressure that the local space charge field formed by the space charge separation can be compared to the applied electric field [11].However,our previous work reported that the diffuse discharge started from the electrode edge and moved toward the downstream region under the effect of airflow,which led to the diffuse discharge in the whole discharge space[12].It was found that there was a filamentary discharge in the central region of the electrode and a diffuse discharge at the electrode edge.This was mainly attributed to the existence of the edge effect [13–15].Most researchers believe that the diffuse discharge at the electrode edge is mainly ascribed to the stronger electric field or the local field enhancement effect [16,17].For example,it was found by Okazaki et al that uniform discharge can be generated by using metal wire meshes as electrodes covered by polyethylene terephthalate,and the discharge was affected by the flatness of the metal wire meshes [18].Other works have also reported that uniform discharge could be obtained under the same wire mesh electrodes [19,20].They pointed out that the generation of uniform discharge was ascribed to the field enhancement which initiated a corona discharge and provided seed electrons.Liu et al found that the floating electrode could improve the uniformity of discharge due to the stronger electric field,which is beneficial for generating the initial discharge,and resulted in sufficient and evenly distributed initial electrons[21].It has been found that a bright ring first appears at the electrode edge and then moves from the edge to the inside,and the result agrees well with a 2D simulation that showed that the electric field at the electrode edge increased and resulted in a more intense discharge[22].Duan et al reported that the pattern shape was consistent with the shape of the electrode.But they did not provide an explanation for the role of the edge effect in pattern formation [13].For the diffuse discharge at the electrode edge,less attention has been paid to the influence of the transverse electric field on discharge[12,23].It has been reported in previous work that the transverse electric field at the electrode edge is responsible for the generation of diffuse discharge[12].Therefore,the effects of the transverse electric field on discharge are investigated.

    In this paper,an asymmetric electrode geometry is proposed to enhance the transverse electric field.To further explore the influence of the mechanisms of the transverse electric field on discharge,the electrical characteristics,discharge images,movement characteristics of the micro-discharge channel,and optical emission spectra (OES) are analyzed.

    2.Experimental setup

    Figure 1.Schematic diagram of the experimental setup.(a) Symmetric geometry where the high-voltage and grounded electrodes are aligned,(b) asymmetric geometry where the transverse electric field is in the same direction as the airflow,(c) asymmetric geometry where the transverse electric field is opposite to the airflow.

    A schematic diagram of the experimental setup is presented in figure 1.The working gas is room air.The electrode geometry(a)is a symmetric electrode structure.Electrode geometries (b) and(c) are the asymmetric electrode structures in order to obtain a transverse electric field component.The only difference between them is the direction of the transverse electric field.The highvoltage and grounded electrodes are in the x-y plane,and the electrodes are only staggered in the y-direction.Therefore,the transverse electric field component is along the y-direction,and the longitudinal electric field component is along the z-direction.A transparent indium tin oxide (ITO) and an aluminum foil,covered by quartz plates with a thickness of 2 mm,are used as the high-voltage and grounded electrodes,respectively.The aligned area of the high-voltage and grounded electrodes is 100 × 50 mm2.The displacement distance between the highvoltage and grounded electrodes is 20 mm.The discharge gaps are varied from 4–6 mm.A repetitive nanosecond pulse power supply applied to the high-voltage electrode can generate pulse voltage with a rising time of 40 ns and a full width at half maximum of 200 ns,which is used to ignite the discharge.The voltage regulator connected to the nanosecond pulse power supply can be used to adjust the voltage amplitude input to the DBD circuit.The pulse repetition frequencies (PRFs) of the nanosecond pulse power supply are adjusted at the fixed values of 100,300,600,1000,and 1200 Hz.In all cases,the discharges are excited by the same voltage waveform and amplitude of 31.5 kV.The applied voltage and total current are measured by a high-voltage probe(Tektronix P6015 A,bandwidth 75 MHz)and a current probe (Pearson 6600,bandwidth 120 MHz),respectively.And a digital oscilloscope (Tektronix DPO4104,bandwidth 1 GHz,sampling frequency 5 GS s-1)is used to record and store these waveforms.The OES are directly measured by the spectrometer equipped with a CCD detector (Acton Research Spectrum 2500i),which is calibrated with a Hg lamp.The entrance slit opening of the monochromator is 500 μm.The resolution of the spectrum measurement system is 0.013 nm.The grating groove is 2400 lines/mm with a glancing wavelength at 240 nm.The integration time is 20 ms.The distance between the optical fiber and the DBD reactor is fixed at 35 mm.The collimator is placed on the fiber head to ensure that the OES can be measured at a local area and realized spatial resolution.Location 1 at the central region of the electrode is selected to represent the aligned region,location 2 where the fiber is tangent to the electrode edge represents the misaligned region,and location 3 indicates the electrode edge.The reason will be given in section 3.1.3.The airflow system consists of an air pump serving as the source of airflow and a self-designed wind tunnel used for introducing airflow to the discharge gap.A pitot tube connected with a pressure differential gauge(Testo512,200 hPa)is placed at the central position of the wind tunnel outlet.The airflow velocity can be changed and displayed on the pressure differential gauge by adjusting the valve on the wind tunnel,and the valve positions corresponding to the different airflow velocities are marked in advance.Airflow parallel to the electrode along y-direction ranges from 0–40 m s-1.Accordingly,the transverse electric field of the electrode geometry (b) is in the same direction as the airflow and the transverse electric field of the electrode geometry (c) is opposite to the airflow.Top-view and side-view discharge images are captured by a digital camera(Nikon D7000) and a high speed motion analyzer (Phantom V1610).The side-view discharge images are taken in the x-direction.

    Figure 2.The applied voltage and total current waveforms of different electrode structures under a 5 mm discharge gap and a 1000 Hz PRF.(1) Electrode geometry (a),(2) electrode geometry (b),(3) electrode geometry (c).

    Figure 3.Side-view images of the discharge of different electrode structures.The exposure time is 1/1000 s.(1) Electrode geometry (a),(2) electrode geometry (b),(3) electrode geometry (c).

    3.Results and discussion

    3.1.Comparison of the discharges between the different electrode structures in static air

    3.1.1.Electrical characteristics.When the discharge gap is 5 mm and the PRF is 1000 Hz,the applied voltage and total current waveforms of different electrode structures are presented in figure 2.It can be seen that the applied voltage and total current waveforms of different electrode structures are almost the same.A sudden drop at the rising front of the applied voltage is about 20 ns.It is usually considered a sign of gas breakdown,which is ascribed to the sudden increase of the discharge current and insufficient output power of the nanosecond supply when the gas breaks down[24].The total current peak occurring at around 48 ns is about 166 A.The total current is composed of the conduction current,displacement current,and electromagnetic interference [25].The displacement current and electromagnetic interference exhibit little difference due to a less misaligned region compared with the aligned region.And the conduction current is mainly contributed by the discharge of the aligned region,and accordingly,the conduction current is almost identical due to the same aligned regions between the different electrode structures.

    Figure 4.Simulated spatial distributions of the electric field intensity of different electrode structures.(1) Electrode geometry (a),(2) electrode geometry (c).

    3.1.2.Discharge images.Although the electrical characteristics are almost the same for the different electrode structures,there is a difference in the side-view discharge images given in figure 3 under the same experimental conditions.The DBD reactor is in static air.The exposure time is 1/1000 s,which is equal to one discharge cycle.The edge locations of the high-voltage and grounded electrodes are marked with red lines and green lines,respectively.It can be seen that the discharge at the electrode edge is more diffuse than that at the electrode center for the symmetric electrode geometry,and the diffuse discharge manifests in the misaligned region for the asymmetric electrode geometries.A filamentary discharge is composed of numerous short-lived micro-discharge channels,which are randomly distributed in space and time.It is found that the diffuse discharge has a more uniform emission or luminescence without any bright spots or an irregular distribution pattern compared with the filamentary discharge,or it means that no bright filaments can be distinguished.The difference between the filamentary discharge and the diffuse discharge is visually distinguishable.In general,for the AC discharge,the current waveform can be used to categorize the different regimes of DBD[18].In this case the current waveform of the filamentary discharge shows multi-current peaks,while the current waveform of the uniform discharge presents a current peak.However,for the nanosecond pulsed DBD,the different discharge states cannot be distinguished by using electrical properties,because the current waveform shows only one current peak for both the filamentary discharge and the diffuse discharge.The generation of the diffuse discharge is attributed to the more uniform distribution of the space charges induced by the transverse electric field.For the symmetric electrode geometry,on the one hand,the electric field at the electrode edge is stronger than that in the central region [17,22],which can generate more electron avalanches.On the other hand,it has been reported that there is a transverse electric field component at the electrode edge [12,23].The spatial ions in the microdischarge channels at the electrode edge have a transverse movement under the effect of the transverse electric field,and drag the surface electrons through an electrostatic coupling force,which leads to the redistribution of the space charges.The combined effects,which are the simultaneous generation of more electron avalanches and the redistribution of the space charges,result in the spatial overlap of the electron avalanches and the smoothness of the space charge field gradient.Thus the excessive development of single electron avalanches into the streamer can be inhibited,which is conducive to the generation of diffuse discharge [26].The reasons for the presence of the diffuse discharge in the misaligned region of the asymmetric electrode geometries are analogous.

    To certify that the generation of the diffuse discharge is caused by the transverse electric field,the spatial electric field intensity distributions of electrode geometries (a) and (c) are simulated by using Ansoft Maxwell software,as shown in figure 4.This method has been used in previous works[21,27,28].In this model,all the simulation parameters are the same as those in the experiment.The high-voltage and grounded electrodes are applied with voltages of 31.5 kV and 0 V,respectively.It can be seen that the electric field intensity distributions at the electrode edges and misaligned regions are not uniform,and there is a transverse electric field component in these regions.Although it is difficult to determine the quantitative intensities of the electric field and the transverse electric field component,the relative intensity between them can be qualitatively compared.They correspond to the phenomena of the discharge images in figure 3,which confirm our analysis of the experimental results.

    Figure 5.Temporal evolution of the discharge with a 1/1000 s exposure time.PRF=1000 Hz,discharge gap=6 mm.1 T represents the first discharge cycle,and the time interval is 1 T=1 ms.(1)Electrode geometry(a),(2)electrode geometry(b),(3)electrode geometry(c).

    In addition,it is concluded from figure 3 that the length of the diffuse discharge region along the direction of the transverse electric field for the asymmetric electrode geometries is longer than that for the symmetric electrode geometry.The lengths of the diffuse discharge region are integrated directly into the captions of figure 3.The length of the diffuse discharge region is considered to be related to the moving velocity of the space charges.Therefore,the moving characteristics of the micro-discharge channel are studied by using the discharge temporal evolution images taken by a high speed motion analyzer.The results are shown in figure 5 with a 6 mm discharge gap.The exposure time is 1/1000 s and the sample rate is 1000 frames per second (fps) for the 1000 Hz PRF.The starting time is 1 T,representing the first discharge.The time interval of 1 T is equal to 1 ms.The cycle corresponding to each discharge image is marked.The two micro-discharge channels marked with red squares are selected and traced frame by frame in figure 5,which are both located at the high-voltage electrode edge.According to the invariance of the relative position between the two microdischarge channels,the two micro-discharge channels marked with red squares in each cycle are considered to be identical.Therefore,the positions of the two micro-discharge channels in the sequential discharge images can be identified.First,the temporal evolution of the discharge images show that the bright micro-discharge filaments manifest in the aligned region,while there is no diffuse discharge in the misaligned region,which is consistent with figure 8 in section 3.2.The diffuse discharge is not obvious,and the relatively weak luminescence intensity of the diffuse discharge cannot be photographed under a larger discharge gap,which leads to no diffuse emission.Then it can be found from figure 5 that the micro-discharge channel of the symmetric electrode geometry is stationary,and the micro-discharge channel moves along the direction of the transverse electric field for the asymmetric electrode geometries.It is found that the movement of the micro-discharge channel is the migration from a point to another in the adjacent discharge cycle [12].For the symmetric electrode geometry,the residual charges generated by the previous discharge always maintain the same position,and the subsequent micro-discharge channel remains at the same location where the residual charges occurred due to the memory effect[5,29,30].The memory effect refers to the influence of the residual charges that survived the previous discharge on the subsequent discharge [31].Although the lifetime of the micro-discharge channel was on the order of 100 ns or less [32,33],it took about 30 min to remove the surface charges [34].The space charges that survived the previous discharge were abundant when the next discharge came [29].Thus the micro-discharge channel appears to be stationary.For the asymmetric electrode geometries,the transverse electric field in the misaligned region as the driving force breaks the equilibrium of the forces between the microdischarge channels,and the repulsive force on the space charges of the aligned region compels the movement along the direction of the transverse electric field.Accordingly,the memory effect is broken.When the next pulse arrives,the micro-discharge channel occurs at the new position where the space charges migrate.The space charges are almost stationary due to the memory effect of the symmetric electrode geometry.However,the space charges move directionally under the effect of the transverse electric field force for the asymmetric electrode geometries.Consequently,there is a longer distance of transverse migration of the space charges during the same time interval,which results in the increase of the length of the diffuse discharge region along the direction of the transverse electric field.

    In summary,the transverse electric field is responsible for the generation of the diffuse discharge.The length of the diffuse discharge region along the direction of the transverse electric field is related to the migration of the space charges induced by the transverse electric field.

    3.1.3.Optical emission spectra.To estimate the plasma gas temperature and electron density,the OES of the three locations are measured with a 5 mm discharge gap and a 1000 Hz PRF.It can be concluded from figure 3 that the discharges at the electrode edges and the misaligned electrode regions are more diffuse than those in the central region.Therefore,the three locations shown in the experimental setup are selected to measure the OES.In order to reduce errors,each group of data is recorded three times under the same conditions.According to a simplified collision-radiation model,the relative intensity ratio of 371.1 and 380.5 nm(I371.1nm/I380.5nm) can be used to calculate the average electron density ne[35–37].The relationship between

    Figure 6.The OES at three locations with a 5 mm discharge gap and a 1000 Hz PRF.

    where C0,C1,and C2are constant,n is the degree of ionization,and Ng=2.5 × 1019cm–3is the number density of the neutral particles at atmospheric pressure.According to the OES shown in figure 6,the range of I371.1nm/I380.5nmat locations 1,2,and 3 is 0.109 ± 0.003,0.107 ± 0.001,and 0.110 ± 0.001,respectively.This is similar to the previous work where the value of I371.1nm/I380.5nmis about 0.118,and the calculated neis about 1.72 × 1011cm–3[38].It can be concluded that the electron density at the electrode edge is slightly higher.When the electric field at the electrode edge is stronger,electrons collide and ionize neutral particles to generate more electrons.

    Due to a small energy gap between the rotational levels,the equilibrium between the translational motion and the rotational motion is easily achieved by frequent collisions between the heavy particles at atmospheric pressure [39].The plasma gas temperature is approximately equal to the rotational temperature of N2.By comparing the measured spectrum of N2(C3∏u→B3∏g,0 -2)with the simulated spectrum provided by the Specair software,the rotational temperature of N2can be obtained [40,41],as shown in figure 7.It is noted that the rotational temperature of N2at different locations is almost the same.This indicates that the plasma gas temperature is nearly constant in the discharge space.

    3.2.Comparison of discharges under different discharge gaps in static air

    Figure 7.The measured and simulated OES of N2(C3∏u→B3∏g,0 -2 )with a 5 mm discharge gap and a 1000 Hz PRF.(1)Location 1,(2) location 2,(3) location 3.

    In order to verify the effects of the transverse electric field on the discharge,the discharges under different discharge gaps are studied.Figure 8 shows side-view discharge images of different electrode structures when the discharge gaps are fixed at 4,5,and 6 mm.The PRF is 1000 Hz,and the corresponding exposure time is 1/1000 s.It can be seen that the diffuse discharge area in the misaligned region increases with the decrease of the discharge gaps.The lengths of the diffuse discharge region are integrated directly into the captions of figure 8.Based on the above analysis,the diffuse discharge area depends on the characteristics of the movement of the space charges during the same time interval.Therefore,the moving velocity of the micro-discharge channel is calculated by u=L/Δt,where L is the length of the path of the movement obtained by the images of the evolution of the discharge,and the time interval Δt between adjacent discharge images is fixed at 20 T.The time interval 1 T is equal to 1 ms.Accordingly,the calculated moving velocity of the micro-discharge channel is the time-average velocity per 20 discharge cycles.The initial positions of the tracked microdischarge channel under different discharge gaps are invariant,starting at the electrode edge.Since it is known from figure 5 that the micro-discharge channel of electrode geometry (a) is stationary,and the characteristics of the movement of the micro-discharge channel in the electrode geometries (b) and (c) are the same,only the moving velocities of the the micro-discharge channel for the electrode geometry(c)are calculated,as shown in figure 9.It is shown that the moving velocity of the micro-discharge channel increases with the decrease of discharge gaps.When the applied voltage is kept consistent,the magnitude of the electric field depends on the discharge gap.A higher applied electric field can generate many extra micro-discharge channels [42],and the magnitude of the transverse electric field component in the misaligned region increases in a smaller discharge gap.Consequently,the effect of the transverse electric field force on the space charges in the micro-discharge channel is enhanced,which leads to a larger moving velocity of the space charges.The result is that the space charges have a longer transverse migration distance.The result is similar to the previous work where the length and the propagation speed of filaments increased with VP[43].A more uniform space charge distribution can be formed in a larger region,and a stronger longitudinal electric field component can maintain the discharge over a longer distance.It is also proved that the diffuse discharge area in the misaligned region lies on the moving velocity of the space charges.

    Figure 8.Side-view images of the discharges of different electrode structures under different discharge gaps and a PRF of 1000 Hz.The exposure time is 1/1000 s.(1) Electrode geometry (a),(2) electrode geometry (b),(3) electrode geometry (c).

    Figure 9.The variation of the moving velocity of the microdischarge channel for the electrode geometry (c) under different discharge gaps and a PRF of 1000 Hz.

    At the same time,it is summarized in figure 9 that the moving velocity of the micro-discharge channel decreases with time.The transverse electric field only acts on the space charges in the misaligned region,which are repelled by the space charges from the aligned region.Thus their movement will be hindered.When the space charges accumulate continuously in the aligned region,the obstacles to the space charges in the misaligned region increase.Therefore,the moving velocity shows a downward trend.

    3.3.Comparison of discharges under different PRFs in static air

    Figure 10.Side-view images of the discharges of different electrode structures with a 5 mm discharge gap under different PRFs.(1)Electrode geometry (a),(2) electrode geometry (b),(3) electrode geometry (c).

    The effects of the transverse electric field on discharge can also be illustrated by comparing the discharges under different PRFs.Figure 10 shows side-view images of the discharge of different electrode structures with a 5 mm discharge gap.The values of the PRF are selected as,i.e.100,600,and 1000 Hz,and the exposure time corresponds to 1/100,1/640,and 1/1000 s,respectively.It can be observed that the diffuse discharge area of the misaligned region at a low PRF is larger than that at a high PRF.The lengths of the diffuse discharge region are integrated directly into the captions of figure 10.The PRF is an important parameter of a nanosecond pulsed DBD,and its effects on the discharge characteristics have been extensively investigated[38,44,45].The memory effect presented by the residual particles is the characteristic of repeated pulse breakdown,and it can affect the development of the next breakdown[26,44].On the one hand,the memory effect is more obvious,since the time interval of two consecutive pulses is shorter under a higher PRF[46].The space charges almost remain at the position where the previous discharge occurred during the shorter time,which causes the next discharge to appear at the location of the previous discharge.When the PRF is low,the space charges have enough time to spread during the longer time interval of two consecutive pulses.The distribution of the space charges is more uniform in a larger region where the space charges spread.On the other hand,based on the above analysis,the diffuse discharge area depends on the moving characteristics of the space charges during the same time interval.In the same way,only the moving velocities of the micro-discharge channel for the electrode geometry(c)under different PRFs are calculated,as presented in figure 11.It is concluded that the moving velocity of the micro-discharge channel increases with the decrease of PRFs.Under the experimental condition of the discharge gap and the applied voltage fixed consistently,the magnitude of the transverse electric field component is constant.Accordingly,the transverse electric field force on the space charges is unchanged.With such a low PRF,the space charges survived the previous discharge decay during the longer time interval of two consecutive pulses[26].Besides,the space charges have enough time to spread due to the density gradient and repulsion.When the next pulse arrives,the space charge density is lower.The result is that the hindering effect of the space charges is weakened under the same transverse electric field force.Accordingly,the effect of the transverse electric field on the movement of the space charges is not obvious at a high PRF.The result is that the space charges have a larger moving velocity at a low PRF.Taking these aspects into consideration,the space charges can be uniformly distributed over a larger region at a lower PRF.Also,it is indicated that the moving velocity of the space charges determines the length of the diffuse discharge region along the direction of the transverse electric field.Similarly,the moving velocity of the micro-discharge channel decreases with the discharge cycles.

    Figure 11.Variation of the moving velocity of the micro-discharge channel for the electrode geometry (c) with a 5 mm discharge gap under different PRFs.

    Figure 12.Top-view images of the discharges of different electrode structures under different airflow velocities.(1)Electrode geometry(a),(2) electrode geometry (b),(3) electrode geometry (c).The PRF is 1000 Hz,and the exposure time is 1/100 s.

    3.4.Discharge state transition in airflow

    Top-view images of the discharges of different electrode structures are shown in figure 12.The discharge gap is 4 mm,the PRF is 1000 Hz,and the exposure time is 1/100 s.The edge locations of the high-voltage and grounded electrodes are marked with red lines and green lines,respectively.The direction of the arrow indicates the direction of the airflow.From our previous works,it is known that the discharge state can be changed from a filamentary discharge to a diffuse discharge after introducing airflow [12,47].The generation of the diffuse discharge starts from the electrode edge at the airflow inlet due to the edge effect,which was observed in our previous experiment [12].Subsequently,the diffuse discharge at the electrode edge moves toward the downstream region under the action of the airflow.Meanwhile,airflow can also affect the distribution of the space charges.Ions with heavier mass leave in the volume,while electrons are accelerated by the external applied electric field and accumulate on the dielectric surface.This is because the drift velocity of electrons is larger than that of ions[9,48].Airflow first applies a drag force to the space ions because the airflow velocity in the boundary layer is almost zero,and the surface electrons move under the action of electrostatic coupling with the space ions [49].Therefore,the space charges in the upstream region move toward the downstream region,resulting in the decrease of the space charge density in the upstream region.Although the space charges in the downstream region also leave the region,it is supplemented by the space charges from the upstream region.The density remains basically unchanged.This has been reported in a previous work [50],which can verify that the PMT signal in the downstream region was earlier and stronger than that in the upstream region.Consequently,the diffuse discharge in the upstream region first transforms into a dense filamentary discharge,and the downstream region still manifests as a diffuse discharge.As the airflow increases further,the dense filamentary discharge extends continuously to the downstream region.The dense filamentary discharge is mainly ascribed to the relatively low space charge density.The lower space charge density is caused by the space charges moving toward the downstream region.Therefore,the length of the dense filamentary discharge in the upstream region can be regarded as the distance that the space charges move toward the downstream region within one discharge cycle.According to the formula u=L/t,where L is the length of the dense filamentary discharge region,and the time t is 1 ms.Tables 1 and 2 show the length of the dense filamentary discharge region and the moving velocity of the space charges for the different electrode structures,respectively.It can be seen that the moving velocity of the space charges is less than the corresponding airflow velocity.This phenomenon is consistent with our previous research using an alternatingvoltage source [49].This is mainly because the surface electrons need to overcome the shallow trap of the dielectric as well as the friction between the surface electrons and the dielectric plate.It is reasonable to believe that the reduction of the space charge density is the cause of the dense filamentary discharge.

    Table 1.Length of the dense filamentary discharge region for different electrode structures under different airflow velocities.

    Table 2.Moving velocity of the space charges for different electrode structures under different airflow velocities.

    As presented in figure 12,compared with the dense filamentary discharge region of the symmetric electrode geometry (a),the dense filamentary discharge region at the airflow inlet of the asymmetric electrode geometry (b) is larger.However,at a 5 m s-1airflow velocity,there is no dense filamentary discharge in the whole discharge space for the asymmetric electrode geometry (c),and the dense filamentary discharge region is smaller than that of symmetric electrode geometry.In other words,the electrode geometry(c)is conducive to improving the uniformity of the discharge.The effect of the transverse electric field on the movement of the space charges acts on the duration of the pulses,and the movement of the space charges induced by airflow occurs during the time interval of two consecutive pulses.When the airflow blows into the discharge gap in the same direction as the transverse electric field,the transverse electric field and the airflow cause the space charges to move in the same direction,which induces a larger displacement so that the space charges move toward the downstream region.The space charge density in the larger region is reduced,accordingly,and the dense filamentary discharge region at the airflow inlet is larger.However,for the asymmetric electrode geometry (c) the direction of the transverse electric field is opposite to the airflow,and the effect of the airflow on the movement of the space charges can be partly counteracted by the transverse electric field.The space charges do not easily move toward the downstream region,which results in the smaller dense filamentary discharge region at the airflow inlet.The above phenomenon can further illustrate the effect of the transverse electric field on the discharge by comparing the dense filamentary discharge region between the different electrode structures after introducing airflow.

    4.Conclusions

    The effects of the transverse electric field on a nanosecond pulsed DBD are studied in this paper.The discharge in the misaligned region is more diffuse than that in the central region for the asymmetric electrode geometry.The reason is that the space charges have transverse migration and are redistributed induced by the transverse electric field.Compared with the diffuse discharge at the edge for the symmetric electrode geometry,the diffuse discharge area is larger for the asymmetric electrode geometries.It can be explained that the micro-discharge channel is almost stationary due to the memory effect for the symmetric electrode geometry,however,the micro-discharge channel of asymmetric electrode geometries moves under the effect of the transverse electric field.In addition,the length of the diffuse discharge region along the direction of the transverse electric field increases with the decrease of discharge gaps and PRFs,which is consistent with the variation of the moving velocity of the micro-discharge channel.It is noted that the space charges can be uniformly distributed over a larger region under a faster moving velocity of the micro-discharge channel,which results in a larger area of diffuse discharge.

    When the airflow is consistent with the transverse electric field of the asymmetric electrode geometry is introduced into the discharge gap,the area of the dense filamentary discharge at airflow inlet is larger than that of the symmetric electrode geometry.Nevertheless,the area of dense filamentary discharge is smaller for the asymmetric electrode geometry and the transverse electric field is opposite to the airflow.The main reason for the observed phenomena is the redistribution of the space charges by the movement of the micro-discharge channel due to the transverse electric field.

    Acknowledgments

    This work is supported by National Natural Science Foundation of China (No.51437002).The authors would like to thank Professor Zhenfeng Ding for his technical assistance in the use of the high speed motion analyzer.

    猜你喜歡
    春生智慧
    曹春生作品
    占春生作品
    曹春生
    不認賬
    雜文選刊(2014年12期)2014-11-17 03:53:48
    有智慧的羊
    智慧派
    智慧決定成敗
    智慧往前沖,統(tǒng)計百分百(1)
    智慧往前沖,統(tǒng)計百分百(2)
    智慧往前沖,統(tǒng)計百分百(3)
    精品亚洲乱码少妇综合久久| 日韩一本色道免费dvd| 女人久久www免费人成看片| 日韩人妻高清精品专区| 国产亚洲最大av| 一本久久精品| 欧美精品高潮呻吟av久久| 久久久国产精品麻豆| 人成视频在线观看免费观看| 七月丁香在线播放| 久久 成人 亚洲| 日韩亚洲欧美综合| 国产成人精品一,二区| 久久精品人人爽人人爽视色| 久久精品国产自在天天线| 中文字幕最新亚洲高清| 青春草国产在线视频| 一区二区日韩欧美中文字幕 | 欧美日本中文国产一区发布| 国产一区二区三区av在线| √禁漫天堂资源中文www| 日韩免费高清中文字幕av| 免费看av在线观看网站| 国产精品99久久久久久久久| 午夜av观看不卡| 91久久精品国产一区二区三区| 老司机影院成人| 久久久a久久爽久久v久久| 精品酒店卫生间| 国产一区二区三区av在线| 制服人妻中文乱码| 免费黄色在线免费观看| 色5月婷婷丁香| 亚洲精品aⅴ在线观看| 国产成人免费观看mmmm| av电影中文网址| 成年av动漫网址| 人妻制服诱惑在线中文字幕| 亚洲伊人久久精品综合| 日韩欧美精品免费久久| 国产亚洲av片在线观看秒播厂| 性高湖久久久久久久久免费观看| 欧美日韩视频精品一区| 亚洲欧美日韩卡通动漫| av线在线观看网站| 夜夜看夜夜爽夜夜摸| 国产伦精品一区二区三区视频9| 午夜免费鲁丝| 99热这里只有精品一区| videosex国产| 美女国产视频在线观看| 午夜老司机福利剧场| 内地一区二区视频在线| 99精国产麻豆久久婷婷| 国产欧美另类精品又又久久亚洲欧美| 午夜福利影视在线免费观看| 亚洲三级黄色毛片| 久久久亚洲精品成人影院| 精品久久国产蜜桃| 国产乱来视频区| 午夜免费男女啪啪视频观看| 男人操女人黄网站| 王馨瑶露胸无遮挡在线观看| 精品99又大又爽又粗少妇毛片| 欧美激情极品国产一区二区三区 | 天堂8中文在线网| 国产69精品久久久久777片| 欧美老熟妇乱子伦牲交| 免费观看a级毛片全部| 日本av手机在线免费观看| 妹子高潮喷水视频| 国产白丝娇喘喷水9色精品| 日韩中字成人| 国语对白做爰xxxⅹ性视频网站| 十八禁高潮呻吟视频| 亚洲成人av在线免费| 国产一区二区在线观看av| 精品国产国语对白av| 丝瓜视频免费看黄片| videosex国产| 在线观看三级黄色| 少妇人妻精品综合一区二区| 伦精品一区二区三区| 欧美日韩在线观看h| 国产成人av激情在线播放 | 日本黄大片高清| 在线观看三级黄色| 下体分泌物呈黄色| 亚洲无线观看免费| 26uuu在线亚洲综合色| 高清不卡的av网站| 国产高清三级在线| 日韩一本色道免费dvd| 久久ye,这里只有精品| 一级黄片播放器| 色5月婷婷丁香| 国产熟女午夜一区二区三区 | 精品熟女少妇av免费看| 国产精品三级大全| 另类精品久久| 大陆偷拍与自拍| 日韩大片免费观看网站| 黄色视频在线播放观看不卡| 欧美+日韩+精品| 免费观看a级毛片全部| 美女福利国产在线| 男女无遮挡免费网站观看| av网站免费在线观看视频| a级毛片免费高清观看在线播放| 亚洲精品第二区| 伦精品一区二区三区| 精品国产乱码久久久久久小说| 好男人视频免费观看在线| 极品少妇高潮喷水抽搐| 22中文网久久字幕| 免费观看无遮挡的男女| 精品人妻一区二区三区麻豆| 国产永久视频网站| 亚洲精品久久午夜乱码| tube8黄色片| 亚洲av日韩在线播放| 精品久久久久久久久av| 晚上一个人看的免费电影| videosex国产| √禁漫天堂资源中文www| 国产视频内射| 欧美日韩综合久久久久久| 国产精品国产三级国产av玫瑰| 日本欧美国产在线视频| 午夜福利,免费看| 久久久久精品久久久久真实原创| 国产老妇伦熟女老妇高清| 国产午夜精品一二区理论片| 少妇被粗大猛烈的视频| 国产精品嫩草影院av在线观看| av专区在线播放| 国产黄频视频在线观看| 亚洲成人手机| 日韩欧美一区视频在线观看| 一区在线观看完整版| 男男h啪啪无遮挡| 老女人水多毛片| 国产伦精品一区二区三区视频9| 一区二区三区免费毛片| 成人午夜精彩视频在线观看| 久久精品国产自在天天线| av电影中文网址| 99久久精品一区二区三区| 在线看a的网站| 丰满乱子伦码专区| 国产精品嫩草影院av在线观看| 黑人巨大精品欧美一区二区蜜桃 | 婷婷成人精品国产| 在线精品无人区一区二区三| 成年美女黄网站色视频大全免费 | 亚洲精品成人av观看孕妇| 午夜影院在线不卡| 亚洲成色77777| 亚洲精品日韩av片在线观看| 久久影院123| 久久久久精品性色| 国产亚洲精品第一综合不卡 | 两个人的视频大全免费| www.色视频.com| 少妇猛男粗大的猛烈进出视频| 啦啦啦啦在线视频资源| 国产午夜精品一二区理论片| 91午夜精品亚洲一区二区三区| 国产精品欧美亚洲77777| 亚洲国产色片| 赤兔流量卡办理| 欧美人与性动交α欧美精品济南到 | 亚洲久久久国产精品| 免费播放大片免费观看视频在线观看| 在现免费观看毛片| 伊人久久精品亚洲午夜| 国产精品免费大片| 飞空精品影院首页| 久久久久久久久久久久大奶| 91久久精品国产一区二区成人| 免费黄频网站在线观看国产| 午夜av观看不卡| 插阴视频在线观看视频| 人妻一区二区av| 少妇人妻久久综合中文| 国产日韩欧美视频二区| 国产av码专区亚洲av| 多毛熟女@视频| 亚洲精品国产av蜜桃| 亚洲av日韩在线播放| 青青草视频在线视频观看| av黄色大香蕉| 91aial.com中文字幕在线观看| 午夜91福利影院| 国产精品国产av在线观看| 视频区图区小说| 国产精品嫩草影院av在线观看| 亚洲,欧美,日韩| 精品一区二区免费观看| 三级国产精品片| 久久久久久久久久人人人人人人| 丰满迷人的少妇在线观看| 国产在线一区二区三区精| 91久久精品电影网| 丝袜美足系列| 蜜桃在线观看..| 亚洲欧洲国产日韩| a级毛色黄片| 国产av码专区亚洲av| 中国国产av一级| 精品人妻一区二区三区麻豆| 视频在线观看一区二区三区| 一本大道久久a久久精品| 色5月婷婷丁香| 中国三级夫妇交换| 国产精品国产三级国产专区5o| 国产高清不卡午夜福利| 久久久久久久精品精品| 一区二区日韩欧美中文字幕 | 精品久久蜜臀av无| 丝袜在线中文字幕| freevideosex欧美| 国产老妇伦熟女老妇高清| 日本黄色日本黄色录像| 国产精品免费大片| 国产高清三级在线| 草草在线视频免费看| 熟女人妻精品中文字幕| 九九在线视频观看精品| 免费观看性生交大片5| 天天影视国产精品| 草草在线视频免费看| 亚洲精品成人av观看孕妇| 色网站视频免费| 黑人欧美特级aaaaaa片| 美女中出高潮动态图| 观看av在线不卡| 中文乱码字字幕精品一区二区三区| 伊人久久精品亚洲午夜| 在线观看国产h片| 日日啪夜夜爽| 欧美人与性动交α欧美精品济南到 | 久久热精品热| 亚洲精品国产av成人精品| 中国美白少妇内射xxxbb| 日韩三级伦理在线观看| 黑人巨大精品欧美一区二区蜜桃 | 久久午夜综合久久蜜桃| 一级毛片aaaaaa免费看小| 美女主播在线视频| 国产69精品久久久久777片| 国产免费福利视频在线观看| 性色avwww在线观看| 久久av网站| 亚洲av福利一区| 久久久久久久国产电影| 国产 一区精品| 久久久久网色| 插逼视频在线观看| xxxhd国产人妻xxx| 蜜桃在线观看..| 大码成人一级视频| 精品亚洲成a人片在线观看| 国产精品一国产av| 午夜激情福利司机影院| 成年美女黄网站色视频大全免费 | 亚洲欧美一区二区三区黑人 | 国产亚洲最大av| 又黄又爽又刺激的免费视频.| 精品久久久久久久久亚洲| 人妻少妇偷人精品九色| 精品熟女少妇av免费看| 精品久久久精品久久久| 五月开心婷婷网| 欧美少妇被猛烈插入视频| 亚洲少妇的诱惑av| 国产精品秋霞免费鲁丝片| 啦啦啦啦在线视频资源| 超色免费av| 国产亚洲精品久久久com| 一个人看视频在线观看www免费| 久久久国产一区二区| 高清在线视频一区二区三区| 视频区图区小说| 久久久精品区二区三区| 国产成人精品无人区| av卡一久久| 欧美+日韩+精品| 97在线视频观看| 国产成人一区二区在线| 色哟哟·www| av在线app专区| 日韩中字成人| 晚上一个人看的免费电影| 高清黄色对白视频在线免费看| 一区二区三区四区激情视频| 人人妻人人爽人人添夜夜欢视频| 丰满乱子伦码专区| 成人午夜精彩视频在线观看| 男女边吃奶边做爰视频| 精品一品国产午夜福利视频| 精品午夜福利在线看| 自线自在国产av| 日本av免费视频播放| 九色亚洲精品在线播放| av网站免费在线观看视频| 一本久久精品| 国产熟女欧美一区二区| 一区二区三区免费毛片| 国产男女超爽视频在线观看| 大片免费播放器 马上看| 另类亚洲欧美激情| 激情五月婷婷亚洲| 亚洲人成网站在线观看播放| a级毛片免费高清观看在线播放| 免费少妇av软件| 人人妻人人澡人人看| 久久精品久久久久久久性| 亚洲激情五月婷婷啪啪| 国产 精品1| 欧美日韩一区二区视频在线观看视频在线| 一区二区三区乱码不卡18| 欧美激情极品国产一区二区三区 | 另类精品久久| 日本av免费视频播放| 99热全是精品| 久久久午夜欧美精品| 五月天丁香电影| 国产精品久久久久久久电影| 美女cb高潮喷水在线观看| 久久精品人人爽人人爽视色| 大片免费播放器 马上看| 黄色一级大片看看| av在线播放精品| 99国产精品免费福利视频| 亚洲欧洲精品一区二区精品久久久 | 国产男女超爽视频在线观看| videossex国产| 国产在线免费精品| 在线观看免费高清a一片| 日本黄色日本黄色录像| 国精品久久久久久国模美| 亚洲在久久综合| 男女边吃奶边做爰视频| 久久毛片免费看一区二区三区| 国产男人的电影天堂91| 美女中出高潮动态图| av又黄又爽大尺度在线免费看| 伊人久久国产一区二区| 精品久久久久久电影网| 91成人精品电影| 美女国产视频在线观看| av在线老鸭窝| 国产男人的电影天堂91| 亚洲丝袜综合中文字幕| 激情五月婷婷亚洲| 日韩成人伦理影院| 人人妻人人添人人爽欧美一区卜| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 高清午夜精品一区二区三区| 婷婷色综合www| 九九在线视频观看精品| 大香蕉97超碰在线| 永久网站在线| 简卡轻食公司| 婷婷色麻豆天堂久久| 国产精品99久久99久久久不卡 | 久久人人爽人人片av| 久久热精品热| 国产日韩欧美视频二区| 欧美日韩在线观看h| 妹子高潮喷水视频| 欧美日韩精品成人综合77777| 亚洲成人一二三区av| 久久99精品国语久久久| 啦啦啦啦在线视频资源| 精品久久久久久久久亚洲| 特大巨黑吊av在线直播| 一区二区av电影网| 精品亚洲成国产av| av免费在线看不卡| av在线app专区| av免费在线看不卡| 成人综合一区亚洲| 久久久精品94久久精品| 国产毛片在线视频| 天美传媒精品一区二区| 性高湖久久久久久久久免费观看| 国模一区二区三区四区视频| 亚洲第一区二区三区不卡| 亚洲综合精品二区| 国产女主播在线喷水免费视频网站| 亚洲av免费高清在线观看| 亚洲欧美精品自产自拍| 久久久久精品性色| 国产精品国产三级专区第一集| 三级国产精品片| 成人手机av| 亚洲精品456在线播放app| 国产 一区精品| 中文字幕亚洲精品专区| 精品一区二区免费观看| 国产成人精品一,二区| 午夜激情久久久久久久| 18+在线观看网站| 青春草亚洲视频在线观看| 国产精品99久久久久久久久| 久久韩国三级中文字幕| 国产有黄有色有爽视频| 男人爽女人下面视频在线观看| 高清黄色对白视频在线免费看| 国产精品一区www在线观看| 日韩av在线免费看完整版不卡| 亚洲人与动物交配视频| 一区二区三区精品91| 性色av一级| 国产精品一国产av| 黑人巨大精品欧美一区二区蜜桃 | 夜夜骑夜夜射夜夜干| 亚洲欧美精品自产自拍| 少妇 在线观看| 日产精品乱码卡一卡2卡三| 亚洲成人一二三区av| 欧美bdsm另类| av黄色大香蕉| 免费久久久久久久精品成人欧美视频 | 国产精品蜜桃在线观看| 亚洲美女视频黄频| 国产成人精品一,二区| 色哟哟·www| 日韩成人伦理影院| 国产成人免费无遮挡视频| 日韩熟女老妇一区二区性免费视频| 欧美日韩国产mv在线观看视频| 制服丝袜香蕉在线| 狂野欧美激情性bbbbbb| 亚洲熟女精品中文字幕| 伊人久久国产一区二区| 日韩av免费高清视频| 精品人妻在线不人妻| 欧美性感艳星| videosex国产| 老司机亚洲免费影院| 久久久亚洲精品成人影院| 韩国高清视频一区二区三区| 亚洲欧洲国产日韩| 久久国产精品男人的天堂亚洲 | 精品人妻熟女av久视频| 国产精品无大码| 亚洲美女搞黄在线观看| 欧美精品国产亚洲| 熟女av电影| 日日爽夜夜爽网站| av.在线天堂| 日韩 亚洲 欧美在线| 中文欧美无线码| 免费黄网站久久成人精品| 亚洲国产精品专区欧美| 欧美一级a爱片免费观看看| 男女高潮啪啪啪动态图| 久久鲁丝午夜福利片| 热re99久久国产66热| 成年美女黄网站色视频大全免费 | 日韩伦理黄色片| 乱码一卡2卡4卡精品| 又粗又硬又长又爽又黄的视频| 热99久久久久精品小说推荐| 日日爽夜夜爽网站| 国产成人精品久久久久久| 亚洲欧美日韩另类电影网站| 国产欧美日韩综合在线一区二区| 少妇精品久久久久久久| 热99国产精品久久久久久7| 午夜福利网站1000一区二区三区| 欧美丝袜亚洲另类| 色婷婷av一区二区三区视频| 精品人妻在线不人妻| 伦理电影免费视频| 性高湖久久久久久久久免费观看| 色婷婷久久久亚洲欧美| 欧美精品一区二区免费开放| 十八禁高潮呻吟视频| 国产精品免费大片| 国产成人freesex在线| 日韩欧美精品免费久久| 亚洲精品乱码久久久久久按摩| 久久精品久久久久久久性| 国产精品秋霞免费鲁丝片| 中文字幕久久专区| 久久久久人妻精品一区果冻| 久久女婷五月综合色啪小说| 熟女人妻精品中文字幕| av福利片在线| 精品一区二区三卡| 天美传媒精品一区二区| 狠狠婷婷综合久久久久久88av| 欧美精品亚洲一区二区| 亚洲熟女精品中文字幕| tube8黄色片| 久久午夜综合久久蜜桃| 老司机影院毛片| 国内精品宾馆在线| 在线观看免费日韩欧美大片 | 精品久久久久久久久av| 亚洲精品美女久久av网站| 欧美日韩国产mv在线观看视频| 国产精品 国内视频| 久久久久久久大尺度免费视频| 一级毛片 在线播放| 热99久久久久精品小说推荐| 亚洲高清免费不卡视频| 精品人妻偷拍中文字幕| 飞空精品影院首页| 亚洲综合色网址| 日日摸夜夜添夜夜添av毛片| 丰满少妇做爰视频| 汤姆久久久久久久影院中文字幕| 亚洲国产精品国产精品| 国产亚洲最大av| 亚洲经典国产精华液单| 亚洲人与动物交配视频| 一区二区av电影网| 午夜视频国产福利| 卡戴珊不雅视频在线播放| 人成视频在线观看免费观看| 中文字幕亚洲精品专区| 国产男女内射视频| www.av在线官网国产| 女性被躁到高潮视频| 成人黄色视频免费在线看| 一级爰片在线观看| 日本av手机在线免费观看| 午夜福利在线观看免费完整高清在| 性色av一级| 成人国语在线视频| 国产成人午夜福利电影在线观看| 国产免费一区二区三区四区乱码| 亚洲av国产av综合av卡| 极品少妇高潮喷水抽搐| 丰满迷人的少妇在线观看| 精品国产露脸久久av麻豆| 麻豆乱淫一区二区| 成人国语在线视频| 99re6热这里在线精品视频| 亚洲四区av| 国产 一区精品| 久久国产精品男人的天堂亚洲 | 如何舔出高潮| 777米奇影视久久| 久久鲁丝午夜福利片| 男人爽女人下面视频在线观看| 婷婷色综合大香蕉| h视频一区二区三区| 亚洲色图综合在线观看| 亚洲av欧美aⅴ国产| 亚洲人成网站在线播| 免费av中文字幕在线| .国产精品久久| 一边摸一边做爽爽视频免费| 亚洲伊人久久精品综合| 中国美白少妇内射xxxbb| 亚洲欧美清纯卡通| 成人黄色视频免费在线看| 精品一区二区三卡| 91精品国产国语对白视频| a级毛片在线看网站| 国产深夜福利视频在线观看| av天堂久久9| 日本免费在线观看一区| 日韩一区二区三区影片| 免费不卡的大黄色大毛片视频在线观看| 亚洲,欧美,日韩| 五月天丁香电影| 久久精品久久久久久噜噜老黄| 两个人免费观看高清视频| 一区二区三区免费毛片| 久久99蜜桃精品久久| 欧美一级a爱片免费观看看| 免费观看无遮挡的男女| 91午夜精品亚洲一区二区三区| 国产乱人偷精品视频| 国产精品女同一区二区软件| 狠狠婷婷综合久久久久久88av| 亚洲精品456在线播放app| 国产男女内射视频| 免费高清在线观看日韩| 国产乱来视频区| 精品午夜福利在线看| 国产成人精品福利久久| 亚洲精品第二区| 国产精品国产三级国产av玫瑰| 国产成人freesex在线| 国产精品99久久99久久久不卡 | 久久久国产欧美日韩av| 欧美性感艳星| 午夜福利网站1000一区二区三区| 国产深夜福利视频在线观看| 欧美性感艳星| 国产女主播在线喷水免费视频网站| 亚洲不卡免费看| 美女xxoo啪啪120秒动态图| 我的老师免费观看完整版| 黑人巨大精品欧美一区二区蜜桃 | 妹子高潮喷水视频| 亚洲内射少妇av| 99精国产麻豆久久婷婷| 嘟嘟电影网在线观看| 午夜久久久在线观看| 亚洲成色77777| 免费看av在线观看网站| 午夜久久久在线观看| 国产av精品麻豆| 如何舔出高潮| 成人影院久久| 日本-黄色视频高清免费观看|