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

    磁路和天線位置對(duì)2cm電子回旋共振離子推力器性能影響的實(shí)驗(yàn)研究*

    2019-12-16 11:39:06夏旭楊涓金逸舟杭觀榮付瑜亮胡展
    物理學(xué)報(bào) 2019年23期
    關(guān)鍵詞:離子源推力器磁路

    夏旭 楊涓? 金逸舟 杭觀榮 付瑜亮 胡展

    1) (西北工業(yè)大學(xué)航天學(xué)院,西安 710072)

    2) (上海空間推進(jìn)研究所,上海 201112)

    磁路和天線位置對(duì)電子回旋共振離子推力器(electron cyclotron resonance ion thruster,ECRIT)的電子加熱、等離子體約束和輸運(yùn)過(guò)程都有很大的影響,進(jìn)而影響離子束流引出和中和器耦合電壓.本文開(kāi)展不同磁路和天線位置下2 cm ECRIT離子束流和耦合電壓變化規(guī)律的實(shí)驗(yàn)研究.通過(guò)比較不同磁路的離子源和中和器的束流引出特性,選出合理的磁路結(jié)構(gòu),再比較不同天線位置對(duì)束流引出的影響.歸納了磁路和天線位置對(duì)ECRIT的性能影響規(guī)律,得到合理的推力器結(jié)構(gòu).實(shí)驗(yàn)結(jié)果表明:功率和流量的增加有助于提高離子引出束流和降低電子引出壓;保持天線空間位置不變,合理的磁路結(jié)構(gòu)能增大電子獲能并減小粒子損失,從而提高引出離子束流并降低耦合電壓;在合理磁路結(jié)構(gòu)條件下,離子源和中和器存在有利于離子引出和降低耦合電壓的合理天線位置.根據(jù)實(shí)驗(yàn)結(jié)果選擇出結(jié)構(gòu)較優(yōu)的中和器和離子源進(jìn)行中和實(shí)驗(yàn).結(jié)果表明:有無(wú)中和器工作時(shí)對(duì)離子源束流引出的影響較?。还β屎土髁繛? W,0.1 sccm (1 sccm = 1 mL/min)的中和器與功率和流量為2 W,0.3 sccm的離子源能良好匹配工作,性能指標(biāo)為離子束流5.3 mA、放電損耗337.5 W/A、推進(jìn)劑利用率24.7%、推力368.6 μN(yùn)、比沖1277.6 s、中和器耦合電壓17.4 V.研究結(jié)果有助于理解推力器工作機(jī)理,并為設(shè)計(jì)和性能優(yōu)化提供參考.

    1 引 言

    近年來(lái)隨著微型航天器技術(shù)的發(fā)展,微型電推力器也獲得蓬勃發(fā)展.這些微型電推力器有脈沖式等離子體推力器[1]、場(chǎng)效應(yīng)發(fā)射離子推力器[2]以及微型離子推力器等.微型離子推力器分為直流放電離子推力器[3-5]、射頻放電離子推力器[6,7]和ECRIT[8-10].其中微型電子回旋共振離子推力器(electron cyclotron resonance ion thruster,ECRIT)具有結(jié)構(gòu)簡(jiǎn)單、工作穩(wěn)定、比沖高、無(wú)陰極燒蝕的特點(diǎn),是未來(lái)值得發(fā)展的一種微型電推力器.微型ECR離子源除了應(yīng)用在航天領(lǐng)域,還能用于材料處理等方面[11,12].

    2 cm ECRIT屬于微型離子推力器,目前國(guó)內(nèi)外已經(jīng)對(duì)這類推力器開(kāi)展了大量研究,并于2014年把氙工質(zhì)穩(wěn)態(tài)工作的ECRIT分別應(yīng)用到HODOYOSHI-4與PROCYON兩顆小衛(wèi)星上[13,14].目前針對(duì)ECRIT進(jìn)行了放電引出實(shí)驗(yàn)、離子源內(nèi)等離子體參數(shù)診斷和數(shù)值模擬等研究工作[15-23].其中,引出實(shí)驗(yàn)方面的工作主要是文獻(xiàn)[24]實(shí)驗(yàn)研究了不同氙氣流量、功率下天線構(gòu)型對(duì)2 cm-μ1型ECRIT離子束流引出的影響,結(jié)果表明環(huán)形天線更有利于束流引出,而且引出的離子束流隨功率和流量的增加而增加.文獻(xiàn)[25]實(shí)驗(yàn)研究了功率和氙氣流量對(duì)雙模式2 cm-μ1型 ECRIT離子束流引出和耦合電壓的影響,結(jié)果表明兩種模式下離子和電子束流與功率和流量成正相關(guān)、中和器耦合電壓隨流量減小而增大.文獻(xiàn)[26]實(shí)驗(yàn)研究了不同腔體長(zhǎng)度對(duì)氬工質(zhì)的2 cm ECRIT離子源束流引出的影響,實(shí)驗(yàn)結(jié)果表明腔體長(zhǎng)度過(guò)短或過(guò)長(zhǎng)會(huì)導(dǎo)致離子束流引出的減小.文獻(xiàn)[27]實(shí)驗(yàn)研究了氙工質(zhì)2 cm 電子回旋共振(ECR)中和器的環(huán)形天線尺寸大小、電子引出板、流量和功率對(duì)電子束流引出的影響,得到了合理的天線和電子引出板結(jié)構(gòu)及引出電子束流隨功率和流量增加而增加的變化規(guī)律.

    目前對(duì)于氙氣或氬氣以及穩(wěn)態(tài)或雙模式工作的2 cm ECRIT研究主要集中在天線構(gòu)型和腔體結(jié)構(gòu)對(duì)束流引出、耦合電壓的影響.天線構(gòu)型的研究目的主要是匹配磁場(chǎng)與微波電場(chǎng),未考慮磁場(chǎng)會(huì)影響共振區(qū)與天線位置的關(guān)系.然而,由于共振區(qū)位置與磁路結(jié)構(gòu)有關(guān),而且?guī)щ娏W拥膿p失與天線位置有關(guān),所以磁路與天線位置是影響等離子體的加熱、壁面損失以及輸運(yùn)過(guò)程的重要因素,也是影響離子束流引出和耦合電壓以及推力器性能的重要因素.為此本文開(kāi)展磁路和天線位置對(duì)2 cm ECRIT離子束流和耦合電壓影響的實(shí)驗(yàn)研究,以獲得合理的離子源和中和器結(jié)構(gòu)參數(shù).研究結(jié)果有助于理解2 cm ECRIT的工作過(guò)程,為提升其性能提供參考.

    2 2 cm ECRIT結(jié)構(gòu)及電子獲能的影響因素分析

    2.1 2 cm ECRIT結(jié)構(gòu)

    2 cm ECRIT由離子源和中和器組成,二者均利用ECR原理加熱電子,產(chǎn)生的高能電子與中性氣體電離并維持等離子體.二者采用相同的放電室結(jié)構(gòu),放電室都主要由磁路、環(huán)形天線和圓柱腔體組成,其中磁路由內(nèi)外環(huán)形永磁體及背部磁軛構(gòu)成,在磁軛位于內(nèi)外磁環(huán)間環(huán)向均布進(jìn)氣孔,如圖1中虛線包圍的區(qū)域.微波能量通過(guò)環(huán)形天線饋入放電室以產(chǎn)生ECR等離子體[28].其中離子源通過(guò)雙柵極引出離子束形成推力,雙柵極由屏柵和加速柵組成,屏柵電位為正,加速柵電位為負(fù),柵極尺寸如表1所列.中和器在外界離子束流正電位作用下引出電子中和離子束,以保證推進(jìn)系統(tǒng)處于電平衡狀態(tài).中和器的電子引出板均布6個(gè)小孔,孔徑為2 mm,分布半徑為5 mm.

    圖1 2 cm ECRIT離子源放電室示意圖Fig.1.Schematic diagram of configuration for discharge chamber of 2 cm ECRIT ion source.

    表1 柵極結(jié)構(gòu)Table 1.Grid geometry.

    天線位置參數(shù)由L1和L2定義,其中L1為天線環(huán)形段上表面與屏柵下表面之間的距離,L2為內(nèi)磁環(huán)上表面與天線環(huán)形段下表面之間的距離,H1表示外磁環(huán)高度,W1表示外磁環(huán)寬度,H2表示內(nèi)磁環(huán)高度,W2表示內(nèi)磁環(huán)寬度.

    2.2 2 cm ECRIT磁路與天線位置對(duì)電子獲能的影響分析

    放電室內(nèi)電子獲得的微波能量主要發(fā)生在ECR區(qū),獲能公式[29]為

    其中 ?B/?s 為磁感應(yīng)強(qiáng)度B沿磁力線的梯度,E⊥為垂直于磁力線方向的微波電場(chǎng)強(qiáng)度,v//為電子平行磁場(chǎng)的速度,e 為電子電量.由此可見(jiàn)影響電子獲能的主要因素為靜磁場(chǎng)和微波電場(chǎng),為此設(shè)計(jì)四種不同的放電室磁路結(jié)構(gòu)如表2所列.它們與電子引出板可以構(gòu)成四種具有不同磁場(chǎng)分布規(guī)律的中和器,與屏柵和加速柵也構(gòu)成四種具有不同磁場(chǎng)分布規(guī)律的離子源.天線的位置是會(huì)影響微波電場(chǎng)分布,通過(guò)調(diào)節(jié)天線參數(shù)L1和L2來(lái)產(chǎn)生不同的微波電場(chǎng)分布規(guī)律.因而實(shí)驗(yàn)將用具有不同天線位置和磁路結(jié)構(gòu)的離子源和中和器開(kāi)展離子束流和耦合電壓的變化規(guī)律研究.

    表2 四種放電室的磁路結(jié)構(gòu)參數(shù)Table 2.The magnetic circuit structure parameters of four discharge chambers.

    當(dāng)L1= 2.5 mm與L2= 0.7 mm時(shí),采用COMSOL軟件計(jì)算四種放電室內(nèi)靜磁場(chǎng)和空間分布,計(jì)算結(jié)果如圖2所示.靜磁場(chǎng)和微波電場(chǎng)計(jì)算時(shí)假設(shè)放電室工作環(huán)境為空載條件,即為大氣環(huán)境,不存在等離子體.計(jì)算靜磁場(chǎng)時(shí),參數(shù)設(shè)置:永磁體(釤鈷合金)的相對(duì)磁導(dǎo)率為1.36;磁軛(軟鐵)的相對(duì)磁導(dǎo)率為4000;空氣的相對(duì)磁導(dǎo)率為1.靜磁場(chǎng)的求解域?yàn)槲锢砟P偷乃袇^(qū)域,邊界條件僅需設(shè)置外邊界為零磁標(biāo)量勢(shì)邊界.從圖2中可以看出1號(hào)放電室ECR區(qū)位于天線環(huán)形段下表面,2號(hào)放電室ECR區(qū)與天線環(huán)形段相交,3,4號(hào)放電室ECR區(qū)與天線環(huán)形段上表面的距離逐漸增大.

    針對(duì)4號(hào)放電室,固定L2= 0.7 mm,取L1分別為2.1,2.5,2.9 mm,采用COMSOL軟件計(jì)算出微波電場(chǎng)分布規(guī)律,如圖3所示.計(jì)算微波電場(chǎng)時(shí),參數(shù)設(shè)置:壁面(不銹鋼)的相對(duì)介電常數(shù)為1,電導(dǎo)率為4×106;天線(銅)的相對(duì)介電常數(shù)為1,電導(dǎo)率為6×107;空氣的相對(duì)介電常數(shù)為1,電導(dǎo)率為0 S/m.微波電場(chǎng)的求解域?yàn)榉烹娛覅^(qū)域,邊界條件僅需設(shè)置微波輸入端口邊界,輸入端口類型設(shè)置為同軸,輸入功率設(shè)為2 W.從圖3中可以發(fā)現(xiàn)L1的細(xì)微變化對(duì)微波電場(chǎng)強(qiáng)度E的影響較小,并根據(jù) (1) 式可知L1的增加對(duì)電子獲能影響較小.但是從磁力線分布可以發(fā)現(xiàn)隨著L1增加,更多磁力線與腔體壁面相交從而造成更多的粒子損失.

    基于COMSOL計(jì)算得到的靜磁場(chǎng)與微波電場(chǎng)分布,采用Matlab軟件編程計(jì)算不同磁路下電子獲能參數(shù).由于只有E⊥和 ?B/?s 影響電子獲能,所以這里定義電子獲能指標(biāo)為其分布如圖4所示.從圖4中可以看出獲能主要集中在天線的環(huán)形段附近.雖然1,2,3號(hào)放電室內(nèi)存在極小區(qū)域的獲能指標(biāo)最大值超過(guò)了5×106V2/(T·m),但4號(hào)放電室的電子獲能整體分布更均勻.雖然3號(hào)放電室ECR區(qū)更靠近天線環(huán)形段,其獲能指標(biāo)高于4號(hào)放電室.但由于高能電子沿磁力線遷移,根據(jù)ECR區(qū)與天線位置關(guān)系可估計(jì)損失大小為1號(hào)放電室 > 2號(hào)放電室 > 3號(hào)放電室 > 4號(hào)放電室.所以可推測(cè)4號(hào)放電室對(duì)電子加熱效果最佳.

    圖2 放電室內(nèi)磁場(chǎng)分布 (a) 1號(hào)放電室;(b) 2號(hào)放電室;(c) 3號(hào)放電室;(d) 4號(hào)放電室Fig.2.Distribution of magnetic flux density inside of the discharge chambers:(a) Chamber 1;(b) chamber 2;(c) chamber 3;(d) chamber 4.

    圖3 不同天線位置下4號(hào)放電室內(nèi)微波電場(chǎng)分布 (a) L1= 2.1 mm,L2= 0.7 mm;(b) L1= 2.5 mm,L2= 0.7 mm;(c) L1=2.9 mm,L2= 0.7 mmFig.3.Distribution of microwave electric field intensity inside of discharge chamber 4:(a) L1= 2.1 mm,L2= 0.7 mm;(b) L1=2.5 mm,L2= 0.7 mm;(c) L1= 2.9 mm,L2= 0.7 mm.

    圖4 不同放電室在ECR區(qū)的電子獲能指標(biāo)分布 (a) 1號(hào);(b) 2號(hào);(c) 3號(hào);(d) 4號(hào)Fig.4.Distribution of electron heating index in ECR layer of different discharge chambers:(a) Chamber 1;(b) chamber 2;(c) chamber 3;(d) chamber 4.

    2.3 實(shí)驗(yàn)系統(tǒng)

    推力器實(shí)驗(yàn)系統(tǒng)如圖5所示,由真空艙、柵極電源、中和器電源、氣體傳輸線路、微波源和微波傳輸線路組成.兩套微波源分別輸出4.2 GHz的微波能并通過(guò)環(huán)形器、同軸線纜、微波隔直器輸入到離子源和中和器,環(huán)形器反射波出口處的功率計(jì)用于測(cè)量反射微波功率.微波隔直器用于隔離直流電,氣路隔直器用于推力器與外界的電絕緣.離子源在中和器協(xié)同工作下進(jìn)行束流引出時(shí),中和器引出的電子電流應(yīng)該等于離子源引出的離子電流,離子羽流才被完全中和,此時(shí)施加在中和器殼體上的電勢(shì)為負(fù),地與中和器電勢(shì)之差為耦合電壓Up1.實(shí)驗(yàn)時(shí)屏柵電壓為1500 V,加速柵電壓為—350 V,引出的離子束流Ib為屏柵電流Is減去加速柵電流Ia.離子源束流單獨(dú)引出的實(shí)驗(yàn)連接線路同圖5,此時(shí)中和器不工作.但是離子源正前方5 cm處安裝接地的金屬板,金屬板在高速離子束作用下發(fā)射二次電子從而中和離子束流.中和器束流單獨(dú)引出實(shí)驗(yàn)的連接線路同圖5,此時(shí)離子源不工作.但中和器正前方5 mm處安放金屬板與中和器電源正極連接,中和器殼體接地,金屬板與中和器之間的電勢(shì)差為電子引出電壓Up2.實(shí)驗(yàn)時(shí),金屬板收集到的電流為引出的電子束流.所有實(shí)驗(yàn)都采用氙工質(zhì).實(shí)驗(yàn)中易存在的測(cè)量誤差以及減小誤差的措施如下.

    圖5 實(shí)驗(yàn)系統(tǒng)圖Fig.5.Schematic of the experimental system.

    1)微波功率的測(cè)量.由于調(diào)節(jié)微波增大或減小時(shí),個(gè)人習(xí)慣會(huì)導(dǎo)致讀數(shù)偏高或偏低,實(shí)驗(yàn)前需對(duì)微波源進(jìn)行多次標(biāo)定,取算術(shù)平均值來(lái)減小誤差.實(shí)驗(yàn)時(shí),由于測(cè)量束流引出時(shí),等離子體內(nèi)部存在波動(dòng)性,這會(huì)導(dǎo)致反射功率出現(xiàn)0.5%的波動(dòng),雖然波動(dòng)很小,但仍需要會(huì)記錄多組數(shù)據(jù),并取算術(shù)平均值來(lái)減小誤差.同時(shí),為了防止微波源長(zhǎng)時(shí)間工作后微波傳輸發(fā)生變化,第二天實(shí)驗(yàn)前仍會(huì)進(jìn)行重新標(biāo)定.

    2)引出離子束流的測(cè)量.同上可知,引出束流時(shí)內(nèi)部等離子體存在波動(dòng),引出束流的測(cè)量值會(huì)出現(xiàn)1%左右的波動(dòng),采取記錄多組數(shù)據(jù)后取算術(shù)平均值來(lái)減小誤差.

    3 實(shí)驗(yàn)結(jié)果與討論

    3.1 不同磁路和天線位置的離子源束流引出實(shí)驗(yàn)

    保持天線位置參數(shù)不變,取L1= 2.5 mm,L2=0.7 mm,采用同樣結(jié)構(gòu)的柵極結(jié)構(gòu),利用表2給出的四種放電室結(jié)構(gòu)形成四種離子源用于實(shí)驗(yàn)研究磁路結(jié)構(gòu)對(duì)離子源性能的影響.離子源獨(dú)立工作時(shí),入射功率Pi為2—3 W、氣體流量q 為0.1—0.3 sccm(1 sccm = 1 mL/min),不同磁路下離子源離子束流引出的實(shí)驗(yàn)結(jié)果如圖6所示.可知不同功率下,引出的離子束流與推進(jìn)劑流量呈正相關(guān),1號(hào)離子源的束流明顯低于其他三種離子源.功率2 W、流量增加時(shí),4號(hào)離子源束流逐漸大于其他三種離子源,2,3號(hào)離子源的離子束流逐漸大于1號(hào)離子源.功率3 W時(shí),流量增加,1號(hào)離子源的離子束流最低,3,4號(hào)離子源的離子束流都高于2號(hào)離子源,3號(hào)離子源的離子束流逐漸大于4號(hào)離子源.

    圖6 固定天線位置下四種離子源引出的離子束流Fig.6.Ion beam current of four ion sources at fixed antenna position.

    在2.2節(jié)中可知4號(hào)源的電子獲能較優(yōu),而且?guī)щ娏W訐p失最小,故4號(hào)源引出效果優(yōu)于1,2,3號(hào)源.但也存在意外情況:當(dāng)功率為3 W,流量大于0.2 sccm時(shí),3號(hào)源引出束流卻高于4號(hào)源.雖然與4號(hào)源相比,3號(hào)源ECR區(qū)整體電子獲能分布較差,帶電粒子損失較大,但3號(hào)源在天線環(huán)形段附近的ECR區(qū)電子獲能高于4號(hào)源,相同流量下且ECR區(qū)附近中性原子密度高于4號(hào)源.故功率為3 W,流量大于0.2 sccm時(shí),3號(hào)源電離產(chǎn)生的離子遠(yuǎn)多于損失的離子,從而其引出束流高于4號(hào)源.

    圖7 不同的天線位置下4號(hào)離子源的離子束流引出實(shí)驗(yàn)結(jié)果 (a) L1= 2.1 mm,L2= 0.7 mm;(b) L1= 2.5 mm,L2=0.7 mm;(c) L1= 2.9 mm,L2= 0.7 mmFig.7.Ion beam current of ion source 4 at different antenna position:(a) L1= 2.1 mm,L2= 0.7 mm;(b) L1=2.5 mm,L2= 0.7 mm;(c) L1= 2.9 mm,L2= 0.7 mm.

    針對(duì)功率為2 W時(shí),束流引出性能較好的4號(hào)離子源,固定L2= 0.7 mm,實(shí)驗(yàn)研究L1對(duì)引出束流的影響.圖7為其實(shí)驗(yàn)結(jié)果,可以看出,引出離子束流與流量、微波輸入功率成正相關(guān).在功率和流量變化范圍內(nèi),4號(hào)離子源的離子束流隨L1先增加后減小,L1= 2.5 mm時(shí)束流最高.其原因是天線位置參數(shù)L1代表了離子源內(nèi)部放電容積,L1過(guò)短,狹小的放電空間使高能電子與內(nèi)壁面碰撞損失速率遠(yuǎn)大于與中性原子碰撞電離速率,造成大部分能量損失,導(dǎo)致束流降低;另外,根據(jù)前述分析,L1增加,電子在遷移擴(kuò)散到柵極的過(guò)程中損失更多,從而導(dǎo)致引出的離子束流減小.文獻(xiàn)[26]的實(shí)驗(yàn)結(jié)果也證實(shí)了過(guò)大或過(guò)小的放電容積會(huì)影響束流引出.由此可見(jiàn),L1= 2.5 mm時(shí),引出的離子束流最大.

    綜上可知,功率、流量、磁路和天線位置都會(huì)影響離子源的束流引出.引出束流隨功率、流量的增加而增加.磁路是影響源內(nèi)電子獲能高低和等離子體分布的重要因素,而在相同功率、流量和磁路下,天線位置是影響帶電粒子的損失的重要因素.故改進(jìn)磁路和天線位置是提高離子源性能最直接的方式.

    3.2 中和器獨(dú)立工作時(shí)電子引出電壓的變化規(guī)律實(shí)驗(yàn)

    根據(jù)離子源實(shí)驗(yàn)情況,確定中和器給定的引出電子電流為4 mA.在此條件下以電子引出電壓Up2最低為目標(biāo),通過(guò)實(shí)驗(yàn)研究,確定出合理的放電室磁路結(jié)構(gòu)、天線位置以及工作參數(shù).確定此實(shí)驗(yàn)?zāi)繕?biāo)是因?yàn)橹泻推鳘?dú)立工作時(shí)Up2越低,則當(dāng)中和器與離子源協(xié)同工作時(shí),耦合電壓Up1越低,壁面鞘層電勢(shì)差越低,由此能緩解中和器壁面的離子濺射并提升其可靠性與壽命[28].

    保持天線位置不變,取L1= 3.5 mm,L2=0.7 mm,采用同樣結(jié)構(gòu)的電子引出孔板,利用表2給出的四種放電室結(jié)構(gòu)形成四種中和器用于實(shí)驗(yàn)研究磁路結(jié)構(gòu)對(duì)中和器Up2的影響,實(shí)驗(yàn)結(jié)果如圖8所示.可以看出,中和器的Up2都隨Ie的增加而增加.在相同流量下,功率越高,中和器引出相同的Ie所需Up2越低,反映出功率較高產(chǎn)生的等離子體密度高,中和器Up2低.在相同功率下,流量增加,中和器引出相同Ie所需Up2略微降低.文獻(xiàn)[27]的實(shí)驗(yàn)結(jié)果也證實(shí)了功率和流量與引出電壓呈負(fù)相關(guān).在微波輸入功率和流量的全部范圍內(nèi),2號(hào)中和器引出相同Ie的Up2均低于其他中和器.與1,2號(hào)中和器相比,3,4號(hào)中和器的ECR區(qū)距離電子引出板較近,當(dāng)?shù)入x子體向引出板擴(kuò)散時(shí),部分高能電子未發(fā)生電離碰撞便被引出,導(dǎo)致放電室內(nèi)等離子體密度較低;同時(shí)3,4號(hào)中和器距離進(jìn)氣孔較遠(yuǎn),在相同流量下ECR區(qū)附近的中性原子密度低,導(dǎo)致電離頻率和等離子體密度低,因而引出相同的Ie所需Up2高.而1號(hào)中和器獲能效果不如2號(hào)中和器,且粒子損失較大,故在相同功率、流量下,2號(hào)中和器內(nèi)能產(chǎn)生較高的等離子體密度,從而引出相同的電子束流所需的電子引出電壓低.

    圖8 電子引出電壓隨電子束流的變化 (a) 1號(hào)中和器;(b) 2號(hào)中和器;(c) 3號(hào)中和器;(d) 4號(hào)中和器Fig.8.Voltage for electron extraction vs.electron beam current from different neutralizers:(a) No.1;(b) No.2;(c) No.3;(d) No.4.

    根據(jù)圖8,在L2= 0.7 mm和微波輸入功率1 W條件下,選擇電子引出電壓Up2較低的2號(hào)中和器,進(jìn)行天線位置參數(shù)L1和流量工作參數(shù)對(duì)Up2影響的實(shí)驗(yàn)研究,實(shí)驗(yàn)結(jié)果如圖9所示.從圖9可以看出,L1相同時(shí)2號(hào)中和器引出相同Ie所需的Up2隨著流量增加而略微降低.原因是L1相同條件下,流量增加即放電室內(nèi)中性原子增多,從而電子與原子碰撞頻率增大,生成的等離子體密度增大,故引出相同Ie所需的Up2降低;但當(dāng)流量大于0.2 sccm時(shí),由于微波功率為1 W,ECR區(qū)加熱的高能電子與中性原子碰撞電離趨于飽和,故引出相同Ie所需的Up2變化不明顯.當(dāng)流量相同時(shí),隨著L1增大,引出相同Ie所需Up2降低.分析原因是L1減小時(shí),相同流量下,L1越小放電室內(nèi)氣體壓強(qiáng)越大,電子與原子之間發(fā)生的碰撞頻率更高,從而破壞了電子的ECR加熱過(guò)程,導(dǎo)致等離子體密度降低,引出的電子束流受限.但當(dāng)推進(jìn)劑流量為0.1 sccm和引出的電子束流為4 mA時(shí),Up2值分別為58,162.9和27.4 V,即隨L1先增加后減小.其原因是L1尺寸的增加會(huì)導(dǎo)致電子在壁面上的損失增加,但又有利于降低電子-原子碰撞頻率,使電子能獲得足夠的能量發(fā)生電離碰撞.故當(dāng)L1增加至4.3 mm時(shí),電子與原子的碰撞頻率合理,電子在ECR區(qū)能夠充分獲得微波能量,形成較高密度的等離子體從而降低電子引出電壓.

    圖9 不同L1尺寸下2號(hào)中和器在Pi= 1 W時(shí)的電子引出結(jié)果Fig.9.Electron extraction results of neutralizer 2 with different L1 at Pi= 1 W.

    綜上可知,功率、流量、磁路和天線位置都會(huì)影響中和器的電子引出電壓.引出壓隨功率、流量的增加而減小.磁路對(duì)中和器的影響與離子源不同,因?yàn)橹泻推鲀?nèi)部等離子體會(huì)受引出板電壓的影響,ECR區(qū)與引出板距離過(guò)近,會(huì)導(dǎo)致高能電子未電離就損耗,從而影響中和器的引出電壓.天線位置會(huì)影響帶電粒子的壁面損失和電離頻率,只有當(dāng)電子與原子的碰撞頻率和粒子損失都合理時(shí),中和器內(nèi)等離子體密度高,從而降低電子的引出電壓.

    3.3 離子源和中和器協(xié)同工作性能實(shí)驗(yàn)

    根據(jù)離子源和中和器的磁路和天線位置優(yōu)化結(jié)果,選擇L1= 2.5 mm,L2= 0.7 mm的4號(hào)離子源,L1= 4.3 mm,L2= 0.7 mm的2號(hào)中和器進(jìn)行中和實(shí)驗(yàn)研究.中和器工作參數(shù)固定為輸入功率1 W、流量0.1 sccm,圖10為離子源在2 W和3 W以及不同流量下完全中和時(shí)引出的束流實(shí)驗(yàn)結(jié)果,圖11為離子束完全中和時(shí)中和器耦合電壓實(shí)驗(yàn)結(jié)果.對(duì)比圖10和圖7 (b),可以發(fā)現(xiàn)離子束中和時(shí)的離子源引出束流有所提高,但低流量下差別不大.文獻(xiàn)[30]實(shí)驗(yàn)結(jié)果也證實(shí)了有無(wú)中和器工作對(duì)離子源引出束流影響甚微.圖11表明離子束流完全中和時(shí),中和器所需的耦合電壓Up1為12.0—17.4 V,比單獨(dú)中和器引出相同電流所需的引出電壓Up2低,這是因?yàn)檎妱?shì)的離子束使電子容易引出,從而降低了中和器所需的耦合電壓.

    圖10 中和時(shí)引出離子束流變化規(guī)律Fig.10.Ion beam current at the neutralization state.

    圖11 中和時(shí)中和器耦合電壓變化規(guī)律Fig.11.The coupling voltage of neutralizer at the neutralization state.

    3.4 推力器最優(yōu)結(jié)構(gòu)性能計(jì)算

    假設(shè)氙原子只進(jìn)行一次電離,2 cm ECRIT的推力T、放電損耗εb、推進(jìn)劑利用率ηm和比沖Isp由下式計(jì)算[28]:

    其中Ub,Pi,Pr,M,和g 分別是離子加速電壓、微波入射功率、反射功率、氙原子質(zhì)量、Xe氣質(zhì)量流量和重力加速度;為修正參數(shù),根據(jù)經(jīng)驗(yàn)取0.98(假離子束流發(fā)散半角為5度,含有5%的二價(jià)離子).根據(jù)3.3節(jié)實(shí)驗(yàn)結(jié)果,計(jì)算最優(yōu)結(jié)構(gòu)的中和器和離子源協(xié)同工作下的性能參數(shù).即,在中和器Pi=1 W,q = 0.1 sccm時(shí),離子源在Pi= 2 W,q =0.2 sccm條件下,Ib= 4.53 mA,計(jì)算得到T=315.3 μN(yùn),εb= 407.1 W/A,ηm= 31.63%,Isp=1639.1 s;離子源在Pi= 2 W,q = 0.3 sccm條件下,Ib= 5.3 mA,計(jì)算得到T= 368.6 μN(yùn),εb=337.5 W/A,ηm= 24.7%,Isp= 1277.6 s.

    4 結(jié) 論

    實(shí)驗(yàn)研究了磁路結(jié)構(gòu)以及天線位置對(duì)2 cm ECRIT離子源引出束流和中和器耦合電壓的影響,并獲得了結(jié)構(gòu)最優(yōu)的離子源和中和器,實(shí)驗(yàn)分析了中和器協(xié)同工作條件下的離子源引出束流特征,獲得了結(jié)構(gòu)最優(yōu)2 cm ECRIT的性能,結(jié)論如下.

    1) 離子源ECR區(qū)位置影響離子束流的引出.當(dāng)ECR區(qū)高于天線環(huán)形段上表面,電子獲能高并且?guī)щ娏W訐p失少,離子源能夠引出最大的離子束流.對(duì)于中和器,若ECR區(qū)低于天線環(huán)形段表面,高能電子在天線表面的損失過(guò)大,形成的等離子體密度低,從而使電子引出效果不佳;若ECR區(qū)高于天線環(huán)形段表面,ECR區(qū)附近的中性原子密度低以及部分高能電子未與中性原子發(fā)生電離碰撞便被引出,使得等離子體密度低從而使電子引出效果差.

    2) 固定的磁路結(jié)構(gòu)與天線參數(shù)L2= 0.7 mm時(shí),離子源存在一個(gè)最佳天線位置參數(shù)L1=2.5 mm,使得離子束流引出效果最優(yōu).L1過(guò)短會(huì)導(dǎo)致引出的離子束流減小,L1過(guò)長(zhǎng)會(huì)使得離子損失增加從而減小引出的離子束流.同樣,中和器也存在一個(gè)最佳參數(shù)L1= 4.3 mm使得電子束流引出效果更好.

    3) 隨功率增加,離子源引出束流增加、中和器電子引出電壓降低.在不同功率下,中和器以及離子源引出束流與流量呈正相關(guān).優(yōu)化的中和器與離子源協(xié)同工作實(shí)驗(yàn)表明離子源引出束流略高于無(wú)中和器情況,并且中和器引出相同的電子束流時(shí)所需的電壓降低.

    4) 中和器輸入功率和流量為1 W,0.1 sccm,離子源輸入功率和流量為2 W,0.3 sccm時(shí),推力器引出離子束流5.3 mA,推力和比沖分別為368.6 μN(yùn)和1277.6 s,推進(jìn)劑放電損耗和利用率分別為337.5 W/A和24.7%,中和器所需的耦合電壓為17.4 V.

    猜你喜歡
    離子源推力器磁路
    基于分布磁路法的非晶合金高速磁浮直線電機(jī)懸浮力計(jì)算
    單組元推力器倒置安裝多余物控制技術(shù)驗(yàn)證
    單電極介質(zhì)阻擋放電離子源的消電子技術(shù)研究
    基于等效磁路法的永磁同步電機(jī)特性分析
    用于小行星探測(cè)的離子推力器技術(shù)研究
    離子推力器和霍爾推力器的異同
    太空探索(2015年6期)2015-07-12 12:48:42
    基于ANSYS Workbench的微型揚(yáng)聲器磁路優(yōu)化分析
    Changes in feed intake,nutrient digestion,plasma metabolites,and oxidative stress parameters in dairy cows with subacute ruminal acidosis and its regulation with pelleted beet pulp
    高頻離子源形成的HD+離子束中H3+比例實(shí)測(cè)研究
    用于補(bǔ)償中子孔隙度測(cè)井的可控中子源技術(shù)
    亚洲,欧美精品.| 宅男免费午夜| 精品久久久久久久末码| 国产爱豆传媒在线观看| 免费看日本二区| 国产亚洲av嫩草精品影院| 黄片小视频在线播放| 久久精品91蜜桃| 国产精品1区2区在线观看.| 成人鲁丝片一二三区免费| 特级一级黄色大片| 国产视频内射| 久久久成人免费电影| 日韩人妻高清精品专区| 久久国产精品人妻蜜桃| 精品日产1卡2卡| 日本撒尿小便嘘嘘汇集6| 在现免费观看毛片| 性插视频无遮挡在线免费观看| 免费看美女性在线毛片视频| 欧美一区二区国产精品久久精品| 亚洲av免费在线观看| 一区福利在线观看| 亚洲欧美日韩无卡精品| 亚洲欧美激情综合另类| 99久久精品国产亚洲精品| 亚洲av美国av| 又黄又爽又免费观看的视频| 中文字幕av在线有码专区| 给我免费播放毛片高清在线观看| 99在线人妻在线中文字幕| 动漫黄色视频在线观看| 成人国产综合亚洲| 伊人久久精品亚洲午夜| 亚洲一区二区三区不卡视频| 好看av亚洲va欧美ⅴa在| 国产欧美日韩精品一区二区| 高清在线国产一区| 老司机午夜福利在线观看视频| 俺也久久电影网| 最近中文字幕高清免费大全6 | 欧美日韩综合久久久久久 | 淫秽高清视频在线观看| 一区福利在线观看| 日韩欧美 国产精品| 国产精品99久久久久久久久| 美女被艹到高潮喷水动态| 欧美另类亚洲清纯唯美| 能在线免费观看的黄片| 成人无遮挡网站| 欧美zozozo另类| АⅤ资源中文在线天堂| 婷婷色综合大香蕉| 日韩欧美国产一区二区入口| 国内精品美女久久久久久| 永久网站在线| 欧美日韩黄片免| 久久精品人妻少妇| 亚洲av成人不卡在线观看播放网| 欧美在线黄色| 少妇人妻一区二区三区视频| 99久久久亚洲精品蜜臀av| 搞女人的毛片| 日韩欧美国产在线观看| 蜜桃亚洲精品一区二区三区| 免费观看的影片在线观看| 夜夜爽天天搞| 亚洲精品色激情综合| 草草在线视频免费看| 精品一区二区免费观看| www日本黄色视频网| 五月伊人婷婷丁香| 久久国产乱子免费精品| 国产真实伦视频高清在线观看 | 成人永久免费在线观看视频| 欧洲精品卡2卡3卡4卡5卡区| 一进一出抽搐动态| 淫秽高清视频在线观看| 国产精品三级大全| 国产三级在线视频| 天天一区二区日本电影三级| 99国产极品粉嫩在线观看| 一区二区三区激情视频| 国产精品国产高清国产av| 亚洲色图av天堂| 又紧又爽又黄一区二区| 午夜福利在线观看免费完整高清在 | 人人妻,人人澡人人爽秒播| 午夜福利18| 男人舔女人下体高潮全视频| 国产极品精品免费视频能看的| 婷婷亚洲欧美| 精品久久久久久久人妻蜜臀av| 高清在线国产一区| 国产精品亚洲av一区麻豆| 日韩有码中文字幕| 免费观看的影片在线观看| 国产精品久久久久久人妻精品电影| 精品久久久久久久久av| 亚洲av免费在线观看| 日韩欧美一区二区三区在线观看| 精品久久久久久久久亚洲 | 国产一区二区三区视频了| av在线蜜桃| 人人妻,人人澡人人爽秒播| 丁香欧美五月| 一进一出好大好爽视频| 亚洲自拍偷在线| 成人三级黄色视频| 他把我摸到了高潮在线观看| 亚洲激情在线av| 色综合站精品国产| 综合色av麻豆| 丝袜美腿在线中文| 一级作爱视频免费观看| 欧美xxxx黑人xx丫x性爽| 亚洲av第一区精品v没综合| 欧美成人性av电影在线观看| 久久久久久久精品吃奶| 男女之事视频高清在线观看| 五月玫瑰六月丁香| 男女下面进入的视频免费午夜| 国产伦精品一区二区三区视频9| 久久这里只有精品中国| 麻豆一二三区av精品| 免费人成在线观看视频色| 亚洲国产精品久久男人天堂| 亚洲avbb在线观看| 国产成人欧美在线观看| 国产成人av教育| 老司机午夜福利在线观看视频| 一本精品99久久精品77| 可以在线观看毛片的网站| 欧美激情国产日韩精品一区| 欧美区成人在线视频| 国产精品99久久久久久久久| 99久久精品一区二区三区| 精品国内亚洲2022精品成人| 最近最新中文字幕大全电影3| 日韩大尺度精品在线看网址| 看免费av毛片| 精品久久久久久久人妻蜜臀av| 天天一区二区日本电影三级| 在线a可以看的网站| 在线观看66精品国产| .国产精品久久| 成人av一区二区三区在线看| 欧美又色又爽又黄视频| av国产免费在线观看| 久久精品影院6| 丰满人妻一区二区三区视频av| 亚洲国产精品成人综合色| 97超级碰碰碰精品色视频在线观看| 窝窝影院91人妻| 内射极品少妇av片p| 国产黄片美女视频| 日日摸夜夜添夜夜添小说| 99久久九九国产精品国产免费| 人妻丰满熟妇av一区二区三区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 可以在线观看的亚洲视频| 老司机午夜福利在线观看视频| 免费看日本二区| 国产成人啪精品午夜网站| 久久国产精品人妻蜜桃| 国产精品人妻久久久久久| 欧美不卡视频在线免费观看| 三级国产精品欧美在线观看| 亚洲欧美日韩无卡精品| 亚洲av美国av| 少妇裸体淫交视频免费看高清| 久久久久国内视频| 国产精品嫩草影院av在线观看 | 国产一区二区在线观看日韩| 成人永久免费在线观看视频| 91九色精品人成在线观看| 91狼人影院| 黄色日韩在线| 国产精品久久久久久久电影| 欧美三级亚洲精品| 两个人的视频大全免费| 国产爱豆传媒在线观看| 国产av在哪里看| 精华霜和精华液先用哪个| 男女下面进入的视频免费午夜| 亚洲精品色激情综合| 中文字幕熟女人妻在线| 久久久久久大精品| 夜夜爽天天搞| 亚洲人与动物交配视频| 蜜桃亚洲精品一区二区三区| 久久伊人香网站| 天堂网av新在线| 亚洲欧美日韩东京热| 亚洲经典国产精华液单 | 性插视频无遮挡在线免费观看| 国产精品嫩草影院av在线观看 | 十八禁网站免费在线| 婷婷色综合大香蕉| 亚洲电影在线观看av| 午夜福利高清视频| 国产av在哪里看| 亚洲成人久久性| 久久精品91蜜桃| 国产精品1区2区在线观看.| 日韩免费av在线播放| 亚洲国产精品999在线| www.色视频.com| 少妇的逼好多水| 69人妻影院| 麻豆国产av国片精品| 男人狂女人下面高潮的视频| 欧美一级a爱片免费观看看| 在线播放无遮挡| 免费观看精品视频网站| 国产高清三级在线| 波多野结衣巨乳人妻| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品久久久久久成人av| 99热这里只有是精品在线观看 | 中文字幕av成人在线电影| 村上凉子中文字幕在线| 亚洲精品色激情综合| 狠狠狠狠99中文字幕| 国产美女午夜福利| 精品人妻一区二区三区麻豆 | av在线蜜桃| 97人妻精品一区二区三区麻豆| 国产精品久久久久久亚洲av鲁大| 日韩中字成人| 十八禁网站免费在线| 变态另类丝袜制服| 噜噜噜噜噜久久久久久91| 免费大片18禁| 日本 欧美在线| 欧美日韩福利视频一区二区| 中出人妻视频一区二区| 丰满人妻熟妇乱又伦精品不卡| 欧美日本视频| 日本熟妇午夜| 在线观看免费视频日本深夜| 亚洲美女黄片视频| 丰满人妻一区二区三区视频av| 18禁在线播放成人免费| 色精品久久人妻99蜜桃| 精品熟女少妇八av免费久了| 亚洲第一欧美日韩一区二区三区| 窝窝影院91人妻| 成年女人看的毛片在线观看| 久久99热这里只有精品18| 亚洲av电影在线进入| 一边摸一边抽搐一进一小说| 一级av片app| 免费一级毛片在线播放高清视频| 成人永久免费在线观看视频| 精品福利观看| 成年女人看的毛片在线观看| 五月伊人婷婷丁香| 床上黄色一级片| 丝袜美腿在线中文| 在线看三级毛片| av在线老鸭窝| 丁香六月欧美| 黄色一级大片看看| 国产男靠女视频免费网站| 五月玫瑰六月丁香| 最近最新免费中文字幕在线| 最新在线观看一区二区三区| 国产av麻豆久久久久久久| 中文亚洲av片在线观看爽| 色噜噜av男人的天堂激情| 国产精品精品国产色婷婷| 午夜精品在线福利| 在线看三级毛片| av国产免费在线观看| 国产日本99.免费观看| 亚洲最大成人手机在线| 美女cb高潮喷水在线观看| 尤物成人国产欧美一区二区三区| 日韩精品中文字幕看吧| 亚洲国产精品999在线| av在线天堂中文字幕| 99精品久久久久人妻精品| 99热这里只有是精品50| 麻豆成人av在线观看| 少妇人妻一区二区三区视频| 国产91精品成人一区二区三区| 成年人黄色毛片网站| 啦啦啦观看免费观看视频高清| 一区二区三区激情视频| 麻豆久久精品国产亚洲av| 亚洲真实伦在线观看| 国产久久久一区二区三区| 别揉我奶头~嗯~啊~动态视频| 蜜桃久久精品国产亚洲av| 91狼人影院| 亚洲专区中文字幕在线| 丰满人妻一区二区三区视频av| 99国产精品一区二区蜜桃av| 国产精品亚洲av一区麻豆| 国产伦精品一区二区三区四那| 国产精品久久久久久亚洲av鲁大| 美女高潮喷水抽搐中文字幕| 亚洲aⅴ乱码一区二区在线播放| 桃红色精品国产亚洲av| 热99re8久久精品国产| 亚洲av中文字字幕乱码综合| 中出人妻视频一区二区| 韩国av一区二区三区四区| av在线观看视频网站免费| 丰满的人妻完整版| 久久99热这里只有精品18| 一级毛片久久久久久久久女| 国产日本99.免费观看| 在线播放国产精品三级| 精品福利观看| 亚洲美女视频黄频| 两个人视频免费观看高清| 麻豆av噜噜一区二区三区| 色综合婷婷激情| 国产乱人视频| 精品久久久久久久久av| 日韩亚洲欧美综合| 啦啦啦观看免费观看视频高清| 三级国产精品欧美在线观看| 亚州av有码| 午夜精品在线福利| 国产大屁股一区二区在线视频| 三级国产精品欧美在线观看| 日韩成人在线观看一区二区三区| 国产不卡一卡二| 美女高潮喷水抽搐中文字幕| 欧美一区二区亚洲| 给我免费播放毛片高清在线观看| 桃色一区二区三区在线观看| 国内精品久久久久精免费| 国产伦一二天堂av在线观看| 永久网站在线| 久久人人爽人人爽人人片va | 午夜福利视频1000在线观看| 国产色爽女视频免费观看| 成人特级黄色片久久久久久久| 亚洲,欧美精品.| 99在线视频只有这里精品首页| 欧美潮喷喷水| 少妇人妻一区二区三区视频| 午夜亚洲福利在线播放| 级片在线观看| 在线观看66精品国产| 亚洲,欧美精品.| 国产真实乱freesex| 赤兔流量卡办理| 在线观看一区二区三区| 一级作爱视频免费观看| 18+在线观看网站| 国产爱豆传媒在线观看| 女生性感内裤真人,穿戴方法视频| 麻豆国产av国片精品| 国产aⅴ精品一区二区三区波| 欧美日本视频| 人人妻,人人澡人人爽秒播| 久久精品综合一区二区三区| 简卡轻食公司| 日日摸夜夜添夜夜添小说| 国产高清视频在线播放一区| 久久精品影院6| 成人性生交大片免费视频hd| 99精品在免费线老司机午夜| 1000部很黄的大片| 久久久久久国产a免费观看| 久99久视频精品免费| 成人美女网站在线观看视频| 中文字幕av成人在线电影| 色吧在线观看| 成人性生交大片免费视频hd| 国产精品人妻久久久久久| 深夜a级毛片| 老鸭窝网址在线观看| 国产色爽女视频免费观看| 亚洲 国产 在线| 日本a在线网址| 两个人视频免费观看高清| 一进一出抽搐动态| 亚洲国产日韩欧美精品在线观看| 无遮挡黄片免费观看| 色播亚洲综合网| www.熟女人妻精品国产| 伦理电影大哥的女人| 黄色女人牲交| 成年人黄色毛片网站| 少妇的逼水好多| 日韩欧美精品v在线| 一边摸一边抽搐一进一小说| 亚洲第一电影网av| 两性午夜刺激爽爽歪歪视频在线观看| 91av网一区二区| 日韩欧美精品免费久久 | 久久国产乱子伦精品免费另类| 欧美激情在线99| 欧美另类亚洲清纯唯美| 看十八女毛片水多多多| 最新中文字幕久久久久| 国内毛片毛片毛片毛片毛片| 91麻豆av在线| 国产黄片美女视频| 日韩欧美三级三区| 国产人妻一区二区三区在| 看片在线看免费视频| 此物有八面人人有两片| 两个人视频免费观看高清| 久久久久精品国产欧美久久久| 麻豆久久精品国产亚洲av| 天美传媒精品一区二区| 99热这里只有是精品在线观看 | 久久久久久久久大av| 欧美成人性av电影在线观看| 99热这里只有精品一区| 欧美区成人在线视频| 九色成人免费人妻av| 国产v大片淫在线免费观看| 又紧又爽又黄一区二区| 国产精品三级大全| 中出人妻视频一区二区| 999久久久精品免费观看国产| 综合色av麻豆| 国产精品久久电影中文字幕| 两个人的视频大全免费| 日韩大尺度精品在线看网址| 一区二区三区高清视频在线| 最近中文字幕高清免费大全6 | 少妇裸体淫交视频免费看高清| 性欧美人与动物交配| 熟妇人妻久久中文字幕3abv| 亚洲国产欧洲综合997久久,| 给我免费播放毛片高清在线观看| 色综合婷婷激情| 亚洲美女黄片视频| 亚洲av一区综合| 丰满人妻熟妇乱又伦精品不卡| 国产老妇女一区| 国产精品自产拍在线观看55亚洲| 日本三级黄在线观看| 国产爱豆传媒在线观看| 亚洲av熟女| 高潮久久久久久久久久久不卡| 免费看a级黄色片| 夜夜爽天天搞| 成人精品一区二区免费| 国产不卡一卡二| 51国产日韩欧美| 免费无遮挡裸体视频| 深夜精品福利| 免费看光身美女| 欧美精品国产亚洲| 成年免费大片在线观看| 国产午夜福利久久久久久| 免费搜索国产男女视频| 韩国av一区二区三区四区| 九色成人免费人妻av| 99久久成人亚洲精品观看| 91久久精品电影网| www.www免费av| 欧美激情久久久久久爽电影| 制服丝袜大香蕉在线| 午夜日韩欧美国产| 亚洲精品色激情综合| 色av中文字幕| 波多野结衣高清无吗| 99riav亚洲国产免费| 91狼人影院| 久久久久久久久大av| 日本撒尿小便嘘嘘汇集6| 色哟哟哟哟哟哟| 岛国在线免费视频观看| av专区在线播放| 日韩中字成人| 国产麻豆成人av免费视频| 婷婷精品国产亚洲av| 中文字幕熟女人妻在线| 国产在线精品亚洲第一网站| 日韩成人在线观看一区二区三区| 毛片女人毛片| 成人国产综合亚洲| 午夜福利免费观看在线| 午夜久久久久精精品| 亚洲精品一区av在线观看| 一本一本综合久久| 成人国产综合亚洲| 欧美精品国产亚洲| 日本在线视频免费播放| 国产男靠女视频免费网站| 久久午夜亚洲精品久久| 又黄又爽又免费观看的视频| 成人特级av手机在线观看| 老女人水多毛片| 村上凉子中文字幕在线| 国产欧美日韩一区二区精品| 日韩精品中文字幕看吧| 久久久久久大精品| 国内精品久久久久精免费| 九色国产91popny在线| 美女大奶头视频| 桃色一区二区三区在线观看| 动漫黄色视频在线观看| 欧美性猛交╳xxx乱大交人| 亚洲人成网站高清观看| 男人舔女人下体高潮全视频| 最新在线观看一区二区三区| 国产免费一级a男人的天堂| 欧美在线一区亚洲| 亚洲成av人片在线播放无| 美女高潮喷水抽搐中文字幕| 亚洲av电影不卡..在线观看| 欧美日本亚洲视频在线播放| 久久久国产成人免费| 国产乱人伦免费视频| 少妇裸体淫交视频免费看高清| 免费人成在线观看视频色| 久久婷婷人人爽人人干人人爱| 国产日本99.免费观看| 成人毛片a级毛片在线播放| 亚洲精品久久国产高清桃花| 亚洲va日本ⅴa欧美va伊人久久| 亚洲最大成人av| 啦啦啦观看免费观看视频高清| 国产伦精品一区二区三区视频9| 首页视频小说图片口味搜索| 大型黄色视频在线免费观看| 国产精华一区二区三区| 波野结衣二区三区在线| 狠狠狠狠99中文字幕| 国产精品不卡视频一区二区 | 一级作爱视频免费观看| 国产极品精品免费视频能看的| 最近最新免费中文字幕在线| 国产探花极品一区二区| 嫩草影院精品99| 精品久久久久久久人妻蜜臀av| 亚洲18禁久久av| 日韩欧美国产在线观看| 搡女人真爽免费视频火全软件 | 小蜜桃在线观看免费完整版高清| 国产主播在线观看一区二区| 无人区码免费观看不卡| 不卡一级毛片| 亚洲黑人精品在线| 国产亚洲精品综合一区在线观看| 亚州av有码| 免费黄网站久久成人精品 | 午夜精品在线福利| 麻豆成人av在线观看| 亚洲精华国产精华精| 黄色配什么色好看| 精品免费久久久久久久清纯| 色5月婷婷丁香| 欧美性猛交黑人性爽| 国产欧美日韩精品一区二区| 日本精品一区二区三区蜜桃| 最近最新中文字幕大全电影3| 日韩欧美一区二区三区在线观看| 两人在一起打扑克的视频| 国产精品久久久久久亚洲av鲁大| 成人性生交大片免费视频hd| 1000部很黄的大片| 亚洲最大成人手机在线| 日本黄色视频三级网站网址| 最近最新免费中文字幕在线| av福利片在线观看| 18禁裸乳无遮挡免费网站照片| 婷婷精品国产亚洲av| 免费搜索国产男女视频| 婷婷亚洲欧美| 啪啪无遮挡十八禁网站| 小说图片视频综合网站| 欧美xxxx黑人xx丫x性爽| 国产精品98久久久久久宅男小说| 亚洲国产精品成人综合色| 757午夜福利合集在线观看| 午夜亚洲福利在线播放| 亚洲,欧美,日韩| 制服丝袜大香蕉在线| 国产免费男女视频| 一夜夜www| 757午夜福利合集在线观看| 91在线精品国自产拍蜜月| 国产免费男女视频| 嫩草影院新地址| 99国产精品一区二区三区| 性色av乱码一区二区三区2| 欧美绝顶高潮抽搐喷水| 噜噜噜噜噜久久久久久91| 欧美日本视频| 俄罗斯特黄特色一大片| 韩国av一区二区三区四区| av在线天堂中文字幕| 国产成+人综合+亚洲专区| 在线观看美女被高潮喷水网站 | 国产高潮美女av| 精品人妻熟女av久视频| 久久热精品热| 亚洲人成伊人成综合网2020| 国产av在哪里看| 久久草成人影院| av欧美777| 十八禁国产超污无遮挡网站| 亚洲熟妇中文字幕五十中出| 日韩大尺度精品在线看网址| 色综合站精品国产| netflix在线观看网站| 偷拍熟女少妇极品色| 韩国av一区二区三区四区| 日本免费a在线| 99久久九九国产精品国产免费| 嫩草影院入口| 亚洲第一区二区三区不卡|