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

    Influence of anode temperature on ignition performance of the IRIT4-2D iodine-fueled radio frequency ion thruster

    2022-02-15 11:08:44XinghuaZHANG張興華ZhenhuaZHANG張振華ShaoxiaJIA賈少霞TingJIN金婷JinghuaYANG楊景華LongLI李龍FangfangLIU劉芳芳YongCAI蔡勇andJianCAI蔡建
    Plasma Science and Technology 2022年1期
    關(guān)鍵詞:李龍興華振華

    Xinghua ZHANG(張興華),Zhenhua ZHANG(張振華),Shaoxia JIA(賈少霞), Ting JIN (金婷), Jinghua YANG (楊景華), Long LI (李龍),Fangfang LIU (劉芳芳), Yong CAI (蔡勇) and Jian CAI (蔡建)

    1 Institute of Microelectronics,Chinese Academy of Sciences,Beijing 100029,People’s Republic of China

    2 Suzhou Nano-Space Dynamics Technology Co.Ltd, Suzhou 215000, People’s Republic of China

    Abstract This paper reports the ignition performance of the iodine-fueled radio frequency (RF) ion thruster (IRIT) at different anode temperatures (Ta0).The experimental results show that the anode temperature plays important role on the ignition process of the IRIT.There were two characteristic temperatures related to the anode: the minimum ignition temperature (Ti0) and the stable ignition temperature (Tis), which were much lower than the pipeline temperature and the storage tank temperature.At Ta0 < Ti 0 ,it failed to discharge.When Ti0 ≤ Ta0

    Keywords:radio frequency ion thruster,iodine propellant,ignition,deposition,sublimation,clog

    Nomenclature

    Iceelectron current emitted by the neutralizer/cathode

    Iscscreen grid current

    ˙maiodine flow rate into the discharge chamber

    ˙mdepiodine flow rate deposited in the anode

    ˙msubiodine flow rate sublimated from the anode

    ˙mtiodine flow rate out from the storage tank

    pvcpressure of the vacuum chamber

    RFnet input radio frequency power

    RFinincident radio frequency power

    RFrereflected radio frequency power

    Taanode/thruster temperature

    Ta0Taat ignition

    Ti0minimum ignition temperature

    Tisstable ignition temperature

    Ttanktemperature of the storage tank

    Tlinetemperature of the pipeline

    Uacacceleration grid voltage

    Uscscreen grid voltage

    VSvalve switching state

    1.Introduction

    The electric propulsion fueled by solid iodine propellant has become a research highlight in the past decades.Iodine is a feasible alternative to xenon to be used as a propellant for electric propulsion systems.There are obvious advantages using iodine instead of xenon:iodine is close to xenon in the propulsion performance, its density is 3 times of xenon(4.9 g cm?3@ room temperature versus 1.6 g cm?3@ room temperature and 2000 psi) and the price is 1/10, and the low pressure propellant feed system can be used, which is expected to make the propulsion system lighter, smaller and lower cost.Since the concept of iodine-fueled electric propulsion was put forward by Dressler et al [1] and Tverdokhlebov et al [2], it has been widely studied.Szabo et al[3–5], Kamhawi et al [6] and Smith et al [7] performed tests on the BHT series iodine Hall thrusters including BHT-200-I,BHT-600-I, BHT-1000 and BHT-8000.Liu et al [8] studied the ignition of the iodine cusped field thruster.Tsay et al[9, 10] studied the propulsion performance of the iodinefueled mini radio frequency (RF) ion thruster (IRIT) BIT3.Experiments were conducted on an iodine-fueled radio-frequency ion-thruster RIT-10 by Holste et al [11].Martínez et al [12] presented the development and testing of the NPT30-I2 ion gridded thruster.Yang et al [13] studied the performance of an iodine radio ion thruster IRIT4.The response of the magnetically enhanced plasma thruster has been studied by Manente et al [14, 15].In the above studies,the two propellants of xenon and iodine have shown comparable performances.

    One of the possible issues that can arise in an iodine propulsion system is that the vapor can redeposit to form a clog/blockage if a cold spot develops in the propellant lines,which makes the ignition characteristics of the iodine electric thruster at a lower temperature quite different from that of xenon.In general, to keep iodine from redepositing somewhere else within the system, the propellant lines should be kept significantly warmer than the solid iodine that is subliming to produce the iodine vapor[16].But for the thruster/anode, the situation is not exactly the same as the propellant lines of the propellant feed systems.On one hand, the anode acts as a cold sink on the lines, and that is a worst thermal case during the thruster ignition period as described by Samples et al [17].On the other hand, the thruster/anode is the propellant exit to vacuum and the low-pressure section,where it is relatively difficult to completely clog and it is possible for the thruster to ignite at a lower temperature.In the experiments of Paganucci et al [18], the HT100 Hall thruster was mounted on a plate with installed heaters to ensure that a temperature higher than 80 °C at the anode to avoid iodine condensation.But Sazbo et al [4] pointed out that the anode gas distributor does not have to be heated during ordinary thruster operation.

    To the best of our knowledge,no systematic study about the ignition characteristics of thrusters at different temperatures has been reported so far.This work will focus on the experimental study of the ignition characteristics of the IRIT4-2D IRIT [14] at different initial temperatures,especially at the temperatures lower than the propellant lines and the storage tank, and try to disclose the effect of iodine deposition-sublimation properties on the propulsion performance.The start-up temperature of iodine fueled thruster will be determined through the research, which is helpful to the thrust control strategy under different temperatures.

    2.Experimental procedure

    2.1.System composition

    The IRIT4-2D RF ion thruster (figure 1(a)) is a mission enabling, iodine-fueled ion propulsion system scheduled for launch in 2021, and it is a fully integrated propulsion system[19].The nominal operating point of the IRIT4-2D is at 150 W input power and 40 000 Ns total impulse, which produces 2 mN thrust and 2000 s specific impulse.The principle of the IRIT is shown in figure 1(b).Solid iodine propellant stored in the low pressure storage tank is heated to sublimate to a gaseous state, which is transferred into the thruster through the pipeline and then ionized to generate plasma by the coil-induced electromagnetic field, finally the ions are accelerated and extracted by the grid electrode to generate the thrust.Ions are neutralized by electrons emitted by neutralizers.

    The IRIT4-2D has a modular design, and include two sets of ion thrusters (anode), two sets of neutralizers (cathode), one iodine feed module as well as passive thermal management,and the power processing unit(PPU)integrated with the controller unit.The thruster has a diameter of 3.5 cm,which is an improved version of the 4 cm diameter thruster[13].By regulating the temperature of the storage tank,a flow rate range of 30–500 μg s?1can be achieved with the error less than±5%.The PPU module provides input power for the thruster, neutralizer and the iodine feed module, and realizes parameter monitoring and closed-loop control of thrust; the radio frequency is 2 MHz, and the DC voltage can reach 1800 V.

    The IRIT4-2D system used in the experiments in this work has passed the qualification level environmental test.Only one thruster, one neutralizer, a set of pipelines and one matcher were used in the experiment.For the laboratory purpose, the thruster, the RF power matcher, and the iodine feed module were placed in the vacuum chamber, while the PPU was located on the atmospheric side for ease of access to it (figure 2(a)).RIGOL DP832A and ITECH IT6942A DC power supply were used to supplied power the PPU with the input voltage of 28 ± 0.01 V.A computer with the special software was used for experiment control and data acquisition.

    Figure 1.The IRIT4-2D engineering prototype.(a) Photo of the prototype, (b) system block diagram.

    The experiment was carried out in a two-stage vacuum system, as described by Yang et al [13].The size of the square vacuum chamber in which the thruster is installed is 0.45 × 0.45 × 0.45 m3, and the size of the primary vacuum cylinder chamber is ?0.8 × 1.6 m (figure 2(a)).The vacuum chamber pressure (pvc) was measured with a ReBorn precision ionization gauge, model ZJ-12/CF35, with the measuring range from 1 × 10?1to 5 × 10?8Pa(1 × 10?3to 5 × 10?10mbar).This gauge was located at a flange on the side of the cylinder chamber approximately 1.2 m downstream from the thruster exit.The background neutral pressure could reach 5 × 10?5Pa (5 × 10?7mbar), andpvcis about 2 × 10?4Pa (2 × 10?6mbar) while the thruster is operational at the flow rate of 1 sccm xenon or the mass flow rate of 100 μg s?1iodine.

    2.2.Experimental methods

    In the experiment, only the anode/thruster temperature at ignition (Ta0) was changed, while the parameters such as the iodine flow rate out from the storage tank (˙mt), the incident radio frequency power (RFin), the screen grid voltage (Usc),the electron current emitted by cathode/neutralizer (Ice) and valve opening time remained unchanged.By monitoring the changes of screen grid current (Isc), the reflected radio frequency power (RFre),RFin,pvcand other parameters, the discharge performance of the thruster under differentTa0was analyzed, and the influence mechanism of iodine deposition and re-sublimation on thruster ignition was revealed.It was expected to obtain the minimum ignition temperature of the thruster (Tis) and the discharge over a wide range of temperature changes, which provides a basis for the design of thruster preheating method.

    The selection ofTa0was considered as follows.Obviously, ifTa0was the same as the temperature of the pipeline (Tline), the path in the thruster was unblocked and it could be ignited normally.Therefore,Tlinecould be used as the maximum temperature of the experiment.Meanwhile,IRIT4 could not ignite at room temperature(RT)according to our previous experiments, soRTcould be used as the lowest temperature of the experiment.Based on the above considerations,the range ofTa0in the experiment was set betweenRTandTline.Therefore,the experiments were divided into the following four cases.

    Figure 2.Experimental setup.(a) Schematic of the setup, (b) schematic of the thruster and the iodine feed module.

    Case 1.Ignition at low temperatures.Ta0was slightly higher thanRT.

    Case 2.Ignition at medium temperatures.Ta0was higher than the temperature in Case 1.Iodine deposit in the anode is not too serious and it could be ignited.However, the upper and lower limits of the temperature were unknown, and multiple attempts were made from lower temperature to higher one.

    Case 3.Ignition at high temperatures.Ta0was higher than the upper limit temperature in Case 2, and the flow path in the anode was completely unblocked.By analyzing thepvccurve at differentTa0, the lower limit of temperature in this case could be obtained.But in fact, it was found that this value could be predicted from the experimental results of Case 2.

    Case 4.Continuous firing at the high temperature.When the thruster is unblocked, ignites for a long time, and discharge performance was investigated during the process of the anode/thruster temperature (Ta) approaching the thermal equilibrium temperature of the thruster.

    Although multiple experiments were carried out, only one of the typical results was given in each case in this paper.

    For each of the above cases, a complete ignition test process includes the following 6 steps (figure 3).

    Step 1.Preheating the iodine feed module.Special heaters were used to heat the storage tank and the pipeline respectively.Closed loop control was adopted, and the temperature control accuracy reached ±0.1 °C.According to our previous flow test results [13], the temperature of the storage tank (Ttank) andTlinewere set to 63.8 °C and 85 °C respectively, with ˙mtof 100 μg s?1representing the nominal set point.Under this condition, the flow path inside the pipeline was completely unblocked.

    Figure 3.Steps of a typical ignition test.

    Step 2.Preheating the thruster (figure 2(b)).The anode was heated by RF power.The discharge chamber was directly heated by the RF coil,and tube A was heated by the discharge chamber.A PT1000 thermometer was installed on tube A to measure its temperature in real time.RFinwas set to about 23 W, andTarose slowly fromRT.RF power was turned off whenTaarrived at the target temperature ofTa0.WhenTawas slightly lower thanTa0,RF power was turned on manually again.Despite the measurement accuracy of PT1000 was as high as±0.1°C,the control error ofTa0was about±2°C due to the open-loop control.

    Step 3.pvctest.The valve was installed on the pipeline and close to the outlet of the storage tank (figure 2(b)).When the valve was opened, the iodine vapor in the storage tank flew into the anode through the pipeline, and then diffused from the anode outlet to the vacuum chamber.After the valve was closed, the iodine vapor flow was cut off and no longer flew into the pipeline and the anode.Whenpvcin the vacuum chamber decreased to about 2 × 10?4Pa (2 × 10?6mbar),the valve was opened and closed after about 30 s.At the rated flow of ˙ma= 100 μg s?1,pvcwas greater than 1 × 10?3Pa(1 × 10?5mbar).According to the relationship betweenpvcand ˙ma(see appendix A), it can be inferred that there was deposition in tube A of the anode or not.If tube A was clear,pvcrose rapidly at the moment of valve opening; moreover,pvcdecreased rapidly at the moment of valve closing.If thepvcrose slowly when the valve was opened,or it fell after the valve closing, it indicated that there was deposition or blockage in tube A.After closing the valve,the next step was proceeded until thepvcwas reduced to about 2 × 10?4Pa(2 × 10?6mbar) again.

    Step 4.Ignition.Whenpvcdropped to about 2 × 10?4Pa(2 × 10?6mbar) andTareachesTa0planned in step 2, the valve was opened again.Ignition was tried about 10 s later,whileRFinremained at 23 W,Ice,Uscand the acceleration grid voltage(Uac)were set to 30 mA,1400 ± 10 V and 200 ± 2 V respectively.A successful ignition or not was judged by the change ofIsc.An instantaneously increase ofIscindicated that a successful ignition and step 5 was carried out.Isc= 0 meant failure to ignition, and ignition was tried once more to eliminate the influence of random factors.After three unsuccessful ignition attempts within 5 s, step 6 would be carried out.

    Step 5.Continuous firing.Different from the closed-loop control mode of constant thrust on orbit, the thrust was not actively controlled during the discharge of the experiment in this work.Duration of the firing was set to be large than 600 s.When the thruster was operated under the rated parameters,Iscshould peak at about 32 mA instantaneously and remain constant, andRFreshould drop to below 0.2 W instantaneously.IfIscandRFrewere very different from the abovementioned values,Tahad a significant impact on the discharge process.

    Step 6.Shutting down and cooling down.The power supply of the screen grid, the accelerator grid, the RF generator, and the neutralizer were powered off in sequence, and the valve was closed finally.The thruster was cooled to room temperature, and the iodine feed module was kept warm or cool as needed.

    In each experiment, all the key parameters were measured in real time by the PPU with integrated controller and logged automatically through a special control software at a frequency of 1 Hz, includingIsc,Usc,Uac,RFin,RFre,Ice,Ta,Ttank,Tline, and the valve switching state (VS), etc.The net input radio frequency power of the thruster is denoted asRF,which is equal to the difference betweenRFinandRFre.VS= 0 means the valve is closed, andVS= 1 denotes it is opened.pvcmeasured by the ionization gauge was recorded by the vacuum chamber control software at a frequency of 1 Hz.

    3.Results and discussion

    3.1.Case 1: Ignition at low temperatures

    Figure 4.Failure of ignition at Ta0 = 27.1 °C.(a) Is c ,(b) RFi n ,RFr e,RF, (c) pvc andVS, (d)Ta.

    Figure 4 shows the time histories ofIsc,RF,RFin,RFre,pvc,VSandTawhen the ignition try was atTa0= 27.1 °C.From 0 to 38 s, only small oscillations inIscwere observed (notice that there was an initial value of about 3 mA on theIsccurve),and the ignition failed.RF,RFinandRFrewere also almost unchanged.pvcindicated 1.8 × 10?4Pa before opening the valve.When the valve was opened,pvcincreased rapidly,reaching 5 × 10?4Pa (5 × 10?6mbar) after 22 s, and then increased slowly.In particular, when the valve was closed at 42 s,pvcdid not drop immediately, but continued to rise slowly.This indicates that iodine has deposited inside the thruster and a clog was produced.The iodine in the clog served as the second source of iodine sublimation.As it was elevated in temperature, although the iodine vapor flow from the storage tank was cut off after the valve was closed, the deposited iodine continued to volatilize and enter the vacuum chamber at a lower flow rate, which is similar to the phenomenon in the glass tube described by Polzin [16].

    Tareached 28.9 °C at 117 s, and the valve was opened again at this time.pvcincreased only by 0.3 × 10?4Pa in 15 s,and the ignition failed again.After the valve was closed,pvcremained unchanged for 50 s.It was until 185 s thatpcbegan to drop rapidly, and it decreased to 3 × 10?4Pa(3 × 10?6mbar) in 20 s, which indicated that the clog was cleared and the flow rate reduced to zero.

    It can be seen that,the iodine vapor flowing from the tank through the pipeline will deposit inside the thruster and form a blockage whenTais too low.This causes ˙mato be much smaller thanm˙t,and it cannot be discharged under the rated ignition parameters ofRF,UscandIce,etc.This shows that there is a minimum ignition temperatureTi0for the iodine thruster.IfTa0is lower thanTi0due to the external environment of the satellite in orbit, the thruster cannot work, and it needs to be actively heated.The value ofTi0will be obtained in section 3.2.

    On the other hand,because the outlet of the thruster is in a high vacuum zone,the thruster has a certain degree of selfclearing ability; even at a lower temperature, all the iodine deposited inside the thruster will volatilize into the vacuum after a certain period of time.

    A simple model of the flow rate in the anode was established to analyze the above phenomena.˙macan be expressed by the following formula:

    where, ˙mdepis the iodine flow rate deposited in the anode,˙msubis the iodine flow rate sublimated from the anode.˙mtis a constant in this article, both ˙mdepand ˙msubare related to the thruster temperature,ambient air pressure,area,etc,which are given by equations (2)–(4) [17].

    where,pvapis the vapour pressure of the solid iodine,pis the partial pressure of iodine in the gas phase,Tais the absolute temperatureTa,Ris the gas constant,αis a constant,Mis the molecular weight, A, B, and C are constants.It should be noted that although the same symbols ofpandαare used in equations (2) and (3), they may be different.

    It can be seen from equations (2) and (3) that ˙mdepincreases with the increasingp,while ˙msubdecreases.For the thruster,it is obvious thatpin tube A is much larger than that in the discharge chamber, which leads to easier to deposit in tube A than the latter.Therefore, equation (1) is mainly used to describe the situation in tube A.

    Figure 5.Pressure in the vacuum chamber at Ta of 35.8 °C.(a) pvc andVS, (b)Ta.

    From equations (2) and (3), the lowerTa,the largerm˙depand the smallerm˙sub.During the progress of ignition in Case 1,m˙dep>m˙sub,thusm˙a

    According to equation (3), sublimation will continue as long aspis less thanpvap.According to formula(4),thepvapof iodine at room temperature is about 94 Pa,while thepvcin the vacuum chamber is about 1 × 10?4–1 × 10?3Pa(1 × 10?6–1 × 10?5mbar) during the experiment.After the valve is closed, the deposited iodine will continue to volatilize, that is, the vacuum environment gives the thruster a certain dredging ability.

    3.2.Case 2: Ignition at medium temperatures

    Prior to the ignition test at a higher temperature, a dedicatedpvctest was done, and the results are presented in figure 5.Notice that valve opening time is taken as 0 time in the figure.The thruster was heated to 35.8°C and the valve was opened for 30 s.pvcrose to 1.0 × 10?3Pa(1.0 × 10?5mbar)in 14 s and reached 1.2 × 10?3Pa (1.2 × 10?5mbar) in 30 s.Its rising slope and maximum value are significantly higher than those ofTa= 27.1°C.After the valve was closed at 30 s,pvcdid not drop immediately, but still rose slowly.pvcpeaked at 1.3 × 10?3Pa (1.3 × 10?5mbar) at 47 s, and then dropped rapidly afterwards.It is shown that there was still iodine deposition in the thruster atTa= 35.8 °C, but the clogging was less than that ofTa= 27.1 °C.

    Figure 6.Discharge with clog-clearing at Ta0 = 38.2 °C: (a) Is c ,(b)RFi n, RFr e ,RF, (c) pvc andVS, (d)Ta.

    Figure 7.Pressure in the vacuum chamber at Ta0 = 50.9 °C.(a) pvc andVS, (b)Ta.

    Results ofTa= 38.2°C are shown in figure 6.The valve is opened 9 s before ignition.pvcrose quickly to 1.1 × 10?3Pa (1.1 × 10?5mbar), and dropped rapidly when it was discharged (appendix A), as illustrated in figure 6(c).Significant variation of rising-falling-stable was observed in figure 6(a) which plotsIscversus time t during discharge of 850 s.After successful discharge,Iscincreased to 27 mA almost instantly, and then linearly increased, peaking at 36 mA at 240 s.Since the test was not fully automated,manual adjustment was used to tuneRFinandUsc.Isc= 17 mA at the beginning corresponded toUsc= 800 V,and it arrived at 27 mA instantly after Uswas adjusted to 1400 V at 10 s.From 240 to 320 s,Iscdecreased from the maximum value of 36 to 33 mA.Iscstabilized to 33 mA after 320 s firing.Under the condition thatUscand other parameters remain unchanged,variation ofIscis mainly determined bym˙aandRFfor the radio frequency ion thruster.Iscincreases with the increasingm˙aorRF, on the contrary, a decrease inIscmeans a decrease of them.

    Oscillations also appeared on theRFin,RFreandRFcurve within 310 s after discharge, ranging from 22.3 to 24.2 W,with a maximum change of 1.9 W, and finally stabilized to 23.2 W.In particular dramatic oscillations occurred inRFrecurve, showing the decrease–increase–decrease–stable trend,as shown in figure 6(b).The discharge process can be divided into 4 stages according to the change ofRFre.Stage I:0–110 s,RFredecreased from 0.8 to 0.0 W; Stage II: 110–240 s,RFreincreased again to 1.2 W;Stage III:240–310 s,RFredecreased to 0.1 W; Stage IV: After 310 s,RFrestabilized to 0.1 W.Here, the change ofRFrewas mainly attributed to the change ofm˙a(appendix B).The above mentioned four stages are discussed in detail below combined with the relationship betweenm˙aandRFre.

    Phase I.At the moment of ignition 0 s, the iodine flow rate reached the discharge condition.Since a part of the iodine was redeposited in the thruster,m˙a

    Phase II.RFreincreased again (from 0 to 1.1 W), whileIsccontinued to increase by 6 mA (from 30 to 36 mA), andTaincreased from 40.6 °C to 47.8 °C.SinceUscand other parameters remained unchanged in the experiment, the increase ofIscwas mainly caused by the increase of ˙maand/orRF.From appendix C,an increase of 2 W(from 22 to 24 W)inRFcontributes about 2 mA toIscnearIsc= 30 mA.This confirms that the additional 4 mA of the increase inIscis caused by the increase inm˙a,that is, the increase inRFrein phase II corresponded to the increase inm˙a.The reason for the increasingm˙ais that asTafurther increased,m˙sincreased,and the clogged flow path gradually became clearing.At this time,m˙t-m˙dep+m˙sub>m˙t.

    Phase III.RFreandIscdecreased simultaneously withTaincreasing from 47.8°C to 52.8°C.This indicates that ˙mawas gradually decreasing because the iodine deposited in the thruster is gradually sublimating and exhausting, although it was still greater thanm˙t.

    Phase IV.Except thatTacontinued to increase (from 52.8°C to 76.9 °C),RFreandIsctended to stabilize.The iodine deposited inside the thruster has been fully utilized, and the flow remained steady atm˙a=m˙t.

    In short, it can be ignited whenTa0reachesTi0.But there are large oscillations atIscandRFre,which takes about a few minutes to stabilize.The thrust change is contributed to the change of ˙maduring the process of the iodine deposition and clearing inside the thruster.WhenTais low, the ˙mais lower than ˙mtdue to the deposition of iodine,resulting in a smallIscand a largeRFre.AsTarises,the deposited iodine serves as the second sublimation source, which may causem˙a=m˙t-m˙dep+m˙sub>m˙t,and bothIscandRFrewill reach larger values.

    Although figure 6 shows the result ofTa0= 38.2 °C, the value ofTi0of the thruster can be obtained by a detailed test.The specific operation method is as follows: heating the thruster to keepTarising continuously and slowly,try to ignite once when the temperature rises by 1°C,until the discharge is successful; The lowest value ofTaobtained upon successful discharge isTi0.Finally, theTi0of the IRIT4-2D is 36 °C,which is much lower thanTlineof 85 °C andTtankof 63.8 °C.For the condition ofTa≥Ti0, although the thrust fluctuates dramatically with manual open-loop controls,the stable thrust output can be easily obtained if feedback linking thrust to RF power is adopted, which is one of the advantages of the RF ion thruster.

    Also, it can be predicted that whenTa0is greater than a certain temperatureTis, there will be no clog, and a stable discharge curve can be obtained.WhenIscreaches the maximum value in Phase II,the thruster is completely unblocked,and the correspondingTavalue at this time isTis.It can be obtained from figure 6 thatIscpeaked at 36 mA at 240 s(figure 6(a)),andTawas 47.8°C(figure 6(d))at this time.So it can be inferred thatTis= 47.8 °C.

    3.3.Case 3: Ignition at high temperatures

    Tawas set to be slightly higher than 47.8 °C according to theTisinferred from section 3.2.Similar to the previous section,a specialpvctest was done firstly to get the change ofpvcwhen the valve was opened and closed.Figure 7 shows thepvccurve ofTa= 50.9 °C.After the valve was opened,pvcrose rapidly, reaching 1.0 × 10?3Pa (1.0 × 10?5mbar) at 15 s and 1.2 × 10?3Pa(1.2 × 10?5mbar)at 34 s.The valve was closed at 34 s andpvcdrops immediately.It took only 5 s

    (from 34 to 39 s) forpvcto drop from 1.2 × 10?3Pa(1.2 × 10?5mbar)to 4.0 × 10?4Pa(4.0 × 10?6mbar).The ups and downs ofpvcwere completely synchronized with the opening and closing of the valve.This indicates that there was no iodine deposit or blockage inside the thruster, and the iodine flow path was completely cleared.This conforms with the prediction of the temperature required for no iodine deposits in the previous section.

    Figure 8 shows the result of discharge atTa0= 48.4 °C.The valve was opened 10 s before discharge.pvcrose rapidly when the valve was opened,and then fell rapidly immediately after the discharge-the trend was the same as that of Case 2.Tacontinued to rise from 48.4 °C to 84.9 °C within 900 s.Disregarding the two steps ofIscin the initial 22 s due to the manual adjustment ofUscandRF,IscandRFquickly stabilized to the maximum values of 32.0 mA and 22.3 W,respectively, whileRFrestabilized to the minimum value of 0.1 W.

    It can be seen that the discharge curve is stable whenTa0≥Tis, andTiscan be called the steady discharge temperature.For the thruster used in this work,Tisis 48.4°C,just as inferred in section 3.2.AlthoughTisis 12°C higher thanTi0,it is still much lower thanTlineof 85 °C andTtankof 63.7 °C.There are two reasons for this:first,the outlet of the thruster is directly facing the vacuum, which leads to the low back pressure and the large sublimation flow rate; second, the larger internal diameter and the shorter flow channel of the thruster than the pipelines result in the smaller flow resistance.It should be noted that the valueTisis related to the specific internal structure of the thruster, and it will be different for different thruster structures.Stable discharge can be achieved whenTais much lower thanTline, which is of great significance for reducing the preheating power consumption of the thruster and shortening its preheating time.

    Figure 8.Normal operation at Ta0 = 48.4 °C.(a) Is c ,(b) RFi n, RFr e,RF, (c) pvc andVS, (d)Ta.

    3.4.Case 4: Continuous firing at the high temperature

    Figure 9.Continuous firing at Ta0 = 48.4 °C.(a) Is c ,(b) RFi n ,RFr e,RF, (c) pvc andVS, (d)Ta.

    For the experiment corresponding to figure 8 in Case 3, the thruster was kept firing continuously for 4000 s to test its long time scale stability with the large increases ofTa,and the result is shown in figure 9.AlthoughTaincreased from 48.4°C to 142.7 °C, the screen grid current and RF power remained stable, withIsc,RFandRFrebeing 32 mA, 23.5 W and 0.1 W,respectively.It can be seen that during the hourly continuous operation of the thruster,the thrust performance is hardly affected byTaeven if the temperature changes up to 94.3 °C.That is to say, whenTais high enough, the internal iodine path of the thruster is completely unblocked; the thermal throttling effect of the internal flow path of the thruster can be ignored so as not to reduce the flow rate, and the thrust can remain stable for a long time under the condition of a large temperature rise.

    4.Conclusions

    In this paper, the experimental results of the IRIT4-2D IRIT ignited at different thruster/anode temperatures are presented.It is found that the inherent cold deposition characteristics of iodine make its ignition characteristics different from that of xenon.The main conclusions are as follows:

    (1) There is a lowest ignition temperatureTi0for the iodine thruster.AtTa

    (2) WhenTarises toTi0(36 °C), it can be ignited.ButIscandRFreoscillate dramatically in the first few minutes,which is mainly caused by the flow changes of the iodine vapor during the process of iodine deposition-toclearing in the thruster.At the beginning,m˙a

    (3) WhenTais higher thanTis(47.8 °C), the discharge process is stable.Tisis much lower thanTlineor even lower thanTtank,which can be attributed to the low back pressure and low flow resistance of the thruster facing the vacuum.

    (4) During the continuous operation of the thruster, the thrust performance is hardly affected by the temperature of the thruster even if the temperature changes up to about 100 °C.

    Notice that the values ofTi0andTismay vary for different internal structure of the thruster, but the ignition process will be similar.Research is helpful for reducing the preheating power consumption of the thruster, shortening the preheating time and designing the thrust control strategy.

    Acknowledgments

    The financial support from National Natural Science Foundation of China (No.11805265) and Key Laboratory of Micro-Satellites, Chinese Academy of Sciences (No.KFKT201903) is gratefully acknowledged.

    Figure B1.Variation of RF power with Ttank (˙ma) at Isc = 31.5 mA.

    Appendix A

    Relation betweenpvcand ˙maare briefly introduced in this section.

    Since the actual iodine flow rate exiting from the thruster is proportional topvcmeasured by the ionization gauge, the real-time ˙macan be obtained from thepvccurve.

    Rapidly decreasingpvcafter the successful ignition does not mean that the propellant flow rate is reduced, which is caused by the working principle of the ionization gauge ZJ12.The hot cathode in the ionization gauge emits electrons to ionize gas molecules, then the ions are collected by the collector, and the gas pressure is measured by the collected ion current.When the thruster is successfully discharged,most of the iodine is ionized and accelerated to the velocity of about 20 km s?1.Particles ofI+/I2/I/I2+/I2+etc impacts on the wall surface of the vacuum chamber,and a complicated sputtering reaction occurred.Only a small amount of iodine arrive at the site of the ionization gauge, resulting in that the measured pressure is very low.(The same phenomenon occurs on xenon in our experiments.)

    Appendix B

    Since a fixed frequency of 2 MHz and a fixed matching network were applied in the system,RFreis mostly affected by the change ofm˙a.This is illustrated by figure B1,which pl otsRF,RFinandRFreagainstTtank(˙ma)while keeping bothIscandUscconstant, wherem˙a=m˙tat high enoughTa.There is a downward-upward trend with the increasingTtank(˙ma), with the minimumRFreappears atTtank= 64.3 °C (˙ma= 110 μg s?1)which is close to the nominal flow rate of the system as designed.

    Appendix C

    A special test was arranged to get the contribution ofRFtoIsc,and the result is shown in figure C1.(Notice thatIscwas time-averaged in figure C1(a).) Keepingm˙a=m˙tat high enoughTa,theIsc?RFcurve was obtained by adjustingRFin.It can be seen that, nearIsc= 30 mA, an increase of 2 W(from 22 to 24 W) inRFcontributes about 2 mA toIsc.

    Figure C1.Isc changes with RF of non-blocking discharge.(a)Is c, (b)RF, (c) Isc?RF.

    ORCID iDs

    猜你喜歡
    李龍興華振華
    家住西安
    毛焰藝術(shù)風(fēng)格中的自我表達(dá)
    興華市林湖鄉(xiāng):村企聯(lián)建共走振興路
    高明的診斷
    攝影作品欣賞
    金沙江文藝(2019年7期)2019-07-29 01:57:06
    高明的診斷
    故事會(huì)(2019年8期)2019-04-20 12:48:04
    Cole-Hopf Transformation Based Lattice Boltzmann Model for One-dimensional Burgers’Equation?
    馬興華攝影作品欣賞
    金沙江文藝(2017年4期)2017-03-31 07:35:16
    “杯”慘
    我想用鮮花修補(bǔ)已不再湛藍(lán)的天空
    椰城(2014年6期)2014-11-17 08:43:58
    天天一区二区日本电影三级| 国产亚洲欧美98| 听说在线观看完整版免费高清| 午夜影院日韩av| 精品不卡国产一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 97碰自拍视频| 精品99又大又爽又粗少妇毛片 | 麻豆成人午夜福利视频| 亚洲欧美日韩无卡精品| 99国产精品一区二区三区| 特大巨黑吊av在线直播| 国产成人aa在线观看| 欧美成人性av电影在线观看| 99国产综合亚洲精品| 国产成年人精品一区二区| 精品久久久久久,| 婷婷精品国产亚洲av在线| 99精品在免费线老司机午夜| 成人性生交大片免费视频hd| АⅤ资源中文在线天堂| 桃色一区二区三区在线观看| 久久精品影院6| 在线a可以看的网站| 亚洲精品乱码久久久v下载方式 | 国产高潮美女av| 一区二区三区国产精品乱码| 操出白浆在线播放| 国内毛片毛片毛片毛片毛片| 日本与韩国留学比较| 国产精品亚洲一级av第二区| 看片在线看免费视频| 人妻久久中文字幕网| 成年人黄色毛片网站| 国产三级黄色录像| 国产精品av视频在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 精品久久久久久久久久久久久| 午夜福利成人在线免费观看| 舔av片在线| 少妇丰满av| 日韩免费av在线播放| 老司机在亚洲福利影院| 高清在线国产一区| av女优亚洲男人天堂 | 色精品久久人妻99蜜桃| 999久久久国产精品视频| 婷婷精品国产亚洲av| 波多野结衣巨乳人妻| 一级作爱视频免费观看| 久9热在线精品视频| 69av精品久久久久久| 观看美女的网站| 19禁男女啪啪无遮挡网站| 99久久精品热视频| 日韩精品青青久久久久久| 男人和女人高潮做爰伦理| 久久婷婷人人爽人人干人人爱| 黄色日韩在线| 99久国产av精品| 91av网站免费观看| 又黄又粗又硬又大视频| 亚洲成人久久爱视频| 免费观看人在逋| 在线观看美女被高潮喷水网站 | 老司机深夜福利视频在线观看| 精品午夜福利视频在线观看一区| xxxwww97欧美| 桃红色精品国产亚洲av| 久久香蕉国产精品| 网址你懂的国产日韩在线| 又爽又黄无遮挡网站| 成在线人永久免费视频| 亚洲专区国产一区二区| 99国产精品99久久久久| 亚洲第一欧美日韩一区二区三区| 琪琪午夜伦伦电影理论片6080| 草草在线视频免费看| 国产乱人视频| 黄色成人免费大全| 国产成人av激情在线播放| 国产成人aa在线观看| 久久久久久国产a免费观看| av福利片在线观看| 国产日本99.免费观看| 麻豆国产av国片精品| 午夜福利视频1000在线观看| 免费看十八禁软件| 精品久久久久久久末码| 久久九九热精品免费| 精品熟女少妇八av免费久了| 国模一区二区三区四区视频 | 国产91精品成人一区二区三区| 国产精品久久久久久久电影 | 日韩 欧美 亚洲 中文字幕| 日韩人妻高清精品专区| 国内精品久久久久精免费| 窝窝影院91人妻| 又爽又黄无遮挡网站| 欧美高清成人免费视频www| xxx96com| 亚洲精品乱码久久久v下载方式 | 69av精品久久久久久| 日韩av在线大香蕉| 舔av片在线| 全区人妻精品视频| 午夜日韩欧美国产| www国产在线视频色| 亚洲九九香蕉| 少妇丰满av| 成年女人看的毛片在线观看| 小蜜桃在线观看免费完整版高清| 天堂网av新在线| 国产激情久久老熟女| 欧美另类亚洲清纯唯美| 久久国产精品人妻蜜桃| 老司机午夜福利在线观看视频| 成人鲁丝片一二三区免费| 欧美成人一区二区免费高清观看 | 淫妇啪啪啪对白视频| 中文字幕熟女人妻在线| av欧美777| 亚洲中文字幕一区二区三区有码在线看 | 免费大片18禁| 99在线人妻在线中文字幕| 18禁国产床啪视频网站| 在线视频色国产色| 男女下面进入的视频免费午夜| 亚洲美女黄片视频| www.自偷自拍.com| 亚洲 国产 在线| 国内精品美女久久久久久| 在线a可以看的网站| 特大巨黑吊av在线直播| 久久精品91无色码中文字幕| 午夜福利视频1000在线观看| 99久久成人亚洲精品观看| 一a级毛片在线观看| 国产午夜精品久久久久久| 最近最新免费中文字幕在线| 国产亚洲欧美98| 丁香六月欧美| 精品国产超薄肉色丝袜足j| 国产aⅴ精品一区二区三区波| 欧美成人一区二区免费高清观看 | 脱女人内裤的视频| 欧美日韩精品网址| 色在线成人网| 亚洲精品国产精品久久久不卡| 丝袜人妻中文字幕| 欧美极品一区二区三区四区| 看免费av毛片| 精品不卡国产一区二区三区| 村上凉子中文字幕在线| 床上黄色一级片| 18禁裸乳无遮挡免费网站照片| 99久国产av精品| 国产成人精品久久二区二区免费| 久久99热这里只有精品18| 日本一二三区视频观看| 久久久久久久久中文| 美女午夜性视频免费| 国产成人精品久久二区二区91| 中文在线观看免费www的网站| 精品一区二区三区视频在线观看免费| 偷拍熟女少妇极品色| 国产精品 欧美亚洲| 小说图片视频综合网站| 成人国产一区最新在线观看| 岛国在线免费视频观看| 综合色av麻豆| 久99久视频精品免费| 怎么达到女性高潮| 国产主播在线观看一区二区| 亚洲欧美激情综合另类| 亚洲精品在线美女| 狠狠狠狠99中文字幕| 国产综合懂色| 91av网一区二区| 亚洲五月婷婷丁香| 国产精品99久久久久久久久| 99在线人妻在线中文字幕| 国产激情久久老熟女| 亚洲熟女毛片儿| 嫩草影院精品99| 欧美国产日韩亚洲一区| 波多野结衣高清无吗| 母亲3免费完整高清在线观看| 啦啦啦韩国在线观看视频| 亚洲成人中文字幕在线播放| 一进一出抽搐动态| 一区二区三区激情视频| 91在线精品国自产拍蜜月 | 日韩欧美在线乱码| 欧美黑人欧美精品刺激| 亚洲国产中文字幕在线视频| 亚洲精品一区av在线观看| 母亲3免费完整高清在线观看| 日日夜夜操网爽| 国产又色又爽无遮挡免费看| 亚洲成人久久爱视频| 欧美不卡视频在线免费观看| 嫩草影院精品99| 久久久久精品国产欧美久久久| www.自偷自拍.com| 久久午夜亚洲精品久久| 在线十欧美十亚洲十日本专区| 天堂动漫精品| 亚洲第一电影网av| 亚洲成人中文字幕在线播放| 非洲黑人性xxxx精品又粗又长| 一二三四在线观看免费中文在| 麻豆国产97在线/欧美| 人妻久久中文字幕网| 国产高潮美女av| 日韩av在线大香蕉| 成人永久免费在线观看视频| 久久精品综合一区二区三区| 男女那种视频在线观看| 久久精品aⅴ一区二区三区四区| 亚洲人与动物交配视频| 999久久久精品免费观看国产| 午夜亚洲福利在线播放| 亚洲av中文字字幕乱码综合| 久久久久久久久久黄片| 日本在线视频免费播放| 99国产极品粉嫩在线观看| www.www免费av| 91在线观看av| 国产精品九九99| 色尼玛亚洲综合影院| 国产精品影院久久| 午夜精品久久久久久毛片777| 老熟妇乱子伦视频在线观看| 精品国产乱码久久久久久男人| 国产成人啪精品午夜网站| 国产精品久久久久久亚洲av鲁大| 十八禁人妻一区二区| 又黄又粗又硬又大视频| 国产精品爽爽va在线观看网站| 久久午夜综合久久蜜桃| 在线视频色国产色| 日韩有码中文字幕| 色哟哟哟哟哟哟| 免费观看人在逋| 免费在线观看日本一区| 国产精品一及| 精品久久久久久成人av| 一级a爱片免费观看的视频| 在线观看日韩欧美| 999久久久精品免费观看国产| 亚洲精品中文字幕一二三四区| 美女被艹到高潮喷水动态| 九九久久精品国产亚洲av麻豆 | 欧美高清成人免费视频www| 日韩中文字幕欧美一区二区| 熟妇人妻久久中文字幕3abv| 亚洲成人免费电影在线观看| 久久精品夜夜夜夜夜久久蜜豆| 欧美性猛交╳xxx乱大交人| 美女cb高潮喷水在线观看 | 日韩精品青青久久久久久| 亚洲精品国产精品久久久不卡| 久久久色成人| 高潮久久久久久久久久久不卡| ponron亚洲| 俄罗斯特黄特色一大片| 国产单亲对白刺激| 国产精品一区二区三区四区免费观看 | 国内精品一区二区在线观看| 每晚都被弄得嗷嗷叫到高潮| 欧美黑人欧美精品刺激| 麻豆国产av国片精品| 精品一区二区三区av网在线观看| 老汉色∧v一级毛片| 黄片小视频在线播放| 真人一进一出gif抽搐免费| 国产伦精品一区二区三区四那| 12—13女人毛片做爰片一| 成年女人永久免费观看视频| 99热6这里只有精品| 久久精品亚洲精品国产色婷小说| 少妇人妻一区二区三区视频| 国产精品 国内视频| 成人永久免费在线观看视频| 国产精品精品国产色婷婷| 99热精品在线国产| 天堂√8在线中文| 国产精品永久免费网站| 婷婷六月久久综合丁香| 国产成人啪精品午夜网站| 亚洲美女黄片视频| 国产美女午夜福利| 大型黄色视频在线免费观看| 在线观看免费午夜福利视频| 老熟妇仑乱视频hdxx| 午夜福利视频1000在线观看| 亚洲中文字幕日韩| 日本一二三区视频观看| 夜夜看夜夜爽夜夜摸| 天天躁日日操中文字幕| 一个人看的www免费观看视频| 黄频高清免费视频| 观看免费一级毛片| 日本免费一区二区三区高清不卡| 欧美不卡视频在线免费观看| 日本黄色视频三级网站网址| 久久这里只有精品中国| 两个人视频免费观看高清| 久久九九热精品免费| 天堂网av新在线| 嫩草影视91久久| 欧美乱妇无乱码| 又紧又爽又黄一区二区| 在线观看午夜福利视频| 黄色成人免费大全| 18美女黄网站色大片免费观看| 不卡av一区二区三区| 国产精品av久久久久免费| 丁香六月欧美| 成人特级av手机在线观看| 人妻丰满熟妇av一区二区三区| 午夜a级毛片| 国产v大片淫在线免费观看| 人人妻人人澡欧美一区二区| 国产97色在线日韩免费| 无遮挡黄片免费观看| 国产日本99.免费观看| 伊人久久大香线蕉亚洲五| av女优亚洲男人天堂 | 天堂动漫精品| 黄色丝袜av网址大全| 欧美在线黄色| 国产麻豆成人av免费视频| 国产精品一区二区三区四区免费观看 | 国内毛片毛片毛片毛片毛片| 国产精品亚洲美女久久久| 嫩草影院精品99| 色综合婷婷激情| 久久精品综合一区二区三区| 啦啦啦观看免费观看视频高清| 宅男免费午夜| 久久久精品大字幕| 亚洲国产欧洲综合997久久,| 国内精品久久久久精免费| 亚洲,欧美精品.| 级片在线观看| 欧美黑人欧美精品刺激| 久久午夜综合久久蜜桃| 国产高清视频在线观看网站| 日日夜夜操网爽| 亚洲熟妇熟女久久| 99久久精品一区二区三区| 长腿黑丝高跟| 波多野结衣巨乳人妻| 哪里可以看免费的av片| 看片在线看免费视频| 国产三级黄色录像| 美女被艹到高潮喷水动态| 波多野结衣高清作品| 亚洲在线观看片| 男女床上黄色一级片免费看| 国产激情欧美一区二区| 床上黄色一级片| av黄色大香蕉| 国产乱人视频| 天天躁日日操中文字幕| 精品国内亚洲2022精品成人| 亚洲国产精品sss在线观看| 男插女下体视频免费在线播放| 亚洲专区字幕在线| 19禁男女啪啪无遮挡网站| 亚洲欧美激情综合另类| 国产精品久久久久久久电影 | 日本免费一区二区三区高清不卡| 日韩成人在线观看一区二区三区| 欧美日韩福利视频一区二区| 999久久久精品免费观看国产| 91av网站免费观看| 在线国产一区二区在线| 深夜精品福利| 日韩欧美国产在线观看| 日本 av在线| 欧美黄色淫秽网站| 人人妻,人人澡人人爽秒播| 国产精华一区二区三区| 天天一区二区日本电影三级| 成人特级av手机在线观看| 淫妇啪啪啪对白视频| 搡老岳熟女国产| 两性夫妻黄色片| 欧美不卡视频在线免费观看| 国产精品日韩av在线免费观看| 伦理电影免费视频| 最好的美女福利视频网| 给我免费播放毛片高清在线观看| 国语自产精品视频在线第100页| 欧美日韩乱码在线| 小蜜桃在线观看免费完整版高清| 国产蜜桃级精品一区二区三区| 亚洲欧美激情综合另类| 国产精品亚洲一级av第二区| 我要搜黄色片| 又黄又粗又硬又大视频| 少妇裸体淫交视频免费看高清| 精品国内亚洲2022精品成人| 国产av在哪里看| 亚洲 国产 在线| 91久久精品国产一区二区成人 | 制服丝袜大香蕉在线| 国产又黄又爽又无遮挡在线| 波多野结衣高清无吗| 母亲3免费完整高清在线观看| 脱女人内裤的视频| 久久久成人免费电影| 亚洲一区高清亚洲精品| 男人舔女人下体高潮全视频| 国产aⅴ精品一区二区三区波| 五月玫瑰六月丁香| 国产1区2区3区精品| 成在线人永久免费视频| 婷婷亚洲欧美| www.www免费av| 色哟哟哟哟哟哟| 亚洲人与动物交配视频| 一区二区三区高清视频在线| 亚洲国产欧美一区二区综合| 制服人妻中文乱码| 日本三级黄在线观看| 国产高潮美女av| 国产精品久久视频播放| 国产av麻豆久久久久久久| 丰满的人妻完整版| 中出人妻视频一区二区| 亚洲 国产 在线| 成人无遮挡网站| 中文资源天堂在线| 岛国视频午夜一区免费看| xxx96com| 一个人免费在线观看的高清视频| 又黄又爽又免费观看的视频| 国产亚洲av嫩草精品影院| 国产精品综合久久久久久久免费| 成人鲁丝片一二三区免费| 国内揄拍国产精品人妻在线| 欧美一级a爱片免费观看看| 国产午夜精品久久久久久| 欧美乱码精品一区二区三区| 国产精品 国内视频| 午夜福利在线观看免费完整高清在 | 嫩草影视91久久| 日本黄色片子视频| 国产蜜桃级精品一区二区三区| 国产成人av激情在线播放| 久久久国产精品麻豆| 免费看光身美女| 女同久久另类99精品国产91| 免费看十八禁软件| 午夜视频精品福利| 国产美女午夜福利| 国内精品美女久久久久久| 久久中文字幕人妻熟女| 岛国在线免费视频观看| 9191精品国产免费久久| 国产精品久久久人人做人人爽| 国产69精品久久久久777片 | 桃色一区二区三区在线观看| 动漫黄色视频在线观看| 欧美zozozo另类| 香蕉久久夜色| 叶爱在线成人免费视频播放| 日本成人三级电影网站| bbb黄色大片| 啦啦啦免费观看视频1| 观看免费一级毛片| 琪琪午夜伦伦电影理论片6080| 免费在线观看视频国产中文字幕亚洲| 18禁裸乳无遮挡免费网站照片| 日韩免费av在线播放| 老司机午夜福利在线观看视频| 久久精品国产清高在天天线| 制服人妻中文乱码| 欧美一区二区国产精品久久精品| 成人国产综合亚洲| 亚洲国产色片| 无遮挡黄片免费观看| 精品久久久久久久久久久久久| 国产精品九九99| 亚洲最大成人中文| 久久久色成人| 亚洲欧美日韩无卡精品| 男女那种视频在线观看| 国产精品av久久久久免费| 一二三四在线观看免费中文在| 亚洲人成电影免费在线| 色噜噜av男人的天堂激情| 国产成人av激情在线播放| 俄罗斯特黄特色一大片| 久久香蕉精品热| 老司机午夜十八禁免费视频| 黄色女人牲交| 亚洲中文日韩欧美视频| 最近最新免费中文字幕在线| 大型黄色视频在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产欧美人成| 老司机在亚洲福利影院| 18禁黄网站禁片午夜丰满| 一a级毛片在线观看| 免费观看精品视频网站| 欧美成人性av电影在线观看| 999久久久精品免费观看国产| 成年免费大片在线观看| 黄片小视频在线播放| 在线观看舔阴道视频| 亚洲精品久久国产高清桃花| 国产精品久久久久久亚洲av鲁大| 国内久久婷婷六月综合欲色啪| 琪琪午夜伦伦电影理论片6080| 欧美不卡视频在线免费观看| 日韩有码中文字幕| 热99在线观看视频| 波多野结衣巨乳人妻| 老鸭窝网址在线观看| 日韩人妻高清精品专区| 日本在线视频免费播放| 国产精品爽爽va在线观看网站| 国语自产精品视频在线第100页| 精品熟女少妇八av免费久了| 午夜a级毛片| 99久久精品一区二区三区| 91久久精品国产一区二区成人 | 美女免费视频网站| 免费在线观看日本一区| 久久久久久久精品吃奶| 国内少妇人妻偷人精品xxx网站 | 国产成人福利小说| 成年版毛片免费区| 久9热在线精品视频| 神马国产精品三级电影在线观看| av黄色大香蕉| 国产精品女同一区二区软件 | 亚洲欧美激情综合另类| 国产91精品成人一区二区三区| 人人妻人人看人人澡| 一边摸一边抽搐一进一小说| 免费看a级黄色片| 最近最新中文字幕大全电影3| 亚洲九九香蕉| 亚洲中文字幕一区二区三区有码在线看 | 精品国产亚洲在线| 免费看日本二区| 亚洲av成人av| 老司机在亚洲福利影院| 精品无人区乱码1区二区| 91麻豆精品激情在线观看国产| 欧美高清成人免费视频www| 丝袜人妻中文字幕| 在线视频色国产色| 欧美最黄视频在线播放免费| 99国产精品一区二区蜜桃av| 1000部很黄的大片| 亚洲第一欧美日韩一区二区三区| 91麻豆精品激情在线观看国产| 色综合婷婷激情| 十八禁人妻一区二区| 在线视频色国产色| 日日夜夜操网爽| 中亚洲国语对白在线视频| 99国产极品粉嫩在线观看| 两个人看的免费小视频| 日日干狠狠操夜夜爽| 日韩人妻高清精品专区| 国产人伦9x9x在线观看| 国产黄色小视频在线观看| 日韩欧美 国产精品| 久久伊人香网站| 久久中文字幕人妻熟女| 麻豆久久精品国产亚洲av| 亚洲国产中文字幕在线视频| 舔av片在线| 亚洲中文字幕一区二区三区有码在线看 | 久久久久久九九精品二区国产| 午夜激情福利司机影院| 国产精品久久久人人做人人爽| 香蕉久久夜色| 国产精品久久视频播放| 午夜福利18| 欧美最黄视频在线播放免费| 中亚洲国语对白在线视频| 久久人人精品亚洲av| 啦啦啦韩国在线观看视频| 中亚洲国语对白在线视频| 久久人人精品亚洲av| 99精品在免费线老司机午夜| 中国美女看黄片| 国产精品久久久av美女十八| 国产1区2区3区精品| 高清毛片免费观看视频网站| 一级a爱片免费观看的视频| 免费一级毛片在线播放高清视频| 99riav亚洲国产免费| 五月伊人婷婷丁香| 中文字幕av在线有码专区| 欧美成人性av电影在线观看| 久久欧美精品欧美久久欧美| 国产成人av激情在线播放| 听说在线观看完整版免费高清| 999精品在线视频| 欧美成人免费av一区二区三区| 成年女人永久免费观看视频| 日韩中文字幕欧美一区二区| 亚洲av电影不卡..在线观看| 成在线人永久免费视频| 99久久久亚洲精品蜜臀av|