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

    同軸槍脈沖放電等離子體輸運過程中密度變化的實驗研究

    2017-08-01 17:16:04楊亮張俊龍閆慧杰滑躍任春生
    物理學(xué)報 2017年5期
    關(guān)鍵詞:氬氣同軸電離

    楊亮 張俊龍 閆慧杰 滑躍 任春生

    (大連理工大學(xué)物理與光電工程學(xué)院,三束材料改性教育部重點實驗室,大連 116024)

    同軸槍脈沖放電等離子體輸運過程中密度變化的實驗研究

    楊亮 張俊龍 閆慧杰 滑躍 任春生?

    (大連理工大學(xué)物理與光電工程學(xué)院,三束材料改性教育部重點實驗室,大連 116024)

    (2016年8月27日收到;2016年12月7日收到修改稿)

    同軸槍放電等離子體具有密度高、輸運速度快等特點,在核物理、航天工程等領(lǐng)域具有廣闊的應(yīng)用前景,已成為國際前沿研究熱點.同軸槍中的等離子體密度是反映其應(yīng)用特性的重要參數(shù)之一,因此等離子體密度在輸運過程中的變化對理論研究和實際應(yīng)用都具有重要意義.利用發(fā)射光譜法測量了Hβ譜線的Stark展寬,從而計算出同軸槍放電等離子體密度在輸運過程中的變化.結(jié)果顯示,當(dāng)電源注入能量為1.08 kJ、同軸槍內(nèi)空氣氣壓為4.0 Pa時,等離子體密度在輸運過程中不斷增加;相同能量注入條件下,當(dāng)同軸槍內(nèi)空氣氣壓增加至10 Pa時,等離子體密度在輸運過程中出現(xiàn)了先增加后減小的趨勢;當(dāng)電源的注入能量達(dá)到7.68 kJ時,等離子體密度在10 Pa氣壓條件下輸運時也出現(xiàn)了一直增加的現(xiàn)象.此外,當(dāng)同軸槍內(nèi)的工作氣體變?yōu)闅鍤鈺r,在注入能量為1.08 kJ、槍內(nèi)氣壓4.0 Pa條件下,等離子體密度在輸運過程中一直減小.

    同軸槍,等離子體密度,Stark展寬

    1 引 言

    同軸槍通常是由兩個同軸電極組成(一個實心圓柱陽極和一個空心圓柱陰極),它通過兩電極間瞬時的強(qiáng)放電產(chǎn)生一個磁化的高密度等離子體團(tuán),該等離子體團(tuán)在電流流過同軸槍時產(chǎn)生的洛倫茲力作用下被加速,很容易達(dá)到幾千米每秒甚至幾十千米每秒的速度噴射出去[1?3].由于其高密度、高輸運速度等特點[1,4?10],已經(jīng)被廣泛應(yīng)用于聚變裝置[2,11,12],等離子體空間推進(jìn)[13?16]、磁約束核聚變中的再燃[17,18]、塵埃粒子加速[1,4,19?21]以及空間微小碎片超高速撞擊的實驗?zāi)M研究[3,22]中.

    在同軸槍放電等離子體的研究中,等離子體密度與輸運速度是非常重要的研究內(nèi)容.在我們前期開展的研究工作中,對同軸槍放電等離子體的軸向輸運速度及其影響因素進(jìn)行了研究[23],而放電等離子體密度及其隨輸運過程的變化是另一個重要參數(shù),不僅決定著放電等離子體應(yīng)用的物理、化學(xué)過程,而且在很大程度上決定著其應(yīng)用特性.所以,研究同軸槍放電中的等離子體密度及變化具有非常重要的意義.

    目前等離子體密度診斷所采用的方法主要有微波法、Thomson散射法、Langmuir探針法以及光譜法[1,4,24?31].微波法和Thomson散射法[25,26]的操作較為復(fù)雜,實際中通常采用Langmuir探針法和光譜法.Langmuir探針法需要測量等離子體的伏安特性[27,28],但由于探針直接與等離子體相接觸,很容易對等離子體造成干擾,而且探針暴露于強(qiáng)電離環(huán)境下也很容易損壞(尤其對于高密度、高能量密度等離子體環(huán)境).相比于探針法,發(fā)射光譜技術(shù)診斷等離子體密度具有操作簡單、選擇性好、靈敏度和準(zhǔn)確度高以及對等離子體本身無干擾等特點,此外,該技術(shù)不僅能實現(xiàn)原位診斷,而且能夠得到時間及空間分辨信息,所以被廣泛應(yīng)用于等離子體的診斷中.Ticos等[1,4]利用Dα的Stark展寬效應(yīng)測量了同軸槍放電等離子體塵埃粒子加速中的等離子體密度,Witherspoon等[31]測量了裝置放電時的Hβ譜線寬度進(jìn)而得到了等離子體密度.然而,在以往諸多的實驗研究中,對等離子體密度的測量往往集中于某一特定位置,而對同軸槍放電等離子體密度在輸運過程中的變化特性的研究相對較少.事實上,由于雪犁模型(等離子體在輸運過程中,不僅會繼續(xù)電離與其接觸的中性粒子,而且會推著電離后的粒子一起加速)[32?34]的提出,等離子體的密度在輸運過程中是變化的,而且該變化能夠決定其裝置的應(yīng)用特性,所以研究同軸槍放電等離子體密度在輸運過程中的變化是十分必要的.基于此,本文利用Hβ譜線Stark展寬效應(yīng)研究了等離子體密度及其在輸運過程中的變化特性,對不同放電條件(不同的放電氣體和壓強(qiáng))下等離子體密度及其在輸運過程中的變化特性也做了比較詳細(xì)的對比研究.

    2 實驗裝置

    實驗研究裝置如圖1所示,主要由真空腔、供電系統(tǒng)、同軸放電槍、輸運通道以及觀察窗口組成.真空腔(圖1左側(cè)球形裝置)為外徑55 cm的空心圓柱體,通過一個組合真空機(jī)組得到本底真空水平為10?3Pa.供電系統(tǒng)通過一個幅值可調(diào)電源給240μF的電容器充電,再通過點火開關(guān)將幾千伏(最大值10 kV)的電壓供給同軸放電槍,可以得到瞬時幾千甚至幾萬安培峰值的放電電流.同軸放電槍由兩個電極組成,一個是同軸槍陽極(直徑2.0 cm,長24 cm的實心銅棒),另一個是同軸槍陰極(內(nèi)徑8 cm,厚度0.3 cm,長24 cm的不銹鋼空心圓柱),陽極與陰極之間用一厚3.0 cm、長4.0 cm的尼龍介質(zhì)隔離,所以放電時等離子體的實際加速距離為20 cm.將同軸放電槍置入輸運通道(內(nèi)徑14.6 cm,長85 cm的不銹鋼空心圓柱)內(nèi),中間用一層3.0 cm厚的聚四氟乙烯材料隔離,該輸運通道與真空腔連接在一起.同時,在輸運通道處設(shè)置一長24 cm、寬10 cm的石英玻璃觀察窗,同軸放電槍輸出端面與觀察窗口的距離約為6.0 cm.觀察窗口設(shè)有A,B,C三個光譜采集點,用于測量等離子體在輸運過程中的密度變化.同時,在該三點處采用快速響應(yīng)高量子效率的光電倍增管(P30A-05)測量放電等離子體發(fā)光信號,以光信號的時間差以及兩個光電倍增管的距離來計算放電等離子體的軸向輸運速度.三個采集點與同軸槍輸出端面的距離分別為8,17,26 cm.采用Acton SpectraPro 2500i型光柵光譜儀測量Hβ譜線并計算其展寬,所選光柵的密度為1200 line/mm.由于同軸槍脈沖放電時間短,產(chǎn)生的等離子體輸運速度很快,所以測量時將光譜儀的積分時間設(shè)置為3.0 s,以保證可以采集到等離子體的發(fā)光光譜.

    圖1 (網(wǎng)刊彩色)實驗裝置示意圖Fig.1.(color online)Experimental setup.

    電信號測量回路及原理如圖2所示,電壓與電流信號分別用Tektronix高壓探頭與Pearson電流探頭測得并輸入高存儲容量數(shù)字示波器(Tektronix DPO 4104)進(jìn)行記錄與分析.高壓探頭置于同軸槍高壓陽極端,Pearson電流探頭置于接地陰極端.實驗過程中,首先需要給電容器兩端充電,充電后的電容器儲存大量的電能,而后通過高壓開關(guān)閉合瞬時供給同軸放電槍使槍內(nèi)氣體擊穿,產(chǎn)生放電等離子體.由于流過同軸槍陽極的軸向電流感生出很強(qiáng)的方位角磁場B,該磁場與放電時的徑向電流J相互作用,形成軸向的洛倫茲力F(F=J×B),該洛倫茲力能使等離子體沿著軸向運動,從而形成沿著軸向輸運的等離子體.

    圖2 (網(wǎng)刊彩色)實驗測量原理圖Fig.2.(color online)Schematic of experimental measurement principle.

    3 結(jié)果與討論

    3.1 放電特性

    選用電源注入能量為1.08 kJ(電容器兩端電壓達(dá)到3.0 kV),空氣氣壓4.0 Pa作為放電條件,可得到圖3所示電壓、電流及光電流波形圖.從圖中可以看出,當(dāng)電源注入能量為1.08 kJ時,放電電流達(dá)到8.5 kA,電流波形很明顯地呈現(xiàn)出多次放電特性,多次放電特性能夠同時從等離子體光電流信號中得到證實.對比電流波形與光電流波形可以看出,電流波形的每一個波峰與波谷分別對應(yīng)一個光信號的波峰,且光信號在經(jīng)過A,B,C三點時依次出現(xiàn)一個峰值,這就證明確實有多次放電發(fā)生且等離子體在運動.多次放電的原因是電源電路阻抗不匹配,即電源電路中不同電容、電感、電阻回路產(chǎn)生的RLC振蕩相互影響,最后在電極間形成多次擊穿而引起多次放電[21,23,24].

    觀察光電流波形,發(fā)現(xiàn)其信號具有一定的寬度,這是因為同軸槍內(nèi)噴射出的等離子體有一定的體積(通常稱等離子體團(tuán)[6,7,12])且在輸運時會在一定的時間段內(nèi)一直出現(xiàn)在光信號接收口的前端.此外,對應(yīng)一次放電等離子體的光信號均出現(xiàn)兩個波峰,通常認(rèn)為Blow-by不穩(wěn)定性是造成該現(xiàn)象的原因[31,35].Blow-by不穩(wěn)定性來源于徑向磁壓分布的不均勻性,如果放電腔室內(nèi)的放電電流層最初為徑向均勻的密度和厚度分布,它沿著內(nèi)電極的加速度會高于沿外電極的加速度(靠近內(nèi)電極的磁場較強(qiáng),受到的洛倫茲力較大,而靠近外電極的磁場相對較弱,受到洛倫茲力較小).該加速度的差異會使載流等離子體產(chǎn)生向后的傾斜,傾斜的電流層會產(chǎn)生一個大的軸向電流分量,進(jìn)而形成一個徑向的洛倫茲力.該洛倫茲力使越來越多的等離子體向低速區(qū)轉(zhuǎn)移,便會在靠近內(nèi)電極處出現(xiàn)一塊質(zhì)量較小但速度很快的等離子體團(tuán),而靠近外電極處出現(xiàn)一塊質(zhì)量較大但速度較慢的等離子體團(tuán),正如圖3中A,B處測量的光電流信號中首先出現(xiàn)一個寬度較小的波峰,而后出現(xiàn)一個寬度較大的波峰.

    從A,B,C處測量的光電流信號還可以看出光信號的寬度在輸運過程中逐漸變窄,這說明等離子體團(tuán)在輸運過程中逐漸壓縮.另外,三個光電流的時間間隔分別為9.3μs和11.7μs,A與B,B與C之間的距離均為9 cm,由此可以計算出等離子體團(tuán)在A與B間的輸運速度約為9.7 km/s,而在B與C間的輸運速度約為7.7 km/s.可以看出等離子體團(tuán)在輸運過程中是減速的.導(dǎo)致等離子體團(tuán)速度減慢的原因主要是雪犁電離模型[32?34].從雪犁模型的基本假設(shè)可知,初始時刻等離子體受洛倫茲力作用開始加速運動,但在輸運過程中,等離子體中的高能電子會電離部分中性氣體分子并拖拽被電離的中性氣體分子一起運動,致使等離子體的質(zhì)量在輸運過程中不斷增大.根據(jù)動量守恒定律,被拖拽電離粒子運動的動量與等離子體減少的動量是相等的,所以等離子體軸向速度在輸運過程中逐漸減慢.事實上,同軸槍放電等離子體均會經(jīng)過一個先加速后減速的過程[36?38],因為在等離子體最初的加速過程中,等離子體的質(zhì)量m和速度v都是變化量,與洛倫茲力F滿足沖量定理,即(t為時間),在電流上升階段,洛倫茲力單調(diào)增加,等離子體的質(zhì)量和速度都隨之增大;在電流達(dá)到峰值并開始下降時,洛倫茲力也按照類似的趨勢達(dá)到極大并開始衰減,而等離子體質(zhì)量又是單調(diào)增大的,因此洛倫茲力F不足以維持質(zhì)量遞增的等離子體高速運動只能通過來反向補(bǔ)償,因此等離子體的速度v開始減小.一般該過程發(fā)生在同軸槍內(nèi)部,而在同軸槍外部,由于洛倫茲力消失,輸運過程中的等離子體速度應(yīng)該是逐漸減小的.

    圖3 (網(wǎng)刊彩色)同軸槍在注入能量為1.08 kJ,空氣氣壓為4.0 Pa放電條件下的電壓、電流及光電流波形圖Fig.3.(color online)Voltage,current and photomultiplier acquired signal waveform for coaxial gun discharge in air with applied power energy of 1.08 kJ and pressure of 4.0 Pa.

    圖4(a)所示為電源注入能量為1.08 kJ,空氣氣壓為10 Pa條件下的電壓、電流波形圖.從圖中可以看出,與相同注入能量下4.0 Pa時的空氣放電相比,空氣氣壓為10 Pa時產(chǎn)生的放電電流更大,其峰值可達(dá)10 kA以上.圖4(b)為電源注入能量為1.08 kJ、氬氣氣壓為4.0 Pa時的放電電壓、電流波形圖,與空氣放電相比,氬氣放電電流變大,放電擊穿電壓變低.圖4中的放電特性為下一步解釋等離子體團(tuán)在不同工作氣體條件下等離子體密度及其隨輸運距離的變化特性提供了依據(jù).

    圖4 (網(wǎng)刊彩色)同軸槍在(a)空氣,注入能量為1.08 kJ,氣壓為10 Pa和(b)氬氣,注入能量為1.08 kJ,氣壓為4.0 Pa放電條件下的電壓、電流波形圖Fig.4.(color online)Voltage and current waveforms for coaxial gun discharge in(a)air with applied power energy of 1.08 kJ and pressure of 10 Pa and(b)argon with applied power energy of 1.08 kJ and pressure of 4.0 Pa.

    3.2 光譜學(xué)診斷

    譜線的展寬一般有自然展寬(?λnat)、共振展寬(?λrev)、van der Waals展寬(?λvan)、Stark展寬(?λStark)、多普勒展寬(?λDoppler)和儀器展寬(?λInstrument)[39,40].前四種展寬為洛倫茲線型,多普勒展寬為高斯線型,而儀器展寬的線型取決于儀器本身.實驗使用的光譜儀通過測量得出的是高斯線型,所以高斯線型的總展寬(?λG)由儀器展寬和多普勒展寬組成[29,39,40],由此,可以得出,

    式中的儀器展寬可由He-Ne激光器中心波長為632.80 nm的譜線測出,實驗中測量的儀器展寬為0.0989 nm.多普勒展寬計算公式為

    式中λ0為對應(yīng)的吸收峰值,Tg為氣體溫度(單位:K),A為氣體的相對分子質(zhì)量[29,30,39].

    由于儀器展寬屬于高斯展寬(?λG),而Stark展寬屬于洛倫茲展寬(?λL),在實驗中所采集到的Hβ譜線,既有均勻展寬(洛倫茲展寬),也有非均勻展寬(高斯展寬),所以實驗所得到譜線的線型為佛克脫線型.在由均勻展寬和非均勻展寬共同導(dǎo)致的佛克脫線型中,洛倫茲線型展寬、高斯線型展寬和總的佛克脫線型展寬之間的關(guān)系為[40]

    式中的?λv為佛克脫線型展寬寬度.在該實驗條件下(通常是幾帕水平),自然展寬、共振展寬和van der Waals展寬可以忽略不計,所以洛倫茲線型展寬主要為Stark展寬,可將洛倫茲展寬近似當(dāng)作Stark展寬,具體值為

    而等離子體密度ne(單位:cm?3)與Stark展寬寬度的關(guān)系為[30,39]

    (4)式中,?λv可以通過擬合實驗所測Hβ譜線得到.而對于高斯展寬,可以通過(1)式得出.但是,在一些文獻(xiàn)中[1,4],即使同軸槍放電電流達(dá)200 kA,離子溫度為2.8 eV時(此時的多普勒展寬約為0.06 nm),在計算等離子體密度時多普勒展寬也未被計算在內(nèi).本實驗中,同軸槍放電電流在10 kA左右,產(chǎn)生的離子溫度小于2.8 eV,而氣體溫度低于離子溫度,所以,本實驗的多普勒展寬應(yīng)小于0.06 nm,該展寬的寬度遠(yuǎn)小于本實驗中的最小Stark展寬0.2186 nm,通過(1)式的計算關(guān)系可以得出,該展寬的寬度對高斯展寬寬度的影響很小.

    因此,多普勒展寬在本實驗中亦可以忽略,即高斯展寬就是儀器展寬.

    該實驗在觀察窗口B處測量的光譜如圖5所示,其中可以明顯看到486.13 nm處的Hβ譜線,該展寬為佛克脫線型展寬.圖6是將圖5中的Hβ譜線展開,同時對其進(jìn)行線型擬合,得到了展寬為0.2949 nm的佛克脫線型展寬寬度,高斯展寬為0.0989 nm,這時由(4)式可得Stark展寬為0.2617 nm,進(jìn)而計算出B處的等離子體密度為1.5×1015cm?3.

    圖5 同軸槍在注入能量為1.08 kJ、空氣氣壓為4.0 Pa放電條件下觀察窗口B處的Hβ光譜譜線Fig.5.Spectrum of Hβ-region measured at the position B for coaxial gun discharge in air with applied power energy of 1.08 kJ and pressure of 4.0 Pa.

    圖6 同軸槍在注入能量為1.08 kJ、空氣氣壓為4.0 Pa放電條件下觀察窗口B處Hβ譜線的擬合Fig.6.Fitting of Hβ-line measured at the position B for coaxial gun discharge in air with applied power energy of 1.08 kJ and pressure of 4.0 Pa.

    3.3 注入能量為1.08 kJ、氣壓為4.0 Pa空氣條件下等離子體團(tuán)在輸運過程中的密度變化

    圖7所示為等離子體團(tuán)在注入能量為1.08 kJ、氣壓為4.0 Pa空氣條件下輸運經(jīng)過A,B,C三點時測量的Hβ譜線展寬.對該展寬進(jìn)行歸一化處理,歸一化后的佛克脫線型展寬寬度分別為0.2567,0.2949,0.3132 nm,經(jīng)過(4)式得到Stark展寬寬度分別為0.2186,0.2617,0.2820 nm,從而計算出A,B,C三點處的等離子體團(tuán)密度分別為1.14×1015,1.5×1015,1.68×1015cm?3.由此可以看出,等離子體團(tuán)的密度在輸運過程中會逐漸增大.由雪犁模型可知[32?34],等離子體團(tuán)在輸運過程中會將一部分能量轉(zhuǎn)移到與其接觸的中性粒子薄層中,從而將中性粒子電離,并推著被電離的粒子運動.此外,由于等離子體團(tuán)的速度非常高,且剛被電離的帶電粒子在軸向方向并沒有定向速度,后面的高速等離子體團(tuán)將其加速的過程中也會存在一種堆積效應(yīng),使帶電粒子聚集,產(chǎn)生碰撞等離子體激波[38],該激波沿軸向傳播使等離子體團(tuán)逐漸被壓縮,等離子體團(tuán)密度在輸運過程中升高,這是一個既電離又壓縮的過程.但是,等離子體團(tuán)在電離前面的中性粒子時會不斷地消耗其能量,且等離子體團(tuán)的能量在輸運過程中也會因自身的淬滅而不斷降低,即等離子體團(tuán)在輸運過程中的能量會逐漸減小,其電離能力也會逐漸減弱.這也就解釋了為什么上述實驗中等離子體團(tuán)在B—C段增加的密度小于其在A—B段增加的密度,即密度增長變緩慢的原因.

    事實上,圖3中的光電流波形能很好地驗證等離子體團(tuán)在輸運過程中存在的壓縮效應(yīng),并且能夠解釋等離子體團(tuán)在輸運過程中減速的原因.圖中的光電流波形是同一等離子體團(tuán)經(jīng)過A,B,C三點時光電倍增管測量得到的光電流波形.通過比較這三處的光電流波形可看出,光電流的波峰在輸運過程中逐漸變窄變高.光電流的波峰變窄說明等離子體團(tuán)的體積縮小(等離子體團(tuán)被壓縮),強(qiáng)度增強(qiáng)說明等離子體團(tuán)的聚集程度升高,這應(yīng)該是一個明顯的聚集過程,且該過程可能也是等離子體碰撞激波形成的過程.此外,根據(jù)雪犁模型描述,等離子體團(tuán)在輸運過程中的質(zhì)量會增加,而圖3中的光電流可以明顯看出等離子體團(tuán)在輸運過程中逐漸被壓縮,進(jìn)而使等離子體團(tuán)的體積逐漸減小,因此,等離子體團(tuán)密度增加.同時,由于等離子體團(tuán)質(zhì)量增大且被壓縮的緣故,可能會形成等離子體碰撞激波,從而使等離子體團(tuán)在輸運過程中的速度逐漸減小.

    圖7 同軸槍在注入能量為1.08 kJ、空氣氣壓為4.0 Pa放電條件下觀察窗口A,B,C處測得的Hβ譜線展寬Fig.7.Broadenings of Hβ-lines measured at the positions A,B,C for coaxial gun discharge in air with applied power energy of 1.08 kJ and pressure of 4.0 Pa.

    3.4 注入能量為1.08 kJ、氣壓為10 Pa空氣條件下等離子體團(tuán)在輸運過程中的密度變化

    圖8所示為同軸槍內(nèi)空氣氣壓為10 Pa、注入能量為1.08 kJ時,在A,B,C三點處測量得到的佛克脫線型展寬(歸一化后),該展寬的寬度分別為0.3299,0.4001,0.3476 nm,Stark展寬寬度分別為0.3002,0.3757,0.3195 nm,計算出等離子體密度分別為1.84×1015,2.57×1015,2.02×1015cm?3.可以看出,該條件下放電產(chǎn)生的等離子體團(tuán)密度大于氣壓為4.0 Pa時的等離子體團(tuán)密度,并且該等離子體團(tuán)的密度在輸運過程中并不是一直增加的,而是先增加后減小的.事實上,當(dāng)氣壓升高后,同軸槍內(nèi)的粒子數(shù)密度增加,則相同電壓條件下放電會更強(qiáng)(圖4(a)),電離的粒子數(shù)會增加,等離子體團(tuán)的密度高于4.0 Pa時的等離子體團(tuán)密度.當(dāng)?shù)入x子體團(tuán)在洛倫茲力作用下向前運動時,等離子體團(tuán)就會一邊使氣體繼續(xù)電離,使等離子體團(tuán)密度提高,一邊不斷地消耗其能量,使其電離能力減弱.這意味著在一定條件下,當(dāng)?shù)入x子體團(tuán)的能量消耗到一定程度,即電離率低于淬滅率時,等離子體團(tuán)的密度便開始下降.通過該分析,不難解釋等離子體團(tuán)在10 Pa條件下輸運過程中出現(xiàn)的現(xiàn)象.首先,在等離子體團(tuán)從A運動到B的過程中,等離子體團(tuán)密度增加,說明在該過程中等離子體團(tuán)的能量足夠高,可以繼續(xù)電離中性粒子.但當(dāng)?shù)入x子體團(tuán)從B運動到C時,等離子體團(tuán)密度減小,表明等離子體團(tuán)的電離率已經(jīng)小于其淬滅率.這意味著氣壓的增大會使等離子體團(tuán)在前進(jìn)的過程中能量消耗得更快.因為背景氣體粒子數(shù)密度的提高,使等離子體團(tuán)在經(jīng)過相同距離時需要電離更多的中性粒子,即等離子體團(tuán)的能量消耗更多,盡管在4.0 Pa時等離子體團(tuán)在輸運過程中密度在不斷增加,但是等離子體團(tuán)的密度不可能一直增加下去,只是相對氣壓為10 Pa時等離子體團(tuán)密度增加所持續(xù)的距離更長,而且恰好超過了觀察窗口的距離.

    圖8 同軸槍在注入能量為1.08 kJ,空氣氣壓為10 Pa放電條件下觀察窗口A,B,C處測得的Hβ譜線展寬Fig.8.Broadenings of Hβ-lines measured at the positions A,B,C for coaxial gun discharge in air with applied power energy of 1.08 kJ and pressure of 10 Pa.

    基于上述分析,可將等離子體團(tuán)的輸運過程分為兩個階段.在第一階段,等離子體團(tuán)的能量高,在輸運過程中會不斷電離中性粒子,使等離子體團(tuán)的密度逐漸升高,但與此同時,等離子體團(tuán)的能量在逐漸下降,即電離能力逐漸減弱.在第二階段,等離子體團(tuán)的能量已經(jīng)降低到其電離率小于淬滅率的水平,等離子體團(tuán)的密度會逐漸降低.第一階段與第二階段的臨界點是等離子體團(tuán)密度最大值出現(xiàn)的時刻.

    3.5 注入能量為7.68 kJ、氣壓為10 Pa空氣條件下等離子體團(tuán)在輸運過程中的密度變化

    為了進(jìn)一步分析上述同軸槍脈沖放電等離子體團(tuán)在輸運過程中的密度變化,實驗中將電源的注入能量增大到7.68 kJ(電容器兩端電壓達(dá)到8.0 kV),測量了空氣氣壓10 Pa條件下,等離子體團(tuán)輸運經(jīng)過A,B,C三點時的密度變化趨勢,如圖9所示.在A,B,C三點處測量得到的佛克脫線型展寬(歸一化后)寬度分別為0.49406,0.60367,0.66076 nm,并通過(4)式分別得到了0.4743,0.5875,0.646 nm的Stark展寬,從而計算出等離子體密度分別為3.65×1015,5.03×1015,5.81×1015cm?3.從該密度的變化趨勢可以看出,當(dāng)電源的注入能量增長至7.68 kJ時,等離子體團(tuán)的密度大于在1.08 kJ時的密度,且在輸運過A,B,C三點時又一次呈現(xiàn)出一直增加的趨勢.這是因為當(dāng)電源注入能量升高后,同軸槍內(nèi)的放電會更強(qiáng),使電離的粒子數(shù)增加,同時使等離子體團(tuán)具有更高的能量.該等離子體團(tuán)在輸運過程中也會發(fā)生雪犁電離,使等離子體團(tuán)的密度升高,能量不斷被損耗.但是,由于該等離子體團(tuán)初始能量高于注入能量為1.08 kJ時產(chǎn)生的等離子體團(tuán)能量,所以在輸運相同的距離時,即經(jīng)過A—B段后,仍具有很高的能量(很強(qiáng)的電離能力),導(dǎo)致該等離子體團(tuán)在經(jīng)過B—C段時能夠繼續(xù)電離中性粒子使等離子體團(tuán)的密度繼續(xù)增加.不過,等離子體團(tuán)在B—C段增加的密度小于其在A—B段增加的密度,即出現(xiàn)了密度增長變緩慢的趨勢,這與注入能量1.08 kJ、氣壓4.0 Pa時等離子體團(tuán)的密度變化趨勢是相同的,也說明等離子體團(tuán)的密度不能一直增加下去.從該結(jié)果可以看出,當(dāng)?shù)入x子體團(tuán)的初始能量足夠高時,等離子體團(tuán)的密度不僅會增加,電離能力也會增強(qiáng),并且該電離能力會維持更長的距離,使等離子體團(tuán)的密度在輸運經(jīng)過更長的距離時一直增加,也就是說,等離子體團(tuán)的密度及持續(xù)電離的能力都取決于等離子體團(tuán)最初所攜帶的能量.但是無論等離子體團(tuán)最初所攜帶的能量有多大,只要在輸運足夠遠(yuǎn)的距離后,便會達(dá)到電離率小于淬滅率的水平而密度開始減小.

    圖9 同軸槍在注入能量為7.68 kJ、空氣氣壓為10 Pa放電條件下觀察窗口A,B,C處測得的Hβ譜線展寬Fig.9.Broadenings of Hβ-lines measured at the positions A,B,C for coaxial gun discharge in air with applied power energy of 7.68 kJ and pressure of 10 Pa.

    3.6 注入能量為1.08 kJ、氣壓為4.0 Pa氬氣條件下等離子體團(tuán)在輸運過程中的密度變化

    當(dāng)同軸槍在氬氣氣壓4.0 Pa、注入能量1.08 kJ條件下放電時,等離子體團(tuán)密度在輸運過程中的變化趨勢如圖10所示.從圖中所示佛克脫線型展寬可以看出,在相同氣壓4.0 Pa條件下,等離子體團(tuán)在輸運過程中的密度變化趨勢與空氣中明顯不同,即在空氣中密度逐漸增加,而在氬氣中逐漸減小.這意味著氬氣放電產(chǎn)生的等離子體在距離同軸槍很近的距離就發(fā)生了聚集現(xiàn)象.引起該現(xiàn)象的原因可能是氣體電離能的差異以及不同種類氣體放電中的反應(yīng).首先,從圖3和圖4(b)中的放電情況可知,相同條件下氬氣的放電強(qiáng)于空氣,這能很好地解釋為什么距離同軸槍最近的A點所測得的氬氣等離子體團(tuán)密度大于相同氣壓條件4.0 Pa時空氣中的等離子體團(tuán)密度(氬氣等離子體團(tuán)密度為1.75×1015cm?3,空氣等離子體團(tuán)密度為1.14×1015cm?3).此外,氬原子的電離能為15.7 eV,而作為空氣中主要成分的氮分子和氧分子的電離能分別為15.58 eV和12.22 eV[41],說明氬氣的電離能高于空氣的電離能.對于一定的等離子體電子能量分布函數(shù)而言,兩者(氬氣、空氣)電離能的差異可能會使等離子體團(tuán)不能繼續(xù)電離氬氣,但卻能使空氣繼續(xù)電離,從而導(dǎo)致等離子體團(tuán)密度在氬氣輸運中逐漸降低,而在空氣輸運中逐漸增大.也可能是由于氬氣電離能高,等離子體團(tuán)在輸運過程中更容易達(dá)到電離率小于淬滅率的水平,進(jìn)而使等離子體團(tuán)密度在輸運較短距離后便開始降低.反之,由于空氣的電離能低,輸運過程中消耗的能量較氬氣中低,所以需要更長的距離才能達(dá)到電離率小于淬滅率的水平,因而等離子體團(tuán)密度在更長的輸運距離中一直增加,本實驗中該距離超過了觀察窗口的距離.引起空氣等離子體團(tuán)密度在窗口觀察距離中一直增加的另一個原因也可能是解離反應(yīng)[42?44],e+N2→2N+e,e+O2→2O+e.

    圖10 同軸槍在注入能量為1.08 kJ、氬氣氣壓為4.0 Pa放電條件下觀察窗口A,B,C處測得的Hβ譜線展寬Fig.10.Broadenings of Hβ-lines measured at the positions A,B,C for coaxial gun discharge in argon with applied power energy of 1.08 kJ and pressure of 4.0 Pa.

    該反應(yīng)能夠使粒子數(shù)增加,且解離后氮、氧原子的電離能(分別為14.5和13.6 eV)[41]低于氬原子的電離能.因此,相比較氬氣等離子體團(tuán)輸運過程而言,空氣等離子體團(tuán)不僅仍具有較高的電離能力,而且因為電離了更多數(shù)量的原子而使密度不斷增加.

    4 結(jié) 論

    本文主要應(yīng)用Hβ譜線的Stark展寬法測量了同軸槍脈沖放電等離子體團(tuán)在輸運過程中的密度變化,并研究了不同參數(shù)下該密度的變化趨勢.結(jié)果發(fā)現(xiàn),等離子體團(tuán)在空氣中輸運時密度會隨著氣壓條件的不同而有差異,在低氣壓時密度一直增加,在高氣壓時密度先增加后減小.同時發(fā)現(xiàn),在高氣壓下如果提高電源注入能量,等離子體團(tuán)密度在輸運過程中也會一直增加.此外,當(dāng)氣壓與放電條件相同時,等離子體團(tuán)密度在氬氣中一直減小.這些現(xiàn)象出現(xiàn)的原因可歸結(jié)為高速運動的等離子體團(tuán)既能繼續(xù)電離前方的中性粒子,又要消耗其自身能量的特性.不同的氣體密度、種類以及等離子體團(tuán)最初攜帶的能量都會導(dǎo)致輸運過程中能量損耗程度的不同,從而使不同參數(shù)下等離子體團(tuán)密度最大值出現(xiàn)的位置有了很大的差異.該實驗結(jié)果中聚集現(xiàn)象的出現(xiàn)表明了等離子體團(tuán)在輸運過程中的物理特性,這些研究加深了我們對同軸槍放電等離子體輸運過程特性的認(rèn)識與了解.

    [1]Ticos C M,Wang Z H,Wurden G A,Kline J L,Montgomery D S 2008Phys.Plasmas15 103701

    [2]Turchi P J,Roderick N F,Degnan J H,Frese M H,Amdahl D J 2008IEEE Tran.Plasma Sci.36 92

    [3]Gao Z X,Feng C H,Yang X Z,Huang J G 2012Acta Phys.Sin.61 145201(in Chinese)[高著秀,馮春華,楊宣宗,黃建國2012物理學(xué)報61 145201]

    [4]Ticos C M,Wang Z H,Wurden G A,Kline J L,Montgomery D S,Dorf L A,Shukla P K 2008Phys.Rev.Lett.100 155002

    [5]Asai T,Itagaki H,Numasawa H,Terashima Y,Hirano Y,Hirose A 2010Rev.Sci.Instrum.81 10E119

    [6]Case A,Messer S,Brockington S,Wu L,Witherspoon F D,Elton R 2013Phys.Plasmas20 012704

    [7]Case A,Messer S,Bomgardner R,Witherspoon F D 2010Phys.Plasmas17 053503

    [8]Hsu S C,Merritt E C,Moser A L,Awe T J,Brockington S,Davis J S,Adams C S,Case A,Cassibry J T,Dunn J P,Gilmore M A,Messer S J,Witherspoon F D 2012Phys.Plasmas19 123514

    [9]Ticos C M,Wang Z H,Wurden G A 2011IEEE Trans.Plasma Sci.39 2088

    [10]Messer S,Case A,Bomgardner R,Phillips M,Witherspoon F D 2009Phys.Plasmas16 064502

    [11]Woodru ffS,Hill D N,Stallard B W,Bulmer R,Cohen B,Holcomb C T,Hooper E B,McLean H S,Moller J,Wood R D 2003Phys.Rev.Lett.90 095001

    [12]Kikuchi Y,Nakanishi R,Nakatsuka M,Fukumoto N,Nagataet M 2010IEEE Trans.Plasma Sci.38 232

    [13]Cassibry J T,Thio Y C,Markusic T E,Wu S T 2006J.Propul.Power.22 628

    [14]Cassibry J T 2008IEEE Trans.Plasma Sci.36 2180

    [15]Markusic T E,Thio Y C,Cassibry J T 2002Proceedings of the 38th AIAA Joint Propulsion ConferenceIndianapolis,Indiana,July 7–10,2002,p1

    [16]McNab I R 2009IEEE Trans.Magn.45 381

    [17]Voronin A V,Gusev V K,Petrov Y V,Sakharov N V,Abramova K B,Sklyarova E M,Tolstyakov S Y 2005Nucl Fusion45 1039

    [18]Voronin A V,Gusev V K,Petrov Y V,Mukhin E E,Tolstyakov S Y,Kurskiev G S,Kochergin M M,Hellblom K G 2008Nukleonika53 103

    [19]Ticos C M,Wang Z H,Wurden G A 2008IEEE Trans.Plasma Sci.36 2770

    [20]Wang Z H,Ticos C M,Wurden G A 2007Phys.Plasmas14 103701

    [21]Gao Z X,Huang J G,Han J W,Yang X Z,Feng C H 2010Spacecraft Environment Engineering27 285(in Chinese)[高著秀,黃建國,韓建偉,楊宣宗,馮春華2010航天器環(huán)境工程27 285]

    [22]Han J W,Zhang Z L,Huang J G,Li X Y,Chen Z F,Quan R H,Li H W 2006Spacecraft Environment Engineering23 205(in Chinese)[韓建偉,張振龍,黃建國,李小銀,陳趙峰,全榮輝,李宏偉 2006航天器環(huán)境工程 23 205]

    [23]Yang L,Yan H J,Zhang J L,Hua Y,Ren C S 2014High Voltage Engineering40 13(in Chinese)[楊亮,閆慧杰,張俊龍,滑躍,任春生2014高電壓技術(shù)40 13]

    [24]Wiechula J,Hock C,Iberler M,Manegold T,Sch?nlein A,Jacoby J 2015Phys.Plasmas22 043516

    [25]Wu Y,Bai S B,Wang J Y,Chen J P,Ni X W 2007Optoelectronic Technology27 49(in Chinese)[吳瑩, 白順波,王俊彥,陳建平,倪曉武2007光電子技術(shù)27 49]

    [26]Warner K,Hieftje G M 2002Spectrochim Acta B57 201

    [27]Clements R M 1978J.Vac.Sci.Technol.15 193

    [28]Yuan F Y,Lu W Q,Lin G Q 2009Chinese Journal of Vacuum Science and Technology29 509(in Chinese)[袁方園,陸文琪,林國強(qiáng)2009真空科學(xué)與技術(shù)學(xué)報29 509]

    [29]Li S,Liu Z W,Chen Q,Liu F P,Wang Z D,Yang L Z 2012Spectrosc.Spectr.Anal.32 33(in Chinese)[李森,劉忠偉,陳強(qiáng),劉福平,王正鐸,楊麗珍 2012光譜學(xué)與光譜分析32 33]

    [30]Qian M Y,Ren C S,Wang D Z,Zhang J L,Wei G D 2010J.Appl.Phys.107 063303

    [31]Witherspoon F D,Case A,Messer S J,Bomgardner R,Phillips M W,Brockington S,Elton R 2009Rev.Sci.Instrum.80 083506

    [32]Hart P J 1962Phys.Fluids5 38

    [33]Hart P J 1964J.Appl.Phys.35 3425

    [34]Zhang J L,Yang L,Yan H J,Hua Y,Ren C S 2015Acta Phys.Sin.64 075201(in Chinese)[張俊龍,楊亮,閆慧杰,滑躍,任春生2015物理學(xué)報64 075201]

    [35]Cassibry J T,Thio Y C,Wu S T 2006Phys.Plasmas13 053101

    [36]Huang J G,Han J W,Li H W,Cai M H,Li X Y,Zhang Z L,Chen Z F,Wang L,Yang X Z,Feng C H 2009Chinese Science Bulletin54 150(in Chinese)[黃建國,韓建偉,李宏偉,蔡明輝,李小銀,張振龍,陳趙峰,王龍,楊宣宗,馮春華2009科學(xué)通報54 150]

    [37]Chen Z F 2006M.S.Thesis(Beijing:Center for Space Science and Applied Research,Chinese Academy of Sciences)(in Chinese)[陳趙峰 2006碩士學(xué)位論文 (北京:中國科學(xué)院空間科學(xué)與應(yīng)用研究中心)]

    [38]Kubo H,Kawashima N,Itoh T 1971Plasma Physics13 131

    [39]Qian M Y 2011 Ph.D.Dissertation(Dalian:Dalian University of Technology)(in Chinese)[錢沐楊 2011博士學(xué)位論文(大連:大連理工大學(xué))]

    [40]Wu L F 2010 M.S.Thesis(Wuhan:Wuhan Institute of Technology)(in Chinese)[吳利峰 2010碩士學(xué)位論文(武漢:武漢工程大學(xué))]

    [41]RaizerY P 1991GasDischargePhysics(Berlin:Springer-Verlag)pp58–59

    [42]Roy S,Singh K P,Gaitonde D V 2007Proceedings of the 45th AIAA Aerospace Sciences Meeting and ExhibitReno,Nevada,January 8–11,p1

    [43]Kossyi I A,Kostinsky A Y,Matveyev A A,Silakov V P 1992Plasma Sources Sci.Technol.1 207

    [44]Singh K P,Roya S 2007J.Appl.Phys.101 123308

    PACS:52.80.Vp,94.20.Fg,32.60.+i DOI:10.7498/aps.66.055203

    Experimental study on coaxial gun pulse discharge plasma density change in transport process

    Yang Liang Zhang Jun-Long Yan Hui-Jie Hua Yue Ren Chun-Sheng?

    (Key Laboratory of Materials Modification by Laser,Ion and Electron Beams of the Ministry of Education,School of Physics and
    Optoelectronic Technology,Dalian University of Technology,Dalian 116024,China)

    27 August 2016;revised manuscript

    7 December 2016)

    Coaxial gun discharge plasma with high density and velocity has a number of potential applications in fusion energy,plasma refueling,disruption mitigation in tokamaks,plasma space propulsion,acceleration of dust particles to hypervelocity etc.,and thus has become an important research topic in fields of nuclear physics and aerospace engineering.In this paper,we report the experimental investigation on electrical and transport characteristics of coaxial gun discharge plasma.Based on electrical and optical diagnoses,the discharge voltage,discharge current and axial velocity of plasma transport are measured.Meanwhile,the emission spectrum technology is employed to measure the Stark broadening of Hβspectral line and then plasma density is calculated.The experimental results show that the discharges in the coaxial gun present a feature of multiple discharges and blow-by instability phenomena are observed by photomultiplier acquired signals.In addition,the plasma velocity and density in the transport process are not constant.It is found that the axial plasma velocity in the transport process decreases due to mass gain caused by the snowplow model and the change tendency of plasma density in the transport process is dependent on various settings.A systematic study has been carried out for exploring plasma density change in transport process,and different experimental parameters are adopted in order to further analyze the physical mechanism of plasma density change in transport process.When the air pressure in the coaxial gun is changed from 4.0 Pa to 10 Pa,for 1.08 kJ applied power energy,an obvious difference appears in transport properties of plasma density,i.e.,plasma density increases gradually in 4.0 Pa air while it increasesfirst and then decreases in 10 Pa air.However,the plasma density increases continually in air pressure of 10 Pa when the power energy is increased to 7.68 kJ.Moreover,when the working gas is replaced with argon and discharge setting is 4.0 Pa pressure and 1.08 kJ applied power energy,the plasma density decreases continually in the transport process.The distinct behaviors,as analyzed,are mainly caused by plasma energy transformation difference in the transport process.As it is known,the plasma movement at high velocity in coaxial guns can ionize neutral particles and consume its energy,which results in the increasing plasma density and the decreasing electron and ion temperatures in the transport process.Then,a maximum density is present in the transport process when the electron and ion temperatures are lower than that at which gas ionization occurs.The axial location of maximum density changes with applied power energy,working gas pressure and species,which means that plasma energy transformation and density change properties in transport process strongly rely on different external parameters.The study provides some insight into how to better apply the coaxial gun discharge plasma to practical engineering field.

    coaxial gun,plasma density,Stark broadening

    PACS:52.80.Vp,94.20.Fg,32.60.+i

    10.7498/aps.66.055203

    ?通信作者.E-mail:rchsh@dlut.edu.cn

    ?Corresponding author.E-mail:rchsh@dlut.edu.cn

    猜你喜歡
    氬氣同軸電離
    電離與離子反應(yīng)高考探源
    示范快堆主容器內(nèi)氬氣空間數(shù)值模擬
    同軸單元的先驅(qū)者 Tannoy(天朗)
    企業(yè)車間氬氣泄漏模擬
    水的電離平衡問題解析
    BPRT同軸機(jī)組在高爐上的應(yīng)用設(shè)計
    變速箱輸入軸內(nèi)孔同軸度檢具設(shè)計
    如何復(fù)習(xí)“水的電離”
    歸類總結(jié)促進(jìn)H2O電離的反應(yīng)
    氬氣的純化方法及純化裝置
    低溫與特氣(2014年4期)2014-03-30 02:09:09
    国产亚洲av片在线观看秒播厂| 在线观看人妻少妇| 日本wwww免费看| 国产免费现黄频在线看| 亚洲第一av免费看| 久久久久久久大尺度免费视频| √禁漫天堂资源中文www| 成年美女黄网站色视频大全免费 | 欧美+日韩+精品| 成人影院久久| 日韩制服骚丝袜av| freevideosex欧美| 熟女电影av网| 国产乱人偷精品视频| 午夜福利视频精品| 男女边吃奶边做爰视频| 成人毛片a级毛片在线播放| 午夜福利,免费看| 日韩精品免费视频一区二区三区 | 久久人妻熟女aⅴ| 又黄又爽又刺激的免费视频.| 一级毛片我不卡| av免费在线看不卡| 成年美女黄网站色视频大全免费 | 欧美人与性动交α欧美精品济南到 | 国产亚洲欧美精品永久| 永久免费av网站大全| 午夜影院在线不卡| 美女cb高潮喷水在线观看| 国产高清不卡午夜福利| 91精品国产九色| 日韩一区二区三区影片| 最近的中文字幕免费完整| 色婷婷久久久亚洲欧美| 精品卡一卡二卡四卡免费| 91在线精品国自产拍蜜月| 久久人人爽人人爽人人片va| 少妇被粗大的猛进出69影院 | 欧美人与性动交α欧美精品济南到 | 在线观看免费视频网站a站| 内地一区二区视频在线| 2021少妇久久久久久久久久久| 中文欧美无线码| 欧美三级亚洲精品| 国产高清不卡午夜福利| 欧美+日韩+精品| 成人二区视频| 国产精品久久久久久av不卡| 色吧在线观看| 亚洲国产精品一区二区三区在线| 成人手机av| 亚洲久久久国产精品| 亚洲av不卡在线观看| videossex国产| 人妻夜夜爽99麻豆av| 简卡轻食公司| 亚洲精品日韩在线中文字幕| 日韩一区二区三区影片| 搡女人真爽免费视频火全软件| 久久狼人影院| 成人无遮挡网站| 亚洲欧美成人综合另类久久久| www.av在线官网国产| 美女视频免费永久观看网站| 午夜免费鲁丝| 国产一区二区三区av在线| 青春草视频在线免费观看| 欧美国产精品一级二级三级| 亚洲av二区三区四区| 亚洲不卡免费看| 大片免费播放器 马上看| 欧美精品国产亚洲| 久久99一区二区三区| 国产高清不卡午夜福利| 蜜桃在线观看..| av黄色大香蕉| 久久精品国产鲁丝片午夜精品| 中文字幕亚洲精品专区| 一区二区三区乱码不卡18| 一边亲一边摸免费视频| 最新的欧美精品一区二区| 久久亚洲国产成人精品v| 亚洲少妇的诱惑av| 国产69精品久久久久777片| 春色校园在线视频观看| 校园人妻丝袜中文字幕| 午夜免费鲁丝| 午夜福利视频精品| 久久午夜福利片| 大香蕉97超碰在线| 丝袜脚勾引网站| 99热这里只有是精品在线观看| 少妇猛男粗大的猛烈进出视频| 中国国产av一级| 99九九线精品视频在线观看视频| 午夜免费鲁丝| 少妇猛男粗大的猛烈进出视频| 一区二区三区免费毛片| 久久久久人妻精品一区果冻| 中国国产av一级| 91精品伊人久久大香线蕉| 国精品久久久久久国模美| 久久久精品94久久精品| 多毛熟女@视频| 久久99热这里只频精品6学生| 久久国产亚洲av麻豆专区| 亚洲精品国产av成人精品| 久久久久网色| 校园人妻丝袜中文字幕| 精品久久久久久久久av| 久久婷婷青草| 亚洲av免费高清在线观看| 91在线精品国自产拍蜜月| 精品一区在线观看国产| 最新中文字幕久久久久| 亚洲精品日韩在线中文字幕| 欧美日韩一区二区视频在线观看视频在线| 国产毛片在线视频| 成人综合一区亚洲| 蜜桃久久精品国产亚洲av| 熟女人妻精品中文字幕| 婷婷色综合www| 色婷婷久久久亚洲欧美| 啦啦啦视频在线资源免费观看| 啦啦啦中文免费视频观看日本| 亚洲成人手机| 日本欧美国产在线视频| 久久久久国产精品人妻一区二区| 国产一级毛片在线| 丝瓜视频免费看黄片| 在线观看人妻少妇| 精品视频人人做人人爽| 国产成人精品在线电影| 日本黄色日本黄色录像| 岛国毛片在线播放| 国产成人freesex在线| 午夜福利,免费看| 国产精品.久久久| 亚洲久久久国产精品| 国产白丝娇喘喷水9色精品| 三级国产精品欧美在线观看| 熟女电影av网| 午夜福利在线观看免费完整高清在| 亚洲一区二区三区欧美精品| 免费播放大片免费观看视频在线观看| 99热国产这里只有精品6| 寂寞人妻少妇视频99o| 国产成人精品一,二区| 少妇猛男粗大的猛烈进出视频| 在线观看免费日韩欧美大片 | kizo精华| 男人爽女人下面视频在线观看| 国产伦精品一区二区三区视频9| 国产精品不卡视频一区二区| 最近手机中文字幕大全| 在线观看免费日韩欧美大片 | 99国产综合亚洲精品| 精品人妻熟女av久视频| 极品人妻少妇av视频| 国产精品国产av在线观看| 麻豆精品久久久久久蜜桃| 亚洲婷婷狠狠爱综合网| 啦啦啦在线观看免费高清www| 最近最新中文字幕免费大全7| 色吧在线观看| 又粗又硬又长又爽又黄的视频| 精品久久国产蜜桃| 日本av手机在线免费观看| 一本久久精品| 99热全是精品| 建设人人有责人人尽责人人享有的| 免费日韩欧美在线观看| 国产av国产精品国产| 一本久久精品| 久久久久国产精品人妻一区二区| 亚洲av综合色区一区| 少妇熟女欧美另类| 建设人人有责人人尽责人人享有的| 中文字幕免费在线视频6| 女的被弄到高潮叫床怎么办| 久久精品国产亚洲网站| 在线亚洲精品国产二区图片欧美 | 久久99精品国语久久久| 一区在线观看完整版| 色吧在线观看| 五月伊人婷婷丁香| 蜜臀久久99精品久久宅男| 国产免费一级a男人的天堂| 久久国产精品大桥未久av| 国产在线免费精品| 欧美激情国产日韩精品一区| 97超碰精品成人国产| 日韩欧美一区视频在线观看| 精品亚洲成a人片在线观看| 久久久久久人妻| 欧美日韩亚洲高清精品| 亚洲欧美日韩另类电影网站| 老司机影院成人| 春色校园在线视频观看| videossex国产| 精品少妇内射三级| 国产成人av激情在线播放 | 啦啦啦视频在线资源免费观看| 国产亚洲欧美精品永久| 黄片播放在线免费| 国产一区二区在线观看日韩| 亚洲精品视频女| 热99久久久久精品小说推荐| 久久久久久久大尺度免费视频| 一区二区三区四区激情视频| 久久热精品热| 97超碰精品成人国产| 婷婷色av中文字幕| 人妻少妇偷人精品九色| 一级黄片播放器| 欧美3d第一页| 蜜桃久久精品国产亚洲av| 涩涩av久久男人的天堂| 亚洲三级黄色毛片| 日韩精品免费视频一区二区三区 | 一本—道久久a久久精品蜜桃钙片| 亚洲国产欧美在线一区| 久久久精品94久久精品| 视频中文字幕在线观看| 欧美日韩av久久| 蜜桃国产av成人99| 欧美日韩亚洲高清精品| 日韩一区二区视频免费看| 欧美日韩综合久久久久久| 亚洲精品,欧美精品| 又粗又硬又长又爽又黄的视频| 美女福利国产在线| 国产乱来视频区| 亚洲三级黄色毛片| 伦理电影大哥的女人| 国产精品99久久久久久久久| 国产免费现黄频在线看| 亚洲精品日韩在线中文字幕| 又粗又硬又长又爽又黄的视频| 能在线免费看毛片的网站| 飞空精品影院首页| 午夜免费观看性视频| 不卡视频在线观看欧美| 国产免费一级a男人的天堂| 久久久久久伊人网av| 国产男女内射视频| 久久国产亚洲av麻豆专区| 久久久久网色| 丁香六月天网| 国产成人精品福利久久| 熟妇人妻不卡中文字幕| 久久久精品区二区三区| 女的被弄到高潮叫床怎么办| 伊人久久国产一区二区| 精品人妻在线不人妻| 黄色配什么色好看| 成年美女黄网站色视频大全免费 | 日日爽夜夜爽网站| 一本一本综合久久| a级毛色黄片| 国产色爽女视频免费观看| 国产午夜精品一二区理论片| 精品酒店卫生间| av播播在线观看一区| 一级片'在线观看视频| 亚洲欧美清纯卡通| 99热这里只有精品一区| 99国产精品免费福利视频| 汤姆久久久久久久影院中文字幕| 亚洲欧洲日产国产| 男人爽女人下面视频在线观看| 国产日韩欧美在线精品| av专区在线播放| 久久鲁丝午夜福利片| av播播在线观看一区| 激情五月婷婷亚洲| 欧美xxxx性猛交bbbb| 青春草视频在线免费观看| 午夜福利,免费看| 黑丝袜美女国产一区| 少妇被粗大猛烈的视频| 自拍欧美九色日韩亚洲蝌蚪91| 制服人妻中文乱码| 王馨瑶露胸无遮挡在线观看| 黑人巨大精品欧美一区二区蜜桃 | 成人午夜精彩视频在线观看| 色网站视频免费| 日韩人妻高清精品专区| 少妇被粗大猛烈的视频| 丰满乱子伦码专区| 啦啦啦视频在线资源免费观看| 大香蕉97超碰在线| 亚洲国产欧美日韩在线播放| 精品一区二区三卡| 五月开心婷婷网| 满18在线观看网站| 国产熟女午夜一区二区三区 | 国产一级毛片在线| 亚洲三级黄色毛片| 女人久久www免费人成看片| 久久久久久久久久久久大奶| 亚洲精品日韩av片在线观看| 亚洲av福利一区| 国精品久久久久久国模美| 黄色视频在线播放观看不卡| 丝瓜视频免费看黄片| 熟女av电影| 天堂中文最新版在线下载| 日本午夜av视频| 欧美国产精品一级二级三级| 日韩欧美精品免费久久| 免费看av在线观看网站| 美女中出高潮动态图| 色网站视频免费| 国产精品偷伦视频观看了| 自拍欧美九色日韩亚洲蝌蚪91| 少妇猛男粗大的猛烈进出视频| 久久久精品区二区三区| 日本欧美视频一区| 熟女电影av网| 日韩中文字幕视频在线看片| av不卡在线播放| 91成人精品电影| 两个人的视频大全免费| 夫妻午夜视频| 少妇 在线观看| 麻豆成人av视频| a级毛片在线看网站| 国产一区亚洲一区在线观看| 日韩成人av中文字幕在线观看| 免费观看在线日韩| 欧美激情国产日韩精品一区| 亚洲成色77777| 久久久午夜欧美精品| 国产精品国产三级专区第一集| 亚洲国产精品专区欧美| xxx大片免费视频| 人人妻人人添人人爽欧美一区卜| 各种免费的搞黄视频| 久久精品久久精品一区二区三区| 少妇 在线观看| 国产爽快片一区二区三区| 日韩不卡一区二区三区视频在线| 国产一区二区三区av在线| 国产亚洲欧美精品永久| 观看美女的网站| 亚洲婷婷狠狠爱综合网| 自线自在国产av| 国产成人精品福利久久| 亚洲av男天堂| 国产极品粉嫩免费观看在线 | 国产在线免费精品| 亚洲一区二区三区欧美精品| 国产欧美另类精品又又久久亚洲欧美| 国产片特级美女逼逼视频| 国产色婷婷99| 国产免费一级a男人的天堂| 亚洲综合色网址| 99久久精品国产国产毛片| 亚洲丝袜综合中文字幕| 亚洲精品乱码久久久久久按摩| 精品久久国产蜜桃| 在线免费观看不下载黄p国产| 久久久欧美国产精品| 蜜臀久久99精品久久宅男| 成人18禁高潮啪啪吃奶动态图 | 18禁观看日本| 人体艺术视频欧美日本| 久久久精品94久久精品| 国产精品久久久久久久久免| 九草在线视频观看| 色婷婷av一区二区三区视频| 国产成人a∨麻豆精品| 美女视频免费永久观看网站| 看十八女毛片水多多多| 亚洲av.av天堂| 最近手机中文字幕大全| 十八禁高潮呻吟视频| 欧美日韩一区二区视频在线观看视频在线| 色吧在线观看| 国产精品嫩草影院av在线观看| 天天影视国产精品| 日本av免费视频播放| 国产精品久久久久成人av| 久久久欧美国产精品| 91aial.com中文字幕在线观看| 91久久精品国产一区二区成人| 飞空精品影院首页| 黄色毛片三级朝国网站| 日韩精品免费视频一区二区三区 | 精品卡一卡二卡四卡免费| 婷婷色av中文字幕| 久久久久国产精品人妻一区二区| 一区二区日韩欧美中文字幕 | 熟女人妻精品中文字幕| av视频免费观看在线观看| 亚洲一区二区三区欧美精品| 91久久精品国产一区二区三区| 99热这里只有是精品在线观看| 精品人妻熟女毛片av久久网站| 如日韩欧美国产精品一区二区三区 | 亚洲少妇的诱惑av| 中文字幕免费在线视频6| 欧美+日韩+精品| 男女免费视频国产| 国产熟女欧美一区二区| 欧美最新免费一区二区三区| 精品午夜福利在线看| 91精品国产九色| 久久这里有精品视频免费| 日韩精品免费视频一区二区三区 | 人妻系列 视频| 亚洲精华国产精华液的使用体验| 午夜日本视频在线| 黑人猛操日本美女一级片| 国产精品久久久久久精品电影小说| 国产免费又黄又爽又色| 男女免费视频国产| 99热国产这里只有精品6| 国产精品一区二区在线不卡| 亚洲av综合色区一区| 99视频精品全部免费 在线| 国产片内射在线| 在线观看人妻少妇| 大话2 男鬼变身卡| videossex国产| 国产成人午夜福利电影在线观看| 日韩大片免费观看网站| 成年av动漫网址| 一区二区三区精品91| av电影中文网址| 亚洲精品美女久久av网站| 大码成人一级视频| 一边摸一边做爽爽视频免费| 国产精品一区二区在线观看99| 国产国语露脸激情在线看| 久久久久久久精品精品| 亚洲国产欧美日韩在线播放| 国产在线免费精品| 26uuu在线亚洲综合色| 久久毛片免费看一区二区三区| 国产精品久久久久久久电影| 免费久久久久久久精品成人欧美视频 | 精品视频人人做人人爽| 日本午夜av视频| 亚洲婷婷狠狠爱综合网| 久久婷婷青草| 日韩中字成人| 最近手机中文字幕大全| 久久精品熟女亚洲av麻豆精品| 一本一本综合久久| 大香蕉97超碰在线| 日韩电影二区| 国产精品成人在线| 国产免费一区二区三区四区乱码| 在线观看人妻少妇| 亚洲美女黄色视频免费看| 午夜福利视频精品| 亚洲精品久久成人aⅴ小说 | 丰满乱子伦码专区| videos熟女内射| 蜜臀久久99精品久久宅男| 在线观看免费高清a一片| 国产午夜精品一二区理论片| 又大又黄又爽视频免费| 亚洲欧美中文字幕日韩二区| 午夜免费观看性视频| 男人操女人黄网站| 嫩草影院入口| 国产乱人偷精品视频| 国产精品成人在线| 亚洲高清免费不卡视频| 亚洲国产日韩一区二区| 色哟哟·www| 久久久a久久爽久久v久久| 亚洲色图 男人天堂 中文字幕 | 中文天堂在线官网| 一级毛片电影观看| 99久久人妻综合| 99视频精品全部免费 在线| 自拍欧美九色日韩亚洲蝌蚪91| 精品一区二区三区视频在线| 亚洲在久久综合| 美女国产视频在线观看| 在线 av 中文字幕| 欧美激情国产日韩精品一区| 国产成人免费观看mmmm| 99热全是精品| 99精国产麻豆久久婷婷| 亚洲情色 制服丝袜| 国产乱人偷精品视频| 亚洲精品久久久久久婷婷小说| 王馨瑶露胸无遮挡在线观看| 亚洲欧美成人精品一区二区| 亚洲国产精品专区欧美| 性色av一级| 满18在线观看网站| 精品少妇内射三级| 成人午夜精彩视频在线观看| 免费播放大片免费观看视频在线观看| 国产成人aa在线观看| 五月开心婷婷网| xxxhd国产人妻xxx| 欧美97在线视频| 亚洲美女视频黄频| 建设人人有责人人尽责人人享有的| 日韩成人伦理影院| 精品一区二区免费观看| 欧美精品亚洲一区二区| 看免费成人av毛片| 日韩不卡一区二区三区视频在线| 在线观看美女被高潮喷水网站| 熟女人妻精品中文字幕| 如何舔出高潮| 久久亚洲国产成人精品v| 亚洲精品aⅴ在线观看| 99久久中文字幕三级久久日本| 丰满迷人的少妇在线观看| 精品一区二区三区视频在线| 乱码一卡2卡4卡精品| 性高湖久久久久久久久免费观看| 大香蕉久久成人网| 男女国产视频网站| 国产 精品1| 欧美97在线视频| 99国产综合亚洲精品| 男女啪啪激烈高潮av片| 最近最新中文字幕免费大全7| 国产国拍精品亚洲av在线观看| 天美传媒精品一区二区| 最近的中文字幕免费完整| 成人综合一区亚洲| 在线看a的网站| 亚洲国产欧美在线一区| 老司机亚洲免费影院| 999精品在线视频| 黄色一级大片看看| 七月丁香在线播放| 国产精品久久久久成人av| 国产精品国产av在线观看| 亚洲一级一片aⅴ在线观看| 纯流量卡能插随身wifi吗| 国产黄频视频在线观看| 人妻系列 视频| 亚洲精品视频女| 老司机影院成人| 天堂中文最新版在线下载| 午夜91福利影院| 人妻制服诱惑在线中文字幕| 制服人妻中文乱码| 日韩精品免费视频一区二区三区 | 婷婷色av中文字幕| 肉色欧美久久久久久久蜜桃| 国产av国产精品国产| 国产探花极品一区二区| 欧美性感艳星| 国产精品国产av在线观看| 免费av中文字幕在线| 成人综合一区亚洲| 五月玫瑰六月丁香| 欧美丝袜亚洲另类| 热re99久久精品国产66热6| 视频在线观看一区二区三区| 少妇 在线观看| 边亲边吃奶的免费视频| 国内精品宾馆在线| 天堂中文最新版在线下载| 国产在视频线精品| 老司机影院毛片| 能在线免费看毛片的网站| 亚洲av二区三区四区| 久久精品久久精品一区二区三区| 午夜福利,免费看| av天堂久久9| 久久av网站| 久久青草综合色| 十八禁网站网址无遮挡| 中国国产av一级| 热99久久久久精品小说推荐| 男女高潮啪啪啪动态图| 赤兔流量卡办理| 黑人猛操日本美女一级片| 赤兔流量卡办理| 久久影院123| 国产伦理片在线播放av一区| 多毛熟女@视频| 免费黄频网站在线观看国产| 又大又黄又爽视频免费| 亚洲成人一二三区av| 亚洲精品自拍成人| 色婷婷久久久亚洲欧美| 狠狠婷婷综合久久久久久88av| 国产av国产精品国产| 人人妻人人爽人人添夜夜欢视频| 在线观看美女被高潮喷水网站| 午夜视频国产福利| 亚洲国产欧美在线一区| 欧美97在线视频| 日韩大片免费观看网站| 乱人伦中国视频| 一区二区av电影网| 国产一级毛片在线| 久久99蜜桃精品久久| 欧美3d第一页| 如日韩欧美国产精品一区二区三区 | 黄色怎么调成土黄色| www.av在线官网国产| 久久韩国三级中文字幕| 色网站视频免费| 插阴视频在线观看视频| 美女脱内裤让男人舔精品视频| 免费观看的影片在线观看| 午夜福利网站1000一区二区三区| 在线 av 中文字幕| 男女边摸边吃奶|