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

    超燃沖壓發(fā)動機模態(tài)轉(zhuǎn)換及推力突變實驗研究

    2021-03-20 03:42:24顧洪斌周芮旭李忠朋
    實驗流體力學 2021年1期
    關(guān)鍵詞:凹腔來流馬赫數(shù)

    連 歡, 顧洪斌,*, 周芮旭,2, 李 拓,2, 李忠朋,2

    1. 中國科學院力學研究所 高溫氣體動力學國家重點實驗室, 北京 100190;2. 中國科學院大學, 北京 101408

    0 引 言

    雙模態(tài)沖壓發(fā)動機具有幾何尺寸固定、單通道的特點,可有效減少重量和迎風面積,擴展了吸氣式?jīng)_壓發(fā)動機的工作范圍和飛行包線。自20世紀60年代Builder[1]提出雙模態(tài)沖壓發(fā)動機概念后,各國開展了大量理論研究、地面試驗及飛行試驗[2-6],通過隨飛行馬赫數(shù)變化調(diào)整供油規(guī)律進行燃燒室內(nèi)釋熱分布調(diào)控,已實現(xiàn)了雙模態(tài)沖壓發(fā)動機性能包線及穩(wěn)定工作包線內(nèi)發(fā)動機總壓恢復系數(shù)、燃燒室流動形式及發(fā)動機推力的平穩(wěn)過渡。然而亞燃及超燃工作模態(tài)下的發(fā)動機氣動特性、燃燒模式卻存在明顯差異[7],導致氣動性能失穩(wěn)、發(fā)動機推力突變的瞬態(tài)流動-燃燒機制仍不清楚,限制了飛行軌跡適應(yīng)性釋熱控制方法的復雜程度,影響飛行軌跡設(shè)計。因此,需要針對雙模態(tài)沖壓發(fā)動機燃燒室模型開展模態(tài)轉(zhuǎn)換及推力突變機制研究,發(fā)展預(yù)測及主動控制方法以擴展雙模態(tài)沖壓發(fā)動機穩(wěn)定工作包線。

    固定幾何構(gòu)型的雙模態(tài)沖壓發(fā)動機宏觀特性由熱力學基本定律描述[2],主要包括絕熱壓縮、等壓釋熱及膨脹釋熱物理過程。燃燒室通常由等截面隔離段及擴張型燃燒段組成,隔離段的主要作用是防止預(yù)燃激波串向上游傳播引發(fā)進氣道不啟動,燃燒室截面設(shè)計應(yīng)保證燃燒室出口流速為聲速或超聲速,采用擴張型尾噴管從而不受傳統(tǒng)亞燃沖壓發(fā)動機固定幾何喉道尾噴管對氣動性能要求的限制。據(jù)此設(shè)計的雙模態(tài)沖壓發(fā)動機可同時具有由理論沿程馬赫數(shù)表征的亞燃和超燃2種工作模態(tài),模態(tài)轉(zhuǎn)換過程可由H-K圖描述[8],主要由隔離段預(yù)燃激波串狀態(tài)定義亞燃和超燃工作模態(tài)[9-13]。為簡化雙模態(tài)沖壓發(fā)動機內(nèi)流形式,開展模態(tài)轉(zhuǎn)換及推力突變過程研究,本文基于單邊擴張的雙模態(tài)沖壓發(fā)動機燃燒室模型,采用凹腔回流區(qū)穩(wěn)焰及固定質(zhì)量流量的壁面橫向射流燃料噴注方式。

    從氣動特性角度分析,非燃燒狀態(tài)下,雙模態(tài)沖壓發(fā)動機亞燃及超燃工作模態(tài)內(nèi)流道流速均為超聲速。燃燒狀態(tài)下,亞燃工作模態(tài)通過燃燒釋熱形成反壓,當燃燒釋熱反壓超過邊界層抗反壓能力后向上游傳播,在隔離段中形成預(yù)燃激波串,使燃燒室入口來流減速至亞聲速,并在亞聲速氣流中組織燃燒;隨著飛行馬赫數(shù)增加,將來流減速至亞聲速并組織燃燒造成的總壓損失無法保證足夠的總壓恢復系數(shù)及發(fā)動機正推力,隔離段中預(yù)燃激波串的強度逐漸衰減,結(jié)構(gòu)形式由邊界層附著擬正激波變化至斜激波至耗散消失,燃燒室入口來流為超聲速,在超聲速氣流中組織燃燒。因此,雙模態(tài)沖壓發(fā)動機燃燒室內(nèi)涉及的主要流動現(xiàn)象可總結(jié)為:1) 燃燒段釋熱反壓上行及隔離段預(yù)燃激波串結(jié)構(gòu)演變[14-15];2) 燃燒室邊界層抗反壓及流動分離[16-19]。以上2種流動現(xiàn)象對發(fā)動機推力的貢獻相反。在亞燃及超燃工作模態(tài),以及通過控制供油位置使釋熱分布后移實現(xiàn)的亞燃至超燃模態(tài)平穩(wěn)過渡過程中,釋熱反壓與邊界層抗反壓作用保持動態(tài)平衡;而雙模態(tài)沖壓發(fā)動機模態(tài)轉(zhuǎn)換中導致推力突變的流動特征尚不明確,由此構(gòu)成本文主要研究內(nèi)容。

    國內(nèi)外實現(xiàn)模態(tài)轉(zhuǎn)換的地面試驗研究主要包括:冷流狀態(tài)下,通過機械/空氣節(jié)流產(chǎn)生反壓以獲得氣動特性匹配規(guī)律[15,20-23];燃燒狀態(tài)下,由于燃燒上行反壓和預(yù)燃激波串主要與釋熱量相關(guān),大量研究通過調(diào)節(jié)燃料的質(zhì)量流量和當量比實現(xiàn)亞燃和超燃工作模態(tài)轉(zhuǎn)換[24-28],其過程伴隨推力滯環(huán)等非線性特征[29];邊界層抗反壓能力由邊界層熱流輸運特性決定,主要受壁面熱流及溫度影響,國內(nèi)外研究人員通過被動改變壁面溫度實現(xiàn)模態(tài)轉(zhuǎn)換[30],對支板凹腔穩(wěn)焰燃燒室構(gòu)型模態(tài)轉(zhuǎn)換規(guī)律及時間序列定量方法已基本掌握[31-38]。過去,與動態(tài)飛行軌跡相關(guān)的研究僅通過改變來流總溫實現(xiàn)模態(tài)轉(zhuǎn)換及推力突變[7,31,39-42];Meng等[43]建設(shè)了來流連續(xù)可變的地面直連燃燒實驗臺,通過連續(xù)改變來流總溫總壓以及飛行馬赫數(shù),模擬了真實的動態(tài)飛行軌跡,獲得了模態(tài)平穩(wěn)過渡的基本規(guī)律。本文在Meng等[43]的研究基礎(chǔ)上,開展來流連續(xù)變化的地面直連燃燒試驗,對碳氫燃料雙模態(tài)沖壓發(fā)動機模態(tài)轉(zhuǎn)換及推力突變現(xiàn)象進行研究,重點討論碳氫燃料雙模態(tài)沖壓發(fā)動機模態(tài)轉(zhuǎn)換及推力突變過程中預(yù)燃激波串強度及流動特征變化規(guī)律。

    值得注意的是,國內(nèi)外地面試驗研究還包括通過改變噴注角度和噴注位置實現(xiàn)模態(tài)轉(zhuǎn)換[42,44-47],這涉及到橫向射流混合層演變強耦合物理化學過程,橫向射流混合層對發(fā)動機內(nèi)流氣動特性的影響主要包括:1) 橫向射流混合層氣動節(jié)流特性[30];2) 橫向射流混合層與燃燒釋熱的強耦合特性[48]。橫向射流自身產(chǎn)生的反壓對來流具有減速效果,表現(xiàn)為氣動節(jié)流特性。橫向射流混合層與燃燒釋熱的強耦合特性主要包括以下物理內(nèi)涵:1) 橫向射流氣相分子擴散及多相湍流輸運特性決定混合層時空演變,影響混合速率及混合效率;2) 混合層的當量比空間分布決定燃燒特性,影響點火、火焰?zhèn)鞑ヌ匦约搬専崧史植迹?) 燃燒釋熱反壓造成的預(yù)燃激波串與混合層界面強干擾,改變氣相分子擴散及多相湍流輸運過程,影響混合層時空演變。因此,本文關(guān)注碳氫燃料混合層演變特征及運動激波串與混合層界面干擾對雙模態(tài)沖壓發(fā)動機模態(tài)轉(zhuǎn)換及推力突變的影響規(guī)律。

    綜上所述,為闡明亞燃及超燃工作模態(tài)轉(zhuǎn)換過程中導致推力突變的瞬態(tài)流動-燃燒過程,本文開展來流馬赫數(shù)連續(xù)變化的推力突變模態(tài)轉(zhuǎn)換直連實驗研究,重點關(guān)注碳氫燃料雙模態(tài)沖壓發(fā)動機模態(tài)轉(zhuǎn)換及推力突變過程中瞬態(tài)流動特征,討論運動激波串與混合層界面干擾以及穩(wěn)焰模式火焰動力學與流動特征之間的相互作用規(guī)律。

    1 連續(xù)變馬赫數(shù)直連實驗臺

    1.1 連續(xù)變馬赫數(shù)直連實驗臺主體

    連續(xù)變馬赫數(shù)直連試驗臺主體由流量動態(tài)調(diào)節(jié)氣源、可調(diào)加熱器、可調(diào)噴管以及實驗?zāi)P腿紵医M成[43]。通過壓力調(diào)節(jié)閥實現(xiàn)空氣、氫氣、氧氣流量動態(tài)調(diào)節(jié),采用環(huán)形多孔噴注中心燃燒的燒氫補氧方式對來流進行加熱,通過氣源流量動態(tài)調(diào)節(jié)實現(xiàn)加熱器內(nèi)模擬總溫總壓實時可調(diào),可調(diào)噴管型面按照高低馬赫數(shù)雙設(shè)計點非對稱方法設(shè)計,由機電伺服機構(gòu)驅(qū)動改變喉道面積實現(xiàn)模擬飛行馬赫數(shù)實時調(diào)節(jié)??烧{(diào)加熱器和可調(diào)噴管均采用水冷熱防護。流量動態(tài)調(diào)節(jié)氣源、可調(diào)加熱器、可調(diào)噴管通過PLC工業(yè)控制系統(tǒng)進行精確時序控制,接受實驗控制臺操作指令,并根據(jù)指令執(zhí)行流量調(diào)節(jié)、伺服機構(gòu)運動等相應(yīng)動作,同時具有實驗臺主體氣源壓力、直連臺雙螺柱安裝型推力傳感器數(shù)據(jù)采集和處理能力。實驗臺模擬最大氣體流量2.5 kg/s,最高總溫1900 K,最高總壓4 MPa,最長工作時間30 s,模擬飛行馬赫數(shù)范圍4.5~6.5。

    1.2 診斷技術(shù)

    直連實驗臺推力傳感器選擇Omega LC203-500雙螺柱安裝稱重傳感器,采樣頻率為1 kHz,安裝于直連實驗臺前端滑動架與地面支架之間,發(fā)動機推力增量由加熱器和發(fā)動機共同工作時實測臺架推力與加熱器單獨工作時臺架推力的差值表示。低頻壓力測點布置在實驗段上壁面中心線,共60個。測壓接頭連接內(nèi)徑1 mm的銅質(zhì)導管至量程為0~690 kPa的壓力模塊處,導管長度約為20 cm;實驗前對測量零點進行了標定,實驗中忽略導管中的壓力損失。通過DTC Initium系統(tǒng)進行數(shù)據(jù)采集,采樣頻率為300 Hz。同時分別在隔離段入口、燃燒室入口上壁面、凹腔底部及凹腔后緣布設(shè)高頻壓力傳感器。傳感器采用Kulite XTL-190(M)高頻壓力傳感器,量程0.7 MPa,采樣頻率為100 kHz。燃燒場測量實驗采用Phantom v1612高速相機,并在鏡頭前部安裝中心波長為430 nm、帶寬±10 nm的定制窄帶濾光片(峰值透過率為0.882)記錄燃燒場CH(A)*自由基的高速攝影圖像,相機快門曝光時間為20 μs,幀頻為6000幀/s,分辨率1280 pixel×800 pixel,記錄時間為10 s。實驗采用常規(guī)Z型紋影布局。利用氙燈冷光源及Phantom v2640高速相機拍攝紋影圖像,相機快門曝光時間為2 μs,幀頻為6000幀/s,分辨率1920 pixel×700 pixel,記錄時間為10 s,與對CH(A)*高速攝像的Phantom v1612相機通過DG645-#1延時發(fā)生器進行同步。低頻壓力測量系統(tǒng)、高頻壓力測量系統(tǒng)、CH(A)*以及紋影高速攝影相機均經(jīng)DG645-#2與直連臺主體PLC工業(yè)控制系統(tǒng)共同按照實驗控制臺操作指令觸發(fā)以實現(xiàn)設(shè)備運行與測量時序同步。

    實驗?zāi)P腿紵覟榻孛?0 mm×40 mm的等截面隔離段和單邊擴張雙凹腔燃燒室模型,其中隔離段長度為337 mm,燃燒室模型下壁面擴張角為2°,凹腔位于隔離段入口下游447 mm處,長度為65 mm,深度為17 mm,傾角為22.5°。凹腔上游60.5和9.5 mm處分別安裝煤油噴塊和引導氫噴塊。此外,凹腔底部安裝火花塞輔助點火。實驗段燃燒室模型及壓力傳感器布設(shè)方式如圖1所示,其中pref,p1,p2,p3分別對應(yīng)隔離段入口、隔離段出口、燃燒室入口、第二凹腔前緣低頻壓力測點位置,紅色框線為石英觀察窗口安裝位置。

    圖1 雙模態(tài)模型燃燒室Fig. 1 Dual-mode model combustor

    2 來流連續(xù)變化的模態(tài)轉(zhuǎn)換及推力突變實驗

    來流連續(xù)變化的地面試驗?zāi)M飛行馬赫數(shù)為5.0~6.0,對應(yīng)燃燒室入口馬赫數(shù)為2.5~3.0。為簡化混合層耦合特性及燃料供應(yīng)系統(tǒng)對實驗結(jié)果的影響,在來流變化過程中,采用固定質(zhì)量流量和固定位置壁面垂直噴注RP-3航空煤油,燃料當量比和通量比僅隨來流馬赫數(shù)變化。模擬飛行高度等實驗工況(ABCD)見圖2。本文討論2種加速上行實驗工況(ABC、ADC)以及2種減速下行實驗工況(CBA、CDA)的發(fā)動機性能、非定常流動-燃燒特性及影響規(guī)律。實驗參數(shù)如表1所示。

    圖2 試驗工況Fig.2 Experimental conditions

    表1 試驗工況Table 1 Experimental parameters

    實驗工況A、B、C分別對應(yīng)固定馬赫數(shù)工作狀態(tài)Ma=5.0、5.6和6.0,固定來流條件下的發(fā)動機沿程壓力、沿程馬赫數(shù)、推力曲線等基本工作特性,以及流動、混合層演變、火焰動力學特征用作開展來流連續(xù)變化模態(tài)轉(zhuǎn)換及推力突變過程的研究基準。以實驗工況ABC為例,表示在10 s穩(wěn)定實驗時間內(nèi),加熱氣體模擬總溫總壓按照AB~BC段線形變化,燒氫補氧加熱器各路氣體(空氣、氫氣、氧氣)流量在AB~BC固定實驗狀態(tài)理論計算值基礎(chǔ)上線形擬合構(gòu)成分段函數(shù)。為模擬飛行馬赫數(shù)在AB~BC段線性變化,拉瓦爾噴管型面及喉道高度由伺服機構(gòu)在AB~BC按線性擬合分段函數(shù)作動,通過加熱氣體總溫總壓調(diào)節(jié)及拉瓦爾噴管型面變化的精確時序配合實現(xiàn)加減速來流連續(xù)變化的實驗工況。典型10 s加速上行實驗工況ABC及ADC的加熱器實驗總壓曲線如圖3所示,加熱器5.0~10.0 s形成穩(wěn)定冷態(tài)空氣,10.0 s時加熱器點火通過燒氫補氧形成高焓空氣,14.5~24.5 s為模擬動態(tài)飛行路徑實驗,各路氣源流量以及噴管型面根據(jù)模擬路徑需求進行動態(tài)調(diào)節(jié),發(fā)動機于14.5 s通過火花點火器和引導氫點火。

    圖3 可調(diào)加熱器總壓Fig. 3 Total pressure of the transient operation heater

    3 結(jié)果及討論

    3.1 直連實驗推力曲線

    2種加速上行實驗工況(ABC、ADC)以及2種減速下行實驗工況(CBA、CDA)的直連實驗燃燒室推力曲線(臺架實測推力與相應(yīng)工況加熱器推力差值)如圖4所示。

    推力傳感器采樣率為1 kHz,最小時間解析為1 ms,由圖4(a)可見ABC與ADC上行路徑均發(fā)生推力突變現(xiàn)象,其中ABC路徑推力突變發(fā)生于t4,ABC=18.986 s,ADC路徑推力突變發(fā)生于t3,ADC=21.809 s。圖4(b) 表示減速下行工況CBA及CDA,推力持續(xù)增加,且下行路徑對推力特性具有明顯影響。

    圖4 推力傳感器測量數(shù)據(jù)Fig.4 Thrust measurement during simulated acceleration and deceleration

    本文研究所用單邊擴張燃燒室熱壅塞煤油當量比上限約為0.4,實驗時煤油質(zhì)量流量不變,當量比根據(jù)來流變化改變約為0.2~0.3,引導氫當量比約為0.1,因此,在引導氫與煤油共同燃燒階段(14.5~16.0 s),引導氫釋熱量對推力特性影響較明顯,推力高于煤油穩(wěn)定燃燒階段(16.0~24.5 s)。在不同飛行路徑條件下火花點火器頻率為15 Hz,單脈沖能量約為50 J,通過點火器電壓電流實時測量,確保點火功率穩(wěn)定。同時,引導氫噴注壓力實時測量,確保引導氫流量穩(wěn)定。ABC、ADC和CBA、CDA飛行路徑的主要區(qū)別是在煤油穩(wěn)定燃燒階段(16.0~24.5 s)由加熱器和噴管聯(lián)動實現(xiàn)動態(tài)飛行軌跡的飛行動壓差異。

    3.2 時間離散沿程馬赫數(shù)計算

    時間離散質(zhì)量加權(quán)沿程馬赫數(shù)計算方法描述如下:首先,對模擬飛行馬赫數(shù)變化曲線進行時間以及馬赫數(shù)離散,得到各個時刻隔離段入口所對應(yīng)的馬赫數(shù);其次,將各個時刻隔離段入口馬赫數(shù)作為初始條件代入定常一維質(zhì)量加權(quán)馬赫數(shù)計算程序中,得到該時刻的馬赫數(shù)沿程分布;最后,將各個時刻的沿程馬赫數(shù)沿時間軸排列,得到馬赫數(shù)時空分布圖。本文針對連續(xù)可變飛行馬赫數(shù)5.0~6.0,采用300個時間離散點(同沿程壓力傳感器采樣率)和三次Hermite內(nèi)插值方法,將計算得到的100個空間離散沿程壓力數(shù)據(jù)進行時間離散的沿程馬赫數(shù)計算。質(zhì)量加權(quán)沿程馬赫數(shù)的時空分布演變?nèi)鐖D5所示。

    圖5(a)~(c)分別表示固定馬赫數(shù)工作狀態(tài)A、B和C(Ma=5.0,5.6和6.0)。14.5 s對應(yīng)非燃燒狀態(tài),燃燒室內(nèi)為超聲速氣流;14.5~16.0 s為輔助點火引導氫與煤油共同燃燒階段,形成預(yù)燃激波對來流進行減速,在雙凹腔內(nèi)以亞聲速組織燃燒;16.0~19.5 s為煤油穩(wěn)定燃燒階段,A、B工況均為亞燃工作模態(tài),在第一凹腔內(nèi)以亞聲速組織燃燒,C工況為超燃工作模態(tài),全流道超聲速組織燃燒。

    圖5(d)表示加速上行ABC路徑。14.5~16.0 s,輔助點火引導氫與煤油共同燃燒形成預(yù)燃激波,燃燒室入口為亞聲速;從16.0 s起為煤油穩(wěn)定燃燒階段,在第一凹腔內(nèi)以亞聲速組織燃燒,至17.050 s左右第一凹腔內(nèi)燃燒狀態(tài)發(fā)生亞燃-超燃工作模態(tài)轉(zhuǎn)換,第一凹腔內(nèi)以超聲速組織燃燒,在第二凹腔內(nèi)穩(wěn)焰以亞聲速組織燃燒,18.960 s左右第二凹腔燃燒狀態(tài)發(fā)生工作模態(tài)轉(zhuǎn)換。由上節(jié)可知推力突變發(fā)生于18.986 s,說明推力突變可能與雙凹腔燃燒組織方式相關(guān)。

    圖5(e)表示加速上行ADC路徑。ADC路徑14.5~16.0 s,輔助點火引導氫與煤油共同燃燒形成預(yù)燃激波,燃燒室入口為亞聲速; 從16.0 s起為煤油穩(wěn)定燃燒階段,在第一凹腔內(nèi)以亞聲速組織燃燒,至21.750 s左右第一凹腔內(nèi)燃燒狀態(tài)發(fā)生亞燃-超燃工作模態(tài)轉(zhuǎn)換,第一凹腔及第二凹腔均以超聲速組織燃燒。由上節(jié)可知推力突變發(fā)生于21.809 s。

    圖5(f)和(g)表示減速下行工況CBA及CDA,全程為亞燃工作模態(tài)。

    圖5 質(zhì)量加權(quán)馬赫數(shù)Fig. 5 Mass weighted average Ma during simulated acceleration and deceleration

    3.3 模態(tài)轉(zhuǎn)換及推力突變過程的瞬態(tài)流動特征

    針對ABC路徑第一、第二凹腔模態(tài)轉(zhuǎn)換以及推力突變時刻,在17.052~17.061 s、18.965~18.967 s、18.986 s之后時段內(nèi)燃燒室入口至第一凹腔的紋影以及CH*自發(fā)光的高速攝影圖像如圖6所示。其中,17.052~17.061 s對應(yīng)由沿程馬赫數(shù)定義的亞燃燃燒模態(tài),動態(tài)流動現(xiàn)象為預(yù)燃激波串由擬正激波向斜激波變化,強度減弱,隨著來流加速過程,預(yù)燃激波串持續(xù)后移并伴隨小范圍振蕩。預(yù)燃激波向下游移動造成運動激波與混合層界面干擾,降低了橫向射流穿透深度,削弱了橫向射流氣動節(jié)流作用,同時改變了混合效率。如圖6(b)中CH*自發(fā)光圖像所示,釋熱區(qū)面積隨之減少,穩(wěn)焰模式由尾跡凹腔穩(wěn)焰逐漸過渡至剪切層穩(wěn)焰。18.965~18.967 s對應(yīng)由沿程馬赫數(shù)定義的超燃燃燒模態(tài),可見第二凹腔上游維持一定的燃燒反壓;隨來流馬赫數(shù)持續(xù)增大至18.986 s時不足以維持燃燒反壓,發(fā)生推力突變,推力突變后燃燒仍然穩(wěn)定,穩(wěn)焰模式為剪切層穩(wěn)焰。18.986 s之后流動特征均與圖6中18.986 s流動特征相同。

    圖6 (a) ABC加速上行軌跡紋影;(b) ABC加速上行軌跡CH*自發(fā)光Fig. 6 (a) Schlieren imaging during acceleration ABC; (b) CH* chemiluminescence during acceleration ABC

    針對ADC路徑模態(tài)轉(zhuǎn)換及推力突變時刻,在21.750~21.759 s時段以及21.809 s內(nèi)燃燒室入口至第一凹腔的紋影以及CH*自發(fā)光的高速攝影圖像如圖7所示。與圖6中ABC路徑相同,預(yù)燃激波串由擬正激波向斜激波變化,強度減弱,隨著來流加速過程,預(yù)燃激波串持續(xù)后移并伴隨小范圍振蕩。預(yù)燃激波向下游移動造成運動激波與混合層界面干擾,降低了橫向射流穿透深度,且與圖6中ABC路徑相同。如圖7中CH*自發(fā)光圖像所示,假設(shè)CH*自發(fā)光強度與釋熱量正相關(guān),則在加速軌跡的亞燃-超燃模態(tài)轉(zhuǎn)換過程中,釋熱量逐漸減少。21.758~21.809 s可見超燃模態(tài)轉(zhuǎn)換及推力突變對應(yīng)的流動特征;21.809 s推力突變之后的流動特征均與推力突變時刻相同。

    圖7 (a) ADC加速上行軌跡紋影;(b) ADC加速上行軌跡CH*自發(fā)光Fig.7 (a) Schlieren imaging during acceleration ADC; (b) CH* chemiluminescence during acceleration ADC

    ABC與ADC路徑模態(tài)轉(zhuǎn)換及推力突變過程具有相同的瞬態(tài)流動發(fā)展過程,根據(jù)圖6及7所示ADC與ABC路徑的共性動態(tài)紋影圖像特征,將來流變化導致的亞燃-超燃模態(tài)轉(zhuǎn)換及推力突變現(xiàn)象中的瞬態(tài)流動過程總結(jié)為以下4個階段:

    第I階段:隔離段預(yù)燃激波階段,以亞聲速組織燃燒,尾跡穩(wěn)焰及凹腔穩(wěn)焰,加速上行ABC與ADC路徑亞燃-超燃模態(tài)轉(zhuǎn)換均發(fā)生于第I階段;

    第II階段:運動激波與混合層界面干擾階段,以亞聲速組織燃燒,氣動反壓能力逐漸減弱,尾跡穩(wěn)焰及剪切層穩(wěn)焰;

    第III階段:燃燒反壓維持階段,亞聲速組織燃燒向超聲速組織燃燒過渡,剪切層穩(wěn)焰;

    第IV階段:超聲速核心流與附面層階段,超聲速組織燃燒,剪切層穩(wěn)焰,加速上行ABC與ADC路徑推力突變均發(fā)生于第IV階段。

    減速下行工況CBA及CDA,全程為亞燃工作模態(tài),流動特征及穩(wěn)焰模式與圖7第I階段隔離段預(yù)燃激波亞燃模態(tài)完全相同。

    3.4 模態(tài)轉(zhuǎn)換及推力突變過程流動機制初步討論

    雙模態(tài)沖壓發(fā)動機燃燒分區(qū)規(guī)律通常在固定來流總壓條件下由當量比和總溫界定,射流通量影響暫不討論。相應(yīng)的燃燒組織方式分別對應(yīng)幾種穩(wěn)焰模式,主要包括尾跡穩(wěn)焰、凹腔穩(wěn)焰、剪切層穩(wěn)焰以及幾種穩(wěn)焰模式之間的振蕩模式。本文實驗?zāi)P腿紵业娜紵謪^(qū)規(guī)律如圖8所示,其中熱壅塞邊界由Rayleigh分析方法估計給出,實驗中變總壓條件下的A、B和C工況位于亞燃與超燃理論工作邊界,分別對應(yīng)尾跡穩(wěn)焰、凹腔穩(wěn)焰、剪切層穩(wěn)焰模式。

    圖8 A、B和C工況燃燒模式Fig. 8 Combustion regime of case A, B and C

    釋熱總量和亞燃與超燃理論工作邊界是本文討論流動機制的基礎(chǔ),首先在絕熱無黏假設(shè)條件下討論隔離段預(yù)燃激波演變特性。隔離段出口及入口壓比p1/pref可表示為[2]:

    (1)

    超聲速核心流與隔離段入口截面積比值可表示為:

    (2)

    理論上,超聲速核心流與隔離段入口截面積比值與壓比為拋物線關(guān)系,如圖9所示。圖9(a)和(b)顯示了路徑ABC及ADC亞燃到超燃模態(tài)轉(zhuǎn)換過程中,Ma=5.0時預(yù)燃激波尚未形成,隔離段預(yù)燃激波形成后由正激波向斜激波強度衰減至耗散消失。圖9(c)和(d)則顯示了路徑CBA和CDA亞燃模態(tài),Ma=6.0時預(yù)燃激波尚未形成,建立預(yù)燃激波后激波強度先發(fā)生衰減,之后逐漸向正激波強度過渡。首次衰減是引導氫停止導致釋熱減少造成的。

    該分析方法給出了超聲速核心流截面與預(yù)燃激波強度的相對關(guān)系,可以用于說明加速及減速過程中隔離段預(yù)燃激波結(jié)構(gòu)的演化規(guī)律。根據(jù)上述分析可知,釋熱總量與內(nèi)流道匹配是模態(tài)轉(zhuǎn)換及推力變化過程的根本,隔離段預(yù)燃激波是主要的流動特征。加速上行ABC和ADC路徑發(fā)生模態(tài)轉(zhuǎn)換和推力突變的根本原因是亞燃與超燃理論工作邊界釋熱總量差別;減速下行CBA和CDA未發(fā)生模態(tài)轉(zhuǎn)換及推力突變是由于點火引導氫釋熱導致釋熱量差別不大。

    由于上述分析方法不考慮熱流邊界層作用,故此處通過加速上行ABC和ADC路徑沿程壓力與隔離段入口靜壓比的時空演變關(guān)系,討論模態(tài)轉(zhuǎn)換及推力突變過程瞬態(tài)流動特征,如圖10所示。實驗段燃燒室模型及壓力傳感器布設(shè)方式見圖1,其中pref,p1,p2,p3分別對應(yīng)隔離段入口、隔離段出口、燃燒室入口第一凹腔前緣、第二凹腔前緣低頻壓力測點位置。

    加速上行ABC及ADC路徑連續(xù)變化來流引起的亞燃-超燃模態(tài)轉(zhuǎn)換發(fā)生在第I階段(ABC:17.050 s,ADC:21.750 s)。從圖6和7可以觀察到主導流動特征是預(yù)燃激波結(jié)構(gòu)演變。相應(yīng)的,圖10(a)和(b)定量表示亞燃-超燃模態(tài)轉(zhuǎn)換過程中由靜壓比表征的預(yù)燃激波串強度逐漸減弱,且隔離段預(yù)燃激波串位置后移。且圖10(a)所示為ABC路徑雙凹腔燃燒組織影響模態(tài)轉(zhuǎn)換規(guī)律,由圖10(b)可見ADC路徑僅在第一凹腔組織燃燒。推力突變發(fā)生于第IV階段(ABC:18.986 s,ADC:21.809 s),圖6和7表示推力突變時內(nèi)流道已觀察不到明顯的預(yù)燃激波結(jié)構(gòu)。相應(yīng)的,圖10(a)和(b) 定量表示推力突變發(fā)生于加速上行ABC和ADC路徑沿程壓比均大幅度降低之后。

    燃燒室入口第二凹腔前緣、第二凹腔前緣低頻壓力測點位置與隔離段入口靜壓比(p2/pref,p3/pref)截面如圖11所示,所描述及結(jié)論與圖10一致。

    圖9 動態(tài)飛行路徑超聲速核心流截面與預(yù)燃激波強度分析Fig. 9 Impulse function analysis during simulated acceleration and deceleration

    圖10 (a) ABC沿程壓比; (b) ADC沿程壓比Fig. 10 (a) Pressure ratio during acceleration ABC; (b) Pressure ratio during acceleration ADC

    圖11 (a) ABC及ADC路徑截面壓比p2/pref; (b) ABC及ADC路徑截面壓比p3/prefFig. 11 (a) p2/pref cutoff with route ABC and ADC; (b) p3/pref cutoff with route ABC and ADC

    釋熱總量與內(nèi)流道匹配是模態(tài)轉(zhuǎn)換及推力變化過程的根本,主導流動特征是隔離段預(yù)燃激波強度演變特性。然而燃料橫向射流氣動節(jié)流以及釋熱反壓在隔離段預(yù)燃激波削弱耗散之后,仍然可對來流進行減速并維持推力。

    動態(tài)飛行軌跡相關(guān)氣動熱及燃燒熱積分效應(yīng),可改變熱流邊界層特性以及發(fā)動機內(nèi)流道抗反壓能力,造成亞燃與超燃工作邊界變化。熱壁面模態(tài)擾動及湍流輸運機制復雜[54-55],基于目前研究的局限性,變動壓條件下釋熱反壓與邊界層抗反壓的動態(tài)失穩(wěn)可能是導致瞬態(tài)推力突變的流動機制。如圖12所示,動態(tài)飛行軌跡氣動熱和燃燒熱帶來的熱流邊界層效應(yīng)可以改變抗反壓能力,造成亞燃與超燃工作邊界的變化,因此在進行飛發(fā)一體軌跡優(yōu)化設(shè)計時應(yīng)予以關(guān)注。目前,由于缺乏反應(yīng)流邊界層演化的有效定量實驗測量手段以及壁面熱流時間積分效應(yīng)邊界層演化理論預(yù)測方法[56-57],故尚無法提供足夠的實驗證據(jù)和理論解釋動態(tài)飛行軌跡帶來的氣動熱及燃燒熱問題對亞燃與超燃工作邊界的定量影響。未來可參考熱湍流邊界層研究,通過定量化紋影技術(shù)以及壁面溫度、熱流測量技術(shù)輔助認知反應(yīng)流邊界層演化規(guī)律。

    圖12 熱流邊界層影響工作邊界示意圖Fig. 12 Illustration of mode transition shift due to boundary layer disturbance

    4 結(jié) 論

    本文基于單邊擴張的雙模態(tài)沖壓發(fā)動機燃燒室模型,采用雙凹腔回流區(qū)穩(wěn)焰及固定質(zhì)量流量的壁面橫向射流燃料噴注方式,開展來流連續(xù)變化飛行馬赫數(shù)5.0~6.0加速上行和減速下行的地面直連試驗。結(jié)論如下:

    1) 通過地面實驗觀察到加速上行過程中來流變化可導致亞燃-超燃工作模態(tài)轉(zhuǎn)換及推力突變,且沖壓發(fā)動機推力性能受飛行路徑強烈影響。

    2) 加速上行模態(tài)轉(zhuǎn)換及推力突變過程的瞬態(tài)流動特征包括以下4個階段:第I階段,隔離段預(yù)燃激波階段,以亞聲速組織燃燒,尾跡穩(wěn)焰及凹腔穩(wěn)焰;第II階段,運動激波與混合層界面干擾階段,以亞聲速組織燃燒,尾跡穩(wěn)焰及剪切層穩(wěn)焰;第III階段,燃燒反壓維持階段,亞聲速組織燃燒向超聲速組織燃燒過渡,剪切層穩(wěn)焰;第IV階段,超聲速核心流與附面層階段,超聲速組織燃燒,剪切層穩(wěn)焰。加速上行ABC及ADC路徑亞燃-超燃模態(tài)轉(zhuǎn)換均發(fā)生于第I階段,推力突變均發(fā)生于第IV階段。

    3) 釋熱總量與內(nèi)流道匹配是模態(tài)轉(zhuǎn)換及推力變化過程的根本,主導流動特征是隔離段預(yù)燃激波強度演變特性,然而燃料橫向射流氣動節(jié)流以及釋熱反壓在隔離段預(yù)燃激波削弱耗散之后,仍然可對來流進行減速并維持推力。同時,動態(tài)飛行軌跡相關(guān)氣動熱及燃燒熱積分效應(yīng),可改變熱流邊界層特性以及發(fā)動機內(nèi)流道抗反壓能力,造成亞燃與超燃工作邊界變化。這在飛發(fā)一體軌跡優(yōu)化設(shè)計時應(yīng)予以關(guān)注。然而由于目前缺乏反應(yīng)流邊界層演化的有效定量實驗測量手段以及壁面熱流時間積分效應(yīng)邊界層演化理論預(yù)測方法,故尚無法給出氣動熱及燃燒熱問題對亞燃與超燃工作邊界的定量影響規(guī)律。

    猜你喜歡
    凹腔來流馬赫數(shù)
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    凹腔對高超聲速邊界層穩(wěn)定性的影響
    兩種典型來流條件下風力機尾跡特性的數(shù)值研究
    能源工程(2022年2期)2022-05-23 13:51:48
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    載荷分布對可控擴散葉型性能的影響
    不同來流條件對溢洪道過流能力的影響
    縫翼凹腔擋板氣動性能和降噪效果數(shù)值研究
    壁面噴射當量比對支板凹腔耦合燃燒的影響
    彈發(fā)匹配驗證試驗系統(tǒng)來流快速啟動技術(shù)研究
    柴油機活塞燃燒室凹腔邊緣的非線性有限元強度分析方法
    欧美成人一区二区免费高清观看| 中文字幕av成人在线电影| 国产精品av视频在线免费观看| 国产精品99久久久久久久久| 国产欧美亚洲国产| 97精品久久久久久久久久精品| 午夜亚洲福利在线播放| av免费在线看不卡| 亚洲精品aⅴ在线观看| 亚洲一级一片aⅴ在线观看| 天堂中文最新版在线下载 | .国产精品久久| 一级av片app| 黄色怎么调成土黄色| 久久久久精品久久久久真实原创| 国产毛片a区久久久久| 好男人在线观看高清免费视频| 大片电影免费在线观看免费| 在线看a的网站| 日韩一区二区视频免费看| 亚洲怡红院男人天堂| 亚洲精品自拍成人| 久久久国产一区二区| 新久久久久国产一级毛片| 欧美亚洲 丝袜 人妻 在线| 最近手机中文字幕大全| 亚洲av不卡在线观看| 久久韩国三级中文字幕| 又大又黄又爽视频免费| 色播亚洲综合网| 久久久久久国产a免费观看| 国产av不卡久久| 中国三级夫妇交换| 亚洲内射少妇av| 亚洲国产欧美人成| 免费看a级黄色片| 国产乱人视频| 国产真实伦视频高清在线观看| 亚洲精品久久久久久婷婷小说| 欧美97在线视频| 天天躁日日操中文字幕| 97人妻精品一区二区三区麻豆| 精品一区二区三卡| 国产男女内射视频| 中文资源天堂在线| 久久久国产一区二区| 成人欧美大片| 国产黄片视频在线免费观看| 身体一侧抽搐| 伊人久久国产一区二区| 一区二区av电影网| 熟女av电影| 99热这里只有精品一区| 国产综合精华液| 精品国产一区二区三区久久久樱花 | 国产精品偷伦视频观看了| 直男gayav资源| 日产精品乱码卡一卡2卡三| 一个人观看的视频www高清免费观看| 秋霞伦理黄片| 黑人高潮一二区| 国产爽快片一区二区三区| 精品熟女少妇av免费看| 久久人人爽人人爽人人片va| 下体分泌物呈黄色| 欧美激情国产日韩精品一区| 亚洲国产精品专区欧美| 国产成人aa在线观看| 欧美日韩综合久久久久久| 新久久久久国产一级毛片| xxx大片免费视频| 精品熟女少妇av免费看| 精品久久久久久久人妻蜜臀av| 亚洲国产精品成人久久小说| 精品一区二区免费观看| 九草在线视频观看| 国产精品一区www在线观看| 我要看日韩黄色一级片| av在线亚洲专区| 一本久久精品| 国产精品秋霞免费鲁丝片| 亚洲无线观看免费| 97精品久久久久久久久久精品| 有码 亚洲区| 国内精品美女久久久久久| 新久久久久国产一级毛片| 美女内射精品一级片tv| 视频区图区小说| 国产高清三级在线| 欧美xxxx性猛交bbbb| 99久久中文字幕三级久久日本| 夫妻性生交免费视频一级片| 欧美丝袜亚洲另类| 爱豆传媒免费全集在线观看| 三级男女做爰猛烈吃奶摸视频| 欧美日韩亚洲高清精品| 国产高清有码在线观看视频| 日产精品乱码卡一卡2卡三| 18+在线观看网站| 成人毛片60女人毛片免费| 国产一区有黄有色的免费视频| 最后的刺客免费高清国语| 久久久久精品久久久久真实原创| 天天一区二区日本电影三级| 午夜精品国产一区二区电影 | 在线看a的网站| 中文字幕久久专区| 亚洲精品第二区| 亚洲精品国产成人久久av| 男女边摸边吃奶| 欧美潮喷喷水| 午夜福利视频1000在线观看| 精品一区二区免费观看| 亚洲电影在线观看av| 一级av片app| 亚洲,一卡二卡三卡| 国产一区二区在线观看日韩| 两个人的视频大全免费| 精品人妻一区二区三区麻豆| 国产黄片视频在线免费观看| 又大又黄又爽视频免费| 国产熟女欧美一区二区| 赤兔流量卡办理| 久久精品熟女亚洲av麻豆精品| 久久久久久国产a免费观看| 十八禁网站网址无遮挡 | 国产精品无大码| 国产一区有黄有色的免费视频| 久久久久精品性色| 久久99热这里只有精品18| .国产精品久久| 国产成人午夜福利电影在线观看| 一级毛片久久久久久久久女| 麻豆久久精品国产亚洲av| 亚洲精品一二三| 可以在线观看毛片的网站| 亚洲不卡免费看| 热99国产精品久久久久久7| 视频中文字幕在线观看| 精品一区二区免费观看| 成人亚洲欧美一区二区av| 别揉我奶头 嗯啊视频| 久久97久久精品| 成人国产av品久久久| 亚洲内射少妇av| 天天躁夜夜躁狠狠久久av| 国产精品一区二区三区四区免费观看| 久久精品国产亚洲网站| 超碰av人人做人人爽久久| 国产成人福利小说| 黄片无遮挡物在线观看| 十八禁网站网址无遮挡 | 亚洲精品第二区| 亚洲自偷自拍三级| 亚洲国产日韩一区二区| 久久99热这里只频精品6学生| 色综合色国产| 国模一区二区三区四区视频| 日本黄大片高清| 毛片一级片免费看久久久久| 人人妻人人爽人人添夜夜欢视频 | 美女xxoo啪啪120秒动态图| 欧美高清成人免费视频www| 亚洲欧洲国产日韩| 免费黄频网站在线观看国产| 久久热精品热| 联通29元200g的流量卡| 97在线人人人人妻| 香蕉精品网在线| 韩国高清视频一区二区三区| 国产黄色视频一区二区在线观看| 自拍偷自拍亚洲精品老妇| 亚洲精品aⅴ在线观看| 欧美激情国产日韩精品一区| av在线老鸭窝| 日本欧美国产在线视频| 午夜亚洲福利在线播放| 国产精品秋霞免费鲁丝片| 肉色欧美久久久久久久蜜桃 | 免费看a级黄色片| 国产黄片美女视频| 国产午夜精品一二区理论片| 亚洲精品国产av蜜桃| 99热这里只有是精品在线观看| 在线观看免费高清a一片| 亚洲内射少妇av| 国产欧美另类精品又又久久亚洲欧美| 高清av免费在线| 国产男人的电影天堂91| 久久99蜜桃精品久久| 精品熟女少妇av免费看| 免费人成在线观看视频色| av黄色大香蕉| 成人美女网站在线观看视频| 秋霞伦理黄片| 80岁老熟妇乱子伦牲交| 99热这里只有精品一区| 91久久精品国产一区二区成人| 久久久久国产网址| 在线免费观看不下载黄p国产| 日韩大片免费观看网站| 成年版毛片免费区| 亚洲成人一二三区av| 精品少妇久久久久久888优播| 看十八女毛片水多多多| 在线免费观看不下载黄p国产| 国产高清国产精品国产三级 | 国产精品99久久久久久久久| 欧美成人a在线观看| 久久国产乱子免费精品| kizo精华| 王馨瑶露胸无遮挡在线观看| 69人妻影院| 国产精品一及| 在线播放无遮挡| 99热6这里只有精品| 国产高清国产精品国产三级 | a级毛片免费高清观看在线播放| 2022亚洲国产成人精品| 两个人的视频大全免费| 亚洲自偷自拍三级| 六月丁香七月| 国产在视频线精品| 热99国产精品久久久久久7| 一本一本综合久久| 建设人人有责人人尽责人人享有的 | 亚洲av.av天堂| 亚洲欧美成人综合另类久久久| 九九在线视频观看精品| 一区二区三区精品91| 国产精品无大码| 国产一区二区三区av在线| 黄片wwwwww| 日韩制服骚丝袜av| 3wmmmm亚洲av在线观看| 一边亲一边摸免费视频| 男女国产视频网站| 久久久精品欧美日韩精品| 国产精品女同一区二区软件| 欧美日韩视频高清一区二区三区二| 成人午夜精彩视频在线观看| 欧美日韩综合久久久久久| 国产亚洲5aaaaa淫片| 欧美bdsm另类| 色哟哟·www| 天天躁日日操中文字幕| 熟妇人妻不卡中文字幕| 午夜免费观看性视频| 国产欧美亚洲国产| freevideosex欧美| 久久精品久久久久久噜噜老黄| 综合色av麻豆| 热99国产精品久久久久久7| 中文字幕亚洲精品专区| 99热这里只有是精品50| 18禁裸乳无遮挡动漫免费视频 | 亚洲一区二区三区欧美精品 | 国产精品国产三级专区第一集| 中文在线观看免费www的网站| 视频中文字幕在线观看| 国产亚洲91精品色在线| 欧美bdsm另类| 街头女战士在线观看网站| 永久网站在线| 插逼视频在线观看| 看免费成人av毛片| 18+在线观看网站| 免费黄色在线免费观看| 亚洲最大成人中文| 另类亚洲欧美激情| 免费少妇av软件| 直男gayav资源| 久久精品久久久久久久性| 久久久久久久大尺度免费视频| 天堂俺去俺来也www色官网| 午夜视频国产福利| 午夜免费观看性视频| 波野结衣二区三区在线| 精品一区二区三卡| 亚洲,欧美,日韩| 18+在线观看网站| 伦理电影大哥的女人| 国产精品国产三级国产专区5o| 国产大屁股一区二区在线视频| 国产精品一区二区性色av| 日韩不卡一区二区三区视频在线| 国产亚洲最大av| 久久综合国产亚洲精品| 高清av免费在线| 亚洲欧美日韩另类电影网站 | 亚洲欧美成人精品一区二区| 中文字幕制服av| 不卡视频在线观看欧美| 成人毛片60女人毛片免费| 少妇丰满av| 最新中文字幕久久久久| 欧美一区二区亚洲| 精品一区二区三区视频在线| 欧美xxxx黑人xx丫x性爽| 成人美女网站在线观看视频| 欧美最新免费一区二区三区| 精品熟女少妇av免费看| 又粗又硬又长又爽又黄的视频| 色婷婷久久久亚洲欧美| 久久精品国产亚洲av天美| 成人一区二区视频在线观看| 国产熟女欧美一区二区| 久久久久久九九精品二区国产| 亚洲精品亚洲一区二区| 国国产精品蜜臀av免费| 卡戴珊不雅视频在线播放| 特大巨黑吊av在线直播| a级一级毛片免费在线观看| 欧美日韩精品成人综合77777| 美女xxoo啪啪120秒动态图| 国产黄色视频一区二区在线观看| 18禁动态无遮挡网站| 久久影院123| 亚洲精品乱久久久久久| 日韩av免费高清视频| 欧美国产精品一级二级三级 | 亚洲人成网站在线观看播放| 观看免费一级毛片| 一级爰片在线观看| 欧美高清成人免费视频www| 在线天堂最新版资源| 日韩人妻高清精品专区| 久久国产乱子免费精品| 各种免费的搞黄视频| 激情 狠狠 欧美| 卡戴珊不雅视频在线播放| 亚洲欧美清纯卡通| 天天躁日日操中文字幕| 在现免费观看毛片| 一级二级三级毛片免费看| 国产成人a∨麻豆精品| 国产成人精品婷婷| 成人特级av手机在线观看| 久久久精品欧美日韩精品| 欧美日韩视频精品一区| 亚洲国产精品专区欧美| 一区二区三区四区激情视频| 在现免费观看毛片| 欧美日韩亚洲高清精品| 秋霞伦理黄片| 在线a可以看的网站| 深夜a级毛片| 汤姆久久久久久久影院中文字幕| 2022亚洲国产成人精品| 99热国产这里只有精品6| 亚洲av男天堂| 精品一区在线观看国产| 嫩草影院入口| 亚洲国产色片| 日本黄色片子视频| 人妻少妇偷人精品九色| 性色av一级| 精品99又大又爽又粗少妇毛片| 国产人妻一区二区三区在| 日韩一区二区三区影片| 欧美日韩亚洲高清精品| 国产亚洲av片在线观看秒播厂| 午夜精品一区二区三区免费看| 成人午夜精彩视频在线观看| 亚洲,欧美,日韩| 老司机影院毛片| 男女啪啪激烈高潮av片| 精品人妻视频免费看| 亚洲国产精品专区欧美| 男女边吃奶边做爰视频| 听说在线观看完整版免费高清| 日韩大片免费观看网站| 91久久精品国产一区二区三区| 久久精品久久精品一区二区三区| 亚洲高清免费不卡视频| 国产精品99久久99久久久不卡 | 国产黄频视频在线观看| 七月丁香在线播放| av一本久久久久| 毛片一级片免费看久久久久| 国产人妻一区二区三区在| 国产91av在线免费观看| 中文字幕久久专区| 国产探花极品一区二区| 一级毛片黄色毛片免费观看视频| 国产黄频视频在线观看| 亚洲综合色惰| 欧美一区二区亚洲| 亚洲最大成人手机在线| 性色av一级| 国产男女内射视频| 国产男女内射视频| 免费观看性生交大片5| 国产亚洲精品久久久com| 青春草国产在线视频| 欧美成人一区二区免费高清观看| 黄色怎么调成土黄色| 少妇 在线观看| 夜夜爽夜夜爽视频| 国产视频内射| 欧美成人午夜免费资源| 国产免费一级a男人的天堂| 国产亚洲91精品色在线| 少妇熟女欧美另类| av福利片在线观看| 2018国产大陆天天弄谢| 免费观看性生交大片5| 亚洲精品国产成人久久av| 可以在线观看毛片的网站| 亚洲国产欧美在线一区| 日日啪夜夜爽| 在线看a的网站| 午夜福利视频精品| 91在线精品国自产拍蜜月| av卡一久久| 干丝袜人妻中文字幕| 亚洲丝袜综合中文字幕| 国内精品美女久久久久久| 黄色欧美视频在线观看| 大又大粗又爽又黄少妇毛片口| 十八禁网站网址无遮挡 | 日韩av在线免费看完整版不卡| 激情 狠狠 欧美| 国产亚洲91精品色在线| 日韩强制内射视频| 国产高清有码在线观看视频| 久热久热在线精品观看| 成年免费大片在线观看| 午夜老司机福利剧场| 久久久久性生活片| 日韩视频在线欧美| 涩涩av久久男人的天堂| 国产成人一区二区在线| 国产老妇女一区| 日日啪夜夜撸| 中国三级夫妇交换| 久久综合国产亚洲精品| 男插女下体视频免费在线播放| 日本欧美国产在线视频| 国产精品久久久久久精品古装| 特级一级黄色大片| 亚洲精华国产精华液的使用体验| 婷婷色麻豆天堂久久| 亚洲精品色激情综合| 成人国产av品久久久| 亚洲欧美日韩无卡精品| 久久ye,这里只有精品| 51国产日韩欧美| 欧美+日韩+精品| 只有这里有精品99| 一区二区三区精品91| 天美传媒精品一区二区| 性色av一级| 人妻一区二区av| 亚洲成人精品中文字幕电影| 国产精品av视频在线免费观看| 国产亚洲午夜精品一区二区久久 | 国产精品人妻久久久影院| 国精品久久久久久国模美| 最近手机中文字幕大全| av在线app专区| 亚洲色图av天堂| 少妇丰满av| 午夜福利在线观看免费完整高清在| 日韩伦理黄色片| 美女xxoo啪啪120秒动态图| 亚洲精品中文字幕在线视频 | 人妻夜夜爽99麻豆av| 永久免费av网站大全| 久久99精品国语久久久| 大片免费播放器 马上看| 欧美97在线视频| 亚洲欧美一区二区三区黑人 | 99热这里只有是精品50| 国产精品偷伦视频观看了| 18禁裸乳无遮挡免费网站照片| 国产又色又爽无遮挡免| eeuss影院久久| 热99国产精品久久久久久7| 99re6热这里在线精品视频| 最近最新中文字幕免费大全7| 久久精品夜色国产| 网址你懂的国产日韩在线| av一本久久久久| 高清在线视频一区二区三区| 男人舔奶头视频| 91久久精品国产一区二区成人| 成人综合一区亚洲| 少妇丰满av| 秋霞在线观看毛片| 国产av码专区亚洲av| 中国国产av一级| 男女那种视频在线观看| 中文欧美无线码| 天天躁日日操中文字幕| 国产午夜精品一二区理论片| 99久久精品热视频| 久久精品国产亚洲网站| 插逼视频在线观看| 亚洲av一区综合| 日本一本二区三区精品| 深夜a级毛片| 亚洲国产色片| 国产白丝娇喘喷水9色精品| 丰满乱子伦码专区| 亚洲精品乱码久久久v下载方式| 男女啪啪激烈高潮av片| 欧美xxxx黑人xx丫x性爽| 人人妻人人爽人人添夜夜欢视频 | 18禁在线无遮挡免费观看视频| 人妻 亚洲 视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品99又大又爽又粗少妇毛片| 精品熟女少妇av免费看| 国产av国产精品国产| 欧美高清成人免费视频www| 久久久久久久久大av| 中文字幕免费在线视频6| 日韩人妻高清精品专区| 国模一区二区三区四区视频| 黄色配什么色好看| 高清av免费在线| 新久久久久国产一级毛片| 亚洲成人中文字幕在线播放| 国产 精品1| 亚洲第一区二区三区不卡| 欧美激情在线99| 国产又色又爽无遮挡免| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品久久久久久精品古装| 亚洲av日韩在线播放| 亚洲国产精品成人综合色| 人人妻人人看人人澡| 亚洲精品色激情综合| 中文字幕免费在线视频6| 日日啪夜夜撸| 中文乱码字字幕精品一区二区三区| www.av在线官网国产| 久久久精品欧美日韩精品| 国产日韩欧美在线精品| 欧美日韩视频精品一区| 熟女人妻精品中文字幕| av天堂中文字幕网| 国产综合精华液| 亚洲在线观看片| 免费观看a级毛片全部| 国产精品久久久久久精品电影小说 | 国产乱人视频| 亚洲欧美清纯卡通| 亚洲成色77777| 久久韩国三级中文字幕| 人妻制服诱惑在线中文字幕| 啦啦啦中文免费视频观看日本| 精品酒店卫生间| 精品久久久精品久久久| 寂寞人妻少妇视频99o| 人妻少妇偷人精品九色| 午夜福利在线在线| 蜜桃亚洲精品一区二区三区| 日韩伦理黄色片| 真实男女啪啪啪动态图| 自拍欧美九色日韩亚洲蝌蚪91 | 久久久精品94久久精品| 亚洲在久久综合| 91久久精品国产一区二区成人| 久久女婷五月综合色啪小说 | av线在线观看网站| 我的女老师完整版在线观看| 久久久a久久爽久久v久久| 18禁裸乳无遮挡动漫免费视频 | 国产男人的电影天堂91| 大又大粗又爽又黄少妇毛片口| 国产男人的电影天堂91| 又黄又爽又刺激的免费视频.| 大码成人一级视频| 精品久久久噜噜| 综合色av麻豆| 天天一区二区日本电影三级| av在线播放精品| 一区二区三区四区激情视频| 亚洲性久久影院| 嫩草影院新地址| 国产又色又爽无遮挡免| 国产精品精品国产色婷婷| 午夜免费观看性视频| 三级经典国产精品| 亚洲婷婷狠狠爱综合网| 精品国产一区二区三区久久久樱花 | 男女那种视频在线观看| 男女边吃奶边做爰视频| 日韩成人伦理影院| 国产高清有码在线观看视频| 黄色一级大片看看| 最近中文字幕高清免费大全6| 伊人久久精品亚洲午夜| 性插视频无遮挡在线免费观看| 成年版毛片免费区| 亚洲精品影视一区二区三区av| 久久女婷五月综合色啪小说 | 成人二区视频| 在线观看av片永久免费下载| 久久精品夜色国产| 只有这里有精品99| 亚洲成人一二三区av| 精品人妻视频免费看| 特大巨黑吊av在线直播| 国产真实伦视频高清在线观看| 免费看av在线观看网站| 精品人妻视频免费看| 中文字幕制服av| 最新中文字幕久久久久| 直男gayav资源| 亚洲精品乱久久久久久| 国产免费又黄又爽又色| 综合色av麻豆|