• <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ǔ)
    亚洲激情五月婷婷啪啪| av女优亚洲男人天堂| 99热这里只有是精品在线观看| 成人美女网站在线观看视频| 亚洲电影在线观看av| 日韩不卡一区二区三区视频在线| 日本爱情动作片www.在线观看| 麻豆成人午夜福利视频| 日韩视频在线欧美| 久久久久久九九精品二区国产| 99久久中文字幕三级久久日本| 在线观看免费高清a一片| 婷婷色av中文字幕| 99国产精品免费福利视频| 国产黄色免费在线视频| 欧美性感艳星| 极品教师在线视频| 亚洲中文av在线| 亚洲熟女精品中文字幕| 亚洲aⅴ乱码一区二区在线播放| 免费看av在线观看网站| 亚洲国产精品999| 久久久亚洲精品成人影院| 777米奇影视久久| 午夜福利影视在线免费观看| 欧美日韩亚洲高清精品| 尾随美女入室| 亚洲精品日韩av片在线观看| 91久久精品国产一区二区成人| 干丝袜人妻中文字幕| 国产精品人妻久久久久久| 两个人的视频大全免费| 高清av免费在线| 一区二区三区乱码不卡18| 蜜桃久久精品国产亚洲av| 日韩,欧美,国产一区二区三区| 亚洲人成网站在线观看播放| 男女啪啪激烈高潮av片| 欧美少妇被猛烈插入视频| 在线免费十八禁| 国产色婷婷99| 久久热精品热| 丰满迷人的少妇在线观看| 国产中年淑女户外野战色| av专区在线播放| 国产av码专区亚洲av| 黄片wwwwww| 国产黄片视频在线免费观看| 蜜桃亚洲精品一区二区三区| 免费大片黄手机在线观看| 亚洲成人中文字幕在线播放| av卡一久久| 欧美极品一区二区三区四区| 久久久久久久久久人人人人人人| 日韩中文字幕视频在线看片 | 国产精品99久久久久久久久| 国产一区二区三区av在线| 欧美一区二区亚洲| 国产一级毛片在线| 国产男女超爽视频在线观看| xxx大片免费视频| 亚洲av免费高清在线观看| 99热这里只有是精品50| 亚洲av日韩在线播放| 国产男人的电影天堂91| 下体分泌物呈黄色| 日韩欧美一区视频在线观看 | 免费久久久久久久精品成人欧美视频 | 成年美女黄网站色视频大全免费 | 久久国产乱子免费精品| 好男人视频免费观看在线| 久久精品国产亚洲av天美| 国产成人精品福利久久| 国产精品一及| 国精品久久久久久国模美| 五月天丁香电影| 久久久久久久国产电影| 国产午夜精品一二区理论片| av在线观看视频网站免费| 欧美性感艳星| 精品99又大又爽又粗少妇毛片| 成人免费观看视频高清| 日韩三级伦理在线观看| 下体分泌物呈黄色| 色5月婷婷丁香| 街头女战士在线观看网站| 天堂8中文在线网| 免费看日本二区| 亚洲欧美日韩东京热| av国产免费在线观看| 国产成人免费无遮挡视频| 成年av动漫网址| 成人亚洲欧美一区二区av| 99热网站在线观看| 在线免费十八禁| 国产精品国产av在线观看| 午夜福利高清视频| 人人妻人人添人人爽欧美一区卜 | 2022亚洲国产成人精品| 啦啦啦中文免费视频观看日本| 国产v大片淫在线免费观看| 国产在线免费精品| 九九久久精品国产亚洲av麻豆| 亚洲精品国产av成人精品| 精品久久久噜噜| 国产爱豆传媒在线观看| 日本av手机在线免费观看| 精品99又大又爽又粗少妇毛片| 老司机影院毛片| 国产精品国产三级国产av玫瑰| 国产亚洲精品久久久com| 国产成人精品一,二区| 成人毛片60女人毛片免费| 一级毛片电影观看| 超碰97精品在线观看| 亚洲欧美日韩东京热| 国产精品av视频在线免费观看| av国产久精品久网站免费入址| 国产一区有黄有色的免费视频| 三级国产精品欧美在线观看| 在线观看美女被高潮喷水网站| 国产v大片淫在线免费观看| 亚洲成色77777| 一本—道久久a久久精品蜜桃钙片| 国产精品嫩草影院av在线观看| 亚洲丝袜综合中文字幕| 在线观看免费日韩欧美大片 | 舔av片在线| 男人狂女人下面高潮的视频| 人人妻人人澡人人爽人人夜夜| 丰满人妻一区二区三区视频av| 下体分泌物呈黄色| 欧美+日韩+精品| 日韩不卡一区二区三区视频在线| 99久久精品国产国产毛片| 亚洲丝袜综合中文字幕| 观看av在线不卡| 亚洲婷婷狠狠爱综合网| 国产日韩欧美在线精品| 91aial.com中文字幕在线观看| 亚洲aⅴ乱码一区二区在线播放| 免费黄频网站在线观看国产| 自拍偷自拍亚洲精品老妇| 国产高清有码在线观看视频| 黄色欧美视频在线观看| 国产 一区精品| 亚洲高清免费不卡视频| 国产精品人妻久久久久久| 最黄视频免费看| 欧美日韩国产mv在线观看视频 | 黄色怎么调成土黄色| 老熟女久久久| 成人毛片60女人毛片免费| 成人漫画全彩无遮挡| 中文字幕制服av| 亚洲成色77777| 国产成人精品婷婷| 男人爽女人下面视频在线观看| 伊人久久国产一区二区| 国产永久视频网站| 各种免费的搞黄视频| 91精品国产国语对白视频| 美女视频免费永久观看网站| 永久网站在线| 偷拍熟女少妇极品色| 久久久久久人妻| 永久网站在线| 99九九线精品视频在线观看视频| 精品亚洲成国产av| 国产精品无大码| 亚洲国产最新在线播放| 人妻 亚洲 视频| 亚洲熟女精品中文字幕| 黄色视频在线播放观看不卡| 国产在线视频一区二区| 国产精品三级大全| 内地一区二区视频在线| 交换朋友夫妻互换小说| 好男人视频免费观看在线| 午夜免费观看性视频| 一区在线观看完整版| 能在线免费看毛片的网站| 80岁老熟妇乱子伦牲交| 七月丁香在线播放| 免费看av在线观看网站| 男的添女的下面高潮视频| 亚洲精品久久午夜乱码| 久热这里只有精品99| 26uuu在线亚洲综合色| 国产精品久久久久成人av| 欧美日本视频| 亚洲av成人精品一二三区| av女优亚洲男人天堂| 99国产精品免费福利视频| 免费不卡的大黄色大毛片视频在线观看| 国产黄色视频一区二区在线观看| 亚洲成人av在线免费| 亚洲国产欧美人成| 精品亚洲成国产av| 五月伊人婷婷丁香| 国产美女午夜福利| 精品国产露脸久久av麻豆| 国产精品欧美亚洲77777| 国产欧美日韩一区二区三区在线 | 国产熟女欧美一区二区| 18禁在线播放成人免费| 成人综合一区亚洲| 亚洲国产精品专区欧美| 三级国产精品片| 欧美区成人在线视频| 狂野欧美激情性xxxx在线观看| 成年av动漫网址| 午夜福利高清视频| 99久久综合免费| 午夜视频国产福利| 97超视频在线观看视频| 成年免费大片在线观看| 午夜福利高清视频| 大香蕉97超碰在线| 性高湖久久久久久久久免费观看| 免费久久久久久久精品成人欧美视频 | 尤物成人国产欧美一区二区三区| 久久精品国产亚洲av天美| 国产精品99久久久久久久久| 两个人的视频大全免费| 久久久久视频综合| 国语对白做爰xxxⅹ性视频网站| 国产爽快片一区二区三区| 国产成人freesex在线| 成人黄色视频免费在线看| 亚洲精品久久午夜乱码| 超碰97精品在线观看| 国产在线视频一区二区| 五月伊人婷婷丁香| 精品少妇久久久久久888优播| 欧美3d第一页| 91狼人影院| 街头女战士在线观看网站| 久久人人爽av亚洲精品天堂 | 天堂中文最新版在线下载| 亚洲国产精品专区欧美| 日韩免费高清中文字幕av| 在线观看人妻少妇| 国产成人精品久久久久久| 自拍欧美九色日韩亚洲蝌蚪91 | 人妻系列 视频| 久久亚洲国产成人精品v| 亚洲欧美日韩无卡精品| 青青草视频在线视频观看| 成年av动漫网址| 欧美+日韩+精品| 成人午夜精彩视频在线观看| 91久久精品电影网| av一本久久久久| 日韩一区二区三区影片| 高清视频免费观看一区二区| 街头女战士在线观看网站| 久久97久久精品| 免费黄频网站在线观看国产| 我的老师免费观看完整版| 日韩强制内射视频| 国产精品伦人一区二区| 亚洲,欧美,日韩| 3wmmmm亚洲av在线观看| 性色avwww在线观看| 成人午夜精彩视频在线观看| 国产黄频视频在线观看| 久久99热这里只频精品6学生| 国产精品国产三级国产专区5o| 免费看光身美女| 日韩制服骚丝袜av| 免费观看的影片在线观看| videos熟女内射| 精品人妻偷拍中文字幕| a级毛片免费高清观看在线播放| 国产精品久久久久成人av| 成人一区二区视频在线观看| 男女啪啪激烈高潮av片| 天堂8中文在线网| 国产爽快片一区二区三区| 欧美区成人在线视频| 麻豆国产97在线/欧美| 97精品久久久久久久久久精品| 亚洲av在线观看美女高潮| 国产有黄有色有爽视频| 欧美少妇被猛烈插入视频| 香蕉精品网在线| 亚洲成色77777| 国产成人精品婷婷| 亚洲精品456在线播放app| 亚洲色图av天堂| 黑人高潮一二区| av女优亚洲男人天堂| 欧美日本视频| 蜜臀久久99精品久久宅男| 国产高清有码在线观看视频| 国产大屁股一区二区在线视频| 成人毛片60女人毛片免费| 青春草亚洲视频在线观看| 欧美日韩亚洲高清精品| 草草在线视频免费看| 青春草亚洲视频在线观看| 国产亚洲一区二区精品| 久久久久国产网址| 一本—道久久a久久精品蜜桃钙片| 下体分泌物呈黄色| 欧美日本视频| 少妇人妻久久综合中文| 国产精品熟女久久久久浪| av视频免费观看在线观看| 国内揄拍国产精品人妻在线| 最近手机中文字幕大全| 97超碰精品成人国产| 在线观看一区二区三区| 国产精品熟女久久久久浪| 国产伦在线观看视频一区| 亚州av有码| 美女视频免费永久观看网站| 亚洲国产高清在线一区二区三| 亚洲精品国产成人久久av| 久热久热在线精品观看| 国产av精品麻豆| 中文精品一卡2卡3卡4更新| 中文乱码字字幕精品一区二区三区| 日韩在线高清观看一区二区三区| 久久精品久久久久久噜噜老黄| av女优亚洲男人天堂| 伊人久久精品亚洲午夜| 亚洲自偷自拍三级| 久久久久久久大尺度免费视频| 亚洲欧美一区二区三区黑人 | 亚洲人与动物交配视频| 日本欧美视频一区| 欧美最新免费一区二区三区| 亚洲天堂av无毛| 六月丁香七月| 嘟嘟电影网在线观看| 91狼人影院| 在线观看免费日韩欧美大片 | 赤兔流量卡办理| 在线看a的网站| 女人十人毛片免费观看3o分钟| 日韩欧美一区视频在线观看 | 日韩av免费高清视频| 久久精品夜色国产| 久久这里有精品视频免费| 欧美xxⅹ黑人| 色吧在线观看| 热99国产精品久久久久久7| 在线观看免费视频网站a站| 你懂的网址亚洲精品在线观看| 九九久久精品国产亚洲av麻豆| 免费观看无遮挡的男女| 久久精品熟女亚洲av麻豆精品| av免费观看日本| 人妻少妇偷人精品九色| 国产成人freesex在线| 欧美国产精品一级二级三级 | 欧美变态另类bdsm刘玥| 又大又黄又爽视频免费| 精品人妻熟女av久视频| 国产毛片在线视频| 日韩大片免费观看网站| 国产精品国产三级国产专区5o| 成人亚洲欧美一区二区av| 99九九线精品视频在线观看视频| videos熟女内射| 日本与韩国留学比较| 插阴视频在线观看视频| 另类亚洲欧美激情| 国产 一区精品| 美女cb高潮喷水在线观看| 久久久欧美国产精品| 2022亚洲国产成人精品| 99九九线精品视频在线观看视频| 国产一级毛片在线| 成人毛片60女人毛片免费| 卡戴珊不雅视频在线播放| 日韩成人伦理影院| 我的女老师完整版在线观看| 日韩一区二区三区影片| 肉色欧美久久久久久久蜜桃| 中国国产av一级| 国产色婷婷99| 男女下面进入的视频免费午夜| 最近2019中文字幕mv第一页| 欧美区成人在线视频| 天美传媒精品一区二区| 亚洲av成人精品一区久久| 一区二区三区精品91| 春色校园在线视频观看| 十八禁网站网址无遮挡 | 各种免费的搞黄视频| 一级av片app| 纵有疾风起免费观看全集完整版| 美女福利国产在线 | 精品人妻视频免费看| 毛片女人毛片| 女人久久www免费人成看片| 舔av片在线| 日本欧美视频一区| 精华霜和精华液先用哪个| 性色avwww在线观看| 欧美+日韩+精品| 国国产精品蜜臀av免费| 国产精品麻豆人妻色哟哟久久| 午夜福利在线在线| 天天躁日日操中文字幕| 亚洲国产精品999| 国产精品人妻久久久影院| 久热久热在线精品观看| 在线观看av片永久免费下载| 18禁在线无遮挡免费观看视频| 国产毛片在线视频| 国产伦精品一区二区三区视频9| 欧美xxxx性猛交bbbb| 国产精品一区二区在线观看99| 偷拍熟女少妇极品色| 夜夜爽夜夜爽视频| 久久久久精品久久久久真实原创| 在现免费观看毛片| 国产在线免费精品| 中国国产av一级| 亚洲第一av免费看| 欧美日韩在线观看h| 亚洲欧美清纯卡通| 国产伦理片在线播放av一区| 日韩三级伦理在线观看| 国产欧美日韩一区二区三区在线 | 卡戴珊不雅视频在线播放| 在线观看免费视频网站a站| 日日撸夜夜添| 麻豆精品久久久久久蜜桃| 1000部很黄的大片| 婷婷色综合www| 亚洲成人av在线免费| 亚洲,一卡二卡三卡| 国产精品99久久久久久久久| 国产成人精品婷婷| 久久精品国产亚洲av天美| 香蕉精品网在线| 色网站视频免费| 一二三四中文在线观看免费高清| 人妻一区二区av| 国产伦在线观看视频一区| 欧美变态另类bdsm刘玥| 日本色播在线视频| 亚洲精品成人av观看孕妇| 男女免费视频国产| 久久午夜福利片| 久久精品国产亚洲av涩爱| 国国产精品蜜臀av免费| 伊人久久国产一区二区| 在线观看免费日韩欧美大片 | 九草在线视频观看| 人人妻人人添人人爽欧美一区卜 | 一区在线观看完整版| 亚洲精品第二区| 天天躁夜夜躁狠狠久久av| 99久久精品国产国产毛片| 高清午夜精品一区二区三区| 六月丁香七月| 五月开心婷婷网| 黄色怎么调成土黄色| 免费看不卡的av| 2018国产大陆天天弄谢| 国产精品麻豆人妻色哟哟久久| 精品久久久久久久久亚洲| 亚洲精品456在线播放app| 免费高清在线观看视频在线观看| 国产亚洲精品久久久com| 一个人看的www免费观看视频| 欧美一区二区亚洲| 久久婷婷青草| 亚州av有码| av在线老鸭窝| 精品久久久久久电影网| 91精品国产九色| 亚洲精品自拍成人| a 毛片基地| 纵有疾风起免费观看全集完整版| 在线 av 中文字幕| 亚洲欧美精品自产自拍| 婷婷色综合www| av.在线天堂| 亚洲欧美清纯卡通| 人体艺术视频欧美日本| 一区二区三区免费毛片| 老司机影院毛片| av国产久精品久网站免费入址| 少妇高潮的动态图| 超碰av人人做人人爽久久| 老女人水多毛片| 中国三级夫妇交换| 亚洲av在线观看美女高潮| 精华霜和精华液先用哪个| 韩国高清视频一区二区三区| 国产黄色免费在线视频| 亚洲欧美中文字幕日韩二区| 制服丝袜香蕉在线| 2021少妇久久久久久久久久久| 老司机影院毛片| 女人十人毛片免费观看3o分钟| 看免费成人av毛片| 日本爱情动作片www.在线观看| 尾随美女入室| 激情 狠狠 欧美| 久久这里有精品视频免费| 国产亚洲午夜精品一区二区久久| 丝袜喷水一区| 精品久久久久久电影网| 国产在线一区二区三区精| 久久久成人免费电影| 亚洲美女搞黄在线观看| 国产伦精品一区二区三区视频9| 国产欧美日韩一区二区三区在线 | 免费看日本二区| 国产极品天堂在线| av女优亚洲男人天堂| 全区人妻精品视频| 亚洲精品国产av成人精品| 国产av码专区亚洲av| 国产精品久久久久久av不卡| 女人久久www免费人成看片| 亚洲高清免费不卡视频| 大陆偷拍与自拍| 性色avwww在线观看| 欧美日韩一区二区视频在线观看视频在线| 国语对白做爰xxxⅹ性视频网站| 亚洲国产高清在线一区二区三| 99久国产av精品国产电影| 天堂中文最新版在线下载| 日本欧美国产在线视频| 成人毛片a级毛片在线播放| 亚洲国产色片| av网站免费在线观看视频| 亚洲欧美中文字幕日韩二区| 边亲边吃奶的免费视频| 久久久久精品性色| 亚洲欧美精品自产自拍| 亚洲精品国产av成人精品| 天堂中文最新版在线下载| av不卡在线播放| 午夜福利在线观看免费完整高清在| 我要看黄色一级片免费的| 在线看a的网站| 女性被躁到高潮视频| 国产免费福利视频在线观看| 国产精品久久久久成人av| 欧美日韩一区二区视频在线观看视频在线| 99久久综合免费| 国产精品久久久久久久电影| 最后的刺客免费高清国语| 插逼视频在线观看| 日韩欧美精品免费久久| 国产熟女欧美一区二区| 国产亚洲精品久久久com| 午夜福利高清视频| 午夜激情福利司机影院| 国产亚洲5aaaaa淫片| 免费看av在线观看网站| 最新中文字幕久久久久| 中文字幕人妻熟人妻熟丝袜美| 七月丁香在线播放| 又黄又爽又刺激的免费视频.| 国产亚洲91精品色在线| 一级毛片aaaaaa免费看小| 欧美少妇被猛烈插入视频| 中国美白少妇内射xxxbb| 青青草视频在线视频观看| 国产成人精品福利久久| 91精品国产国语对白视频| 国产精品.久久久| 最近手机中文字幕大全| 91久久精品电影网| 国产精品一二三区在线看| 免费观看a级毛片全部| 欧美高清性xxxxhd video| 亚洲国产精品专区欧美| 日韩视频在线欧美| 在线观看国产h片| 成人高潮视频无遮挡免费网站| 王馨瑶露胸无遮挡在线观看| 免费观看在线日韩| 狂野欧美白嫩少妇大欣赏| 免费观看a级毛片全部| 日日啪夜夜爽| 啦啦啦在线观看免费高清www| av视频免费观看在线观看| 久久综合国产亚洲精品| h日本视频在线播放| 久久国内精品自在自线图片| 1000部很黄的大片| 亚洲欧美日韩另类电影网站 | 狂野欧美激情性xxxx在线观看| 国产毛片在线视频| 国产免费视频播放在线视频| 国产一区二区三区综合在线观看 | 老师上课跳d突然被开到最大视频| 欧美精品国产亚洲| 亚洲性久久影院| 精品久久久久久久末码| 国产精品一区二区性色av| 美女高潮的动态| 色婷婷av一区二区三区视频| 亚州av有码| av专区在线播放| 午夜免费鲁丝| 夜夜看夜夜爽夜夜摸| 秋霞在线观看毛片| 最黄视频免费看| tube8黄色片| 国产精品三级大全| 在线精品无人区一区二区三 | 看免费成人av毛片|