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

    Numerical study of the grid erosion of field emission electric propulsion

    2021-10-31 08:14:32ChangLU鹿暢YangLUO羅楊GuangqingXIA夏廣慶HuiGAO高輝andNuoXU許諾
    Plasma Science and Technology 2021年10期
    關(guān)鍵詞:許諾

    Chang LU (鹿暢), Yang LUO (羅楊), Guangqing XIA (夏廣慶),Hui GAO (高輝) and Nuo XU (許諾)

    1 State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology,Dalian 116024, People’s Republic of China

    2 Key Laboratory of Advanced Technology for Aerospace Vehicles of Liaoning Province, Dalian University of Technology, Dalian 116024, People’s Republic of China

    3 Beijing Institute of Mechanical Equipment, Beijing 100854, People’s Republic of China

    Abstract In this paper,the immersed finite element particle-in-cell Monte Carlo collision(IFE-PIC-MCC)model is used to study the cause of the grid erosion in field emission electric propulsion(FEEP).The simulation results show that charge exchange (CEX) ions are the main cause of the grid erosion, while beam ions basically do not impinge on the grid.The CEX ions are mainly generated near the grid.Some of these CEX ions return to the upper surface, lower surface and notch side surface of the grid.The impact angle of CEX ions returning to the upper and side surfaces is large, but their energy is low, while the impact angle of CEX ions returning to the lower surface is small, but their energy is high.The above characteristics lead to the similar erosion rates of these three surfaces.

    Keywords: IFE-PIC-MCC, FEEP, grid erosion, CEX ions

    1.Introduction

    With the development of manufacturing technology to micro/nanoscale, the size of satellites also becomes smaller and smaller.As a result, micro or nano satellites have decome more predominant[1].With the advantages of small size,low power consumption, short development cycle, formation networking and low cost, micro or nano satellites play an important role in scientific research, national defense, commercial sector, etc [2–4].However, their limited size and power supply demand low power, small volume and small mass.In addition, other space missions (gravitational wave detection, Aurora communication, etc) require higher thrust accuracy.For example, the Laser Interferometer Space Antenna (LISA) plan needs 5–30 μN thrust, thrust accuracy≤0.1 μN, and thrust noise≤0.1 μN/√Hz [5].Conventional chemical and cold gas thrusters struggle to meet the thrust requirements of the LISA plan.Field emission electric propulsion(FEEP)using liquid metal as a working fluid is one of the excellent solutions to realize such space missions [6].Its structure is very compact, and is expected to achieve high specific impulse, high efficiency, and high thrust accuracy.

    The thrust generation of FEEP is completely different from that of traditional electric propulsion, which uses gas ionization and accelerated plasma jets.FEEP uses liquid metal as a working medium to avoid the use of high-pressure tanks and working medium transport [7].Meanwhile, the energy needed to produce ions is much lower than the ionization energy of most inert gases, which further avoids the loss of working medium and efficiency caused by the low ionization rate.Moreover, FEEP can also avoid the performance degradation caused by the size effect, which commonly occurs in micro propulsion.

    However, the above characteristics of FEEP also bring some problems.The ionization in FEEP is based on the field emission, which requires a very high voltage to be applied between the grid and the emitter.Due to the high conductivity of metals (106–107S m?1), the very high surface tension of liquid metal makes the working voltage of FEEP generally above 5 kV[8–10].This makes the specific impulse of FEEP very high,generally between 4000 s and 8000 s.On the other hand, it also results in the high kinetic energy of the ions impacting the grid, which then causes serious grid erosion.Furthermore, in the experiment, the serious grid erosion is accompanied by a large amount of refluxed charge exchange(CEX) ions, which further aggravates the grid erosion.However, most research has mainly focused on the optimization of FEEP grid parameters.For example, the Centrospazio space research laboratory has done a lot of research on narrow slit FEEP, including micro thrust, plume characteristic,life,etc[11–13].Genoveseet alstudied the launch performance of narrow slot FEEP and found that there is an exponential relationship between the launch current and the launch voltage, and the thrust has a linear relationship with the power [11].Marcuccioet alcalculated the efficiency,thrust and specific impulse of FEEP, and proposed that the performance of FEEP mainly depends on the emission electrical parameters and the geometric parameters of FEEP slit[12];Andrenucciet alcarried out the numerical simulation of FEEP, mainly used to find the best geometric parameters of FEEP,so as to improve the performance of FEEP[13].It can be seen that the mechanism of FEEP grid erosion is still not well studied.

    Therefore, the grid erosion characteristics of FEEP are studied in this paper.The mechanism of grid erosion is analyzed by the three-dimensional immersed finite element particle-in-cell Monte Carlo collision (IFE-PIC-MCC) model[14–17].In section 2, the simulation model and parameter settings based on the IFE-PIC-MCC method are given.Section 3 presents the simulation results and discussion.Finally, the conclusion is given in section 4.

    2.Simulation model

    The working principle of FEEP is as follows: in the working state,the liquid metal propellant after heating is continuously transported to the emitter tip through the micro channel on the emitter surface under capillary action.Meanwhile, a highvoltage electric field of 6–10 kV is applied between the absorber and emitter.According to the principle of field emission and field evaporation,metal ions form a Taylor cone under the action of surface tension and electric field force,and finally form the emission plume.

    In view of the above process, we propose the following assumptions.First, ignoring the formation process of the Taylor cone, ions are emitted directly from the tip.Because the Taylor cone is in the nanometer scale, which is much smaller than the grid size, the initial position of all ions is located at the same point of the emitter tip.The initial velocity of ions agrees with Maxwell distribution.Second, since most of the electrons are absorbed by the emitter in the field emission process, it is assumed that there are no electrons between the emitter and the grid.However, in the quasi neutral region downstream of the grid, electrons are assumed to conform to the Boltzmann distribution.Third, because the flow rate of neutral atoms is relatively small,usually 1/100 of the ion flow rate, and the macro velocity of neutral atoms(about 30–40 m s?1)is far less than its thermal velocity(about 300 m s?1), it can be considered that the macro velocity of neutral atoms is close to 0, and the distribution of neutral atoms is relatively uniform.In addition, only CEX collisions between ions and atoms are considered, and other collisions are ignored.

    Based on the above assumptions, the calculation process of our model is as follows.First, the motion of ions is calculated to obtain the beam ion distribution.Then, with the uniformly distributed atoms, the CEX collisions are calculated to obtain the CEX ion distribution.Finally, the ions impinging on the grid surface are collected to analyze the erosion mechanism.Next,the specific calculation methods of ion motion, CEX collision and grid erosion are introduced.

    2.1.Ion motion and field solution

    The movement of charged particles in the simulation region will generate an electromagnetic field.The electromagnetic field formed by ion motion is given by Maxwell equations.Since the magnetic field produced by the ions has little influence on their movements, the electromagnetic field can be simplified as a problem of the electrostatic field.Hence,the Maxwell equations are simplified as follows:

    where Φ is the potential, E is the electric field,niandneare the ion charge density and electron charge density, and ε0is the dielectric constant of the vacuum.Note that in the region between the emitter and the grid,ne= 0; while the electron density downstream of grid is obtained by

    wheren∞is the ion density downstream of ion optics,eis the elementary charge quantity,kis the Boltzmann constant, Φ∞is the plasma potential downstream of ion optics, andTe∞is the temperature of the electron downstream of ion optics.

    Then, the force of the ion withqin the electromagnetic field can be obtained from equation (3),

    where F is the electric field force,andqis the charge quantity of ions.

    In our model, the immersed finite element (IFE) method[14–17] is used to solve the electric field.The IFE method developed in recent years is a promising numerical tool for solving isotropic/anisotropic Poisson equations in structured meshes.The IFE method uses the slice basis function technique in the element divided by the interface,so that the mesh division does not depend on the interface.The other elements,which are not divided by the interface, use the traditional finite element basis functions.Therefore, the IFE method can effectively solve electromagnetic fields with complex boundary conditions in structured meshes.Moreover, compared with the traditional finite element method, the IFE method has few changes, so it is easy to program.

    After obtaining the electric field, the motion of ions can be solved by Newton’s second law,

    wheremiis the ion mass,viis the velocity of ions,and x is the position of ions.Then, equation (4) is solved by Boris’s leapfrog method [18].

    2.2.Charge exchange collision

    The low-speed CEX ions accelerate under the high voltage of the grid system, and some of them will bombard the surface of the grid, causing erosion.Assuming that the neutral atoms are uniformly distributed in the calculation area and obey the law of ideal gas, the number density of neutral atoms can be obtained from equation (5) [19],

    wherennis the atom density,Iiis the beam currents, η is the ratio of atom to ion flow rate,vnis the atom velocity,andAsis the outlet area of atoms.In the region between the emitter and the grid,Asis the area of the grid slot.In the region downstream of the grid,Asis the exit of the calculation region.

    Since the velocity of ions is much higher than that of atoms,in the calculation,the relative velocity between ions and atoms is replaced by the velocity of ions.The cross-sectionsσT(vi) of CEX collisions between atoms and ions is obtained as follows.

    wherek1andk2are the collision section coefficients,k1=?1.492×10?10s andk2=2.6997×10?9m [19],viis the mode of ion velocity,andσT(vi)is the cross-section of CEX collisions.The collision frequencyμiis calculated by eq u ation (7),

    wherenn(xi) is the atom density atxi.Then the collision probabilityPican be calculated according to equation (8),

    whereΔtis the time step, andPiis the collision probability.

    The null collision technique is used to deal with the particle collision in the Monte Carlo method [19].First, the maximum collision frequencyμmof all ions is calculated,

    whereμmis the maximum collision frequency of all ions.

    Then the maximum value of the collision probabilityPnullcan be obtained by using the maximum value of the collision frequency,

    wherePnullis the maximum value of the collision probability.

    The total number of ions in the calculation domain is expressed byNtot,and the collision frequencies μ of all ions in one step can be calculated by randomly selectingNtotPnullions.Selecting a random numberRbetween 0 and 1, ifμ μm>R,the collision occurs.

    2.3.Grid erosion

    The number of atoms bombarded by an ion from the material surface is called the sputtering yield,which is expressed byY.The sputtering yield is related to the energy of ion bombardment on the material, the angle between the impact trajectory and the normal of the material surface, and the physical properties of the material atoms.The erosion depth per unit time of a point on the grid surface impacted by incident ions is called the erosion rate,which is expressed byRE.The erosion depth of the grid can be obtained by multiplying the erosion rate by the working time.On a small areaS, it is assumed that ions uniformly bombard the surface of the grid.Assuming thatNions with unit charge bombard the small areaSin unit time, the current density formed on the grid surface can be calculated by equation (11),

    whereJis the current density,Sis a small area, andNis the number of ions impactingS.According to the definition of sputtering yield, the number of atoms bombarded by the incident ions in the areaSper unit time can be obtained from equation (12),

    whereNais the number of the bombarded atoms,andYis the sputtering yield.

    If the mass of the material bombarded from the grid is expressed asM,the volumeVof the material can be obtained from equation (13),

    whereVis the volume of the whole bombarded atoms,ρgis the density of the gird material, andMgis the atom mass of the gird material.In this way, we can get the following expression according to the definition of erosion rate,

    whereREis the erosion rate.

    2.4.Sputtering yield

    The grid material is molybdenum.According to the experimental data of molybdenum, the formula of sputtering yieldYis [20],

    whereEXeis the ion impact energy, and

    The sputtering yield is not only affected by the incident energy, but also related to the incident angleθi,hence

    whenθi≤0.689,equation (18) is used to correctY,

    whenθi> 0.689,equation (19) is used to correctY,

    2.5.Parameter setting

    The working medium is indium.The relative atomic mass of indium is 114.5, which is about 1.91×10?25kg.The radius of the emitter tip is about 10 μm and the length of the emitter extending into the simulation domain is 0.05 mm.The emission current of the emitter is 100 μA.The ratio η is 2:3.In the region between the emitter and the grid,Asis about 1.04 mm2,while in the region downstream of the grid,Asis about 3 mm2.The velocity of the atom is set to 1000 ms–1.In our simulation,the neutral atoms are assumed to be uniformly distributed[19, 21]; hence the atom density is set to be 1.6×1018m?3between the emitter and the grid and 5.5×1017m?3downstream of the grid,according to the atom flow rate.The emitter voltage is set to 6 kV,and the grid voltage is set to ?8 kV.The input parameters of the grid geometry are shown in table 1.

    Table 1.Parameters of the FEEP grid model (normalized).

    As shown in figure 1, the simulation model is simplified as a 60×20×50 cuboid with a mesh size of 0.05 mm.As for the field boundary condition,the surface withZ=Zmaxis set as the Dirichlet boundary,and the potential value is equal to 0, while the other surfaces are set as zero Newman boundary condition.As for the particle boundary conditions,theZ=ZminandZ=Zmaxsurfaces are set as absorption boundaries,while the others are set as symmetric boundaries.In addition, when the ion hits the grid, it is considered to be absorbed by the grid, so the ion is deleted in the simulation.

    Finally, the simulation results in this paper are all normalized,and the reference parameters can be found in reference[22].

    3.Simulation results and analysis

    Figures 2–4 show the potential, ion and CEX ion density distributions, respectively.Because the distance between the grid and the emitter is very close,a large potential gradient is formed,while the potential gradient downstream of the grid is relatively small, and the equipotential surface is basically parallel to the grid.Hence, between the grid and the emitter,the beam ions will be accelerated by the strong electric field.While,in the downstream of the grid,the beam ion is mainly affected by the vertical downward electric field, which will cause the gradually increase of the beam divergence angle.

    Figure 1.The geometry model of the FEEP grid.

    Figure 2.Potential distribution.

    Figure 3.Beam ion density distribution.

    Figure 4.Density distribution of beam ions and CEX ions.

    Figure 5.Comparison of density distribution between beam ion and CEX ion.

    Figure 6.Velocity phase distribution of beam ions and CEX ions.

    In figure 4, the green to red parts are mainly beam ions,and the blue parts are mainly CEX ions.From the density distribution of beam ions and CEX ions, it can be seen that the beam ions can be extracted and accelerated well, and no beam ions impinge directly on the grid.However, the distribution of CEX ions is relatively wide, with a large number of CEX ions impacting the grid slot and the end faces.

    The density distributions of beam ions and CEX ions are compared in figure 5.In order to show the difference between the distribution of CEX ions and beam ions,the upper part of figure 5 contains both CEX ions and beam ions, while the lower part contains only beam ions.CEX ions are mainly concentrated near the grid.This is mainly due to the high velocity and high density of ions near the grid,which leads to obvious CEX collisions.However, in the downstream of the grid,the ion velocity and density decrease significantly so that the probability of CEX collision decreases.

    Figure 6 shows the velocity phase distribution comparison of beam ions and CEX ions.Consistent with the analysis in figure 1, the beam ions undergo a rapid acceleration process and then slow deceleration, but the velocity of the beam ions is always along the positive direction of thez-axis.However,the velocity of most CEX ions is along the negative direction of thez-axis, which indicates that most CEX ions will return to the grid and cause grid erosion.

    Figure 7.Distribution of CEX ions, which cause the grid erosion.

    Figure 8.Location distribution of CEX ions causing grid erosion.

    Figures 7 and 8 show the distribution of CEX ions,which cause the grid erosion, and their generation positions.In figures 7 and 8, thez-coordinate means the generation positions of CEX ions on thez-axis.It can be seen that these CEX ions mainly come from the position betweenz=20–70.This is consistent with the CEX ion distribution shown in figures 5 and 6.That is, the CEX ions causing the grid erosion are mainly generated between the grid and the emitter, and near the downstream of the grid.

    Figures 9–12 show the CEX ion angle, energy, current density and erosion rate distribution of the upper surface,respectively.It can be seen that their impact angles on the upper surface is large, which is basically between 80° and 90°.The impact energy is about 1000 eV.The current density distribution shows that CEX ions are mainly concentrated near the grid notch.Due to the relatively small impact angle near the notch, the final erosion rate distribution on the upper face is relatively uniform, which is concentrated at 3.9×10?4m/1000 h.

    Figure 9.Impact angular distribution of the CEX ions impacting the upper surface.

    Figure 10.Energy distribution of the CEX ions impacting the upper surface.

    Figure 11.Current density distribution of the CEX ions impacting the upper surface.

    Figure 12.Erosion rate distribution of the upper surface.

    Figure 13.Impact angular distribution of the CEX ions impacting the lower surface.

    Figure 14.Energy distribution of the CEX ions impacting the lower surface.

    Figure 15.Current density distribution of the CEX ions impacting the lower surface.

    Figures 13–16 show the CEX ion angle, energy, current density and erosion rate distribution of the lower surface.Comparing with figures 8–11,it can be seen that the CEX ion angle of the lower surface is very small, which is basically between 10° and 50°.This is because these CEX ions will first backflow from the downstream of the grid into the grid slot, then be rebounded by the emitter’s extremely high potential, and finally impinge on the lower surface.Meanwhile,because of the obvious radial electric field between the grid and the emitter, the impact angle of these CEX ions on the lower surface is very small.However,their impact energy is relatively high, about 2000 eV, and the current density distribution shows that CEX ions are also concentrated near the grid notch.Due to the small impact angle,the erosion rate distribution of the lower surface is similar to that of the upper surface, which is about 3.9×10?4m/1000 h.

    Figures 17–20 show the CEX ion angle, energy, current density and erosion rate distribution on the side surface of the notch.It can be seen that the CEX ion angle on the side surface is also large,which is basically between 70°and 90°.However, the impact energy of these CEX ions is relatively low, about 600 eV.The current density distribution shows that the CEX ions on the side surface are concentrated in the middle of the notch.Because of the large impact angle, the erosion rate distribution on the side surface of notch is also about 3.9×10?4m/1000 h.

    Figure 16.Erosion rate distribution of the lower surface.

    Figure 17.Impact angular distribution of the CEX ions impacting the side surface.

    Figure 18.Energy distribution of the CEX ions impacting the side surface.

    Figure 19.Current density distribution of the CEX ions impacting the side surface.

    Figure 20.Erosion rate distribution of the side surface.

    In summary, the grid erosion is mainly caused by the returned CEX ions.These CEX ions are mainly generated betweenz=20–70, that is, near the grid.These CEX ions return to the upper surface, lower surface and notch side surface of the grid under the influence of the grid negative potential.The impact angle of these CEX ions returning to the upper and side surfaces is large (about 70°–90°), but their energy is low(about 600–1000 eV),while the impact angle of CEX ions returning to the lower surface is very small (about 10°–50°),but their energy is high(about 2000 eV).The above characteristics lead to the similar erosion rates of these three surfaces (about 3.9×10?4m/1000 h).

    4.Conclusion

    In this paper, the three-dimensional simulation model for the FEEP grid is established by the IFE-PIC-MCC method.The erosion mechanism of the FEEP grid was analyzed by this model.The simulation results show that the grid erosion is mainly caused by the returned CEX ions,while the beam ions make little contribution to the grid erosion.The velocity and density of the beam ions near the grid are high, so the CEX collision frequency near the grid is high, leading to CEX ions mainly generated near the grid.Some of these CEX ions return to the upper surface,lower surface and notch side surface of the grid.The impact angle of CEX ions returning to the upper and side surfaces is large,but their energy is low,while the impact angle of CEX ions returning to the lower surface is small,but their energy is high.The above characteristics lead to the similar erosion rates of these three surfaces.

    Acknowledgments

    This work was supported by National Natural Science Foundation of China (No.11675040), the Fundamental Research Funds for the Central Universities of China (Nos.DUT20LAB203 and DUT21GJ206), the Key Research and Development Project of Liaoning Province (No.2020JH2/10500003), and the Open Fund for Tianqin Project Key Laboratory of Ministry of Education(in preparation)(TQRF-2020-2.5).

    猜你喜歡
    許諾
    許諾作品
    古道上的“雷鋒”
    開心果——朱小柯
    誤闖拉瑪大沙漠(一)
    情 話
    許諾
    參花(下)(2020年4期)2020-04-16 12:49:04
    現(xiàn)實與未來需要雙重許諾
    商周刊(2018年10期)2018-06-06 03:04:09
    許諾(外三首)
    天津詩人(2017年2期)2017-11-29 01:24:15
    網(wǎng)絡(luò)許諾銷售行為地域范圍的認定
    淺析英美法系合同中許諾拘束力的基礎(chǔ)
    考比视频在线观看| 国产精品.久久久| 亚洲 欧美一区二区三区| 亚洲五月婷婷丁香| 午夜影院在线不卡| 老司机在亚洲福利影院| 免费久久久久久久精品成人欧美视频| 免费看十八禁软件| 欧美变态另类bdsm刘玥| 久久久久久免费高清国产稀缺| 久久青草综合色| 亚洲欧美一区二区三区黑人| 国产日韩欧美在线精品| 色播在线永久视频| 亚洲国产中文字幕在线视频| 日本精品一区二区三区蜜桃| 午夜影院在线不卡| 亚洲国产成人一精品久久久| 精品福利观看| 九色亚洲精品在线播放| 大片免费播放器 马上看| 免费不卡黄色视频| 国产成人精品无人区| 一个人免费在线观看的高清视频 | 最黄视频免费看| 国内毛片毛片毛片毛片毛片| 国产精品.久久久| 一二三四社区在线视频社区8| 在线观看人妻少妇| 天天躁狠狠躁夜夜躁狠狠躁| 欧美+亚洲+日韩+国产| 午夜视频精品福利| 黄色 视频免费看| 国产一级毛片在线| 久久国产精品影院| 午夜免费鲁丝| 纵有疾风起免费观看全集完整版| 日韩熟女老妇一区二区性免费视频| 国产91精品成人一区二区三区 | 在线观看免费高清a一片| 午夜免费鲁丝| 中文字幕av电影在线播放| 午夜免费鲁丝| 99久久99久久久精品蜜桃| 午夜福利视频精品| 国产成人av激情在线播放| 亚洲精品国产色婷婷电影| 国产亚洲午夜精品一区二区久久| 欧美中文综合在线视频| 成人影院久久| 黑人欧美特级aaaaaa片| 伦理电影免费视频| 亚洲 欧美一区二区三区| 一区二区三区乱码不卡18| 久久精品亚洲熟妇少妇任你| 亚洲精品成人av观看孕妇| 999久久久精品免费观看国产| 女人被躁到高潮嗷嗷叫费观| 好男人电影高清在线观看| 国产精品麻豆人妻色哟哟久久| 少妇精品久久久久久久| 波多野结衣av一区二区av| 欧美+亚洲+日韩+国产| 久久国产精品影院| 丝袜在线中文字幕| 91成人精品电影| 亚洲国产毛片av蜜桃av| 丁香六月天网| 啦啦啦免费观看视频1| 男人添女人高潮全过程视频| 黑人巨大精品欧美一区二区蜜桃| av天堂在线播放| 日韩视频在线欧美| 高清视频免费观看一区二区| 美女高潮到喷水免费观看| 成人亚洲精品一区在线观看| 宅男免费午夜| 日韩视频一区二区在线观看| 精品免费久久久久久久清纯 | 国产成人影院久久av| 久久精品亚洲av国产电影网| 久久香蕉激情| 最新的欧美精品一区二区| 大陆偷拍与自拍| 亚洲一码二码三码区别大吗| 一级片免费观看大全| 色精品久久人妻99蜜桃| 美女大奶头黄色视频| 黑丝袜美女国产一区| 国产欧美日韩精品亚洲av| 国产精品久久久人人做人人爽| 女人被躁到高潮嗷嗷叫费观| 嫩草影视91久久| 亚洲全国av大片| 国产在线一区二区三区精| 乱人伦中国视频| 精品亚洲成a人片在线观看| 精品人妻1区二区| 免费黄频网站在线观看国产| 亚洲精品在线美女| 免费在线观看视频国产中文字幕亚洲 | videos熟女内射| 9色porny在线观看| 性高湖久久久久久久久免费观看| 一区二区日韩欧美中文字幕| 狠狠婷婷综合久久久久久88av| 国产色视频综合| av国产精品久久久久影院| 淫妇啪啪啪对白视频 | 少妇猛男粗大的猛烈进出视频| 国产男女内射视频| 黄色毛片三级朝国网站| 婷婷成人精品国产| 在线精品无人区一区二区三| 18在线观看网站| 一级片免费观看大全| 成人国产av品久久久| 亚洲国产欧美在线一区| 美国免费a级毛片| 嫩草影视91久久| 精品人妻熟女毛片av久久网站| 国产福利在线免费观看视频| 国产一区二区在线观看av| 免费在线观看日本一区| 国产精品久久久久久精品古装| 欧美日韩亚洲高清精品| 国产精品国产av在线观看| 亚洲av日韩精品久久久久久密| 日本91视频免费播放| 国产精品久久久人人做人人爽| 亚洲七黄色美女视频| 美女主播在线视频| 狂野欧美激情性xxxx| 色播在线永久视频| 老司机靠b影院| 青春草亚洲视频在线观看| 99国产精品一区二区三区| 亚洲 国产 在线| 国产成人系列免费观看| 中文字幕另类日韩欧美亚洲嫩草| 99热网站在线观看| 中文精品一卡2卡3卡4更新| 丰满少妇做爰视频| 丁香六月天网| 国产日韩欧美在线精品| 深夜精品福利| 男女边摸边吃奶| 亚洲欧美成人综合另类久久久| 欧美一级毛片孕妇| 丁香六月欧美| 亚洲欧美激情在线| 亚洲精品一二三| 日韩制服丝袜自拍偷拍| 18禁观看日本| 欧美黑人欧美精品刺激| 美女高潮喷水抽搐中文字幕| 亚洲人成电影观看| 最近最新免费中文字幕在线| 黑人猛操日本美女一级片| 高潮久久久久久久久久久不卡| 精品高清国产在线一区| 亚洲全国av大片| 精品少妇黑人巨大在线播放| 交换朋友夫妻互换小说| 国产男人的电影天堂91| 一边摸一边抽搐一进一出视频| tube8黄色片| 91av网站免费观看| 免费久久久久久久精品成人欧美视频| 少妇 在线观看| av欧美777| 男女免费视频国产| 各种免费的搞黄视频| 一本大道久久a久久精品| 午夜精品国产一区二区电影| 成在线人永久免费视频| 操美女的视频在线观看| 欧美性长视频在线观看| 搡老岳熟女国产| 男女无遮挡免费网站观看| 中文字幕人妻丝袜制服| 国产区一区二久久| 色综合欧美亚洲国产小说| 在线观看免费高清a一片| 欧美日韩亚洲高清精品| 午夜福利乱码中文字幕| 亚洲五月婷婷丁香| 免费一级毛片在线播放高清视频 | 亚洲欧美清纯卡通| 亚洲精品一区蜜桃| 久久久久精品人妻al黑| 久久久国产成人免费| 成人18禁高潮啪啪吃奶动态图| 午夜两性在线视频| av有码第一页| 人妻久久中文字幕网| 性色av一级| 欧美日韩中文字幕国产精品一区二区三区 | av超薄肉色丝袜交足视频| 国产熟女午夜一区二区三区| 热re99久久国产66热| 99国产极品粉嫩在线观看| 老汉色av国产亚洲站长工具| 久久久精品国产亚洲av高清涩受| 热re99久久国产66热| 日韩中文字幕欧美一区二区| 五月开心婷婷网| 成年人黄色毛片网站| 国产亚洲欧美精品永久| 国产精品国产三级国产专区5o| 国产亚洲av高清不卡| 黄片小视频在线播放| 国产精品成人在线| 亚洲国产欧美网| 亚洲伊人色综图| 国内毛片毛片毛片毛片毛片| 亚洲国产中文字幕在线视频| 久热爱精品视频在线9| 精品人妻在线不人妻| 国产成人精品久久二区二区91| 亚洲欧洲日产国产| 老汉色∧v一级毛片| 国产一卡二卡三卡精品| 99热国产这里只有精品6| 淫妇啪啪啪对白视频 | 精品国产一区二区三区四区第35| 中国美女看黄片| 国产日韩一区二区三区精品不卡| 久久久久久人人人人人| 久久精品国产亚洲av高清一级| 大片电影免费在线观看免费| 国产日韩欧美视频二区| av天堂久久9| 美女视频免费永久观看网站| 亚洲精品国产精品久久久不卡| 日韩制服骚丝袜av| 老司机影院成人| avwww免费| 视频在线观看一区二区三区| 亚洲精品久久午夜乱码| 国产av一区二区精品久久| av线在线观看网站| tube8黄色片| 亚洲黑人精品在线| 天堂8中文在线网| 99久久人妻综合| 51午夜福利影视在线观看| 国产亚洲午夜精品一区二区久久| 丰满饥渴人妻一区二区三| 日韩,欧美,国产一区二区三区| 国产成人精品无人区| 国产成人免费观看mmmm| 国产精品国产三级国产专区5o| 日韩三级视频一区二区三区| 亚洲av男天堂| 久久人人爽人人片av| 一本一本久久a久久精品综合妖精| 精品熟女少妇八av免费久了| 亚洲va日本ⅴa欧美va伊人久久 | 90打野战视频偷拍视频| www.精华液| 另类精品久久| 亚洲精品国产精品久久久不卡| 操美女的视频在线观看| 久久久久精品国产欧美久久久 | 亚洲国产精品999| 如日韩欧美国产精品一区二区三区| 久久 成人 亚洲| www.自偷自拍.com| 国产伦人伦偷精品视频| 亚洲情色 制服丝袜| www.自偷自拍.com| 日韩人妻精品一区2区三区| 黑人欧美特级aaaaaa片| 精品亚洲成国产av| 中亚洲国语对白在线视频| 看免费av毛片| 精品人妻1区二区| 成年美女黄网站色视频大全免费| 9色porny在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 在线观看免费视频网站a站| a在线观看视频网站| 国产av国产精品国产| 极品少妇高潮喷水抽搐| 又大又爽又粗| 一本综合久久免费| 国产精品国产av在线观看| 波多野结衣一区麻豆| av在线播放精品| 好男人电影高清在线观看| 免费观看人在逋| 欧美亚洲 丝袜 人妻 在线| 啦啦啦视频在线资源免费观看| 丁香六月天网| 亚洲av电影在线进入| 777米奇影视久久| 亚洲中文字幕日韩| 99国产精品免费福利视频| 精品一区二区三区四区五区乱码| 国产精品99久久99久久久不卡| 亚洲av成人不卡在线观看播放网 | 成人三级做爰电影| 久久久久久久大尺度免费视频| 国产有黄有色有爽视频| 亚洲精品国产色婷婷电影| 久久久久久人人人人人| 亚洲国产欧美一区二区综合| 脱女人内裤的视频| av在线app专区| a级毛片在线看网站| 我要看黄色一级片免费的| 啦啦啦 在线观看视频| av天堂久久9| 精品视频人人做人人爽| 一级片'在线观看视频| 老熟女久久久| 黄色片一级片一级黄色片| 黄网站色视频无遮挡免费观看| 69av精品久久久久久 | 亚洲全国av大片| 欧美精品一区二区免费开放| 黄网站色视频无遮挡免费观看| 美女大奶头黄色视频| 我要看黄色一级片免费的| 90打野战视频偷拍视频| 少妇猛男粗大的猛烈进出视频| 日本精品一区二区三区蜜桃| 久久综合国产亚洲精品| 亚洲精品一二三| www.精华液| 免费少妇av软件| 国产成人一区二区三区免费视频网站| 91字幕亚洲| 成人黄色视频免费在线看| 欧美成人午夜精品| 亚洲激情五月婷婷啪啪| 亚洲午夜精品一区,二区,三区| 免费女性裸体啪啪无遮挡网站| 亚洲国产欧美网| 久久久水蜜桃国产精品网| 国产精品欧美亚洲77777| 免费黄频网站在线观看国产| 精品久久久久久电影网| 91字幕亚洲| 十八禁人妻一区二区| 一本久久精品| 欧美变态另类bdsm刘玥| 男女国产视频网站| 少妇精品久久久久久久| 日本撒尿小便嘘嘘汇集6| 国产精品国产三级国产专区5o| 老鸭窝网址在线观看| 热re99久久国产66热| 啦啦啦 在线观看视频| av天堂久久9| 久久毛片免费看一区二区三区| 久久这里只有精品19| 我要看黄色一级片免费的| 最近中文字幕2019免费版| 亚洲av电影在线观看一区二区三区| 嫩草影视91久久| 日韩有码中文字幕| 又紧又爽又黄一区二区| 亚洲av日韩在线播放| 又黄又粗又硬又大视频| 欧美激情极品国产一区二区三区| 欧美另类亚洲清纯唯美| 国产又爽黄色视频| 人妻一区二区av| 最新的欧美精品一区二区| 日韩国内少妇激情av| 身体一侧抽搐| 亚洲av熟女| 一级毛片高清免费大全| videosex国产| 亚洲人成网站高清观看| 亚洲18禁久久av| 日本在线视频免费播放| 亚洲全国av大片| 国产成人aa在线观看| 亚洲av五月六月丁香网| 久久中文看片网| 日本在线视频免费播放| 国产高清视频在线播放一区| 成人手机av| 免费无遮挡裸体视频| 日本免费a在线| 看片在线看免费视频| 久久精品影院6| 正在播放国产对白刺激| 美女免费视频网站| 亚洲国产中文字幕在线视频| 97人妻精品一区二区三区麻豆| tocl精华| 麻豆成人av在线观看| 岛国在线免费视频观看| 国产午夜精品论理片| 丝袜美腿诱惑在线| 欧美大码av| 亚洲国产看品久久| 老司机午夜福利在线观看视频| 亚洲 国产 在线| 美女 人体艺术 gogo| 欧美成人午夜精品| 一进一出好大好爽视频| 九色国产91popny在线| 国产一区在线观看成人免费| 嫁个100分男人电影在线观看| 全区人妻精品视频| 午夜福利在线观看吧| 日日爽夜夜爽网站| 国产欧美日韩精品亚洲av| 亚洲国产欧美人成| 成人三级做爰电影| 国产精品久久久av美女十八| 日本在线视频免费播放| 国产高清激情床上av| 99在线视频只有这里精品首页| 国产精品综合久久久久久久免费| 1024手机看黄色片| 婷婷精品国产亚洲av在线| 麻豆成人午夜福利视频| 老汉色∧v一级毛片| 国产精品美女特级片免费视频播放器 | 成人一区二区视频在线观看| 在线免费观看的www视频| 精品久久久久久成人av| 久久精品国产清高在天天线| 美女扒开内裤让男人捅视频| 老熟妇仑乱视频hdxx| 18禁黄网站禁片午夜丰满| 黑人巨大精品欧美一区二区mp4| 91成年电影在线观看| 欧美日韩国产亚洲二区| 精品久久久久久久毛片微露脸| 久久精品亚洲精品国产色婷小说| 久久久精品欧美日韩精品| 日韩大尺度精品在线看网址| 亚洲av日韩精品久久久久久密| 成人欧美大片| 国产蜜桃级精品一区二区三区| xxx96com| 午夜福利欧美成人| 久久久久精品国产欧美久久久| 一区福利在线观看| 在线观看免费午夜福利视频| 久久国产乱子伦精品免费另类| 99久久国产精品久久久| 一级a爱片免费观看的视频| 国产97色在线日韩免费| 国产三级在线视频| 午夜福利成人在线免费观看| 国产又色又爽无遮挡免费看| 久久热在线av| 国产主播在线观看一区二区| 色尼玛亚洲综合影院| 99热只有精品国产| 成人18禁在线播放| 久久久水蜜桃国产精品网| 亚洲成av人片在线播放无| 久久久久久大精品| 成人高潮视频无遮挡免费网站| 久久久久国内视频| 久热爱精品视频在线9| 久久午夜亚洲精品久久| 亚洲一区二区三区色噜噜| 国产伦人伦偷精品视频| 欧美中文综合在线视频| 亚洲,欧美精品.| 国产探花在线观看一区二区| 国产精品精品国产色婷婷| 国内精品一区二区在线观看| 日韩av在线大香蕉| 午夜两性在线视频| 欧美乱妇无乱码| 国产午夜福利久久久久久| 久久天躁狠狠躁夜夜2o2o| 亚洲专区字幕在线| 亚洲精品粉嫩美女一区| 黄色丝袜av网址大全| 国产精品综合久久久久久久免费| 我要搜黄色片| 99热这里只有是精品50| 老熟妇乱子伦视频在线观看| 国产精品一及| 欧美大码av| 日本五十路高清| 国产精品久久视频播放| 亚洲欧美精品综合久久99| 亚洲精品在线美女| 国产一区二区在线av高清观看| 亚洲男人的天堂狠狠| 欧美一区二区精品小视频在线| 亚洲成av人片免费观看| 又黄又粗又硬又大视频| 日本精品一区二区三区蜜桃| 91麻豆av在线| 老司机深夜福利视频在线观看| 亚洲男人天堂网一区| 欧美成人免费av一区二区三区| 精品国产超薄肉色丝袜足j| 99精品欧美一区二区三区四区| 久久午夜综合久久蜜桃| 少妇人妻一区二区三区视频| 毛片女人毛片| 特级一级黄色大片| 中文字幕熟女人妻在线| 午夜久久久久精精品| 亚洲欧美日韩无卡精品| 欧美日韩一级在线毛片| 好男人电影高清在线观看| 99热6这里只有精品| 12—13女人毛片做爰片一| 欧美日韩福利视频一区二区| 波多野结衣高清无吗| 亚洲成a人片在线一区二区| 亚洲成人久久性| 午夜精品在线福利| 久久久国产欧美日韩av| www国产在线视频色| 日韩 欧美 亚洲 中文字幕| 特大巨黑吊av在线直播| 长腿黑丝高跟| 亚洲av中文字字幕乱码综合| 国产一区二区三区在线臀色熟女| 搞女人的毛片| 欧美成狂野欧美在线观看| 欧美一区二区精品小视频在线| 亚洲欧美精品综合久久99| 亚洲国产看品久久| svipshipincom国产片| 色综合站精品国产| 国产精品美女特级片免费视频播放器 | 免费观看精品视频网站| 哪里可以看免费的av片| 亚洲中文日韩欧美视频| 在线观看美女被高潮喷水网站 | 99热6这里只有精品| 天天躁夜夜躁狠狠躁躁| 夜夜夜夜夜久久久久| 啦啦啦免费观看视频1| 午夜成年电影在线免费观看| 在线免费观看的www视频| 91九色精品人成在线观看| 性色av乱码一区二区三区2| 日本黄色视频三级网站网址| 最近最新免费中文字幕在线| 午夜两性在线视频| 久久久国产成人免费| 亚洲精品在线美女| 亚洲人成电影免费在线| 人人妻人人澡欧美一区二区| 一级作爱视频免费观看| 黄色 视频免费看| 成人三级做爰电影| 国产蜜桃级精品一区二区三区| 日本一本二区三区精品| 久久久久久久午夜电影| 女人爽到高潮嗷嗷叫在线视频| 日本a在线网址| 久久欧美精品欧美久久欧美| 丰满人妻熟妇乱又伦精品不卡| 午夜福利免费观看在线| 久久香蕉精品热| 欧美在线一区亚洲| 亚洲熟女毛片儿| 桃色一区二区三区在线观看| 日韩欧美精品v在线| 日韩大尺度精品在线看网址| 欧美在线一区亚洲| 特级一级黄色大片| 麻豆成人午夜福利视频| 久久精品aⅴ一区二区三区四区| 黄色片一级片一级黄色片| av欧美777| 色综合婷婷激情| 国产精品九九99| 在线观看66精品国产| 久久九九热精品免费| 国产爱豆传媒在线观看 | 不卡一级毛片| 午夜激情福利司机影院| 他把我摸到了高潮在线观看| 一卡2卡三卡四卡精品乱码亚洲| 欧美一区二区国产精品久久精品 | √禁漫天堂资源中文www| 90打野战视频偷拍视频| 亚洲人成网站高清观看| 两个人视频免费观看高清| 18禁观看日本| 精品久久蜜臀av无| 国产精品九九99| 午夜a级毛片| 成人18禁在线播放| 两个人的视频大全免费| 国产成人一区二区三区免费视频网站| 亚洲精品一区av在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲自拍偷在线| 亚洲 欧美一区二区三区| 国产精品av久久久久免费| 国产精品1区2区在线观看.| 亚洲精品久久国产高清桃花| 欧美一区二区精品小视频在线| 成人18禁在线播放| 十八禁网站免费在线| 好看av亚洲va欧美ⅴa在| 一级片免费观看大全| 亚洲精品中文字幕在线视频| 久久久国产成人免费| 日本黄大片高清| 国产精品久久电影中文字幕| 国产精品一区二区三区四区久久| 精品久久久久久久久久久久久|