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

    The dynamics of a nanosecond gas discharge development with an extended slot cathode in argon

    2020-12-02 08:11:16ASHURBEKOVIMINOVSHAKHSINOVZAKARYAEVAandRABADANOV
    Plasma Science and Technology 2020年12期

    N A ASHURBEKOV , K O IMINOV, G S SHAKHSINOV,M Z ZAKARYAEVA,2 and K M RABADANOV

    1 Dagestan State University, Faculty of Physics, Makhachkala 367000, Russia

    2 Institute of Physics, DFRC of RAS, Makhachkala 367015, Russia

    Abstract

    Keywords: nanosecond gas discharge, numerical simulation, argon, excited atoms

    1.Introduction

    In recent years,plasma technologies based on low-energy ion flux sources have undergone significant development for various precision technological applications, for example, in atomic/molecular-layer surface etching of micro- and nanoelectronic materials [1-9].For such technological applications, various plasma systems for the formation and generation of ion fluxes with controlled ion energy are investigated.One of the effective sources of low-energy ion fluxes consists of a beam plasma having the spatial configuration of a‘plasma sheet’whose surface functions as a wide-aperture ion emitter in the transverse direction [4, 10, 11].For the formation of such ion flux sources, extended ribbon electron beams are required, particularly in terms of a high-voltage transverse nanosecond discharge generated by an extended slot cathode.In such a gas-discharge system,a ribbon electron beam is generated as a result of the emission of electrons from a cold cathode under the influence of bombardment by highenergy ions and atoms, as well as their acceleration in the region of the cathode potential drop(CPD)[12,13].By using ribbon electron beams behind the mesh anode in the region beside the collector, it has been shown that a beam plasma is created in the form of a thin ‘plasma sheet’ acting as the emitting surface of a wide-aperture beam formed by accelerated ions [11].Such a constructive use of a nanosecond discharge with a slot cathode requires a more detailed study of the dynamics of formation and the space-time distribution of the main parameters of the discharge.It is necessary to optimize the conditions for the generation of ribbon electron beam and to control the average electron energy in the‘plasma sheet’,which in its turn will make possible to control the energy of ion flows.

    Figure 1.Block diagram of experimental setup: (1) monochromator-spectrograph; (2) objective; (3) chamber; (4) A, C-anode, cathode;(5) high-speed Princeton Instruments PI-MAX3 ICCD camera; (6) Tektronix AFG 3022B synchronisation generator; (7) oscilloscope;(8) PC; (9) gas inlet and evacuation system; (10) HVPG.

    The present work is aimed at studying the dynamics of the formation of the main parameters of a high-voltage nanosecond discharge generated by a slot cathode in argon at gas pressures from 1 to 10 Torr at various values of an external electric field applied to the electrodes.The work studies the features of the spatiotemporal distribution of the main parameters of the discharge.

    2.Experimental part and results

    The experimental setup consists of a high voltage nanosecond pulse generator, synchronisation generator, discharge chamber, gas inflow-evacuation and pressure monitoring systems,as well as the diagnostic systems of discharge electrical and optical characteristics [13] (figure 1).The discharge chamber is assembled as a quartz tube with an internal diameter of 5 cm having a system of aluminium electrodes located inside at a distance of 0.6 cm apart.The anode consists of a flat plate 5 cm long,2 cm wide and 0.5 cm thick,while the cathode has the shape of a cylindrical rod measuring 5 cm in length and 1.2 cm in diameter,with a lengthwise cavity cut in the form of a slot with a width of 0.2 cm and depth of 0.6 cm [14].A Tektronix TDS 2024B oscilloscope was used for measuring the electrical characteristics of the discharge.In order to study the spatiotemporal dynamics of nanosecond discharge development, we used a PI-MAX3 Princeton Instruments high-speed camera equipped with an optical shutter having an adjustable exposure time of 2 ns and a time resolution between frames of up to 3 ns.The depth of field during highspeed photography of the discharge development was about 50 mm, when using a Nikon 200 mm lens with an f/32 fnumber.

    The PI-MAX3 system was synchronised with discharge current using a Tektronix AFG 3022B pulse generator.

    The radiation of the discharge was projected onto the input of a high-speed photorecorder PI-MAX3 ICCD by means of a photographic objective; next, a time-lapse recording was carried out at various time intervals between the triggering of the thyratron and the launch time of the ICCD camera.The exposure time for frame-by-frame photorecording was selected in the range from 5 to 20 ns,depending on the specific formation conditions for a nanosecond discharge.The signals from the PI-MAX3 camera were sent to a personal computer with subsequent fixation and processing of the optical images.For measurements, a pulseperiodic discharge mode with a frequency of 50 Hz was applied.On one hand, such a discharge repetition frequency greatly increases the intensification of the glow due to accumulation of an optical signal in the recording system.On the other hand, this frequency allows complete recombination of the plasma and the destruction of metastable particles in the interval between individual current pulses (for time≈2×10-2s) to take place.

    Figure 2.Optical patterns of radiation from a nanosecond discharge with a gap cathode in argon at U0=1 kV and p=1 Torr.(a) t=150 ns, (b) t=200 ns, (c) t=300 ns, (d) t=450 ns.

    Figure 3.Optical pattern of integral radiation (a) and pulses of combustion voltage and discharge current in argon (b) at U0=1 kV and p=1 Torr.

    Table 1.Estimated discharge parameters.

    In order to provide additional verification of the results of numerical simulation of excited atoms density, a nuanced investigation of panoramic discharge emission spectra depending on the amplitude of voltage pulses and gas pressure was carried out.For these purposes, an MS 7504i monochromator/spectrograph from SOL instruments Ltd equipped with a four-position turret for automated change of diffraction gratings (200, 1200, 1800, 2400 lines mm?1) was used.The recorded wavelengths are in the range of 267-1200 nm, the inverse linear dispersion of the device (for a 1200 lines mm?1grating) is 1 nm mm?1, and the spectral resolution is 0.037 nm (with 10 μm slit).The monochromator/spectrograph is equipped with a stepping motor for rotating diffraction gratings with the wavelength setting accuracy of about 0.05 nm, and a Hamamatsu HS 102H S 9840 CCD camera for recording the discharge emission spectrum.

    Figure 2 presents the characteristic optical patterns of the discharge at different moments of time from the beginning of a high-voltage voltage pulse (U0) being applied to the discharge gap.The dynamics of the spatiotemporal distribution of optical radiation in the discharge gap shows that the intensity of the optical radiation of the discharge initially increases, reaches its maximum value approximately 300 ns after the application of a voltage pulse to the electrodes, then subsequently decreases, at which point the discharge goes to the final stage.At gas pressures of about 1 Torr, the experiment demonstrates the characteristic spatial distribution of the optical radiation discharge.At all stages in the development of the nanosecond discharge, no discharge is observed inside the cathode cavity:the characteristic glow only appears in the gap between the electrodes.In this case, the radiation intensity is lower in the centre of the gap against the cathode cavity, leading to the formation of a dark region reaching almost to the surface of the anode.The narrow glow layer observed along the cathode surface,including in the cavity,is apparently due to the reflection of radiation from the polished cathode surface.Due to insufficient discharge radiation intensity at the initial stages, it was not possible to experimentally determine the propagation of the ionisation front from the anode surface or the overlap of the discharge gap.

    Figure 4.The cross section of the electrode system (a) and the geometry of the region of gas-discharge processes simulation (b).

    Figure 3 depicts the optical representation of the discharge integral radiation along with the corresponding discharge current-voltage characteristics.At a gas pressure of 5 Torr, it can be seen that the discharge penetrates into the cathode cavity, a dense plasma is formed in the cavity and a radiation intensity maximum appears at the exit of the cathode slot (figure 3(a)).

    The half-height voltage pulse duration is approximately 100 ns relative to a total discharge duration of 500 ns.The discharge current, which reaches its maximum value 200 ns after the start of electric gas breakdown, is equal to about 60 A (figure 3(b)).Using current-voltage characteristics, the approximate discharge parameters obtained at various stages in the development of the nanosecond discharge are as follows:JD-discharge current density,E/N-reduced electric field strength,udr-electron drift velocity and ne-electron density.These estimated values are given in table 1.

    The concentration of free electrons in the discharge was calculated on the basis of plasma conductivity.Since the mobility of ions is known to be two orders of magnitude lower than that of electrons, any ionic contribution to the electric current can be discounted.On this basis, for the density of the discharge current,the formulaJD=eneudr,can be written, where e is the electron charge,neis the concentration of free electrons andudris the electron drift velocity.TheJD=I/Scurrent density was calculated by the experimental values of the I current intensity and the S crosssectional area of the discharge.The electron drift velocity was determined from the plots on the basis of the E/N reduced electric field strength, where N is the concentration of gas atoms [15, 16].When determining the values of E in the plasma column, it was considered that, in a discharge with a hollow cathode,as well as in an anomalous discharge,almost all the voltage applied to the gap decreases in the region of the CPD.For estimates, the Up~UD/5 values of the voltage drop across the plasma column were taken [17].

    3.Numerical model for the dynamics of ionisation process

    Numerical modelling of the problem was carried out in the COMSOL Multiphysics Plasma Module.The ionization coefficients, excitation coefficients and transport coefficients of particles are the required inputs to the drift-diffusion equations, but they are dependent upon the electron energy distribution function (EDF), which is not known beforehand.The EDF was calculated as follows:the slow part of the EDF(the EDF of a slow group of electrons) was calculated as a Maxwell function, and a Gaussian function was used to find the fast part of the EDF(the EDF of a fast group of electrons),where the Gaussian function being 10% of the Maxwell function.Further, both functions were stitched.After that an interpolation table was formed for the EDF and this table was used as inputs of the Plasma Model.The geometry of the modelling area was selected on the basis of the actual dimensions of the discharge chamber; the electrode system used for the experimental studies was as described in the previous paragraph.The geometric model in the program was set to axisymmetric (figure 4).For the cylindrical coordinate system()φr z, ,(dependence onφis not considered), axial symmetry is specified relative to r=0.Further, COMSOL Multiphysics allows for reflecting the results obtained relative to the boundary of axial symmetry in the model.

    The aim of numerical modelling is to study the dynamics of the formation and development of a nanosecond gas with an extended slot cathode and flat anode, including the formation of the spatio-temporal structure of the discharge and the distribution of the main discharge parameters depending on the gas pressure and amplitude of the voltage pulses.The results of the modelling support a deeper understanding of the physics of the high-voltage nanosecond discharge operation,as well as helping to explain experimentally-observed features of this discharge in argon.

    In the field of modelling, a self-consistent system is solved, consisting of equations for the electric field, density and average momentum of electrons,along with equations for electron energy and transport of heavy particles.

    In order to numerically solve the system of equations,thene(t=0,x,y)=ni(t=0,x,y)=n0initial conditions were chosen withn0~108cm?3presenting the initial background plasma density determined by the residual processes in the gas-discharge system under the frequency-periodic plasma generation regime [18].

    The boundary conditions were selected as follows:

    (1)V=0 is the potential on the walls of the cathode;

    (2)V=U(t) is the potential on the anode wall, set from the experimentally measured value of the voltage pulse at the anode;

    (3) Boundary conditions on the dielectric walls of the discharge chamber:

    whereεris the relative dielectric constant;dsis the thickness of the dielectric layer;Vrefis the reference potential determined by the charge captured on the surface of the dielectric;σsis defined as follows:σs=n×Ji+n×Je,Jiis the normal component of the total density for the wall ion current,Jeis the normal component of the total density for the wall electron current,which was measured experimentally on the electrodes for these discharge conditions.

    In order to set the shape and parameters of the electric field potential,the U(t)voltage pulse measured experimentally on the electrodes for the corresponding discharge conditions was digitised.A ballast resistance was introduced into the calculation model that took the characteristics of the discharge circuit into account in accordance withU=U0?Ip×Rb,equation using the value selected by the coincidence of the calculated discharge current pulse with the experimentally measured voltage pulse.

    For numerical modelling, a triangular grid is built in the modelling area and the computational cycle is started with a given time step.The values of the electric field potential and the density of charged particles at each internal nodal point of the grid are calculated.In this case, elastic collisions of electrons with argon atoms,processes of electronic excitation,direct and stepwise ionisation of argon atoms, electronic deactivation of metastable states of atoms, as well as the Penning ionisation are taken into account.Here,the destruction of charged particles due to recombination in the volume and on the surface of the boundary material as a result of adhesion processes is also taken into account, along with secondary electron emission from the cathode surface having a coefficient of 0.02.

    In numerical modelling, the grid step varies within(1-9)×10-3cm, while the time step is taken to be equal to 10-12s[19].The iterative process used for refining the values of the calculated parameters stops when the specified accuracy is achieved.For the calculation results presented in this paper, the relative accuracy comprises 10-4.

    4.Results of numerical modelling for the dynamics of ionisation processes in argon

    A numerical modelling for the distribution of the electric field potential in the discharge gap and inside the cathode cavity was carried out during the development of ionisation.Under similar conditions, the formation and development dynamics of the ionisation wave front was calculated along with the spatiotemporal dynamics for the distribution of the charged particle and excited argon atom density, as well as for the energy of electrons along the discharge gap in a gas pressure range from 1 to 10 Torr and voltage pulse of 1 kV amplitude applied to the electrodes.The initial concentration of charged particles is taken to be equal to 108cm-3and uniformly distributed both in the discharge gap and inside the cathode cavity, corresponding in order of magnitude to the residual charge in pulse-periodic mode of gas ionisation.

    4.1.Calculation results at gas pressure of 1 Torr

    The results of numerical modelling for ionisation processes in a high-voltage nanosecond discharge generated by a gap cathode and obtained at a gas pressure of 1 Torr and a pulse amplitude of 1 kV are provided below.Figure 5 presents the dynamics of the spatiotemporal distribution of the density of charged particles in the discharge gap.As can be seen from the figure,the density of charged particles is redistributed due to electrons being pulled to the anode at the initial stages in the formation of electric gas breakdown.In this case, the electron concentration in the cathode cavity and around the cathode decreases, then gradually increases near the anode surface; in addition, the concentration of charged particles mainly increases in regions located closer to the edges of the flat anode.When a concentration of 5×108cm-3is reached in these regions after about 50 ns, two ionisation waves,directed towards the cathode,are formed(figure 5).Estimates indicate the propagation velocity of ionisation waves to be around 2×107cm s?1.After approximately 150 ns, a solid ionisation front forms at the anode surface between highconcentration regions and advances toward the cathode(figure 5(b)).As a result, a characteristic discharge structure is formed having an uneven distribution of plasma density and a maximum electron concentration of 3×1011cm-3(figure 5(d)).It should be noted that the plasma structure is formed at a distance of about 0.2 cm from the cathode surface;no penetration is observed into the cathode cavity with an approximate CPD region length of lc~0.2 cm (figure 5(d)).

    Figure 5.The results of numerical simulation of the dynamics of the spatio-temporal distribution of electron density in the discharge gap atU0=1 kV and p=1 Torr.(a) t=80 ns, (b) t=150 ns, (c) t=200 ns, (d) t=400 ns.

    Figure 6.The results of numerical simulation of the distribution of the electric field potential(a)and the distribution of the average electron energy (b) in the center of the discharge gap at U0=1 kV and p=1 Torr.

    The characteristic distribution features of the intensity of the optical discharge radiation in the experimentally-obtained optical patterns(figure 2)are in qualitative agreement with the results of numerical modelling (figure 5).However, as compared to the results of numerical modelling, the beginning of the plasma column in optical patterns turns out to be closer to the cathode surface; here, the lc~0.2 cm value can be considered as an upper estimate.At a gas pressure of p=1 Torr,a classical high-voltage hollow-cathode discharge pattern is realised.If the cathode cavity from the end is completely open, i.e.the h slot width in the cathode coincides with the size of the exit hole and no enough space is provided between the walls of the slot to accommodate the CPD areas(h<2lc),then no transmission of the electric field occurs in the cathode cavity.

    For an abnormal discharge in argon,an estimate of the lcvalue can be provided from a formula obtained by extrapolating the experimental data [20]

    At a gas pressure of p=1 Torr, lc=0.122 cm is obtained from formula (1).Since the cathode slot width h=0.2 cm under the conditions of this work,h<0.244 cm;consequently, no discharge formed by the cathodic potential drop will be observed inside the cathode cavity.This is also confirmed by the results of modelling the distribution of the electric field potential along the centre of the discharge gap(figure 6(a)), which predicts that the electric field at the entrance of the cathode cavity will first increase,then achieve its maximum value and finally drop.During the entire period of the formation of electric breakdown of the gas,the electric field hardly penetrates into the cathode cavity (figure 6(a)).

    Although some of the ions penetrate the cathode cavity,their density inside the cathode cavity is insufficient for a substantial change in the distribution of the electric field inside the cavity.These ions then reach the inner surface of the cathode cavity to cause secondary electron emission.Since there is practically no accelerating electric field inside the cathode cavity, secondary electrons diffuse from the cavity towards the exit for some time.At the exit, these electrons enter the cathode region and are accelerated to high energy.The mean free path of an electron with respect to inelastic processes comprises λ=1/(N·σ)=0.17 cm,where N=3.3×1016cm-3and σ=1.8×10-16cm2[21].Consequently, lc<λ and the electrons in collisionless mode pass the CPD region to gain ε=eUc~900 eV energy

    Figure 7.The results of numerical simulation of the distribution of electron density(a)and the density of excited atoms(b)in the center of the discharge gap at U0=1 kV and p=1 Torr.

    Figure 8.The results of numerical simulation of the dynamics of the spatio-temporal distribution of electron density in the discharge gap atU0=1 kV and p=5 Torr.(a) t=100 ns, (b) t=110 ns, (c) t=130 ns, (d) t=400 ns.

    (figure 6(a)).The results of numerical modelling confirm the high electron energy at the exit from the cathode cavity(figure 6(b)).An average electron energy of up to 80 eV designates the presence of a sufficient number of high-energy electrons in the discharge.Under the same conditions,the results of numerical modelling presented in [22] indicate the presence of electrons with energies of above 1000 eV in the discharge gap.The path length of accelerated electrons can be estimated according to the formulawhere εi=26 eV is the energy of formation of an argon ion pair.These estimates demonstrate that the accelerated electrons freely reach the anode surface at p=1 Torr.The low electron concentration in the discharge gap (figure 7(a)) is explained by the absence of the hollow cathode effect in the discharge,resulting in a dramatic increase in the electron concentration and discharge current.Under the considered conditions, the results of numerical modelling also showed the concentration of excited atoms to be comparable with the concentration of free electrons (figure 7(b)).

    This fact also confirms the presence of high-energy electrons in the discharge, producing more than ten acts of electronic excitation of atoms.The results suggest that highvoltage nanosecond discharge generated by a slot cathode at gas pressures of about p=1 Torr is applicable for the effective generation of a high-energy electron beam.

    4.2.Calculation results at gas pressure of 5 Tor r

    Figure 9.The results of numerical simulation of the distribution of the electric field potential(a)and the distribution of the average electron energy (b) in the center of the discharge gap at U0=1 kV and p=5 Torr.

    Figure 10.The results of numerical simulation of the distribution of electron density(a)and the density of excited atoms(b)in the center of the discharge gap at U0=1 kV and p=5 Torr.

    According to the numerical modelling, the dynamics of discharge formation at a gas pressure of 5 Torr remains the same as at p=1 Torr in the initial stages.Then, approximately 100 ns after the application of a voltage pulse to the electrodes from the centre of the flat anode, a third ionisation wave is formed, qualitatively altering the consequent discharge development dynamics along with the parameters of the formed plasma.An ionisation wave propagating at a speed of about 1.7×107cm s?1toward the cathode (figure 8(b))approaches the entrance to the cathode cavity, penetrating it when the electron concentration reaches ~1012cm-3(figure 8(c)).An estimate by formula (1) at a gas pressure of 5 Torr gives the value lc=0.024 cm; therefore, h>2×lc=0.048 cm and the cathode cavity is penetrated by the electric field (figure 9(a)),forming plasma behind it.A dense plasma with a concentration of about 1012cm?3shields the applied electric field, which is mainly localised in the cathode layer.Further,ions are accelerated in the cathode layer,penetrating into the cathode cavity to enhance the electric field.An electric field penetrating the cathode cavity also accelerates the initial electrons inside the cathode cavity with the subsequent generation of a quasi-neutral plasma accelerating the penetration of plasma from the discharge gap into the cavity of the cathode.As a result, a CPD region is formed in the cavity where ions are accelerated causing secondary electron emission.Electrons emitted from the side surfaces of the cavity pass through the CPD region without collisions (lc=0.024 cm<λ=0.030 cm)and are accelerated to energies of ε=eUc~650 eV (figure 9(a)).This is reflected in the reversed field from the opposite side, making oscillatory movements in the cathode cavity and repeatedly returning to the region of negative glow.

    As a result of such oscillations of electrons inside the cathode cavity, a multiple increase in gas ionisation occurs,while areas of negative glow at the side surfaces of the slot overlap and the hollow cathode effect is fully manifested.The presence of high energy electrons in the cathode cavity is also confirmed by the results of numerical modelling(figure 9(b)).Instead of oscillating between the walls of the cathode cavity,electrons accelerated inside the cavity in the cathode move towards the anode under the influence of an external electric field, resulting in an increase in their concentration at the cavity exit.At the same time, due to energy relaxation,accelerated electrons at the exit move to the energy region of the maximum values for excitation and ionisation crosssections.These two factors lead to the formation of a region having the maximum concentration of excited atoms and charged particles at the cathode cavity exit(figures 8(d),10(a)and (b)).The electron concentration in this plasma structure of ne~1015cm?3is an order of magnitude greater than that obtaining in the discharge gap (figures 8(d) and 10(a)).The concentration of excited atoms at the maximum reaches the value of 1014cm-3(figure 10(b)).

    The results of numerical modelling appear to be in good agreement with those obtained experimentally (figure 3(a)),where the bright elliptical glow region at the cathode cavity exit is visible in the optical picture of the discharge.The concentration of free electrons in the discharge gap obtained in the numerical experiment (ne~1014cm-3) is also consistent with the estimate of the electron concentration using conductivity (ne~8×1013cm-3) and given in table 1.

    Figure 11.The results of numerical simulation of the dynamics of the spatio-temporal distribution of electron density in the discharge gap at U0=1 kV and p=10 Torr.(a) t=50 ns, (b) t=85 ns, (c) t=150 ns, (d) t=400 ns.

    Figure 12.The results of numerical simulation of the distribution of the electric field potential(a)and the distribution of the average electron energy (b) in the center of the discharge gap at U0=1 kV and p=10 Torr.

    Some accelerated electrons (Λ~0.68 cm) reach the anode to provide ionisation both in the entire gap and on the surface of the anode,as well as compensating for ionic losses at the cathode.

    The current density of the anomalous glow discharge in argon can be estimated using the approximation formula[20]:

    An estimation of current density according to formula(2)for the conditions of table 1 at t=200 ns gives the value JAD~0.36 A cm?2,which is almost 50 times lower than the current density obtained from the experimental data for the considered discharge (table 1).Such a large difference in current density indicates that a transverse discharge generated by a slot cathode is fundamentally different from an anomalous discharge and is of independent interest for studying this type of discharge.The obtained results demonstrate the applicability of the studied discharge for obtaining high current densities at relatively low voltages applied to the electrodes.The formation of an extended bright plasma column is particularly appropriate as an extended optical radiation source for various applied problems.

    4.3.Calculation results at gas pressure of 10 Torr

    Figures 11-13 present the numerical modelling results of the development dynamics of ionisation processes,along with the spatial distribution of the main plasma discharge parameters obtained for a gas pressure of 10 Torr.The calculations indicate that the results obtained for a gas pressure of 10 Torr are qualitatively similar to those obtained for a gas pressure of 5 Torr.Indeed,even at a gas pressure of 10 Torr,two extreme ionisation waves are initially formed followed by a third ionisation wave in the centre of the discharge gap(figures 11(a) and (b)).

    Figure 13.The results of numerical simulation of the distribution of electron density(a)and the density of excited atoms(b)in the center of the discharge gap at U0=1 kV and p=10 Torr.

    At p=10 Torr,ionisation waves are formed earlier than those at p=5 Torr; here, the concentration of charged particles in the discharge gap grows faster,while the propagation speed of the ionisation wave decreases.When the density of charged particles reaches 1012cm-3, the plasma begins to penetrate into the cathode cavity and the plasma column contracts to the centre of the discharge gap (figure 11(c)).Next, the cathode layer forms; this is followed by the production of charged particles inside the cathode cavity.The electrons emitted from the side surfaces of the cavity in the cathode pass the CPD region without collisions(for this case,lc=0.012 cm <λ=0.015 cm) and are accelerated to energies of ε=eUc~400 eV(figure 12(a)).Further,accelerated electrons result in overlapping areas of negative glow at the side surfaces of the cathode slot and the hollow cathode effect is manifested.The generation of accelerated electrons in the cathode cavity is confirmed by the numerical modelling results (figure 12(b)).The average electron energy at the bottom of the cathode slot of around 32 eV (figure 12(b)) is almost three times less than that at p=5 Torr (figure 9(b)).

    The development of ionisation processes in the discharge also leads to the formation of a plasma structure at the exit of the cathode cavity with a maximum concentration of charged particles (figure 11(d)), promoting the self-limiting discharge current effect as discovered and studied in [23].

    The electron concentration in this plasma structure equal to ne~1015cm-3is an order of magnitude greater than in the discharge gap (figure 13(a)).Under similar conditions, the concentration of excited argon atoms behaves in a complex way.Initially, the atom electronic excitation processes proceed mainly in the discharge gap between the electrodes;however, after 120 ns, the density of excited atoms reaches a maximum of~5×1013cm-3at the entrance to the cathode cavity (figure 13(b)).Further, an increase is observed in the concentration of excited atoms inside the cathode cavity,especially at the cavity base; after about 200 ns, their characteristic spatial distribution is established with a sharp density peak of~1014cm-3at the cavity base.By this time, the density of excited atoms inside the cathode cavity and in the discharge gap drops to ~5×1012cm-3.

    It is interesting to note that the spatiotemporal distribution of the density of excited atoms appears to be appropriate for qualitative confirmation of the presence of high-energy electrons in the discharge,as well as for estimating the energy of these electrons and their density distribution in the discharge gap.At p=1 Torr, the maximum concentration of free electrons in the discharge comprises~1011cm-3, with the concentration of excited atoms reaching ~4×1011cm-3at the cathode cavity exit followed by a slight decrease in the gap between the electrodes (figures 6(a) and (b)).This indicates that the generation of high-energy electrons occurs at the cathode cavity exit and that the energy of these electrons is sufficient for carrying out multiple acts of atomic excitation and ionisation in the discharge gap.

    At a pressure of p=5 Torr,the maximum concentration of free electrons in the discharge is~1015cm-3, with the concentration of excited atoms equal to ~2×1014cm-3distributed over the entire length of the discharge with a small maximum at the exit from the cathode cavity(figures 8(a)and(b)).Therefore, in this case, accelerated electrons are generated both inside the cathode cavity and in the gap; nevertheless, the energy of these electrons is smaller as compared to p=1 Torr conditions.At a gas pressure of p=10 Torr,the concentration of free electrons turns out to be an order of magnitude higher than the concentration of excited atoms at the base of the cathode cavity and more than two orders of magnitude higher than in the rest of the discharge gap.Therefore, in this particular case, accelerated electrons are only generated inside the cathode cavity; this is due to that they have an insufficient quantity in the discharge gap and an energy of even less than that at p=5 Torr (figure 12(b)).Indeed,the mean free path of accelerated electrons comprises

    Λ~0.23 cm;attainment of the middle of the discharge gap is not observed for these electrons due to a relaxation of the electron energy in inelastic collisions with atoms.

    5.Conclusion

    The experimental study and numerical modelling of the formation and development dynamics of a high-voltage transverse nanosecond discharge generated by a slot cathode in argon indicate that, at gas pressures from 1 to 10 Torr,ionisation waves formed near the anode surface are propagated to the cathode at a speed of about 107cm s?1with the propagation velocity of ionisation wave decreasing with an increase in the gas pressure.At a pressure of p=1 Torr, a classical version of high-voltage discharge is realised with no penetration of plasma into the cathode cavity or hollow cathode effect observed to take place.The CPD region is formed between the cathode and the plasma column and electrons, which pass this region without collisions and are accelerated and freely reach the surface of the anode.An increase in gas pressure to 5 Torr leads to a penetration of plasma into the cathode cavity with the formation of a CPD region in the cavity,where ions are accelerated,subsequently causing secondary electron emission.The electrons emitted from the side surfaces of the cavity pass through the CPD region without collisions to oscillate inside the cathode cavity.This leads to overlapping areas of negative glow at the side surfaces of the cathode slot and a manifestation of the hollow cathode effect.

    At p=10 Torr, the modelling results qualitatively coincide with the results at p=5 Torr;however,in this case,no accelerated electrons are observed in the gap between the electrodes due to a relaxation in their energy, which occurs inside the cathode cavity as well as at the exit therefrom.In both cases, a single plasma structure is formed at the exit of the cathode cavity with the concentration of charged particles being an order of magnitude higher than in the rest of the gap,producing the self-limiting discharge current effect.

    The results of numerical modelling are in good agreement with those obtained over the course of experimental studies.

    Acknowledgments

    This work was supported by a grant from the Russian Foundation for Basic Research No.19-32-90179 and state assignment FZNZ-2020-0002.

    ORCID iDs

    N A ASHURBEKOV https://orcid.org/0000-0003-3045-9800

    av网站在线播放免费| 午夜福利一区二区在线看| 亚洲九九香蕉| 菩萨蛮人人尽说江南好唐韦庄| 色视频在线一区二区三区| 色老头精品视频在线观看| 99久久人妻综合| 高清视频免费观看一区二区| 黄片大片在线免费观看| 欧美久久黑人一区二区| 妹子高潮喷水视频| 亚洲,欧美精品.| 99国产精品一区二区蜜桃av | 亚洲免费av在线视频| 看免费av毛片| 超碰97精品在线观看| 久久天堂一区二区三区四区| 亚洲精华国产精华精| 成年女人毛片免费观看观看9 | 9热在线视频观看99| 一个人免费看片子| 日本a在线网址| 国产成人欧美在线观看 | 999久久久精品免费观看国产| 日韩精品免费视频一区二区三区| 老司机靠b影院| 久久毛片免费看一区二区三区| 亚洲一码二码三码区别大吗| 国产亚洲午夜精品一区二区久久| 久久国产精品影院| 免费在线观看完整版高清| 久久精品91无色码中文字幕| 丰满迷人的少妇在线观看| 黄色成人免费大全| 宅男免费午夜| 夜夜爽天天搞| 欧美大码av| av免费在线观看网站| 日本欧美视频一区| 电影成人av| 亚洲精品国产一区二区精华液| 亚洲专区字幕在线| 亚洲国产看品久久| 一边摸一边抽搐一进一小说 | 一边摸一边做爽爽视频免费| 少妇被粗大的猛进出69影院| 色精品久久人妻99蜜桃| 国产精品欧美亚洲77777| www日本在线高清视频| 无遮挡黄片免费观看| 99国产综合亚洲精品| 极品人妻少妇av视频| 亚洲中文日韩欧美视频| 在线观看人妻少妇| 人人妻人人澡人人爽人人夜夜| 中文字幕人妻熟女乱码| 午夜日韩欧美国产| 少妇精品久久久久久久| 中文字幕另类日韩欧美亚洲嫩草| 日本vs欧美在线观看视频| 最近最新中文字幕大全电影3 | 99国产极品粉嫩在线观看| av在线播放免费不卡| 美女高潮喷水抽搐中文字幕| 亚洲av电影在线进入| 免费看十八禁软件| 国产精品美女特级片免费视频播放器 | 亚洲精品美女久久av网站| 欧美黑人精品巨大| 国产一区二区三区在线臀色熟女 | 在线观看免费视频网站a站| 欧美黑人欧美精品刺激| 午夜精品久久久久久毛片777| 我的亚洲天堂| 99国产极品粉嫩在线观看| 亚洲九九香蕉| av网站在线播放免费| 手机成人av网站| 国产精品香港三级国产av潘金莲| a级毛片在线看网站| 国产精品欧美亚洲77777| 国产精品亚洲一级av第二区| 深夜精品福利| 一区二区av电影网| 欧美 亚洲 国产 日韩一| 在线观看一区二区三区激情| 69av精品久久久久久 | 日本黄色日本黄色录像| 亚洲成a人片在线一区二区| 91九色精品人成在线观看| 在线观看人妻少妇| 一级,二级,三级黄色视频| 成人影院久久| 国产精品国产av在线观看| 欧美激情高清一区二区三区| 亚洲欧洲日产国产| 国产一区二区三区视频了| 欧美激情极品国产一区二区三区| 国产单亲对白刺激| 成年版毛片免费区| www.999成人在线观看| 欧美日韩视频精品一区| 9191精品国产免费久久| avwww免费| av网站免费在线观看视频| 国产高清videossex| 亚洲三区欧美一区| av有码第一页| 亚洲免费av在线视频| 欧美性长视频在线观看| 欧美变态另类bdsm刘玥| 日本一区二区免费在线视频| av不卡在线播放| 国产不卡一卡二| 一级,二级,三级黄色视频| 久久影院123| 久久久久久免费高清国产稀缺| 啦啦啦 在线观看视频| 99热国产这里只有精品6| 亚洲av成人一区二区三| 国产主播在线观看一区二区| 制服诱惑二区| 丝袜美足系列| 久久国产精品男人的天堂亚洲| 少妇裸体淫交视频免费看高清 | 蜜桃国产av成人99| 日韩欧美三级三区| 国产区一区二久久| 亚洲欧美日韩另类电影网站| 亚洲中文av在线| aaaaa片日本免费| 成人精品一区二区免费| 丁香六月天网| 最新的欧美精品一区二区| xxxhd国产人妻xxx| 午夜福利视频在线观看免费| 午夜福利在线观看吧| 欧美大码av| 999精品在线视频| videos熟女内射| 久久影院123| 国产成人免费无遮挡视频| 日韩一区二区三区影片| a级片在线免费高清观看视频| 久久精品成人免费网站| 91大片在线观看| 动漫黄色视频在线观看| 80岁老熟妇乱子伦牲交| 国产黄色免费在线视频| 黑人猛操日本美女一级片| 欧美精品啪啪一区二区三区| 新久久久久国产一级毛片| 国产欧美日韩精品亚洲av| 肉色欧美久久久久久久蜜桃| 夜夜爽天天搞| 日本撒尿小便嘘嘘汇集6| av天堂久久9| 精品久久久精品久久久| 美女视频免费永久观看网站| 两个人免费观看高清视频| 精品国产一区二区三区久久久樱花| 欧美日韩视频精品一区| 黄色丝袜av网址大全| 成年人午夜在线观看视频| 久久婷婷成人综合色麻豆| 国产不卡一卡二| 99国产综合亚洲精品| 中文字幕人妻熟女乱码| 我的亚洲天堂| 国产亚洲一区二区精品| 亚洲av电影在线进入| 亚洲综合色网址| 亚洲三区欧美一区| 国产一区二区三区在线臀色熟女 | 1024视频免费在线观看| 亚洲成人手机| av视频免费观看在线观看| 黑丝袜美女国产一区| 欧美在线一区亚洲| 日韩精品免费视频一区二区三区| 深夜精品福利| 69av精品久久久久久 | 麻豆成人av在线观看| 精品一区二区三卡| 日本a在线网址| 一本—道久久a久久精品蜜桃钙片| 真人做人爱边吃奶动态| 女人被躁到高潮嗷嗷叫费观| 国产欧美日韩一区二区三| av网站免费在线观看视频| 国产精品亚洲av一区麻豆| 真人做人爱边吃奶动态| 大型黄色视频在线免费观看| 十八禁网站网址无遮挡| 亚洲avbb在线观看| 国产精品影院久久| 久久精品91无色码中文字幕| 99国产精品99久久久久| 亚洲男人天堂网一区| 精品久久久久久电影网| 国产无遮挡羞羞视频在线观看| 欧美日韩亚洲综合一区二区三区_| 国产成人啪精品午夜网站| 999久久久国产精品视频| 亚洲精品中文字幕一二三四区 | 国产有黄有色有爽视频| 久久影院123| 男女午夜视频在线观看| 日本av免费视频播放| 欧美日韩国产mv在线观看视频| 欧美黑人精品巨大| 欧美黄色淫秽网站| 动漫黄色视频在线观看| 亚洲国产欧美一区二区综合| 国产精品电影一区二区三区 | cao死你这个sao货| 亚洲国产欧美网| 国产不卡av网站在线观看| 99久久精品国产亚洲精品| 在线观看66精品国产| 女性生殖器流出的白浆| 淫妇啪啪啪对白视频| 欧美av亚洲av综合av国产av| e午夜精品久久久久久久| 啦啦啦中文免费视频观看日本| 国产精品久久久av美女十八| 18禁黄网站禁片午夜丰满| 一二三四在线观看免费中文在| 成人免费观看视频高清| 淫妇啪啪啪对白视频| 久久国产精品男人的天堂亚洲| 午夜精品久久久久久毛片777| 欧美激情久久久久久爽电影 | 9热在线视频观看99| 国产亚洲欧美在线一区二区| 50天的宝宝边吃奶边哭怎么回事| 久久毛片免费看一区二区三区| 操美女的视频在线观看| 国产高清视频在线播放一区| 丝袜人妻中文字幕| 成年版毛片免费区| 亚洲欧美一区二区三区黑人| 亚洲av成人一区二区三| 大香蕉久久网| 高潮久久久久久久久久久不卡| 天天影视国产精品| 亚洲免费av在线视频| 婷婷成人精品国产| 免费人妻精品一区二区三区视频| av免费在线观看网站| 黑丝袜美女国产一区| 国产在线精品亚洲第一网站| 久久久欧美国产精品| 亚洲成国产人片在线观看| 亚洲精品国产区一区二| 最新的欧美精品一区二区| 又黄又粗又硬又大视频| 90打野战视频偷拍视频| 人人澡人人妻人| 精品视频人人做人人爽| 夜夜夜夜夜久久久久| 国产精品久久久久久精品古装| 精品福利永久在线观看| 女性被躁到高潮视频| av有码第一页| 亚洲色图综合在线观看| 99国产精品免费福利视频| 国产亚洲午夜精品一区二区久久| 精品一区二区三卡| 精品欧美一区二区三区在线| 精品亚洲乱码少妇综合久久| 18禁黄网站禁片午夜丰满| 日韩中文字幕欧美一区二区| 国产精品久久久久久人妻精品电影 | 国产精品久久久av美女十八| 国产97色在线日韩免费| 韩国精品一区二区三区| 久久久国产成人免费| 黑人巨大精品欧美一区二区mp4| 91av网站免费观看| 天天躁日日躁夜夜躁夜夜| 啦啦啦 在线观看视频| 国产人伦9x9x在线观看| 考比视频在线观看| 国产激情久久老熟女| 97在线人人人人妻| 日本vs欧美在线观看视频| 欧美日韩亚洲高清精品| 两性夫妻黄色片| 亚洲欧美精品综合一区二区三区| 亚洲av日韩在线播放| 黑人操中国人逼视频| 国产欧美亚洲国产| 黄色怎么调成土黄色| 人妻 亚洲 视频| 日韩大片免费观看网站| 男女午夜视频在线观看| e午夜精品久久久久久久| 亚洲全国av大片| 久久午夜亚洲精品久久| 精品国产亚洲在线| 久久性视频一级片| 国产精品免费视频内射| av网站免费在线观看视频| 亚洲自偷自拍图片 自拍| 性高湖久久久久久久久免费观看| 国产免费视频播放在线视频| 欧美日本中文国产一区发布| 色视频在线一区二区三区| 他把我摸到了高潮在线观看 | 亚洲熟妇熟女久久| 我的亚洲天堂| 99re6热这里在线精品视频| 欧美国产精品一级二级三级| 国产av一区二区精品久久| 大码成人一级视频| 欧美精品高潮呻吟av久久| 热99re8久久精品国产| 12—13女人毛片做爰片一| 最新的欧美精品一区二区| 老司机在亚洲福利影院| av超薄肉色丝袜交足视频| 色婷婷久久久亚洲欧美| 精品国产乱码久久久久久男人| 日日夜夜操网爽| 精品一品国产午夜福利视频| 久久精品亚洲av国产电影网| 我的亚洲天堂| 精品人妻1区二区| 性色av乱码一区二区三区2| 日韩中文字幕视频在线看片| 在线观看免费视频日本深夜| av有码第一页| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲五月色婷婷综合| 在线天堂中文资源库| 18禁裸乳无遮挡动漫免费视频| 黑人巨大精品欧美一区二区mp4| 国产在线免费精品| av又黄又爽大尺度在线免费看| 激情视频va一区二区三区| 国产一区二区三区在线臀色熟女 | 国产男女超爽视频在线观看| 精品第一国产精品| 宅男免费午夜| 淫妇啪啪啪对白视频| a级片在线免费高清观看视频| www.熟女人妻精品国产| 两个人免费观看高清视频| 美女扒开内裤让男人捅视频| 国产淫语在线视频| 久久久精品94久久精品| 热99久久久久精品小说推荐| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品久久成人aⅴ小说| av又黄又爽大尺度在线免费看| 久久久国产欧美日韩av| 自拍欧美九色日韩亚洲蝌蚪91| 国产日韩欧美在线精品| 午夜久久久在线观看| av欧美777| 免费久久久久久久精品成人欧美视频| 久热这里只有精品99| 国产极品粉嫩免费观看在线| 欧美+亚洲+日韩+国产| 一进一出抽搐动态| 97在线人人人人妻| av又黄又爽大尺度在线免费看| 国产精品麻豆人妻色哟哟久久| 国产单亲对白刺激| 久久精品熟女亚洲av麻豆精品| 欧美日韩亚洲国产一区二区在线观看 | 久久久久久人人人人人| 人人妻人人澡人人看| 自线自在国产av| 日本一区二区免费在线视频| 97在线人人人人妻| 老司机午夜福利在线观看视频 | 超色免费av| 一区二区三区激情视频| 别揉我奶头~嗯~啊~动态视频| 免费在线观看影片大全网站| 操出白浆在线播放| 国产男靠女视频免费网站| e午夜精品久久久久久久| 男女午夜视频在线观看| 欧美成人免费av一区二区三区 | 日本精品一区二区三区蜜桃| 国精品久久久久久国模美| 久久av网站| 人人妻人人添人人爽欧美一区卜| 欧美日本中文国产一区发布| 午夜日韩欧美国产| 国产免费福利视频在线观看| 欧美精品亚洲一区二区| 一区二区三区精品91| 国产无遮挡羞羞视频在线观看| 久久久国产成人免费| 三上悠亚av全集在线观看| 精品国产乱码久久久久久小说| 久久人人97超碰香蕉20202| 日韩一卡2卡3卡4卡2021年| 精品人妻1区二区| 一本色道久久久久久精品综合| 国产日韩欧美在线精品| av一本久久久久| 国产欧美日韩一区二区三| 久久热在线av| 国产三级黄色录像| netflix在线观看网站| 91成年电影在线观看| 黄色丝袜av网址大全| 免费在线观看完整版高清| 亚洲精品久久午夜乱码| 免费高清在线观看日韩| 免费在线观看视频国产中文字幕亚洲| 99国产精品一区二区蜜桃av | 欧美黄色片欧美黄色片| 日日摸夜夜添夜夜添小说| 国产一区二区三区在线臀色熟女 | 美女福利国产在线| 老司机影院毛片| 高清欧美精品videossex| 亚洲av第一区精品v没综合| 下体分泌物呈黄色| 黄片大片在线免费观看| 天天躁夜夜躁狠狠躁躁| 久久久久久亚洲精品国产蜜桃av| 美女午夜性视频免费| 国产精品欧美亚洲77777| 欧美老熟妇乱子伦牲交| 中文字幕av电影在线播放| 国产精品自产拍在线观看55亚洲 | 免费少妇av软件| 欧美国产精品一级二级三级| 亚洲精品美女久久久久99蜜臀| 精品高清国产在线一区| 免费观看人在逋| 欧美日韩精品网址| 男人操女人黄网站| 乱人伦中国视频| av又黄又爽大尺度在线免费看| 久热这里只有精品99| 国产精品一区二区免费欧美| 无人区码免费观看不卡 | 欧美国产精品一级二级三级| 男女床上黄色一级片免费看| 啦啦啦视频在线资源免费观看| 丝袜美腿诱惑在线| 激情在线观看视频在线高清 | 亚洲欧美日韩高清在线视频 | 97人妻天天添夜夜摸| 国产精品偷伦视频观看了| 在线亚洲精品国产二区图片欧美| 十分钟在线观看高清视频www| 精品一品国产午夜福利视频| 久久久久久久久免费视频了| 国产三级黄色录像| 国产精品熟女久久久久浪| 久久精品成人免费网站| 女人被躁到高潮嗷嗷叫费观| 亚洲精品乱久久久久久| kizo精华| 99国产精品99久久久久| 男女下面插进去视频免费观看| 亚洲成人免费电影在线观看| 黄色片一级片一级黄色片| 亚洲国产毛片av蜜桃av| 黄色丝袜av网址大全| 在线天堂中文资源库| 亚洲专区字幕在线| 亚洲av国产av综合av卡| 12—13女人毛片做爰片一| 国产区一区二久久| 久久天堂一区二区三区四区| 少妇猛男粗大的猛烈进出视频| 日韩欧美免费精品| 亚洲专区国产一区二区| 日本一区二区免费在线视频| 精品国产一区二区三区四区第35| 亚洲情色 制服丝袜| 欧美日韩中文字幕国产精品一区二区三区 | 精品亚洲成国产av| 午夜成年电影在线免费观看| 啦啦啦 在线观看视频| 国产亚洲精品一区二区www | 国产精品一区二区在线不卡| 久久久国产精品麻豆| 国产精品熟女久久久久浪| 大香蕉久久成人网| 男女午夜视频在线观看| 国产高清国产精品国产三级| 精品久久久久久久毛片微露脸| 曰老女人黄片| 欧美日韩视频精品一区| 丝袜喷水一区| 最新美女视频免费是黄的| 成人亚洲精品一区在线观看| 欧美日韩精品网址| 日韩有码中文字幕| 精品一区二区三卡| 久久青草综合色| 一夜夜www| 亚洲欧美激情在线| 日本撒尿小便嘘嘘汇集6| 两个人免费观看高清视频| 国精品久久久久久国模美| 亚洲av第一区精品v没综合| 精品国产亚洲在线| 一本色道久久久久久精品综合| 国产一卡二卡三卡精品| 精品卡一卡二卡四卡免费| 电影成人av| 在线亚洲精品国产二区图片欧美| 欧美日韩亚洲国产一区二区在线观看 | av在线播放免费不卡| 免费在线观看日本一区| 肉色欧美久久久久久久蜜桃| 亚洲精品在线美女| 国产精品亚洲av一区麻豆| 黄片小视频在线播放| 丰满少妇做爰视频| 黄片大片在线免费观看| 无遮挡黄片免费观看| 国产aⅴ精品一区二区三区波| 视频区欧美日本亚洲| 97人妻天天添夜夜摸| 国产极品粉嫩免费观看在线| 久久精品aⅴ一区二区三区四区| 日日爽夜夜爽网站| 母亲3免费完整高清在线观看| 999久久久国产精品视频| 国产激情久久老熟女| 国产三级黄色录像| 黄色成人免费大全| 老司机午夜十八禁免费视频| 国产免费av片在线观看野外av| 热re99久久国产66热| 久久精品91无色码中文字幕| 亚洲伊人色综图| 国产高清国产精品国产三级| 叶爱在线成人免费视频播放| 动漫黄色视频在线观看| 日韩中文字幕视频在线看片| 老汉色av国产亚洲站长工具| a级毛片黄视频| 女警被强在线播放| 国产免费av片在线观看野外av| 大片免费播放器 马上看| 国产精品 国内视频| 日韩欧美国产一区二区入口| 国产精品一区二区在线不卡| 久久人妻福利社区极品人妻图片| 国产一区二区三区视频了| 777米奇影视久久| 亚洲精品粉嫩美女一区| 国产精品国产高清国产av | 久久国产亚洲av麻豆专区| 欧美精品人与动牲交sv欧美| 乱人伦中国视频| 最新美女视频免费是黄的| 肉色欧美久久久久久久蜜桃| 亚洲国产欧美日韩在线播放| 国产精品久久电影中文字幕 | a级毛片在线看网站| 人人妻,人人澡人人爽秒播| 成人影院久久| 国产伦人伦偷精品视频| 12—13女人毛片做爰片一| 国产精品免费大片| 精品久久久久久久毛片微露脸| 国产日韩欧美视频二区| 中文字幕色久视频| 国产午夜精品久久久久久| 黄色视频在线播放观看不卡| 一区二区av电影网| 亚洲第一青青草原| 真人做人爱边吃奶动态| 18在线观看网站| 国产成人系列免费观看| 日日摸夜夜添夜夜添小说| 1024视频免费在线观看| 免费黄频网站在线观看国产| 亚洲欧美日韩高清在线视频 | 国产高清视频在线播放一区| 成人18禁在线播放| 国产av一区二区精品久久| 欧美乱码精品一区二区三区| 国产av又大| 男女无遮挡免费网站观看| 午夜日韩欧美国产| 国产色视频综合| 日韩精品免费视频一区二区三区| 亚洲国产欧美网| 日本av免费视频播放| 一区福利在线观看| 午夜日韩欧美国产| 香蕉国产在线看| 一区福利在线观看| 美国免费a级毛片| 后天国语完整版免费观看| 午夜福利欧美成人| 高潮久久久久久久久久久不卡| 啪啪无遮挡十八禁网站| 免费看a级黄色片| av国产精品久久久久影院| 母亲3免费完整高清在线观看| 新久久久久国产一级毛片| av国产精品久久久久影院| 亚洲av国产av综合av卡| 亚洲国产毛片av蜜桃av| 少妇猛男粗大的猛烈进出视频| 嫁个100分男人电影在线观看| 国产不卡一卡二| 搡老熟女国产l中国老女人|