• <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).

    一本—道久久a久久精品蜜桃钙片| 亚洲一区二区三区欧美精品| 黄色视频不卡| 叶爱在线成人免费视频播放| 午夜精品久久久久久毛片777| 成人影院久久| 侵犯人妻中文字幕一二三四区| 亚洲少妇的诱惑av| 老汉色av国产亚洲站长工具| 少妇裸体淫交视频免费看高清 | 成人永久免费在线观看视频 | 亚洲国产av影院在线观看| 久久久久久久久免费视频了| 亚洲专区国产一区二区| 伊人久久大香线蕉亚洲五| 好男人电影高清在线观看| 人妻 亚洲 视频| 国产亚洲av高清不卡| www.999成人在线观看| 色在线成人网| 丁香欧美五月| 午夜福利欧美成人| 搡老岳熟女国产| 丝袜喷水一区| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲免费av在线视频| 久久国产精品男人的天堂亚洲| 精品久久久久久电影网| 自线自在国产av| 免费少妇av软件| 久久中文看片网| 一本一本久久a久久精品综合妖精| 久久精品熟女亚洲av麻豆精品| 国产一卡二卡三卡精品| 国产精品.久久久| 免费女性裸体啪啪无遮挡网站| 亚洲一区二区三区欧美精品| 操美女的视频在线观看| 精品第一国产精品| 欧美性长视频在线观看| 欧美日韩黄片免| 亚洲 国产 在线| 日韩视频在线欧美| 老熟妇仑乱视频hdxx| 啦啦啦视频在线资源免费观看| av国产精品久久久久影院| 色婷婷久久久亚洲欧美| 999精品在线视频| 蜜桃国产av成人99| 嫩草影视91久久| a级片在线免费高清观看视频| 成人亚洲精品一区在线观看| 国产精品久久久人人做人人爽| 亚洲伊人色综图| 成人av一区二区三区在线看| 伊人久久大香线蕉亚洲五| 午夜福利乱码中文字幕| 亚洲熟妇熟女久久| 高清黄色对白视频在线免费看| 国产色视频综合| 色婷婷av一区二区三区视频| 久久热在线av| 一本久久精品| 91精品国产国语对白视频| 精品高清国产在线一区| 飞空精品影院首页| 欧美午夜高清在线| 啪啪无遮挡十八禁网站| 老司机午夜福利在线观看视频 | 99久久99久久久精品蜜桃| 麻豆成人av在线观看| 欧美av亚洲av综合av国产av| videosex国产| 一级,二级,三级黄色视频| 国产成人av教育| 午夜久久久在线观看| 中文字幕高清在线视频| 亚洲三区欧美一区| 亚洲美女黄片视频| 免费少妇av软件| 国产一区二区三区在线臀色熟女 | 免费少妇av软件| 国产精品.久久久| 亚洲国产成人一精品久久久| av在线播放免费不卡| 亚洲性夜色夜夜综合| 亚洲天堂av无毛| 午夜久久久在线观看| 丝袜美足系列| 一区二区三区激情视频| 嫩草影视91久久| 亚洲久久久国产精品| 中文字幕另类日韩欧美亚洲嫩草| 肉色欧美久久久久久久蜜桃| 午夜精品久久久久久毛片777| 在线av久久热| 黑人欧美特级aaaaaa片| 国产在视频线精品| 在线观看免费视频日本深夜| 欧美av亚洲av综合av国产av| 美女高潮喷水抽搐中文字幕| 日韩欧美免费精品| 欧美成狂野欧美在线观看| 久久精品国产99精品国产亚洲性色 | 成人手机av| 久久久久国产一级毛片高清牌| 丝袜美足系列| 在线观看一区二区三区激情| 十八禁人妻一区二区| 国产日韩一区二区三区精品不卡| 久久人人97超碰香蕉20202| 国产精品自产拍在线观看55亚洲 | 成人永久免费在线观看视频 | 精品一区二区三区av网在线观看 | 欧美亚洲 丝袜 人妻 在线| 久久久久久久久免费视频了| 丁香六月天网| 超色免费av| 久久青草综合色| 91字幕亚洲| 俄罗斯特黄特色一大片| 久久ye,这里只有精品| 999久久久精品免费观看国产| 美女高潮喷水抽搐中文字幕| 首页视频小说图片口味搜索| 日韩一区二区三区影片| 午夜福利,免费看| 91精品三级在线观看| 女同久久另类99精品国产91| 国产精品1区2区在线观看. | 日日摸夜夜添夜夜添小说| 黄色成人免费大全| 乱人伦中国视频| 精品人妻熟女毛片av久久网站| 一级,二级,三级黄色视频| √禁漫天堂资源中文www| 久久 成人 亚洲| 老熟妇乱子伦视频在线观看| 国产精品秋霞免费鲁丝片| 后天国语完整版免费观看| 首页视频小说图片口味搜索| 国产高清国产精品国产三级| 免费在线观看视频国产中文字幕亚洲| 婷婷成人精品国产| 久久精品国产亚洲av香蕉五月 | 国产亚洲精品第一综合不卡| 欧美黑人精品巨大| 黄色视频,在线免费观看| 成年人黄色毛片网站| 2018国产大陆天天弄谢| 一本综合久久免费| 悠悠久久av| 日韩欧美三级三区| 国内毛片毛片毛片毛片毛片| 国产亚洲午夜精品一区二区久久| 久久久久久免费高清国产稀缺| av有码第一页| 午夜激情久久久久久久| 18禁观看日本| 满18在线观看网站| 啦啦啦中文免费视频观看日本| 国产在线一区二区三区精| 亚洲av日韩在线播放| 91大片在线观看| 亚洲国产成人一精品久久久| cao死你这个sao货| 精品高清国产在线一区| 亚洲视频免费观看视频| 精品国产乱码久久久久久小说| 国产在视频线精品| 他把我摸到了高潮在线观看 | videos熟女内射| 精品福利永久在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲第一av免费看| 一区在线观看完整版| 黑人操中国人逼视频| 精品国产一区二区三区四区第35| 美女主播在线视频| 中文字幕精品免费在线观看视频| 国产不卡av网站在线观看| bbb黄色大片| 一本综合久久免费| 亚洲熟妇熟女久久| 法律面前人人平等表现在哪些方面| 亚洲精品在线美女| 亚洲一码二码三码区别大吗| 在线观看免费高清a一片| 国产日韩欧美在线精品| 久久午夜亚洲精品久久| 国产精品国产高清国产av | 国产精品成人在线| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品熟女久久久久浪| 久久天堂一区二区三区四区| 18禁黄网站禁片午夜丰满| 五月天丁香电影| 亚洲久久久国产精品| 成人国语在线视频| 国产亚洲av高清不卡| 一个人免费在线观看的高清视频| 成人18禁高潮啪啪吃奶动态图| 99国产极品粉嫩在线观看| 久久久久久久久久久久大奶| 国产99久久九九免费精品| 丁香欧美五月| 操出白浆在线播放| 中文亚洲av片在线观看爽 | 女人精品久久久久毛片| 久久久久久免费高清国产稀缺| 久久久国产成人免费| 色尼玛亚洲综合影院| 欧美国产精品va在线观看不卡| 国产亚洲精品第一综合不卡| 久久中文字幕一级| e午夜精品久久久久久久| 国产免费视频播放在线视频| 好男人电影高清在线观看| 国产一区二区三区在线臀色熟女 | 亚洲精品在线观看二区| 国产av精品麻豆| 久久久国产一区二区| 成年版毛片免费区| 中文字幕人妻丝袜一区二区| 最黄视频免费看| 成人永久免费在线观看视频 | 正在播放国产对白刺激| 亚洲精品成人av观看孕妇| 日本精品一区二区三区蜜桃| 别揉我奶头~嗯~啊~动态视频| 在线观看舔阴道视频| 国产又爽黄色视频| 如日韩欧美国产精品一区二区三区| 国产高清国产精品国产三级| 麻豆乱淫一区二区| 亚洲精品美女久久av网站| 午夜老司机福利片| 久久热在线av| 精品亚洲成a人片在线观看| 亚洲欧美一区二区三区黑人| 搡老乐熟女国产| 精品亚洲乱码少妇综合久久| 亚洲黑人精品在线| 国产在线精品亚洲第一网站| 欧美日韩亚洲高清精品| 一个人免费在线观看的高清视频| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久精品电影小说| 国产一卡二卡三卡精品| 免费观看a级毛片全部| 精品国产乱子伦一区二区三区| 亚洲国产成人一精品久久久| 国产成人精品在线电影| 人人妻人人澡人人看| 丰满少妇做爰视频| 欧美 日韩 精品 国产| 欧美激情极品国产一区二区三区| 亚洲国产欧美一区二区综合| 亚洲第一青青草原| av网站免费在线观看视频| 日本一区二区免费在线视频| 日日夜夜操网爽| 欧美+亚洲+日韩+国产| 国产男女超爽视频在线观看| 国产精品亚洲av一区麻豆| 欧美人与性动交α欧美软件| 久久狼人影院| 中国美女看黄片| 青青草视频在线视频观看| 夜夜夜夜夜久久久久| 女人精品久久久久毛片| 男女高潮啪啪啪动态图| 老司机亚洲免费影院| 欧美精品高潮呻吟av久久| av有码第一页| 日本av免费视频播放| 男女高潮啪啪啪动态图| 国产精品亚洲一级av第二区| 久久av网站| 最近最新免费中文字幕在线| 国产在线免费精品| 黄色片一级片一级黄色片| 成年人午夜在线观看视频| 亚洲,欧美精品.| 国产区一区二久久| 女同久久另类99精品国产91| 91成年电影在线观看| 久久久国产一区二区| 久久久精品区二区三区| 90打野战视频偷拍视频| 桃花免费在线播放| 757午夜福利合集在线观看| av天堂在线播放| 亚洲av国产av综合av卡| 99久久99久久久精品蜜桃| 他把我摸到了高潮在线观看 | 50天的宝宝边吃奶边哭怎么回事| 男人舔女人的私密视频| 亚洲av第一区精品v没综合| 国产日韩一区二区三区精品不卡| 自线自在国产av| 男女下面插进去视频免费观看| 淫妇啪啪啪对白视频| 日韩欧美国产一区二区入口| 中文字幕人妻丝袜制服| 国产激情久久老熟女| 天堂俺去俺来也www色官网| 高清视频免费观看一区二区| 69精品国产乱码久久久| 无限看片的www在线观看| 精品人妻在线不人妻| 中文欧美无线码| 美国免费a级毛片| 宅男免费午夜| 黄色片一级片一级黄色片| 激情在线观看视频在线高清 | 天天添夜夜摸| 精品午夜福利视频在线观看一区 | 少妇被粗大的猛进出69影院| 精品国产乱码久久久久久小说| 亚洲国产欧美在线一区| 婷婷丁香在线五月| 五月天丁香电影| a级片在线免费高清观看视频| 国产色视频综合| 亚洲午夜理论影院| 国产精品久久久人人做人人爽| 99久久国产精品久久久| 高清欧美精品videossex| 精品一品国产午夜福利视频| 亚洲av第一区精品v没综合| 日韩中文字幕视频在线看片| 精品少妇一区二区三区视频日本电影| 欧美黑人精品巨大| 亚洲熟妇熟女久久| 精品少妇一区二区三区视频日本电影| avwww免费| 黑人操中国人逼视频| 啦啦啦免费观看视频1| 欧美精品人与动牲交sv欧美| 色综合婷婷激情| 91麻豆精品激情在线观看国产 | netflix在线观看网站| h视频一区二区三区| 久久精品成人免费网站| 久久热在线av| 又紧又爽又黄一区二区| 中文字幕精品免费在线观看视频| 国产无遮挡羞羞视频在线观看| 国产在线视频一区二区| 在线av久久热| 日韩中文字幕视频在线看片| 人人澡人人妻人| 免费少妇av软件| 性少妇av在线| 国产亚洲一区二区精品| 母亲3免费完整高清在线观看| 搡老熟女国产l中国老女人| 国产一卡二卡三卡精品| 99国产综合亚洲精品| 色视频在线一区二区三区| 免费黄频网站在线观看国产| 亚洲黑人精品在线| 日日爽夜夜爽网站| 97在线人人人人妻| 在线观看免费日韩欧美大片| 亚洲专区字幕在线| 99re6热这里在线精品视频| 精品欧美一区二区三区在线| 12—13女人毛片做爰片一| cao死你这个sao货| 美女福利国产在线| 久热这里只有精品99| 久久久国产成人免费| 一区二区三区乱码不卡18| 欧美黄色片欧美黄色片| e午夜精品久久久久久久| 国产精品亚洲一级av第二区| 99精国产麻豆久久婷婷| 久久热在线av| 久久国产精品男人的天堂亚洲| 精品福利永久在线观看| 午夜福利视频精品| 亚洲av美国av| 狂野欧美激情性xxxx| 日韩制服丝袜自拍偷拍| 亚洲一区二区三区欧美精品| 香蕉久久夜色| 老汉色av国产亚洲站长工具| 99久久国产精品久久久| 精品国产亚洲在线| 国产成人精品无人区| 免费在线观看日本一区| 亚洲国产看品久久| 一区在线观看完整版| 精品欧美一区二区三区在线| 国产av又大| 热re99久久国产66热| 美女午夜性视频免费| 亚洲熟女精品中文字幕| 五月天丁香电影| 欧美日韩中文字幕国产精品一区二区三区 | 丁香欧美五月| 久久这里只有精品19| 建设人人有责人人尽责人人享有的| 女性生殖器流出的白浆| 啦啦啦在线免费观看视频4| 五月开心婷婷网| 大片免费播放器 马上看| 91成年电影在线观看| 美女福利国产在线| 欧美激情高清一区二区三区| 手机成人av网站| 一级片免费观看大全| 在线 av 中文字幕| 啪啪无遮挡十八禁网站| 亚洲欧美日韩另类电影网站| tocl精华| 亚洲专区中文字幕在线| 夜夜夜夜夜久久久久| 日韩成人在线观看一区二区三区| 久久久国产一区二区| 国产高清激情床上av| √禁漫天堂资源中文www| 国产免费av片在线观看野外av| 制服人妻中文乱码| 宅男免费午夜| 午夜免费鲁丝| 久久国产精品大桥未久av| 高清欧美精品videossex| 亚洲黑人精品在线| 另类精品久久| 日本撒尿小便嘘嘘汇集6| 亚洲国产毛片av蜜桃av| 免费av中文字幕在线| 人妻一区二区av| 国产精品电影一区二区三区 | 最近最新中文字幕大全免费视频| 国产成人啪精品午夜网站| 国产激情久久老熟女| 女人高潮潮喷娇喘18禁视频| 男人操女人黄网站| 亚洲精品国产区一区二| 国产精品偷伦视频观看了| 欧美成人午夜精品| 捣出白浆h1v1| 男女高潮啪啪啪动态图| 丝瓜视频免费看黄片| 久久久久国内视频| 夜夜骑夜夜射夜夜干| 亚洲av日韩精品久久久久久密| 欧美不卡视频在线免费观看| 神马国产精品三级电影在线观看| 亚洲精品粉嫩美女一区| 综合色av麻豆| 国产精品永久免费网站| 精品国产亚洲在线| 久久中文字幕一级| 久久人妻av系列| 97碰自拍视频| 亚洲精品中文字幕一二三四区| 久久久久精品国产欧美久久久| 人妻夜夜爽99麻豆av| 搡老妇女老女人老熟妇| 91老司机精品| 国内揄拍国产精品人妻在线| av女优亚洲男人天堂 | 真人做人爱边吃奶动态| 国产又色又爽无遮挡免费看| 成人av在线播放网站| 免费在线观看亚洲国产| 一进一出抽搐动态| 欧美日韩精品网址| 一二三四社区在线视频社区8| 熟女人妻精品中文字幕| 欧美不卡视频在线免费观看| 十八禁网站免费在线| 国产激情久久老熟女| av天堂在线播放| 亚洲欧美日韩高清专用| 婷婷亚洲欧美| 亚洲性夜色夜夜综合| 观看免费一级毛片| 一级a爱片免费观看的视频| 午夜激情欧美在线| 国产视频一区二区在线看| 久久中文字幕人妻熟女| 国产成人av激情在线播放| 国产精品av视频在线免费观看| 老司机午夜十八禁免费视频| 三级国产精品欧美在线观看 | 亚洲午夜理论影院| 亚洲18禁久久av| www.精华液| 少妇人妻一区二区三区视频| 国产成人一区二区三区免费视频网站| 国产极品精品免费视频能看的| 免费搜索国产男女视频| 国产精品99久久久久久久久| 久久久久久大精品| 国产高清videossex| 人人妻,人人澡人人爽秒播| 操出白浆在线播放| 巨乳人妻的诱惑在线观看| 国产综合懂色| 视频区欧美日本亚洲| 欧美日韩瑟瑟在线播放| 后天国语完整版免费观看| 精品人妻1区二区| 亚洲男人的天堂狠狠| 亚洲熟妇中文字幕五十中出| 国产亚洲精品一区二区www| 欧美黑人巨大hd| 黄色 视频免费看| 性欧美人与动物交配| 一进一出好大好爽视频| 亚洲av片天天在线观看| 久久精品91无色码中文字幕| 国产三级中文精品| 成人18禁在线播放| 午夜福利免费观看在线| 免费无遮挡裸体视频| 亚洲欧美日韩高清专用| av国产免费在线观看| 国产精品一区二区免费欧美| 国产精品精品国产色婷婷| 欧美日韩综合久久久久久 | 亚洲成人免费电影在线观看| 中文字幕精品亚洲无线码一区| 国产真人三级小视频在线观看| 亚洲av成人精品一区久久| 国产精品女同一区二区软件 | av片东京热男人的天堂| 一进一出好大好爽视频| 国产亚洲精品av在线| 国产精品久久视频播放| 好男人电影高清在线观看| 成熟少妇高潮喷水视频| 国产一区二区在线av高清观看| 岛国在线观看网站| 精品熟女少妇八av免费久了| 亚洲成人中文字幕在线播放| 亚洲欧美日韩高清在线视频| 国产高清有码在线观看视频| 不卡一级毛片| 99久久综合精品五月天人人| 国产精品免费一区二区三区在线| 韩国av一区二区三区四区| 一本久久中文字幕| 色在线成人网| 午夜精品一区二区三区免费看| 精品免费久久久久久久清纯| 每晚都被弄得嗷嗷叫到高潮| 成熟少妇高潮喷水视频| 国产一区二区激情短视频| 曰老女人黄片| 18禁黄网站禁片午夜丰满| 久久这里只有精品中国| 一级毛片高清免费大全| 69av精品久久久久久| 国产精品综合久久久久久久免费| 欧美日韩亚洲国产一区二区在线观看| 久久香蕉国产精品| av天堂中文字幕网| 99久久精品国产亚洲精品| 久久婷婷人人爽人人干人人爱| 小蜜桃在线观看免费完整版高清| 免费看十八禁软件| 色在线成人网| 国产探花在线观看一区二区| 啪啪无遮挡十八禁网站| 欧美+亚洲+日韩+国产| 国产高清videossex| 国产精品女同一区二区软件 | 国产精品亚洲av一区麻豆| 亚洲av美国av| 亚洲国产精品合色在线| 黄色片一级片一级黄色片| 九九久久精品国产亚洲av麻豆 | 老汉色av国产亚洲站长工具| 色老头精品视频在线观看| 亚洲人成电影免费在线| 色综合亚洲欧美另类图片| 午夜日韩欧美国产| 一本久久中文字幕| 五月伊人婷婷丁香| 日韩欧美精品v在线| 久久久精品大字幕| 免费人成视频x8x8入口观看| 18禁观看日本| 精品久久久久久久久久久久久| 天天添夜夜摸| 国产一区二区三区视频了| 成熟少妇高潮喷水视频| 亚洲av电影在线进入| 午夜两性在线视频| 色综合欧美亚洲国产小说| 国产一区二区激情短视频| 在线观看免费视频日本深夜| 又爽又黄无遮挡网站| 亚洲精品国产精品久久久不卡| 又紧又爽又黄一区二区| 亚洲精品456在线播放app | 亚洲va日本ⅴa欧美va伊人久久| 久久性视频一级片| 成人国产一区最新在线观看| 99热6这里只有精品| 村上凉子中文字幕在线| 色av中文字幕| 亚洲激情在线av| 国内精品久久久久精免费| 欧美性猛交黑人性爽| 国产伦在线观看视频一区| 久久九九热精品免费| 亚洲国产看品久久|