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

    Matching characteristics of magnetic field configuration and chamfered channel wall in a magnetically shielded Hall thruster

    2021-10-31 08:14:56ZhaoyuWANG王昭宇HongLI李鴻ChaoZHONG鐘超YanlinHU扈延林YongjieDING丁永杰LiqiuWEI魏立秋andDarenYU于達(dá)仁
    Plasma Science and Technology 2021年10期

    Zhaoyu WANG (王昭宇), Hong LI (李鴻),2, Chao ZHONG (鐘超),Yanlin HU (扈延林), Yongjie DING (丁永杰),2, Liqiu WEI (魏立秋),2 and Daren YU (于達(dá)仁),2

    1 Lab of Plasma Propulsion, Harbin Institute of Technology, Harbin 150001, People’s Republic of China

    2 Key Laboratory of Aerospace Plasma Propulsion, Ministry of Industry and Information Technology,Harbin 150001, People’s Republic of China

    3 Beijing Institute of Control Engineering, Beijing 100190, People’s Republic of China

    Abstract To date,the selection of the magnetic field line used to match the chamfered inner and outer channel walls in a magnetically shielded Hall thruster has not been quantitatively studied.Hence, an experimental study was conducted on a 1.35 kW magnetically shielded Hall thruster with a xenon propellant.Different magnetic field lines were chosen, and corresponding tangentially matched channel walls were manufactured and utilized.The results demonstrate that high performance and a qualified anti-sputtering effect cannot be achieved simultaneously.When the magnetic field lines that match the chamfered wall have a strength at the channel centerline of less than 12% of the maximum field strength, the channel wall can be adequately protected from ion sputtering.When the magnetic field lines have a strength ratio of 12%–20%,the thruster performance is high.These findings provide the first significant quantitative design reference for the match between the magnetic field line and chamfered channel wall in magnetically shielded Hall thrusters.

    Keywords: magnetically shielded Hall thruster, magnetic field line, chamfered channel wall,tangent matching

    1.Introduction

    With the rapid development of space science and technology in recent years [1], space missions and research have shown significant changes in their long-term characteristics.Commercial aerospace is rapidly developing and has attracted significant attention.Commercial and civil-oriented low-orbit communications constellations, Internet of Things constellations, and other space missions require long-term operations in orbit[2].Excellent development opportunities and adverse application-related challenges both contribute to demanding requirements for the working life of Hall thrusters [3–5], which are expected to achieve a cumulative working time of over 150 000 h [6].Therefore,developing long-life Hall thrusters will be imperative in the future.

    The most critical factor that influences the lifespan of a Hall thruster is the sputtering erosion of the discharge channel wall caused by high-energy ions[7–9].In addition to maintaining and participating in the discharge process in the Hall thruster[10–14],another important role of the channel is to protect the magnetic circuit.When the channel wall is eroded through, high-energy ions in the acceleration zone bombard and erode the magnetic circuit.Once the components of the magnetic circuit are eroded and damaged,further degradation or overheating may occur,and the magnetic field distribution may deviate from the designed configuration and may not be able to maintain a stable discharge.This would ultimately lead to thruster failure.

    NASA’s Jet Propulsion Laboratory (JPL) conducted life tests for BPT-4000 in 2012[15,16].They found that the wall of BPT-4000 remained unchanged from 5600 to 10 400 h.Since then, the magnetic shielding theory of anti-sputtering has been refined based on simulations and experiments.The magnetic field lines of the magnetically shielded field have been made concave toward the anode and almost parallel to the channel wall.In this magnetic field configuration,the high potential and low electron temperature near the anode can be maintained at the wall by isothermal and equipotential magnetic field lines [17–19].The electric field parallel to the direction of the magnetic fieldE‖is negligible,and the electric field is almost perpendicular to the magnetic field.The potential near the wall is as high as the potential of the channel center.Moreover, the wall has a low electron temperature and low potential drop of the sheath, and the electric field along the wall being directed to the channel center effectively reduces the probability of ion bombardment on the wall [20].Compared with those of an unshielded magnetic field,the position of the maximum magnetic field strength and the ionization zone constrained by the magnetic field both move downstream, moving the acceleration zone further outside the channel.Because high-energy ions are produced outside the channel owing to the outward shift of the acceleration zone,wall erosion caused by high-energy ions can be avoided.It has been found that in a magnetically shielded Hall thruster, a wall profile that is parallel to the magnetic field lines can protect the wall from ion bombardment[21, 22].A chamfered wall in a 50 kW high-power magnetically shielded Hall thruster reduced the interaction intensity between the plasma and wall [23, 24], and the wall was protected from ion bombardment.The chamfered wall in a 1.35 kW magnetically shielded Hall thruster underwent almost no changes after an experiment [25].However, there are no institutional quantitative design references, and the knowledge regarding magnetic field line selection is insufficient.In this study, the matching characteristics of the magnetic field and wall chamfer were investigated by choosing different magnetic field lines, and the first quantitative design reference for the wall chamfer was obtained.

    To ensure that the thruster can prevent the wall from sputtering, the electron temperature at the wall must be low and the wall cannot hinder the magnetic field lines.By studying the discharge characteristics of the thruster considering different matched walls, the change regularity in performance and the anti-sputtering effect with the change in magnetic field lines can be determined.Finally, a design reference for magnetic field line selection for wall profiles in magnetically shielded thrusters with high performance can be established.This will serve as a basis for the wall design of magnetically shielded Hall thrusters.

    The remainder of this paper is organized as follows.Section 2 introduces the experimental methodology, and section 3 analyzes the change in the anti-sputtering effect and discharge characteristics under the wall profiles matching different magnetic field lines and the mechanism.The final section summarizes the paper.

    2.Research methodology and schemes

    2.1.Experimental design of matched wall profiles for antisputtering effect

    For a magnetically shielded Hall thruster, when the magnetic field and the wall are matched, the wall does not hinder the magnetic field lines and is equipotential.The magnetic field lines graze the wall, and the intersection points of the magnetic field lines and wall are the start and end positions of the chamfer, as shown in figure 1.In this study, different matched wall profiles were designed according to the magnetic field lines.

    A 1.35 kW magnetically shielded Hall thruster with inner and outer diameters of 70 mm and 100 mm,respectively,was independently developed by the Harbin Institute of Technology, as shown in figure 2.One inner and one outer magnetic coil were used for the magnetic field, providing satisfactory circumferential symmetry.The operation flow rate range of the thruster is 50–56 sccm,and the power ranges from 1.3 to 1.5 kW.The characteristics of the magnetic field can be expressed in detail using the FEMM two-dimensional axisymmetric finite element analysis software program.

    To quantitatively and conveniently designate and distinguish the magnetic field line used to match the chamfered channel wall, the characteristic parameteris defined as

    whereBcis the magnetic field strength at the intersection of the magnetic field line and the channel centerline,andBc,maxis the maximum magnetic field strength on the channel centerline.is dimensionless and describes the ratio of the magnetic field strength.Magnetic field lines withvalues of 36%, 28%, 20%, 12%, and 5% were selected using FEMM,as shown in figures 3(a) and (b), and the corresponding tangent wall profiles were designed and referred to as Cases 1,2,3,4,and 5,respectively.A straight channel,denoted as Case 0,was designed for comparison.A wall matched with a smallerhas a larger flow area and a larger chamfer.Table 1 lists the details of the chamfer size;the radial length represents the radial distance between the exit of the channel and the start of the chamfer, and θ is the angle between the chamfered wall and axial axis.

    Table 1.Chamfer size specifications.

    Table 2.Energy spectrum analysis results.

    The discharge characteristics of the thruster were analyzed using different magnetic field lines under constant flow rate and constant power conditions.The wall states were observed, and the wall profiles were measured after approximately 17 h of ignition.The relationship between the deposition and erosion rates can be determined by comparing the wall colors before and after the experiment.Before the experiment, the inner and outer walls were the original white of the boron nitride ceramic.When the walls were black after the experiment,the erosion rate of the ions on the wall was negligible during discharge.These black deposits consist of various different products formed by discharging, including sputtering of the thruster shell and vacuum chamber, as well as evaporation of the insulation layer of the magnetic coil due to the high temperature.If the thruster discharges normally when the wall has black depositions, the life of the thruster can be significantly extended.When the wall was white after the experiment,the sputtering rate was non-negligible, and the white area was the sputtering zone.The existence of a sputtering zone indicates that the life of the thruster is threatened.

    2.2.Experimental equipment

    The experiments were performed at the Plasma Propulsion Laboratory of the Harbin Institute of Technology.The vacuum chamber has a size of ?1.5 m×4 m and is evacuated by two oil-diffusion pumps, three mechanical booster pumps, and one rotary pump.The ultimate vacuum pressure reached 5×10?3Pa.When the thruster operated at a volume flow rate of 50 sccm,the vacuum pressure was below 5×10?3Pa.The uncertainty of the facility pressure was below 10%[26,27].The uncertainty of the gas flow controller was approximately 1%,and the purity of the propellant was 99.9995%.

    The thrust (T) was measured using a three-wire torsion pendulum device [28], which was a spring–mass–damper system [29].The thruster was placed on a suspended pallet,and the mass of the counterweight was adjusted to level the pallet.The generated thrust rotated the pallet, causing the angle of the laser incident on the mirror to change.After the laser was reflected, the position on the scale moved.The propulsion produced by the thruster can be calibrated based on the laser displacement produced by the weight of the known mass as a reference.Based on the discharge current(Id), thrust (T), discharge voltage (Ud), and anode flow (ma.),the thruster anode efficiency (ηa) can be calculated as follows [30]:

    A Faraday probe was used to measure the change in the ion beam current.The probe was mounted on a circular arc scanning rail, and the center of rotation coincided with the center of the thruster exit plane.The distance between the Faraday collector and center of rotation was 300 mm,the diameter of the collector was 9.5 mm, and the gap between the collector and the guard ring was 0.5 mm,which was less than 10 times the Debye length[31,32].Figure 4 shows the structure of the Faraday probe.Parallel capacitor and resistor connections were adopted to eliminate noise (R=1 kΩ,C=3 μF).

    Figure 1.Matching of the magnetic field lines and wall.

    Figure 2.1.35 kW magnetically shielded Hall thruster.

    Figure 3.Selected magnetic field lines with different in (a) axial distribution and (b) planar distribution .

    Figure 4.Structure of the Faraday probe.

    Figure 5.Measurement of wall profile using CMM.

    Figure 6.Diagram of the placement of tantalum sheet.

    Figure 7.Inner channel walls after the experiments.

    During the measurements, both collector and shield electrodes were provided with a bias voltage of ?24 V.The total ion current (Ii) was calculated as

    wherelis the distance between the probe and thruster center.

    The plume divergence angle was obtained by determining the ion current density distribution (θj( )) in the plume area.With the center of the thruster as the center, the plume divergence half-angle (α) was calculated as follows:

    The thruster wall profile measurement was commissioned by Harbin Ultra Precision Equipment Engineering and Technology Center Co., Ltd.The measurements were performed using a coordinate measuring machine (CMM) (Leitz, Germany), as shown in figure 5.The CMM model was a PMMCUtra12107.The contours of the four circumferential positions,i.e.0°, 90°, 180°, and 270°, at both the inner and outer walls were measured before and after the experiment (the circumferential azimuth of the cathode was defined as 0°).Then, the change in wall profile was determined by calculating the average difference between the values before and after the experiment.Unshielded Hall thrusters, such as the SPT-100, PPS1350, and H6US, have a minimum prophase erosion rate of 10 μm h?1[33], while those of magnetically shielded thrusters are 2–3 orders of magnitude lower; the H6MS and MSHT-600 have maximum erosion rates of 0.1 and 0.08 μm h?1, respectively[34].The measurement accuracy of the CMM is of the submicron order.To accurately determine the change in wall profile,the thruster was ignited for approximately 17 h in each case.

    A tantalum sheet was used to collect the deposits during the experiment.The sheet was placed at the thruster exit plane,located outside the outer pole and radially from the channel centerline, as shown in figure 6.The deposition spectrum was analyzed using a high-resolution fully automatic field-emission scanning electron microscope (FE-SEMSU5000).

    3.Experimental results and analysis

    3.1.Anti-sputtering effects

    By observing the wall color after a long ignition time under standard conditions (300 V, 50 sccm) in figure 7, we analyzed the thruster’s ability to protect the wall against ion sputtering.For Cases 4 and 5,it can be seen that when the walls match the magnetic field lines withvalues of 12% and 5%, black depositions appear on the chamfered walls, and the thruster is qualified with the anti-sputtering effect.With a reduction in chamfer area,part of the chamfer in Case 3 appears gray,and the anti-sputtering effect is poor.Prominent sputtering bands are observed on the straight section in the channels for Cases 1 and 2;this indicates that erosion could not be prevented.In addition,a 2–3 mm white sputtering band was observed in Case 0.

    The energy spectrum analysis results for the depositions are presented in table 2.The weight percentage of C in the deposition was 62.81%.The main composition of the depositions indicates that the depositions are conductive.

    The differences in the wall profile after the experiments were determined using the CMM, as shown in figure 8.For Cases 1 and 3,the relative change in the outer wall is only 1/3 of that in the inner wall,indicating that the outer wall has a better anti-sputtering effect; moreover, the outer wall in Case 5 is thickened after the experiment.The wall profile changes the most in Case 0;the erosion of the inner wall significantly decreases with increasing chamfer area.The prophase erosion rate for Case 0 reached 2.5 μm h?1.For the unshielded thruster, the erosion rate in the straight channel was greater than 10 μm h?1, which indicates that the application of the magnetically shielded magnetic field significantly reduced the wall erosion.The erosion of the inner wall in Case 5 was too small to be measured, the average difference in the profile near the exit was approximately 1.8 μm, and the prophase erosion rate was estimated to be 0.106 μm h?1, which is the same as the value determined by JPL.Figure 8 shows that the erosion in Case 0 continuously varied along the axial direction and that there was a dramatic change in erosion rate when the matched wall chamfer could not maintain the anode potential and low electron temperature.For the inner wall,drastic drops occurred near the beginning of the chamfer,indicating that the erosion was severe downstream at the beginning of the chamfer, as shown in figure 8.For Cases 1 and 3, substantial changes were observed at 23.6 mm and 22.8 mm, respectively.

    Figure 8.Changes in profiles of (a) inner and (b) outer channel walls after the experiments.

    Figure 9.Sketch of electron motion in the matching area of magnetic field and chamfered wall.

    Figure 10.Atom density along the channel centerline.

    When the chamfer area is reduced,although the matching magnetic field line strength is sufficient to restrain electrons(as indicated by the green line in figure 9), there exist magnetic field lines with lowervalues that can also restrain electrons, causing the field line with the lowerto be hindered by the wall(as indicated by the yellow line in figure 9).Electrons restrained by yellow lines tend to hit the chamfer directly.In this case, the wall potential is reduced, and the electron temperature increases with theof the line.The decrease in chamfer potential increases the energy that ion gained through acceleration from the anode to the chamfer.The wall potential with a reduced chamfer failed to repel high-energy ion bombardment at the wall.In addition,because the electrons confined by the magnetic field line hit the wall, an electronic current parallel to the magnetic direction was generated, which further led to ion movement toward the wall and weakened the anti-sputtering effect.From the foregoing analysis,we found that reducingenhances the anti-sputtering effect.In this experiment, the thruster met the requirement for protection against ion sputtering when the wall matched the magnetic field lines withvalues under 12%.

    3.2.Discharge characteristics

    The channel is the boundary of the discharge chamber, and the chamfer area has a significant effect on the gas flow and density distribution.In addition, the channel material affects electron conduction.High-energy ions are ejected from the channel exit, which develop the force.Furthermore, the chamfer significantly affects the strong interaction between the plasma and wall and determines the energy loss.For different matching cases,the potential distribution of the wall and its vicinity determines the beam utilization.In this study,the influences of the chamfer on the neutral gas distribution,electron conduction, interaction between ions and walls, and spatial potential were analyzed to elucidate the relationship between the chamfer and discharge characteristics.

    The ratio of the mean free path of the xenon atom to the channel characteristic length in the Hall thruster is greater than 10, which satisfies the condition that the Knudsen coefficient exceeds 10.Therefore, the free molecular flow model can be used to simulate anode distribution [35].The atom density distributions of the channel centerline in the steady-state for the different cases are shown in figure 10.These distributions demonstrate that the area of the matched chamfer wall increases and that the atom density in the channel decreases with decreasingof the field line.The atomic density in Case 0 was much higher than that in the other cases.

    The ionization process of atoms in the channel can be expressed as

    wherenais the anode density,neis the electron density,σiis the ionization cross-section,veis the velocity of the electrons,and〈σive〉 is the mean ionization rate.Because the maximum magnetic field strength exists outside the channel in the magnetically shielded thruster, the ionization zone is at the channel exit,and the acceleration zone is outside the channel.As theof the field line decreases, a low atom density can reduce the ionization intensity.

    3.2.1.Results under the constant power condition.Figure 11 shows the performance characteristics of the thruster under a constant power of 1.5 kW (300 V).In all cases apart from Case 0,the anode flow rate decreases with a decrease in theB? of the field line.Case 3 exhibits a clear advantage in terms of performance.

    Figure 11.(a)Volume flow rate and thrust and(b)anode efficiency and specific impulse under different chamfered wall cases at a constant discharge power of 1.5 kW.

    Figure 12.(a)Discharge current and thrust and(b)anode efficiency and specific impulse under different chamfered wall cases at a constant anode flow rate of 50 sccm.

    Figure 13.Ion current density distribution under different chamfered wall cases.

    Figure 14.Normalized plume divergence angle.

    The flow rates in Cases 3 and 5 were 54 sccm and 52.5 sccm,respectively.Compared with that in Case 3,the flow rate in Case 5 decreased by approximately 2.8%;however, the thrust dropped from 97.2 to 91.4 mN,i.e.a reduction of 6.0%.The decrease in thrust is not proportional to the decrease in flow rate,and it can be inferred that the ionization intensity decreased from Case 3 to Case 5.Importantly, although the ionization intensity decreased, the flow rate declined continuously from Case 3 to Case 5 as well.This represents the deposition that accelerates the electron conduction to the anode and increases the electron current in the anode current in Case 5, which further contributes to the low current efficiency,low efficiency,and low specific impulse in Case 5.

    The energy loss at the wall is another important factor affecting the performance of the thruster.Case 0 is a straight channel, which has the highest atom density in the channel but stops the ions from scattering at the channel exit and causes energy loss.Therefore, the thrust and efficiency in Case 0 are 2.5% and 3.6% less than those in Case 1,respectively.

    The flow rate decreased by 2.7% from Case 1 to Case 3,whereas the thrust increased by 5.4%.The reduction in energy loss significantly improves performance and increases propellant utilization, contributing to the improvement in performance in Case 3.

    3.2.2.Results under the constant flow rate condition.Further experiments were performed at a constant flow rate of 50 sccm.The discharge voltage was 300 V.As shown in figure 12,the current increased gradually with decreasingof the field line (except for Case 0), while the thrust first increased and subsequently decreased.Case 4 had a high specific impulse, and Case 3 exhibited a high efficiency.Therefore, high performance is achieved when the wall matches the magnetic field lines with aof 12%–20%.

    Compared with those in Case 1, the discharge current in Case 0 was 4.8%higher,while the thrust was approximately 1%lower.Because many ions are lost at the wall,the low propellant utilization contributes to the low thrust in Case 0;this is because the energy loss at the wall is relative to the anti-sputtering effect,and Case 0 cannot prevent ions from sputtering.From Case 1 to Case 3, the discharge current changed from 4.15 to 4.32 A,an increase of 4%; however, the thrust increased from 80.4 to 87.53 mN, i.e.by 8.9%, indicating that the ions lost at the wall in Case 1 form the thrust in Cases 2 and 3.As theB?of the magnetic field lines decreases from Case 0 to Case 5,the energy loss gradually decreases.Therefore, the wall can protect itself from ion sputtering when it matches the magnetic field lines with a lowB,? which is consistent with the results in section 3.1.

    The discharge current continues to increase from Case 1 to Case 5.Since the depositions are conductive,the electrons accelerated by the depositions maintain a high discharge current.The discharge current increases by 2.9%from Case 3 to Case 4, whereas the thrust in Case 4 is only 0.3% higher than that in Case 3;we can infer from the experimental results that the increase in the electron current contributes to the increase in the discharge current and that the decrease in energy loss contributes to the increase in thrust.The rising gradient in the discharge current from Case 4 to Case 5 is lower than that from Case 3 to Case 4.Because the depositions in Case 4 are sufficiently firm to reduce the surface resistance, the increased deposition in Case 5 cannot significantly increase the discharge current, and the electron loss at the wall reduces the ionization intensity in Case 5;this causes the thrust to decrease in Case 5.

    A Faraday probe can be used to obtain the ion current and beam profile in different cases under constant flow rate conditions,as shown in figure 13.Case 3 has the highest peak ion current density, and the ion current obtained using the integral along the curve is 4.3 A.The ion current in Case 5 decreased by 2.33% compared with that in Case 3, which is consistent with the regularity of the thrust shown in figure 12.The thrust decreased by 2.58% from Case 3 to Case 5.The decrease in the ion current from Case 3 to Case 5 is attributed to the decrease in atom density (approximately 5%, as seen from figure 10) and the shortened electron resident time(approximately 10%, represented by the increase in electron current).The decrease in the current efficiency explains the decrease in performance.The ion current increased by 10.2%from Case 0 to Case 3, while the thrust increased by 9.8%;this is attributed to the reduced ion loss at the wall as the high potential wall repels the ion bombardment.

    After further calculation, the divergence angle of the plume was obtained, as shown in figure 14.The divergence angle first exhibits a steep drop and then a slight rise.Because of the outward shifts of the magnetic field, the acceleration zone is located outside the channel, and the flow area of the channel exit determines the plume divergence.When the wall chamfer is large,as in Cases 4 and 5,it is difficult for the wall to restrain the beam from diverging, and the ions accelerate outward at the large flow area; thus, the plume is divergent.With a decrease in chamfer area,the radial divergence of ions is hindered by the high potential wall, and the plume divergence is slightly reduced, as in Case 3.However, when the chamfer area is further reduced, the magnetic field lines are interrupted by the wall, and the wall can no longer maintain a high potential and low electron temperature.The electric field is parallel to the local magnetic field, and the ions acquire radial acceleration near the wall and significantly diverge along the radial direction in Cases 0, 1, and 2.Thus,the large loss of energy and large plume divergence result in poor performance, as shown in figure 12.Accordingly, the thrust and efficiency are low in Cases 0, 1, and 2.

    4.Conclusion

    An experimental study was performed to quantitatively select the magnetic field line to match the chamfered channel wall tangentially in a magnetically shielded Hall thruster.For this purpose, a dimensionless characteristic parameterwas defined,i.e.the ratio of the field strength at the intersection of the selected magnetic field line and channel centerline to the maximum magnetic field strength along the centerline.

    A reduction inbrings the selected magnetic field line closer to the anode; consequently, the channel wall is better protected from ion sputtering.However, more deposition accumulates on the wall surface, degrading the wall insulation; thus, electron cross-field transport is enhanced through the wall,and the current efficiency deteriorates.In addition,asdecreases,the chamfered area of the channel wall increases.On one hand,the atom density in the channel decreases, and the ionization intensity declines;conversely,the wall energy loss caused by the ion bombardment is reduced.Considering these effects, it is inferred that asdecreases in the range of 20%–36%, the reduction in the wall energy loss is the primary factor controlling the change in thruster performance.Asdecreases further below 12%, the decline in both the current efficiency and ionization intensity governs the change in thruster performance.As a result, the thruster performance first increases and then decreases asdecreases.

    To quantitatively achieve high performance, the selected magnetic field line should have aof 12%–20%.Moreover,to ensure a satisfactory anti-sputtering effect at the wall, the selected magnetic field line should have aof less than 12%.

    Acknowledgments

    This work was funded by National Natural Science Foundation of China (Nos.52076054 and 51736003), Civil Aerospace Technology Pre-research Project (No.D03015), and Defense Industrial Technology Development Program (No.JCKY2019603B005).

    亚洲国产色片| 成年人免费黄色播放视频| 人妻少妇偷人精品九色| 观看美女的网站| 国产毛片在线视频| 97精品久久久久久久久久精品| 少妇的逼好多水| 免费观看a级毛片全部| 久久久久久久久久久免费av| 亚洲精品一区蜜桃| 婷婷色麻豆天堂久久| 国产亚洲一区二区精品| 最近2019中文字幕mv第一页| 欧美人与性动交α欧美软件 | 人妻少妇偷人精品九色| 人人澡人人妻人| 亚洲av中文av极速乱| 青春草视频在线免费观看| 欧美精品亚洲一区二区| 久久久久精品久久久久真实原创| 久久女婷五月综合色啪小说| 欧美 亚洲 国产 日韩一| 免费在线观看黄色视频的| 深夜精品福利| 麻豆精品久久久久久蜜桃| 丰满乱子伦码专区| 熟女电影av网| videosex国产| 欧美激情极品国产一区二区三区 | 黄色一级大片看看| 国产精品久久久久久久久免| 久久热在线av| av在线观看视频网站免费| 成年女人在线观看亚洲视频| 在现免费观看毛片| 中国国产av一级| 国产精品久久久久久久电影| 午夜精品国产一区二区电影| 精品第一国产精品| 爱豆传媒免费全集在线观看| 观看av在线不卡| 日韩一本色道免费dvd| 久久国产精品大桥未久av| 国产成人午夜福利电影在线观看| 久久精品国产a三级三级三级| 久久精品夜色国产| 美女xxoo啪啪120秒动态图| xxx大片免费视频| 视频区图区小说| 寂寞人妻少妇视频99o| 国产淫语在线视频| 超碰97精品在线观看| 一区二区三区乱码不卡18| 欧美国产精品一级二级三级| www.色视频.com| 18在线观看网站| 成人毛片60女人毛片免费| 伦精品一区二区三区| 少妇精品久久久久久久| 下体分泌物呈黄色| 午夜老司机福利剧场| 成人影院久久| 搡女人真爽免费视频火全软件| 国产成人精品无人区| 成年人免费黄色播放视频| 考比视频在线观看| 丝瓜视频免费看黄片| 侵犯人妻中文字幕一二三四区| 午夜福利视频精品| 精品久久国产蜜桃| av福利片在线| 免费高清在线观看视频在线观看| 亚洲人成网站在线观看播放| 飞空精品影院首页| 精品久久久精品久久久| 日韩三级伦理在线观看| 汤姆久久久久久久影院中文字幕| 欧美日本中文国产一区发布| 国产色爽女视频免费观看| 少妇人妻精品综合一区二区| 国产成人a∨麻豆精品| 欧美少妇被猛烈插入视频| 亚洲精品,欧美精品| 国产一区二区激情短视频 | 晚上一个人看的免费电影| 久久久久人妻精品一区果冻| 伊人亚洲综合成人网| 国产在线免费精品| 成年动漫av网址| 久久久久国产精品人妻一区二区| 国产成人精品久久久久久| 欧美+日韩+精品| 日韩一区二区三区影片| videossex国产| 国产高清国产精品国产三级| 哪个播放器可以免费观看大片| 捣出白浆h1v1| 夜夜爽夜夜爽视频| 国产乱人偷精品视频| 亚洲成人一二三区av| 国产老妇伦熟女老妇高清| 天美传媒精品一区二区| 一本—道久久a久久精品蜜桃钙片| 香蕉国产在线看| 伊人亚洲综合成人网| 色94色欧美一区二区| 80岁老熟妇乱子伦牲交| 乱人伦中国视频| 久久人人97超碰香蕉20202| a级毛色黄片| 男人爽女人下面视频在线观看| 国产免费一区二区三区四区乱码| 国产在视频线精品| 亚洲欧洲日产国产| 99热网站在线观看| 国产一区二区三区av在线| 老司机影院毛片| 国产又色又爽无遮挡免| 18禁观看日本| 欧美亚洲日本最大视频资源| 中文字幕制服av| 三上悠亚av全集在线观看| 中文字幕亚洲精品专区| 欧美日韩视频高清一区二区三区二| 国产欧美另类精品又又久久亚洲欧美| 日本av免费视频播放| 蜜桃在线观看..| 乱码一卡2卡4卡精品| 少妇人妻精品综合一区二区| 高清欧美精品videossex| 亚洲 欧美一区二区三区| 欧美性感艳星| 亚洲欧美成人综合另类久久久| 日本与韩国留学比较| 超色免费av| 美女大奶头黄色视频| 成人毛片60女人毛片免费| 一级片免费观看大全| 国产精品 国内视频| 满18在线观看网站| 午夜福利网站1000一区二区三区| 肉色欧美久久久久久久蜜桃| 国产一区二区激情短视频 | 久久久久国产精品人妻一区二区| 久久人人爽人人爽人人片va| 日韩制服骚丝袜av| 久久99精品国语久久久| 日本-黄色视频高清免费观看| a级片在线免费高清观看视频| 国产一区有黄有色的免费视频| 自线自在国产av| 欧美3d第一页| 亚洲精品第二区| 夜夜骑夜夜射夜夜干| 亚洲婷婷狠狠爱综合网| 啦啦啦中文免费视频观看日本| 视频中文字幕在线观看| 久久久久久久精品精品| av又黄又爽大尺度在线免费看| 视频在线观看一区二区三区| 欧美日韩成人在线一区二区| 蜜臀久久99精品久久宅男| 18在线观看网站| 亚洲激情五月婷婷啪啪| 日韩人妻精品一区2区三区| 22中文网久久字幕| 男女边摸边吃奶| 日韩av在线免费看完整版不卡| 只有这里有精品99| 欧美精品一区二区免费开放| 国产精品成人在线| 日产精品乱码卡一卡2卡三| 美女大奶头黄色视频| 久久99一区二区三区| 寂寞人妻少妇视频99o| 中文字幕人妻丝袜制服| av.在线天堂| 纵有疾风起免费观看全集完整版| 纵有疾风起免费观看全集完整版| 狠狠婷婷综合久久久久久88av| 亚洲人成77777在线视频| 亚洲第一av免费看| 成年av动漫网址| 18禁观看日本| 精品亚洲乱码少妇综合久久| 成人二区视频| 久久精品夜色国产| 国产精品一二三区在线看| 国产成人aa在线观看| 欧美 亚洲 国产 日韩一| 卡戴珊不雅视频在线播放| 亚洲精品国产色婷婷电影| 国国产精品蜜臀av免费| 多毛熟女@视频| 亚洲成国产人片在线观看| 亚洲精品色激情综合| 日本猛色少妇xxxxx猛交久久| av天堂久久9| 一二三四中文在线观看免费高清| 精品少妇内射三级| av不卡在线播放| 国产精品99久久99久久久不卡 | 99视频精品全部免费 在线| 亚洲久久久国产精品| 亚洲精品一二三| 九草在线视频观看| 国产精品.久久久| 伦理电影大哥的女人| 亚洲中文av在线| 国产一区二区激情短视频 | 免费播放大片免费观看视频在线观看| 亚洲美女视频黄频| 我的女老师完整版在线观看| 日韩中文字幕视频在线看片| 久久久久视频综合| 国产亚洲最大av| 欧美性感艳星| 9色porny在线观看| av在线播放精品| 大片免费播放器 马上看| 夜夜骑夜夜射夜夜干| 亚洲精品日本国产第一区| 欧美日韩av久久| 男人爽女人下面视频在线观看| 欧美 亚洲 国产 日韩一| 欧美人与性动交α欧美软件 | 1024视频免费在线观看| 成人毛片60女人毛片免费| 精品人妻一区二区三区麻豆| 日日摸夜夜添夜夜爱| 精品久久国产蜜桃| 看免费av毛片| 久久综合国产亚洲精品| 国精品久久久久久国模美| 久久韩国三级中文字幕| 人妻人人澡人人爽人人| 又黄又粗又硬又大视频| 国产亚洲精品久久久com| 丝袜脚勾引网站| 麻豆精品久久久久久蜜桃| 九九爱精品视频在线观看| 日韩视频在线欧美| 欧美激情国产日韩精品一区| 欧美最新免费一区二区三区| 亚洲一码二码三码区别大吗| 亚洲av中文av极速乱| 国产精品一国产av| 精品一区二区三区四区五区乱码 | 亚洲人成网站在线观看播放| 免费在线观看完整版高清| 在线天堂最新版资源| 这个男人来自地球电影免费观看 | 咕卡用的链子| 国产老妇伦熟女老妇高清| 热re99久久精品国产66热6| 精品99又大又爽又粗少妇毛片| 亚洲精品久久成人aⅴ小说| 亚洲人成77777在线视频| 男女免费视频国产| 婷婷色综合大香蕉| 九色成人免费人妻av| 久久青草综合色| 在线 av 中文字幕| 久久久久久人妻| 80岁老熟妇乱子伦牲交| 看免费成人av毛片| 亚洲精品自拍成人| 精品一区二区免费观看| 久久综合国产亚洲精品| 色婷婷久久久亚洲欧美| 午夜激情av网站| 熟女人妻精品中文字幕| 国产精品久久久久久精品电影小说| videos熟女内射| 狠狠婷婷综合久久久久久88av| 大香蕉久久成人网| 十分钟在线观看高清视频www| 日韩av免费高清视频| 美女福利国产在线| 人妻系列 视频| 亚洲国产看品久久| 国产免费现黄频在线看| 亚洲精品视频女| 亚洲国产精品专区欧美| 亚洲国产精品一区三区| 夜夜骑夜夜射夜夜干| 飞空精品影院首页| 精品亚洲成国产av| 尾随美女入室| 成年美女黄网站色视频大全免费| 免费在线观看黄色视频的| 国产精品一区二区在线不卡| 国产老妇伦熟女老妇高清| 九九爱精品视频在线观看| 国产成人午夜福利电影在线观看| 免费观看av网站的网址| 欧美bdsm另类| 卡戴珊不雅视频在线播放| 麻豆乱淫一区二区| h视频一区二区三区| 亚洲av中文av极速乱| 黑人欧美特级aaaaaa片| 9热在线视频观看99| 国产乱人偷精品视频| 亚洲国产日韩一区二区| 午夜日本视频在线| 一级片'在线观看视频| 咕卡用的链子| 欧美亚洲日本最大视频资源| 99国产精品免费福利视频| 久久99精品国语久久久| 大片免费播放器 马上看| 久久久久久久精品精品| 亚洲精品第二区| 人妻 亚洲 视频| 精品人妻偷拍中文字幕| 少妇的丰满在线观看| 热re99久久国产66热| www日本在线高清视频| tube8黄色片| 久久精品久久久久久久性| 亚洲国产精品国产精品| 亚洲av福利一区| 亚洲色图 男人天堂 中文字幕 | 国产欧美日韩一区二区三区在线| 亚洲国产av影院在线观看| av在线观看视频网站免费| 高清黄色对白视频在线免费看| av有码第一页| 免费看不卡的av| 精品人妻在线不人妻| 亚洲美女搞黄在线观看| 水蜜桃什么品种好| 少妇猛男粗大的猛烈进出视频| 女的被弄到高潮叫床怎么办| 超色免费av| 精品一区二区三区视频在线| 人体艺术视频欧美日本| 久久 成人 亚洲| kizo精华| 色哟哟·www| 国产日韩欧美亚洲二区| 国产国语露脸激情在线看| 熟女人妻精品中文字幕| 男女免费视频国产| 一本大道久久a久久精品| 少妇高潮的动态图| 日韩成人av中文字幕在线观看| 国产色爽女视频免费观看| 国产欧美日韩综合在线一区二区| 另类精品久久| a级毛片在线看网站| 亚洲精品自拍成人| 晚上一个人看的免费电影| 伦精品一区二区三区| 少妇人妻精品综合一区二区| 男女下面插进去视频免费观看 | 日本-黄色视频高清免费观看| 看十八女毛片水多多多| 看免费成人av毛片| 精品少妇内射三级| 一级毛片电影观看| 搡女人真爽免费视频火全软件| 丰满少妇做爰视频| 久久女婷五月综合色啪小说| 亚洲色图综合在线观看| 99国产精品免费福利视频| 国产片内射在线| 亚洲国产看品久久| 五月玫瑰六月丁香| videos熟女内射| 久久久久精品久久久久真实原创| 国产精品久久久久久久久免| 国产 一区精品| 韩国av在线不卡| 亚洲av欧美aⅴ国产| 高清欧美精品videossex| 亚洲av电影在线进入| 亚洲精品av麻豆狂野| 欧美激情国产日韩精品一区| 青青草视频在线视频观看| 午夜免费男女啪啪视频观看| 国产视频首页在线观看| 午夜免费男女啪啪视频观看| 国产视频首页在线观看| 亚洲欧美成人精品一区二区| 伦理电影免费视频| 热99久久久久精品小说推荐| 激情五月婷婷亚洲| 少妇被粗大的猛进出69影院 | 丝瓜视频免费看黄片| 欧美人与性动交α欧美精品济南到 | 大香蕉久久成人网| 韩国av在线不卡| 国产黄色免费在线视频| 少妇人妻精品综合一区二区| 亚洲精品久久久久久婷婷小说| 亚洲人与动物交配视频| 汤姆久久久久久久影院中文字幕| 韩国av在线不卡| 国产极品天堂在线| 建设人人有责人人尽责人人享有的| 日韩av不卡免费在线播放| 人妻少妇偷人精品九色| 亚洲av中文av极速乱| 午夜激情av网站| 精品福利永久在线观看| 9热在线视频观看99| 男女下面插进去视频免费观看 | 国产精品无大码| 国产乱来视频区| 日韩一本色道免费dvd| 日韩 亚洲 欧美在线| 精品人妻偷拍中文字幕| 日韩欧美精品免费久久| 欧美+日韩+精品| 欧美性感艳星| 亚洲久久久国产精品| 欧美日韩视频高清一区二区三区二| 宅男免费午夜| 久久99热这里只频精品6学生| 久久鲁丝午夜福利片| 亚洲色图综合在线观看| 久久女婷五月综合色啪小说| videossex国产| 一级毛片黄色毛片免费观看视频| 国产精品嫩草影院av在线观看| 满18在线观看网站| 欧美精品亚洲一区二区| 2022亚洲国产成人精品| 热re99久久精品国产66热6| 一区在线观看完整版| 中文字幕亚洲精品专区| 久久这里只有精品19| 乱码一卡2卡4卡精品| xxx大片免费视频| 精品一品国产午夜福利视频| 9色porny在线观看| 性高湖久久久久久久久免费观看| 一本久久精品| 亚洲av欧美aⅴ国产| 国产精品国产三级国产专区5o| 亚洲色图 男人天堂 中文字幕 | 欧美精品一区二区免费开放| 午夜精品国产一区二区电影| 免费观看a级毛片全部| 狠狠精品人妻久久久久久综合| 亚洲国产av新网站| 国产精品久久久av美女十八| 黑人猛操日本美女一级片| 久久av网站| 亚洲,一卡二卡三卡| 久久女婷五月综合色啪小说| 男人爽女人下面视频在线观看| 22中文网久久字幕| 另类精品久久| 亚洲国产av新网站| 成年女人在线观看亚洲视频| 观看美女的网站| 精品酒店卫生间| 国产免费一区二区三区四区乱码| 永久网站在线| 一边摸一边做爽爽视频免费| 精品国产乱码久久久久久小说| 黑人猛操日本美女一级片| 欧美另类一区| 人妻一区二区av| 中文字幕av电影在线播放| 大陆偷拍与自拍| 国产精品嫩草影院av在线观看| 婷婷色麻豆天堂久久| 黄色 视频免费看| 中文天堂在线官网| 观看美女的网站| 777米奇影视久久| 大片免费播放器 马上看| 伊人久久国产一区二区| 日韩熟女老妇一区二区性免费视频| 婷婷色综合大香蕉| 日韩成人av中文字幕在线观看| 欧美丝袜亚洲另类| 熟女av电影| 国产 一区精品| 免费观看无遮挡的男女| 亚洲精品国产av蜜桃| 亚洲欧美精品自产自拍| 青春草亚洲视频在线观看| 亚洲国产成人一精品久久久| 你懂的网址亚洲精品在线观看| 久久影院123| 亚洲图色成人| 一区二区三区精品91| 国产精品久久久久成人av| 在线观看免费视频网站a站| 少妇人妻 视频| 少妇熟女欧美另类| 夫妻性生交免费视频一级片| 亚洲国产精品一区三区| 午夜影院在线不卡| 久久精品国产自在天天线| 黄色视频在线播放观看不卡| 波野结衣二区三区在线| 亚洲四区av| 在线观看免费日韩欧美大片| 亚洲精品一区蜜桃| 一级a做视频免费观看| 亚洲,欧美,日韩| 成人无遮挡网站| 永久免费av网站大全| 妹子高潮喷水视频| av免费在线看不卡| 少妇被粗大的猛进出69影院 | 18禁动态无遮挡网站| 爱豆传媒免费全集在线观看| 日韩成人伦理影院| 午夜老司机福利剧场| 日韩一本色道免费dvd| 精品99又大又爽又粗少妇毛片| 又粗又硬又长又爽又黄的视频| 亚洲久久久国产精品| 国产亚洲午夜精品一区二区久久| 欧美+日韩+精品| 最近手机中文字幕大全| 大香蕉久久成人网| 一个人免费看片子| 青春草亚洲视频在线观看| 久久久久久久久久人人人人人人| 最近2019中文字幕mv第一页| 99热国产这里只有精品6| 国产乱人偷精品视频| 国国产精品蜜臀av免费| 黄色一级大片看看| 久久久久网色| 日本色播在线视频| 夜夜爽夜夜爽视频| 国产乱来视频区| 精品一区二区三区四区五区乱码 | 99热网站在线观看| 一个人免费看片子| 视频中文字幕在线观看| 一级毛片黄色毛片免费观看视频| 国产精品一二三区在线看| 精品一区二区免费观看| 国产又色又爽无遮挡免| 18禁在线无遮挡免费观看视频| 母亲3免费完整高清在线观看 | 日韩电影二区| 啦啦啦视频在线资源免费观看| 国产成人精品无人区| 国产精品久久久久久久久免| 久久久国产一区二区| 亚洲中文av在线| 美女国产高潮福利片在线看| 99香蕉大伊视频| 夜夜爽夜夜爽视频| 久久 成人 亚洲| a级毛片在线看网站| 丝瓜视频免费看黄片| 国产视频首页在线观看| 国产av码专区亚洲av| 亚洲精品成人av观看孕妇| 日日撸夜夜添| 亚洲久久久国产精品| www日本在线高清视频| 国产熟女午夜一区二区三区| 夜夜爽夜夜爽视频| 在线观看免费高清a一片| 1024视频免费在线观看| 久久精品国产鲁丝片午夜精品| 国产精品一区二区在线不卡| 国产av码专区亚洲av| 免费看av在线观看网站| 久久久国产一区二区| videosex国产| 久久久久国产精品人妻一区二区| 亚洲国产av影院在线观看| 国产精品国产av在线观看| 美国免费a级毛片| 亚洲色图 男人天堂 中文字幕 | 我的女老师完整版在线观看| 亚洲,欧美精品.| 欧美国产精品va在线观看不卡| av免费观看日本| 日韩电影二区| 精品卡一卡二卡四卡免费| 在线观看免费日韩欧美大片| 巨乳人妻的诱惑在线观看| 91在线精品国自产拍蜜月| 精品人妻偷拍中文字幕| 精品亚洲成a人片在线观看| 久久综合国产亚洲精品| 边亲边吃奶的免费视频| 久久精品国产a三级三级三级| 成人漫画全彩无遮挡| 免费女性裸体啪啪无遮挡网站| 国产无遮挡羞羞视频在线观看| 九色成人免费人妻av| 国产男女超爽视频在线观看| 国产亚洲一区二区精品| 大码成人一级视频| 五月开心婷婷网| 久久精品夜色国产| 亚洲精品乱久久久久久| 少妇高潮的动态图| 99热这里只有是精品在线观看| 又粗又硬又长又爽又黄的视频| 久久精品人人爽人人爽视色| 99热全是精品| 中文精品一卡2卡3卡4更新| 女人久久www免费人成看片| 亚洲精品成人av观看孕妇| 在线观看三级黄色| 亚洲人成77777在线视频| 欧美丝袜亚洲另类| 午夜激情久久久久久久|