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

    Electron sheaths near a positively biased plate subjected to a weak electron beam

    2023-03-15 00:54:16QiLIU劉琦ZianWEI衛(wèi)子安JinxiuMA馬錦秀KaiyangYI弋開陽WanTAO陶萬andZhengyuanLI李政元
    Plasma Science and Technology 2023年3期
    關(guān)鍵詞:錦秀劉琦開陽

    Qi LIU (劉琦), Zian WEI (衛(wèi)子安), Jinxiu MA (馬錦秀),*, Kaiyang YI (弋開陽),Wan TAO (陶萬) and Zhengyuan LI (李政元)

    1 School of Physical Sciences, University of Science and Technology of China, Hefei 230026, People's Republic of China

    2 Department of Plasma Physics and Fusion Engineering,University of Science and Technology of China,Hefei 230026, People's Republic of China

    3 Southwestern Institute of Physics, Chengdu 610041, People's Republic of China

    Abstract Electron sheaths have previously only been measured near a positively biased small electrode,in which a potential dip was often observed.In this paper,we present an experimental study on the electron sheath near a stainless steel plate in the presence of a weak electron beam.It is shown that the electron beam,though its density is much lower than that of the background plasma,will substantially alter the sheath structure,i.e.,it causes the disappearance of the potential dip when the beam energy just exceeds the ionization potential of the neutral gas but later enhances the dip for higher energies.It is also shown that proper biases on the plate and chamber wall are the key to the formation of the electron sheath and the dip.For a fixed plate bias but with different electron beam energy,the measured thickness of the ion-free Child-Langmuir sheath agrees well with that of the theoretical model.

    Keywords: electron sheath, potential dip, electron-beam-plasma system

    1.Introduction

    Sheaths are non-neutral regions in plasma boundaries connecting to a wall or an electrode and have been widely studied since the 1920s [1-3].The sheaths are usually classified as electron and ion sheaths according to potential distributions[4].The ion sheaths are more commonly encountered both in laboratories and in nature.However, electron sheaths usually occur near a small electrode that is biased positively with respect to the plasma potential [5-7].Under the same experimental conditions, the electron sheaths are often much thinner than the ion sheaths,which makes the measurement of the electron sheaths more difficult [8-12].

    Previous experiments have shown that the potential distribution of the electron sheath was not monotonic, but often accompanied by a potential dip between the sheath and bulk plasma [5, 13-15].The sheath with the potential dip was also known as a double sheath,which is a kind of double layer since it consists of two regions of opposite space charge.Studies of the double sheath have attracted much interest over the past decade.There were also some studies on the sheath structures being influenced by injecting extra directional particles into the boundary of the electrode.For example, Fujita et al studied the V-shaped potential distribution near a grid due to electron beam injection [16].Bailung et al studied the effect of an ion beam on the electron sheath near a small plate in an ion-beam-plasma system[14].However, it is not clear how the electron sheath could form near an electrode immersed in an electron-beam-plasma system and how the electron beam affects the sheath structure.

    Past observations of the electron sheath were only made near a sufficiently-small positively-biased electrode, which was usually covered with a ceramic coating on the back.The existence of the previously measured potential dip was once explained by an ‘ion pumping mechanism’ [13, 17, 18].Because of the ionization process within the sheath, the generated ions will be trapped in the dip, but the generated electrons will be accelerated toward the electrode.Thus, the accumulation of positive charge will prevent the formation of the stationary sheath.To explain the observed stationary dip,Forest and Hershkowitz[13]once measured the equipotential contours near the electrode and found a potential drop near the edge of the electrode because the back side of the electrode was electrically floated at a negative potential.Thus,they proposed the removal (pumping) of the trapped ions along the surface to the edge of the electrode, which is, in fact, a 2D process.However, in many low-pressure discharges, the collision mean-free-path for the ionization is much longer (about the size of the discharge chamber) than the sheath thickness.Thus, the ion pumping might not be applicable to the low pressure discharges.

    Recently, Baalrud et al proposed a new model of global nonambipolar flow to explain under what conditions a monotonic electron sheath, a monotonic ion sheath, or a nonmonotonic(double)sheath could form[5]for the simplest case of the plasma consisting of a wall and a planar auxiliary electrode.They found that the area of the auxiliary electrode is a key parameter in determining which kind of sheath could form near it.If the area ratio of the auxiliary electrode to the wall is greater than 1.7μ, whereμ= (and meand miare the electron and ion mass, respectively, a monotonic ion sheath will form near the auxiliary electrode regardless of its bias.If the area ratio is less than μ, a monotonic electron sheath can form when the auxiliary electrode is biased at least a few volts more positive than the plasma.If the area ratio is between μ and 1.7μ, a double sheath (potential dip) will form near the auxiliary electrode.The act of the potential dip is to reflect some random electron flux back to the bulk plasma in order to keep the global balance of the electron and ion currents to the wall and to the auxiliary electrode.This model was verified in hot filament and hollow cathode discharges in recent experiments[5, 6, 19, 20].

    Another difficulty in the study of the electron sheath is that the positive bias on the electron-collecting electrode will significantly alter the plasma parameters, and in particular,raise the plasma potential.The bias of the electrode is relative to grounded chamber walls, which means that the potential drop from the electrode to the bulk plasma is a variable.However, the effect of the potential drop was not taken into account in previous experiments.

    In this work, we experimentally investigate the electron sheath formation near a positively biased plate immersed in a weak electron-beam-plasma system.It will be shown that the presence of the electron beam,though weak,will appreciably modify the electron sheath structure.The experimental arrangement is described in section 2.The formation of the electron-beam-plasma system and the measured plasma parameters are given in section 3.The electron sheath structures are presented in section 4.The main conclusions are summarized in section 5.

    2.Experimental setup

    The experiment was carried out in a double plasma device,which is schematically shown in figure 1.The device consists of two identical cylindrical chambers of 500 mm length and 500 mm diameter insulated from each other.The base pressure is p=2×10-4Pa and the working pressure is p=2.5×10-2Pa with Ar gas injected.Two separation grids,SG1 and SG2, with 60% geometric transmittance, were placed in parallel between the chambers.The plasma-generating part(in the source chamber)consists of an anode mesh(AM) and a hot cathode composed of 12 tungsten wires (hot filament, HF) of 120 mm length and 0.2 mm diameter.The plasma was produced with the discharge between the AM and HF.The discharge voltage and current are 30 V and 48 mA,respectively.The bulk plasma parameters in each chamber were diagnosed with two identical stainless steel planar Langmuir probes (10 mm diameter and 0.5 mm thickness,located 140 mm to the separation grids).The typical plasma parameters in the experimental chamber are electron density ne=(1-2)×1014m-3and temperature Te=3-4 eV.A stainless steel circular plate (80 mm in diameter, positively biased)was placed in the center of the experimental chamber for the purpose of the electron sheath formation.An axially movable emissive probe was used to measure the plasma potential distribution in the upstream side of the plate.

    Figure 1.Schematic diagram of the experimental setup.AM: anode mesh, HF: hot filaments, PP: planar probe, EP: emissive probe.

    The diagnostics using Langmuir probes are typical invasive methods,which were previously applied to diagnose some sheath-related properties,for example,the thickness of a collisional pre-sheath[21]and the expansion of a sheath[22].The size of the cold Langmuir probes is usually about (or even larger than) the sheath thickness.Thus, the interference of the cold probes with the sheath results in significant measurement deviation from reality [23].Hot (emissive)probe [24] diagnosis has become the most common method for the sheath diagnosis.The fine structured emissive probes(with sizes much smaller than the dimension of the sheath)were widely used in the studies of both the ion and electron sheaths [5,14,25-28].The emissive probe used in this work consists of a probe tip made of tungsten wire 0.025 mm in diameter and 5 mm in length (positioned with its length parallel to the surface of the plate) and lead wires insulated and supported within glass capillary tubes of 0.32 mm in diameter.The details of the probe structure can be found elsewhere [25].The resistance of the emissive probe is 1.5 Ω(including lead wires and connections)and the heating current is 230-240 mA.The spatial resolution of the probe movement is 0.02 mm.The plasma potential is measured using the inflection-point method from the probe I-V trace in the zeroemission limit [24].

    Table 1.Dependence of the plasma parameters on the bias of the plate, Vpl.The source and experimental chamber walls as well as SG1 and SG2 are grounded.

    3.Formation of electron-beam-plasma system

    Previous experiments showed that an ion-beam-plasma system with a large proportion of beam ions can be easily obtained in the double-plasma device [26, 29].To keep the experiment as stable as possible, both SG2 and the experimental chamber were grounded.To form an electron-beamplasma system, a necessary condition is that the potential of the experimental region must be higher than that of the source region.Two ways were used to control the electron beam energy by adjusting the potential difference between the experimental and the source regions.

    First, we investigated the influence of the plate bias, Vpl,on the plasma parameters while keeping the source chamber and SG1 grounded.A positive bias Vpl=5-25 V with respect to the ground was applied to the plate.The electron density ne,temperature Teand the plasma potentials φSand φEin the bulk plasmas of the source and experimental regions, measured with the planar Langmuir probes, together with the calculated Debye lengths λD,are listed in table 1.It is shown that the plasma parameters in the source region remain almost unchanged.In the experimental region,the changes in ne,Te,and λDare very small, but the plasma potential φEincreases with the increase of Vpl.It should be noted that the increase of φEis not linear with Vpl.When Vplexceeds 30 V,φEis raised close to the bias of the plate.

    It can be seen from table 1 that the plasma potential is higher in the experimental region than that in the source region.Thus, the electrons gain energy traversing from the source region through SG1 and SG2 into the experimental region,forming an electron beam.Since the plasma is weakly collisional in this device (collision mean-free-path is comparable to the size of the chambers),the energy of the electron beam, Eb, is approximately equal to the potential difference between the experimental and source region, i.e.,Eb≈e(φE-φS).Figure 2 shows the dependence of the potential drops φE-φSand Vpl-φEon the plate bias Vpl.It is shown that the electron beam energy increases with Vpl.Furthermore, the bias on the plate is higher than the plasma potential φE.Therefore, the electron sheath forms near the plate and the potential drop across the sheath increases with the increase of Vpl.

    Second, we keep the plate bias fixed at Vpl=20 V and change the negative bias,VS,on the source chamber wall and SG1(SG1 is connected to the source chamber)while keeping the experimental chamber and SG2 grounded.In this case,the measured plasma parameters are listed in table 2.It can be seen that the changes in ne,Te,and λDin the source chamber are very small except for the VS=0 case, and the changes in the experimental chamber are not appreciable, especially when VS≤-10 V.With the decrease of VSfrom 0 to-80 V,the plasma potential in the experimental region,φE,shows an initial increase of about 4 V and then remains almost unchanged after VS≤-10 V.However, the plasma potential in the source region, φS, decreases with VS, but φSis always about 4 V higher than VSexcept for VS=0 case.Figure 3 shows the dependence of the potential drops φE-φSand Vpl-φEon the bias VS.It is shown that φE-φS(i.e., the electron beam energy) decreases almost linearly with the increase of VS.The potential drop Vpl-φEonly shows a very slight decrease when VSchanges from -80 to -20 V and then a small increase after VS>-20 V.

    From the above results it can be seen that the potential difference between the experimental and source regions accelerates the electrons to form a beam.The beam electrons ionize the neutral gas to create a background plasma.Thus,the electron-beam-plasma system establishes in the experimental region.However, since the plasma potential in the source region is higher than the bias on SG1, most of the electrons from the source region are reflected by the potential barrier of SG1.Only a few electrons enter into the experimental region.Therefore, the electron beam is weak.This is manifested by the current-voltage trace of a spherical probe of 8 mm diameter inserted into the experimental region (about 100 mm from the SG2)as shown in figure 4(a)for the case of VS=-40 V and Vpl=20 V.The main knee position(inflection point) is near the spherical probe bias Vsp~15 V,corresponding to the bulk plasma potential of the experimental region.There is a second knee at Vsp~-39 V as shown in the magnified portion of the I-V trace in the inset of figure 4(a), which corresponds to the contribution from the electron beam.The difference between the two knee points(~44 V) corresponds to the electron beam energy.The soobtained beam energy is close to the plasma potential difference between the experimental and source regions(~39 V)measured with the planar probes.Note that there exist slight differences in the measured results among the spherical and planar probes and that the plasma potentials in the source and experimental regions are measured not with a single planar probe.The collected current from the electron beam to the spherical probe can be estimated by the current difference between the two tangents near the second knee[30]as shown in the inset of figure 4(a).The inferred beam-to-background density ratio, nb/ne=(vte/vb)(Ib/Ie) (vb, vte, Ib, and Ieare the speeds and the probe collected currents of the beam and thermal electrons,respectively)versus the potential difference φE-φSis shown in figure 4(b).It shows only a slight increase in nb/newhen the beam energy increases from 20 to 80 eV.nb/neis only about ~0.2%for most cases.Though the concentration of the electron beam is rather low, it has an appreciable effect on the sheath structure,as will be shown in the next section.

    Figure 2.Variations in the potential drop from the plate to the bulk plasma Vpl-φE in the experimental region and the plasma potential difference φE-φS between the experimental and source regions with respect to the bias on the plate Vpl.Both chamber walls and SG1 and SG2 are grounded.

    Figure 3.Variations in the potential drop from the plate to the bulk plasma Vpl-φE in the experimental region and the plasma potential difference φE-φS between the experimental and source regions with respect to the bias on the source chamber wall VS.The plate bias is Vpl=20 V.

    Table 2.The dependence of the plasma parameters on the bias of the source chamber wall and SG1, VS.The experimental chamber and SG2 are grounded and the plate is biased at 20 V with respect to the ground.

    4.Potential distribution of the electron sheaths

    For the case when both the source and experimental chambers are grounded, the potential distributions with respect to the distance from the plate for different plate biases are shown in figure 5.It is shown that the electron sheath structures are nearly monotonic in this case.As the plate bias increases,both the sheath potential drop and the electron beam energy increase (see also figure 2 and table 1), and thus the electron sheath thickness increases.

    For a given bias on the plate, Vpl=20 V, and different negative biases on the source chamber, the potential distributions with respect to the distance from the plate are shown in figure 6.As shown in the figure, when VS=0, the electron beam energy is low (see figure 3 and table 2) and there exists an obvious potential dip structure near the plate.As VSdecreases to -20 V, the dip structure weakens and almost disappears.With a further decrease in VS, the dip structure reappears and gradually enhances.

    For convenience of discussion, we define a few characteristic quantities as schematically shown in figure 7.We define the position of the potential minimum as the dip position,xdip,and the corresponding potential difference from the bulk plasma to the dip as φdip.The position where the potential starts to deviate from the bulk plasma is defined as the sheath edge, xs.The inflection point, xinfl, at which the charge neutrality satisfies (d2φ/dx2=-ρ/?0=0, where ρ is the net charge density and ?0is the vacuum dielectric constant),is defined as the interface between the electron and ion sheaths.The region between xinfland the plate is a negative charge layer(electron sheath)and that between xsand xinflis a positive charge layer(ion sheath).In the case when there is no potential dip structure,the inflection point and the sheath edge should coincide, i.e., xinfl=xs.Because for this kind of discharge the ion temperature (close to room temperature) is much lower than the electron temperature, all ions are essentially reflected by the positive sheath potential with respect to the bulk plasma potential.Hence, the point in the region between the dip and the plate where the potential is equal to the bulk plasma potential, xCL,can be defined as the edge of the electron Child-Langmuir sheath since the ion density can be neglected in the region between x=xCLand x=0.

    Figure 4.(a)Current-voltage trace of the 8 mm diameter spherical probe.The inset is a magnified portion of the trace near the second knee.(b)Beam-to-background electron density ratio,nb/ne,inferred from the second and main knees of the I-V trace,versus the bias on the source chamber wall.

    Figure 5.Potential distributions near the positively biased plate in the upstream side for different bias voltages on the plate.x=0 is the position of the plate.The source and experimental chambers and SG1 and SG2 are grounded.

    Figure 6.Potential distributions near the positively biased plate(Vpl=20 V) in the upstream side for different bias voltages on the source chamber and SG1.x=0 is the position of the plate.The experimental chamber and SG2 are grounded.The lines are the polynomial fits to the experimental data.

    Using the analog to the electron-free Child-Langmuir law, the ion-free Child-Langmuir law can similarly be obtained.The relationship between the thickness of the Child-Langmuir sheath, wCL, and the sheath potential drop from the plate bias Vplto the Child-Langmuir sheath-edge potential (equal to the bulk plasma potential of the experimental region φE) can be written as [31]

    where λDis the electron Debye length evaluated at the sheath edge.Figure 8 presents a comparison of the thickness of the Child-Langmuir sheath, |xCL|, obtained from the measured potential profiles,with that evaluated from equation(1),wCL,using the bulk plasma parameters obtained in section 3, versus the potential difference φE-φS.Good agreement between the measured and the theoretically evaluated Child-Langmuir sheath thicknesses is observed.However, the experimental Child-Langmuir sheath thickness is systematically larger by about 1 mm than the theoretical thickness.

    Figure 7.Schematic drawing of the electron sheath with the potential dip showing the characteristic points.xs is the edge of the sheath with the bulk plasma, xinfl is the inflection point separating the ion and electron sheaths, xdip is the point of the potential minimum, xCL is the edge of the ion-free Child-Langmuir sheath, and φdip is the potential difference between the bulk plasma and the potential minimum.

    Figure 8.Comparison of the measured (|xCL|, bullets) and the theoretically evaluated (wCL, closed circles) ion-free Child-Langmuir sheath thickness versus the potential difference φE-φS.

    In the presence of the potential dip, the variations of the depth φdipand the width xCL-xsof the dip together with the whole thickness of the electron sheath |xinfl| with respect to the electron beam energy, ~e(φE-φS), are shown in figure 9.It is shown that with the increase of the beam energy,the dip depth suddenly drops to almost zero and then slowly increases and saturates at about 0.2 V.The dip width and the whole electron sheath thickness show similar variations.

    Figure 9.Variations in the depth φdip (a), the width xCL-xs (b) of the potential dip,and the whole thickness of the electron sheath|xinfl|(c) with respect to the potential difference φE-φS.

    To explain the experimental results, we make use of the balance between the loss and generation of the charged particles [5].For simplicity, we only consider the balance in the experimental chamber.Since the potential of the bulk plasma is higher than the grounded chamber wall and the bias on the planar plate is higher than the bulk plasma potential, the ions are mainly lost to the chamber wall and the background electrons are mainly lost to the plate because they are reflected by the corresponding potential barrier.The beam electrons, on the other hand, can be lost both to the chamber wall and to the plate because of their higher energy to penetrate the barrier except for the case when VS=0.The charged particle generation in the experimental chamber is due to the ionization of the neutral argon gas by both the thermal electrons and the electron beam.In the presence of the potential dip near the plate, the balance of the electrons and ions can be written in a simplified form

    where G represents the generation of the charged particles,Apland Aware the areas of the plate and chamber wall, respectively, nband vbare the density and velocity of the electron beam,is the mean speed of the background electrons,is the ion-acoustic speed,and α accounts for the correction of the ion loss flux due to the presence of the electron beam.The quasi-neutrality condition requires ni=ne+nb.From equations (2) and (3), and noting Apl?Awandwe obtain the depth of the dip

    Equation (4) shows that the presence of the electron beam deepens the potential dip.Note that when VS=0 or-5 V (the first two data in figure 9), the energy of the electron beam is lower than the ionization potential of the Ar gas(15.6 eV),so that no loss of the beam electrons due to the ionization occurs.In this case, the fraction of the beam electrons is relatively high and hence the dip is relatively deep.When VS=-10 V,the energy of the electron beam is Eb~17.9 eV,which is higher than the ionization potential, and hence the impact ionization by the beam electrons takes place,which is evidenced by the sudden increase in the background electron density (from 0.7×1013to 1.3×1013m-3) as observed in table 2.In this case, the beam electrons are substantially consumed and nb/neis substantially lower than the case when VS>-10 V.This explains the sudden decrease in φdiptill Eb~17.9 eV.The ionization process is further enhanced until VS=-20 V (the background electron density is further increased from 1.3×1013to 1.8×1013m-3),after which the ionization seems to be saturated (the background electron density is around 2×1013m-3,see table 2).This can also be seen in figure 4(b) where nb/neis about 0.2%.However, the beam energy increases in this process.Thus,the increase in vbcan explain the slow increase in φdipafter Eb>17.9 eV.

    5.Discussion and conclusion

    In this paper, the formation of the electron-beam-plasma system and its influence on the electron sheath near a positively biased plate were studied in the experimental chamber of a double-plasma device.It was shown that the large-area(about the cross section of the device) and low-energy electron beam can be introduced by means of adjusting the biases on the plate and on the source chamber wall.Though the concentration of the electron beam is small,it has a significant effect on the electron sheath.

    For the case when both the source and experimental chambers are grounded,raising the bias on the plate results in the increase both in the electron sheath potential drop and in the electron beam energy.In this case, the electron sheaths near the plate are nearly monotonic and the potential dip structure is not obvious except for in the case when the bias is rather low (~5 V) when the dip structure tends to appear.

    For a given bias on the plate,raising the negative bias VSon the source chamber wall results in a nearly linear increase in the electron beam energy but the sheath potential drop remains almost fixed.The presence of the electron beam enhances the dip structure, i.e., the increase in the beam energy causes the deepening and widening of the potential dip except for the initial rapid drop in the dip depth when the beam energy increases from below to above the ionization potential of the argon gas.The latter effect is attributed to the rapid loss of the beam electrons to the background electrons in the ionization process.Furthermore, the experimentally measured variation of the thickness of the ion-free Child-Langmuir sheath agrees well with that of the theoretical model.

    It was previously shown that the formation of the dip structure near an electron sheath is the result of the selfadjustment of the plasma to satisfy the global balance between the fluxes of the electron and ion losses[5,6].In the presence of the electron beam, the total electron loss flux increases.Thus, it is expected that a deeper potential dip forms self-consistently to repel more thermal electrons from being lost to the plate in order to maintain the global current balance.Therefore, our experimental results showing that the dip depth increases with the beam energy are in agreement with the theoretical expectation.

    The electron sheath exists widely both in laboratory plasmas and in applications, such as that near the Langmuir probe and the surface of a spacecraft.The existence of the electron beam significantly modifies the sheath properties.Thus, the investigation carried out in this work may have practical applications in interpreting the possible deformation of the collected currents of the Langmuir probe and the surface charging of the spacecraft when they are in the plasma containing the electron beam.In addition, it may also have potential applications in plasma-material processing using the electron beam.

    Acknowledgments

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

    猜你喜歡
    錦秀劉琦開陽
    巧思維切入,妙歸納總結(jié)
    ———一道平面向量考題的探究
    咕咕叫的肚皮
    初心引航,構(gòu)建“雙減”新樣態(tài)
    現(xiàn)實版《桃姐》:雇主一家合力贍養(yǎng)百歲保姆
    百姓生活(2018年10期)2018-11-05 06:12:22
    現(xiàn)實版《桃姐》:雇主一家合力贍養(yǎng)百歲保姆
    新傳奇(2018年47期)2018-08-09 09:41:54
    GLOBAL EXISTENCE OF CLASSICAL SOLUTIONS TO THE HYPERBOLIC GEOMETRY FLOW WITH TIME-DEPENDENT DISSIPATION?
    開陽十里畫廊
    你幫我養(yǎng)小,我為你養(yǎng)老
    Application of Class Activities in English Teaching
    仙非青梅君非馬
    飛魔幻A(2015年3期)2015-05-14 09:46:33
    中文字幕精品亚洲无线码一区| 一进一出抽搐动态| 一个人免费在线观看电影 | 国产成年人精品一区二区| 午夜免费观看网址| 男女下面进入的视频免费午夜| 美女免费视频网站| 中文字幕人成人乱码亚洲影| 久久人人精品亚洲av| 午夜影院日韩av| 亚洲全国av大片| 三级国产精品欧美在线观看 | av天堂在线播放| 成人国产一区最新在线观看| 日韩免费av在线播放| 日本黄色视频三级网站网址| 麻豆成人av在线观看| 亚洲专区国产一区二区| av福利片在线观看| 国产精品99久久99久久久不卡| 欧美不卡视频在线免费观看 | 最近最新中文字幕大全免费视频| 日本免费a在线| 日本免费a在线| 成人午夜高清在线视频| 狠狠狠狠99中文字幕| 欧美成人性av电影在线观看| 精品国产乱码久久久久久男人| 听说在线观看完整版免费高清| 亚洲无线在线观看| 国产精品日韩av在线免费观看| 亚洲av日韩精品久久久久久密| 国产一区二区激情短视频| 国产午夜精品久久久久久| 深夜精品福利| 国产91精品成人一区二区三区| 欧美黑人欧美精品刺激| 十八禁网站免费在线| 变态另类丝袜制服| 老汉色av国产亚洲站长工具| 日韩欧美三级三区| 亚洲激情在线av| 亚洲激情在线av| www.999成人在线观看| 国产精品九九99| a级毛片a级免费在线| 亚洲性夜色夜夜综合| 搡老妇女老女人老熟妇| 国产亚洲精品久久久久5区| 18禁美女被吸乳视频| 91成年电影在线观看| 久久久水蜜桃国产精品网| 国产欧美日韩一区二区三| 此物有八面人人有两片| 欧美日韩瑟瑟在线播放| 国产精品国产高清国产av| 亚洲av五月六月丁香网| 国产精品一区二区精品视频观看| 国产精品一区二区精品视频观看| 午夜两性在线视频| 老司机靠b影院| 91字幕亚洲| 人人妻,人人澡人人爽秒播| 亚洲成人精品中文字幕电影| 国产男靠女视频免费网站| 精品久久久久久,| 午夜a级毛片| 亚洲av美国av| 久久精品综合一区二区三区| 国产精品av久久久久免费| 老鸭窝网址在线观看| 久久午夜亚洲精品久久| 午夜老司机福利片| 免费看日本二区| 亚洲 欧美一区二区三区| 国产精品电影一区二区三区| 久久久久免费精品人妻一区二区| 日韩精品中文字幕看吧| 在线观看66精品国产| 国产一区二区三区视频了| 国产三级中文精品| 99re在线观看精品视频| 精品乱码久久久久久99久播| 岛国视频午夜一区免费看| 色综合站精品国产| 日本a在线网址| 免费看十八禁软件| 狂野欧美白嫩少妇大欣赏| 桃色一区二区三区在线观看| 99热这里只有精品一区 | 女人被狂操c到高潮| 在线播放国产精品三级| 午夜精品一区二区三区免费看| 亚洲国产日韩欧美精品在线观看 | 亚洲成a人片在线一区二区| 成熟少妇高潮喷水视频| 国产精品一区二区免费欧美| 亚洲中文日韩欧美视频| 国产精品98久久久久久宅男小说| 亚洲中文字幕一区二区三区有码在线看 | 亚洲av电影在线进入| 国产成人精品无人区| 日本一二三区视频观看| 久久这里只有精品19| 欧美极品一区二区三区四区| 最近最新中文字幕大全电影3| 在线国产一区二区在线| 亚洲av五月六月丁香网| 欧美一区二区精品小视频在线| 日韩成人在线观看一区二区三区| 午夜免费观看网址| 舔av片在线| 国产精品 国内视频| 在线观看免费午夜福利视频| 首页视频小说图片口味搜索| 级片在线观看| 成人一区二区视频在线观看| 欧美绝顶高潮抽搐喷水| 久久香蕉精品热| 最好的美女福利视频网| 一进一出好大好爽视频| 精品久久久久久久毛片微露脸| 99久久无色码亚洲精品果冻| 99热6这里只有精品| 国产精品一区二区免费欧美| 欧美不卡视频在线免费观看 | 国产精品亚洲av一区麻豆| 人妻丰满熟妇av一区二区三区| 午夜激情av网站| 亚洲精品一区av在线观看| 国产精品一及| 一夜夜www| 黄色片一级片一级黄色片| 中文在线观看免费www的网站 | 岛国在线免费视频观看| 免费在线观看黄色视频的| 天堂影院成人在线观看| 国产午夜精品论理片| 亚洲专区中文字幕在线| 国产精品久久久久久精品电影| 日本a在线网址| 在线看三级毛片| 久久久精品大字幕| 国产亚洲av高清不卡| 国产在线观看jvid| 国产精品爽爽va在线观看网站| 禁无遮挡网站| 特大巨黑吊av在线直播| 日韩欧美在线乱码| 黄色 视频免费看| 久久精品aⅴ一区二区三区四区| 国产精品自产拍在线观看55亚洲| 99在线视频只有这里精品首页| 国产人伦9x9x在线观看| www日本在线高清视频| 999久久久国产精品视频| 久久天堂一区二区三区四区| 欧美日韩亚洲综合一区二区三区_| 精品熟女少妇八av免费久了| 国产aⅴ精品一区二区三区波| 777久久人妻少妇嫩草av网站| 国产在线观看jvid| 亚洲狠狠婷婷综合久久图片| 在线观看午夜福利视频| 欧美丝袜亚洲另类 | 亚洲av电影不卡..在线观看| 亚洲成av人片在线播放无| 两个人的视频大全免费| 久久这里只有精品中国| 日本免费a在线| 国产麻豆成人av免费视频| 怎么达到女性高潮| 成人特级黄色片久久久久久久| 久久精品夜夜夜夜夜久久蜜豆 | 日本在线视频免费播放| 成人高潮视频无遮挡免费网站| 99国产极品粉嫩在线观看| 久久久水蜜桃国产精品网| 午夜福利视频1000在线观看| 成人18禁在线播放| 国产精品av视频在线免费观看| 亚洲国产欧洲综合997久久,| av超薄肉色丝袜交足视频| 精品久久蜜臀av无| 一本大道久久a久久精品| 大型黄色视频在线免费观看| 国产高清视频在线观看网站| 久久久久九九精品影院| 高清在线国产一区| 亚洲色图av天堂| av在线播放免费不卡| 色播亚洲综合网| 成熟少妇高潮喷水视频| 欧美性长视频在线观看| 在线视频色国产色| 99国产综合亚洲精品| 国产成人一区二区三区免费视频网站| 免费观看精品视频网站| 国产一区二区三区视频了| 免费电影在线观看免费观看| 久久天躁狠狠躁夜夜2o2o| 成人亚洲精品av一区二区| 亚洲va日本ⅴa欧美va伊人久久| 无遮挡黄片免费观看| 制服丝袜大香蕉在线| 夜夜夜夜夜久久久久| 亚洲av片天天在线观看| 男人舔女人的私密视频| av有码第一页| 亚洲精品久久国产高清桃花| 成人18禁在线播放| 在线观看66精品国产| 国产蜜桃级精品一区二区三区| 国产成人精品无人区| 亚洲精品国产一区二区精华液| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品自产拍在线观看55亚洲| 精品无人区乱码1区二区| 成人欧美大片| 悠悠久久av| 一二三四社区在线视频社区8| tocl精华| 国产激情欧美一区二区| 麻豆久久精品国产亚洲av| 波多野结衣高清作品| 午夜福利视频1000在线观看| 国产精品一区二区免费欧美| 国产一级毛片七仙女欲春2| 亚洲精品美女久久久久99蜜臀| 高清毛片免费观看视频网站| 国产片内射在线| 欧美成人午夜精品| 免费电影在线观看免费观看| 一本久久中文字幕| 午夜久久久久精精品| 一二三四在线观看免费中文在| 国产主播在线观看一区二区| 中亚洲国语对白在线视频| 国产av一区在线观看免费| 变态另类成人亚洲欧美熟女| 少妇被粗大的猛进出69影院| 麻豆av在线久日| 国内久久婷婷六月综合欲色啪| 亚洲一区中文字幕在线| 久久亚洲精品不卡| 亚洲国产欧美一区二区综合| 久久久久精品国产欧美久久久| 国产亚洲av高清不卡| 精品乱码久久久久久99久播| 后天国语完整版免费观看| 国产99久久九九免费精品| 人妻夜夜爽99麻豆av| 亚洲熟妇熟女久久| 国产成人一区二区三区免费视频网站| 欧美日本视频| 在线视频色国产色| 人成视频在线观看免费观看| √禁漫天堂资源中文www| 又紧又爽又黄一区二区| 免费在线观看完整版高清| 久久午夜亚洲精品久久| 国产成人精品久久二区二区免费| av天堂在线播放| 国产高清有码在线观看视频 | 精品一区二区三区视频在线观看免费| 亚洲国产欧洲综合997久久,| 女警被强在线播放| 9191精品国产免费久久| 五月伊人婷婷丁香| 午夜福利在线观看吧| 欧美大码av| 12—13女人毛片做爰片一| 亚洲无线在线观看| 国产单亲对白刺激| 老熟妇乱子伦视频在线观看| 香蕉av资源在线| 国产av又大| 国产亚洲精品综合一区在线观看 | 色综合欧美亚洲国产小说| 国语自产精品视频在线第100页| 精品国产美女av久久久久小说| 久久九九热精品免费| 成人精品一区二区免费| 午夜福利18| 一个人观看的视频www高清免费观看 | 亚洲欧美一区二区三区黑人| 成人亚洲精品av一区二区| 三级男女做爰猛烈吃奶摸视频| 五月玫瑰六月丁香| 国产成人av教育| 欧美性长视频在线观看| 久久精品国产亚洲av香蕉五月| 午夜福利视频1000在线观看| 黄色毛片三级朝国网站| 久久人人精品亚洲av| 亚洲国产精品成人综合色| 91老司机精品| 久久精品人妻少妇| 日韩欧美免费精品| 日本三级黄在线观看| 99久久久亚洲精品蜜臀av| or卡值多少钱| 99久久综合精品五月天人人| 三级男女做爰猛烈吃奶摸视频| 波多野结衣高清作品| 久久久久性生活片| 国产精品乱码一区二三区的特点| 久久婷婷人人爽人人干人人爱| 国产精华一区二区三区| 日日夜夜操网爽| 午夜影院日韩av| a级毛片在线看网站| 精品国内亚洲2022精品成人| 午夜免费成人在线视频| 亚洲成av人片在线播放无| 无遮挡黄片免费观看| 精品久久久久久久末码| 一级毛片女人18水好多| 午夜福利在线在线| 最近最新中文字幕大全免费视频| 中文字幕最新亚洲高清| 我的老师免费观看完整版| 亚洲九九香蕉| 高清在线国产一区| 男人舔女人的私密视频| 性色av乱码一区二区三区2| 在线播放国产精品三级| 成人av一区二区三区在线看| 亚洲九九香蕉| 一本精品99久久精品77| 嫩草影视91久久| 18禁美女被吸乳视频| 制服人妻中文乱码| 国产精品综合久久久久久久免费| 国产精品爽爽va在线观看网站| 欧美国产日韩亚洲一区| 两个人的视频大全免费| 男女之事视频高清在线观看| 久久午夜综合久久蜜桃| 婷婷精品国产亚洲av在线| 国产精品 国内视频| 成人一区二区视频在线观看| 亚洲欧美日韩无卡精品| 一边摸一边抽搐一进一小说| 一级毛片精品| 少妇人妻一区二区三区视频| a在线观看视频网站| 亚洲av成人一区二区三| 国产av一区二区精品久久| 成人av在线播放网站| 超碰成人久久| 国产精品一区二区三区四区久久| 亚洲成a人片在线一区二区| 日本成人三级电影网站| 国产三级黄色录像| 国产成人啪精品午夜网站| 大型黄色视频在线免费观看| 色播亚洲综合网| 亚洲精品av麻豆狂野| 亚洲精品久久成人aⅴ小说| 国产精品美女特级片免费视频播放器 | 成人18禁高潮啪啪吃奶动态图| 亚洲精品在线观看二区| 久久精品成人免费网站| 久久久久久久久中文| 少妇被粗大的猛进出69影院| 中文字幕精品亚洲无线码一区| 午夜激情av网站| or卡值多少钱| 制服人妻中文乱码| 免费看a级黄色片| 国产精华一区二区三区| 国产亚洲精品一区二区www| 18禁国产床啪视频网站| 国产日本99.免费观看| 黑人巨大精品欧美一区二区mp4| 亚洲av日韩精品久久久久久密| 婷婷丁香在线五月| 老司机深夜福利视频在线观看| 99re在线观看精品视频| 男插女下体视频免费在线播放| 精品久久久久久久久久久久久| 每晚都被弄得嗷嗷叫到高潮| 日日干狠狠操夜夜爽| 日韩大尺度精品在线看网址| av国产免费在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 色老头精品视频在线观看| 床上黄色一级片| 久久伊人香网站| 国产男靠女视频免费网站| 久久九九热精品免费| 老司机深夜福利视频在线观看| 久久久水蜜桃国产精品网| 色尼玛亚洲综合影院| 久久欧美精品欧美久久欧美| 在线观看午夜福利视频| 久久精品国产综合久久久| 免费在线观看完整版高清| 国产v大片淫在线免费观看| 欧美一区二区精品小视频在线| 老鸭窝网址在线观看| 精品久久久久久久人妻蜜臀av| 国产精品久久久人人做人人爽| 亚洲免费av在线视频| 成人特级黄色片久久久久久久| 搡老妇女老女人老熟妇| 午夜精品一区二区三区免费看| 搡老岳熟女国产| 女警被强在线播放| 全区人妻精品视频| 岛国在线免费视频观看| 久久久久久久久中文| 88av欧美| 欧洲精品卡2卡3卡4卡5卡区| 此物有八面人人有两片| 国产成人精品久久二区二区91| 国产精品电影一区二区三区| 国产蜜桃级精品一区二区三区| 97超级碰碰碰精品色视频在线观看| 妹子高潮喷水视频| cao死你这个sao货| 男男h啪啪无遮挡| 美女高潮喷水抽搐中文字幕| 欧美成人免费av一区二区三区| 国产精品久久久久久亚洲av鲁大| 禁无遮挡网站| 成年版毛片免费区| 香蕉丝袜av| 免费在线观看日本一区| 一二三四在线观看免费中文在| 亚洲精品美女久久久久99蜜臀| 亚洲国产中文字幕在线视频| 97碰自拍视频| 桃红色精品国产亚洲av| 亚洲av五月六月丁香网| 久久中文字幕一级| 日日夜夜操网爽| 精品乱码久久久久久99久播| www国产在线视频色| 久久精品人妻少妇| 久久九九热精品免费| 国产精品一区二区精品视频观看| 国产黄色小视频在线观看| 久久 成人 亚洲| 黄片小视频在线播放| 亚洲精品国产一区二区精华液| 亚洲乱码一区二区免费版| 欧美日韩福利视频一区二区| 女同久久另类99精品国产91| 免费看美女性在线毛片视频| 日日摸夜夜添夜夜添小说| 精品国内亚洲2022精品成人| 丰满的人妻完整版| 婷婷丁香在线五月| 国产精品久久久久久亚洲av鲁大| 欧美色视频一区免费| 动漫黄色视频在线观看| 午夜亚洲福利在线播放| 亚洲成av人片在线播放无| 一个人免费在线观看的高清视频| 成人欧美大片| 19禁男女啪啪无遮挡网站| 国内精品久久久久精免费| 国产私拍福利视频在线观看| svipshipincom国产片| 欧美在线黄色| 精品欧美一区二区三区在线| 国产亚洲精品第一综合不卡| 午夜老司机福利片| 久久精品人妻少妇| 99精品久久久久人妻精品| 亚洲精品在线美女| 99riav亚洲国产免费| 欧美 亚洲 国产 日韩一| 淫秽高清视频在线观看| 丁香六月欧美| 亚洲18禁久久av| 国产aⅴ精品一区二区三区波| 国产精品av视频在线免费观看| 亚洲欧美精品综合久久99| 午夜福利18| 国产视频内射| 嫩草影视91久久| 久99久视频精品免费| 亚洲国产高清在线一区二区三| 午夜福利在线观看吧| 国产黄片美女视频| 日本a在线网址| 欧美 亚洲 国产 日韩一| 欧美日韩亚洲综合一区二区三区_| 色综合婷婷激情| 九九热线精品视视频播放| 窝窝影院91人妻| 床上黄色一级片| 亚洲精品av麻豆狂野| 变态另类丝袜制服| 中文字幕久久专区| www.999成人在线观看| 国产精品一区二区免费欧美| 美女黄网站色视频| 国产又黄又爽又无遮挡在线| 香蕉av资源在线| 国产片内射在线| 日韩欧美一区二区三区在线观看| 欧美中文综合在线视频| 美女午夜性视频免费| 色综合欧美亚洲国产小说| 国产一区在线观看成人免费| 亚洲欧美精品综合久久99| 亚洲成人中文字幕在线播放| 狂野欧美激情性xxxx| 久久人人精品亚洲av| 少妇的丰满在线观看| 国产精品久久久av美女十八| 国产单亲对白刺激| 丁香欧美五月| 婷婷丁香在线五月| 男女下面进入的视频免费午夜| 色av中文字幕| 色噜噜av男人的天堂激情| 亚洲欧美精品综合久久99| 一a级毛片在线观看| 欧美中文综合在线视频| 婷婷六月久久综合丁香| 黄色丝袜av网址大全| 亚洲av美国av| 啦啦啦观看免费观看视频高清| 国产亚洲欧美在线一区二区| 中文字幕高清在线视频| 亚洲性夜色夜夜综合| 亚洲av片天天在线观看| 免费观看人在逋| 免费无遮挡裸体视频| 1024视频免费在线观看| 亚洲精品在线美女| 久99久视频精品免费| 婷婷精品国产亚洲av在线| 欧美3d第一页| 日韩精品免费视频一区二区三区| 久久伊人香网站| 男女床上黄色一级片免费看| 成人一区二区视频在线观看| 国产精品精品国产色婷婷| 少妇熟女aⅴ在线视频| 麻豆av在线久日| 人妻夜夜爽99麻豆av| bbb黄色大片| 国产三级在线视频| 校园春色视频在线观看| 一进一出好大好爽视频| 九色国产91popny在线| 99在线人妻在线中文字幕| 99久久99久久久精品蜜桃| 国产一区二区激情短视频| 久久 成人 亚洲| 国内毛片毛片毛片毛片毛片| 好看av亚洲va欧美ⅴa在| 亚洲av电影在线进入| 窝窝影院91人妻| 国产精品精品国产色婷婷| 亚洲欧美一区二区三区黑人| 精品福利观看| 亚洲国产欧美人成| 精品久久久久久久久久久久久| 久久精品国产亚洲av香蕉五月| 欧美久久黑人一区二区| 麻豆国产97在线/欧美 | 91九色精品人成在线观看| 国产伦人伦偷精品视频| 欧美最黄视频在线播放免费| 视频区欧美日本亚洲| 欧美大码av| 欧美一区二区国产精品久久精品 | 欧美日本视频| av视频在线观看入口| 日日摸夜夜添夜夜添小说| 亚洲性夜色夜夜综合| 日韩大尺度精品在线看网址| 国产精品av久久久久免费| 国产精品国产高清国产av| 亚洲九九香蕉| 白带黄色成豆腐渣| 一a级毛片在线观看| 首页视频小说图片口味搜索| 两个人免费观看高清视频| 免费在线观看日本一区| 成在线人永久免费视频| 欧美日韩精品网址| 麻豆成人午夜福利视频| 成人国产一区最新在线观看| 嫁个100分男人电影在线观看| 欧美一区二区精品小视频在线| 又爽又黄无遮挡网站| 国产精品久久久久久精品电影| 亚洲欧美精品综合久久99| 美女高潮喷水抽搐中文字幕| 国产成人精品久久二区二区免费| 19禁男女啪啪无遮挡网站| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲aⅴ乱码一区二区在线播放 | 丰满人妻一区二区三区视频av | 国产探花在线观看一区二区| 九色成人免费人妻av| 黄色a级毛片大全视频| 久久久久久大精品| 国产97色在线日韩免费| 国产欧美日韩一区二区三| 91老司机精品| 国产真人三级小视频在线观看| 别揉我奶头~嗯~啊~动态视频| 国产精品98久久久久久宅男小说| 亚洲精华国产精华精| 一夜夜www| 法律面前人人平等表现在哪些方面|