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

    無中和器射頻離子推力器原理研究

    2016-02-13 07:43:37夏廣慶徐宗琦王鵬朱雨陳茂林
    中國空間科學(xué)技術(shù) 2016年1期
    關(guān)鍵詞:電負(fù)性推力器柵極

    夏廣慶,徐宗琦,王鵬,朱雨,陳茂林

    1.大連理工大學(xué) 工業(yè)裝備結(jié)構(gòu)分析國家重點(diǎn)實(shí)驗(yàn)室,大連116024 2.中國科學(xué)院深圳先進(jìn)技術(shù)研究院,深圳518055 3.西北工業(yè)大學(xué) 燃燒、熱結(jié)構(gòu)與內(nèi)流場重點(diǎn)實(shí)驗(yàn)室,西安710072

    無中和器射頻離子推力器原理研究

    夏廣慶1,*,徐宗琦1,王鵬1,朱雨2,陳茂林3

    1.大連理工大學(xué) 工業(yè)裝備結(jié)構(gòu)分析國家重點(diǎn)實(shí)驗(yàn)室,大連116024 2.中國科學(xué)院深圳先進(jìn)技術(shù)研究院,深圳518055 3.西北工業(yè)大學(xué) 燃燒、熱結(jié)構(gòu)與內(nèi)流場重點(diǎn)實(shí)驗(yàn)室,西安710072

    在電推力器的廣泛應(yīng)用中,大部分都采用加速正離子的方式產(chǎn)生推力,且需要安裝中和器發(fā)射電子對噴射出的正離子進(jìn)行中和,否則會導(dǎo)致航天器自充電,對其通信及電子器件造成損害。因此,中和器的性能成為制約電推力器工作狀態(tài)和壽命的重要因素。為了克服該缺點(diǎn),介紹了一種基于同時加速正、負(fù)離子的無中和器射頻離子推力器,闡述了其結(jié)構(gòu)組成及推進(jìn)原理,分析了離子-離子的產(chǎn)生、正負(fù)離子的加速過程,指出了關(guān)鍵技術(shù)包括電負(fù)性工質(zhì)氣體放電特性、磁場過濾電子束縛效能,以及可周期性交替加速正、負(fù)離子的柵極偏置電壓加載方式。分析表明,該推力器在低軌航天器及深空探測器中具有潛在的應(yīng)用前景。

    電推進(jìn);負(fù)離子;射頻離子推力器;電負(fù)性氣體;離子-離子推力器;磁過濾;交變電壓

    隨著電推進(jìn)技術(shù)的發(fā)展,電推力器已經(jīng)越來越廣泛地應(yīng)用于不同種類的衛(wèi)星上。1998年以前,空間科學(xué)探測衛(wèi)星上電推力器的應(yīng)用率約為3%,相應(yīng)的電推進(jìn)裝置有200種左右[1]。如今,約有20%的衛(wèi)星使用電推進(jìn)系統(tǒng),美國、俄羅斯及歐洲等國已將電推力器廣泛用于衛(wèi)星和空間飛行器,中國也已經(jīng)在實(shí)踐9號衛(wèi)星上成功驗(yàn)證和應(yīng)用霍爾推力器及離子推力器,未來將發(fā)展更大規(guī)模的電推進(jìn)應(yīng)用平臺。

    通常認(rèn)為,按不同離子加速方式,電推力器主要分為三種:電熱式、靜電式和電磁式。

    電熱式推力器主要有電阻加熱推力器和電弧加熱推力器,電阻加熱推力器是利用電能加熱工質(zhì),使其汽化分解,經(jīng)噴管膨脹,加速噴出產(chǎn)生推力,工作方式與化學(xué)推進(jìn)基本相同。電弧加熱推力器是靠兩個電極之間加電壓后形成的電弧加熱工質(zhì),經(jīng)擴(kuò)張噴管高速噴出產(chǎn)生推力。電熱式推力器的特點(diǎn)是能夠獲得較大推力,但是比沖相對較低。

    靜電式推力器除了傳統(tǒng)的考夫曼離子推力器外,還有微波離子推力器和射頻離子推力器,其主要是利用微波或射頻的能量使易于電離的工質(zhì)電離,形成電子和正離子,正離子在靜電場的作用下加速噴出產(chǎn)生推力。由于需要柵極系統(tǒng)提供靜電場,因此也可稱柵極離子推力器。靜電式推力器的特點(diǎn)是比沖大、推力效率高。

    電磁式推力器主要有霍爾推力器和磁等離子體推力器,其主要是利用電能使工質(zhì)形成等離子體,在洛侖茲力的作用下加速噴出產(chǎn)生推力。電磁式推力器的特點(diǎn)是離子束集中、推力較大。

    柵極離子推力器和霍爾推力器的共同點(diǎn)是都能夠獲得較高比沖[2-3],且為了使噴射到空間中的羽流呈電中性,都必須安裝中和器發(fā)出電子與正離子中和。因此,中和器的穩(wěn)定性、使用壽命及功率損耗成為制約衛(wèi)星總體工作性能的重要因素。另外,中和器占據(jù)推力器結(jié)構(gòu)空間,增大推力器總體質(zhì)量,從而增加發(fā)射成本。因此,人們普遍認(rèn)為中和器是電推力器的一個軟肋[3-4],同時也在考慮研制一種基于可以分別加速正離子和負(fù)離子并最終使其中和為中性粒子的無中和器射頻離子推力器。

    1 推力器研究現(xiàn)狀

    關(guān)于無中和器射頻離子推力器的理論研究及相關(guān)試驗(yàn)工作已在國外開展,美國喬治亞理工學(xué)院設(shè)計出一種名為馬歇爾離子-離子推力器(Marshall’s Ion-io N Thruster,MINT),以電負(fù)性氣體作為工質(zhì)氣體,交替加速正負(fù)離子產(chǎn)生推力,并通過試驗(yàn),獲得了理想的結(jié)果[5]。法國國家科學(xué)研究中心也根據(jù)無需中和器這一目標(biāo)研究出一種電負(fù)性離子推力器,經(jīng)過試驗(yàn)發(fā)現(xiàn),電負(fù)性工質(zhì)氣體具有較高的電離度,可以獲得較高的等離子體密度。當(dāng)柵極加速電壓為500 V時,可以達(dá)到與傳統(tǒng)柵極推力器相近的推力及比沖[6]。美國海軍研究實(shí)驗(yàn)室通過試驗(yàn)證實(shí)了在柵極板上加載正弦偏置電壓能夠從離子-離子等離子體中引出正負(fù)離子[7]。美國休斯頓大學(xué)研究學(xué)者建立了時變流體模型用以研究射頻電壓與直流電壓對離子-離子等離子體動力學(xué)特征的影響并發(fā)現(xiàn)其鞘層結(jié)構(gòu)與傳統(tǒng)的柵極離子推力器相比具有很大不同,但是對鞘層的空間尺度及預(yù)鞘層是否存在影響尚未進(jìn)行詳細(xì)的研究[8-9]。

    電負(fù)性等離子體是一種含有正離子、負(fù)離子及電子的等離子體,其電負(fù)性程度可以通過負(fù)離子密度與電子密度的比值α衡量,即α= n-/ne。電子溫度與負(fù)離子溫度的比值γ也是描述電負(fù)性等離子體性質(zhì)的重要參數(shù)之一,即γ=Te/T-。由于負(fù)離子比正離子質(zhì)量大,因此電負(fù)性等離子體有很多特征參數(shù)不同于電正性等離子體,比如玻姆速度、鞘層厚度及懸浮電位等[10-11]。

    有研究表明,無中和器射頻離子推力器能量傳輸效率高達(dá)90%,推力可達(dá)20 m N/k W,比沖在30 000 m/s左右。但是,為實(shí)現(xiàn)激發(fā)后產(chǎn)生電負(fù)性等離子體,該射頻離子推力器需要使用電負(fù)性氣體作為工質(zhì)氣體。

    2 推力器工作性能及原理

    2.1 推力器的推力和比沖

    對于航天推力器,其推力可表示為

    式中:vex為離子噴出速度;dm/dt為推力器工質(zhì)總質(zhì)量的變化率,推力大小可以根據(jù)任務(wù)需要控制。

    對式(1)進(jìn)行積分,得

    式中:Δv為航天器速度變化增量;m0和mf分別為航天器初始質(zhì)量和最終質(zhì)量。因此,Δm= m0-mf為推力器工質(zhì)質(zhì)量變化量。式(2)為火箭理想速度公式,又稱為齊奧爾科夫斯基公式。

    比沖可表示為

    式中:g0為重力加速度常數(shù)。從式(2)和式(3)中可以看出,比沖影響推進(jìn)工質(zhì)消耗:速度增量相同情況下,比沖越大,推進(jìn)工質(zhì)消耗越少。電推力器具有較高比沖,因此可以使用較少推進(jìn)工質(zhì),節(jié)約發(fā)射成本。

    無中和器射頻離子推力器與傳統(tǒng)的柵極離子推力器推力和比沖產(chǎn)生方式相同。柵極離子推力器主要依靠引出由電離室產(chǎn)生的離子進(jìn)行靜電加速來工作,通過在多個柵極之間施加電偏壓實(shí)現(xiàn)對離子的加速,并產(chǎn)生推力。對于離子束,離子噴出速度vex與離子運(yùn)動速度vb相等,離子總質(zhì)量隨時間的變化率dm/dt可表示為離子流量Γi。

    根據(jù)動能定理,忽略離子能量損失,則

    式中:e為基本電荷,U0為柵極之間的電勢差,Mi為離子質(zhì)量。于是,推力可表示為

    式中:Ai為柵極有效面積,忽略等離子體發(fā)散及能量損失;ns為邊界處等離子體密度;Te為等離子體溫度。當(dāng)推力達(dá)到最大值時,根據(jù)離子電流密度的連續(xù)性方程和能量方程,求解泊松方程可得電極最大電流[12]式中:ε0為真空介電常數(shù);d為柵極板間距離。假設(shè)等離子體密度均勻,當(dāng)?shù)入x子體流達(dá)到最大值時,有Γi=JCL/e,則最大推力

    從式(7)中可以看出,柵極離子推力器的最大推力與柵極有效面積、柵極間距離及加速電壓有關(guān),柵極有效面積越大,柵極板間距離越小,加速電壓越大,推力器產(chǎn)生的推力越大。同時,也可以看出離子質(zhì)量對推力沒有影響,但會影響推力器的比沖及工質(zhì)氣體的損耗。

    2.2 推力器的工作原理

    無中和器射頻離子推力器是一種根據(jù)電負(fù)性等離子體推進(jìn)(Plasma propulsion with Electronegative GASES,PEGASES)原理探索研究的新型靜電式電推力器,故也稱為電負(fù)性氣體等離子體推力器。它的主要特點(diǎn)是基于射頻耦合激勵的方式,使用電負(fù)性氣體和常規(guī)電推進(jìn)氣體的混合工質(zhì)放電,在推力器放電室內(nèi)同時產(chǎn)生正離子、負(fù)離子和電子,通過徑向磁場的磁障屏蔽效應(yīng)過濾電子,使得正、負(fù)離子均可到達(dá)柵極加速系統(tǒng)。在柵極系統(tǒng)交變加速電壓的作用下,正離子與負(fù)離子周期性噴出,進(jìn)而產(chǎn)生推力,最后二者復(fù)合為中性粒子,保證羽流宏觀上呈電中性,因此不需要中和器,從而可以大大減輕結(jié)構(gòu)質(zhì)量、提高壽命和穩(wěn)定性,并且可減少帶電離子的回流和對航天器的濺射作用。圖1所示為推力器工作過程簡圖。

    圖1 無中和器射頻離子推力器工作過程簡圖Fig.1 Working process of the RIT without neutralizer

    無中和器射頻離子推力器工作過程主要分為三步:1)等離子體的產(chǎn)生;2)離子-離子的形成;3)正負(fù)離子的加速。

    (1)等離子體的產(chǎn)生

    等離子體主要是由射頻功率源通過匹配器將能量耦合給天線,天線激發(fā)電離中性工質(zhì)氣體而產(chǎn)生放電,有關(guān)射頻放電理論在此不詳細(xì)闡述[5]。該過程關(guān)鍵在于射頻功率源的能量能否最大程度地傳遞給中性工質(zhì)氣體,即產(chǎn)生高密度的等離子體。通過電磁場的作用將能量最終耦合給等離子體主要有3種方式:容性耦合(CCP)、感應(yīng)耦合(ICP)和波耦合。在容性耦合模式中,通過隨射頻電壓變化而改變的振蕩鞘層使電子獲得能量并激發(fā)放電。在感應(yīng)耦合模式中,通過射頻天線的電流在等離子體中產(chǎn)生感應(yīng)電場,電子從電場中獲得能量并激發(fā)放電,電場向等離子體內(nèi)部傳播一段距離直至場強(qiáng)衰減為零,這段距離就是所謂的趨膚深度。在波耦合模式中,外加磁場使電磁波通過朗道阻尼作用將能量傳給電子,然后激發(fā)放電。

    無中和器射頻離子推力器采用的等離子源是感應(yīng)耦合等離子體源,射頻頻率為4 MHz,射頻功率為200 W,兩個柵極板放在等離子體出口處,柵極板由不銹鋼制成且厚度為0.8 mm,離子通過率為60%,柵極孔直徑為2 mm,柵極板間距為2 mm。磁過濾裝置提供徑向變化范圍為0~0.034 T的磁場,天線纏繞在石英放電室外并集成封裝,這樣能夠減少功率損失,增大ICP耦合效率[13-14],圖2為其整體結(jié)構(gòu)簡圖。首先使用傳輸線將射頻功率源與匹配器連接,然后能量通過匹配器傳輸?shù)缴漕l天線,射頻天線通過感應(yīng)耦合作用加熱放電室中等離子體使其溫度不斷升高,直至發(fā)生電離。ICP放電實(shí)際上削弱了CCP放電,從而可以忽略等離子體射頻電勢能。在這個系統(tǒng)中,功率效率能夠達(dá)到90%以上,因此可以有效地應(yīng)用于空間推進(jìn)系統(tǒng)。

    圖2 無中和器射頻離子推力器整體結(jié)構(gòu)簡圖Fig.2 Structure of the RIT without neutralizer

    電負(fù)性工質(zhì)氣體轉(zhuǎn)變?yōu)樨?fù)離子通常需要兩個過程:電離和復(fù)合。中性氣體分子首先被電離為正離子和電子,然后電子與未被電離的中性粒子復(fù)合成負(fù)離子,這個過程必須使用電負(fù)性氣體作為工質(zhì)氣體才能實(shí)現(xiàn),且在等離子體內(nèi)部及低逸出功金屬表面均能夠發(fā)生。例如,負(fù)氫離子能夠在金屬銫表面產(chǎn)生[15-16]。

    在電負(fù)性氣體試驗(yàn)中發(fā)現(xiàn),負(fù)離子能夠通過電子的離解俘獲碰撞產(chǎn)生。在這個過程中,電子能量過高而不能直接吸附于中性粒子,因此電子與中性粒子碰撞后會使其電離,或者通過自身的振動及旋轉(zhuǎn)運(yùn)動消耗能量,然后與中性粒子復(fù)合形成負(fù)離子。對于大部分電負(fù)性氣體,高溫電子與中性氣體分子碰撞使其電離為正離子,低溫電子與中性氣體結(jié)合為負(fù)離子。負(fù)離子的密度與電子密度有關(guān),被束縛的電子越多,負(fù)離子密度越高,電負(fù)性程度越大,因此使用磁場過濾裝置使電子被束縛在放電室中可有效地產(chǎn)生負(fù)離子。

    (2)離子-離子的形成

    這個工作過程的主要目的是將等離子體分成兩個區(qū)域:1)等離子體中心區(qū)域;2)下游離子-離子等離子體區(qū)域。在等離子體中心區(qū)域,電子溫度很高,熱運(yùn)動劇烈,與中性氣體分子碰撞后會使其有效地電離。在下游離子-離子等離子體區(qū)域,電子密度幾乎可以忽略,等離子體的能量主要指的是離子能量。綜上所述,為了能使這兩個區(qū)域形成,需要控制等離子體中電子的溫度。通常情況下,在低溫等離子體中,電子溫度Te與電離平衡 (電子產(chǎn)生和消失)、工質(zhì)氣體種類及氣體壓強(qiáng)p與離子有效擴(kuò)散特征尺寸L乘積有關(guān)[11]。電子加熱的3種射頻耦合方式如前文所述,放電功率的值對電子溫度影響較小。在負(fù)離子源中,電子通過安裝在離子引出孔前邊的磁過濾裝置 (局部橫向磁場)后被冷卻,能量降低而更易與中性粒子結(jié)合。這項(xiàng)技術(shù)被應(yīng)用于PEGASES推力器,其中的磁過濾裝置由磁場方向垂直于放電室軸向的永磁體構(gòu)成,電子的冷卻程度與磁感應(yīng)強(qiáng)度和方向有關(guān)。試驗(yàn)中發(fā)現(xiàn),在磁感應(yīng)強(qiáng)度最大處電子溫度最低[14],同時,在磁感應(yīng)強(qiáng)度最大處附近存在離子-離子等離子體[17],且離子數(shù)密度是電子數(shù)密度的三倍左右。在局部區(qū)域中,離子密度保持很高,這說明離子在穿過磁場過程中不受磁場影響。RF功率僅為120 W時,離子-離子等離子體區(qū)域中的離子密度就能達(dá)到3×1017m-3。

    (3)正負(fù)離子的加速

    無中和器射頻離子推力器屬于柵極離子推力器的一種,其離子加速機(jī)制與傳統(tǒng)的柵極離子推力器相同,但是其不同點(diǎn)在于使用同一個柵極交替加速正負(fù)離子。為了實(shí)現(xiàn)正負(fù)離子的交替加速,需要在第一個柵極板上施加正弦偏置電壓,第二個柵極板接地[18]。為了引出并加速正負(fù)離子,需要使推力器壁面懸浮而非接地,并且兩個柵極板之間的電場方向在一個偏置周期內(nèi)改變。柵極系統(tǒng)中,正負(fù)離子運(yùn)動軌跡如圖3所示。

    在任何離子源中,加速柵極的設(shè)計必須與加速電壓及等離子體參數(shù)(密度和溫度)相匹配,這樣可以獲得性質(zhì)較好的等離子體羽流,即離子束具有較為合適的正負(fù)離子流量和擴(kuò)散角,從而能夠滿足推力要求,其推力可以通過式(5)表示。對于單純的離子-離子等離子體,進(jìn)入鞘層的離子具有相同的玻姆速度,即vB=此時離子溫度可用電子溫度表示。推力器的推力效率為

    圖3 正負(fù)離子運(yùn)動軌跡簡圖Fig.3 Motion trajectories of the positive and negative ions

    式中:QT為總質(zhì)量流率;Pin為輸入功率。

    3 關(guān)鍵技術(shù)

    相比于其他柵極離子推力器,無中和器射頻離子推力器的關(guān)鍵技術(shù)主要在于以下幾個方面:電負(fù)性工質(zhì)氣體的選擇、徑向磁場過濾裝置的作用及柵極偏置電壓的加載方式。下面就這幾個方面進(jìn)行簡要分析。

    3.1 電負(fù)性工質(zhì)氣體的選擇

    無中和器射頻離子推力器最主要的特點(diǎn)在于必須使用電負(fù)性氣體作為工質(zhì)氣體。電負(fù)性氣體是具有極易與自由電子結(jié)合形成負(fù)離子并保持穩(wěn)定狀態(tài)能力的氣體。電負(fù)性氣體是由非金屬元素組成的單質(zhì)或化合物,非金屬性較強(qiáng)。其通常是由鹵族元素組成的單質(zhì)分子氣體,比如氯氣、碘蒸汽,或者含有鹵族元素的化合物分子氣體,比如SF6和CF4。另外,電負(fù)性氣體還包括O2和H2等非鹵族元素氣體。因此,選擇電負(fù)性工質(zhì)氣體時應(yīng)該首先考慮上述氣體,或者幾種氣體的混合物。然而,在這些電負(fù)性氣體中,I2被看作是最有發(fā)展前景的工質(zhì),這是因?yàn)镮2是一種電負(fù)性極強(qiáng)的物質(zhì),且在衛(wèi)星發(fā)射到太空過程中極易以固態(tài)形式儲存,很大程度上節(jié)省了推進(jìn)劑所占空間[19]。

    分子的電離和復(fù)合率系數(shù)通常受電子溫度的影響,反應(yīng)過程為[20]

    由式(9)、式(10)和式(11)可以看出,中性氣體分子首先與高溫電子發(fā)生電離碰撞產(chǎn)生正離子和電子,低溫電子與未被電離的中性氣體分子結(jié)合成負(fù)離子,此時在放電室上游同時存在正、負(fù)離子及電子,為了維持電離-復(fù)合反應(yīng)的不斷進(jìn)行,必須通過磁場約束電子使其留在放電室中,從而提高負(fù)離子的密度,即增大電負(fù)性程度。

    產(chǎn)生電負(fù)性等離子體的關(guān)鍵在于控制電子能量分布,即高溫電子與低溫電子的數(shù)量,從而使發(fā)生電離與復(fù)合的中性氣體分子具有合適的比例。

    國外已有試驗(yàn)表明,在壓強(qiáng)為0.13 Pa、射頻功率為100 W的條件下,采用SF6作為工質(zhì)氣體產(chǎn)生正、負(fù)離子的速率相當(dāng),并具有較高的復(fù)合率[20]。

    3.2 徑向磁場過濾裝置的作用

    在傳統(tǒng)推力器內(nèi)部的低溫等離子體中,電子的作用是引發(fā)中性氣體放電及中和噴出的正離子,然而在無中和器射頻離子推力器中,電子不僅被用于引發(fā)中性氣體放電,而且用于與未被電離的中性粒子復(fù)合為負(fù)離子[21]。在多數(shù)情況下,離子-離子等離子體更適合用于低溫等離子體設(shè)備,離子-離子可以通過脈沖放電或者磁場過濾裝置形成[22-23]。在脈沖放電中,余輝里的電子能量較低,可以與中性離子復(fù)合為負(fù)離子。當(dāng)使用磁場過濾裝置時,負(fù)離子的形成區(qū)域能夠被限定在一定空間范圍內(nèi)。通過磁場的作用,電子被約束在上游電離區(qū)域,下游只有可與中性粒子復(fù)合的少量低溫電子。

    徑向磁場過濾裝置對離子-離子的形成具有很大影響,徑向磁場可以由永磁體或電磁鐵提供,其作用主要體現(xiàn)在以下3個方面:

    1)約束等離子體源中的電子,確保工質(zhì)氣體能夠有效電離并維持放電。

    2)減少擴(kuò)散到放電室下游離子-離子區(qū)的電子數(shù)量,使其沿著磁力線運(yùn)動并最終附著于器壁。

    3)降低通過磁場的電子的溫度,使其更易與氣體分子結(jié)合形成負(fù)離子。

    雖然這些理論已經(jīng)建立,但是電子冷卻機(jī)理、場強(qiáng)梯度分布方式等問題還沒有完全分析清楚,因此采用徑向磁場過濾裝置產(chǎn)生離子-離子等離子體的方法還需深入研究,關(guān)鍵是如何提高加速柵極上游的磁場對電子的束縛能力,即阻止電子被加速噴出的能力。

    3.3 柵極偏置電壓的加載方式

    無中和器射頻離子推力器柵極板上通常加載正弦偏置電壓,電場方向隨時間交替變化,以使正負(fù)離子能夠分別被加速。但是在一個周期中,被加速的正負(fù)離子的總量不一定相同,因此需要根據(jù)正負(fù)離子的加速效率,對二者的加速時間做出適當(dāng)調(diào)整,以保證噴射出的正負(fù)離子總量相同,進(jìn)而能夠完全復(fù)合為中性粒子,使羽流呈電中性。通過進(jìn)一步研究發(fā)現(xiàn),方向交替變化的方波電壓也能實(shí)現(xiàn)正負(fù)離子加速。通過粒子網(wǎng)格法(Particle In Cell,PIC)仿真發(fā)現(xiàn),從等離子體源中引出的正負(fù)離子流量比例由電負(fù)性程度決定,負(fù)離子密度遠(yuǎn)高于電子密度[24]。當(dāng)負(fù)離子從等離子體引出時,也能檢測到同時被引出的部分電子。在這種情況下,負(fù)電荷加速周期應(yīng)該比相應(yīng)的正電荷加速周期短,以此來補(bǔ)償電荷密度和電荷流量,電荷的加速周期可以根據(jù)需要做適當(dāng)?shù)膬?yōu)化。這種周期的不對稱性可以推廣到僅加速正電荷和電子的系統(tǒng)中。對一個等離子體加速系統(tǒng)進(jìn)行優(yōu)化時,正負(fù)離子的加速電壓大小及周期可以根據(jù)正負(fù)離子溫度的不同作適當(dāng)調(diào)整,以保證加速后的正負(fù)離子電荷量相等,使羽流最終呈電中性。國外已有研究學(xué)者采用一維PIC方法對于在柵極板上加載不同波形偏置電壓時的鞘層結(jié)構(gòu)進(jìn)行了模擬[25],相同的問題也被采用聯(lián)立離子連續(xù)性方程、動量方程及泊松方程求解的方法研究[26],結(jié)果表明鞘層結(jié)構(gòu)與偏置電壓頻率及離子與柵極板碰撞頻率密切相關(guān),離子運(yùn)動速度遠(yuǎn)高于熱運(yùn)動的平均速度,因此說明當(dāng)偏置電壓模式不同時,形成的鞘層都能對離子有明顯的作用。

    柵極板上加載電壓的頻率范圍應(yīng)考慮下列因素:1)離子-離子等離子體自身回旋頻率和離子通過柵極板所需時間;2)單束離子電荷在下游空間產(chǎn)生的電勢大小。這種分析模型建立在一維基礎(chǔ)上,完全忽略二次電子散射及其他外部的空間電荷中和源,最終推測出加載電壓頻率應(yīng)該在幾兆赫茲區(qū)域內(nèi)[27]。

    另有研究發(fā)現(xiàn),利用射頻自偏壓原理也能實(shí)現(xiàn)柵極加速等離子體。這種方法是根據(jù)離子和電子對時變、電場的響應(yīng)時間不同,通過形成于電極前不斷變化的空間電荷鞘層調(diào)節(jié)電壓,并最終使施加在電極板上的射頻偏置電壓與等離子體相匹配。實(shí)現(xiàn)這種加速方法完全不需要磁場及中和器,僅需要控制等離子體的流動參數(shù)。在施加于柵極間的射頻振蕩電場的作用下,離子能夠連續(xù)不斷地被加速,而電子卻周期性地被引出。離子能量由場強(qiáng)的均值決定。電子的加速過程具有自洽性,因此不能通過一般方法控制電子能量的空間分布[28]。已經(jīng)證實(shí),擴(kuò)散到下游空間區(qū)域內(nèi)的等離子體流呈電中性。由于這種加速方法的特殊性,邊界條件對整個等離子體流都會產(chǎn)生影響,引起其他粒子輸運(yùn)過程產(chǎn)生,如等離子體的非雙極擴(kuò)散[29]。射頻自偏壓加速等離子體的關(guān)鍵在于粒子流能否視為非均質(zhì)等離子體,即離子與電子的流量相等,或者視為傳統(tǒng)的等離子體束流,即正離子流被低溫均質(zhì)電子云包圍[30-33]。因此,可以通過試驗(yàn)測量并比較離子流與電子流在下游射頻自偏壓系統(tǒng)中角分率的大小來判斷離子與電子流量是否相等。

    法國巴黎綜合理工學(xué)院根據(jù)射頻自偏壓原理設(shè)計出一種名為“海王星”的柵極加速等離子體系統(tǒng),并在試驗(yàn)中使用馬赫探針測量等離子體流中離子與電子數(shù)。經(jīng)過試驗(yàn)發(fā)現(xiàn),在此系統(tǒng)中電子與離子流量幾乎相等,等離子體流速在10 000~35 000 m/s之間,并隨著射頻電壓的增大而增大,射頻電壓約為375 V,等離子體流中心處的離子流量約為1016s-1·m-2[34],且等離子體流具有很強(qiáng)的方向性,因而能夠準(zhǔn)確判斷出推力器羽流方向[35]。然而,采用馬赫探針方法測量粒子流量具有一定誤差,不能準(zhǔn)確說明離子數(shù)與電子數(shù)完全相等,因此還需進(jìn)行更深一步地理論分析與試驗(yàn)驗(yàn)證。

    4 結(jié)束語

    與傳統(tǒng)的柵極離子推力器相比,無中和器射頻離子推力器具有很大的發(fā)展空間,其特點(diǎn)在于使用同一個等離子體源分別加速正負(fù)離子,這樣可以不需要外加中和器就能實(shí)現(xiàn)等離子體羽流最終呈電中性,因此具有優(yōu)越性能。雖然將這種推力器最終應(yīng)用到空間飛行器上還需經(jīng)過深入的理論研究和試驗(yàn)分析,但是已經(jīng)證實(shí)在磁過濾裝置的下游存在離子-離子等離子體,且中性粒子電離程度已經(jīng)滿足柵極離子推力器粒子電離程度的要求。離子-離子等離子體源能夠產(chǎn)生足夠高密度的正負(fù)離子,從而無中和器射頻離子推力器能夠在500 V加速電壓下工作,具有與傳統(tǒng)柵極推力器相近的推力和比沖。

    總之,無中和器射頻離子推力器具有潛在的應(yīng)用前景。其結(jié)構(gòu)組成相對簡單,無需中和器,因此可以避免粒子對傳統(tǒng)中和器熱陰極材料污染的制約,從而有效將電推力器的應(yīng)用從地球同步軌道拓展至低地球軌道,且具有較好的推力及比沖性能,能夠滿足未來空間科學(xué)探測任務(wù)的需求。

    References)

    [1] SANCHEZ M M,POLLARD J E.Spacecraft electric propulsion-an overview[J].Journal of Propulsion and Power,1998,14(5):688-699.

    [2] KAUFAMAN H R.Technology of closed drift thrusters AIAA-1985-0023[R].Reston:AIAA,1985.

    [3] GOEBEL D M,KATZ I.Fundamentals of electric propulsion:ion and Hall thrusters[M].New Jersey: Wiley,2008:89-90.

    [4] NISHIYAMA K,KUNINAKA H.Discussion on performance history and operations of hayabusa ion engines[J].Transactions of the Japanese Society for Artificial Intelligence,Aerospace Technology Japan, 2012,10:1-8.

    [5] NATALIE R S,THOMAS M L,MITCHELL L R W, et al.Design and preliminary testing plan of electronegative ion thruster,AIAA-2014-3425[R].Reston:AIAA,2014. [6] AANESLAND A,BREDIN J,POPELIER L,et al. Direction for the future-successive acceleration of positive and negative ions applied to space propulsion[R].CERN Yellow Report,2013:575-584.

    [7] WALTON S G,LEONARDT D,FEMSLER R F,et al.Extraction of positive and negative ions from electronbeam generated plasmas[J].Applied Physics Letters, 2002,81(6):987-989.

    [8] MIDHA V,ECONOMOU D J.Dynamics of ion-ion plasmas under radio frequency bias[J].Journal of Applied Physics,2001,90(3):1102-1114.

    [9] MIDHA V,RAMAMURTHI B,ECONOMOU D J. Time evolution of an ion-ion plasma after the application of a direct current bias[J].Journal of Applied Physics, 2002,91(10):6383-6387.

    [10] CHABERT P,BRAITHWAITE N S J.Physics of radio frequency plasma[M].Cambridge:Cambridge University,2011.

    [11] LIEBERMAN M A,LICHTENBERG A J.Principles of plasma discharges and materials processing[M]. New York:Wiley-Interscience,2005.

    [12] HARPER J M,CUOMO J J,KAUFMAN H R. Technology and applications of broad-beam ion sources used in sputtering[J].Journal of Vacuum Science and Technology,1982,21(3):737.

    [13] GODYAK V A.Electrical and plasma parameters of ICP with high coupling efficiency[J].Plasma Sources Science Technology,2011,20:025004.

    [14] AANESLAND A,BREDIN J,CHABERT P,et al. Electron energy distribution function and plasma parameters across magnetic filters[J].Applied Physics Letters,2012,100(4):044102.

    [15] BELCHENKO Y I,DIMOV G I,DUDNIKOV V G.A powerful injector of neutrals with a surface plasma source of negative ions[J].Nuclear Fusion, 1974,14(1):113.

    [16] BACAL M.Physics aspects of negative ion sources[J]. Nuclear Fusion,2006,46(6):250-259.

    [17] BREDIN J,CHABERT P,AANESLAND A.Langmiur probe analysis of highly electronegative plasma[J]. Applied Physics Letters,2013,102(15):154107.

    [18] AANESLAND A,CHABERT P,IRZYK M,et al. Method and device for forming a plasma beam:US 2013/0300288 A1[P].2013-11-14.

    [19] SZABO J,ROBIN M,PAINTAL S,et al.Iodine propellant space propulsion[C]∥33rd International Electric Propulsion Conference,Washington,2013.

    [20] RAPP D,BRIGLIA D D.Total cross sections for ionization and attachment in gases by electron impact.Ⅱ.Negative-ion formation[J].Journal of Chemical Physics,1965,43:1480-1486.

    [21] GOEBEL D M,KATZ I.Fundamentals of electric propulsion:ion and hall thrusters[M].New Jersey: Wiley,2008:202-204.

    [22] POPELIER L,AANESLAND A,CHABERT P.Response of an ion-ion plasma to dc biased electrodes[J].Journal of Physics D:Applied Physics,2011,44(31):315203.

    [23] HEMSWORTH R S,INOUE T.Positive and negative ion sources for magnetic fusion[J].IEEE Transactions on Plasma Science,2005,33(6):1799-1813.

    [24] OUDINI N,RAIMBAULT J L,CHABERT P,et al.Particle-in-cell simulation of an electronegative plasma under direct current bias studied in a large range of electronegativity[J].Physics Plasma,2013, 20(4):043501.

    [25] MEIGE A,LERAY G,RAIMBAULT J,et al. Sheath and presheath in ion-ion plasmas via particle in cell simulation[J].Applied Physics Letters,2008, 92(6):061501.

    [26] MIDHA V,ECONOMOU J.Dynamics of ion-ion plasmas under radio frequency bias[J].Journal of Applied Physics,90(3):1102.

    [27] AANESLAND A.Low temperature plasma exposed to external magnetic and electric fields[D].Habilitation a Diriger des Recherches:Universite Pierre et Marie Curie,2013.

    [28] RAFALSKYI D,AANESLAND A.Coincident ion acceleration and electron extraction for space propulsion using the self-bias formed on a set of RF biased grids bounding a plasma source[J].Journal of Physics D: Applied Physics,2014,47(49):495203.

    [29] LAFLEUR T,AANESLAND A.Ambipolar and nonambipolar diffusion in an rf plasma source containing a magnetic filter[J].Physics of Plasmas,2014,21(6): 063510.

    [30] SOLOSHENKO I A.Physics of ion beam plasma and problems of intensive ion beam transportation [J].Review of Scientific Instruments,1996,67(4): 1646-1652.

    [31] DUDIN S V,ZYKOV A V,FARENIK V I.Low energy ion beam space charge neutralization[J].Review of Scientific Instruments,1994,65(4):1451.

    [32] RAFALSKYI D,AANESLAND A.Hysteresis effects in the formation of a neutralizing beam plasma at low ion energy[J].Europhysics Letter,2013,104(3): 35004.

    [33] HUMPHRIES S JR.Intense ion-beam neutralization in free space[J].Applied Physics Letters,1978,32(12): 792.

    [34] RAFALSKYI D,AANESLAND A.Plasma acceleration using a radio frequency self-bias effect[J].Physics of Plasmas,2015,22(6):063502.

    [35] RAFALSKYI D,AANESLAND A.Plasma propulsion using novel concepts[C]∥32nd International Conference on Phenomena in Ionized Gases,Romania,July 26-31, 2015.

    (編輯:車曉玲)

    Research on the principle of RF ion thruster without a neutralizer

    XIA Guangqing1,*,XU Zongqi1,WANG Peng1,ZHU Yu2,CHEN Maolin3
    1.State Key Laboratory of Structural Analysis for Industrial Equipment,Dalian University of Technology,Dalian 116024,China 2.Shenzhen Institutes of Advanced Technology,Chinese Academy of Sciences,Shenzhen 518055,China 3.Science and Technology on Combustion,Internal Flow and Thermo-Structure Laboratory,Northwestern Polytechnical University,Xi'an 710072,China

    In the wide application of electric thrusters,most of them accelerate positive ions to generate thrust,and need to install the neutralizer which can emit electrons to neutralize the positive ions.Otherwise it will cause the spacecraft self-charging,and the communication equipment and electronic devices will be damaged.Therefore,the neutralizer performance becomes an important factor restricting the working state of the electric thrusters.In order to overcome the shortcoming,a kind of RF ion thruster without a neutralizer based on accel-erating positive and negative ions alternately was introduced.The structural components and the propulsion principles were described.The generation of ion-ions and accelerating process of positive and negative ions were discussed.The key techniques including electronegative gas discharge characteristics and magnetic field filter the electrons efficiency,as well as the applied bias voltage way for grids in which the positive and negative ions could be accelerated alternately and periodically were pointed.The analysis results indicate this thruster has potential application prospect in the low orbit spacecraft and deep space probes.

    electric propulsion;negative ion;RF ion thruster(RIT);electronegative gas; ion-ion thruster;magnetic filter;alternating voltage

    V439+.1

    :A

    10.3780/j.issn.1000-758X.2016.0004

    2015-11-26;

    :2015-12-30;錄用日期:2016-01-18;< class="emphasis_bold">網(wǎng)絡(luò)出版時間

    時間:2016-02-24 13:26:25

    http:∥www.cnki.net/kcms/detail/11.1859.V.20160224.1326.003.html

    裝備預(yù)研共用技術(shù)基金;中國博士后科學(xué)基金(2013M541230);中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金[DUT15ZD(G)01]

    *通迅作者:夏廣慶(1979-),男,副教授,gq.xia@dlut.edu.cn,主要研究方向?yàn)殡娡七M(jìn)

    夏廣慶,徐宗琦,王鵬,等.無中和器射頻離子推力器原理研究[J].中國空間科學(xué)技術(shù),2016,36(1):1-8.

    XIANG G Q,XU Z Q,WANG P,et al.Research on the principle of RF ion thruster without a neutralizer[J].Chinese Space Science and Technology,2016,36(1):1-8(in Chinese).

    http:∥zgkj.cast.cn

    猜你喜歡
    電負(fù)性推力器柵極
    元素電負(fù)性的常見考查
    單組元推力器倒置安裝多余物控制技術(shù)驗(yàn)證
    離子推力器三柵極組件熱形變仿真分析及試驗(yàn)研究
    真空與低溫(2022年2期)2022-03-30 07:11:22
    元素不同氧化態(tài)電負(fù)性的一種簡單標(biāo)度方法
    安徽化工(2020年6期)2020-12-21 09:49:04
    關(guān)于電負(fù)性的一些思考
    柵極液壓成型專用設(shè)備的研制
    超原子理論計算基團(tuán)電負(fù)性的研究
    用于小行星探測的離子推力器技術(shù)研究
    IGBT柵極驅(qū)動電阻的選擇
    離子推力器和霍爾推力器的異同
    太空探索(2015年6期)2015-07-12 12:48:42
    99九九线精品视频在线观看视频| 91在线精品国自产拍蜜月| 欧美日韩精品成人综合77777| 免费观看在线日韩| 我的老师免费观看完整版| 日本wwww免费看| 国产黄频视频在线观看| 亚洲电影在线观看av| 国产男女超爽视频在线观看| 久久久国产精品麻豆| 女人久久www免费人成看片| 午夜福利,免费看| 欧美高清成人免费视频www| 国产av一区二区精品久久| 久久婷婷青草| 欧美精品高潮呻吟av久久| 久久久a久久爽久久v久久| 99热这里只有是精品50| 中文字幕久久专区| 十八禁高潮呻吟视频 | av有码第一页| 国产精品一区二区三区四区免费观看| 岛国毛片在线播放| 性高湖久久久久久久久免费观看| 人人妻人人澡人人爽人人夜夜| 国产精品一区二区在线观看99| 三级经典国产精品| 老司机影院毛片| 色网站视频免费| 天堂中文最新版在线下载| 我的老师免费观看完整版| 精品人妻偷拍中文字幕| 久久人妻熟女aⅴ| 国内少妇人妻偷人精品xxx网站| 亚洲av日韩在线播放| 制服丝袜香蕉在线| 国产免费福利视频在线观看| 中文天堂在线官网| 国内少妇人妻偷人精品xxx网站| 国产极品天堂在线| 成人黄色视频免费在线看| 最新中文字幕久久久久| 草草在线视频免费看| 99久久中文字幕三级久久日本| 18禁裸乳无遮挡动漫免费视频| 夫妻午夜视频| 久久影院123| 天美传媒精品一区二区| 热re99久久国产66热| 黑人猛操日本美女一级片| 七月丁香在线播放| 少妇丰满av| 亚洲av成人精品一二三区| xxx大片免费视频| 欧美少妇被猛烈插入视频| 人人妻人人添人人爽欧美一区卜| 国产成人a∨麻豆精品| a 毛片基地| 晚上一个人看的免费电影| 丁香六月天网| 两个人免费观看高清视频 | 九色成人免费人妻av| 国产综合精华液| 国产成人a∨麻豆精品| 亚洲国产精品成人久久小说| 天堂俺去俺来也www色官网| 久热久热在线精品观看| 中文字幕免费在线视频6| 久久国产精品男人的天堂亚洲 | 18禁在线播放成人免费| 黄色毛片三级朝国网站 | 欧美日韩综合久久久久久| 在线看a的网站| 深夜a级毛片| 日韩中文字幕视频在线看片| 成人18禁高潮啪啪吃奶动态图 | 日韩av免费高清视频| 偷拍熟女少妇极品色| 久久午夜福利片| 国产精品嫩草影院av在线观看| 日韩一区二区三区影片| 国产日韩欧美亚洲二区| 80岁老熟妇乱子伦牲交| 欧美日本中文国产一区发布| 成年美女黄网站色视频大全免费 | 午夜av观看不卡| 亚洲av.av天堂| 丰满乱子伦码专区| 在线观看免费高清a一片| 久久久久网色| 亚洲精品456在线播放app| 亚洲av成人精品一区久久| a级一级毛片免费在线观看| 五月开心婷婷网| 交换朋友夫妻互换小说| av福利片在线| 日韩av免费高清视频| 久久影院123| 99久国产av精品国产电影| 99热6这里只有精品| 伊人久久国产一区二区| 人妻 亚洲 视频| 国产伦理片在线播放av一区| 欧美区成人在线视频| 自线自在国产av| 午夜精品国产一区二区电影| 欧美97在线视频| 亚洲欧美日韩东京热| 在线看a的网站| 日韩中字成人| 色哟哟·www| 在线观看免费高清a一片| 韩国av在线不卡| 又爽又黄a免费视频| 欧美高清成人免费视频www| 黑丝袜美女国产一区| 男女边摸边吃奶| 色婷婷av一区二区三区视频| 99九九线精品视频在线观看视频| 午夜91福利影院| 麻豆乱淫一区二区| 中文在线观看免费www的网站| a级片在线免费高清观看视频| 黄色日韩在线| 亚洲精品亚洲一区二区| 色婷婷av一区二区三区视频| 偷拍熟女少妇极品色| av不卡在线播放| 女的被弄到高潮叫床怎么办| 另类亚洲欧美激情| 欧美日韩综合久久久久久| 国产一区二区在线观看日韩| 2018国产大陆天天弄谢| 精品少妇内射三级| 亚洲精品第二区| 老熟女久久久| 久久av网站| 国内少妇人妻偷人精品xxx网站| 人人妻人人澡人人看| 一级片'在线观看视频| 久久97久久精品| 国产精品一区www在线观看| 久久精品国产亚洲av天美| kizo精华| av有码第一页| 男女边摸边吃奶| 各种免费的搞黄视频| 久久国产亚洲av麻豆专区| 少妇熟女欧美另类| 中文在线观看免费www的网站| 一个人看视频在线观看www免费| 欧美精品人与动牲交sv欧美| 久久久久久伊人网av| 国产精品女同一区二区软件| 美女cb高潮喷水在线观看| 欧美 日韩 精品 国产| 蜜桃在线观看..| 99精国产麻豆久久婷婷| 一级,二级,三级黄色视频| 我的女老师完整版在线观看| 亚洲,一卡二卡三卡| 免费在线观看成人毛片| 成年女人在线观看亚洲视频| 美女脱内裤让男人舔精品视频| 精品久久久久久久久亚洲| 伦理电影免费视频| 老司机影院成人| 老女人水多毛片| 伦理电影大哥的女人| 最新的欧美精品一区二区| 国产欧美另类精品又又久久亚洲欧美| 色视频在线一区二区三区| 午夜激情福利司机影院| 黄色欧美视频在线观看| 国产av精品麻豆| 五月伊人婷婷丁香| 久久人人爽av亚洲精品天堂| 人妻少妇偷人精品九色| 国产日韩欧美亚洲二区| 久久久久久久久久久免费av| 在线观看免费日韩欧美大片 | 在线观看三级黄色| av国产精品久久久久影院| 少妇熟女欧美另类| 国产一级毛片在线| 日日啪夜夜爽| 免费人成在线观看视频色| 中文资源天堂在线| 国产永久视频网站| 一边亲一边摸免费视频| 大片免费播放器 马上看| 成人二区视频| 国产在线男女| 日韩欧美 国产精品| 日日撸夜夜添| 91久久精品国产一区二区成人| 国产精品国产三级国产专区5o| 精品视频人人做人人爽| av.在线天堂| 亚洲欧洲国产日韩| 色视频在线一区二区三区| 麻豆乱淫一区二区| 人人妻人人看人人澡| 日韩人妻高清精品专区| 2018国产大陆天天弄谢| 国产亚洲精品久久久com| 亚洲不卡免费看| 九草在线视频观看| 免费高清在线观看视频在线观看| 亚洲自偷自拍三级| 一级毛片我不卡| 亚洲av在线观看美女高潮| 岛国毛片在线播放| 国产一区二区三区av在线| 亚洲av电影在线观看一区二区三区| 简卡轻食公司| 亚洲国产毛片av蜜桃av| 高清午夜精品一区二区三区| 少妇人妻久久综合中文| 热re99久久国产66热| 自线自在国产av| h日本视频在线播放| 成人国产麻豆网| av国产精品久久久久影院| 国产精品一区www在线观看| 精品一品国产午夜福利视频| 三级国产精品欧美在线观看| 丰满人妻一区二区三区视频av| 美女cb高潮喷水在线观看| 精品亚洲成国产av| av福利片在线观看| 亚洲美女视频黄频| 大片免费播放器 马上看| 免费黄网站久久成人精品| 国产在线视频一区二区| 观看av在线不卡| 91久久精品国产一区二区成人| 久久久久久人妻| 国产成人精品福利久久| 一区在线观看完整版| 日本av手机在线免费观看| 一级av片app| 国产探花极品一区二区| 少妇人妻一区二区三区视频| 久久97久久精品| 国产欧美日韩一区二区三区在线 | 女性生殖器流出的白浆| 日韩中字成人| 男男h啪啪无遮挡| 久久久久久人妻| 欧美最新免费一区二区三区| 亚洲怡红院男人天堂| av卡一久久| 在线观看人妻少妇| 国产成人精品久久久久久| xxx大片免费视频| 久久久久精品性色| 欧美变态另类bdsm刘玥| 久久免费观看电影| 黑人猛操日本美女一级片| 久久人人爽人人片av| 国产av一区二区精品久久| 久久99精品国语久久久| 一级二级三级毛片免费看| 婷婷色av中文字幕| 色5月婷婷丁香| 欧美亚洲 丝袜 人妻 在线| 中国国产av一级| 超碰97精品在线观看| 免费大片黄手机在线观看| 成人国产麻豆网| 能在线免费看毛片的网站| 国产免费又黄又爽又色| 成人毛片60女人毛片免费| 激情五月婷婷亚洲| 亚洲精品日本国产第一区| 欧美变态另类bdsm刘玥| 一区二区三区精品91| 久久久久久久精品精品| 久热久热在线精品观看| 美女脱内裤让男人舔精品视频| 桃花免费在线播放| 日本欧美国产在线视频| 久热这里只有精品99| 国产熟女欧美一区二区| 亚洲欧美成人综合另类久久久| 精品国产国语对白av| 一个人免费看片子| 纵有疾风起免费观看全集完整版| 丰满乱子伦码专区| 超碰97精品在线观看| 久久久久久久久久久免费av| 国产精品国产三级国产av玫瑰| 亚洲人与动物交配视频| 欧美激情国产日韩精品一区| 欧美最新免费一区二区三区| 亚洲一级一片aⅴ在线观看| 国产成人aa在线观看| 一个人看视频在线观看www免费| www.色视频.com| 国产精品99久久99久久久不卡 | 久久久久久久国产电影| 高清视频免费观看一区二区| 日本wwww免费看| 亚洲精品第二区| av天堂中文字幕网| 多毛熟女@视频| 久久久久久久大尺度免费视频| 熟女电影av网| 久久影院123| 亚洲国产精品国产精品| 国产永久视频网站| 少妇猛男粗大的猛烈进出视频| 久久久a久久爽久久v久久| 国产精品福利在线免费观看| 春色校园在线视频观看| a级一级毛片免费在线观看| 国产极品天堂在线| av卡一久久| 精品一品国产午夜福利视频| 蜜臀久久99精品久久宅男| 亚洲高清免费不卡视频| 亚洲av成人精品一区久久| 日日啪夜夜撸| 精品卡一卡二卡四卡免费| 免费看不卡的av| 亚洲av男天堂| 最近的中文字幕免费完整| 国产一区二区三区综合在线观看 | 免费看不卡的av| 免费高清在线观看视频在线观看| 日本av手机在线免费观看| 极品教师在线视频| 99视频精品全部免费 在线| 一级毛片电影观看| 国产精品久久久久久久久免| 久久国产精品大桥未久av | 国产精品麻豆人妻色哟哟久久| 91久久精品电影网| 国产日韩一区二区三区精品不卡 | 嫩草影院入口| 中文字幕久久专区| 午夜av观看不卡| 性色avwww在线观看| 成人午夜精彩视频在线观看| 亚洲熟女精品中文字幕| 性色avwww在线观看| 男女无遮挡免费网站观看| 国产精品成人在线| av卡一久久| 男女无遮挡免费网站观看| 国产精品成人在线| 国产免费一区二区三区四区乱码| 国产日韩欧美在线精品| 日韩av在线免费看完整版不卡| 校园人妻丝袜中文字幕| 97在线人人人人妻| 久久久久精品性色| 国产在线男女| 最新中文字幕久久久久| 亚洲欧美中文字幕日韩二区| 伦精品一区二区三区| 国产 一区精品| 国内精品宾馆在线| 久久久精品免费免费高清| 一区在线观看完整版| 美女福利国产在线| 国产精品麻豆人妻色哟哟久久| 自拍偷自拍亚洲精品老妇| 亚洲高清免费不卡视频| 久久这里有精品视频免费| 深夜a级毛片| 久久国内精品自在自线图片| 少妇被粗大的猛进出69影院 | 我要看黄色一级片免费的| 国产午夜精品一二区理论片| 3wmmmm亚洲av在线观看| 99久国产av精品国产电影| 日韩成人伦理影院| 亚洲综合色惰| 久久99热6这里只有精品| 成人亚洲精品一区在线观看| 国产成人精品婷婷| 性色av一级| 欧美日韩一区二区视频在线观看视频在线| 一级毛片电影观看| 国产91av在线免费观看| 亚洲精品aⅴ在线观看| 嫩草影院新地址| 成人综合一区亚洲| 日本午夜av视频| 国产淫片久久久久久久久| 国产一区二区三区综合在线观看 | 成人午夜精彩视频在线观看| 欧美日韩精品成人综合77777| 高清欧美精品videossex| 日日摸夜夜添夜夜添av毛片| 99久久精品热视频| 久久精品夜色国产| 99视频精品全部免费 在线| 久久久久人妻精品一区果冻| 精品亚洲乱码少妇综合久久| 在线观看av片永久免费下载| 中文字幕av电影在线播放| 免费人成在线观看视频色| 国产欧美日韩综合在线一区二区 | 成人午夜精彩视频在线观看| 亚洲av成人精品一二三区| 亚洲人与动物交配视频| 欧美最新免费一区二区三区| av天堂久久9| 国产亚洲av片在线观看秒播厂| 国产精品成人在线| 国内精品宾馆在线| 国产精品.久久久| 亚洲欧洲日产国产| 亚洲精品乱码久久久v下载方式| 国产永久视频网站| videossex国产| 久久亚洲国产成人精品v| 亚州av有码| 三上悠亚av全集在线观看 | 国产淫片久久久久久久久| 免费大片黄手机在线观看| 亚洲内射少妇av| 久久久久久人妻| 国产极品天堂在线| 久久国内精品自在自线图片| 在线精品无人区一区二区三| 久久人人爽人人爽人人片va| 内地一区二区视频在线| 在线观看av片永久免费下载| 亚洲高清免费不卡视频| 麻豆成人av视频| 久久久国产一区二区| 亚洲精品一区蜜桃| 深夜a级毛片| 亚洲欧美一区二区三区国产| 高清午夜精品一区二区三区| 免费观看性生交大片5| 赤兔流量卡办理| 欧美精品国产亚洲| 蜜桃久久精品国产亚洲av| 只有这里有精品99| 大香蕉97超碰在线| 三级经典国产精品| 久久人人爽av亚洲精品天堂| 亚洲三级黄色毛片| 国产成人91sexporn| 国产一级毛片在线| 亚洲精品日韩在线中文字幕| 亚洲av免费高清在线观看| 国产白丝娇喘喷水9色精品| av一本久久久久| 中文字幕久久专区| 99久久中文字幕三级久久日本| 男人狂女人下面高潮的视频| 亚洲精品自拍成人| 99久久精品国产国产毛片| 女人久久www免费人成看片| 亚洲精品国产色婷婷电影| 亚洲三级黄色毛片| 亚洲av男天堂| 亚洲人成网站在线播| 免费人妻精品一区二区三区视频| 国产淫片久久久久久久久| 99视频精品全部免费 在线| 日日爽夜夜爽网站| 亚洲无线观看免费| √禁漫天堂资源中文www| 国产色爽女视频免费观看| 视频区图区小说| 建设人人有责人人尽责人人享有的| 国产精品一区二区性色av| 午夜福利影视在线免费观看| 国产精品女同一区二区软件| 三级国产精品片| 欧美精品人与动牲交sv欧美| 欧美日韩国产mv在线观看视频| 欧美激情极品国产一区二区三区 | 亚洲欧美日韩东京热| 亚洲欧洲日产国产| 免费黄频网站在线观看国产| 搡女人真爽免费视频火全软件| 亚洲成人一二三区av| 午夜日本视频在线| 丝袜脚勾引网站| 日本黄大片高清| 一区二区三区乱码不卡18| 性高湖久久久久久久久免费观看| 久久人人爽av亚洲精品天堂| 国产成人午夜福利电影在线观看| 大片电影免费在线观看免费| 中文在线观看免费www的网站| 搡女人真爽免费视频火全软件| 久久久久久久国产电影| 欧美日韩视频精品一区| 国产乱来视频区| 街头女战士在线观看网站| 少妇 在线观看| 男的添女的下面高潮视频| 国产黄色视频一区二区在线观看| 成人毛片60女人毛片免费| av女优亚洲男人天堂| 人妻夜夜爽99麻豆av| 国产精品女同一区二区软件| 午夜免费男女啪啪视频观看| 日本猛色少妇xxxxx猛交久久| 99热这里只有是精品50| 99久久中文字幕三级久久日本| 精品亚洲成a人片在线观看| 日本vs欧美在线观看视频 | 免费观看a级毛片全部| 一级黄片播放器| 男女边摸边吃奶| 最近的中文字幕免费完整| 22中文网久久字幕| 国产黄色视频一区二区在线观看| 我的老师免费观看完整版| 久久久久网色| 国产成人精品无人区| 欧美日韩av久久| 日本-黄色视频高清免费观看| 午夜福利影视在线免费观看| 全区人妻精品视频| 欧美日韩国产mv在线观看视频| 欧美日韩综合久久久久久| 99久久综合免费| 三级经典国产精品| 91精品国产国语对白视频| 看十八女毛片水多多多| 男男h啪啪无遮挡| 伦理电影大哥的女人| 日日摸夜夜添夜夜爱| 人人妻人人澡人人爽人人夜夜| 好男人视频免费观看在线| 亚洲成色77777| 99久久综合免费| 99精国产麻豆久久婷婷| 亚洲欧美清纯卡通| 欧美xxⅹ黑人| 久久久久精品性色| 亚洲av中文av极速乱| 建设人人有责人人尽责人人享有的| 日本与韩国留学比较| 在线免费观看不下载黄p国产| 亚洲精品视频女| 少妇 在线观看| 欧美亚洲 丝袜 人妻 在线| 麻豆乱淫一区二区| 免费人妻精品一区二区三区视频| 麻豆乱淫一区二区| 黄片无遮挡物在线观看| 婷婷色av中文字幕| 五月伊人婷婷丁香| 国产爽快片一区二区三区| 一区二区av电影网| 丁香六月天网| 精品午夜福利在线看| 99久久精品一区二区三区| 丰满乱子伦码专区| 尾随美女入室| 欧美日韩视频高清一区二区三区二| 十八禁网站网址无遮挡 | 十八禁网站网址无遮挡 | 精品酒店卫生间| 热re99久久国产66热| 国产熟女欧美一区二区| 久久久久网色| 日韩av不卡免费在线播放| 免费在线观看成人毛片| 欧美精品一区二区免费开放| 国产精品久久久久久久久免| 又大又黄又爽视频免费| 又爽又黄a免费视频| 18禁裸乳无遮挡动漫免费视频| 男女无遮挡免费网站观看| 国产亚洲最大av| 高清视频免费观看一区二区| av国产精品久久久久影院| 性高湖久久久久久久久免费观看| 大片免费播放器 马上看| 国产精品久久久久久久电影| 亚洲国产精品成人久久小说| 国产精品福利在线免费观看| 中文字幕制服av| 18禁裸乳无遮挡动漫免费视频| 成年女人在线观看亚洲视频| 亚洲精品国产色婷婷电影| 日日摸夜夜添夜夜添av毛片| 亚洲,一卡二卡三卡| 免费av不卡在线播放| 日韩免费高清中文字幕av| av网站免费在线观看视频| 免费看av在线观看网站| 久久久久精品久久久久真实原创| 超碰97精品在线观看| 最黄视频免费看| 中文字幕人妻熟人妻熟丝袜美| 大片电影免费在线观看免费| 国产精品不卡视频一区二区| 天美传媒精品一区二区| 黄色日韩在线| 成年人午夜在线观看视频| 婷婷色麻豆天堂久久| 如何舔出高潮| 日韩大片免费观看网站| 人妻 亚洲 视频| 校园人妻丝袜中文字幕| 久久国产精品大桥未久av | 欧美日韩一区二区视频在线观看视频在线| 国产淫片久久久久久久久| 一区二区三区精品91| 欧美成人午夜免费资源| 中文精品一卡2卡3卡4更新|