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

    亞大氣壓六相交流電弧等離子體射流特性研究:實(shí)驗(yàn)測(cè)量?

    2018-03-26 22:07:00郭恒蘇運(yùn)波李和平曾實(shí)聶秋月李占賢李志輝
    物理學(xué)報(bào) 2018年4期
    關(guān)鍵詞:空腔電弧射流

    郭恒 蘇運(yùn)波 李和平 曾實(shí) 聶秋月 李占賢 李志輝

    1)(清華大學(xué)工程物理系,北京 100084)

    2)(華北理工大學(xué)機(jī)械工程學(xué)院,唐山 063500)

    3)(哈爾濱工業(yè)大學(xué)電氣工程及自動(dòng)化學(xué)院,哈爾濱 150001)

    4)(中國(guó)空氣動(dòng)力研究與發(fā)展中心超高速所,綿陽(yáng) 621000)

    5)(國(guó)家計(jì)算流體力學(xué)實(shí)驗(yàn)室,北京 100191)

    1 引 言

    距海平面20—100 km的臨近空間區(qū)域是人類航空與航天活動(dòng)的過(guò)渡區(qū)域,該區(qū)域內(nèi)相對(duì)稀薄的大氣環(huán)境有效減小了飛行器所受的空氣阻力,使得飛行器容易實(shí)現(xiàn)高超聲速飛行,具有廣闊的商業(yè)和軍事應(yīng)用前景.但同時(shí),臨近空間飛行器在高超聲速飛行過(guò)程中也面臨著復(fù)雜介質(zhì)環(huán)境的強(qiáng)烈影響,特別是其表面激波熱致高密度、強(qiáng)碰撞、非均勻和非平衡(包括熱力學(xué)和化學(xué)非平衡效應(yīng))等離子體的產(chǎn)生[1,2].一方面,能量的注入和等離子體的存在使得飛行器周圍局部區(qū)域的當(dāng)?shù)芈曀僭黾?邊界層特性(如局部流場(chǎng)和壓力場(chǎng)的空間分布、當(dāng)?shù)伛R赫數(shù)的變化以及邊界層分離特性等)與常溫空氣相比將發(fā)生顯著變化,從而導(dǎo)致飛行器飛行阻力的變化[3];另一方面,包覆飛行器表面的等離子體鞘套環(huán)境對(duì)電磁波的吸收、散射和折射,會(huì)引發(fā)“黑障”效應(yīng)[4?6],從而影響飛行器的通信和控制系統(tǒng).因此,開(kāi)展高超聲速飛行器與其周圍復(fù)雜介質(zhì)環(huán)境間相互作用機(jī)理的研究對(duì)于指導(dǎo)飛行器的設(shè)計(jì)具有重要意義.

    目前,有關(guān)飛行器特性的研究主要采用飛行試驗(yàn)、地面風(fēng)洞實(shí)驗(yàn)和數(shù)值模擬方法.飛行試驗(yàn)是直接對(duì)飛行器在空間飛行過(guò)程中的關(guān)鍵參數(shù)進(jìn)行在線診斷,包括美國(guó)的無(wú)線傳輸衰減測(cè)量(Radio At-tenuation Measurements,RAM)、日本的超高速飛行試驗(yàn)(Hypersonic Flight Experiment,HYFLEX)以及俄羅斯(前蘇聯(lián))的載人航空航天器(Manned Aero-Space Vehicle,VKA)等飛行試驗(yàn)項(xiàng)目[7?11].但這種方法代價(jià)高昂,且在這種極端條件下的實(shí)驗(yàn)診斷手段也比較匱乏,一般通過(guò)探針或者電磁波測(cè)量飛行器壁面附近等離子體鞘套的電子數(shù)密度.因此,測(cè)量數(shù)據(jù),特別是公開(kāi)發(fā)表的實(shí)驗(yàn)數(shù)據(jù)十分有限.自20世紀(jì)60年代以來(lái),地面模擬就成了研究空間飛行器與其周圍環(huán)境相互作用過(guò)程的重要手段之一.其中,用來(lái)產(chǎn)生復(fù)雜介質(zhì)環(huán)境的模擬設(shè)備主要為風(fēng)洞(如等離子體風(fēng)洞、激波風(fēng)洞等)和能夠引出高速、強(qiáng)束流的等離子體發(fā)生器[12?14].此類大型地面模擬裝置的建設(shè)和維護(hù)成本,特別是大型風(fēng)洞的建設(shè)費(fèi)用和運(yùn)行成本也同樣非常昂貴,開(kāi)展大量重復(fù)性的實(shí)驗(yàn)仍然具有一定的困難.近年來(lái),隨著氣體放電理論和等離子體源技術(shù)的不斷發(fā)展,涌現(xiàn)出了不少小規(guī)模實(shí)驗(yàn)室模擬裝置,比如螺旋波等離子體源、電弧等離子體源、輝光放電、電感/電容耦合放電等離子體源等模擬裝置[15?18].盡管這些等離子體源的產(chǎn)生方式和產(chǎn)生條件與飛行器實(shí)際飛行過(guò)程中強(qiáng)激波所產(chǎn)生的等離子體環(huán)境并不完全相同,但從模擬飛行器與其周圍復(fù)雜介質(zhì)環(huán)境在某些方面的相互作用過(guò)程,如“黑障”中等離子體-波相互作用過(guò)程的研究依然是十分有意義的.

    從臨近空間飛行器“黑障”問(wèn)題研究的需求來(lái)看,需要在地面模擬實(shí)驗(yàn)中產(chǎn)生具有足夠大體積(如有效試驗(yàn)區(qū)長(zhǎng)度L3=30—50 cm)、足夠厚(如鞘套厚度δ=5—10 cm)以及足夠高的等離子體密度(如電子數(shù)密度ne=108—1012cm?3)的包覆飛行器表面的等離子體鞘套環(huán)境(如圖1所示),進(jìn)而能夠研究等離子體鞘套特性及其與入射電磁波間的相互作用規(guī)律.在大氣壓或亞大氣壓條件下,盡管熱等離子體的電子數(shù)密度相對(duì)輝光放電要高得多,但產(chǎn)生大體積、穩(wěn)定的電弧等離子體通常十分困難,例如在等離子體材料處理應(yīng)用中所產(chǎn)生的直流電弧放電等離子體射流盡管在層流條件下其射流的長(zhǎng)度可以達(dá)到1 m量級(jí)[19],但其徑向尺寸通常在1 cm量級(jí)或更小[20?23];采用電感耦合方式所產(chǎn)生的放電等離子體的直徑盡管相對(duì)直流放電要大得多,但也往往限于幾個(gè)厘米之內(nèi)[20,24],很難產(chǎn)生更大直徑的等離子體射流.日本的Watanabe教授研究組[25,26]在進(jìn)行玻璃粉末的熱等離子體球化處理時(shí)提出了采用多相交流電弧放電裝置產(chǎn)生直徑可達(dá)10 cm、厚度可達(dá)15 cm的片狀高溫等離子體放電區(qū).這種放電裝置從電極結(jié)構(gòu)的設(shè)計(jì)上解決了在大氣壓或亞大氣壓條件下產(chǎn)生大直徑電弧放電等離子體的問(wèn)題,但依然無(wú)法產(chǎn)生用于模擬包覆飛行器表面的、具有足夠長(zhǎng)度的大體積等離子體環(huán)境.在本研究中,我們將雙射流直流電弧放電等離子體發(fā)生器的設(shè)計(jì)理念[27]與上述多相交流電弧放電等離子體發(fā)生器的設(shè)計(jì)理念[25,26]相結(jié)合,建立了能夠在大氣壓和亞大氣壓條件下放電的六相交流電弧等離子體實(shí)驗(yàn)平臺(tái),產(chǎn)生了大體積等離子體電弧射流,并研究了不同工作參數(shù)對(duì)等離子體射流特性的影響規(guī)律,為下一步開(kāi)展臨近空間高超聲速飛行器以及航天器再入大氣環(huán)境飛行過(guò)程中“黑障”問(wèn)題的研究準(zhǔn)備了良好的條件.

    圖1 包覆飛行器的等離子體鞘套示意圖Fig.1.Schematic of the plasma layer surrounding an aerocraft.

    2 實(shí)驗(yàn)平臺(tái)介紹

    2.1 物理設(shè)計(jì)的基本思路

    在大氣壓及亞大氣壓條件下產(chǎn)生包覆飛行器的大體積等離子體鞘套,必須產(chǎn)生大體積的非平衡態(tài)等離子體.而在中高氣壓條件下,由于氣體的臨界擊穿場(chǎng)強(qiáng)較大(如在大氣壓條件下該值通常在105—106V/m量級(jí)[28]),等離子體發(fā)生器的特征尺度(此處以發(fā)生器電極間距表征)較小(通常在毫米—厘米量級(jí)),導(dǎo)致很難產(chǎn)生大體積(如本研究中所需要的直徑在10—20 cm、長(zhǎng)度在30—50 cm以上)的非平衡等離子體射流.研究中我們提出以文獻(xiàn)[25,26]所報(bào)道的多相交流電弧放電等離子體源為原型(如圖2(a)所示),結(jié)合雙射流直流電弧等離子體發(fā)生器的設(shè)計(jì)理念(如圖2(b)所示)[27],通過(guò)增加電極區(qū)水冷約束管和電極上游冷氣體注入噴嘴(如圖2(c)所示)獲得大體積等離子體射流.

    圖2 新型六相交流電弧放電等離子體射流源設(shè)計(jì)思路示意圖 (a)文獻(xiàn)[25,26]中的電弧等離子體源;(b)多射流直流電弧等離子體源;(c)本文提出的新型電弧等離子體源Fig.2.Schematics for the design of the new type sixphase arc discharge plasma jet source:(a)Plasma source presented in Refs.[25,26];(b)multi-jet direct current arc plasma source;(c)new type plasma source proposed in this study.

    具體的物理設(shè)計(jì)思路是如下.

    1)采用文獻(xiàn)[25,26]所報(bào)道的多相交流電弧放電等離子體源獲得較大的放電區(qū)直徑.當(dāng)采用如圖2(a)所示的多電極、多相放電發(fā)生器結(jié)構(gòu)設(shè)計(jì)時(shí),一方面,在保持相鄰電極間的距離為厘米量級(jí)的同時(shí),可以使得所有電極所形成的放電區(qū)直徑隨著電極數(shù)目的增加而顯著增加,通??梢赃_(dá)到數(shù)十厘米量級(jí);而另一方面,由于采用了多相交流電源驅(qū)動(dòng)放電,不僅避免了單相交流放電時(shí)電流過(guò)零所造成的熄弧問(wèn)題,而且可以通過(guò)設(shè)計(jì)不同的電極放電時(shí)序來(lái)調(diào)整放電所形成的高溫區(qū)位置[26].如圖2(a)所示,當(dāng)采用相鄰電極順序放電的連接方式時(shí),將形成類似環(huán)狀的等離子體高溫區(qū),這有利于在鈍體表面形成具有一定厚度的等離子體鞘套.

    2)盡管采用如圖2(a)所示的等離子體發(fā)生器結(jié)構(gòu)設(shè)計(jì)可以獲得具有較大放電區(qū)直徑的等離子體,但弧柱區(qū)的高溫部分電離氣體在發(fā)生器幾何軸線方向均勻膨脹,無(wú)法獲得長(zhǎng)等離子體射流,所形成的等離子體區(qū)厚度較小.為此,Yao等[25]提出將12個(gè)電極分成兩組,分兩層進(jìn)行布置,通過(guò)調(diào)整層間距(給定的層間距為15 cm)從而在放電區(qū)獲得具有一定厚度的等離子體,但該方法依然無(wú)法在電極區(qū)下游的射流區(qū)形成長(zhǎng)而穩(wěn)定的等離子體射流.為此,我們借鑒雙射流直流電弧等離子體發(fā)生器的設(shè)計(jì)思路,即將電極幾何軸線與發(fā)生器幾何軸線間的夾角由圖2(a)的90°變?yōu)?0°,從而通過(guò)電極尖下游區(qū)等離子體射流的相互作用形成具有一定長(zhǎng)度的、總體上沿發(fā)生器幾何軸線流動(dòng)的高溫等離子體射流(如圖2(b)所示)[29].

    3)為了進(jìn)一步約束放電區(qū)高溫部分電離氣體的自由膨脹,我們?cè)趫D2(b)的基礎(chǔ)上增加水冷約束管和電極平面上游的冷氣體注入噴嘴.如圖2(c)所示,增加上述部件后,不僅約束了高溫電弧等離子體的自由膨脹,而且通過(guò)電極平面上游冷氣體的注入,使得更多的冷氣體被電弧加熱后沿發(fā)生器幾何軸線向電極平面的下游流動(dòng),從而在保證放電區(qū)直徑基本不變的情況下有利于形成更長(zhǎng)的等離子體射流[27].

    2.2 六相交流電弧放電等離子體實(shí)驗(yàn)平臺(tái)

    圖3 MPX-2015 (a)示意圖;(b)實(shí)物照片F(xiàn)ig.3.Schematic(a)and picture(b)of MPX-2015.

    按照上述設(shè)計(jì)思路,我們?cè)O(shè)計(jì)了多相交流電弧等離子體實(shí)驗(yàn)平臺(tái)(multiphase gas discharge plasma experimental platform-2015,MPX-2015).圖3(a)和圖3(b)分別為該實(shí)驗(yàn)平臺(tái)的示意圖和實(shí)物圖,包括真空腔體及其內(nèi)部的等離子體發(fā)生器、電源系統(tǒng)和控制系統(tǒng).實(shí)驗(yàn)平臺(tái)總長(zhǎng)1.7 m,總高約2 m.真空腔由兩段組成,其中前段真空腔側(cè)壁安裝有6個(gè)伺服電動(dòng)缸,每個(gè)電動(dòng)缸帶動(dòng)一個(gè)電極做直線運(yùn)動(dòng);前段真空腔的內(nèi)部安裝有一個(gè)內(nèi)徑為100 mm、采用水冷結(jié)構(gòu)的約束管,主要起到對(duì)等離子體射流的約束作用.該約束管為不銹鋼材料,亦可作為與6個(gè)電極相對(duì)應(yīng)的地電極.在前段真空腔的端部中心位置安裝有1個(gè)伺服電動(dòng)缸,用以帶動(dòng)冷氣體進(jìn)氣噴嘴做直線運(yùn)動(dòng).后段真空腔采用雙層水冷結(jié)構(gòu),在其側(cè)壁上開(kāi)有2個(gè)長(zhǎng)方形觀察窗和2個(gè)圓形觀察窗,以及1個(gè)連接真空泵的法蘭.2個(gè)長(zhǎng)方形觀察窗的長(zhǎng)度均為380 mm,靠近真空泵法蘭的長(zhǎng)方形觀察窗寬70 mm,而該窗口對(duì)面的長(zhǎng)方形觀察窗的寬度則為130 mm;兩個(gè)圓形觀察窗直徑均為150 mm.另外,真空腔底部觀察窗也可以更換為安裝用于固定鈍體的軌道臺(tái)架.

    MPX-2015的電源系統(tǒng)主要由6個(gè)弧焊電源和1個(gè)6相隔離變壓器組成,同時(shí)也包括了用于非接觸引弧的高頻引弧器和用于抑制變壓器涌流的軟啟動(dòng)電路.該電源系統(tǒng)以六個(gè)單相交流弧焊變壓器為主功率源,采用六相隔離變壓器將三相交流電變換成為具有60°相位差的六相交流電,為六個(gè)單相交流弧焊電源供電,并采用如圖4所示的線路連接方式將等離子體發(fā)生器的6個(gè)電極與6臺(tái)單相交流弧焊機(jī)連接.放電過(guò)程中,相鄰兩個(gè)電極間所加交流電壓的相位差為60°;每相電流的波形均為標(biāo)準(zhǔn)正弦波,當(dāng)其中1相電流過(guò)零時(shí),其他5相的電流并不會(huì)同時(shí)經(jīng)過(guò)零點(diǎn).于是,盡管當(dāng)電流過(guò)零時(shí)電弧會(huì)熄滅,但之后電流會(huì)持續(xù)增大,并依靠其他5相完成自燃,從而保證電弧的持續(xù)和一定程度上的穩(wěn)定.由于在放電過(guò)程中沒(méi)有中心點(diǎn),因此,每一電極在充當(dāng)功率電極的同時(shí)也在輪流地充當(dāng)?shù)仉姌O.在發(fā)生器穩(wěn)定放電的情況下,每相的電流在70—110 A的范圍內(nèi)變化,兩相之間的工作電壓則維持在24 V左右.

    MPX-2015的控制系統(tǒng)由三菱Q系列可編程邏輯控制器(PLC)、三菱伺服電機(jī)、觸摸屏和上位機(jī)組成.PLC和觸摸屏之間采用RS422數(shù)據(jù)傳輸協(xié)議通訊,PLC和上位機(jī)之間則采用RS232數(shù)據(jù)傳輸協(xié)議通訊.在放電過(guò)程中,觸摸屏通過(guò)PLC主要進(jìn)行系統(tǒng)的啟動(dòng)和監(jiān)控,7個(gè)伺服電機(jī)分別帶動(dòng)等離子體發(fā)生器的6個(gè)電極和1個(gè)冷氣體進(jìn)氣噴嘴運(yùn)動(dòng),并采用2個(gè)QD75D4運(yùn)動(dòng)控制模塊對(duì)7個(gè)電機(jī)進(jìn)行定位運(yùn)動(dòng)控制.上位機(jī)則使用LabVIEW程序,主要用于后續(xù)實(shí)驗(yàn)數(shù)據(jù)的采集.

    2.3 等離子體射流關(guān)鍵參數(shù)的實(shí)驗(yàn)測(cè)量方法概述

    我們采用發(fā)射光譜法來(lái)測(cè)量等離子體射流區(qū)的氣體溫度(Th)、電子激發(fā)溫度(Te)和電子數(shù)密度(ne)這三個(gè)關(guān)鍵參數(shù).其中,采用Boltzmann圖法測(cè)量了等離子體的電子激發(fā)溫度(Te)[30];采用OH(羥基)基團(tuán)在306—310 nm間的發(fā)射光譜測(cè)量了氣體的轉(zhuǎn)動(dòng)溫度(Trot)[31],并假定Th≈Trot;采用等離子體在380—430 nm間的連續(xù)輻射譜測(cè)量了電子數(shù)密度(ne)[32].根據(jù)第1節(jié)所述,包覆飛行器表面的等離子體鞘套環(huán)境需要有足夠高的電子密度(ne=108—1012cm?3).因此,為了定量確定等離子體射流的尺寸,包括最大射流直徑(d)和射流長(zhǎng)度(l),我們以等離子體電子數(shù)密度1010cm?3為界限(記為ne,crt),待每次放電穩(wěn)定之后,通過(guò)對(duì)等離子體射流進(jìn)行徑向和軸向二維掃描測(cè)量,得到所測(cè)電子數(shù)密度大于ne,crt的等離子體射流區(qū)尺寸.

    圖4 MPX-2015電源系統(tǒng)電路連接圖Fig.4.Circuit for the MPX-2015.

    3 等離子體自由射流特性的實(shí)驗(yàn)研究

    3.1 典型實(shí)驗(yàn)工況下的等離子體自由射流特性

    圖5給出了等離子體自由射流形貌(曝光時(shí)間texp=50 ms),實(shí)驗(yàn)條件為:氬-氮混合氣體(氮含量10%)作為等離子體工作氣體,總流量(Qtot)為16.0 slpm(其中從每個(gè)電極環(huán)縫注入的氣體流量(Qe)為1.0 slpm(標(biāo)準(zhǔn)升每分鐘),電極平面上游的冷氣體注入量(Qc)為10.0 slpm)、單相弧電流峰值(I)為80 A,相鄰兩電極間距(le)為40 cm,水冷約束管出口內(nèi)徑(r1)為2.5 cm,真空腔壓強(qiáng)(p)為1000 Pa.為了便于后續(xù)分析,以水冷約束管出口中心作為坐標(biāo)原點(diǎn)(O),以發(fā)生器幾何軸線作為z軸,并將等離子體射流的流動(dòng)方向定義為z軸正方向,與之相垂直的方向定義為徑向(r軸).從圖5可以看到,從水冷約束管噴出的等離子體射流沿流動(dòng)方向逐漸擴(kuò)張.圖6則給出了相同工況下等離子體發(fā)生器軸線上z=10.0 cm處的發(fā)射光譜.利用這些譜線,可以得到等離子體射流關(guān)鍵參數(shù)的空間分布規(guī)律.

    圖5 典型工況下的等離子體自由射流照片F(xiàn)ig.5.Image of the plasma free jet under a typical operating condition.

    圖6 發(fā)生器軸線上z=10.0 cm處的等離子體發(fā)射光譜圖Fig.6.Emission spectrum of the plasma jet at the location z=10.0 cm on the jet axis.

    圖7給出了相同工況下實(shí)驗(yàn)測(cè)量與采用非平衡等離子體模型計(jì)算[33]得到的等離子體自由射流區(qū)氣體溫度、電子激發(fā)溫度和電子數(shù)密度沿發(fā)生器幾何軸線的變化規(guī)律,其中有關(guān)非平衡態(tài)等離子體數(shù)值模擬物理數(shù)學(xué)模型的具體介紹可參見(jiàn)文獻(xiàn)[33].可以看到:1)上述三個(gè)特征參數(shù)(Th,Te和ne)均沿射流流動(dòng)方向在z=0—50.0 cm的范圍內(nèi)逐漸降低,實(shí)驗(yàn)測(cè)量得到的電子激發(fā)溫度和重粒子溫度分別從水冷約束管出口處的7601 K和3600 K降低到了3037 K和1100 K,而電子數(shù)密度則從出口處的1.25×1015cm?3下降到了1.1×1011cm?3;數(shù)值模擬結(jié)果與實(shí)驗(yàn)測(cè)量結(jié)果在定性的變化趨勢(shì)上是完全一致的;在z=50.0 cm處,實(shí)驗(yàn)測(cè)量得到電子溫度和重粒子溫度較數(shù)值模擬結(jié)果分別低22%和26%,實(shí)驗(yàn)測(cè)量得到的電子數(shù)密度與計(jì)算結(jié)果相差約18%,這在一定程度上驗(yàn)證了本文實(shí)驗(yàn)測(cè)量結(jié)果的可靠性;2)相應(yīng)地,等離子體射流區(qū)電子數(shù)密度大于ne,crt(=1010cm?3)的射流直徑亦從水冷約束管出口處的6.0 cm擴(kuò)大到了11.6 cm(如圖8所示).

    圖7 等離子體射流區(qū)氣體溫度、電子激發(fā)溫度和數(shù)密度沿軸線的變化規(guī)律Fig.7.Pro files of Th,Teand nealong the geometrical axis of the plasma jet.

    圖8 等離子體射流直徑沿軸線的變化規(guī)律Fig.8.Variations of the diameter of the plasma jet along the axis.

    3.2 主要工作參數(shù)對(duì)等離子體自由射流特性的影響

    為了進(jìn)一步研究不同放電參數(shù)(如真空腔壓強(qiáng)、電極間距、弧電流、等離子體工作氣體流量等)對(duì)自由射流特性的影響規(guī)律,針對(duì)圖5所對(duì)應(yīng)的典型工況設(shè)計(jì)了如表1所列的16組實(shí)驗(yàn)條件,其中水冷約束管內(nèi)徑以及從電極平面上游注入的冷氣體量均保持不變,即r1=2.5 cm,Qc=10.0 slpm.

    表1 自由射流特性研究實(shí)驗(yàn)工況表Table 1.Operating conditions for the plasma free jets.

    圖9和圖10分別給出了對(duì)應(yīng)于表1不同放電工況下的等離子體射流長(zhǎng)度和最大直徑.其中,從圖5中O點(diǎn)開(kāi)始,軸向方向每隔5.0 cm,徑向方向每隔1.0 cm進(jìn)行等離子體射流光譜空間分布的掃描測(cè)量,射流長(zhǎng)度和直徑均以電子密度約為ne,crt為邊界得到的最大直徑以及長(zhǎng)度定義.實(shí)驗(yàn)測(cè)量結(jié)果表明:1)真空腔壓強(qiáng)(p)會(huì)對(duì)等離子體射流尺度產(chǎn)生顯著的影響,隨著壓強(qiáng)的降低,等離子體射流長(zhǎng)度與直徑均有所增加;2)在弧電流(I)保持不變的情況下,相鄰兩電極間的距離(le)對(duì)放電穩(wěn)定性有明顯的影響,隨著弧電流的降低,放電穩(wěn)定性變差;且真空腔壓強(qiáng)越低,這種變化越明顯;3)其他條件保持不變的情況下,增加弧電流(I)和沿電極環(huán)縫注入的氣流量(Qe)均有利于提高放電的穩(wěn)定性、增大射流尺度,并且也使得擊穿氣體產(chǎn)生放電變得更加容易.在本文所涉及的工況范圍內(nèi),等離子體射流的直徑和長(zhǎng)度分別在10.0—14.0 cm以及42.0—60.0 cm的范圍內(nèi)變化.

    圖9 不同放電工況下的等離子體射流的最大直徑Fig.9.Maximum values of the plasma jet diameter under different operating conditions.

    圖10 不同放電工況下的等離子體射流長(zhǎng)度Fig.10.Lengths of the plasma jet under different operating conditions.

    4 等離子體沖擊射流特性的實(shí)驗(yàn)研究

    4.1 典型實(shí)驗(yàn)工況下的等離子體沖擊射流特性

    圖11給出了典型工況(Qe=1.0 slpm,Qc=10.0 slpm,I=100 A,le=30 mm,r1=4.0 cm,p=1000 Pa)、有鈍體存在條件下氬-氮混合氣體(氮含量10%)放電所產(chǎn)生的等離子體沖擊射流的形貌(texp=20 ms).由于MPX-2015真空腔測(cè)試段長(zhǎng)度的限制,所采用的鈍體長(zhǎng)度為35.0 cm,最大直徑為5.0 cm,鈍體尖端與發(fā)生器水冷約束管出口之間的距離(L1)為8.0 cm.從圖11可以看到,從水冷約束管噴出的等離子體射流在沿流動(dòng)方向緊貼鈍體逐漸擴(kuò)張,整個(gè)鈍體全部被包覆在等離子體射流區(qū)內(nèi).

    圖12給出了相同工況下距離鈍體表面2.0 cm處等離子體射流區(qū)氣體溫度、電子激發(fā)溫度和電子數(shù)密度沿射流流向的變化規(guī)律(以鈍體尖端為坐標(biāo)零點(diǎn)).可以看到:1)上述三個(gè)特征參數(shù)(Th,Te和ne)均沿射流的流動(dòng)方向在z=0—40.0 cm的范圍內(nèi)逐漸降低,其中,電子激發(fā)溫度和重粒子溫度分別從尖端處的6120 K和3050 K降低到了3100 K和1020 K,而電子數(shù)密度則從4.00×1014cm?3下降到了1.40×1011cm?3;2)相應(yīng)地,等離子體射流區(qū)電子數(shù)密度大于ne,crt(=1010cm?3)的等離子體鞘套厚度亦從鈍體尖端處的4.7 cm擴(kuò)大到了7.6 cm(如圖13所示),有效覆蓋長(zhǎng)度可以達(dá)到約33.0 cm,占鈍體總長(zhǎng)度(35.0 cm)的94.2%.

    圖11 典型工況有鈍體存在條件下的等離子體沖擊射流照片F(xiàn)ig.11.Image of the plasma impinging jet upon a bluffbody under a typical operating condition.

    圖12 距鈍體表面2.0 cm處等離子體射流區(qū)氣體溫度、電子激發(fā)溫度和數(shù)密度沿軸線的變化Fig.12.Distributions of Th,Teand nealong the geometrical axis of the plasma jet and 2.0 cm away from the surface of the bluffbody.

    4.2 主要工作參數(shù)對(duì)等離子體沖擊射流特性的影響

    表2列出了電極間距、水冷約束管內(nèi)徑、從電極平面上游注入的冷氣體量以及鈍體半徑(r2)均保持不變的條件下(le=3.0 cm,r1=4.0 cm,Qc=10.0 slpm,r2=2.5 cm)研究真空腔壓強(qiáng)、弧電流、等離子體工作氣體流量對(duì)等離子體沖擊射流特性影響規(guī)律的8組實(shí)驗(yàn)工況.

    圖13 等離子體鞘套厚度沿發(fā)生器幾何軸線的變化規(guī)律(δcrit為等離子體射流有效覆蓋厚度下限)Fig.13.Variations of the plasma layer thickness along the geometrical axis of the plasma jet(δcritis the lower limit of the effective plasma layer thickness).

    表2 沖擊射流特性研究實(shí)驗(yàn)工況表Table 2.Operating conditions for the plasma impinging jets.

    圖14和圖15分別給出了對(duì)應(yīng)于表2中不同工況下所形成的等離子體鞘套有效工作長(zhǎng)度和最大厚度,其中,從圖11中O點(diǎn)開(kāi)始至鈍體根部,軸向方向每隔5.0 cm、徑向方向每隔1.0 cm進(jìn)行等離子體射流的光譜空間分布掃描測(cè)量,射流長(zhǎng)度和直徑均以電子密度約ne,crt為邊界得到的最大直徑以及長(zhǎng)度定義.實(shí)驗(yàn)測(cè)量結(jié)果表明:1)真空腔壓強(qiáng)(p)會(huì)對(duì)鈍體區(qū)的等離子體鞘套厚度產(chǎn)生顯著的影響,隨著壓強(qiáng)的升高,等離子體鞘套的有效工作長(zhǎng)度和厚度均減小;2)相同真空腔壓強(qiáng)條件下,增大弧電流(I)和沿電極環(huán)縫注入的氣流量(Qe)均會(huì)使得所形成的等離子體鞘套長(zhǎng)度和厚度有所增加,同時(shí)也使擊穿氣體產(chǎn)生放電變得更加容易,但相比真空腔壓強(qiáng)的影響要小得多.在本文所研究的參數(shù)范圍內(nèi),等離子體射流包覆長(zhǎng)35.0 cm、半徑2.5 cm的鈍體后所形成的等離子體鞘套的最大厚度和長(zhǎng)度分別在7.2—8.6 cm以及27.0—35.0 cm的范圍內(nèi)變化.

    圖14 不同放電工況下的等離子體鞘套厚度Fig.14.Plasma layer thicknesses surrounding the bluffbody under different operating conditions.

    圖15 不同放電工況下的等離子體鞘套長(zhǎng)度Fig.15.Plasma layer lengths surrounding the bluffbody under different operating conditions.

    5 結(jié)論與展望

    本文以臨近空間高超聲速飛行器以及航天器再入大氣環(huán)境飛行過(guò)程中的“黑障”問(wèn)題為研究背景,建立了六相交流電弧等離子體實(shí)驗(yàn)平臺(tái)(MPX-2015),對(duì)其所產(chǎn)生的等離子體自由射流和有鈍體存在條件下的沖擊射流特性進(jìn)行了初步實(shí)驗(yàn)研究.在實(shí)驗(yàn)參數(shù)范圍內(nèi)得到的主要結(jié)論為:

    1)采用本文所提出的六相交流電弧放電等離子體射流源結(jié)構(gòu)設(shè)計(jì),能夠在中等壓力范圍(500—1000 Pa)內(nèi)獲得比常規(guī)直流電弧放電和射頻感應(yīng)耦合放電等離子體源更大尺寸的等離子體射流;

    2)等離子體工作氣體流量、真空腔壓強(qiáng)、電極間距以及弧電流等因素均會(huì)對(duì)等離子體射流特性產(chǎn)生一定的影響,其中,真空腔壓強(qiáng)對(duì)等離子體射流特性的影響最為顯著,隨著壓強(qiáng)由500 Pa升高到1000 Pa,電子與中性粒子間的碰撞頻率提高,電子數(shù)密度沿軸向迅速衰減,等離子體自由射流長(zhǎng)度和最大直徑減小;在沖擊鈍體條件下所形成的等離子體鞘套的有效工作長(zhǎng)度和厚度均減小;而相同氣壓條件下,提高沿電極環(huán)縫注入的等離子體工作氣體流量(由0.5 slpm增加到1.0 slpm)或放電的弧電流(由80 A增加到100 A)均有利于等離子體鞘套尺寸的增大.

    本文著重研究了當(dāng)放電達(dá)到準(zhǔn)穩(wěn)態(tài)時(shí),不同工況下時(shí)間平均的等離子體射流長(zhǎng)度、直徑以及射流區(qū)氣體溫度、電子激發(fā)溫度和電子數(shù)密度等關(guān)鍵參數(shù),但并未對(duì)等離子體射流隨時(shí)間的演化特性以及時(shí)間平均的等離子體射流速度分布等進(jìn)行深入的研究.今后應(yīng)進(jìn)一步開(kāi)展如下研究工作:設(shè)計(jì)由信號(hào)發(fā)生器控制的放電圖像與光譜信號(hào)采集系統(tǒng)[34?36],研究發(fā)生器的放電特性以及放電過(guò)程中等離子體電弧射流特性隨時(shí)間的變化規(guī)律;設(shè)計(jì)具有不同型線的水冷約束管獲得超聲速的等離子體射流,以期更好地模擬臨近空間飛行器與其周圍復(fù)雜介質(zhì)環(huán)境間的氣動(dòng)熱效應(yīng)和“黑障”問(wèn)題.

    [1]Keidar M,Kim M,Boyd I D 2008J.Spacecr.Rockets45 445

    [2]Morris R A,Bench P M,Golden K E,Sutton E A 199937th Aerospace Sciences Meeting and ExhibitReno,NV,USA,January 11–14,1999 AIAA-99-0630

    [3]Evans J S,Schexnayder Jr C J,Huber P W 1973NASA TN D-7332

    [4]Gillman E D,Foster J E,Blankson I M 2010NASA/TM-2010-216220

    [5]Rybak J P,Churchill R J 1971IEEE Trans.Aerosp.Electron.Syst.7 879

    [6]Mather D E,Pasqual J M,Sillence J P,Lewis P 2005AIAA/CIRA 13th International Space Planes and Hypersonics Systems and Technologies ConferenceCapua,Italy,May 16–20,2005 AIAA-2005-3443

    [7]Akey N D 1971NASA Special Publication252 19

    [8]Watillon P,Berthe P,Chavagnac C 2003AIAA/ICAS International Air and Space Symposium and Exposition:The Next 100 YearsDayton,OH,USA,July 14–17,2003 AIAA-2003-2913

    [9]Shirouzu M,Yamamoto M 1996Space Plane and Hypersonic Systems and Technology ConferenceNorfolk,VA,USA,November 18–22,1996 AIAA-96-4524-CP

    [10]Yanagihara M,Munenaga T 200424th International CongressoftheAeronauticalSciencesYokohama,Japan,August 29–September 3,2004 p2004-7

    [11]Sakurai H,Kobayasi M,Yamazaki I,Shirouzu M,Yamamoto M 1997Acta Astronaut.40 105

    [12]Auweter-Kurtz M,Kurtz H L,Laure S 1996J.Propul.Power12 1053

    [13]Zhao L,Liu X X,Su H S 2015Journal of Telemetry Tracking and Command36 28(in Chinese)[趙良,劉秀祥,蘇漢生2015遙測(cè)遙控36 28]

    [14]Hermann T,L?hle S,Zander F,Fulge H,Fasoulas S 2016J.Thermophys.Heat Transf.30 673

    [15]Lemmer K M 2009Ph.D.Dissertation(Ann Arbor:University of Michigan)

    [16]Yang M,Li X,Xie K,Liu Y,Liu D 2013Phys.Plasmas20 012101

    [17]Shashurin A,Zhuang T,Teel G,Keidar M,Kundrapu M,Loverich J,Beilis I I,Raitses Y 2014J.Spacecr.Rockets51 838

    [18]Zhang B L,Zhuang Z,Li Y W,Wang Y T,Duan P Z,Zhang M K 2017High Voltage Engineering43 3055(in Chinese)[張百靈,莊重,李益文,王宇天,段朋振,張茗柯2017高電壓技術(shù)43 3055]

    [19]Pan W,Zhang W,Zhang W,Wu C 2001Plasma Chem.Plasma Process.21 23

    [20]Fauchais P,Vardelle A 1997IEEE Trans.Plasma Sci.25 1258

    [21]Colombo V,Concetti A,Ghedini E,Rotundo F,Sanibondi P,Boselli M,Dallavalle S,Gherardi M,Nemchinsky V,Vancini M 2012Plasma Chem.Plasma Process.32 411

    [22]Vardelle A,Moreau C,Akedo J,et al.2016J.Therm.Spray Technol.25 1376

    [23]Jin F,Li P,Ge N 2014High Voltage Engineering40 2057(in Chinese)[金鋒,李鵬,葛楠 2014高電壓技術(shù) 40 2057]

    [24]Riaby V A,Masherov P E,Obukhov V A,Savinov V P 2013High Voltage Engineering39 30596

    [25]Yao Y,Hossain M M,Watanabe T,Matsuura T,Funabiki F,Yano T 2008Chem.Eng.J.139 390

    [26]Watanabe T,Liu Y,Tanaka M 2014Plasma Chem.Plasma Process.34 443

    [27]Wang Z,Wu G Q,Ge N,Li H P,Bao C Y 2010IEEE Trans.Plasma Sci.38 2906

    [28]RaizerY P 1991GasDischargePhysics(Berlin:Springer)p136

    [29]Zhang H,Wu G Q,Li H P,Bao C Y 2009IEEE Trans.Plasma Sci.37 1129

    [30]Mahmood S,Shaikh N M,Kalyar M A,Ra fi q M,Piracha N K,Baig M A 2009J.Quant.Spectrosc.Radiat.Transf.110 1840

    [31]de Izarra C 2000J.Phys.D:Appl.Phys.33 1697

    [32]Park J,Henins I,Herrmann H W,Selwyn G S 2000Phys.Plasmas7 3141

    [33]Guo H,Zhang X N,Nie Q Y,Li H P,Zeng S,Li Z H 2018Acta Phys.Sin.67 055201(in Chinese)[郭恒, 張曉寧,聶秋月,李和平,曾實(shí),李志輝 2018物理學(xué)報(bào) 67 055201]

    [34]He L M,Lei J P,Chen Y,Liu X J,Chen G C,Zeng H 2017High Voltage Engineering43 3061(in Chinese)[何立明,雷健平,陳一,劉興建,陳高成,曾昊 2017高電壓技術(shù)43 3061]

    [35]Li X D,Zhang M,Zhu F S,Zhang H,Bo Z 2015High Voltage Engineering41 2022(in Chinese)[李曉東,張明,朱鳳森,張浩,薄拯2015高電壓技術(shù)41 2022]

    [36]Zhang H,He L,Yu J,Qi W,Chen G 2018Plasma Sci.Technol.20 024001

    猜你喜歡
    空腔電弧射流
    深海逃逸艙射流注水均壓過(guò)程仿真分析
    低壓天然氣泄漏射流擴(kuò)散特性研究
    煤氣與熱力(2022年4期)2022-05-23 12:45:00
    故障電弧探測(cè)器與故障電弧保護(hù)裝置在工程中的應(yīng)用分析
    基于邊光滑有限元法的二維復(fù)合彈性空腔聲振特性分析
    2219鋁合金激光電弧復(fù)合焊接及其溫度場(chǎng)的模擬
    空腔參數(shù)對(duì)重力壩穩(wěn)定的影響分析
    航空電氣系統(tǒng)中故障電弧的分析
    電子制作(2017年22期)2017-02-02 07:10:15
    前置污水去油池
    前置污水去油池
    射流齒形噴嘴射流流場(chǎng)與氣動(dòng)聲學(xué)分析
    免费在线观看黄色视频的| 国产精品国产三级国产专区5o| 午夜精品久久久久久毛片777| 捣出白浆h1v1| 欧美亚洲日本最大视频资源| 啦啦啦中文免费视频观看日本| 少妇人妻久久综合中文| 国内毛片毛片毛片毛片毛片| 在线观看免费视频网站a站| 久久国产精品影院| 人人妻人人添人人爽欧美一区卜| 如日韩欧美国产精品一区二区三区| 久久精品人人爽人人爽视色| 天天影视国产精品| a 毛片基地| 丁香六月欧美| 日本一区二区免费在线视频| 美女脱内裤让男人舔精品视频| 在线观看免费视频网站a站| 少妇粗大呻吟视频| 精品国产国语对白av| 成年女人毛片免费观看观看9 | 后天国语完整版免费观看| 一区二区三区激情视频| 男女国产视频网站| 亚洲av日韩在线播放| 国产欧美日韩一区二区三 | 青春草视频在线免费观看| 亚洲精品在线美女| 久久 成人 亚洲| 国产97色在线日韩免费| 麻豆av在线久日| 性高湖久久久久久久久免费观看| 久久精品国产亚洲av香蕉五月 | 久久久久久久久久久久大奶| 美女高潮喷水抽搐中文字幕| avwww免费| 欧美日韩黄片免| 又黄又粗又硬又大视频| 欧美大码av| 成人黄色视频免费在线看| 亚洲国产日韩一区二区| 免费黄频网站在线观看国产| 精品国产超薄肉色丝袜足j| 日日摸夜夜添夜夜添小说| 久久久久视频综合| 母亲3免费完整高清在线观看| 久久亚洲精品不卡| 少妇粗大呻吟视频| 欧美国产精品va在线观看不卡| 正在播放国产对白刺激| 欧美日韩一级在线毛片| 中文字幕人妻丝袜制服| 美女中出高潮动态图| 少妇粗大呻吟视频| 国产日韩欧美视频二区| 国产成人a∨麻豆精品| 欧美日韩国产mv在线观看视频| 午夜91福利影院| 又黄又粗又硬又大视频| 美国免费a级毛片| 久久精品亚洲熟妇少妇任你| 久久99热这里只频精品6学生| 夜夜夜夜夜久久久久| 欧美精品一区二区大全| 日本a在线网址| 日本精品一区二区三区蜜桃| 热99久久久久精品小说推荐| 丝袜喷水一区| 国产在线一区二区三区精| 久久久久久久国产电影| 乱人伦中国视频| 精品亚洲成a人片在线观看| 99国产极品粉嫩在线观看| 男人添女人高潮全过程视频| 久久av网站| 考比视频在线观看| 精品国产一区二区三区四区第35| 亚洲精品美女久久av网站| 性高湖久久久久久久久免费观看| 女人久久www免费人成看片| 18禁国产床啪视频网站| 一个人免费看片子| 一个人免费看片子| 男女无遮挡免费网站观看| 久热这里只有精品99| 一区二区日韩欧美中文字幕| 亚洲,欧美精品.| 女人久久www免费人成看片| 高清视频免费观看一区二区| 蜜桃在线观看..| 久久毛片免费看一区二区三区| 男女床上黄色一级片免费看| 少妇裸体淫交视频免费看高清 | 久久免费观看电影| 91精品三级在线观看| 亚洲精品一二三| 欧美97在线视频| 国产一区有黄有色的免费视频| 在线观看人妻少妇| 美女午夜性视频免费| 精品视频人人做人人爽| 欧美在线黄色| 国产在线免费精品| 亚洲精品粉嫩美女一区| 黄色视频,在线免费观看| 交换朋友夫妻互换小说| 天天添夜夜摸| 视频在线观看一区二区三区| 国产黄色免费在线视频| 啦啦啦免费观看视频1| 久久精品国产a三级三级三级| 精品亚洲成国产av| 久久精品国产a三级三级三级| 久久久国产欧美日韩av| 欧美av亚洲av综合av国产av| 波多野结衣av一区二区av| 国产野战对白在线观看| 一个人免费看片子| 亚洲男人天堂网一区| 人人妻人人爽人人添夜夜欢视频| 啦啦啦 在线观看视频| 丝袜美腿诱惑在线| 欧美日韩成人在线一区二区| 欧美人与性动交α欧美精品济南到| 欧美激情高清一区二区三区| 18在线观看网站| 最近中文字幕2019免费版| 欧美午夜高清在线| 午夜福利免费观看在线| 国产免费福利视频在线观看| 美女大奶头黄色视频| 亚洲熟女精品中文字幕| 亚洲成国产人片在线观看| 高清在线国产一区| 精品少妇内射三级| 美国免费a级毛片| netflix在线观看网站| 又黄又粗又硬又大视频| 啦啦啦啦在线视频资源| 免费黄频网站在线观看国产| 欧美激情 高清一区二区三区| 大香蕉久久网| 男女边摸边吃奶| 麻豆国产av国片精品| 菩萨蛮人人尽说江南好唐韦庄| 国产成人欧美在线观看 | 97在线人人人人妻| 国产91精品成人一区二区三区 | 国产精品国产av在线观看| 亚洲精品久久成人aⅴ小说| 美女中出高潮动态图| 亚洲国产精品999| 美女高潮到喷水免费观看| 久久人妻福利社区极品人妻图片| 狠狠狠狠99中文字幕| a在线观看视频网站| 99热网站在线观看| 黑人欧美特级aaaaaa片| 大片免费播放器 马上看| 老司机影院成人| 超碰97精品在线观看| 91成人精品电影| 美女中出高潮动态图| 欧美黑人欧美精品刺激| 日韩欧美免费精品| 麻豆av在线久日| 这个男人来自地球电影免费观看| 国产精品偷伦视频观看了| 欧美人与性动交α欧美精品济南到| 国产精品秋霞免费鲁丝片| 成人国产av品久久久| 亚洲av电影在线进入| 国产一卡二卡三卡精品| 美女高潮喷水抽搐中文字幕| 成年人免费黄色播放视频| 免费高清在线观看日韩| 狠狠狠狠99中文字幕| 嫁个100分男人电影在线观看| 别揉我奶头~嗯~啊~动态视频 | 飞空精品影院首页| 新久久久久国产一级毛片| 亚洲第一av免费看| 久久毛片免费看一区二区三区| 老汉色av国产亚洲站长工具| 色94色欧美一区二区| 日韩欧美免费精品| 免费看十八禁软件| 视频区图区小说| 91麻豆av在线| 中文字幕最新亚洲高清| 免费一级毛片在线播放高清视频 | 天天添夜夜摸| 欧美精品人与动牲交sv欧美| 超色免费av| 99久久国产精品久久久| 爱豆传媒免费全集在线观看| 欧美另类亚洲清纯唯美| 久久午夜综合久久蜜桃| 50天的宝宝边吃奶边哭怎么回事| 啦啦啦在线免费观看视频4| 无限看片的www在线观看| 欧美激情久久久久久爽电影 | 久久久久久久久久久久大奶| 久久性视频一级片| 后天国语完整版免费观看| 五月天丁香电影| 五月天丁香电影| 亚洲国产欧美日韩在线播放| 免费黄频网站在线观看国产| 少妇被粗大的猛进出69影院| 国产精品.久久久| 美女脱内裤让男人舔精品视频| 久久久久久久久久久久大奶| 欧美日韩一级在线毛片| 欧美黄色淫秽网站| 高清在线国产一区| 真人做人爱边吃奶动态| a级毛片黄视频| 亚洲中文字幕日韩| 成人国产一区最新在线观看| av在线老鸭窝| av视频免费观看在线观看| 91精品三级在线观看| 99国产综合亚洲精品| 肉色欧美久久久久久久蜜桃| 久久精品成人免费网站| 亚洲欧美成人综合另类久久久| 国产成人啪精品午夜网站| 久久 成人 亚洲| 久久精品亚洲熟妇少妇任你| 999精品在线视频| 国产亚洲av片在线观看秒播厂| 欧美黑人精品巨大| 一二三四社区在线视频社区8| 欧美xxⅹ黑人| 国产黄频视频在线观看| 精品久久久久久电影网| 午夜久久久在线观看| 欧美黑人精品巨大| 97在线人人人人妻| h视频一区二区三区| kizo精华| 欧美激情极品国产一区二区三区| 大型av网站在线播放| 国产色视频综合| 老司机午夜福利在线观看视频 | 国产精品 国内视频| 国产成人一区二区三区免费视频网站| 日本猛色少妇xxxxx猛交久久| 啦啦啦在线免费观看视频4| 性色av一级| 乱人伦中国视频| 欧美老熟妇乱子伦牲交| 日韩三级视频一区二区三区| 男女免费视频国产| 精品国产乱码久久久久久小说| 亚洲欧洲日产国产| 日日摸夜夜添夜夜添小说| a 毛片基地| 99国产极品粉嫩在线观看| 精品乱码久久久久久99久播| 少妇裸体淫交视频免费看高清 | 天天躁日日躁夜夜躁夜夜| 一级黄色大片毛片| 乱人伦中国视频| 在线永久观看黄色视频| 日韩视频一区二区在线观看| 国产黄频视频在线观看| 国产男女内射视频| 天堂8中文在线网| 超碰成人久久| 成人国产一区最新在线观看| 久久亚洲精品不卡| 国产精品自产拍在线观看55亚洲 | 侵犯人妻中文字幕一二三四区| 国产日韩欧美视频二区| 国产精品麻豆人妻色哟哟久久| 十八禁网站免费在线| 中文字幕人妻丝袜制服| 午夜精品久久久久久毛片777| 菩萨蛮人人尽说江南好唐韦庄| 丁香六月欧美| 亚洲专区中文字幕在线| 亚洲男人天堂网一区| 国产精品 国内视频| 桃花免费在线播放| 99国产精品99久久久久| 丰满饥渴人妻一区二区三| 日韩精品免费视频一区二区三区| 黄频高清免费视频| 国产极品粉嫩免费观看在线| 黑人欧美特级aaaaaa片| 亚洲va日本ⅴa欧美va伊人久久 | 99国产精品免费福利视频| 久久免费观看电影| 国产亚洲欧美精品永久| 精品国产国语对白av| 亚洲全国av大片| 亚洲第一青青草原| 美女大奶头黄色视频| 在线av久久热| 91九色精品人成在线观看| 精品亚洲成a人片在线观看| 国产精品九九99| 久久精品亚洲av国产电影网| 桃红色精品国产亚洲av| 亚洲av国产av综合av卡| 欧美日韩亚洲国产一区二区在线观看 | 黄片大片在线免费观看| 国产高清videossex| 另类精品久久| 国产有黄有色有爽视频| 久久久国产一区二区| 少妇 在线观看| 欧美黑人欧美精品刺激| 热99re8久久精品国产| 一本久久精品| 久热这里只有精品99| 色精品久久人妻99蜜桃| av国产精品久久久久影院| 男人舔女人的私密视频| 久久久久精品人妻al黑| e午夜精品久久久久久久| av不卡在线播放| 日韩视频在线欧美| 国产精品成人在线| 一区二区三区四区激情视频| 午夜免费成人在线视频| 肉色欧美久久久久久久蜜桃| 国产激情久久老熟女| 国产精品99久久99久久久不卡| 日韩欧美一区二区三区在线观看 | 法律面前人人平等表现在哪些方面 | 欧美亚洲日本最大视频资源| 老汉色av国产亚洲站长工具| 欧美日韩福利视频一区二区| xxxhd国产人妻xxx| 午夜精品国产一区二区电影| 日韩熟女老妇一区二区性免费视频| 高清av免费在线| 欧美日韩亚洲高清精品| 色94色欧美一区二区| 欧美黄色淫秽网站| 不卡av一区二区三区| 精品一区二区三卡| 精品国产乱码久久久久久男人| 色综合欧美亚洲国产小说| 亚洲一区中文字幕在线| 精品人妻在线不人妻| 在线观看免费日韩欧美大片| 久久久久久久久免费视频了| 婷婷成人精品国产| videosex国产| 亚洲欧美成人综合另类久久久| 国产免费一区二区三区四区乱码| 丝袜喷水一区| av超薄肉色丝袜交足视频| 久久久国产一区二区| 黑人巨大精品欧美一区二区mp4| 蜜桃在线观看..| 欧美日韩福利视频一区二区| 一区二区三区乱码不卡18| 国产亚洲精品第一综合不卡| 亚洲av日韩在线播放| 国产片内射在线| 亚洲成人免费电影在线观看| 久久中文看片网| 麻豆国产av国片精品| 美女福利国产在线| av天堂在线播放| 国产伦人伦偷精品视频| 一级毛片电影观看| 女人被躁到高潮嗷嗷叫费观| 久久中文字幕一级| 亚洲avbb在线观看| 国产无遮挡羞羞视频在线观看| 中文字幕人妻丝袜一区二区| 波多野结衣一区麻豆| 亚洲 欧美一区二区三区| 久久久久久人人人人人| 国产一区二区三区在线臀色熟女 | 亚洲伊人久久精品综合| 我要看黄色一级片免费的| 午夜影院在线不卡| 一级,二级,三级黄色视频| 国产成人啪精品午夜网站| 亚洲三区欧美一区| 日韩,欧美,国产一区二区三区| 免费在线观看黄色视频的| 国产免费一区二区三区四区乱码| 精品亚洲成a人片在线观看| 99国产综合亚洲精品| 亚洲第一欧美日韩一区二区三区 | av网站在线播放免费| 亚洲专区中文字幕在线| 汤姆久久久久久久影院中文字幕| 亚洲精品国产一区二区精华液| 亚洲五月色婷婷综合| 国产精品二区激情视频| av福利片在线| 高清欧美精品videossex| 美女中出高潮动态图| 欧美成狂野欧美在线观看| 亚洲成人手机| 天天躁狠狠躁夜夜躁狠狠躁| 大型av网站在线播放| 久久国产精品人妻蜜桃| 日韩制服丝袜自拍偷拍| 国产有黄有色有爽视频| 久久久精品94久久精品| 性色av一级| 黄色视频,在线免费观看| 国产成人啪精品午夜网站| 午夜两性在线视频| 91精品三级在线观看| av片东京热男人的天堂| 九色亚洲精品在线播放| 欧美变态另类bdsm刘玥| 亚洲av男天堂| 欧美激情久久久久久爽电影 | 午夜免费鲁丝| 男女午夜视频在线观看| 一级片免费观看大全| 黑人欧美特级aaaaaa片| 天天影视国产精品| 日韩大码丰满熟妇| 亚洲色图综合在线观看| 欧美精品高潮呻吟av久久| 两性夫妻黄色片| 美国免费a级毛片| 热99久久久久精品小说推荐| 久久香蕉激情| 亚洲免费av在线视频| 亚洲专区国产一区二区| 久久中文字幕一级| 亚洲情色 制服丝袜| 精品人妻1区二区| 亚洲黑人精品在线| 日日摸夜夜添夜夜添小说| 波多野结衣av一区二区av| 十八禁人妻一区二区| 亚洲av成人一区二区三| 一本一本久久a久久精品综合妖精| 国产一区二区在线观看av| 老司机影院成人| 男女无遮挡免费网站观看| 国产国语露脸激情在线看| 亚洲av日韩精品久久久久久密| 视频在线观看一区二区三区| 精品国产国语对白av| 欧美黑人欧美精品刺激| 亚洲综合色网址| 欧美 日韩 精品 国产| 交换朋友夫妻互换小说| 黄色毛片三级朝国网站| 午夜免费观看性视频| 亚洲伊人久久精品综合| 黄网站色视频无遮挡免费观看| 久久人妻福利社区极品人妻图片| 淫妇啪啪啪对白视频 | 女人久久www免费人成看片| 亚洲avbb在线观看| 日本91视频免费播放| 日本av免费视频播放| 亚洲欧美精品综合一区二区三区| 90打野战视频偷拍视频| 搡老乐熟女国产| 欧美97在线视频| 视频区图区小说| 久久久久久免费高清国产稀缺| 天天操日日干夜夜撸| 91精品国产国语对白视频| 国产精品国产av在线观看| 少妇的丰满在线观看| 午夜免费成人在线视频| 欧美日韩成人在线一区二区| 夜夜骑夜夜射夜夜干| 五月开心婷婷网| 黑人欧美特级aaaaaa片| 国产精品香港三级国产av潘金莲| 国产精品秋霞免费鲁丝片| 国产精品一二三区在线看| 黄频高清免费视频| 91精品伊人久久大香线蕉| 黄片小视频在线播放| 欧美另类亚洲清纯唯美| 精品亚洲乱码少妇综合久久| 香蕉国产在线看| 男人爽女人下面视频在线观看| 国产免费现黄频在线看| 中国国产av一级| 久久精品熟女亚洲av麻豆精品| 99热网站在线观看| 亚洲国产成人一精品久久久| 这个男人来自地球电影免费观看| 亚洲va日本ⅴa欧美va伊人久久 | 波多野结衣一区麻豆| 亚洲专区国产一区二区| 久久精品国产亚洲av高清一级| 亚洲国产欧美网| 久久精品人人爽人人爽视色| 亚洲精品国产一区二区精华液| 亚洲国产看品久久| 亚洲精品乱久久久久久| 国产xxxxx性猛交| 99香蕉大伊视频| 精品福利观看| 国产av精品麻豆| 窝窝影院91人妻| 国产精品.久久久| 人人澡人人妻人| 国产精品国产av在线观看| 免费久久久久久久精品成人欧美视频| 精品亚洲成a人片在线观看| 国产极品粉嫩免费观看在线| 久久久精品区二区三区| 男女边摸边吃奶| 男人爽女人下面视频在线观看| 亚洲七黄色美女视频| 欧美日韩国产mv在线观看视频| 日韩三级视频一区二区三区| 黑人欧美特级aaaaaa片| 国产精品久久久人人做人人爽| 99热网站在线观看| av天堂在线播放| 免费在线观看黄色视频的| 久久久久久免费高清国产稀缺| 电影成人av| 国产精品熟女久久久久浪| videosex国产| av国产精品久久久久影院| 深夜精品福利| 亚洲国产精品一区三区| 久久久精品免费免费高清| 自拍欧美九色日韩亚洲蝌蚪91| av天堂在线播放| 国产极品粉嫩免费观看在线| 热re99久久精品国产66热6| 免费看十八禁软件| 欧美激情久久久久久爽电影 | 搡老熟女国产l中国老女人| 亚洲国产日韩一区二区| 男男h啪啪无遮挡| 人人妻人人澡人人爽人人夜夜| 丁香六月欧美| 国产成人精品在线电影| 咕卡用的链子| 大陆偷拍与自拍| av不卡在线播放| 狠狠狠狠99中文字幕| 亚洲自偷自拍图片 自拍| 久久中文字幕一级| 国产在视频线精品| 精品国产一区二区三区久久久樱花| 黑人操中国人逼视频| 成年美女黄网站色视频大全免费| 日日爽夜夜爽网站| 亚洲国产日韩一区二区| 国产1区2区3区精品| 在线天堂中文资源库| 成人18禁高潮啪啪吃奶动态图| 国产精品国产av在线观看| 国产亚洲午夜精品一区二区久久| 久久狼人影院| √禁漫天堂资源中文www| 久久精品aⅴ一区二区三区四区| 国产成人啪精品午夜网站| 久久精品aⅴ一区二区三区四区| 欧美中文综合在线视频| 午夜免费观看性视频| 欧美少妇被猛烈插入视频| 久久毛片免费看一区二区三区| 正在播放国产对白刺激| 亚洲成人国产一区在线观看| 久久国产精品人妻蜜桃| 免费在线观看完整版高清| 精品第一国产精品| 美女高潮喷水抽搐中文字幕| 黄片播放在线免费| 亚洲久久久国产精品| 久久99热这里只频精品6学生| 伊人亚洲综合成人网| av不卡在线播放| 成年美女黄网站色视频大全免费| 亚洲国产精品成人久久小说| 交换朋友夫妻互换小说| 九色亚洲精品在线播放| 亚洲精品自拍成人| 中文字幕高清在线视频| 热99久久久久精品小说推荐| 欧美大码av| 考比视频在线观看| 国产成人免费观看mmmm| 欧美+亚洲+日韩+国产| 菩萨蛮人人尽说江南好唐韦庄| 新久久久久国产一级毛片| 成人国语在线视频| 亚洲精品久久成人aⅴ小说| 丰满迷人的少妇在线观看| www.av在线官网国产| 中文字幕av电影在线播放| 最近中文字幕2019免费版| www日本在线高清视频| 久久久久久久久久久久大奶| 啦啦啦啦在线视频资源| 亚洲人成电影免费在线| 97在线人人人人妻| 色综合欧美亚洲国产小说| 久久久精品免费免费高清| 黑人欧美特级aaaaaa片| 国产欧美日韩精品亚洲av| 啦啦啦视频在线资源免费观看| videosex国产| 夜夜骑夜夜射夜夜干|