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

    On abnormal behaviors of ion beam extracted from electron cyclotron resonance ion thruster driven by rod antenna in cross magnetic field

    2021-08-05 08:29:54JuanLI李娟ShenghuiFU付省輝YurouYANG楊雨柔andZhenfengDING丁振峰
    Plasma Science and Technology 2021年8期
    關鍵詞:李娟

    Juan LI (李娟), Shenghui FU (付省輝), Yurou YANG (楊雨柔) and Zhenfeng DING (丁振峰)

    School of Physics, Dalian University of Technology, Dalian 116024, People’s Republic of China

    Abstract In a 2.45 GHz electron cyclotron resonance(ECR)ion thruster powered with rod antenna under a cross magnetic field, abnormal behaviours such as sudden drop of ion beam current (Ib) and larger increasing-rate of Ib in the high microwave power (Pw) discharges at high gas flow rates were observed.A differential method was proposed to reveal the changes in the radial profiles of gray values extracted from the end-view discharge images.The increasing-rate of Ib with respect to Pw was used to evaluate efficiencies of ion production and transport.Analyses indicate that discharges are dominantly sustained by ordinary wave via electron heating in the electron plasma resonance layer that can shift along the rod-antenna, and extraordinary wave can only ignite a discharge in the ECR layer in the low gas flow rate regime.In terms of the confinement region defined by the magnetic field lines intercepting with the screen grid, the confinement region of the optimized 2.45 GHz cross magnetic field takes the shape of hourglass, enabling the high increasing-rate of Ib with respect to Pw in high power discharges at high gas flow rates.Correlated with the accompanied bright boundary layer appearing in the differentiated image,the sudden drop of Ib in the low gas flow rate regime is attributed to the discharge ignited by the enhanced extraordinary wave in the ECR layer neighbouring the narrowest confinement region,where the produced ions can promptly enter the loss region.

    Keywords: ECR ion thruster, ion beam current, discharge mode transition

    1.Introduction

    Gridded microwave electron cyclotron resonance (ECR) ion sources have advantages of both electron heating and magnetic confinement, and therefore are widely applied in significant fields including electric thrust of satellites.To control orbit and attitude of micro satellites,small size,low discharge power and high-performance ECR ion thrusters were developed [1–15].Among the smallest (2 cm) ECR ion thrusters with various magnetic field configurations [5–7, 14, 15], the highest ion beam current was achieved in the so-called cross magnetic field[14,15],produced by two sets of magnets with magnetizations in two orthogonal directions(axial and radial).One reason for this is the axial and radial magnetic-gradient forces that respectively reduce plasma loss to the chamber walls and drive ions to move towards to the extraction grids[14].In the optimization of the 2.0 cm cross magnetic field ECR ion thruster, ion beam current and plasma impedance were studied as functions of discharge power, gas flow rate and microwave frequency [14, 15].The highest ion beam current achieved at the driving frequency of 1.6 GHz was attributed to the optimized ECR layer that was away from the solid surfaces of the rod antenna and the discharge side wall[14].However, the side-view image of the discharge demonstrated that the most luminous region was not in the expected ECR layer (see figure 2(c) in [14]), but in the vicinity of the rod antenna, implying that the practical discharges excited with the rod antenna in the cross-field configuration remain to be identified.

    In designing a 5 cm ECR ion thruster powered with a rod antenna, we improved the cross magnetic field configuration by using a closed hollow magnet-cylinder consisting of eight arc-shaped magnets and covering the undesired ECR layer with aluminum plate [13].The rod-antenna discharge in the 5 cm cross-magnetic-field ECR ion thruster was failed because an arcing was ignited on the front end of the antenna in case of ion extraction.The problem was solved by replacing the rod-antenna with a circular plate one.We designed a 2 cm 2.45 GHz ECR ion thruster with a cross magnetic field configuration improved in the same way as for the above 5 cm ion source and successfully extracted ions from discharges excited with a rod antenna.

    Plasmas in gridded ion thrusters are floated at the high positive voltage of the screen grid,so that diagnosing them by intrusive instruments such as electric probe is technically difficult.As the remote monitors of ion thrusters, commercially-available cameras have been used to capture discharge images.However, the raw images have been not further analyzed to gain quantified data about the monitored discharges.

    Although ECR ion thrusters have been researched for more than two decades [16], and μ10 ion engines had been employed in space as the main thrust system for Hayabusa[17], ECR ion thrusters are still less developed [18].Optimizations of ECR ion thrusters rely highly on experimental test, because of the diagnostic difficulty and complicated processes occurring in the interactions between electromagnetic waves and anisotropic magnetized ECR plasmas.In contrast, investigations on mechanisms underlying ECR ion thrusters are insufficient.Effects of magnetic configuration on ion beam current from ECR thrusters have been experimentally studied [3, 14, 19–22], but the involved microwave mode transitions are not specifically analyzed.As a result,the dependence of ion beam current on discharge power is not well understood.

    In this paper,we report the observed abnormal behaviors and the characteristic regions of ion beam current against discharge power, propose a differential method to analyze gray values of discharge images, and establish a simplified physical model to describe plasma production and ion transport cross magnetic field, and reveal their effects on ion extraction.

    2.Experimental apparatus and diagnostics

    Figure 1.Schematic diagram of experimental system.MPG:Microwave Power Generator, MPM: Microwave Power Meter,TSCT: Three Stub Coaxial Tuner, DCB: Direct Current Block, SG:Screen Grid,AG:Accelerator Grid.Note that the ion thruster and the vacuum chamber are not plotted to scale.In fact, the former is smaller than one tenth of the latter.

    The experimental system for the 2 cm microwave ECR ion thruster is shown in figure 1.The vacuum chamber has a diameter of 510 mm and a length of 660 mm and can be evacuated to a background pressure of 2.0×10?4Pa by a vacuum pumping system consisting of a turbo molecular pump (1000 l s?1) and a mechanical pump (14 l s?1).Xe gas injected into the discharge chamber of the ion thruster is regulated by mass flow controller (Beijing Horiba Metron,S49 35/MT).2.45 GHz microwave from a microwave power generator is fed into the discharge chamber via a microwave power meter (Bird Technologies 7022), a three-stub coaxial tuner, a direct current block and a molybdenum (Mo) rod antenna.A direct current block is used to isolate the current path between the ion thruster and the grounded vacuum chamber.The discharge chamber of the ion thruster is electrically connected to the screen grid at a positive voltage(+1100 V in experiments), and thus the plasma is floated at the positive voltage of the biased screen grid.A negative voltage of –250 V is applied to the acceleration grid to suppress back-streaming electrons.A grounded shielding-cage is mounted to protect the electrically-floated ion source from the bombardment of secondary electrons emitted from the metal surface of the vacuum chamber under the bombardment of extracted ions.The ion beam current (Ib) is calculated fromIs–Ia,hereIsandIaare the currents collected by the screen and acceleration grids, respectively [14].

    The structure of the ion thruster and the magnetic field configuration are shown in figure 2, where the directions of the black arrows denote those of the magnetization.The sizes of key components in the microwave ECR ion source are as follows: the discharge chamber of the ion thruster has a diameter of 21.0 mm and a length of 14.0 mm; the diameter of the ion extraction region is 18.0 mm; the diameter and the length of the rod Mo antenna are 3.0 mm and 10.0 mm,respectively;the gap between the antenna’s front end and the screen grid is 4.0 mm.

    As shown in figure 2,the magnetic field in the discharge chamber is produced by an annulus magnet at the rear end and a hollow-cylinder magnet composed of eight arc-shaped magnets surrounding the side wall of the discharge chamber.The cross magnetic field was calculated by the AC/DC module of the COMSOL software [23].In optimizing the cross magnetic field,the 2.45 GHz ECR layer was designed to locate at the midway between the rod antenna and the side wall of the discharge chamber.The optimization was performed by gradually adjusting the parameters of two magnets,namely inner radius, outer radius, axial length and remnant magnetization.Under the optimized condition, the achieved magnetic field distribution is shown in figure 2, where the white lines represent the magnetic field lines,and the red one for the location of the ECR layer corresponding to 2.45 GHz(B=0.0875 T).Obviously, the ECR layer zigzags around at the midway between the rod antenna and the side wall of the discharge chamber.

    Figure 2.Magnetic field distribution of in azimuthal cross-section.AB: Aluminum Block, AT: Antenna, AM: Annular Magnet, ASM:Arc Shaped Magnet, CW: Cylindrical Wall, PD: Polytetrafluoroethylene Dielectric, SB: Soft-iron Block, SSB: Stainless Steel Block, FP: Front Plate.The white arrows represent magnetic field lines; the red line denotes the ECR layer; the black arrows indicate the magnetization direction.

    A digital camera (Canon EOS 50D, 22.3×14.9 mm CMOS sensor array consisting of 4752×3168 pixels) was installed to capture end-view images of discharge glows (see figure 1).The lens aperture and the shutter speed of the camera were set at F5.6 and 1/200 s, respectively.The recorded 8 bit RGB (red, green and blue) true-colour images in JPEG (joint photographic experts group) format were converted to 8 bit grayscale (0–255 scale) files by the image processing software Image J using the average method [24],enabling quantitative analyses of radial plasma profiles.Unlike ordinary glow images, optical emissions from discharge chambers of gridded ion sources have to pass through the ion extraction holes of the screen and acceleration grids that actually act as a simple optical collimator.In this case,the spatial distribution of gray values extracted from the pixel array of an end-view image can roughly represent the radial profile of the glow.The intensity of optical emission from Xe plasma depends on Xe atom density, electron density and electron energy [25].In a small ECR ion source operating at low gas pressure,the pressure gradient of injected Xe atoms is not high, so that the local emission intensity is mainly determined by the local electron density (ne) and electron temperature (Te) [26].Therefore, the radial profile of gray values extracted from the end-view glow image can approximately represent that of the discharge intensity(defined roughly asneTe).

    Figure 3.Ion extraction currents as functions of microwave discharge power under different Xe gas flow rates.The dashed-lines denote the intersecting points in 0.3 and 0.35 sccm curves in the high-power regime.

    3.Experimental results and discussion

    3.1.Ion beam current

    In the range of Xe gas flow rate from 0.15 to 0.35 sccm, the obtainedIb(defined in the above section) against the microwave discharge power (Pw=Pin–Pref, herePinandPrefare incident and reflected microwave powers, respectively) is shown in figure 3.At 0.15 sccm,Ibvaries almost linearly withPwin the lowPwregime below 2.82 W.Above this value,Ibenters a new quasi-linear region with a lower curve slope.WhenPwfurther rises to 5.01 W,Ibdrops suddenly, accompanied by the sudden increments in both the incident and reflected microwave powers, indicating that the microwave discharge itself also varies discontinuously.In the highPwregime above 6.66 W, the curve ofIbversusPwresumes the quasi-linearity with a slope almost equal to that in prior to the sudden drop.

    As the gas flow rate increases to 0.20 sccm, the dependence ofIbonPWis qualitatively the same as that at 0.15 sccm, i.e.the segmented linear regimes and the sudden drop ofIb.Comparing the results at 0.15 and 0.2 sccm, one finds that the criticalPWat the sudden drop ofIbdecreases with the increasing gas flow rate (gas pressure).When the flow rate further rises to 0.25 sccm,the curve ofIbversusPWbecomes continuous, i.e.the sudden drop ofIbdisappears.Obviously, the continuous curve for 0.25 sccm consists of two quasi-linear curve segments intersecting at a smallerPWcompared with low gas flow rates.As the gas flow rate enters the high gas flow regime of 0.3–0.35 sccm, the intersecting point of the two quasi-linear curves shifts gradually to an even lowerPW.Comparing all the intersecting points,one can find that they are at an almost equalIb, the reason for which will be analysed in section 3.5.4.

    Figure 4.Variations of total gray value extracted from glow images with microwave power under different Xe gas flow rates.

    In the high power discharges at 0.3–0.35 sccm, the third quasi-linear region, which is beyond the second intersecting point denoted with the vertical dotted line in figure 3,appears with a higher curve slope than that in the secondary quasilinear region.ThePWat the second intersecting point also decreases with the increasing gas flow rate.Compared with the 1.6 GHz ECR ion thruster with the cross magnetic field[14], theIbfrom the 2.45 GHz one exhibits different behaviors in two aspects:(1)Ibis a monotonic increasing function of gas flow rate at equalPw;(2)in the high gas flow regime of 0.3–0.35 sccm, the increasing rate ofIbwithPwin the third linear curve segment is abnormally larger, being contrary to the low efficiency of ion production in high-Pwdischarges of conventional ECR ion thrusters [3, 27–30].

    3.2.Total gray value of discharge image

    To gain more insight into the discharge processes, all glow images corresponding to the measured ion beam currents were recorded simultaneously.The extracted total gray-values using Image J software are shown in figure 4.Comparing the two sets of curves in figures 3 and 4,one can find that,in the low-Pwregime, all the curves of total gray value versusPware similar to those ofIbversusPW, i.e.an intersecting point lies between the two quasi-linear curves.More worthy of noting is that the corresponding intersecting-points in the two sets of curves are at an equalPwfor all gas flow rates.However,during the sudden drop ofIbat 0.15 and 0.2 sccm,the corresponding changes of the total gray-values are not as pronounced, but remain constant and decrease slightly,respectively.In the high gas flow region (0.25–0.35 sccm),the total gray-values are the monotonic increasing functions ofPw, but in the high-Pwregime for 0.3 and 0.35 sccm, the shapes of curves of total gray-value versusPwobviously differ from those ofIbversusPw.

    3.3.Primary analysis of discrepancy between ion extraction and optical emission

    The above discrepancy betweenIband the total gray-value is related to the different processes in ion extraction and optical emission.During the ion transport from the upstream, where ions (electrons) are generated, to the screen grid at the downstream,some of the ions(electrons)will lose on the side wall of the discharge chamber, and accordingly make no contribution toIb.In terms of residence time,such ions can be regarded to have a short lifetime.Essentially,the measuredIbis the current of long-lived ions that can arrive at the screen grid and subsequently pass through the acceleration grid.In contrast, the processes in the optical emission are more complicated.Firstly, ion emission makes less contribution to the total optical emission in the case of low ionization degree,while Xe atom emission dominates.Secondly, ions and electrons are confined in the magnetic field designed for 2.45 GHz ECR, but Xe atoms are not directly affected by magnetic confinement, so that the transport properties of particles dominating ion extraction and optical emission are different.Thirdly, the lifetime of excited Xe atom due to spontaneous emission is on the order of nanosecond, being much shorter than that of ion (electron) lost on the side wall.WhenIbis considerably reduced due to ion loss to the side wall (refer the curves for 0.15 and 0.2 sccm in figure 3), the corresponding total intensity of the optical emission (roughly proportional to the total gray value) might drop slightly or remain constant (refer the curves for 0.15 and 0.2 sccm in figure 4), because electrons (ions) transporting across magnetic field lines to the side wall still make contributions to the optical emission.AsPwrises in a same discharge mode, the plasma profile would evolve continuously.Consequently,ion loss flux to the side wall rises at a rate that is roughly proportional to ion density[31].The sudden drop ofIbin figure 3 indicates that the discharge features such as the radial plasma profile should simultaneously undergo steep variations during an unrevealed discharge transition.

    3.4.Differentiated images between two consecutive discharges

    To validate the above analysis, the radial profiles of the discharges in figure 3 should be investigated beside their total gray-values shown in figure 4.The discharges at 0.15 and 0.35 sccm are selected as the representative cases, where the sudden drop ofIbappears and disappears, respectively.The unprocessed discharge images are shown in figures 5 and 6.In the two cases, one can find out with naked eye that the optical emission intensity increases globally with the increasingPw,but the subtle difference between two adjacent images taken with a smallPwstep is almost indiscernible.Even during the sudden drop ofIb(see figures 5(g) and (h)),the variation in the radial profile of optical emission,which is informative of the discharge transition, is still not obvious to naked eye.

    Figure 5.Images of discharge glow corresponding to those shown in figures 3 and 4 at gas flow rate of 0.15 sccm.These photos are 8 bit RGB true-colour images in JPEG format.The sudden drop of ion beam current occurs in the discharge transition from (g) to (h).

    Figure 6.Images of discharge glows corresponding to those shown in figures 3 and 4 at gas flow rate of 0.35 sccm.These photos are 8 bit RGB true-colour images in JPEG format.

    Figure 7.8 bit grayscale images of differentiated gray values from raw images in figure 5.Refer the text for the procedures to achieve the differentiated gray values.Note that figure 7(g) corresponds to the sudden drop of ion beam current.

    In this situation, quantifying radial profile by software Image J is needed.In the initial attempt, the radial profiles of gray values were extracted along the selected extraction holes in a row that passed through the centre of the grids.The obtained radial profiles were unsatisfying because the difference between two consecutive discharges was found to be irregular, depending on the azimuthal orientation of the selected hole-array due to the azimuthally non-uniform gas injection.The further attempt was made to achieve two-dimensional (2D) profiles including the azimuthal direction,but the difference between the absolute gray-values from two consecutive discharges was also less distinguishable.To remove the above drawbacks,a differential method is proposed to analyse two consecutive images, the procedures of which are as follows: (1) extracting the 2D gray values from the images shown in figures 5 and 6;(2)subtracting the gray values of the image in a lowerPwdischarge from those of the adjacent higherPwdischarge; (3) the negative data from the above differential calculation are artificially set to zero, while the positive ones retain their original values.The differentiated gray values achieved in the above procedures are plotted as 8 bit grayscale images (see figures 7 and 8).

    According to the differential method, the gray values increase with the increasingPwin the bright regions of differentiated images shown in figures 7 and 8.From the primary analysis in section 3.3, one knows that the variation in the gray-value approximately represents that of electron density.At 0.15 sccm, the bright regions across the entire radial domain(see figures 7(a)–(f)and(h)–(l))indicate that electron density increases with the increasingPw.The qualitatively unchanged profiles of the differentiated images in the twoPwregions suggest that the discharge modes do not vary abruptly.In the narrow windows ofPw(see figures 5(h)–(g)and 7(g)), however, the electron density decreases and increases, respectively, in the black central region and the bright boundary layer.This transition of the differentiated image profile is a physically reasonable result of the accompanied sudden drop ofIb.Moreover, the bright boundary layer is also consistent with the drop ofIb, the reasons for which are primarily analyzed in section 3.3 and will be further elucidated in section 3.5.4.

    Two possibilities for the increase of electron(ion)density in the bright boundary layer in figure 7(g) are: (1) electrons(ions)produced in the upstream are more easily transported to the boundary region in the downstream; (2) more power is dissipated to produce electrons(ions)in the boundary region.In either case, the ion loss flux to the side wall increases,leading to the decrement in ion current extracted by the biased grids at the downstream.

    Figure 8.8 bit grayscale images of differentiated gray values from raw images in figure 6.Refer the text for the procedures to achieve differentiated gray values.

    At 0.2 sccm, the differentiated images (not presented in the paper) are qualitatively the same as those shown in figure 7.At 0.35 sccm, the bright boundary layer disappears from the differentiated images (see figure 8), coinciding with the absence of the ion current drop.

    As stated in section 3.1,the curves ofIbversusPwin the high flow rate regime of 0.3–0.35 sccm (see figure 3) can be classified,in terms of the curve slope,into three characteristic regions.The smallest curve slope in the moderatePwregion implies that the plasma density should start to increase in the boundary region, as evidenced by the widest bright region of the differentiated image (see figure 8(g)).

    Despite that the differentiated images offer an insight into the correlation between the ion extraction and plasma production, but the underlying mechanism is unknown.

    3.5.Physical model underlying discharge mode transition

    As stated in section 2, measuring the plasma profile in a small-sized ECR ion thruster is technologically difficult.In this case, numerical simulation is an alternative approach.However, the numerical modeling is also a tough task for low-pressure non-uniform ECR magnetoplasmas in the less collisional regime by 2D self-consistent PIC (particle-in-cell)or kinetic simulation of full-wave microwave fields [32, 33].In previously reported works on the sudden mode transitions of ECR discharges,numerical simulations were not successful in re-producing the experimental results.Taking the above into consideration, we discuss the correlation between ion extraction and plasma production as well as ion transportation by establishing a simplified physical model, and focusing attention on the significant features.

    Figure 9.Microwave electric field(yellow lines)calculated with the radio frequency module of COMSOL software and perpendicular component of microwave electric field (E⊥, illustrated with colour contour) in ECR layer. E⊥is defined as the component of electric field perpendicular to the static magnetic field.The white arrows represent magnetic field lines.

    The power deposition profile in electron heating and the distribution of ECR plasma depend on the properties of the wave modes excited by the rod antenna,such as propagation,cutoff and absorption.These properties, as indicated by the evolutions in figures 3–8, will also vary withPwor plasma density.Starting the analysis in the sequence from low to highPwis an intuitive selection.For the practical discharges,however, the low-Pwdischarge evolves from the high one after the gas breakdown.We therefore first analyze the ECR gas breakdown.

    3.5.1.Electric field and electron heating in low plasma density regime during the early phase of gas breakdown.During the pre-breakdown stage, the microwave electric field in the discharge chamber is less influenced by the charged particles with a very low density,and is almost the same as that in the vacuum case.The electric field distribution in the ECR vacuum chamber, calculated by the radio-frequency module of COMSOL software, is shown in figure 9 with the yellow line.The electric field has two characteristic regions: (1)region I between the front end of the rod antenna and the screen grid, where the electric field is most intensive and the majority of electric field lines are nearly parallel with magnetic field ones; (2) region II between the cylindrical surface of the rod antenna and the side wall of the chamber,where the electric field is relatively lower and the majority of electric field lines are almost vertical to magnetic field ones.

    The gas breakdown is not initiated in region I due to the following reasons: (1) the collisionless heating of cosmicradiation-produced electrons in the oscillating microwave electric field, which is almost in parallel with the magnetic field,is ineffective at low gas flow rate(gas pressure);(2)the gas gap in region I is on the order of 4 mm,being so short that the development of electron avalanche is substantially limited at low gas flow rate.The above hypothesis is also supported by the central dark regions appearing in the end-view discharge images (see figure 5).

    Unlike the electron heating in region I, electrons in region II are accelerated transversely in the microwave electric field that is vertical to the external magnetic field.In the ECR layer of ω=ωce(here ω and ωceare the microwave circular frequency and the electron cyclotron frequency,respectively), electrons rotating on the Larmor orbit are continuously accelerated in each half cycle by the linearlypolarized electric field.The electrons accelerated to a high transverse kinetic energy in the transverse electric field cannot escape rapidly from plasma due to slow transport across magnetic field lines, but have a long residence time.In this situation, the ionization by electrons heated in the so-called ECR layer is more efficient than that in region II, and consequently the gas breakdown is initiated at the ECR layer in region II.

    3.5.2.Exact location of gas breakdown in ECR layer.A new problem arises from the gas breakdown in region II:the ECR layer is distributed in the discharge chamber,so that the exact location of gas breakdown remains to be determined.From the transverse heating of electrons discussed above, such a gas breakdown should take place at one segment of the ECR layer where the vertical component of electric field(E⊥)with respect to the static magnetic field is highest.The calculated distribution ofE⊥over the ECR layer is depicted in a colour contour (see figure 9) plotted together with the electric field lines and magnetic field ones.One can see that the highestE⊥is around the position ofr=4.1 mm andz=9.5 mm.It is worth noting here that this location is also helpful in determining the ignition of the discharge during the sudden drop ofIb.

    3.5.3.Rod-antenna excited wave modes in plasma confined by cross magnetic field.After the gas breakdown, the plasma produced in the ECR zone of the highestE⊥is transported along and across magnetic field(see figures 9 and 10(b)).For the magnetic field designed for 2.45 GHz ECR, the transport coefficient across the magnetic field is much lower than that along the field.In fact, the transient plasma transport during gas breakdown is rather complicated, so do the evolutions of the wave modes in the discharges.A thorough description of the above processes is beyond the scope of this paper.Our focus is set on the discussion of the wave modes in the representative cases.

    Figure 10.Schematic illustrations of plasma density distributions in low (a), moderate (b) and high (c) discharge power regimes.In the three-dimensional space, the discharge chamber is divided by the separatrix surface consisting of the magnetic field lines denoted with cyan colour,into the confinement region(inside)and the loss region one (outside).The cyan-coloured magnetic field lines are those that just intersect the corner between the screen grid and the front plate.The yellow arrows represent the direction of wave propagation.For conciseness, the reflection along the rod antenna is neglected.The ECR ignition in (b) only occurs in the low gas flow rate regime.

    When region II defined above is filled with plasma, the wave excited therein by the rod antenna is named as X-(Extraordinary) mode.In the magnetic field configuration shown in figure 9, the magnetic field around the rod antenna is lower than that at the side wall of the chamber, forming a so-called low-field-side injection for the X-wave [34].In this case, the injected X-wave will sequentially encounter, in the Clemmow-Mullaly-Allis diagram [35], the cutoff layer of right-hand polarized wave and the upper hybrid resonance(UHR) layer, before reaching the ECR layer.In small-sized discharge chambers such as 2 cm diameter, the gap between the cutoff layer and the UHR layer is so narrow that the evanescent X-wave penetrating through the cutoff layer can be absorbed at the UHR layer, producing plasma therein.As the plasma density increases with the progress of gas breakdown, the UHR layer retreats together with the cutoff layer to the low magnetic field side, namely to the rod antenna.Therefore, the plasma produced at the UHR layer also contracts to the rod antenna.

    The magnetic field at the cylindrical surface of the rod antenna is 0.078 T, which is slightly lower than 0.0875 T for the 2.45 GHz ECR From the UHR frequencyone knows that, the UHR layer should occur between the ECR layer and the rod antenna (see figure 9).Owing to the high magnetic field, the plasma density at the UHR layer is accordingly rather low.As the plasma density further increases with the development of the gas breakdown, the UHR layer will move closer to the rod antenna, where a strongly non-uniform plasma sheath is formed in the transverse magnetic field.Owing to the non-uniform plasma,the UHR layer near the rod antenna is narrow.Consequently,the absorption of X-wave in the thin UHR layer with a low density is so weak that one should seek a new approach of wave absorption in the high plasma density regime.

    In the images shown in figures 5 and 6, all the brightest regions are always very close to the rod antenna, indicating that there exist ionization (excitation) layers.In fact, the plasma surrounding the rod antenna is confined by the parallel magnetic field whose force lines intercept almost vertically with the extraction grid(see figures 9 and 10).Consequently,the mobile electrons form conducting channels connecting the extraction grid.In this situation,electric field lines originating from the rod antenna’s cylindrical surface no longer terminate at the side wall as in vacuum case(see figure 9),but first pass through the narrow plasma sheath,and immediately deflect to the extraction grid along the conducting plasma channels surrounding the antenna.In the plasma channels where electric field lines are parallel with magnetic field ones, the electromagnetic wave propagates in the form of ordinary (O)wave.As stated above, the microwave power cannot be effectively absorbed in the UHR layer, so that the O-wave in the conducting plasma channels surrounding the rod antenna(see figure 10(a))makes the dominant contribution to electron heating.

    Like the early phase of gas breakdown when the seeding electron is rarefied,in low plasma-density discharges,O-wave electron heating around the rod antenna is also weak.After the ECR gas breakdown, the plasma density reaches a high level.Taking the 0.15 sccm discharge as an example,the gas breakdown occurs at aPwof 19.64 W, which is already beyond the highPwregion in figure 3.In experiments shown in figure 3,Pwwas reduced from the high value after the gas breakdown to the low one for ion extraction.

    As the electron density around the rod antenna increases to the critical value satisfying the equation of ofωpe=ω(hereωpeandωare electron plasma frequency and angular microwave frequency, respectively), the electron plasma resonance (EPR) layer appears in the conducting plasma channel where the local microwave electric field is maximized, and thus the electron heating is simultaneously enhanced [36].When the bulk plasma is in the overdense regime ofωpe>ω,the EPR layer will always occur around the rod antenna, sustaining the plasma production therein.Owing to the strong absorption at the EPR layer, the wave penetrating the EPR layer is evanescent, so that electron heating in the ECR layer is not the dominant one.

    It is worth noting again that, as evidenced by the central dark regions of all the images in figure 5,the EPR layer does not appear in the short gap between the front end of the rod antenna and the screen grid.

    3.5.4.Significant features of ion extraction.As aforementioned, plasma dynamics in ECR discharge is rather complicated, therefore the discussion below is focused the significant features of ion extraction.

    3.5.4.1.The high slope in curve of Ibversus Pwin the low power regime.Ion beam current is involved with ion production, ion transportation and ion extraction.The effects of the above processes onIbare analyzed in the followings.

    In the lowPwregime, the bright regions appear close to the rod antenna (see figures 5 and 6), indicating that all discharges are sustained by EPR electron heating.Consequently,the efficiency of ion production in the lowPwregime changes slightly.The ion transport to the screen grid in the lowPwregime is efficient, the reasons for which are as follows:(1)the plasma produced in the EPR layer around the rod antenna is magnetically confined and driven by the gradient force of the divergent magnetic field, facilitating the fast transport toward the extraction grid;(2)all magnetic field lines, in both the EPR layer and the adjacent region filled by the plasma transported across the magnetic field lines,intercept with the extraction grid (see figure 10(a)).As a result,almost all ions produced in the upstream can reach the screen grid.

    From the principle of ion extraction optics, one knows that not all ions arriving at the screen grid can be extracted out of the acceleration grid due to insufficient or the excessive ion focusing.In the low plasma density (lowPw) regime, the crossover ion impingement occurs.With increasingPwin the low plasma density regime, the percentage of ions that can pass through the acceleration grid increases.From the above discussion, one knows that, in case that the efficiencies of both ion production and ion transport do not change much,the decreasing crossover ion impingement due to the increasing plasma density with the increasingPwwill lead to the rapid rise ofIb.This deduction is also supported by the fact that the intersecting points of allIb?Pwcurves in the lowPwregime(see figure 3)are at almost an equalIb.This result is reasonable because the ion current limit of the crossover ion impingement is equal for a same extraction grid, and should be independent off plasma production at different gas flow rates (gas pressures).In addition, the shift of the intersecting point to the lowerPwwith the increasing gas flow rate (see figure 3)is also reasonable simply due to the increased plasma density at a higher gas flow rate (gas pressure).

    3.5.4.2.Sudden drop of Ib.When a low power is injected into the small region surrounding the rod antenna, the discharge can be operated in the overdense regime ofωpe>ω.As stated above, the discharge in this regime is characterized by the bright layer surrounding the rod antenna.The presence of the bright layer at all conditions indicates that the EPR electron heating contributes to all discharges.The arising question is that whether the ECR electron heating can coexist with the EPR one in thePwregime.The answer is that ECR electron heating does not exist in thePwregime before the drop ofIb,but occurs to induce the sudden drop ofIb,the involved processes in which are presented below.

    Although the gas breakdown is initiated at the ECR layer when the density of seed electrons is very low, the ECR electron heating does not always outmatch the EPR one as electron density is increased.In the lowPwregime, as the microwave discharge is strongly confined around the rodantenna, the electric field is already intensified in the EPR layer.In contrast,the electric field at the large radius is much low due to the strong wave absorption at the EPR layer, so that the ECR electron heating is weak.

    AsPwincreases,the variations that facilitate the ignition in the ECR layer occur in three aspects.Firstly, the electric field in the ECR layer rises simply due to the increase ofPw.Secondly,the seed electron density in the ECR layer increases as a result of plasma transport across magnetic field lines.Thirdly, the EPR layer shifts upstream as the conducting plasma channels become longer at higherPw,because plasma transport along the magnetic field lines is much fast than that across field lines.Since the EPR layer is a special location where the microwave is absorbed, cut-off and reflected, the microwave cannot propagate deep into the bulk plasma in the downstream, leading to the so-called skin effect (see figure 10(b)).The microwave reflected back to the upstream region enhances the electric field in the ECR layer.

    Owing to the above three variations,the discharge can be ignited in the ECR layer.Despite that the exact prediction of the ignition location relies on the strict dynamic modeling of ECR plasma,but it is near the ECR zone of the highestE⊥in the vacuum case(see figures 9 and 10(b)).The magnetic field lines threading the ECR ignition region intercept the front plate (FP in figure 10(b)) mounted at the side wall in the downstream region.Consequently, many of the ions (electrons) produced in the ECR zone will move downstream along the diverging magnetic field lines,and are finally lost at the front plate, making no contribution toIb.After the ECR discharge is ignited, a partition of the injected microwave power is dissipated in the ECR layer, while that in the EPR layer accordingly decreases.Since ions produced in the ECR layer are less efficient for ion extraction,so that the rise ofPwleads to the drop ofIb.In contrast,when more electrons(ions)are transported to the front plate after the drop ofIb, the intensity of optical emission increases at the boundary region,forming a bright annulus in the differentiated image shown as figure 7(g).

    As the Xe gas flow rate(gas pressure)increases to above 0.2 sccm,two counteracting variations are:(1)the rise in gas pressure reduces the electric field required by the ECR ignition (named as the first effect); (2) the absorption of microwave wave in the EPR layer is increased while the reflected microwave is decreased (named as the second effect).This counteraction is involved in complicated kinetic processes, so that the prediction of the critical power for the ECR ignition is not simple.Comparing the critical powers at the sudden drop ofIbfor 0.15 and 0.2 sccm(see figure 3),one knows that the first effect dominates over the second one in this gas flow rate regime.

    As the Xe gas flow rate(gas pressure)further increases to or above 0.25 sccm, the sudden drop ofIbdisappears.There are two possibilities for this:in one case,the above first effect still dominates the second one, but the ECR electron heating takes its role at an even lowerPwin a continuously manner;in the other case, the second effect surpasses the first one, and the ignition at the ECR layer does not take place.We believe that the second case is more plausible because if ECR ignitions occur at even lowerPwfor 0.25 sccm and above,they will also lead to the remarkable reduction ofIbas in the cases of lower gas flow rates.

    From the above,one can conclude that the ECR electron heating is much more efficient in the gas breakdown where electron density is very low, but is later outmatched by the EPR one in the overdense plasma ofωpe>ω.In the low gas flow rate regime, the ECR electron heating can coexist with the EPR one,igniting its discharge and leading to the sudden drop ofIb.In the high gas flow rate regime, the weak ECR electron heating is unable to ignite its discharge.

    3.5.5.Abnormal high slope in curves of Ibversus Pwat high gas flow rate and discharge power.In figure 10, the cyancoloured magnetic field lines are those that just intersect the corner between the screen grid and the front plate.In the three-dimensional space,these lines form a separatrix surface,beyond which ions (electrons) will lose at the side wall and make no contribution to ion extraction.For sake of concise statement, the regions inside and outside of the cyan colour lines are named as confinement and lose ones, respectively.The confinement region takes the shape of hourglass characterized by two wide ends and a narrow waist.

    In explaining the high slope of theIb?Pwcurve in the lowPwregime,ion production and transportation in the wide confinement region are regarded as the underlying reasons.As the ECR discharge (together with the EPR one in the low gas flow rate regime) is ignited at the narrow waist of the hourglass-shaped confinement region,the transport efficiency of the produced ions to the screen grid is reduced due to the high ion loss along magnetic field lines that intercept with the front plate,leading to a low increasing-rate ofIbwith respect toPwbeyond the lowPwregime.

    If plasmas can be produced at the upstream wide-end of the hourglass-shaped confinement region, the transport efficiency of ions to the screen grid will be elevated again because more ions are confined in magnetic field lines intersecting with the screen grid (see figure 10(c)).In discharges at the high gas flow rates (0.3 and 0.35 sccm),the evolutions ofIbwithPwseem to undergo through all the above transitions.The arising question is why only the discharges at the high gas flow rates can be produced at the upstream wide-end of the hourglass-shaped magnetic field.

    As already discussed, with increasing thePw, the discharges are sustained by electron heating in EPR layer that moves upstream along the rod antenna.In low gas flow rate regime, the produced plasma is characterized by low electron density and high electron temperature.Both of them slow down the upstream shifting of the EPR layer.The reason for this is that,in the case of magnetic confinement,the large electron’s Larmor radius and the high electron-neutral collision frequency at the high electron temperature are factors promoting electron transport across magnetic field lines[37].When ions gain more energy in collision with hightemperature electrons, ion transport across magnetic field is also raised.Therefore,the percentage of ions entering the loss region increases, and the formation of an EPR layer close to the rod antenna requires a higherPw.In discharges at the gas flow rates below 0.25 sccm, the EPR layer can only shift upstream to the waist of the confinement region (see figure 10(b)) in thePwregion of figure 3.

    In the high gas flow rate regime, electron density and electron temperature become higher and lower, respectively.The produced plasma tends to build up around the rod antenna as a result of the reduced transport across magnetic field lines, facilitating the upstream shift of the EPR layer.In this way,the EPR layers in the discharges at the high gas flow rates can move to the upstream wide-end of the confinement region (see figure 10(c)), so that ions produced therein are efficiently transported to the screen grid along the magnetic field lines.Consequently,ion beam currents from the high-Pwdischarges at high gas flow rates increase abnormally at a higher rate with the increasingPw(see figure 3).Moreover,the difference between the two critical powers for the second high-slope regions in 0.3 and 0.35 sccm discharges is consistent with what discussed above: a highPwis required for the EPR layer in the discharge at the lower gas flow rate to enter the upstream wide-end of the confinement region (see figure 10(c)).

    3.5.6.Comparison with ion beam current of 1.6 GHz cross magnetic field ECR ion thruster.On the cross magnetic field ECR ion thruster optimized at 1.6 GHz [14], ion beam currents at the gas flow rates are saturated in the high-Pwdischarges(see figure 4(a)in[14]).The highestIbachieved at thePwof 8.0 W is about 15.2 mA (see figure 5 in [14]).On our cross magnetic field ECR ion thruster optimized at 2.45 GHz,Ibat equalPwstill rises with the increasing gas flow rate up to 0.35 sccm(see figure 3).Since our highestPwin the 0.35 sccm discharge is lower than 8.0 W,Ibis linearlyextrapolated to 8.0 W.The obtained 18.0 mA is 18.4%higher than the counterpart of the 1.6 GHz ion thruster.Considering thatIbis rising rapidly with the increasing gas flow rate (see figure 3), one might expect an even higherIbat 8.0 W.

    For 2 cm ECR ion thruster applied to microsatellites, the comparison between ion beam currents at lowPwand gas flow rate is practical.In the 2.0 W discharge at the gas glow rate of 15 μg s?1(about 0.15 sccm), the obtainedIbfrom 1.6 GHz ion thruster is 4.5 mA(see figure 4(a)in[14]),while the counterpart from 2.45 GHz ion thruster is 11.0 mA (see figure 3).The increment inIbby 144.3% is comparatively much higher.

    Both 1.6 and 2.45 GHz ion thrusters have the cross magnetic field,but the behavior difference between the two is large.In the 1.6 GHz ECR ion thruster with a relatively lower magnetic field, the transverse plasma-confinement was accordingly weaker.This is one of reasons for the smaller ion beam currents at low and high gas flow rates.Another reason lies in the configuration of the confinement region.In the 1.6 GHz ion thruster with a donut (annular in this paper)magnet of a smaller diameter (see figure 1 in [14]), the magnetic field lines are over concentrated and diverged in the upstream and downstream regions, respectively, both leading to the reduction of the corresponding confinement regions.In the low gas flow rate region,the percentage of ions produced in the downstream over-divergent magnetic field is higher,and the obtainedIbis lower.In the high-Pwdischarges at large gas flow rates, more of ions produced in the EPR layer in the upstream narrow confinement-region will be transported to the side wall,and thus the obtainedIbis saturated at a lower level.In comparison, our 2.45 GHz ion thruster is advantageous in terms of efficient ion extraction due to enlarged plasma confinement regions.

    4.Conclusion

    In a 2 cm 2.45 GHz microwave ECR ion thruster driven with rod antenna under cross magnetic field, ion beam current exhibits abnormal behaviors such as sudden drop and high increasing-rate with respect to discharge power at low and high gas flow rates, respectively.The proposed differential method is successful in revealing the variation in the radial profile of the gray values extracted from the discharge images during the sudden drop of ion beam current.The optimized 2.45 GHz cross magnetic field configuration, which has relatively wider confinement regions at the downstream and upstream of the chamber, is advantageous over the existing 1.6 GHz one in term of ion magnetic confinement.Ion beam current from 2.45 GHz ion thruster is higher at equal discharge power.The evolution of ion beam current with discharge power is correlated with the upstream shifting of electron-plasma-resonance layer through different segments of the hourglass-shaped confinement region.The sudden drop of ion beam current takes place as a result of the discharge ignited by the enhanced extraordinary wave in the ECR layer where the confinement region is narrowest.

    Acknowledgments

    This work was partially supported by National Natural Science Foundation of China (No.11975070) and Open Funds for Science and Technology on Vacuum Technology and Physics Laboratory, Lanzhou Institute of Physics (No.zwk1609).

    猜你喜歡
    李娟
    老虎萌萌的幸福生活
    學與玩(2022年2期)2022-05-03 00:08:32
    Therapeutic effect of heat-sensitive moxibustion plus medications for senile osteoporosis and its effect on serum BMP-2 and OPG levels
    那個不愿意幫忙帶孩子的婆婆
    那個不愿意幫忙帶孩子的婆婆
    幸福(2019年22期)2019-09-04 06:50:22
    李娟散文集《這世間所有的白》的哲理蘊藉探析
    The Book of Songs “Like a boat of cypress wood”
    Analysis on China’s real estate bubble by Oct.2016
    商情(2017年26期)2017-07-28 12:30:30
    程序框圖題盤點
    睡眠酒店
    相親
    故事會(2010年20期)2010-10-11 08:37:24
    91麻豆精品激情在线观看国产| 18禁黄网站禁片午夜丰满| 69av精品久久久久久| 国产单亲对白刺激| 亚洲欧美日韩无卡精品| 尾随美女入室| 亚洲成a人片在线一区二区| 俺也久久电影网| 乱人视频在线观看| 97人妻精品一区二区三区麻豆| 亚洲色图av天堂| 啦啦啦观看免费观看视频高清| 亚洲五月天丁香| 波多野结衣高清作品| 亚洲成人中文字幕在线播放| 欧美日韩黄片免| 国产精品亚洲美女久久久| 亚洲在线观看片| 国模一区二区三区四区视频| 亚洲狠狠婷婷综合久久图片| 精品人妻1区二区| 亚洲欧美精品综合久久99| 午夜老司机福利剧场| 亚洲av免费在线观看| 男人舔奶头视频| 一级黄色大片毛片| 啦啦啦韩国在线观看视频| 夜夜夜夜夜久久久久| 变态另类丝袜制服| 国内精品美女久久久久久| 在线国产一区二区在线| 又爽又黄无遮挡网站| 国产亚洲精品久久久com| 国产av不卡久久| av天堂在线播放| 欧美三级亚洲精品| 国产精品国产高清国产av| 看片在线看免费视频| 亚洲人成网站在线播| 亚洲人成伊人成综合网2020| 在线免费观看不下载黄p国产 | 欧美3d第一页| 国产成人av教育| 非洲黑人性xxxx精品又粗又长| 精品一区二区三区视频在线| 综合色av麻豆| 在线观看一区二区三区| 两个人的视频大全免费| 国产探花极品一区二区| 国产在视频线在精品| 国产精品一区www在线观看 | 亚洲av熟女| 91午夜精品亚洲一区二区三区 | 免费观看人在逋| 日本色播在线视频| 1000部很黄的大片| 日韩精品中文字幕看吧| 直男gayav资源| 欧美色视频一区免费| 精品久久久久久久末码| 91久久精品国产一区二区三区| 99热网站在线观看| 日本与韩国留学比较| 熟女电影av网| av在线亚洲专区| 日本三级黄在线观看| 欧美中文日本在线观看视频| 午夜激情福利司机影院| 最近在线观看免费完整版| 久久久久九九精品影院| 一夜夜www| 在线观看美女被高潮喷水网站| 琪琪午夜伦伦电影理论片6080| 99热只有精品国产| 香蕉av资源在线| 97人妻精品一区二区三区麻豆| 中文字幕av成人在线电影| 成年女人看的毛片在线观看| 一级黄色大片毛片| 伦理电影大哥的女人| 成人永久免费在线观看视频| 欧美最黄视频在线播放免费| 精品午夜福利视频在线观看一区| 99久久精品国产国产毛片| 69av精品久久久久久| 亚洲最大成人av| 天堂动漫精品| 国产精品国产三级国产av玫瑰| 久久人人精品亚洲av| 欧美一级a爱片免费观看看| 日本一本二区三区精品| 一进一出抽搐动态| 最后的刺客免费高清国语| 69人妻影院| 长腿黑丝高跟| 又粗又爽又猛毛片免费看| 国产精品人妻久久久影院| 午夜福利欧美成人| 欧美激情国产日韩精品一区| 真人做人爱边吃奶动态| 久久国产精品人妻蜜桃| 黄色丝袜av网址大全| 精品日产1卡2卡| 99热这里只有是精品在线观看| 熟女人妻精品中文字幕| 日本一二三区视频观看| 日本 欧美在线| 亚洲性夜色夜夜综合| 亚洲18禁久久av| 九九久久精品国产亚洲av麻豆| 国产三级中文精品| 中文亚洲av片在线观看爽| 精品福利观看| 国产精品亚洲美女久久久| 九色国产91popny在线| 天美传媒精品一区二区| 男女做爰动态图高潮gif福利片| 男女下面进入的视频免费午夜| 国产精品福利在线免费观看| 国产精品综合久久久久久久免费| 永久网站在线| 1000部很黄的大片| 麻豆久久精品国产亚洲av| 又黄又爽又免费观看的视频| 欧美+日韩+精品| 午夜久久久久精精品| 婷婷丁香在线五月| 成人亚洲精品av一区二区| 女人被狂操c到高潮| 天堂网av新在线| 亚洲无线观看免费| 国产麻豆成人av免费视频| 免费高清视频大片| 床上黄色一级片| 我要看日韩黄色一级片| 久久天躁狠狠躁夜夜2o2o| 男插女下体视频免费在线播放| 99热这里只有是精品在线观看| av在线老鸭窝| 日日摸夜夜添夜夜添小说| 国产精品人妻久久久久久| 免费av不卡在线播放| 欧美+亚洲+日韩+国产| 国产主播在线观看一区二区| 亚洲熟妇熟女久久| 美女xxoo啪啪120秒动态图| 高清毛片免费观看视频网站| 人妻少妇偷人精品九色| 色吧在线观看| 国产黄色小视频在线观看| 俺也久久电影网| x7x7x7水蜜桃| 久久精品国产亚洲av天美| av专区在线播放| 成人永久免费在线观看视频| 亚洲国产高清在线一区二区三| 亚洲精品一区av在线观看| 久久久久久久久大av| 国产精品精品国产色婷婷| 免费人成视频x8x8入口观看| 欧美成人一区二区免费高清观看| 成人av在线播放网站| 亚洲久久久久久中文字幕| 欧美最黄视频在线播放免费| 婷婷亚洲欧美| 国产乱人视频| 一进一出好大好爽视频| 乱人视频在线观看| 久久久色成人| 亚洲国产精品久久男人天堂| 国产精品野战在线观看| 伦理电影大哥的女人| 亚洲真实伦在线观看| 亚洲av.av天堂| 一夜夜www| 日本撒尿小便嘘嘘汇集6| 成人美女网站在线观看视频| 国内精品一区二区在线观看| 一进一出抽搐gif免费好疼| 亚洲国产精品sss在线观看| 日本精品一区二区三区蜜桃| 欧美激情在线99| 午夜激情福利司机影院| 亚洲欧美日韩高清在线视频| 少妇高潮的动态图| 嫩草影视91久久| 黄色女人牲交| 国产一区二区三区在线臀色熟女| 色播亚洲综合网| aaaaa片日本免费| 国产成人a区在线观看| 国产av麻豆久久久久久久| 午夜福利视频1000在线观看| 国内精品久久久久久久电影| 国产大屁股一区二区在线视频| 国产精品国产高清国产av| 男人舔奶头视频| 天堂动漫精品| 成人无遮挡网站| 老司机午夜福利在线观看视频| 赤兔流量卡办理| 成年女人看的毛片在线观看| 欧美bdsm另类| or卡值多少钱| av女优亚洲男人天堂| 最近最新中文字幕大全电影3| 一个人看视频在线观看www免费| 一区二区三区高清视频在线| 国产主播在线观看一区二区| 一夜夜www| 狂野欧美白嫩少妇大欣赏| 国产成人a区在线观看| 欧美日韩综合久久久久久 | 国产精品免费一区二区三区在线| 欧美日韩乱码在线| 国内精品久久久久久久电影| 可以在线观看毛片的网站| 中文字幕高清在线视频| 看免费成人av毛片| 婷婷精品国产亚洲av在线| 美女黄网站色视频| 97热精品久久久久久| 亚洲久久久久久中文字幕| 内射极品少妇av片p| 免费看a级黄色片| 日本精品一区二区三区蜜桃| 俺也久久电影网| 国产av一区在线观看免费| 网址你懂的国产日韩在线| 久久国产乱子免费精品| 日韩强制内射视频| 乱系列少妇在线播放| 久久久久久久亚洲中文字幕| h日本视频在线播放| 午夜福利欧美成人| 一个人免费在线观看电影| 国产美女午夜福利| 国产高清三级在线| 亚洲av五月六月丁香网| 国产av在哪里看| 国产高清视频在线观看网站| 我要搜黄色片| 麻豆精品久久久久久蜜桃| 中文在线观看免费www的网站| 日韩国内少妇激情av| 亚洲国产高清在线一区二区三| 日本 av在线| 国产黄色小视频在线观看| 麻豆精品久久久久久蜜桃| 少妇的逼水好多| 色哟哟哟哟哟哟| 美女高潮喷水抽搐中文字幕| videossex国产| 黄色丝袜av网址大全| 欧美最新免费一区二区三区| 久久精品国产自在天天线| 精品免费久久久久久久清纯| 国产单亲对白刺激| 韩国av一区二区三区四区| 精品人妻1区二区| av在线天堂中文字幕| 国产精品美女特级片免费视频播放器| 联通29元200g的流量卡| 麻豆久久精品国产亚洲av| 亚洲av成人av| 最后的刺客免费高清国语| 久久人妻av系列| 伊人久久精品亚洲午夜| 熟女电影av网| 男女做爰动态图高潮gif福利片| 美女高潮喷水抽搐中文字幕| 午夜老司机福利剧场| 亚洲天堂国产精品一区在线| 三级男女做爰猛烈吃奶摸视频| 我的女老师完整版在线观看| 久久天躁狠狠躁夜夜2o2o| 精品人妻一区二区三区麻豆 | 校园春色视频在线观看| 欧美激情国产日韩精品一区| 久久国产乱子免费精品| 国内毛片毛片毛片毛片毛片| 久久久久性生活片| 国产伦一二天堂av在线观看| 亚洲乱码一区二区免费版| 制服丝袜大香蕉在线| 啦啦啦啦在线视频资源| 国产精品99久久久久久久久| 亚洲成人精品中文字幕电影| 亚洲av第一区精品v没综合| 日韩欧美在线二视频| 国内精品一区二区在线观看| 午夜福利欧美成人| 天天躁日日操中文字幕| 亚洲不卡免费看| 三级男女做爰猛烈吃奶摸视频| 国产蜜桃级精品一区二区三区| 嫩草影院入口| 国产av一区在线观看免费| 亚洲人成网站在线播放欧美日韩| 久久久久久久久久黄片| a在线观看视频网站| 99久久成人亚洲精品观看| 欧美区成人在线视频| 亚洲av美国av| 国产真实乱freesex| 免费av观看视频| 精品人妻视频免费看| 一个人免费在线观看电影| 精品久久久久久成人av| 婷婷精品国产亚洲av| 久久久久久久久大av| 99在线视频只有这里精品首页| 欧美潮喷喷水| 男女啪啪激烈高潮av片| av在线天堂中文字幕| 欧美成人免费av一区二区三区| 免费无遮挡裸体视频| 国内揄拍国产精品人妻在线| 久久久久九九精品影院| 99久久成人亚洲精品观看| 午夜福利在线观看免费完整高清在 | 1000部很黄的大片| 波多野结衣高清作品| 一进一出抽搐动态| 国产精品亚洲美女久久久| 18禁黄网站禁片免费观看直播| 嫩草影院入口| 天堂网av新在线| 久久久久久九九精品二区国产| 久久久久久大精品| 男人狂女人下面高潮的视频| 天天一区二区日本电影三级| 婷婷六月久久综合丁香| 毛片一级片免费看久久久久 | 欧美最新免费一区二区三区| 国产蜜桃级精品一区二区三区| 国产高清有码在线观看视频| 国产久久久一区二区三区| 999久久久精品免费观看国产| 搡女人真爽免费视频火全软件 | 91狼人影院| 久久久久国产精品人妻aⅴ院| 22中文网久久字幕| 欧美日韩中文字幕国产精品一区二区三区| 一级毛片久久久久久久久女| 两人在一起打扑克的视频| 熟妇人妻久久中文字幕3abv| avwww免费| videossex国产| 天堂网av新在线| 日本五十路高清| 日韩av在线大香蕉| 日韩欧美国产在线观看| 身体一侧抽搐| 亚洲性夜色夜夜综合| 成年女人毛片免费观看观看9| 看十八女毛片水多多多| 国产成年人精品一区二区| 非洲黑人性xxxx精品又粗又长| 韩国av一区二区三区四区| 国内揄拍国产精品人妻在线| 免费观看的影片在线观看| 免费av观看视频| 人人妻人人澡欧美一区二区| 校园人妻丝袜中文字幕| 九色国产91popny在线| 老熟妇仑乱视频hdxx| 观看美女的网站| 久久久久精品国产欧美久久久| av在线蜜桃| 美女xxoo啪啪120秒动态图| 久久精品国产鲁丝片午夜精品 | 极品教师在线免费播放| 欧美极品一区二区三区四区| 在线a可以看的网站| 女生性感内裤真人,穿戴方法视频| 韩国av在线不卡| 免费无遮挡裸体视频| 日本精品一区二区三区蜜桃| 国产成人av教育| 白带黄色成豆腐渣| 成年女人永久免费观看视频| 如何舔出高潮| 三级国产精品欧美在线观看| 一区二区三区四区激情视频 | 长腿黑丝高跟| 男人狂女人下面高潮的视频| 在线天堂最新版资源| 久久久精品欧美日韩精品| 免费黄网站久久成人精品| 悠悠久久av| 亚洲黑人精品在线| 成人欧美大片| 精品久久久久久久久亚洲 | 亚洲一区二区三区色噜噜| 极品教师在线免费播放| 欧美国产日韩亚洲一区| 成人毛片a级毛片在线播放| 国产伦人伦偷精品视频| 成人美女网站在线观看视频| 国产 一区精品| 俺也久久电影网| 尾随美女入室| 国产精品,欧美在线| x7x7x7水蜜桃| 中文亚洲av片在线观看爽| 男女之事视频高清在线观看| 日韩欧美在线二视频| 日韩一本色道免费dvd| 色精品久久人妻99蜜桃| 麻豆一二三区av精品| 男人舔女人下体高潮全视频| 成人性生交大片免费视频hd| 黄色一级大片看看| 波多野结衣高清作品| 别揉我奶头~嗯~啊~动态视频| 国产久久久一区二区三区| 精品人妻一区二区三区麻豆 | 女同久久另类99精品国产91| 亚洲国产高清在线一区二区三| 亚洲国产高清在线一区二区三| 午夜a级毛片| 在线观看午夜福利视频| 久久久久久久久中文| 亚洲自拍偷在线| 美女cb高潮喷水在线观看| 国产免费男女视频| 日日撸夜夜添| 美女高潮的动态| 又紧又爽又黄一区二区| 国产一区二区三区在线臀色熟女| 色尼玛亚洲综合影院| 国产高清视频在线播放一区| 69av精品久久久久久| 免费人成在线观看视频色| 久久精品夜夜夜夜夜久久蜜豆| 可以在线观看的亚洲视频| av福利片在线观看| 国产精品爽爽va在线观看网站| 欧美国产日韩亚洲一区| 中文字幕av成人在线电影| 亚洲精品影视一区二区三区av| 小蜜桃在线观看免费完整版高清| 亚洲美女视频黄频| 97超级碰碰碰精品色视频在线观看| 一个人看视频在线观看www免费| 亚洲18禁久久av| 国产精品久久久久久久久免| 日本成人三级电影网站| 久久人人爽人人爽人人片va| 成年女人看的毛片在线观看| 精品乱码久久久久久99久播| 久久99热6这里只有精品| 女同久久另类99精品国产91| 欧美精品国产亚洲| 欧美zozozo另类| 免费一级毛片在线播放高清视频| 99热网站在线观看| 成人二区视频| 成人特级av手机在线观看| 日本免费a在线| 国内久久婷婷六月综合欲色啪| 我的老师免费观看完整版| 久久亚洲真实| 色5月婷婷丁香| 97超视频在线观看视频| 一级av片app| 久久久久久九九精品二区国产| 亚洲国产精品成人综合色| 精品久久久久久久人妻蜜臀av| 国产精品免费一区二区三区在线| 国产v大片淫在线免费观看| 欧美又色又爽又黄视频| 亚洲四区av| 欧美日韩综合久久久久久 | 女的被弄到高潮叫床怎么办 | 亚洲人成伊人成综合网2020| 999久久久精品免费观看国产| 国产伦人伦偷精品视频| 国产亚洲精品av在线| 亚洲内射少妇av| 国产色婷婷99| 国产精品久久久久久精品电影| 日本黄色片子视频| 波多野结衣高清无吗| 成年人黄色毛片网站| 亚洲人成网站在线播放欧美日韩| 夜夜爽天天搞| 亚洲,欧美,日韩| 国国产精品蜜臀av免费| 免费在线观看日本一区| 99热网站在线观看| 久久久久精品国产欧美久久久| 日韩精品有码人妻一区| 精品午夜福利视频在线观看一区| 亚洲乱码一区二区免费版| 免费在线观看成人毛片| 男女啪啪激烈高潮av片| 别揉我奶头~嗯~啊~动态视频| 国产一区二区三区视频了| 一区二区三区免费毛片| 国产午夜精品论理片| 俄罗斯特黄特色一大片| 欧美高清成人免费视频www| 国产一区二区在线av高清观看| 久久久久久久亚洲中文字幕| 91麻豆精品激情在线观看国产| 欧美高清性xxxxhd video| eeuss影院久久| x7x7x7水蜜桃| 久久久久久九九精品二区国产| 人妻少妇偷人精品九色| 成人av一区二区三区在线看| 一本久久中文字幕| 麻豆成人av在线观看| 久久午夜福利片| 欧美最新免费一区二区三区| 中亚洲国语对白在线视频| 99久久九九国产精品国产免费| 91在线精品国自产拍蜜月| 亚洲精华国产精华液的使用体验 | 国产免费男女视频| 亚洲av一区综合| 日本熟妇午夜| 日本a在线网址| 免费看光身美女| 国内毛片毛片毛片毛片毛片| 国产伦一二天堂av在线观看| 18禁在线播放成人免费| 尾随美女入室| 免费看美女性在线毛片视频| 在线观看舔阴道视频| 深爱激情五月婷婷| 亚洲,欧美,日韩| 搡老妇女老女人老熟妇| 国产免费男女视频| 搡老妇女老女人老熟妇| 国产伦一二天堂av在线观看| 日本免费a在线| 日本熟妇午夜| av视频在线观看入口| 老熟妇乱子伦视频在线观看| 欧美激情在线99| 白带黄色成豆腐渣| 亚洲精华国产精华精| 欧美最新免费一区二区三区| 中国美女看黄片| 在线天堂最新版资源| 日本欧美国产在线视频| 国产精品久久久久久久电影| 日韩高清综合在线| 精品不卡国产一区二区三区| 亚洲黑人精品在线| 免费av不卡在线播放| 欧美成人a在线观看| 色综合色国产| 久久午夜福利片| 在线观看免费视频日本深夜| 性色avwww在线观看| 亚洲成人免费电影在线观看| 日本一二三区视频观看| 永久网站在线| 亚洲四区av| 国产女主播在线喷水免费视频网站 | a级一级毛片免费在线观看| 国产一区二区激情短视频| 欧美一级a爱片免费观看看| 国产精品国产三级国产av玫瑰| 观看免费一级毛片| 88av欧美| 成人综合一区亚洲| 国产国拍精品亚洲av在线观看| 成人一区二区视频在线观看| 久久精品综合一区二区三区| a级一级毛片免费在线观看| 久久久久久大精品| 欧美高清成人免费视频www| 偷拍熟女少妇极品色| 熟妇人妻久久中文字幕3abv| 国产精品自产拍在线观看55亚洲| 色5月婷婷丁香| 亚洲七黄色美女视频| 国产精品一区二区三区四区免费观看 | 大型黄色视频在线免费观看| 一进一出抽搐动态| 免费看a级黄色片| 又爽又黄无遮挡网站| 亚洲欧美日韩无卡精品| 欧美日韩亚洲国产一区二区在线观看| 国产蜜桃级精品一区二区三区| 黄色丝袜av网址大全| 亚洲精品色激情综合| 黄色欧美视频在线观看| 久久久久久国产a免费观看| 亚洲中文日韩欧美视频| 99九九线精品视频在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 岛国在线免费视频观看| 婷婷精品国产亚洲av| 久久久久九九精品影院| 国内精品美女久久久久久| 日韩强制内射视频| 精品久久久久久久人妻蜜臀av| 久久久久久久午夜电影| 欧美人与善性xxx| 嫩草影院入口| 亚洲成人免费电影在线观看| 国产欧美日韩一区二区精品| 国产爱豆传媒在线观看| 两人在一起打扑克的视频| 亚洲精品成人久久久久久| 黄色丝袜av网址大全| 国产亚洲91精品色在线| 成人精品一区二区免费| 嫩草影院新地址|