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

    磁流體動(dòng)力加速風(fēng)洞技術(shù)發(fā)展分析

    2019-01-07 06:24:32左光齊玢歐東斌
    航天返回與遙感 2018年6期
    關(guān)鍵詞:研究

    左光 齊玢 歐東斌

    ?

    磁流體動(dòng)力加速風(fēng)洞技術(shù)發(fā)展分析

    左光1齊玢1歐東斌2

    (1 中國空間技術(shù)研究院 載人航天總體部,北京 100094)(2 中國航天空氣動(dòng)力技術(shù)研究院 電弧等離子應(yīng)用裝備北京市重點(diǎn)實(shí)驗(yàn)室,北京 100074)

    針對(duì)高超聲速飛行器地面模擬試驗(yàn)需求,傳統(tǒng)試驗(yàn)方法難以實(shí)現(xiàn)真實(shí)氣體溫度、清潔空氣、大尺度、長時(shí)間、高馬赫數(shù)模擬能力,磁流體動(dòng)力加速風(fēng)洞提供了全新技術(shù)路線。文章歸納了國內(nèi)外磁流體動(dòng)力加速風(fēng)洞研究發(fā)展現(xiàn)狀,介紹了磁流體動(dòng)力加速風(fēng)洞原理。文章對(duì)基于熱電離的磁流體動(dòng)力加速風(fēng)洞方案進(jìn)行了論述,采用高頻等離子發(fā)生器為設(shè)備提供加熱源,從而避免電極燒損所引起的污染問題,控制氣體總溫不超過3 500K,使氣體離子化低于30%,通過兩級(jí)加速,達(dá)到出口馬赫數(shù)為15的模擬環(huán)境。文章進(jìn)一步分析了磁流體動(dòng)力加速風(fēng)洞關(guān)鍵技術(shù)問題。超聲速氣流電離技術(shù)方面,核心問題在于超聲速氣流電離規(guī)律與機(jī)理以及電離種子注入、電子束電離等關(guān)鍵技術(shù);磁流體動(dòng)力加速通道設(shè)計(jì)方面,重點(diǎn)考慮氣流密度及磁感應(yīng)強(qiáng)度等因素的綜合影響以及電極設(shè)計(jì)技術(shù);高超聲速模擬測(cè)試方面,關(guān)鍵技術(shù)包括電磁屏蔽技術(shù)、微波干涉儀技術(shù)、平面激光誘導(dǎo)熒光技術(shù)、高分辨率高性能的光譜測(cè)試技術(shù)等。最后,提出了磁流體動(dòng)力加速風(fēng)洞技術(shù)發(fā)展建議。

    磁流體動(dòng)力 高超聲速 風(fēng)洞 地面模擬試驗(yàn) 天地往返飛行器

    0 引言

    天地往返飛行器是航天器重要的一個(gè)分支,傳統(tǒng)的載人飛船和返回式衛(wèi)星均屬于這一類。新的天地往返系統(tǒng)的發(fā)展趨勢(shì)是水平起飛水平滑翔著陸,這也是未來天地往返領(lǐng)域激烈爭(zhēng)奪的戰(zhàn)場(chǎng)??仗煲惑w的思路已經(jīng)被業(yè)界所重視,在這個(gè)領(lǐng)域航天器和航空器必將融合,界限也越來越模糊,基于組合動(dòng)力的高超聲速飛行器是近年來世界空天技術(shù)研究的前沿?zé)狳c(diǎn)。2004年,美國成功開展了高超聲速飛行器X-43A馬赫數(shù)為7和馬赫數(shù)為10的飛行演示驗(yàn)證,實(shí)現(xiàn)了8s和10s超燃沖壓發(fā)動(dòng)機(jī)推動(dòng)下的高超聲速飛行器自由飛行;2013年,X-51實(shí)現(xiàn)了馬赫數(shù)為4.8~5.1、飛行時(shí)間240s的有動(dòng)力飛行,標(biāo)志著高超聲速飛行器研究取得重大進(jìn)展[1-3]。近年來,國內(nèi)也在該領(lǐng)域取得了長足的進(jìn)步,隨著研究及工程化的深入,采用火箭基組合動(dòng)力(Rocket Based Combined Cycle, RBCC)、渦輪基組合動(dòng)力(Turbine Based Combined Cycle, TBCC)、空氣渦輪沖壓動(dòng)力(Air Turbo Rocket, ATR)以及渦輪火箭沖壓組合動(dòng)力(SABRE/Trijet)等發(fā)動(dòng)機(jī)的高超聲速飛行器方案,有望克服單一動(dòng)力飛行包線局限,以其寬空域和廣速度域具有較好的整體性能,特別是具備馬赫數(shù)為8以上長時(shí)間飛行能力,成為高超聲速飛行器的重要發(fā)展方向,在重復(fù)使用天地往返運(yùn)輸?shù)榷囝I(lǐng)域具有廣闊的應(yīng)用前景[4]。

    進(jìn)行天地往返飛行器氣動(dòng)設(shè)計(jì)時(shí),要依托仿真計(jì)算、地面試驗(yàn)來精確確定飛行器所受的氣動(dòng)載荷,在此基礎(chǔ)上提供數(shù)據(jù)輸入以完成飛行器的彈道設(shè)計(jì)、飛行控制設(shè)計(jì)、防熱結(jié)構(gòu)設(shè)計(jì)等分系統(tǒng)設(shè)計(jì)。天地往返飛行器氣動(dòng)分析若重點(diǎn)采用數(shù)值仿真技術(shù),那么在很多空域及速域缺乏驗(yàn)證與評(píng)估手段,具體表現(xiàn)在以下方面:1)數(shù)值仿真建模方面,網(wǎng)格拓?fù)浼敖诿婢W(wǎng)格對(duì)氣動(dòng)仿真結(jié)果存在較大影響;2)湍流模型方面,對(duì)現(xiàn)有湍流模型在精度上缺乏共同的認(rèn)識(shí),同時(shí),模型經(jīng)驗(yàn)參數(shù)修正也是面臨的難題;3)空間離散格式方面,不同空間離散格式數(shù)值特性差異明顯,對(duì)氣動(dòng)仿真結(jié)果影響明顯;4)限制器方面,不同限制器求解性能具有較大差異,且對(duì)氣動(dòng)仿真結(jié)果具有顯著影響。風(fēng)洞試驗(yàn)?zāi)軌蛴行?duì)氣動(dòng)仿真計(jì)算方法及模型進(jìn)行修正,對(duì)形成適用于全空域全速域的氣動(dòng)仿真方案,提高氣動(dòng)數(shù)值仿真精度,進(jìn)而減小彈道、飛行控制、防熱結(jié)構(gòu)等設(shè)計(jì)誤差具有重要意義。

    上述分析表明,天地往返飛行器的發(fā)展離不開風(fēng)洞試驗(yàn)設(shè)備的支撐,但在空氣動(dòng)力學(xué)地面試驗(yàn)設(shè)備上百年的發(fā)展歷程中,建造用于發(fā)展先進(jìn)天地往返飛行器所必須的、能夠模擬馬赫數(shù)為8以上飛行狀態(tài)的試驗(yàn)設(shè)備還存在技術(shù)瓶頸。為了在地面試驗(yàn)設(shè)備中開展吸氣式推進(jìn)系統(tǒng)的飛行器一體化設(shè)計(jì)性能測(cè)試和耐久力考核試驗(yàn),模擬全尺度條件下高超聲速飛行環(huán)境——包括速度、溫度、壓力、空氣組分和試驗(yàn)時(shí)間——是極其重要的。據(jù)資料顯示,即便美國等航天大國也缺乏真實(shí)氣體溫度、清潔空氣、大尺度、長時(shí)間、高馬赫數(shù)試驗(yàn)?zāi)芰ΑD壳耙呀?jīng)存在的儲(chǔ)能加熱潔凈空氣設(shè)備可模擬的氣流最高馬赫數(shù)約為7;燃燒加熱的污染空氣設(shè)備可模擬的氣流速最高馬赫數(shù)約為8;電弧加熱設(shè)備模擬的氣流最高馬赫數(shù)可以達(dá)到9~10,但試驗(yàn)尺度較小,且電極熔化會(huì)對(duì)空氣造成污染;壓縮加熱(如脈沖激波風(fēng)洞)可提供最高馬赫數(shù)約14~15的試驗(yàn)條件,但試驗(yàn)時(shí)間極短,為毫秒量級(jí)。1994年出版的NASA和美國國防部《國家設(shè)備研究報(bào)告》和1997年的美國國防部報(bào)告《航空試驗(yàn)設(shè)備評(píng)估》中確認(rèn)的研究?jī)?nèi)容揭示了吸氣式高超試驗(yàn)?zāi)芰εc相對(duì)應(yīng)的高超聲速飛行器飛行試驗(yàn)需求之間存在著巨大的差距[5]。

    對(duì)于發(fā)展馬赫數(shù)大于8的真實(shí)溫度的高超聲速風(fēng)洞,至少存在以下數(shù)個(gè)主要技術(shù)瓶頸或障礙:1)來流空氣所需的足夠能量;2)正確的空氣化學(xué)成分;3)模擬吸氣發(fā)動(dòng)機(jī)工作的長的工作時(shí)間;4)可承受長時(shí)間高溫高壓試驗(yàn)氣體的設(shè)備材料與工藝技術(shù)。近些年來,為了解決以上難題,出現(xiàn)了兩種技術(shù)途徑嘗試進(jìn)行技術(shù)突破,1)考慮通過集束能量向試驗(yàn)氣體中注入熱能;2)對(duì)試驗(yàn)氣體進(jìn)行磁流體動(dòng)力(Magneto Hydro Dynamics,MHD)加速,從而使氣流速度進(jìn)一步增加,總壓進(jìn)一步提高[6]。磁流體動(dòng)力加速試驗(yàn)設(shè)備因其獲得更高馬赫數(shù)試驗(yàn)條件的同時(shí)避免在噴管喉道前極高的飛行總溫,減少喉道傳熱設(shè)計(jì)難度,降低氣體離解效應(yīng)而受到國外研究人員的關(guān)注。各國針對(duì)該項(xiàng)技術(shù)開展了大量研究,從而研制能夠長時(shí)間穩(wěn)定工作的高焓高超聲速風(fēng)洞,同時(shí)這也成為超聲速氣流磁流體加速技術(shù)在工程領(lǐng)域的最早應(yīng)用。

    1 磁流體動(dòng)力加速風(fēng)洞發(fā)展現(xiàn)狀

    1.1 國外發(fā)展現(xiàn)狀

    國外磁流體動(dòng)力加速風(fēng)洞領(lǐng)域的研究開展較早且相關(guān)研究機(jī)構(gòu)較多,具有代表性的研究單位包括美國航空航天局蘭利研究中心(NASA LRC)、美國空軍阿諾德工程發(fā)展中心(USAF AEDC)、美國航空航天局馬歇爾空間飛行中心(NASA MSFC)以及俄羅斯中央空氣動(dòng)力研究所(TsAGI)等。

    (1)NASA LRC磁流體動(dòng)力加速器

    NASA LRC是磁流體加速的先驅(qū)并在20世紀(jì)60年代維持了超過十年的相關(guān)研究工作。他們的雙重目標(biāo)是開發(fā)用于行星再入的高超聲速風(fēng)洞并將研究結(jié)果應(yīng)用于空間推進(jìn)。采用的工作介質(zhì)是氮?dú)猓∟2),并添加銫(Cs)作為電離種子。NASA LRC設(shè)計(jì)了三種磁流體動(dòng)力加速器,第一種加速器出口為1cm×1cm、第二種加速器出口為2.54cm×2.54cm、第三種加速器出口為6.35cm×6.35cm,根據(jù)設(shè)計(jì)方案,可以模擬飛行器在53km高度下達(dá)到13km/s的出口速度,但由于設(shè)備限制,最終利用30~36對(duì)電極,達(dá)到9.6km/s的出口速度。

    (2)USAF AEDC磁流體動(dòng)力加速風(fēng)洞原理驗(yàn)證

    USAF AEDC制定了兩個(gè)磁流體動(dòng)力加速風(fēng)洞原理驗(yàn)證計(jì)劃,即“低密度驗(yàn)證計(jì)劃”(Low Density Focus Program,LORHO)和“高密度驗(yàn)證計(jì)劃”(High Density Focus Program,HIRHO)[7-8]。LORHO基于1.2MW電弧加熱器,磁場(chǎng)強(qiáng)度2T,利用117對(duì)電極,出口速度達(dá)到3.9km/s。該計(jì)劃同時(shí)設(shè)計(jì)制造了20MW火箭發(fā)動(dòng)機(jī)驅(qū)動(dòng)引導(dǎo)設(shè)備,但未進(jìn)行試驗(yàn)。HIRHO基于激波風(fēng)洞,磁場(chǎng)強(qiáng)度7.5T,利用46cm長的分段法拉第通道及11對(duì)電極,添加鉀作為電離種子,速度較磁流體加速前增加75%,但由于對(duì)于高超聲速飛行模擬的需求下降而終止。

    (3)USAF AEDC輻射驅(qū)動(dòng)/磁流體動(dòng)力混合加速方案

    磁流體動(dòng)力加速要求試驗(yàn)氣體具備導(dǎo)電性,高溫低壓氣體中加入堿金屬種子可以提高電導(dǎo)率,但僅在氣流靜溫大于2 500K時(shí)有效。然而,氣流靜溫大于2 500K時(shí),將產(chǎn)生較高的熵值,同時(shí)因高溫化學(xué)反應(yīng)而產(chǎn)生氣體雜質(zhì)(如NO等),從而導(dǎo)致模擬條件與真實(shí)飛行條件的差異。另外,加入的堿金屬也會(huì)對(duì)飛行器試驗(yàn)帶來未知干擾效應(yīng)。因此,USAF AEDC提出了一種馬赫數(shù)為8~15模擬能力的混合加速方案[9]。該方案首先利用電子束輻射增加的超聲速氣流內(nèi)能,通過噴管的膨脹使清潔空氣馬赫數(shù)達(dá)到12;進(jìn)一步,利用第二個(gè)連續(xù)操作的電子束或電弧以達(dá)到磁流體動(dòng)力加速所需的電導(dǎo)率,從而通過磁流體動(dòng)力加速使空氣馬赫數(shù)達(dá)到15。該混合方法可以避免電離污染及高溫種子的注入。其概念如圖1所示。

    圖1 輻射驅(qū)動(dòng)/磁流體動(dòng)力混合加速風(fēng)洞概念

    (4)NASA MSFC磁流體動(dòng)力加速風(fēng)洞

    NASA MSFC建造了磁流體動(dòng)力加速風(fēng)洞,并對(duì)磁流體加速器在未來大推力推進(jìn)器中的應(yīng)用進(jìn)行了可行性研究[10],如圖2所示。該風(fēng)洞采用1.5MW疊片電弧加熱,流量130g/s,總壓8.9×105Pa,總溫2 700K;采用2MW磁流體動(dòng)力加速器,通道長度96cm,出口尺寸3.6cm×3.6cm,電極65對(duì),磁場(chǎng)強(qiáng)度2T,電離種子為質(zhì)量分?jǐn)?shù)1%的鈉鉀合金(NaK),具備速度增加150%的潛力,加速器出口馬赫數(shù)為3.52,最高速度3 550m/s,出口靜壓2.8×104Pa,出口靜溫3 000K,運(yùn)行時(shí)間小于1s,并可通過第二噴管繼續(xù)膨脹加速,出口速度取決于試驗(yàn)段背壓。該風(fēng)洞可采用引射氮?dú)庑纬尚∮?×103Pa的背壓條件。

    圖2 NASA MSFC磁流體動(dòng)力加速風(fēng)洞試驗(yàn)設(shè)備示意

    (5)俄羅斯TsAGI磁流體動(dòng)力加速風(fēng)洞

    俄羅斯針對(duì)磁流體動(dòng)力加速風(fēng)洞的研究始于前蘇聯(lián)時(shí)期,研制了大型同軸電弧加熱型磁流體動(dòng)力風(fēng)洞[11],結(jié)構(gòu)示意如圖3所示。該設(shè)備輸入功率1MW,可以把空氣加熱到3 700K;電離種子采用質(zhì)量分?jǐn)?shù)1%的鈉鉀合金(NaK),通過位于第一噴管喉道前部的混合室注入;加速通道中布置有80對(duì)相互獨(dú)立的電極,長0.5m,擴(kuò)張角1°,磁感應(yīng)強(qiáng)度達(dá)2.5T;通道后還連接有第二噴管,使氣流進(jìn)一步加速,加速通道出口速度6.5km/s,第二噴管出口速度8.0km/s。該磁流體動(dòng)力系統(tǒng)采用熱沉設(shè)計(jì),運(yùn)行時(shí)間不大于1s,可滿足空氣動(dòng)力測(cè)量參數(shù)要求。

    研究表明,第二噴管內(nèi)的非平衡效應(yīng)較明顯,但可采用依賴于加速器出口組分濃度的凍結(jié)流模型進(jìn)行較好的預(yù)估,且種子的存在使得氮?dú)夥肿幼杂烧駝?dòng)能趨于減少,從而導(dǎo)致噴管出口速度的增加和靜溫的升高。

    1-電弧加熱器;2-混合室;3-種子注入;4-初級(jí)超音速噴管;5-磁流體加速器;6-二級(jí)噴管;7-試驗(yàn)段。

    (6)國外幾種風(fēng)洞的主要參數(shù)

    國外主要磁流體動(dòng)力加速風(fēng)洞主要參數(shù)如表1所示。

    表1 磁流體動(dòng)力加速風(fēng)洞主要參數(shù)

    Tab.1 Main parameters of MHD acceleration wind tunnel

    1.2 國內(nèi)磁流體動(dòng)力加速風(fēng)洞發(fā)展現(xiàn)狀

    國內(nèi)針對(duì)磁流體動(dòng)力加速風(fēng)洞研究開展較晚,研究?jī)?nèi)容則主要集中于磁流體加速器的仿真分析及實(shí)驗(yàn)驗(yàn)證等方面。

    文獻(xiàn)[12]進(jìn)行了磁流體動(dòng)力加速器通道的性能仿真研究,分段法拉第型加速通道長1m,入口面積0.2m2,來流馬赫數(shù)為2,總溫3 000K,總壓1×106Pa,假設(shè)電離度均勻且電導(dǎo)率為20S/m,當(dāng)磁感應(yīng)強(qiáng)度為6T,在電磁力作用下流體的速度增加18.1%,馬赫數(shù)增加31.9%,總溫增加31.0%,總壓增加32.7%。與普通的加熱方式相比,由磁流體動(dòng)力加速器加入的焓值中有一部分是以推力功的形式加入的,不會(huì)導(dǎo)致熵增。

    文獻(xiàn)[13]初步研制出基于激波風(fēng)洞的超聲速氣流磁流體動(dòng)力技術(shù)實(shí)驗(yàn)系統(tǒng),設(shè)計(jì)了馬赫數(shù)為2的超聲速噴管及實(shí)驗(yàn)段,采用氦氣驅(qū)動(dòng)氬氣,在平衡接觸面運(yùn)行方式下得到高溫氣體,通過在低壓段注入電離種子K2CO3粉末,實(shí)現(xiàn)高溫條件下導(dǎo)電流體的產(chǎn)生。得到了以下結(jié)論:通過合理設(shè)計(jì)激波管高低壓段長度,實(shí)驗(yàn)時(shí)間為9ms左右,出口總溫達(dá)到3 500K以上,馬赫數(shù)為2的超聲速氣流的電導(dǎo)率則可達(dá)40S/m。在此基礎(chǔ)上,文獻(xiàn)[13]還開展了超聲速氣流磁流體加速的初步實(shí)驗(yàn)研究,當(dāng)噴管入口總壓0.704 9MPa、理論平衡溫度8 372.8K、噴管出口馬赫數(shù)為1.5時(shí),在電容充電400V、磁感應(yīng)強(qiáng)度0.5T的條件下,超聲速氣流的電導(dǎo)率約為150S/m、電效率約為28%,采用電參數(shù)測(cè)試方法對(duì)磁流體加速效果進(jìn)行評(píng)估,速度增加約15.7%;研究還表明電導(dǎo)率對(duì)加速通道的電效率及加速效果等有很重要的影響。

    從研究進(jìn)展來看,國內(nèi)磁流體動(dòng)力加速風(fēng)洞尚處于關(guān)鍵技術(shù)研究階段,在磁流體動(dòng)力加速器原理驗(yàn)證及性能評(píng)估等方面取得了一定的研究成果,但研究方向不全面,研究深度不足,與實(shí)際工程應(yīng)用尚有較大差距。

    2 磁流體動(dòng)力加速風(fēng)洞原理

    磁流體動(dòng)力加速風(fēng)洞的基本原理是采用電弧或感應(yīng)耦合方式加熱,通過氣動(dòng)噴管加速混合室生成高溫高壓等離子體,然后利用磁流體加速器進(jìn)行電磁力學(xué)加速,在不改變混合室總壓的條件下,進(jìn)一步提高超聲速氣流的速度和總焓,從而獲得所需的試驗(yàn)氣流。

    20世紀(jì)50年代,在國際航空會(huì)議上提出了磁控氣動(dòng)力學(xué)的概念[14]。在此基礎(chǔ)上,各國針對(duì)超聲速氣流磁流體加速技術(shù)開展了大量研究,并首先在研制高焓高超聲速風(fēng)洞中開展了應(yīng)用[15-16]。

    磁流體動(dòng)力加速是一種電磁加速方式,其基本原理如圖4所示。圖中為電流密度;為磁感應(yīng)強(qiáng)度;為電場(chǎng)強(qiáng)度;為速度,在外加電場(chǎng)和磁場(chǎng)的耦合作用下,注入導(dǎo)電流體的電能轉(zhuǎn)化為動(dòng)能而使氣流速度增加。一般認(rèn)為磁流體主要受到洛倫茲效應(yīng)(力效應(yīng))和焦耳效應(yīng)(熱效應(yīng))兩種效應(yīng)的綜合影響。洛倫茲效應(yīng)表現(xiàn)為超聲速導(dǎo)電氣流在電場(chǎng)作用下與外界形成回路,并在磁場(chǎng)作用下受洛倫茲力而加速;焦耳效應(yīng)則主要表現(xiàn)為能量釋放的消極影響,即由于焦耳熱產(chǎn)生的逆壓梯度而給氣流帶來阻滯作用。理論研究表明,當(dāng)磁感應(yīng)強(qiáng)度為2T時(shí),焦耳熱僅能使氣體靜溫升高6K;實(shí)驗(yàn)研究也發(fā)現(xiàn),由于焦耳耗散大部分能量被用來激發(fā)氣體振動(dòng)能,焦耳加熱造成的影響可以忽略。

    圖4 磁流體動(dòng)力加速原理

    研究磁流體動(dòng)力加速風(fēng)洞時(shí)主要應(yīng)考慮加速通道結(jié)構(gòu)形式及采用的電離方式等方面。

    磁流體加速通道的常見結(jié)構(gòu)形式主要有線性霍爾型、分段法拉第型、連續(xù)斜聯(lián)型及導(dǎo)體斜壁型[17],如圖5所示,為斜角,w為壁面斜角。一般來說,霍爾型加速通道更適合于低密度的氣體,而分段法拉第型加速通道則更適合于高密度的氣體,連續(xù)斜聯(lián)型及導(dǎo)體斜壁型加速通道具備與分段法拉第型加速器相同的優(yōu)勢(shì),同時(shí)降低了系統(tǒng)復(fù)雜性,但其性能一般低于分段法拉第型通道。

    圖5 磁流體加速通道的主要結(jié)構(gòu)形式

    磁流體電離方式則主要包括熱電離及外部電離兩類。當(dāng)氣體溫度較高(高于2 500K)時(shí),采用熱電離方式,在氣體中加入堿金屬(堿金屬鹽)物質(zhì)作為電離種子,利用堿金屬電離電位較低的特點(diǎn), 在相對(duì)較低的溫度下獲得部分電離氣體;當(dāng)氣體溫度較低(低于2 500K)時(shí),采用外部電離方式,其電離需要外界能量的注入。熱電離方式大多基于激波管設(shè)備實(shí)現(xiàn);外部電離方式則包括電子束電離、納秒脈沖放電、微波放電等多種形式。

    3 磁流體動(dòng)力加速風(fēng)洞方案

    為滿足先進(jìn)天地往返飛行器地面試驗(yàn)高馬赫、潔凈空氣、正確化學(xué)成分、長時(shí)間運(yùn)行的要求,本文提出一種基于熱電離的磁流體動(dòng)力加速風(fēng)洞方案。通過加熱設(shè)備將空氣加熱至總溫0=3 000~3 500K,并在混合室加入0.1%~1%堿金屬銫Cs,再通過一次噴管將氣流速度提升至=1.5,隨后通過磁流體動(dòng)力加速通道進(jìn)行加速,磁場(chǎng)強(qiáng)度>5T,使得出口氣流速度達(dá)到=3.5,并利用二次噴管進(jìn)一步膨脹加速,在不小于0.04m2(200mm×200mm)的噴管出口范圍內(nèi),最高達(dá)到出口=15的模擬環(huán)境??傮w方案示意如圖6所示。

    圖6 磁流體動(dòng)力加速風(fēng)洞總體方案示意

    (1)氣體電離方式

    采用高頻等離子發(fā)生器為設(shè)備提供加熱源。高頻等離子體發(fā)生器中的工作氣體因擊穿而導(dǎo)電,并在交變磁場(chǎng)作用下形成電流,產(chǎn)生焦耳熱。條件適當(dāng)時(shí),放電將能維持下去并不斷使送入發(fā)生器的冷工作氣體加熱形成等離子體流。由于這種等離子體發(fā)生器沒有電極燒損所引起的污染問題,提升了等離子體潔凈度。加入電離種子提高氣體電導(dǎo)率,同時(shí),為實(shí)現(xiàn)正確空氣化學(xué)成分的模擬,應(yīng)盡量降低氣體離子化程度。根據(jù)氣動(dòng)熱力學(xué)原理,高溫空氣平衡成分隨溫度的變化有以下的趨勢(shì):隨著壓力的升高,離解反應(yīng)要推遲到較高的溫度才發(fā)生。相反,如果壓力減小,則離解反應(yīng)在較低溫度下就會(huì)發(fā)生。當(dāng)氣體總壓0<1×106Pa,控制氣體總溫0=3 000~3 500K,可使氣體離子化低于30%。

    (2)磁流體動(dòng)力加速通道形式

    為滿足高密度氣流加速需求,采用分段法拉第型加速通道。同時(shí),由于加速通道長時(shí)間運(yùn)行過程中將產(chǎn)生大量熱量,采用水冷循環(huán)對(duì)通道進(jìn)行主動(dòng)冷卻。

    (3)氣流參數(shù)測(cè)試方案

    由于磁流體動(dòng)力加速風(fēng)洞內(nèi)部工作條件惡劣,傳統(tǒng)接觸式測(cè)試分析手段難以應(yīng)用。因此,采用可調(diào)諧二極管激光吸收光譜技術(shù)(Tunable Diode Laser Absorption Spectroscopy,TDLAS),將適當(dāng)頻率的激光經(jīng)準(zhǔn)直后射入試驗(yàn)段的氣流中,并測(cè)定透射光譜線,利用譜線的線型信息隨溫度變化特性實(shí)現(xiàn)溫度測(cè)量,利用吸收率的絕對(duì)強(qiáng)度獲取氣體濃度信息。

    4 磁流體動(dòng)力加速風(fēng)洞關(guān)鍵技術(shù)問題

    磁流體動(dòng)力加速風(fēng)洞與傳統(tǒng)高超聲速風(fēng)洞相比,其優(yōu)勢(shì)在于獲得更高試驗(yàn)氣流馬赫數(shù)的同時(shí)避免噴管喉道前總溫過高,從而避免喉道高溫破壞及試驗(yàn)氣流污染。同時(shí),磁流體動(dòng)力加速風(fēng)洞系統(tǒng)復(fù)雜,超聲速氣流磁流體加速試驗(yàn)存在大量技術(shù)難題,相關(guān)試驗(yàn)主要由美俄等少數(shù)航空航天大國完成。超聲速氣流磁流體加速試驗(yàn)關(guān)鍵技術(shù)問題主要包括超聲速氣流電離技術(shù)、磁流體動(dòng)力加速通道設(shè)計(jì)、高超聲速模擬測(cè)試技術(shù)等方面。

    4.1 超聲速氣流電離技術(shù)

    應(yīng)用于超聲速流動(dòng)的磁流體加速技術(shù),其核心技術(shù)是超聲速氣流電離。難點(diǎn)問題在于如何使氣流具有并保持較好的電導(dǎo)率水平及分布。氣流的電導(dǎo)率與氣體種類、氣動(dòng)參數(shù)、通道尺寸等因素都有關(guān)。一般來說,實(shí)現(xiàn)氣體電離時(shí)所需的能量會(huì)隨著通道體積的增大而增加,需要根據(jù)實(shí)際條件設(shè)計(jì)合理的磁流體通道尺寸。同時(shí),對(duì)于熱電離而言,氣體電導(dǎo)率隨溫度(降低)或壓力(增加)都存在近似的指數(shù)衰減關(guān)系,需要綜合選取氣動(dòng)條件。針對(duì)綜合因素條件下超聲速氣流電離規(guī)律及機(jī)理仍需深入研究。

    此外,針對(duì)不同電離方式在磁流體加速風(fēng)洞的工程應(yīng)用,還需開展相關(guān)關(guān)鍵技術(shù)研究。

    對(duì)于熱電離方式,需重點(diǎn)發(fā)展超聲速氣流條件下的電離種子注入技術(shù)。不同的堿金屬電離種子, 其狀態(tài)不一樣, 使用方法也存在差異。固態(tài)種子可采用燃油增溶溶解方式,液態(tài)種子可采用壓力噴射混合方法。

    對(duì)于外部電離而言,從理論上分析,電子束電離方式是最有效的外部電離方法,但目前電子束電離方式在工程應(yīng)用研究中,也還存在諸如真空加速、透射窗力熱強(qiáng)度設(shè)計(jì)、電子束散射特性等一系列的技術(shù)問題需要解決。

    4.2 磁流體動(dòng)力加速通道設(shè)計(jì)

    磁流體動(dòng)力加速風(fēng)洞設(shè)計(jì)首先需要確定磁流體動(dòng)力加速通道的結(jié)構(gòu)形式。文獻(xiàn)顯示,國外多采用分段法拉第結(jié)構(gòu)形式,這是由工作溫度和壓力決定的。磁流體加速試驗(yàn)時(shí)一般需要提高磁感應(yīng)強(qiáng)度以增強(qiáng)流場(chǎng)的電磁效應(yīng),可以選用合適的電磁體設(shè)備。但對(duì)于法拉第型通道而言,當(dāng)流場(chǎng)密度較低時(shí),磁感應(yīng)強(qiáng)度的增大使得帶電粒子的回旋頻率增大,霍爾系數(shù)放大,通道內(nèi)會(huì)產(chǎn)生高度非均勻的磁流體流動(dòng),從而影響磁流體加速的性能。因此,分段法拉第通道的設(shè)計(jì)應(yīng)綜合考慮氣流密度及磁感應(yīng)強(qiáng)度等因素。

    同時(shí),磁流體加速通道中的電極設(shè)計(jì)也是試驗(yàn)中的關(guān)鍵技術(shù),其設(shè)計(jì)方案還需與具體的電離方式相適應(yīng)。例如,為了減小通道中的氣流擾動(dòng),磁流體動(dòng)力電極的型面應(yīng)設(shè)計(jì)成流線型;對(duì)于熱電離方式而言,需要應(yīng)考慮長時(shí)間運(yùn)行的結(jié)構(gòu)強(qiáng)度及高溫冷卻等問題。

    4.3 高超聲速模擬測(cè)試技術(shù)

    由于熱電離設(shè)備(電弧加熱器或高頻感應(yīng)加熱器)和磁流體動(dòng)力電極上均有外加電場(chǎng)和磁場(chǎng),往往發(fā)生瞬間放電的現(xiàn)象。在強(qiáng)磁場(chǎng)的作用下,易產(chǎn)生強(qiáng)電磁干擾,從而影響各參數(shù)的測(cè)量。往常工程中曾出現(xiàn)未經(jīng)屏蔽的高頻感應(yīng)加熱條件下的壓力傳感器損壞的情況。因此,磁流體動(dòng)力加速風(fēng)洞中參數(shù)測(cè)試時(shí)應(yīng)采用電磁屏蔽技術(shù),避免強(qiáng)電磁環(huán)境對(duì)測(cè)試傳感器的干擾。

    等離子體測(cè)試的重要參數(shù)包括加速通道中的電導(dǎo)率和速度、噴管出口速度和氣體組分濃度等,后續(xù)應(yīng)針對(duì)一系列測(cè)試技術(shù)開展研究:1)研究適用于加速通道中電導(dǎo)率測(cè)量的微波干涉儀技術(shù);2)研究適用于氣流溫度和速度測(cè)量的平面激光誘導(dǎo)熒光技術(shù);3)研究適用于高溫氣流組分、濃度和溫度等參數(shù)測(cè)量的高分辨率、高性能的光譜測(cè)試技術(shù)。

    5 磁流體動(dòng)力加速風(fēng)洞發(fā)展建議

    結(jié)合國內(nèi)外發(fā)展現(xiàn)狀,在后續(xù)研究中,須著重考慮以下幾個(gè)方面。

    (1)提前開展磁流體動(dòng)力加速風(fēng)洞研制規(guī)劃

    依托各類高超聲速飛行器研制任務(wù),我國在等離子加熱設(shè)備及試驗(yàn)技術(shù)上具備了一定的研究基礎(chǔ),形成了電弧加熱設(shè)備、高頻感應(yīng)加熱設(shè)備等應(yīng)用能力。為進(jìn)一步滿足真實(shí)氣體溫度、清潔空氣、大尺度、長時(shí)間、高馬赫數(shù)地面模擬試驗(yàn)需求,我國應(yīng)提前在相關(guān)領(lǐng)域開展研究規(guī)劃。針對(duì)磁流體動(dòng)力加速通道設(shè)計(jì)、超聲速氣流電離技術(shù)及高超聲速模擬測(cè)試技術(shù)三大關(guān)鍵技術(shù)問題,重點(diǎn)突破復(fù)雜電源系統(tǒng)技術(shù)、加熱通道主動(dòng)冷卻技術(shù)、熱電離技術(shù)、高分辨率高性能的光譜測(cè)試技術(shù),為磁流體動(dòng)力加速風(fēng)洞研制奠定技術(shù)基礎(chǔ)。

    (2)注重磁流體動(dòng)力學(xué)原理研究

    磁流體動(dòng)力學(xué)原理研究是開展磁流體動(dòng)力加速風(fēng)洞研制的基礎(chǔ)。磁流體加速涉及電場(chǎng)、磁場(chǎng)及流場(chǎng)之間的復(fù)雜耦合關(guān)系,物理過程復(fù)雜。特別是超聲速氣流電離、磁流體動(dòng)力加速通道設(shè)計(jì)等問題都依賴于磁流體動(dòng)力學(xué)原理研究基礎(chǔ)。目前國內(nèi)外針對(duì)該問題的研究尚不深入,其內(nèi)在機(jī)理尚不完全清晰,使得工程研制的理論依據(jù)不足。因此,磁流體動(dòng)力學(xué)原理研究是后續(xù)面臨的一個(gè)重要方向。一方面,發(fā)展磁流體動(dòng)力學(xué)多場(chǎng)耦合數(shù)值仿真模型,描述磁流體流動(dòng)的演化過程及其機(jī)理;另一方面,利用電弧加熱設(shè)備或高頻感應(yīng)加熱設(shè)備提供高溫實(shí)驗(yàn)氣體,建立磁流體動(dòng)力實(shí)驗(yàn)系統(tǒng),研究不同的通道結(jié)構(gòu)、工作介質(zhì)及電離方式下的流動(dòng)特性,為磁流體動(dòng)力加速風(fēng)洞的研制提供理論基礎(chǔ)。

    (3)拓展磁流體動(dòng)力技術(shù)多領(lǐng)域應(yīng)用

    隨著磁流體動(dòng)力技術(shù)成熟度的提高,除應(yīng)用于磁流體動(dòng)力加速風(fēng)洞外,該技術(shù)逐漸在沖壓發(fā)動(dòng)機(jī)推進(jìn)、高超聲速流動(dòng)控制以及機(jī)載設(shè)備供電等方面表現(xiàn)出廣闊的應(yīng)用前景。沖壓發(fā)動(dòng)機(jī)推進(jìn)方面,采用磁流體能量旁路的磁流體沖壓組合發(fā)動(dòng)機(jī),可獲得更高的飛行馬赫數(shù),提高發(fā)動(dòng)機(jī)性能;高超聲速流動(dòng)控制方面,利用磁流體動(dòng)力對(duì)高超聲速飛行器進(jìn)氣道流動(dòng)進(jìn)行調(diào)控,可改善進(jìn)氣道內(nèi)流場(chǎng)品質(zhì),拓寬發(fā)動(dòng)機(jī)工作范圍;機(jī)載設(shè)備供電方面,利用燃燒室后的高溫燃?xì)膺M(jìn)行磁流體發(fā)電,為長時(shí)間大功率機(jī)載設(shè)備供電提供新的技術(shù)途徑。因此,通過開展磁流體動(dòng)力加速風(fēng)洞的理論研究及工程化,突破關(guān)鍵技術(shù),拓展磁流體動(dòng)力技術(shù)在上述領(lǐng)域應(yīng)用,牽引相關(guān)學(xué)科發(fā)展,推動(dòng)顛覆性創(chuàng)新,具有重要意義。

    6 結(jié)束語

    磁流體動(dòng)力加速風(fēng)洞是拓展高超聲速地面模擬試驗(yàn)?zāi)芰Φ娜录夹g(shù)路線,是解決真實(shí)氣體溫度、清潔空氣、大尺度、長時(shí)間、高馬赫數(shù)地面模擬試驗(yàn)難題的重要手段。國外相關(guān)機(jī)構(gòu)在該領(lǐng)域開展了大量研究工作,完成了相關(guān)原理驗(yàn)證及方案設(shè)計(jì),并開展了很多工程研制工作。本文結(jié)合國內(nèi)外發(fā)展現(xiàn)狀,提出了基于熱電離的磁流體動(dòng)力加速風(fēng)洞方案,分析了超聲速氣流電離技術(shù)、磁流體動(dòng)力加速通道設(shè)計(jì)、高超聲速模擬測(cè)試技術(shù)等磁流體動(dòng)力加速風(fēng)洞關(guān)鍵技術(shù)問題,并闡述了開展磁流體動(dòng)力加速風(fēng)洞技術(shù)路徑規(guī)劃,注重磁流體動(dòng)力學(xué)原理研究,拓展磁流體動(dòng)力加速技術(shù)多領(lǐng)域應(yīng)用的發(fā)展建議。

    [1] PHILLIP J J, JOHN B P. The Hyper-X Launch Vehicle: Challenges and Design Considerations for Hypersonic Flight Testing[C]//AIAA/CIRA 13thInternational Space Planes & Hypersonics Systems & Technologies Conference. Capua, Italy, 2005.

    [2] RONDEAU C M, JORRIS T R. X-51A Scramjet Demonstrator Program: Waverider Ground & Flight Test[R]. Virginia: Armed Services Technical Information Agency Documents, 2013.

    [3] 羅金玲, 周丹, 康宏琳, 等. 典型氣動(dòng)問題試驗(yàn)方法研究的綜述[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2014, 32(2): 600-609. LUO Jinling, ZHOU Dan, KANG Honglin, et al. Summarization of Experimental Methods Associated with Typical Aerodynamic Issues[J]. Acta Aerodynamica Sinica, 2014, 32(2): 600-609. (in Chinese)

    [4] 張時(shí)空, 李江, 秦飛, 等. 兩級(jí)入軌運(yùn)載器RBCC動(dòng)力系統(tǒng)內(nèi)流道設(shè)計(jì)與性能計(jì)算[J]. 推進(jìn)技術(shù), 2015, 36(4): 520-526. ZHANG Shikong, LI Jiang, QIN Fei, et al. Design & Computational Investigation of a RBCC Propulsion System Flowpath for a TSTO Concept Vehicle[J]. Journal of Propulsion Technology, 2015, 36(4): 520-526. (in Chinese)

    [5] 弗蘭克K, 丹E馬倫. 先進(jìn)高超聲速試驗(yàn)設(shè)備[M]. 柳森, 黃訓(xùn)銘, 譯. 北京: 航空工業(yè)出版社, 2015: 23-33. FRANG K L, DAN E M. Advanced Hypersonic Test Facilities[M]. Translated by: LIU Seng, HUANG Xunming. Beijing: Aviation Industry Press, 2015: 23-33. (in Chinese)

    [6] MARION L L, GLOYD S, THOMAS F, et al. A Research Program for the Development of a True-temperature Mach[C]// 37thAerospace Sciences Meeting & Exhibit. Reno, NV, 2000.

    [7] THOMAS R B. The 20MW LORHO MHD Accelerator for Wind Tunnel Drive-Design, Construction &Critique[C]//30thPlasmadynamic & Lasers Conference. Norfolk, VA, 1999.

    [8] PATE S R, SILER L G, STALLINGS D W, et al. Development of an MHD-augmented, High Enthalpy, Shock Tunnel Facility[J] . AIAA Journal, 1974, 12(3): 289-297.

    [9] BEST J T, FETTERHOFF T P, LASTER M L, et al. RDHWT/MARIAH II Hypersonic Wind Tunnel Research Program Update[C]//10th AIAA/NAL-NASDA-ISAS International Space Planes & Hypersonic Systems & Technologies Conference. Kyoto, Japan, 2001.

    [10] LITCHFORD R J, COLE J W. Magnetohydrodynamic Augmented Propulsion Experiment: I. Performance Analysis & Design[C]//33rd Plasmadynamics & Lasers Conference. Maui, Hawaii, 2002.

    [11] ALFEROV V I. Current Status & Potentialities of Wind Tunnels with MHD Acceleration[J]. High Temperature, 2000, 38: 300-313.

    [12] 鄭小梅, 呂浩宇, 徐大軍, 等. MHD加速器模式磁控進(jìn)氣道的優(yōu)化設(shè)計(jì)[J]. 航空學(xué)報(bào), 2010, 31(2): 223-230. ZHENG Xiaomei, LYU Haoyu, XU Dajun, et al. Optimization of Accelerator Mode MHD Controlled Inlet[J]. Acta Aeronautica et Astronautica Sinica, 2010, 31(2): 223-230. (in Chinese)

    [13] 李益文, 李應(yīng)紅, 張百靈, 等. 超聲速氣流磁流體加速初步實(shí)驗(yàn)研究[J]. 力學(xué)學(xué)報(bào), 2012, 44(2): 238-244. LI Yiwen, LI Yinghong, ZHANG Bailing, et al. Preliminary Experimental Investigation on Supersonic Flow Magnetohydrodynamic (MHD) Acceleration[J]. Chinese Journal of Theoretical & Applied Mechanics, 2012, 44(2): 238-244. (in Chinese)

    [14] RESLER E L, SEARS W R. The Prospects for Magneto-aerodynamics[J]. Journal of the Aeronautical Sciences, 1958(25): 235-245.

    [15] 李楠. 磁流體加速的機(jī)理研究和實(shí)驗(yàn)系統(tǒng)設(shè)計(jì)[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2014. LI Nan. Mechanism Study of MHD Acceleration & Experimental System Design[D]. Harbin: Harbin Institute of Technology, 2014. (in Chinese)

    [16] 李益文, 張百靈, 李應(yīng)紅, 等. 磁流體動(dòng)力學(xué)在航空工程中的應(yīng)用與展望[J]. 力學(xué)進(jìn)展, 2017, 47(1): 452-502. LI Yiwen, ZHANG Bailing, LI Yinghong, et al. Review on the Application & Prospect of Magnetohydrodynamics in Aeronautical Engineering[J]. Advances in Mechanics, 2017, 47(1): 452-502. (in Chinese)

    [17] MATTHEW W T, CLARK W H, RON J L. Three-dimensional Numerical Modeling of Magnetohydrodynamic Augmented Propulsion Experiment[C]//46thAIAA Aerospace Sciences Meeting & Exhibit. Reno, Nevada, 2008.

    Research on Development of Magneto-hydro-dynamics Acceleration Wind Tunnel Technology

    ZUO Guang1QI Bin1OU Dongbin2

    (1 Beijing Institute of Manned Space System Engineering, China Academy of Space Technology, Beijing 100094, China)(2 Beijing Key Laboratory of Arc Plasma Application Equipment, China Academy of Aerospace Aerodynamics, Beijing 100074, China)

    According to the requirement of ground simulation test for hypersonic vehicle, the traditional test method is difficult to reflect the real gas temperature, clean air, large-scale, long-term and high Mach number simulation ability, and the magneto-hydro-dynamic acceleration wind tunnel provides a new technology route. Domestic and external development status of the magneto-hydro-dynamic acceleration wind tunnel is reviewed. The principle of magneto-hydro-dynamic acceleration wind tunnel is introduced. A magneto-hydro- dynamic acceleration wind tunnel based on thermoelectric ionization is proposed. High frequency plasma generator is used to provide heating source for the equipment, so as to avoid the pollution caused by electrode burning. The total temperature of the gas is controlled below 3 500 K, so that the gas ionization is less than 30%. Through two-stage acceleration, the simulating environment of 15 Mach at the outlet is achieved. Furthermore, the key technical problems of magneto-hydro-dynamic acceleration wind tunnel are analyzed. In the field of supersonic gas ionization, the key problems lie in the law and mechanism of supersonic gas ionization and the key technologies such as ionization seed injection and electron beam ionization. In the design of magneto-hydro-dynamic acceleration channel, the comprehensive influence of gas flow density and magnetic induction intensity and the design technology of electrodes are mainly considered. The key technology of hypersonic simulation test include electromagnetic shielding technology, microwave interferometer technology, planar laser induced fluorescence technology, high resolution and high performance spectral testing technology, etc. Finally, suggestions are provided for future development of magneto-hydro-dynamic acceleration wind tunnel.

    Magneto-hydro-dynamics; hypersonic; wind tunnel; ground simulation test; re-entry space vehicle

    R857

    A

    1009-8518(2018)06-0001-11

    10.3969/j.issn.1009-8518.2018.06.001

    左光,男,1971年生,畢業(yè)于天津大學(xué)及莫斯科航空學(xué)院,碩士,研究員,中國航天科技集團(tuán)學(xué)術(shù)技術(shù)帶頭人。研究方向?yàn)樘斓赝碉w行器氣動(dòng)設(shè)計(jì)及不同升阻比飛行器EDL技術(shù)。E-mail:lunar_cast@126.com。

    齊玢,男,1986年生,畢業(yè)于北京航空航天大學(xué)航空科學(xué)與工程學(xué)院,博士,工程師。研究方向?yàn)楹教炱骺傮w設(shè)計(jì)。E-mail:qionline@163.com。

    2018-08-25

    (編輯:劉穎)

    猜你喜歡
    研究
    FMS與YBT相關(guān)性的實(shí)證研究
    2020年國內(nèi)翻譯研究述評(píng)
    遼代千人邑研究述論
    視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    關(guān)于遼朝“一國兩制”研究的回顧與思考
    EMA伺服控制系統(tǒng)研究
    基于聲、光、磁、觸摸多功能控制的研究
    電子制作(2018年11期)2018-08-04 03:26:04
    新版C-NCAP側(cè)面碰撞假人損傷研究
    關(guān)于反傾銷會(huì)計(jì)研究的思考
    焊接膜層脫落的攻關(guān)研究
    電子制作(2017年23期)2017-02-02 07:17:19
    在线av久久热| 亚洲电影在线观看av| 国产野战对白在线观看| 黄色视频不卡| 桃红色精品国产亚洲av| 欧美另类亚洲清纯唯美| 757午夜福利合集在线观看| 久久久久亚洲av毛片大全| 亚洲,欧美精品.| 免费人成视频x8x8入口观看| 成人三级做爰电影| 亚洲午夜精品一区,二区,三区| 精品人妻1区二区| АⅤ资源中文在线天堂| 久久九九热精品免费| 午夜福利成人在线免费观看| 激情在线观看视频在线高清| 日本 av在线| 51午夜福利影视在线观看| 精品一区二区三区四区五区乱码| 91在线观看av| 嫩草影视91久久| 国产一区二区在线av高清观看| 成人手机av| www.熟女人妻精品国产| 亚洲成人免费电影在线观看| 免费在线观看完整版高清| 精品久久久久久久末码| 丰满人妻熟妇乱又伦精品不卡| 欧美激情久久久久久爽电影| 久久99热这里只有精品18| 亚洲av第一区精品v没综合| 国产亚洲欧美98| 热re99久久国产66热| 黄色 视频免费看| 在线免费观看的www视频| 1024视频免费在线观看| 黄片播放在线免费| 老司机深夜福利视频在线观看| 久久国产精品影院| 国产精品影院久久| 亚洲第一av免费看| 欧美激情 高清一区二区三区| 亚洲精品粉嫩美女一区| 中文字幕精品亚洲无线码一区 | 一a级毛片在线观看| 久久久精品欧美日韩精品| 久99久视频精品免费| 免费在线观看日本一区| 日韩欧美国产一区二区入口| 日韩有码中文字幕| 国产亚洲精品av在线| 97碰自拍视频| 午夜福利免费观看在线| x7x7x7水蜜桃| 亚洲av成人一区二区三| 99国产极品粉嫩在线观看| 午夜福利免费观看在线| 久久久久国产一级毛片高清牌| 亚洲一区二区三区不卡视频| av欧美777| 满18在线观看网站| 中文字幕最新亚洲高清| 少妇粗大呻吟视频| 日韩三级视频一区二区三区| 一本一本综合久久| 人人妻人人看人人澡| 日韩精品青青久久久久久| 亚洲无线在线观看| 成年女人毛片免费观看观看9| 国产亚洲精品一区二区www| 最好的美女福利视频网| 亚洲国产精品成人综合色| 老司机福利观看| 欧美国产精品va在线观看不卡| 日韩成人在线观看一区二区三区| 99久久无色码亚洲精品果冻| 两个人看的免费小视频| 亚洲av成人不卡在线观看播放网| 黑丝袜美女国产一区| 日韩一卡2卡3卡4卡2021年| 9191精品国产免费久久| 中文字幕人成人乱码亚洲影| 1024视频免费在线观看| 亚洲精品av麻豆狂野| ponron亚洲| 国产精品日韩av在线免费观看| 免费在线观看黄色视频的| 麻豆一二三区av精品| 91字幕亚洲| 国产精品日韩av在线免费观看| 久久久精品国产亚洲av高清涩受| 中文字幕精品免费在线观看视频| 亚洲片人在线观看| 在线观看免费日韩欧美大片| 亚洲av中文字字幕乱码综合 | 性欧美人与动物交配| 老司机深夜福利视频在线观看| 欧美在线黄色| 午夜免费激情av| 日韩大尺度精品在线看网址| av中文乱码字幕在线| 香蕉丝袜av| 老司机福利观看| av免费在线观看网站| 1024香蕉在线观看| 国产成人av教育| 91九色精品人成在线观看| 国产亚洲精品第一综合不卡| 1024视频免费在线观看| 久久国产亚洲av麻豆专区| 国产av在哪里看| 日本免费a在线| 成人亚洲精品av一区二区| 老司机在亚洲福利影院| 哪里可以看免费的av片| 精品无人区乱码1区二区| 韩国av一区二区三区四区| 国产野战对白在线观看| 亚洲精品久久国产高清桃花| 免费一级毛片在线播放高清视频| 丁香欧美五月| 99国产精品一区二区蜜桃av| 亚洲欧美日韩高清在线视频| 亚洲av成人不卡在线观看播放网| 99国产精品一区二区蜜桃av| 50天的宝宝边吃奶边哭怎么回事| 91成人精品电影| 欧美性猛交黑人性爽| 亚洲一区高清亚洲精品| 在线观看免费视频日本深夜| 国产亚洲欧美在线一区二区| 欧美中文综合在线视频| 12—13女人毛片做爰片一| a级毛片a级免费在线| 国产成人欧美在线观看| 又紧又爽又黄一区二区| 色综合站精品国产| 亚洲avbb在线观看| 每晚都被弄得嗷嗷叫到高潮| 一夜夜www| 色尼玛亚洲综合影院| 成人特级黄色片久久久久久久| 国产欧美日韩精品亚洲av| 久久久精品国产亚洲av高清涩受| 国产精品1区2区在线观看.| 黄片小视频在线播放| 久9热在线精品视频| 一级毛片高清免费大全| 国产又色又爽无遮挡免费看| 亚洲成av人片免费观看| 欧美乱色亚洲激情| 久久精品国产清高在天天线| xxxwww97欧美| 手机成人av网站| 亚洲熟妇中文字幕五十中出| 成年女人毛片免费观看观看9| 一级黄色大片毛片| or卡值多少钱| www.精华液| 妹子高潮喷水视频| 男女之事视频高清在线观看| 久久人妻av系列| 午夜视频精品福利| 啪啪无遮挡十八禁网站| 丁香六月欧美| 国产一区在线观看成人免费| 久99久视频精品免费| 99国产精品一区二区三区| 国产一卡二卡三卡精品| 国产单亲对白刺激| 成人av一区二区三区在线看| 色播亚洲综合网| 人妻丰满熟妇av一区二区三区| 2021天堂中文幕一二区在线观 | 欧美成人午夜精品| 午夜福利欧美成人| 亚洲狠狠婷婷综合久久图片| 欧美色视频一区免费| 日本一本二区三区精品| 一级a爱视频在线免费观看| 久久国产精品男人的天堂亚洲| 亚洲无线在线观看| 精品一区二区三区四区五区乱码| 一级毛片女人18水好多| 国产伦在线观看视频一区| 国产三级在线视频| 女性生殖器流出的白浆| 亚洲精华国产精华精| 美女扒开内裤让男人捅视频| 久久精品亚洲精品国产色婷小说| 自线自在国产av| 国产爱豆传媒在线观看 | 波多野结衣高清无吗| www.999成人在线观看| 国产色视频综合| 午夜久久久在线观看| 久久精品亚洲精品国产色婷小说| 国产99白浆流出| 亚洲成人国产一区在线观看| 免费人成视频x8x8入口观看| 99久久综合精品五月天人人| 精品国产国语对白av| 亚洲一区二区三区色噜噜| 亚洲av美国av| 欧美日韩一级在线毛片| 欧美另类亚洲清纯唯美| 日韩精品中文字幕看吧| 午夜福利18| 免费观看人在逋| 国产亚洲欧美98| 可以在线观看毛片的网站| 国产一级毛片七仙女欲春2 | 精品午夜福利视频在线观看一区| 国产精品乱码一区二三区的特点| 韩国精品一区二区三区| 久久精品aⅴ一区二区三区四区| 亚洲精品粉嫩美女一区| 黄片大片在线免费观看| 精品久久久久久成人av| 欧美日本视频| 国产一级毛片七仙女欲春2 | 成年女人毛片免费观看观看9| av电影中文网址| 无遮挡黄片免费观看| a级毛片a级免费在线| 国产成人av教育| 亚洲av片天天在线观看| 亚洲全国av大片| 国产成年人精品一区二区| 中亚洲国语对白在线视频| 欧美丝袜亚洲另类 | 麻豆一二三区av精品| 伊人久久大香线蕉亚洲五| 亚洲性夜色夜夜综合| 成人国语在线视频| 波多野结衣巨乳人妻| 自线自在国产av| 亚洲av成人av| www.999成人在线观看| 91国产中文字幕| bbb黄色大片| 九色国产91popny在线| 亚洲中文字幕一区二区三区有码在线看 | 国产高清有码在线观看视频 | 老鸭窝网址在线观看| 人妻久久中文字幕网| 男男h啪啪无遮挡| 国产精品综合久久久久久久免费| 欧美日韩中文字幕国产精品一区二区三区| 精品熟女少妇八av免费久了| 亚洲中文日韩欧美视频| 此物有八面人人有两片| 免费av毛片视频| 给我免费播放毛片高清在线观看| 国产激情欧美一区二区| ponron亚洲| 88av欧美| av超薄肉色丝袜交足视频| 亚洲国产中文字幕在线视频| 中文亚洲av片在线观看爽| 午夜免费观看网址| 久久久久国内视频| 国产片内射在线| 国产精品99久久99久久久不卡| 久久精品夜夜夜夜夜久久蜜豆 | 国产精品久久久久久亚洲av鲁大| 亚洲性夜色夜夜综合| www国产在线视频色| 美女国产高潮福利片在线看| 久久精品国产亚洲av香蕉五月| 亚洲精品久久国产高清桃花| 国内精品久久久久精免费| 999久久久国产精品视频| 亚洲人成77777在线视频| 成人国语在线视频| 精品福利观看| 欧美黑人欧美精品刺激| 亚洲人成电影免费在线| 一本大道久久a久久精品| 久久精品国产亚洲av香蕉五月| 制服诱惑二区| 精品卡一卡二卡四卡免费| 国产精品久久久久久精品电影 | 他把我摸到了高潮在线观看| 精品久久蜜臀av无| 在线视频色国产色| 欧美黄色淫秽网站| 国内少妇人妻偷人精品xxx网站 | 满18在线观看网站| 超碰成人久久| 国产精品精品国产色婷婷| 国产精品电影一区二区三区| 香蕉av资源在线| 99久久综合精品五月天人人| 成人三级黄色视频| 国产亚洲欧美98| 18禁观看日本| 大香蕉久久成人网| 亚洲三区欧美一区| 亚洲一码二码三码区别大吗| 亚洲最大成人中文| 久久精品国产亚洲av高清一级| 后天国语完整版免费观看| 亚洲国产欧美网| 久久精品影院6| 国产又黄又爽又无遮挡在线| 国产一卡二卡三卡精品| 一本精品99久久精品77| 一进一出好大好爽视频| 级片在线观看| 亚洲片人在线观看| 国产成人啪精品午夜网站| 久久久国产欧美日韩av| 亚洲一区高清亚洲精品| 午夜福利18| 老司机靠b影院| 一个人观看的视频www高清免费观看 | 69av精品久久久久久| 久久精品国产综合久久久| 久久久水蜜桃国产精品网| 久久午夜综合久久蜜桃| 亚洲激情在线av| 亚洲avbb在线观看| 久久精品aⅴ一区二区三区四区| 欧美中文日本在线观看视频| 亚洲成人精品中文字幕电影| 久久久国产成人精品二区| 韩国精品一区二区三区| 99热6这里只有精品| 久久婷婷成人综合色麻豆| 嫁个100分男人电影在线观看| 午夜激情av网站| 亚洲av熟女| 超碰成人久久| 日本a在线网址| 搡老熟女国产l中国老女人| 成人亚洲精品一区在线观看| 亚洲午夜理论影院| 久久中文字幕一级| 欧美绝顶高潮抽搐喷水| 黄频高清免费视频| 亚洲va日本ⅴa欧美va伊人久久| 听说在线观看完整版免费高清| 99国产精品99久久久久| 成人免费观看视频高清| 亚洲国产欧美网| 视频区欧美日本亚洲| 成人欧美大片| 女性生殖器流出的白浆| 亚洲av成人一区二区三| 免费看美女性在线毛片视频| 99国产精品一区二区三区| 久久中文字幕一级| 日韩国内少妇激情av| 91av网站免费观看| 免费在线观看视频国产中文字幕亚洲| 中文字幕另类日韩欧美亚洲嫩草| 又黄又粗又硬又大视频| 满18在线观看网站| 国产精品98久久久久久宅男小说| 午夜福利在线观看吧| 777久久人妻少妇嫩草av网站| 日韩av在线大香蕉| 国语自产精品视频在线第100页| 亚洲黑人精品在线| 变态另类成人亚洲欧美熟女| 国产成人影院久久av| 婷婷亚洲欧美| 亚洲国产欧美一区二区综合| 国产精品一区二区精品视频观看| 欧美成人一区二区免费高清观看 | 一本综合久久免费| 18禁美女被吸乳视频| 无限看片的www在线观看| 日韩欧美免费精品| 欧美 亚洲 国产 日韩一| 免费女性裸体啪啪无遮挡网站| 欧美乱色亚洲激情| 亚洲aⅴ乱码一区二区在线播放 | 国产精品香港三级国产av潘金莲| 国产精品永久免费网站| 免费在线观看完整版高清| 黑人欧美特级aaaaaa片| 一二三四社区在线视频社区8| 精品国产乱子伦一区二区三区| 亚洲久久久国产精品| 我的亚洲天堂| 19禁男女啪啪无遮挡网站| 久久国产亚洲av麻豆专区| 久久久国产精品麻豆| 亚洲国产欧美网| 极品教师在线免费播放| 欧美日韩瑟瑟在线播放| 久久精品91无色码中文字幕| 日本黄色视频三级网站网址| 久久午夜亚洲精品久久| 伊人久久大香线蕉亚洲五| 精品无人区乱码1区二区| 男女午夜视频在线观看| 特大巨黑吊av在线直播 | 国产成人av激情在线播放| 国产欧美日韩精品亚洲av| 日日爽夜夜爽网站| 欧美午夜高清在线| 女人被狂操c到高潮| 久久久久久九九精品二区国产 | 中国美女看黄片| av电影中文网址| 99在线人妻在线中文字幕| 国产成人啪精品午夜网站| 嫩草影视91久久| 天天一区二区日本电影三级| 欧美激情极品国产一区二区三区| 天堂影院成人在线观看| 美女大奶头视频| 在线观看免费午夜福利视频| 侵犯人妻中文字幕一二三四区| 国产精品av久久久久免费| 亚洲 欧美一区二区三区| 精品国产美女av久久久久小说| 首页视频小说图片口味搜索| 午夜福利高清视频| 俺也久久电影网| 久久久久国产精品人妻aⅴ院| 视频区欧美日本亚洲| 亚洲在线自拍视频| 国产高清激情床上av| 久久久久久亚洲精品国产蜜桃av| 亚洲国产欧美日韩在线播放| 俺也久久电影网| 国产黄片美女视频| 亚洲七黄色美女视频| 一区二区三区国产精品乱码| or卡值多少钱| 欧美日韩瑟瑟在线播放| 高清毛片免费观看视频网站| 国产精品综合久久久久久久免费| 最新在线观看一区二区三区| av欧美777| 久久国产精品人妻蜜桃| 亚洲精品粉嫩美女一区| 一本大道久久a久久精品| 亚洲国产精品成人综合色| 亚洲人成网站在线播放欧美日韩| 免费在线观看影片大全网站| 看片在线看免费视频| 香蕉av资源在线| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产毛片av蜜桃av| 黄色片一级片一级黄色片| 黄色 视频免费看| 久久性视频一级片| 一二三四在线观看免费中文在| 亚洲欧美一区二区三区黑人| 看黄色毛片网站| 十八禁网站免费在线| 天堂影院成人在线观看| 国产三级在线视频| 亚洲人成电影免费在线| 亚洲三区欧美一区| 两个人视频免费观看高清| 我的亚洲天堂| 欧美午夜高清在线| 亚洲国产欧洲综合997久久, | 悠悠久久av| 一本一本综合久久| 激情在线观看视频在线高清| 嫩草影院精品99| 成人一区二区视频在线观看| 久久99热这里只有精品18| 亚洲avbb在线观看| 亚洲国产精品sss在线观看| 非洲黑人性xxxx精品又粗又长| 欧美一区二区精品小视频在线| 久久九九热精品免费| 国产亚洲精品久久久久久毛片| 色播亚洲综合网| 婷婷亚洲欧美| 大型av网站在线播放| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av片天天在线观看| 久久久久久久久久黄片| 最近最新中文字幕大全电影3 | 91国产中文字幕| 欧美成人午夜精品| 一级片免费观看大全| 老熟妇乱子伦视频在线观看| 国产亚洲精品久久久久5区| 国产伦在线观看视频一区| 18禁美女被吸乳视频| 两性午夜刺激爽爽歪歪视频在线观看 | 韩国精品一区二区三区| 免费在线观看黄色视频的| 日韩精品青青久久久久久| 91麻豆av在线| 国产精品二区激情视频| 成人永久免费在线观看视频| 18禁裸乳无遮挡免费网站照片 | 亚洲五月婷婷丁香| 午夜免费成人在线视频| 高清在线国产一区| 亚洲无线在线观看| 国产黄片美女视频| 女警被强在线播放| 在线观看一区二区三区| 一级毛片女人18水好多| 亚洲精品中文字幕一二三四区| 一边摸一边做爽爽视频免费| 非洲黑人性xxxx精品又粗又长| 日韩精品免费视频一区二区三区| 国产99久久九九免费精品| 性色av乱码一区二区三区2| 亚洲人成网站在线播放欧美日韩| 免费在线观看成人毛片| 他把我摸到了高潮在线观看| 色综合婷婷激情| 一级片免费观看大全| 久久久久久九九精品二区国产 | 精品福利观看| 亚洲成人久久爱视频| 国产不卡一卡二| 日韩大尺度精品在线看网址| 51午夜福利影视在线观看| 精品一区二区三区视频在线观看免费| 欧美色欧美亚洲另类二区| 婷婷六月久久综合丁香| 亚洲中文日韩欧美视频| 麻豆av在线久日| 日韩精品青青久久久久久| 亚洲国产精品久久男人天堂| 久久久国产精品麻豆| 国产精品一区二区精品视频观看| 国产麻豆成人av免费视频| 国产主播在线观看一区二区| 99国产精品99久久久久| 最近最新免费中文字幕在线| 国产精品影院久久| 久久午夜综合久久蜜桃| 精品国产乱码久久久久久男人| 麻豆成人av在线观看| 国内精品久久久久久久电影| 国产免费av片在线观看野外av| 国产1区2区3区精品| 99久久无色码亚洲精品果冻| 成人国产一区最新在线观看| www.熟女人妻精品国产| 亚洲狠狠婷婷综合久久图片| 熟女电影av网| 欧美中文日本在线观看视频| 一二三四社区在线视频社区8| 国产精品乱码一区二三区的特点| 午夜精品久久久久久毛片777| 国产不卡一卡二| 国产激情偷乱视频一区二区| 日本 欧美在线| 好男人电影高清在线观看| 午夜成年电影在线免费观看| or卡值多少钱| 一区二区三区激情视频| 成年女人毛片免费观看观看9| 免费高清视频大片| 757午夜福利合集在线观看| 欧美色视频一区免费| 久久狼人影院| 麻豆一二三区av精品| 18美女黄网站色大片免费观看| 亚洲精品在线观看二区| 欧美大码av| 男女之事视频高清在线观看| 国产高清视频在线播放一区| 国产精品 欧美亚洲| 久久久精品欧美日韩精品| 亚洲三区欧美一区| 中出人妻视频一区二区| 久久久水蜜桃国产精品网| 亚洲色图av天堂| 免费在线观看亚洲国产| 欧美黄色片欧美黄色片| 欧美日韩瑟瑟在线播放| 国产麻豆成人av免费视频| 久久天堂一区二区三区四区| 无限看片的www在线观看| 欧美人与性动交α欧美精品济南到| 亚洲精品久久成人aⅴ小说| 亚洲av成人一区二区三| 91国产中文字幕| 美女午夜性视频免费| 日日摸夜夜添夜夜添小说| 中国美女看黄片| 波多野结衣av一区二区av| 午夜免费鲁丝| 久久亚洲真实| 最好的美女福利视频网| 国内精品久久久久精免费| 母亲3免费完整高清在线观看| 婷婷精品国产亚洲av在线| 亚洲在线自拍视频| 免费搜索国产男女视频| 女人高潮潮喷娇喘18禁视频| 国产成人av激情在线播放| 国产精品久久久久久人妻精品电影| 久久久久亚洲av毛片大全| 精品久久久久久久久久免费视频| 欧美在线一区亚洲| av中文乱码字幕在线| 国产成人av激情在线播放| 99久久无色码亚洲精品果冻| 亚洲精品一区av在线观看| 日韩大尺度精品在线看网址| 亚洲精华国产精华精| 亚洲成人精品中文字幕电影| 亚洲精品av麻豆狂野| 俺也久久电影网| 午夜福利欧美成人|