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

    A study of grid failure mode drivers and methods for accelerated life testing of a 30cm diameter ion thruster

    2023-11-19 13:32:42MingmingSUN孫明明JianfeiLONG龍建飛JuanjuanCHEN陳娟娟WeiYANG楊威WeilongGUO郭偉龍andXinweiCHEN陳新偉
    Plasma Science and Technology 2023年11期
    關(guān)鍵詞:楊威

    Mingming SUN (孫明明) ,Jianfei LONG (龍建飛) ,Juanjuan CHEN (陳娟娟) ,Wei YANG (楊威) ,Weilong GUO (郭偉龍) and Xinwei CHEN (陳新偉)

    1 Science and Technology on Vacuum Technology and Physics Laboratory,Lanzhou Institute of Physics,Lanzhou 730000,People’s Republic of China

    2 School of Fundamental Physics and Mathematical Sciences,Hangzhou Institute for Advanced Study,University of Chinese Academy of Sciences,Hangzhou 310024,People’s Republic of China

    3 Gravitational Wave Universe Taiji Laboratory,Hangzhou Institute for Advanced Study,University of Chinese Academy of Sciences,Hangzhou 310024,People’s Republic of China

    4 Key Laboratory of Gravitational Wave Precision Measurement of Zhejiang Province,University of Chinese Academy of Sciences,Hangzhou 310024,People’s Republic of China

    5 School of Electrical Engineering,Chongqing University of Science and Technology,Chongqing 404100,People’s Republic of China

    Abstract In view of the high cost caused by the 1:1 lifetime verification test of ion thrusters,the lifetime acceleration test should be considered.This work uses the PIC-MCC (Particle-in-Cell Monte-Carlo Collision) method to analyze the five failure factors that lead to the failure of the accelerator grid of a 30 cm diameter ion thruster under the working mode of 5 kW.Meanwhile,the acceleration stress levels corresponding to different failure factors are obtained.The results show that background pressure has the highest stress level on the grid’s erosion.The accelerator grid aperture’s mass sputtering rate under the rated vacuum degree(1×10-4 Pa)of 5 kW work mode is 8.78 times that of the baseline vacuum degree (1×10-6 Pa),and the mass sputtering rate under worse vacuum degree (5×10-3 Pa) is 5.08 times that of 1×10-4 Pa.Under the influence of the other four failure factors,namely,the voltage of the accelerator grid,upstream plasma density,the screen grid voltage and mass utilization efficiency,the mass sputtering rates of the accelerator grid hole are 2.32,2.67,1.98 and 2.51 times those of the accelerator grid hole under baseline condition,respectively.The ion sputtering results of two 30 cm diameter ion thrusters (both installed with new grids assembly) after working for 1000 h show that the mass sputtering rate of the accelerator grid hole under vacuum conditions of 5×10-3 Pa is 4.54 times that under the condition of 1×10-4 Pa,and the comparison error between simulation results and test results of acceleration stress is about 10%.In the subsequent ion thruster lifetime verification,the working vacuum degree can be adjusted according to the acceleration stress level of background pressure,so as to shorten the test time and reduce the test cost.

    Keywords: ion thruster,failure factors,acceleration stress level

    1.Introduction

    Ion thrusters have been widely used in the field of spacecraft propulsion based on the characteristics of large specific impulse,long life and high efficiency.With the expansion of applications,the working life verification of ion thrusters has become the focus of research [1-3].Taking the current communication,meteorological,reconnaissance and other satellites as examples,the ion thrusters are required to be in orbit for more than 15 years,with a total cumulative operating life of more than 30000 h [4,5].The 1:1 lifetime verification takes at least 5 years to complete the ground demonstration test.Therefore,in terms of financial resources and human resources,a 1:1 lifetime verification test is a huge cost.Moreover,it also seriously restricts the development cycle of ion thrusters.Therefore,the current lifetime verification test cannot meet the needs of the rapid development of spacecraft equipment,and it is necessary to study a new,fast and efficient lifetime evaluation method.The lifetime acceleration test is an effective way to meet the long-life requirements of ion thrusters.

    Over the past decade,many studies focused on the sputtering mechanism of the grids and paid attention to the sputtering on the accelerator grid by CEX (charge exchange)ions.Brophy et al [6] studied the failure mechanism and erosion characteristics of a two-grid system of a 30 cm diameter ion thruster.The research results indicated that there are three failure mechanisms to the grids,which are grid shorting failure due to flaking deposited on the screen grid,electron backstreaming failure,and structural failure by pits-andgrooves erosion,respectively.The vacuum chamber pressure increases the flatness of the accelerator grid mass loss,and the accelerator grid current increased from 30 to 42 mA as the tank pressure increased from 2.4×10-5to 3.4×10-5Torr.Brophy et al further indicated that the structural failure is not the first failure mechanism for the accelerator grid.Nakano[7] studied the relationship between the electron backstreaming failure of a three-grid system and the extraction current of a grid single hole through modeling,and used the model to simulate and predict the sputtering shape of the accelerator grid aperture when electron backstreaming occurred.Xia et al [8] used the PIC-MCC method to model and study the beam-focusing difference between a two-grid system and a three-grid system.It indicates that the sputtering of CEX ions on the accelerator grid hole is the main cause of the electron backstreaming and eventually the failure of a three-grid system.Zhang et al [9] established a simulation model of a three-grid system,and studied the influence of background pressure on the ion energy distribution and ion current collected at different locations of the accelerator grid and the decelerator grid.Sengupta et al [10] tested and studied the sensitives of the flow rate of the anode,the hollow cathode and the neutralizer of the NSTAR ion thruster to the effects of discharge voltage and beam current,but did not study the effects of working parameters on specific failure modes.Chen et al [11] developed a three-dimensional numerical model by PIC-MCC method to simulate the beam ion acceleration,extraction,CEX generation and hitting processes of a three-grid system.In 2008,Noord et al [12] tested the influence of parameters adjustment on grid wear of NEXT ion thruster.The results showed that when the screen grid voltage was increased from 1021 to 1179 V,the sputtering rate of the accelerator grid apertures reached the maximum.Based on this mass sputtering rate,due to the continuous expansion of the accelerator grid apertures,the negative potential barrier in the center of the grid aperture to prevent the electrons from the plume near-field region to the discharge chamber,will gradually decrease,and eventually lead to the electron backstreaming failure of the accelerator grid.These previous studies have well demonstrated that the grid,especially the accelerator grid,is the core component that affects the life of the ion thruster,and the parameter adjustment test results prove that there are failure factors on the grid life and different acceleration stress levels on the grid lifetime.However,most of the studies have not fully analyzed the multiple factors that cause the failure of the accelerator grid and the stress level of various factors on the erosion of the accelerator grid.

    Based on the fact that the accelerator grid is the critical component for a three-grid system,then the key to designing the lifetime acceleration test of the ion thruster is to study the effect of different failure factors on the erosion of the grids,and determine the most important influence factor and its acceleration stress level on the lifetime of the grids.In this work,the PIC-MCC method is used to analyze the effects of different failure factors on a three-grid system of a 30 cm diameter ion thruster(developed by LIP,Lanzhou Institute of Physics)under 5 kW work mode,especially ion sputtering of the accelerator grid.The most important influence factor and its acceleration stress level on the erosion of the grids is determined by calculation,thus providing a reference for the subsequent design of the lifetime acceleration test of ion thrusters,so as to achieve the purpose of greatly reducing the test verification period and the test cost.

    2.Failure factors and setting range

    The purpose of the lifetime acceleration test is to compress the test time from 5 years to 1 year or less without affecting the verification life cycle.The design steps are as follows.Firstly,according to the failure modes of the grids of the thruster,the multiple failure factors are sorted out,and the most important influence factor is identified.Secondly,on the basis of obtaining the most important influence factor,the acceleration stress level is calculated by numerical simulation.According to previous research results[3-7],the main failure modes of the accelerator grid are electron backstreaming failure,structural failure,and shorting failure due to screen grid erosion.The screen grid erosion is not obvious in the lifetime test of ion thrusters.The life test of the 30 cm diameter ion thruster conducted in LIP showed a slight change in the aperture and thickness of the screen grid after 6500 h of operation[13].The structural failure of the accelerator grid is mainly due to the pits-and-grooves erosion caused by CEX ions sputtering,and the pits-and-grooves erosion around individual apertures will eventually wear through the accelerator grid,resulting in short with the screen grid by thin metal sputtering.However,for a three-grid system,the structural failure is not the first failure mechanism for the accelerator grid [6],and the decelerator grid effectively reduces the pits-and-grooves corrosion caused by CEX ions on the accelerator grid,thus greatly reducing the risk of structural failure of the accelerator grid[7,8].Meanwhile,the accelerator grid hole is enlarged due to CEX ions sputtering,eventually leading to potentials in the center of the grid apertures that are not sufficiently negative to prevent electron backstreaming into the discharge chamber.Therefore,electron backstreaming failure is the main failure mode for a three-grid system,and in this work,only the failure factors that can cause electron backstreaming failure are studied.In addition,it should be noted that although the performance of the thruster will change during the lifetime acceleration test,mechanisms of the different failure modes will not change except to accelerate the sputtering of ions and the failure of the grids,so there are no new failure modes introduced.Meanwhile,the acceleration test should select reasonable test conditions and duration according to the stress level so as to avoid the abnormal operation of thrusters for a long time.

    There are many factors that can cause electron backstreaming failure of the accelerator grid[6-14].At present,it can be summarized that the test and the analysis results mainly include background pressure,the diameter of the grid hole,the voltage of the accelerator grid,upstream plasma density,the screen grid voltage,propellant flow rate,mass utilization efficiency,thickness of the accelerator grid and the screen grid,etc.In view of these failure factors,it is necessary to determine the research prerequisite,and carry out equivalent replacement of the failure factors with coupling relationship.First of all,the thruster design is not changed in the analysis,that is,the parameters changed by subjective factors are only used as input boundary conditions.For example,the grid gap is only considered as an input parameter rather than as a failure factor,since the equilibrium hot gap under different powers is a fixed value[14,15].Other parameters,such as the diameter of the grid hole,electric field,and magnetic field,are also used as input boundaries.Secondly,the ion beam anomalies caused by the design defects of grid geometric parameters are not considered,and only the failure caused by normal ion erosion is targeted.Finally,parameters with coupling relationships can be replaced with a single parameter.For example,the plasma density upstream of the grids is affected by the mass flow rate.According to the previous calculation results of the influence of the mass flow rate on the distribution of plasma density [16],the grid erosion caused by the mass flow rate can be converted into the grid erosion caused by upstream plasma density,so as to reduce the number of failure factors.

    After analysis and replacement,there are only five failure factors leading to the electron backstreaming failure,which are background pressure,the screen grid voltage,the accelerator grid voltage,upstream plasma density and mass utilization efficiency (referring to that of the discharge chamber),and these parameters have no coupling relationship with each other and are all single factors.Meanwhile,these factors allhave significant effects on the generation of CEX ions.For example,according to the Child-Langmuir equation,the beam voltage,or the screen grid voltage,can change the average beam density,thus affecting the CEX ion density and changing the sputtering rate of the accelerator grid.The variation of the accelerator grid voltage will change the attraction effect of slow CEX ions,thus affecting the sputtering rate.The extraction beam is proportional to the upstream plasma density,and due to the space charge effect,the upstream plasma density will change the focus of the beam and affect the density and distribution of CEX ions.Increasing the background pressure will enhance the density of residual atoms and the probability of charge exchange collision.The mass utilization efficiency affects the amount of xenon atoms escaping from the discharge chamber through the grid holes.With the decrease of the mass utilization efficiency,the production rate of CEX ion will increase accordingly,so its influence mechanism is similar to the background pressure.

    On these bases,the PIC-MCC method is used to calculate the acceleration stress levels of the five failure factors and the corresponding influence on the accelerator grid corrosion.Table 1 shows the rated working parameters of the 30 cm diameter ion thruster under 5 kW mode.In the table,the hot gap between the screen grid and the accelerator grid is the measurement results when the thruster reaches thermal balance [17].Since the grid gap can be stabilized within 10 min after the thruster is started,the grid gap is taken as a fixed boundary in the calculation.

    The setting points should be selected for different failure factors,and a baseline value is set for comparison of the influence of failure factors.The baseline values are generally selected for conditions that cause low ion sputtering to the grids,and the in-orbit operating conditions can be selected.Then,the rated working condition of the 5 kW mode is compared with it,and the corresponding stress level is calculated as the evaluation basis.In calculation,when one failure factor is changed,the other failure factors remain unchanged.The setting for the failure factors is shown in table 2.It should be noted that the setting points of some failure factors shown in table 2 are far beyond the adjustment range in the actual performance test(the flow rate and voltage adjustment range generally do not exceed ±10% of the rated condition).However,in order to more obviously compare the stress levels of different failure factors,a wider range of parameters is set.

    Table 2. Setting of the failure factors.

    3.Calculation model

    Because the grid apertures are axisymmetric,the twodimensional axisymmetric model can be used to simulate the three-dimensional process of beam extraction,which is shown in figure 1(a).The ions and electrons in the grid acceleration region are treated as charged particles and fluids,respectively.The neutral atoms are also treated as particles.The particle density corresponding to the vacuum degree of the facility is treated as the background neutral density.The PIC method is used to describe the motion of beam ions and CEX ions in the calculation region.The MCC method is used to solve the exchange charge collision process between high-speed beam ions and slow neutral atoms,that is,the generation of CEX ions.As is shown in figure 1(a),the upper,lower and right boundaries of the model are set as Neumann boundaries,that is,=0.It is noted that the left boundary is set as Φ=Vsc+Vp,which means the voltage on the left boundary is the sum of the screen grid voltage and the plasma potential.As is shown in figure 1(b),the upstream of the model is set as the plasma inlet,and the downstream of the model and the surfaces of the three grids are all set as absorb boundary,that is,when the ion moves to these boundaries,the ion information is deleted,and the ion movement is no longer calculated and tracked.The upper and lower boundaries are set as reflect boundary,that is,when the ion moves to these boundaries,the ion is treated by specular reflection and return to the calculation area and re-participate in the calculation.In addition,the upstream ion density shown in table 1 is on the order of 1017m-3,in order to reduce the computational burden of PIC simulation,the number of macroscopic ions ineach time step is entered into the calculation area according to an appropriate ratio of 10000:1 to the calculated ion number.

    Table 3. Parameters setting of the calculation model.

    Table 3 shows the parameters setting of the calculation model,and the simulated work mode is 5 kW.Wherersc,racandrdelare the radii of the screen grid,the accelerator grid and the decelerator grid respectively,tsc,tacandtdelare the thicknesses of the screen grid,the accelerator grid and the decelerator grid,respectively.ds-aandda-dare the gaps between the screen grid and the accelerator grid,the accelerator grid and the deceleration grid,respectively (the values in the table are measured results).Vsc,VaccandVdelare the voltages of the screen grid,the accelerator grid and the decelerator grid,respectively.Tiis ion temperature which is approximated with the temperature of the discharge chamber[18],Teuis electron temperature upstream of the grids,and the electrons in this region are mainly from ionization in the discharge chamber and emission of the hollow cathode.Tedis electron temperature downstream of the grids,and the electrons in this region are mainly from ionization in the plume region and emission of the neutralizer.Vpis the plasma potential,which is set according to the calculation result of the plasma sheath potential [19].

    By ignoring the effect of magnetic fields on the ion beam,the potential can be solved by the Poisson equation,as shown in equation (1),whereε0is the dielectric constant,niandneare ion density and electron density,respectively.

    For a two-dimensional axisymmetric model,the Poisson equation can be expressed as equation (2),whererandzare radial and axial coordinates of the mesh grid.When solving the Poisson equation,the whole calculation region is divided into a large number of equally meshed grids,the size of which is smaller than Debye lengthλd.Therefore,the size of the meshed grid is set as 0.025 mm and the time step is set as 1.5×10-10s in the model after calculation (the minimum Debye length is calculated to be 0.029 mm).The finite difference method and Gauss-Seidel iteration are used to solve the Poisson equation to obtain the potential distribution in the calculation region,and the relaxation factor is introduced to optimize the calculation.The electric field intensity at the ion location is solved by linear interpolation.

    Figure 1.Solution region of the model.(a) Simulation model and (b) boundary setting.

    According to the Newton-Lorentz equation for ions,which is given by

    whereMis ion mass,viandxiare speed and position of ions,respectively,Eis electric field intensity at the location of the ions.In the model,it is assumed that xenon atom velocity follows the Maxwell distribution,and the atom enters the calculation region from the upstream boundary at thermal velocity.The position and velocity of the atom are updated according to the defined number of cycles.Electrons in the acceleration region are treated as a fluid,whose density follows the Boltzmann distribution,which is shown in equation (4).Whereneis the electron density,?is the potential at the location of electrons,ne,ref,Te,refand?refare the electron density and electron temperature and potential of reference points,respectively.In the calculation of electron density upstream and downstream of the accelerator grid,the reference points are selected in the interior of the discharge chamber and in the neutralization area of the plume [8].

    The double charged ions are not considered in the calculation,and the number of ions entering the calculation area from the upstream boundary (which is shown in figure 1) in the unit time step can be expressed as equation(5),whereTeis the electron temperature.It is noted that on the basis of obtaining the time step of each ion and the total time for the beam to reach stability,in each time step,a certain number of ions,ΔN,which enter the calculation region from the left boundary,are accelerated by the grid electric field and extracted from the right boundary.Then,the number of ions extracted from the right boundary and the corresponding calculated value of the beam current is counted in each time step.When the calculated value of the beam current is equal to the designed value of the beam current,it is considered that the extracted beam current reaches the steady state,and the ion density distribution is the result after the time average.

    The initial velocity of ions entering the calculation region is set as Bohm velocity,and the initial position of ions entering the calculation region within each time step is random.As is shown in figure 1(a),the initial position of the ion in directionz(axial direction) is 0,the radial direction (rdirection) and the incidence angle of the ion are random numbers.Therefore,the initial position of the ion entering the calculation region can be expressed asr0=,whereranis a random number (generally in the range of 0-1).

    In the MCC collision model,elastic collisions between ions and atoms are considered to simulate the charge exchange process.In the unit time stepΔt,the collision probabilityPcan be expressed as equation (6),wherexiandviare the position and velocity of the ion,σis the charge exchange cross section [20],which can be expressed asσ=(k1ln Δv+k2)2×10-20(wherek1andk2are all constant,andΔvis the velocity of ions relative to atoms at the time of collision),nis the atomic density at the location of the collision.

    After ignoring the charge exchange process between other neutral particles in the vacuum chamber,such as O2,N2and H2O,and further ignoring the charge exchange process between the Mo atom of the grid and xenon ions [21],it can be considered that the main source of CEX ion is residual xenon atoms in the vacuum chamber and unionized xenon atoms permeating from the discharge chamber to the grids.With the assumption that the residual xenon atomic density is uniformly distributed near the grids,and the xenon atomic number density follows the ideal gas equation,that is,the residual xenon atomic densityn1can be expressed asn1=P/kT(P,Tandkare the pressure,gas temperature and Boltzmann constant,respectively).In the calculation,the xenon gas temperatures in the upstream and downstream of the grids (shown in figure 1(b)) are set to 500 K and 300 K,respectively.The density of unionized xenon atoms permeating into the three grids from the discharge chamber,namelyn2,can be calculated according to the mass utilization efficiency and the number of ions entering the calculation region in each time step by the PIC method.The atomic densityn(xi)shown in equation (6) can be expressed as the sum ofn1andn2,which is shown in equation (7).

    Figure 2.Density and location of CEX ions under different background pressures.(a)1×10-6 Pa,(b)1×10-4 Pa,(c)5×10-3 Pa and(d)1 × 10-2 Pa.

    In the calculation,the number of xenon ions bombarded to the surface of the accelerator grid hole and the decelerator grid hole is counted,respectively,in each time step,and the sputtering yield of Mo material corresponding to the number of the counted ions is calculated.That is,when a single ion with a particular energy bombards the grid surface,it will produce a certain number of atoms or sputtering,and the unit is n/ ion or g/ ion.The amount of sputtering yield is calculated by the preset number of iterations,and sputtering yield is calculated according to the fitting formula given by Rosenberg,Brophy and Yim et al [22-24].

    4.Calculation results and analysis

    4.1.Influence of background pressure

    The background gas in the vacuum facility is mainly composed of xenon.Because the ion beam contains a large number of high-energy particles,gas discharge can still occur in the near and far fields of the plume,resulting in a large number of CEX ions and etching the grid.In addition,a high vacuum degree will cause the neutral gas to re-flux into the discharge chamber to participate in the discharge,and the ratio of electron re-flux to ion beam is required to be controlled within 1% [25].In this work,a total of eight kinds of background pressure states of 1 × 10-6Pa (simulating inorbit environment),1×10-5Pa,1×10-4Pa,5×10-4Pa,1×10-3Pa,5×10-3Pa,1×10-3Pa and 5×10-2Pa are calculated,respectively.The distributions of beam ions,neutral atoms and CEX ions and the sputtering rate of the grid hole are calculated.As is shown in figure 2,due to the emphasis on ion erosion,only the position and density distribution of CEX ions under background pressure of 1×10-6Pa(baseline),1×10-4Pa (5 kW rated condition),5×10-3Pa(common condition)and 1×10-2Pa(ultimate condition)are given in this paper,and density and position distribution of beam ions are not given.

    Figure 3.Influence of background pressure on the sputtering rate of the grid hole.

    According to the simulation results,with the increase of background pressure,the distribution of beam ions is almost not affected,but the density of CEX ion increases obviously.As is shown in figure 2,the CEX ion density between the screen grid and the accelerator is almost unaffected by background pressure,mainly because most of the CEX ions generated in the region between the screen grid and the accelerator grid are extracted to the region between the accelerator grid and the decelerator grid by the influence of the electric field.Since the main source of CEX ions is residual xenon atoms and unionized xenon atoms,and with the increase of background pressure,the residual atoms near the grids will increase significantly,and the probability of charge exchange collision will also increase.Therefore,the CEX ion density increases most significantly in the downstream regions of the grids.As is shown in figure 2,CEX ions bombarded the surfaces of the accelerator grid and the decelerator grid are mainly generated in the region close to the downstream of the calculation region.The CEX ions generated in this region are also the main source of ion sputtering erosion to the accelerator and decelerator grids.

    In the calculation,the mass sputtering rate of CEX ions bombarded to the accelerator grid and the decelerator grid is also counted.Figure 3 shows the relationship between the mass sputtering rate of aperture walls and the background pressure.It is calculated that when the pressures are 5×10-3Pa and 1 × 10-4Pa (5 kW rated condition),the accelerator grid hole sputtering rates are 44.56 times and 8.78 times of that under 1×10-6Pa (baseline),respectively,and the accelerator grid hole sputtering rate under 5×10-3Pa is 5.08 times of that under 1×10-4Pa.

    Figure 4.Density and location of CEX ions under different accelerator grid voltages.(a)-100 V,(b)-300 V,(c)-500 V and(d)-700 V.

    Figure 5.Influence of accelerator grid voltage on the sputtering rate of the grid hole.

    4.2.Influence of the accelerator grid voltage

    When other parameters in table 1 remain unchanged and the accelerator grid voltages are adjusted to -100 V,-300 V,-500 V and -700 V,respectively,the simulation results of position and density distribution of CEX ions are shown in figures 4(a)-(d).As is shown in the figure,with the increase of the accelerator grid voltage,the density and location of CEX ions do not change significantly.The density of CEX ions is still in the order of 1013m-3,and most of the CEX ions that can bombard the grids are mainly concentrated in the region between the accelerator grid and the decelerator grid,and the downstream of the decelerator grid.In addition,according to the distribution of beam ions (not given in the paper),the focusing characteristics of beam ions are normal,and no direct sputtering of beam ions on the grid surface occurs.

    According to the statistical results of the mass sputtering rate of CEX ions bombarding the accelerator grid and the decelerator grid,figure 5 shows that the mass sputtering rate of the aperture of the accelerator grid changes linearly with the increase of the accelerator grid voltage,while the sputtering rate of the decelerator grid hole is almost not affected by the voltage change.When the accelerator grid voltage is -400 V,the sputtering rate of the accelerator grid hole is 2.72×10-15kg s-1,which is 2.32 times the sputtering rate of the baseline voltage of -100 V (1.17×10-15kg s-1).

    4.3.Influence of upstream plasma density

    The upstream plasma density has a significant influence on the focusing characteristics of the ion beam [25].When changing this boundary,it is necessary to investigate the change of ion beam focusing.Therefore,under the condition that other parameters in table 1 remain unchanged,the ion beam focusing and location distribution of CEX ions are given when the upstream plasma density of the screen grid is changed from 1×1017m-3to 8×1017m-3,as is shown in figures 6(a)-(h).

    Figure 6(a)shows that when the upstream plasma density decreases to 1×1017m-3,the ion beam is over-focused and the intercepted current of the decelerator grid is 3.6×10-7A,and beam ions will directly bombard the decelerator grid apertures.Figure 6(h) shows that when the upstream plasma density increases to 8 × 1017m-3,the ion beam is underfocused and beam ions will directly bombard the surface of the accelerator grid.Therefore,when the upstream plasma density is lower than 1 × 1017m-3or higher than 8 × 1017m-3,it is not a normal ion sputtering process.In this paper,the acceleration effect of upstream plasma density on ion sputtering is mainly considered in the normal focusing,that is,both under-focusing and over-focusing will not occur.Figure 7 shows the relationship between the mass sputtering rate and upstream plasma density.The calculation results show that when the upstream plasma density is 6×1017m-3(5 kW rated condition),the sputtering rate of the accelerator grid aperture is 5.74×10-15kg s-1,which is 2.67 times the sputtering rate of the baseline condition of 2×1017m-3(2.15 × 10-15kg s-1).

    4.4.Influence of the screen grid voltage

    Figure 6.Beam focusing and location of CEX ions under different upstream plasma densities.(a) 1.0×1017 m-3,(b) 2.0×1017 m-3,(c)3.0×1017 m-3,(d) 4.0×1017 m-3,(e) 5.0×1017 m-3,(f) 6.0 × 1017 m-3,(g) 7.0×1017 m-3 and (h) 8.0×1017 m-3.

    Figure 7.Influence of upstream plasma density on sputtering rate of the grid hole.

    The perveance of the grids is proportional to the 3/2 power of the screen grid voltage [25].Therefore,the variation of the screen grid voltage will change the average beam density,and then further affect CEX ion density and the mass sputtering rate of the accelerator grid.Figures 8(a)-(d) show ion beam focusing and location of CEX ions when the screen grid voltage changes from 1200 to 2000 V while other parameters in table 1 remain unchanged.It can be concluded that with the increase of the screen grid voltage,the divergence angle of the ion beam increases gradually.Figure 9 shows that when the screen grid voltage is 2000 V (baseline),the mass sputtering rate of the accelerator grid is 1.98 times that of 1200 V(5 kW rated condition).

    4.5.Influence of mass utilization efficiency

    With keeping other parameters in table 1 unchanged,figures 10(a)-(d)show the position and density distribution of CEX ions when the mass utilization efficiency decreases from 90% to 60%.According to the analysis results in figure 10,the beam focusing characteristics do not change significantly,while the density of CEX ions gradually increased with the decrease of mass utilization efficiency.Figure 11 shows that for the accelerator grid erosion,the mass sputtering rate at 70% mass utilization efficiency(5 kW rated condition)is 2.51 times higher than that at 90% mass utilization efficiency(baseline).

    4.6.Acceleration stress level and test verification

    Figure 8.Beam focusing and location of CEX ions under different screen grid voltages.(a)1200 V,(b)1400 V,(c)1800 V and(d)2000 V.

    Figure 9. Influence of the screen grid voltage on the sputtering rate of the grid hole.

    Figure 10.Beam focusing and location of CEX ions under different mass utilization efficiencies.(a) 90%,(b) 80%,(c) 70% and (d) 60%.

    According to the above simulation results,the influences of different failure factors on ion sputtering of the grids are summarized.The stress levels of different failure factors of the accelerator grid are shown in table 4.It can be concluded from table 4 that the acceleration stress level is different,but the background pressure has the highest stress level on the accelerator grid erosion.Therefore,background pressure is the most important influence factor.In addition,considering the adjustment of background pressure in the test is convenient,it is recommended that the background pressure can be taken as an influence factor in the lifetime acceleration test first,and other influence factors can be selected according to the test conditions.

    At present,the special lifetime acceleration verification of adjusting the background pressure in the vacuum chamber has not been carried out,but the simulation results can be verified through the performance test and life test of a 30 cm diameter ion thruster in different vacuum facilities.From 2017 to 2021,two 30 cm diameter ion thrusters(number M04 and M06,with newly installed grids) were,respectively,subjected to 1500 h and 6500 h verification tests in the vacuum facility numbered TS-7 and TS-7A (as shown in figure 12)of LIP[13,19].The TS-7 vacuum facility is 3.8 m in diameter and 8 m in length,which can provide a vacuum level of 5 × 10-3Pa when the thruster is under 5 kW rated condition.The TS-7A vacuum facility has a sub-chamber of 2 m in diameter and 2 m in length,and the main chamber is 4.5 m in diameter and 10 m in length.There are six cryogenic pumps and the total pumping speed is about 260 kl s-1,which can provide a vacuum level close to 1 × 10-4Pa when the thruster is under 5 kW rated condition.

    Figure 11. Influence of mass utilization efficiency on the sputtering rate of the grid hole.

    Figure 12.Test facility of 30 cm diameter ion thruster.(a) TS-7 vacuum facility and (b) TS-7A vacuum facility.

    Figure 13.Diameter variation of the accelerator grid hole under different vacuum degrees after 1000 h.(a) Under 1 × 10-4 Pa,(b) under 5 × 10-3 Pa and (c) comparison of variation.

    Since the grids of the two thrusters are all newly installed grids,the variation of the accelerator grid aperture within 1000 h is taken to verify the simulation results.The erosions of the accelerator grid apertures in TS-7 and TS-7A within 1000 h test are,respectively,shown in figures 13(a) and (b),and the comparison results of the two diameters’variation are shown in figure 13(c).Equation (7) converts the aperture variation into mass sputtering rate,wherevSis the mass sputtering rate,ρis the material density of Mo,r0andr1are the initial value of the radius of the accelerator grid and the measured radius after 1000 h,respectively.According to equation (7),the mass sputtering rates of the accelerator grid holes at 1×10-4Pa and 5×10-3Pa are 1.97×10-14kg s-1and 8.94×10-14kg s-1,respectively.Therefore,the acceleration stress level of the vacuum degree of 5×10-3Pa is 4.54 times that of the vacuum degree of 1×10-4Pa,which is consistent with the simulation result of 5.08 times.

    5.Conclusions

    In this work,the acceleration stress levels of five different failure factors are analyzed under 5 kW operation mode for a 30 cm diameter ion thruster.The simulation results indicate that background pressure is the most important influence factor and has the highest stress level for the mass sputtering rate of the accelerator grid apertures,and the experimental results also show that the higher background pressure does lead to faster erosion of the apertures.It is recommended that the background pressure can be taken as an influential factor for the lifetime acceleration test of ion thrusters.But at the same time,it should be noted that high background pressure will cause a higher probability of charge exchange collision and more xenon atoms to backflow into the discharge chamber,which will greatly increase the CEX ion density,and even the thruster is not working properly.The 30 cm diameter ion thruster was previously tested under background pressure of 7 × 10-3Pa,where beam interruption and grid breakdown occurred frequently.Therefore,in the design of the acceleration lifetime test,it is suggested that reasonable background pressure and test duration should be selected according to the acceleration stress level,which can fully verify the effectiveness of the influence factor,and the accelerator grid will not fail prematurely within the designed test duration.

    Acknowledgments

    This work is supported by Key Laboratory Funds for the Science and Technology on Vacuum Technology and Physics Laboratory,Lanzhou Institute of Physics (Nos.HTKJ2022KL510003 and 6142207210303),Independent project of Hangzhou Institute for Advanced Study (No.2022ZZ01009),and Science and Technology Project Affiliated to the Education Department of Chongqing Municipality(No.KJZD-K202101506).

    ORCID iDs

    猜你喜歡
    楊威
    引熱議!奧運(yùn)冠軍楊威讓女兒在家上學(xué)
    Precisely controlling the twist angle of epitaxial MoS2/graphene heterostructure by AFM tip manipulation
    采用樁基托換進(jìn)行既有鐵路橋加固分析
    河南科技(2022年9期)2022-05-31 00:42:40
    Thermally induced band hybridization in bilayer-bilayer MoS2/WS2 heterostructure?
    楊威:田園城市有多遠(yuǎn)
    商周刊(2018年16期)2018-08-14 01:51:42
    幻覺(jué)
    衛(wèi)國(guó)喋血女英烈
    鐵軍(2014年6期)2014-06-03 00:11:09
    “云”起“威”揚(yáng)冠軍夫妻婚后變身“全能選手”
    女性天地(2012年1期)2012-04-29 00:44:03
    誰(shuí)能拯救你
    嚴(yán)厲打造了他
    身体一侧抽搐| 国产91av在线免费观看| 淫秽高清视频在线观看| 国产 一区精品| 建设人人有责人人尽责人人享有的 | 国产真实乱freesex| 高清午夜精品一区二区三区| 久久国内精品自在自线图片| 亚洲无线观看免费| 亚洲内射少妇av| videos熟女内射| 久久精品人妻少妇| 亚洲一区高清亚洲精品| 中文欧美无线码| 欧美变态另类bdsm刘玥| 噜噜噜噜噜久久久久久91| 波多野结衣巨乳人妻| 久久久久久国产a免费观看| 又粗又爽又猛毛片免费看| 特大巨黑吊av在线直播| 免费电影在线观看免费观看| 小蜜桃在线观看免费完整版高清| 少妇猛男粗大的猛烈进出视频 | or卡值多少钱| 久久亚洲国产成人精品v| 一二三四中文在线观看免费高清| 亚洲av中文字字幕乱码综合| 18禁在线无遮挡免费观看视频| 中文欧美无线码| 91精品一卡2卡3卡4卡| 精品无人区乱码1区二区| 亚洲内射少妇av| 亚洲第一区二区三区不卡| 久久久久久久国产电影| 99久久无色码亚洲精品果冻| 久久久久久九九精品二区国产| 天堂av国产一区二区熟女人妻| 国产精品熟女久久久久浪| 亚洲av男天堂| 亚洲av中文av极速乱| 国产伦一二天堂av在线观看| 亚洲av电影不卡..在线观看| 美女xxoo啪啪120秒动态图| av免费观看日本| 中文欧美无线码| 国产伦精品一区二区三区四那| 嫩草影院入口| 欧美日韩精品成人综合77777| 天天一区二区日本电影三级| 成人鲁丝片一二三区免费| 国产视频首页在线观看| 欧美日韩一区二区视频在线观看视频在线 | 69人妻影院| 亚洲自拍偷在线| 久久精品熟女亚洲av麻豆精品 | 看免费成人av毛片| 超碰av人人做人人爽久久| 少妇的逼水好多| 99在线人妻在线中文字幕| 99热这里只有是精品50| 黄片wwwwww| 又爽又黄无遮挡网站| 免费看a级黄色片| 一区二区三区乱码不卡18| 日韩欧美在线乱码| 亚洲国产欧美在线一区| 又爽又黄无遮挡网站| 在线观看66精品国产| a级毛色黄片| 日本黄大片高清| 美女内射精品一级片tv| 亚洲成色77777| 中文字幕免费在线视频6| 亚洲最大成人av| 夜夜看夜夜爽夜夜摸| 99热这里只有精品一区| 亚洲国产精品成人综合色| 亚洲欧美日韩无卡精品| 波野结衣二区三区在线| 级片在线观看| 亚洲欧美日韩卡通动漫| 亚洲怡红院男人天堂| 国产熟女欧美一区二区| 久久久久久久久久久丰满| 性色avwww在线观看| 国产伦精品一区二区三区视频9| 国产精品一区二区三区四区免费观看| 色视频www国产| 久久久精品94久久精品| 国产乱人视频| 人人妻人人看人人澡| 麻豆成人午夜福利视频| 变态另类丝袜制服| 日日啪夜夜撸| 天堂av国产一区二区熟女人妻| 美女cb高潮喷水在线观看| 日本猛色少妇xxxxx猛交久久| 国产精品国产三级专区第一集| 少妇熟女aⅴ在线视频| 嘟嘟电影网在线观看| 亚洲av一区综合| 亚洲在久久综合| 纵有疾风起免费观看全集完整版 | 欧美一区二区亚洲| 秋霞伦理黄片| 国产成人aa在线观看| 噜噜噜噜噜久久久久久91| 亚洲成人精品中文字幕电影| 国产精品1区2区在线观看.| 国产又黄又爽又无遮挡在线| 国产精品精品国产色婷婷| 日本色播在线视频| 日韩,欧美,国产一区二区三区 | 久久久亚洲精品成人影院| 亚洲欧美成人精品一区二区| 欧美zozozo另类| 九九在线视频观看精品| 国产精品久久久久久精品电影小说 | 久久精品91蜜桃| 狂野欧美白嫩少妇大欣赏| 中文欧美无线码| 国产又黄又爽又无遮挡在线| 国产成人福利小说| 日韩成人伦理影院| 午夜亚洲福利在线播放| 91狼人影院| 欧美3d第一页| 亚洲最大成人手机在线| 亚洲美女视频黄频| 午夜福利网站1000一区二区三区| 亚洲国产日韩欧美精品在线观看| 亚洲在线观看片| 欧美丝袜亚洲另类| 亚洲欧美清纯卡通| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 91午夜精品亚洲一区二区三区| 国产伦一二天堂av在线观看| 国产高清三级在线| 欧美成人精品欧美一级黄| 97热精品久久久久久| 99久久精品国产国产毛片| 国产精品爽爽va在线观看网站| 欧美日本亚洲视频在线播放| 黄色欧美视频在线观看| 又爽又黄无遮挡网站| 亚洲综合精品二区| 超碰av人人做人人爽久久| 99热网站在线观看| 国产精品一区二区在线观看99 | 日韩大片免费观看网站 | 亚洲伊人久久精品综合 | 自拍偷自拍亚洲精品老妇| 白带黄色成豆腐渣| 亚洲人成网站在线观看播放| 免费黄频网站在线观看国产| 欧美精品av麻豆av| 国产精品国产三级专区第一集| 在线精品无人区一区二区三| 精品亚洲成a人片在线观看| 国产精品一二三区在线看| 久久国产亚洲av麻豆专区| 男女边吃奶边做爰视频| 日韩中字成人| 免费观看在线日韩| 97超碰精品成人国产| 激情视频va一区二区三区| 午夜福利影视在线免费观看| 黑人欧美特级aaaaaa片| 五月玫瑰六月丁香| 美女国产视频在线观看| 欧美日韩视频精品一区| av网站免费在线观看视频| 欧美人与善性xxx| 80岁老熟妇乱子伦牲交| 亚洲成av片中文字幕在线观看 | 王馨瑶露胸无遮挡在线观看| 十分钟在线观看高清视频www| 成年美女黄网站色视频大全免费| a级片在线免费高清观看视频| 国产免费一区二区三区四区乱码| 成人国产av品久久久| 欧美精品高潮呻吟av久久| a 毛片基地| 卡戴珊不雅视频在线播放| 亚洲综合色网址| 中文字幕人妻熟女乱码| www.av在线官网国产| 在线观看免费视频网站a站| 免费日韩欧美在线观看| 在线观看三级黄色| 这个男人来自地球电影免费观看 | 精品一区二区三区四区五区乱码 | h视频一区二区三区| 97人妻天天添夜夜摸| 最近中文字幕高清免费大全6| 久久久久久伊人网av| 黄色配什么色好看| 久久久a久久爽久久v久久| 久久精品国产a三级三级三级| 国产精品久久久久久久久免| 国产精品99久久99久久久不卡 | 涩涩av久久男人的天堂| 久久精品夜色国产| 欧美国产精品va在线观看不卡| 国产综合精华液| 男女边吃奶边做爰视频| 极品人妻少妇av视频| 日韩一本色道免费dvd| 国产免费现黄频在线看| 国产在线视频一区二区| 两个人免费观看高清视频| 在线亚洲精品国产二区图片欧美| 性色av一级| 51国产日韩欧美| 国产成人aa在线观看| 黄网站色视频无遮挡免费观看| 亚洲精品国产色婷婷电影| 多毛熟女@视频| 亚洲一级一片aⅴ在线观看| 女人精品久久久久毛片| av在线观看视频网站免费| 90打野战视频偷拍视频| 免费大片18禁| 少妇的丰满在线观看| 亚洲美女黄色视频免费看| 男人操女人黄网站| 免费观看a级毛片全部| 超色免费av| 久久99一区二区三区| 国产精品一区二区在线观看99| 成年人免费黄色播放视频| 赤兔流量卡办理| 日韩精品免费视频一区二区三区 | 久热久热在线精品观看| 91在线精品国自产拍蜜月| 黄色毛片三级朝国网站| 国产精品国产三级国产av玫瑰| 一区二区三区精品91| 精品久久久精品久久久| 人妻一区二区av| 久久久久久伊人网av| 51国产日韩欧美| 国产精品久久久av美女十八| 国产精品99久久99久久久不卡 | 夫妻性生交免费视频一级片| 午夜福利在线观看免费完整高清在| 成人毛片60女人毛片免费| av在线观看视频网站免费| 亚洲精品久久久久久婷婷小说| 亚洲精品aⅴ在线观看| 亚洲一码二码三码区别大吗| 精品福利永久在线观看| 日韩av在线免费看完整版不卡| 草草在线视频免费看| 看非洲黑人一级黄片| 99香蕉大伊视频| av福利片在线| 日韩三级伦理在线观看| 一级毛片 在线播放| 日本色播在线视频| 日产精品乱码卡一卡2卡三| 亚洲国产精品一区二区三区在线| 日本免费在线观看一区| 日本免费在线观看一区| 亚洲一码二码三码区别大吗| 亚洲国产精品一区二区三区在线| 男的添女的下面高潮视频| 成人毛片a级毛片在线播放| 亚洲精品乱码久久久久久按摩| 日韩精品有码人妻一区| 国产成人精品久久久久久| 国产精品久久久久久久电影| 成人亚洲欧美一区二区av| 亚洲中文av在线| 欧美最新免费一区二区三区| 国产高清不卡午夜福利| 在线天堂中文资源库| 黄色配什么色好看| 国产不卡av网站在线观看| 国产日韩欧美在线精品| 丝袜喷水一区| 热99久久久久精品小说推荐| 女的被弄到高潮叫床怎么办| 欧美成人精品欧美一级黄| 性色av一级| 久久久久久人人人人人| 亚洲国产精品国产精品| 国产免费又黄又爽又色| 成年女人在线观看亚洲视频| 在线观看一区二区三区激情| 成人综合一区亚洲| 久久精品熟女亚洲av麻豆精品| 在线观看免费视频网站a站| 欧美激情 高清一区二区三区| 在线精品无人区一区二区三| 极品少妇高潮喷水抽搐| 精品亚洲乱码少妇综合久久| 午夜免费男女啪啪视频观看| 少妇高潮的动态图| 天堂中文最新版在线下载| 母亲3免费完整高清在线观看 | 高清毛片免费看| 亚洲精品一区蜜桃| 青春草视频在线免费观看| 日韩一区二区视频免费看| 18禁国产床啪视频网站| 久久久久国产网址| 久久久久久久久久人人人人人人| 欧美xxⅹ黑人| 香蕉国产在线看| 大片电影免费在线观看免费| 91久久精品国产一区二区三区| 9色porny在线观看| 亚洲成人手机| 国产精品久久久久成人av| 99热这里只有是精品在线观看| 亚洲一码二码三码区别大吗| 一区二区三区精品91| 亚洲国产av影院在线观看| 高清av免费在线| 亚洲,一卡二卡三卡| 最后的刺客免费高清国语| 性色av一级| 美女中出高潮动态图| 欧美激情极品国产一区二区三区 | 黑人猛操日本美女一级片| 99久久综合免费| 老熟女久久久| 中国美白少妇内射xxxbb| 美女脱内裤让男人舔精品视频| 精品人妻偷拍中文字幕| 日韩av在线免费看完整版不卡| 国产极品粉嫩免费观看在线| 国产精品人妻久久久影院| 成人漫画全彩无遮挡| 亚洲国产最新在线播放| 久久人人爽人人片av| 十八禁高潮呻吟视频| 午夜免费男女啪啪视频观看| 美女脱内裤让男人舔精品视频| 视频在线观看一区二区三区| 久久99蜜桃精品久久| 又粗又硬又长又爽又黄的视频| 精品少妇久久久久久888优播| 中文字幕精品免费在线观看视频 | 波野结衣二区三区在线| 午夜福利网站1000一区二区三区| 哪个播放器可以免费观看大片| 中文精品一卡2卡3卡4更新| 日本av手机在线免费观看| 女性被躁到高潮视频| 久久精品国产亚洲av天美| 日韩成人伦理影院| 中文天堂在线官网| 国产男人的电影天堂91| 免费大片18禁| 国产在视频线精品| 亚洲av在线观看美女高潮| 国产成人精品福利久久| 99精国产麻豆久久婷婷| 婷婷色麻豆天堂久久| xxx大片免费视频| 欧美日韩国产mv在线观看视频| 国产日韩欧美亚洲二区| 人人妻人人澡人人爽人人夜夜| 交换朋友夫妻互换小说| 精品福利永久在线观看| 男人操女人黄网站| 免费久久久久久久精品成人欧美视频 | 少妇人妻 视频| 美女xxoo啪啪120秒动态图| 好男人视频免费观看在线| 妹子高潮喷水视频| 18+在线观看网站| 22中文网久久字幕| 中文精品一卡2卡3卡4更新| 久久 成人 亚洲| 激情视频va一区二区三区| 午夜免费男女啪啪视频观看| 国产高清不卡午夜福利| 少妇人妻 视频| 中文字幕av电影在线播放| 啦啦啦啦在线视频资源| 久久精品国产a三级三级三级| 成人综合一区亚洲| 日韩熟女老妇一区二区性免费视频| 观看av在线不卡| 91精品国产国语对白视频| 天天躁夜夜躁狠狠躁躁| 久久久国产欧美日韩av| 国产精品.久久久| 波多野结衣一区麻豆| 久久热在线av| 97精品久久久久久久久久精品| 成人午夜精彩视频在线观看| 亚洲一区二区三区欧美精品| 爱豆传媒免费全集在线观看| 嫩草影院入口| 老熟女久久久| 黄色一级大片看看| 91aial.com中文字幕在线观看| 精品一区在线观看国产| 最近手机中文字幕大全| 狠狠精品人妻久久久久久综合| 成人午夜精彩视频在线观看| 插逼视频在线观看| 欧美xxⅹ黑人| 婷婷色综合大香蕉| 黄片播放在线免费| www日本在线高清视频| 国产精品一区www在线观看| 老熟女久久久| 美女中出高潮动态图| 在线精品无人区一区二区三| 超色免费av| √禁漫天堂资源中文www| 午夜91福利影院| 精品99又大又爽又粗少妇毛片| 久久久国产精品麻豆| 精品一区二区免费观看| 亚洲色图综合在线观看| 在线 av 中文字幕| videos熟女内射| 国产有黄有色有爽视频| 色吧在线观看| 91在线精品国自产拍蜜月| 国产精品久久久av美女十八| 免费观看av网站的网址| 亚洲国产欧美在线一区| 久久韩国三级中文字幕| a级毛色黄片| 国产精品久久久av美女十八| 亚洲欧美一区二区三区国产| 久久久久久久大尺度免费视频| 十八禁网站网址无遮挡| 2022亚洲国产成人精品| 2018国产大陆天天弄谢| 咕卡用的链子| 乱码一卡2卡4卡精品| 色94色欧美一区二区| 免费在线观看黄色视频的| 制服人妻中文乱码| 少妇熟女欧美另类| 久热这里只有精品99| 国产精品一区二区在线观看99| 少妇人妻 视频| 久久国内精品自在自线图片| 精品人妻在线不人妻| 男人操女人黄网站| 精品国产国语对白av| 两个人免费观看高清视频| 波多野结衣一区麻豆| 人人妻人人添人人爽欧美一区卜| 黑人巨大精品欧美一区二区蜜桃 | 王馨瑶露胸无遮挡在线观看| 伊人久久国产一区二区| 婷婷色麻豆天堂久久| 亚洲欧美日韩卡通动漫| av在线app专区| 91午夜精品亚洲一区二区三区| 一二三四中文在线观看免费高清| 在线观看一区二区三区激情| 亚洲一区二区三区欧美精品| 一二三四中文在线观看免费高清| 精品亚洲成a人片在线观看| 男女啪啪激烈高潮av片| 久久久久久久精品精品| tube8黄色片| 国产精品一国产av| videossex国产| 80岁老熟妇乱子伦牲交| 爱豆传媒免费全集在线观看| 国产乱人偷精品视频| 亚洲欧美精品自产自拍| av免费在线看不卡| 亚洲国产色片| 亚洲成av片中文字幕在线观看 | 久久久久国产精品人妻一区二区| kizo精华| 22中文网久久字幕| 观看美女的网站| 91成人精品电影| 午夜日本视频在线| 大片免费播放器 马上看| 人妻 亚洲 视频| 亚洲美女视频黄频| a级毛片黄视频| 亚洲精品美女久久av网站| 精品久久蜜臀av无| 另类精品久久| 国产男女超爽视频在线观看| 五月伊人婷婷丁香| 欧美少妇被猛烈插入视频| 国语对白做爰xxxⅹ性视频网站| 自线自在国产av| 赤兔流量卡办理| 午夜福利网站1000一区二区三区| 国产麻豆69| 亚洲综合精品二区| 亚洲,欧美,日韩| 欧美bdsm另类| 亚洲欧美一区二区三区国产| 男的添女的下面高潮视频| 精品一区二区三区视频在线| 欧美日韩综合久久久久久| 亚洲 欧美一区二区三区| 丝袜喷水一区| 免费人妻精品一区二区三区视频| 卡戴珊不雅视频在线播放| 国产xxxxx性猛交| 免费在线观看黄色视频的| 国产 精品1| 99国产精品免费福利视频| 777米奇影视久久| 在线亚洲精品国产二区图片欧美| 人成视频在线观看免费观看| 国产成人a∨麻豆精品| 80岁老熟妇乱子伦牲交| 亚洲国产av新网站| 欧美成人午夜精品| 天堂中文最新版在线下载| 91在线精品国自产拍蜜月| 丝袜喷水一区| 午夜av观看不卡| 欧美人与性动交α欧美软件 | 草草在线视频免费看| 久久精品夜色国产| 五月开心婷婷网| 男女边吃奶边做爰视频| av线在线观看网站| 国产精品一国产av| www日本在线高清视频| 欧美激情极品国产一区二区三区 | 三上悠亚av全集在线观看| 成人毛片a级毛片在线播放| 精品国产国语对白av| 国产精品不卡视频一区二区| av又黄又爽大尺度在线免费看| 欧美成人精品欧美一级黄| 亚洲一区二区三区欧美精品| av国产久精品久网站免费入址| 人人妻人人爽人人添夜夜欢视频| 91久久精品国产一区二区三区| 久久国产精品男人的天堂亚洲 | 久久久久久人人人人人| 亚洲av国产av综合av卡| 国内精品宾馆在线| 亚洲一码二码三码区别大吗| 丝袜脚勾引网站| 18在线观看网站| 性色av一级| 如日韩欧美国产精品一区二区三区| 韩国av在线不卡| 欧美精品一区二区大全| 美女福利国产在线| 大香蕉久久网| 黑人欧美特级aaaaaa片| 国产成人午夜福利电影在线观看| 午夜激情久久久久久久| 久久久久久久亚洲中文字幕| av片东京热男人的天堂| 99久久人妻综合| 免费不卡的大黄色大毛片视频在线观看| 丝袜喷水一区| 97在线视频观看| 80岁老熟妇乱子伦牲交| 国产精品人妻久久久影院| 97在线人人人人妻| 少妇人妻 视频| 国产亚洲最大av| 欧美变态另类bdsm刘玥| 又粗又硬又长又爽又黄的视频| 涩涩av久久男人的天堂| 搡女人真爽免费视频火全软件| 欧美日韩av久久| 黄色一级大片看看| 大片电影免费在线观看免费| 香蕉精品网在线| 永久网站在线| 久久人人爽av亚洲精品天堂| 男女午夜视频在线观看 | 丰满乱子伦码专区| 国产深夜福利视频在线观看| 视频在线观看一区二区三区| 久久午夜福利片| 午夜视频国产福利| 菩萨蛮人人尽说江南好唐韦庄| 欧美精品亚洲一区二区| 免费av不卡在线播放| 黑人猛操日本美女一级片| 亚洲美女搞黄在线观看| 日韩av在线免费看完整版不卡| 国产成人a∨麻豆精品| 午夜视频国产福利| 亚洲美女视频黄频| 99热6这里只有精品| 日本-黄色视频高清免费观看| 综合色丁香网| 美女中出高潮动态图| 国产成人91sexporn| 国产无遮挡羞羞视频在线观看| 大话2 男鬼变身卡| 久久国产亚洲av麻豆专区| 免费看光身美女| 啦啦啦视频在线资源免费观看| 成人亚洲精品一区在线观看| 99九九在线精品视频| 寂寞人妻少妇视频99o| 国产精品久久久久久精品电影小说| a级毛片在线看网站| 三级国产精品片| 亚洲天堂av无毛| 久久人人97超碰香蕉20202| 国产成人91sexporn| 亚洲国产毛片av蜜桃av|