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

    塞式噴管對旋轉(zhuǎn)爆震起爆、傳播與推進特性的影響

    2022-10-29 11:58:52朱亦圓趙明皓王致程焦中天
    燃燒科學(xué)與技術(shù) 2022年5期
    關(guān)鍵詞:爆震周期性燃燒室

    朱亦圓,王 可,趙明皓,王致程,焦中天,范 瑋

    (1.西北工業(yè)大學(xué)動力與能源學(xué)院,西安 710072;2.西北工業(yè)大學(xué)陜西省航空動力系統(tǒng)熱科學(xué)重點實驗室,西安 710129;3.西安航天動力研究所液體火箭發(fā)動機技術(shù)重點實驗室,西安 710100)

    不同于緩燃波以亞聲速傳播,爆震波通過激波和燃燒波的耦合,以超聲速在未燃反應(yīng)物中傳播[1].正因其傳播速度快、反應(yīng)速率高,爆震燃燒接近等容過程,故基于爆震燃燒的動力裝置具有熱循環(huán)效率高的潛在優(yōu)點[2].旋轉(zhuǎn)爆震發(fā)動機(rotating detonation engine,簡稱RDE)作為一種可行的爆震推進方案,近年來已成為推進領(lǐng)域的研究熱點之一[3].

    早在20 世紀(jì)60 年代,Voitsekhoviskii 在圓盤形實驗裝置中發(fā)現(xiàn)了C2H2/O2的旋轉(zhuǎn)爆震現(xiàn)象[4].隨后,國內(nèi)外學(xué)者圍繞旋轉(zhuǎn)爆震開展了大量研究.為實現(xiàn)旋轉(zhuǎn)爆震的穩(wěn)定傳播,圍繞燃燒室特征尺寸[5-6]、燃料和氧化劑噴注器結(jié)構(gòu)[7-8]、燃料的種類[9-10]等展開了探索.隨著實驗測量技術(shù)和數(shù)值模擬的進步,人們逐漸闡明了旋轉(zhuǎn)爆震燃燒室內(nèi)詳細的流場結(jié)構(gòu)[11]與出口流場特征[12].

    前述大量研究主要聚焦于燃燒室內(nèi)爆震波的傳播過程,對推進性能關(guān)注較少,燃燒室出口一般不安裝尾噴管.然而,尾噴管作為噴氣發(fā)動機的關(guān)鍵部件,使氣流膨脹加速產(chǎn)生推力,對提高RDE 的推進性能同樣必不可少.Ishihara 等[13]研究表明,與不安裝噴管相比,安裝噴管后旋轉(zhuǎn)爆震燃燒室獲得了6%~10%的推力增益.不同結(jié)構(gòu)類型的尾噴管對推進性能的增益并不相同.高劍等[14]開展了不同類型尾噴管(收斂尾噴管、擴張尾噴管和收斂擴張尾噴管)以及不安裝尾噴管下的旋轉(zhuǎn)爆震推進性能實驗.結(jié)果表明,燃燒室穩(wěn)定工作范圍內(nèi),收斂尾噴管的增推效果最明顯,與不安裝噴管時相比,推力增益為34.4%;然而,安裝擴張尾噴管反而會降低推進性能.Goto 等[15]進行了RDE 地面軌道滑跑實驗,整個推進裝置獲得的平均加速度為2.77 m/s2,混合物比沖為144 s.

    安裝噴管會顯著提高推進性能,但尾噴管的收斂段和喉部將導(dǎo)致波系反射等,影響上游燃燒場和流動狀態(tài),嚴(yán)重時會對上游燃燒室內(nèi)爆震波的傳播狀態(tài)產(chǎn)生一定影響.作為與環(huán)形爆震燃燒室匹配的噴管構(gòu)型,塞式噴管自然進入了人們的視野.Jourdaine 等[12]采用三維數(shù)值模擬,研究了帶有塞式噴管的旋轉(zhuǎn)爆震燃燒室詳細內(nèi)流場,結(jié)果表明,噴管喉部反射激波向上游傳播與斜激波相互作用形成碰撞點,改變了原有流場結(jié)構(gòu),并可能會增加壁面的熱負荷.Ross 等[16]研究結(jié)果表明,安裝塞式噴管后,燃燒室內(nèi)存在與爆震波旋轉(zhuǎn)方向相反的壓力波,強度雖不足以產(chǎn)生新的爆震波,但可能影響爆震波的穩(wěn)定傳播.Frolov 等[17]以氫氣/空氣為反應(yīng)物,針對大尺寸RDE 安裝塞式噴管下的工作狀態(tài)和推進性能開展了實驗研究,發(fā)現(xiàn)出口阻塞比會增加燃燒室內(nèi)爆震波數(shù)量;另外,還觀察到一種爆震波縱向傳播的特殊模態(tài).針對無內(nèi)柱燃燒室,一般采用Laval 噴管與燃燒室匹配.Peng 等[18]在安裝不同喉道面積噴管的無內(nèi)柱燃燒室中,觀察到3種不同的傳播模態(tài),分別是單波、雙波和鋸齒波,噴管收縮比較小時的貧油邊界會出現(xiàn)鋸齒波模態(tài),并且隨著噴管收縮比增加,燃燒室穩(wěn)定工作的當(dāng)量比范圍減小.Zhang 等[19]在安裝Laval 噴管的燃燒室中,除觀察到單波和雙波模態(tài),還得到一種新的爆震波波形,即具有兩個主峰的單波模態(tài)(two dominant peak one-wave mode).他們認為,該波形由Laval 噴管收縮段產(chǎn)生的反射激波向上游傳播形成.

    綜上所述,尾噴管的存在,可以提高推進性能,但對上游燃燒室內(nèi)的爆震波傳播特性有較大影響,甚至誘發(fā)不穩(wěn)定的工作模態(tài).然而,已有研究大多聚焦于尾噴管對燃燒室內(nèi)爆震波的影響,針對安裝噴管后不同傳播模態(tài)下的起爆特性和傳播模態(tài)對推進性能影響的研究相對較少.針對上述問題,基于安裝和未安裝塞式噴管的環(huán)形旋轉(zhuǎn)爆震燃燒室開展了實驗研究,燃料和氧化劑分別采用乙烯和富氧空氣,在不同流量和當(dāng)量比下,分析對比尾噴管對上游爆震波傳播模態(tài)的影響,闡明了不同傳播模態(tài)下爆震波的傳播特性、起爆特性和推進性能的變化規(guī)律,可為相關(guān)設(shè)計提供支撐.

    1 實驗裝置

    實驗裝置示意圖見圖1,包括旋轉(zhuǎn)爆震實驗件、供給系統(tǒng)、測量與采集系統(tǒng)和控制系統(tǒng).旋轉(zhuǎn)爆震實驗件由旋轉(zhuǎn)爆震燃燒室和塞式噴管組成,供給系統(tǒng)包括燃料、氧化劑和吹除氮氣供給單元,采用的燃料和氧化劑分別為乙烯和氧氣體積分?jǐn)?shù)為50%的富氧空氣.測量與采集系統(tǒng)負責(zé)壓力、流量和推力數(shù)據(jù)的獲取與記錄,實驗中所有信號指令均通過控制系統(tǒng)觸發(fā).

    圖1 實驗裝置示意Fig.1 Schematic of the test bench

    1.1 旋轉(zhuǎn)爆震燃燒室和塞式噴管

    旋轉(zhuǎn)爆震燃燒室和塞式噴管示意圖見圖2.旋轉(zhuǎn)爆震燃燒室為環(huán)形結(jié)構(gòu),主要由燃燒室外環(huán)和內(nèi)柱構(gòu)成,內(nèi)流道外壁直徑為72 mm,長度為75 mm,內(nèi)柱直徑為34 mm.為與塞式噴管中心體平順連接,內(nèi)柱直徑在燃燒室出口位置漸擴為42 mm,噴管入口面積為2 686 mm2.采用環(huán)縫-噴孔的噴注形式,氧化劑通過喉道寬度為0.4 mm 的環(huán)縫進入燃燒室,燃料經(jīng)過100 個周向均布的噴孔噴注,噴孔直徑為0.6 mm.塞式噴管為完全外膨脹式設(shè)計,由噴管外環(huán)和中心體構(gòu)成,喉道面積為845 mm2,長度為91.2 mm,具體幾何參數(shù)見表1.

    圖2 旋轉(zhuǎn)爆震燃燒室和塞式噴管示意Fig.2 Schematic of the rotating detonation chamber with an aerospike nozzle

    表1 塞式噴管的幾何參數(shù)Tab.1 Geometric parameters of the aerospike nozzle

    1.2 實驗數(shù)據(jù)采集與時序控制

    為監(jiān)測供給流量,在燃料和氧化劑供給管路上各安裝1 臺質(zhì)量流量計(測量精度為±1%).燃燒室和預(yù)爆管上各傳感器的分布位置如圖3 所示.燃燒室頭部布置4 個高頻動態(tài)壓力傳感器p1~p4(SINOCERA CY-YD-205 H,測量精度為±3%),其中p1~p3 周向均布,任意兩個傳感器之間呈120°,p1和 p 3 位于同一軸向位置,距離噴注環(huán)縫出口13 mm,p2 與預(yù)爆管出口位于同一軸向位置,距離噴注環(huán)縫出口20 mm.p4 與p1 位于同一軸線位置,與p1 距離20 mm,具體空間位置分布見圖3.為監(jiān)測預(yù)爆管出口是否已形成充分發(fā)展的爆震波,在預(yù)爆管出口上游布置了2 個動態(tài)壓力傳感器ppre1和ppre2,間距為50 mm.在噴管喉道上游軸向距離15 mm 處,布置壓阻式壓力傳感器pN(KELLER PA-21Y,測量精度為±1%),以測量燃燒室對應(yīng)位置處的靜壓;當(dāng)燃燒室出口未安裝噴管時,該傳感器安裝于p3 的正后方,與p4 位于同一排.此外,分別在燃料和氧化劑的供給腔布置壓阻式壓力傳感器以監(jiān)測噴前壓力.旋轉(zhuǎn)爆震燃燒室實驗件整體置于動架上,動架滑塊與靜架滑軌相接,且整體可在滑軌上自由滑動,實驗件和靜架推力壁之間布置推力傳感器(KISTLER 9331B),用以測量產(chǎn)生的推力.傳感器測得數(shù)據(jù)均由最高采樣頻率為1 MHz 的十六通道采集儀記錄.

    圖3 壓力傳感器分布Fig.3 Distribution of the pressure transducers

    單次實驗的時序控制如圖4 所示.首先,觸發(fā)數(shù)據(jù)采集儀開始記錄,并向燃燒室供給燃料和氧化劑;隨后,向預(yù)爆管內(nèi)填充燃料和氧化劑,填充完成后切斷供給;對燃燒室供給燃料和氧化劑約330 ms 后,初步建立穩(wěn)定的冷態(tài)流場,并觸發(fā)點火能量約為50 mJ的火花塞點燃預(yù)爆管內(nèi)的可燃混氣;燃燒室內(nèi)燃燒過程持續(xù)300 ms 后,關(guān)閉燃燒室供給燃料和氧化劑的通路,并打開氮氣管路電磁閥,吹除燃燒室內(nèi)殘余的未燃氣體與燃燒產(chǎn)物;最后,切斷氮氣供給并停止數(shù)據(jù)采集,完成一次實驗過程.

    圖4 實驗過程時序圖Fig.4 Time sequence of the control signals for the operations

    2 實驗結(jié)果與討論

    2.1 傳播模態(tài)分析

    不論燃燒室出口是否安裝噴管,均開展了流量約為200 g/s 和120 g/s 兩種工況下的實驗,對應(yīng)當(dāng)量比為0.3~1.6.安裝噴管時,出現(xiàn)了5 種燃燒波的傳播模態(tài),分別為穩(wěn)定單波模態(tài)、雙波對撞模態(tài)、轉(zhuǎn)變模態(tài)、周期性振蕩模態(tài)和緩燃模態(tài);未安裝噴管時,僅觀察到穩(wěn)定單波模態(tài)、雙波對撞模態(tài)和緩燃模態(tài),如圖5 所示.其中,穩(wěn)定單波模態(tài)是指燃燒室內(nèi)存在一個爆震波沿燃燒室周向順時針或逆時針旋轉(zhuǎn).雙波對撞模態(tài)中,燃燒室內(nèi)存在兩個相向運動的爆震波,發(fā)生周期性碰撞,碰撞點附近的局部高壓會阻礙燃料和氧化劑進入燃燒室,導(dǎo)致爆震波解耦,后因供氣過程恢復(fù)重新生成兩個反向傳播的爆震波.若二者傳播速度相近,碰撞點位置幾乎固定不變;如果傳播速度略有差異,則碰撞點位置不斷變化.轉(zhuǎn)變模態(tài)是指,前期爆震波的峰值壓力呈現(xiàn)周期性振蕩,后期轉(zhuǎn)變?yōu)榉€(wěn)定單波模態(tài).周期性振蕩模態(tài)是指,在整個工作過程中,均存在明顯的峰值壓力振蕩.緩燃模態(tài)對應(yīng)于無爆震產(chǎn)生的連續(xù)緩燃情形,高頻動態(tài)壓力傳感器未監(jiān)測到明顯的壓力脈動,該模態(tài)下可在噴管外(未安裝噴管時,則是燃燒室外)觀察到較長的尾焰.

    圖5 不同流量下燃燒波模態(tài)隨當(dāng)量比的變化Fig.5 Operating equivalence ratio ranges of the propagation modes at different mass flow rates

    當(dāng)燃燒室出口未安裝噴管、流量為200 g/s 時,僅觀察到穩(wěn)定單波模態(tài);總流量降低至120 g/s 時,傳播模態(tài)以穩(wěn)定單波模態(tài)為主,僅在當(dāng)量比較高時出現(xiàn)雙波對撞和緩燃模態(tài),且高當(dāng)量比下的緩燃邊界寬于有噴管情形.然而,在燃燒室尾部安裝噴管后,不同供給流量、相同當(dāng)量比范圍內(nèi),依次出現(xiàn)了上述5 種不同的傳播模態(tài),轉(zhuǎn)變模態(tài)和周期性振蕩模態(tài)下均會出現(xiàn)爆震波峰值壓力振蕩的現(xiàn)象.可以看出,兩種模態(tài)均在當(dāng)量比較高時出現(xiàn),轉(zhuǎn)變模態(tài)可認為是穩(wěn)定單波模態(tài)向周期性振蕩模態(tài)轉(zhuǎn)變的過渡.不同流量下,出現(xiàn)轉(zhuǎn)變模態(tài)的臨界當(dāng)量比范圍相近(均為1.0 左右),出現(xiàn)緩燃模態(tài)的當(dāng)量比邊界也相近(均為1.4 左右),但出現(xiàn)周期性振蕩模態(tài)的當(dāng)量比下邊界有所降低.

    穩(wěn)定單波模態(tài)和雙波對撞模態(tài)在旋轉(zhuǎn)爆震工作過程中較為常見,已有大量文獻針對二者開展了較為充分的分析[6,20-21],本文對其傳播特性不再展開討論.由圖5 可知,未安裝噴管時,傳播模態(tài)主要為穩(wěn)定單波模態(tài),而轉(zhuǎn)變模態(tài)和周期性振蕩模態(tài)僅在安裝噴管后出現(xiàn),針對二者的分析較少.因此,下面對其進行分析討論.

    2.1.1 轉(zhuǎn)變模態(tài)

    當(dāng)供給流量為192.1 g/s,當(dāng)量比為1.0 時,壓力傳感器測得的p1、p2和pN在60 ms 內(nèi)的壓力波形如圖6 所示.可以看出,約36 ms 前,測得p1和p2的峰值壓力有明顯振蕩,之后峰值壓力趨于穩(wěn)定.因傳感器軸向位置不同,測得的爆震波峰值壓力也略有差別,p1平均峰值壓力約為0.65 MPa,p2平均峰值壓力約為0.89 MPa.此外,測得pN的壓力波形同樣表明,模態(tài)轉(zhuǎn)變之前,燃燒室內(nèi)靜壓存在小幅波動,平均室壓為0.223 MPa;模態(tài)轉(zhuǎn)變之后,燃燒室內(nèi)靜壓趨于穩(wěn)定,平均室壓變?yōu)?.216 MPa.圖7 為模態(tài)轉(zhuǎn)變過程34~38 ms 間的壓力波形.可以看出,p1和p2均在35 ms 左右發(fā)生峰值壓力振蕩,后轉(zhuǎn)變?yōu)榉€(wěn)定的傳播模態(tài),并且爆震波后存在一個反射波波峰,以相同的傳播速度跟隨爆震波同向旋轉(zhuǎn).

    圖6 轉(zhuǎn)變模態(tài)下測得的壓力波形Fig.6 Pressure profiles in the self-adjestive mode

    圖7 模態(tài)轉(zhuǎn)變過程中的壓力波形Fig.7 Enlargement pressure profiles during the transition process

    對圖6 中p1的壓力波形進行頻譜分析,結(jié)果如圖8 所示.可以看出,模態(tài)轉(zhuǎn)變之后的穩(wěn)定傳播過程,僅存在一個主頻8 362 Hz,對應(yīng)的平均傳播速度可由式(1)得到.

    圖8 轉(zhuǎn)變模態(tài)下測得的p1 壓力波形頻譜分析Fig.8 Frequency spectrum of the self-adjustive mode

    式中:rout為燃燒室內(nèi)流道外壁半徑;f 為頻譜分析中的主頻率;n 為波數(shù).由傳感器的周向位置分布和壓力波形可判斷,燃燒室內(nèi)的壓力波均同向傳播.假設(shè)燃燒室內(nèi)存在單波,則可得平均波速為1 891.4 m/s,為理論CJ 速度的88.5%;若燃燒室內(nèi)存在雙波,則單個壓力波的平均傳播速度為945.7 m/s,低于爆震波波后聲速(1 160.4 m/s),即壓力波為弱壓縮波,這與測得的壓力波形上升沿僅約為5 μs 矛盾,故為穩(wěn)定的單波傳播模態(tài).受側(cè)向膨脹[22]、摻混效果[23]和邊界層[24]的影響,其傳播速度通常有一定的速度虧損.前36 ms 峰值壓力存在振蕩時,除主頻8 362 Hz 外,存在7 141 Hz 的次頻,以及二階頻率16 724 Hz(為主頻的2 倍),對應(yīng)的平均傳播速度為1 615.3 m/s,為理論CJ 速度的75.6%,該現(xiàn)象亦出現(xiàn)在周期性振蕩模態(tài).此外,在轉(zhuǎn)變?yōu)榉€(wěn)定單波模態(tài)的前36 ms 內(nèi),峰值壓力存在明顯的振蕩現(xiàn)象,振蕩頻率為946 Hz.

    2.1.2 周期性振蕩模態(tài)

    當(dāng)供給流量為120.5 g/s、當(dāng)量比為1.28 時,壓力傳感器測得的p1、p2和pN在60 ms 內(nèi)壓力波形如圖9 所示.可以看出,p1和p2峰值壓力均存在周期性的振蕩,波動范圍在0.25~1.50 MPa,pN亦存在小幅波動,燃燒室的平均壓力為0.160 MPa,與2.1.1 小節(jié)中模態(tài)未轉(zhuǎn)變前相似.相應(yīng)地在20~30 ms 內(nèi),高頻動態(tài)傳感器的局部原始壓力波形如圖10(a)所示,可以看出,峰值壓力的振蕩具有明顯的周期性.為便于觀察兩個壓力傳感器間的信號差異,將圖10(a)紅色虛線框內(nèi)的壓力波形繪制于同一坐標(biāo)軸下,并分別提取p1和p2壓力波形的峰值壓力包絡(luò)線,如圖10(b)所示.可以看出,兩個傳感器測得的壓力峰值波動類似于準(zhǔn)正弦曲線,雖然壓力傳感器p1 和p2 空間分布位置呈120°,但兩者峰值壓力振蕩相位幾乎相同.

    圖9 周期性振蕩模態(tài)下測得的壓力波形Fig.9 Pressure profiles in the periodic oscillation mode

    圖10 周期性振蕩模態(tài)下的局部壓力波形Fig.10 Enlargement pressure profiles in the the periodic oscillation mode

    對60 ms 內(nèi)p1壓力波形做頻譜分析,結(jié)果如圖11 所示.可以看出,傳播過程中均存在兩個傳播頻率,即主頻8 301 Hz(對應(yīng)平均波速為1 877.6 m/s,為理論CJ 速度的88.3%)和次頻7 004 Hz(對應(yīng)的平均波速為1 584.3 m/s,為理論CJ 速度的74.5%).圖11左側(cè)時頻圖表明,主頻8 301 Hz 和次頻7 004 Hz 并不連續(xù).原因在于,爆震波壓力峰值存在同相振蕩,在振蕩波峰附近,爆震波的傳播速度較高,頻率約為8 301 Hz;在振蕩波谷附近,爆震波強度減弱,傳播速度降低,頻率降至約7 004 Hz.此外,峰值壓力周期性振蕩的頻率為1 304 Hz,均為1 000 Hz 左右的中高頻振蕩.

    圖11 周期性振蕩模態(tài)下測得的p1 壓力波形頻譜分析Fig.11 Frequency spectrum of the periodic oscillation mode

    2.2 不同模態(tài)下的起爆特性

    由前文分析可知,不同傳播模態(tài)下,爆震波的傳播特性不同.未安裝噴管條件下,供給流量為200.8 g/s、當(dāng)量比為1.08 時,預(yù)爆管上測得的ppre1和ppre2峰值壓力之間的時間差Δtpre為0.019 ms,由兩個傳感器之間的距離可得壓力波經(jīng)過的平均傳播速度為2 631 m/s,與估算的理論CJ 速度相近,可認為預(yù)爆管中已得到充分發(fā)展的爆震波.爆震波由預(yù)爆管傳出后,進入物理空間更大的旋轉(zhuǎn)爆震燃燒室,傳播通道突擴,可能發(fā)生解耦,轉(zhuǎn)變?yōu)榫徣疾ê颓皩?dǎo)激波向四周傳播,若未發(fā)生解耦,爆震波的強度亦會有所減弱.穩(wěn)定單波模態(tài)的起始過程中,測得的ppre2和p3壓力波形如圖12 所示.可以看出,預(yù)爆管內(nèi)爆震波傳入燃燒室后,測得p3有一個很高的壓力峰值(如圖12 約5 ms 處,約4.4 MPa).原因在于,點火之前旋轉(zhuǎn)爆震燃燒室內(nèi)已經(jīng)填充滿燃料和氧化劑,預(yù)爆管出口射流點燃可燃混氣并形成爆震波,爆震波的高度接近整個燃燒室長度,壓力傳感器p3 靠近燃燒室頭部,受側(cè)向膨脹的影響很小,故測得的壓力峰值明顯高于建立穩(wěn)定旋轉(zhuǎn)爆震后的情形.隨后,燃料和氧化劑的填充過程受阻,爆震波無法維持,解耦成為緩燃波和前導(dǎo)激波.待燃燒室壓力下降,填充過程恢復(fù),局部熱點會觸發(fā)穩(wěn)定爆震過程,如圖12 中紅色虛線框內(nèi)放大圖.將預(yù)爆管出口壓力傳感器ppre2處得到峰值壓力時刻與燃燒室內(nèi)再次起爆時刻的時間間隔Δt 定義為起爆時間,燃燒室內(nèi)最終形成穩(wěn)定爆震的時間Δts定義為建立穩(wěn)定爆震的時間,可見起爆時間與建立穩(wěn)定爆震的時間相等,為2.19 ms.

    圖12 未安裝噴管下穩(wěn)定單波模態(tài)起始過程的壓力波形Fig.12 Pressure profiles of the initiation process in the single-wave mode without nozzle

    安裝噴管條件下,供給流量為194.1 g/s、當(dāng)量比為0.98 時,穩(wěn)定單波模態(tài)的起始過程如圖13 所示.與未安裝噴管工況的起爆過程相同,預(yù)爆管出口射流點燃燃燒室內(nèi)的可燃混氣并形成爆震波,產(chǎn)生一道強度較高的爆震波并影響填充過程,爆震波解耦為緩燃波.然而,之后緩燃波轉(zhuǎn)變?yōu)楸鸬倪^程中,噴管的節(jié)流作用影響了燃料和氧化劑的噴注過程,在爆震形成的初期,爆震波的峰值壓力存在幾個周期的振蕩之后立即轉(zhuǎn)變?yōu)榉€(wěn)定旋轉(zhuǎn)的單波模態(tài)(見圖13 內(nèi)紅色虛線框中的放大部分),該過程與轉(zhuǎn)變模態(tài)中模態(tài)的轉(zhuǎn)變過程類似.因此,建立穩(wěn)定傳播爆震波的時間 Δts大于未安裝噴管的條件,圖中起爆時間為5.53 ms,建立穩(wěn)定爆震的時間為7.59 ms.

    圖13 安裝噴管下穩(wěn)定單波模態(tài)起始過程的壓力波形Fig.13 Pressure profiles of the initiation process in the single-wave mode with nozzle

    與周期性振蕩模態(tài)起始過程相同,轉(zhuǎn)變模態(tài)亦存在峰值壓力的周期性振蕩.因此,以周期性振蕩的起爆過程為例,分析其起爆特性.安裝噴管條件下,供給流量為205.6 g/s、當(dāng)量比為1.15 時,周期性振蕩模態(tài)起始過程如圖14 所示.可以看出,起爆過程與未安裝噴管時的穩(wěn)定單波模態(tài)不同,緩燃轉(zhuǎn)變?yōu)楸鸬某跗冢鸩ǚ逯祲毫创嬖谝欢ㄕ袷幏?12~19 ms 時),約22 ms 之后峰值壓力出現(xiàn)了明顯且穩(wěn)定的周期性振蕩,表明峰值壓力的振蕩在起爆過程已形成.由圖14 可知,起爆時間為3.39 ms,與未安裝噴管工況的起爆時間相近.然而,由于傳播過程中未形成穩(wěn)定的旋轉(zhuǎn)爆震,此模態(tài)下建立穩(wěn)定爆震的時間不作統(tǒng)計.

    圖14 安裝噴管下周期性振蕩模態(tài)起始過程的壓力波形Fig.14 Pressure profiles of the initiation process in the periodic oscillation mode with nozzle

    2.3 不同模態(tài)下的推進特性

    穩(wěn)定單波模態(tài)下,流量約為200 g/s 時,安裝(流量為191.6 g/s,當(dāng)量比為1.00)與未安裝噴管(流量為200.8 g/s,當(dāng)量比為1.08)時的燃燒室室壓與推力曲線如圖15 所示.可以看出,在點火之前,燃料和氧化劑的供給會產(chǎn)生一定的冷態(tài)推力,安裝噴管時冷態(tài)推力(約40 N)略大于未安裝噴管工況(約34 N).點火之后,燃燒室內(nèi)室壓和推力信號均有明顯的過沖現(xiàn)象.原因在于,與穩(wěn)定工作時爆震波只在燃燒室頭部旋轉(zhuǎn)的情形相比,預(yù)爆管初始點燃了充滿整個燃燒室的可燃氣體,爆震波高度大、強度高,增壓效果明顯,產(chǎn)生的高溫已燃氣體膨脹并向出口排出,使初始室內(nèi)靜壓和推力出現(xiàn)極大值.未安裝噴管時室壓無明顯上升,平均室壓為0.094 MPa;然而,安裝噴管后,室壓較點火前明顯上升,平均室壓為0.215 MPa.取燃燒過程結(jié)束前約100 ms 內(nèi)的平均推力(見圖15(a))作為穩(wěn)定工作下的推力,未安裝噴管時為99.4 N,安裝噴管后為200.7 N,推力增益為101.9%.

    圖15 穩(wěn)定單波模態(tài)下測得的燃燒室壓力和推力曲線Fig.15 Static pressure and thrust in the single-wave mode

    由前文可知,轉(zhuǎn)變模態(tài)發(fā)生轉(zhuǎn)變過程前和周期性振蕩模態(tài)中,爆震波的峰值壓力均存在明顯的振蕩過程.轉(zhuǎn)變模態(tài)中,發(fā)生轉(zhuǎn)變后傳播模態(tài)變?yōu)榉€(wěn)定的單波,針對轉(zhuǎn)變點前后的變化可對比分析周期性振蕩和穩(wěn)定傳播過程之間燃燒室室壓和推力的差異.轉(zhuǎn)變模態(tài)下(流量為192.1 g/s,當(dāng)量比為1.02),燃燒室室壓和推力曲線如圖16 所示.可以看出,模態(tài)轉(zhuǎn)變前,pN的幅值變化范圍較大,使噴管入口氣流的非穩(wěn)態(tài)特性增強;轉(zhuǎn)變之后,盡管pN變得更穩(wěn)定,但平均室壓有所下降,由0.223 MPa 變?yōu)?.216 MPa,下降了3.1%.推力曲線上模態(tài)轉(zhuǎn)變過程不明顯,但轉(zhuǎn)變前后推力仍有一定差異,分別取約100 ms 內(nèi)的平均推力,模態(tài)轉(zhuǎn)變前為216.3 N,模態(tài)轉(zhuǎn)變后為199.8 N,下降了7.6%.周期性振蕩模態(tài)下(供給流量為205.6 g/s,當(dāng)量比為1.14),燃燒室室壓和推力曲線如圖17 所示.與轉(zhuǎn)變模態(tài)轉(zhuǎn)變前相似,室壓在整個工作過程中均存在一定的波動,平均室壓為0.237 MPa,對應(yīng)的平均推力約為229.7 N.需要注意,與圖15(b)相比,圖16 和圖17 推力曲線明顯變粗,表明爆震波峰值壓力出現(xiàn)振蕩時,對應(yīng)的推力曲線振蕩加劇.盡管轉(zhuǎn)變模態(tài)下爆震波最終以穩(wěn)定單波模態(tài)傳播,燃燒室壓力的波動不明顯,但推力曲線的振蕩幅值并未減弱,說明峰值壓力周期性振蕩引起的推力振蕩無法快速消除.

    圖16 轉(zhuǎn)變模態(tài)下的燃燒室室壓和推力曲線Fig.16 Static pressure and thrust in the self-adjestive mode

    圖17 周期性振蕩模態(tài)下的燃燒室室壓和推力曲線Fig.17 Static pressure and thrust in the periodic oscillation mode

    盡管周期性振蕩過程中室壓和推力均略大于穩(wěn)定單波過程,但亦會造成噴管入口氣流的非穩(wěn)態(tài)特性增大,導(dǎo)致推力下降.此外,周期性振蕩模態(tài)下,爆震波峰值壓力和燃燒室室壓存在周期性振蕩,使推力振蕩加劇,在實際工作中,應(yīng)避免產(chǎn)生此模態(tài).

    3 結(jié)論

    針對燃燒室出口安裝塞式噴管與否,開展了200 g/s 和120 g/s 兩種流量下的旋轉(zhuǎn)爆震實驗研究,得到了不同當(dāng)量比下的旋轉(zhuǎn)爆震波傳播模態(tài)、起爆特性和推進性能,主要結(jié)論如下:

    (1) 兩種流量下,燃燒室出口未安裝噴管時,傳播模態(tài)以穩(wěn)定單波模態(tài)為主;安裝噴管時,隨著當(dāng)量比增加,出現(xiàn)了5 種不同的傳播模態(tài),分別為穩(wěn)定單波模態(tài)、雙波對撞模態(tài)、轉(zhuǎn)變模態(tài)、周期性振蕩模態(tài)和緩燃模態(tài).

    (2) 轉(zhuǎn)變模態(tài)和周期性振蕩模態(tài)下,均觀察到峰值壓力的周期振蕩現(xiàn)象,周向均布的3 個傳感器測得的峰值壓力振蕩幾乎同相,振蕩頻率約為1 000 Hz;轉(zhuǎn)變模態(tài)下,在轉(zhuǎn)變點之后發(fā)展為穩(wěn)定單波模態(tài).

    (3) 與未安裝噴管的穩(wěn)定單波模態(tài)起爆過程相比,安裝噴管后,由于噴管的節(jié)流作用,穩(wěn)定單波模態(tài)的起爆過程存在峰值壓力振蕩向無振蕩的轉(zhuǎn)變,與轉(zhuǎn)變模態(tài)中的模態(tài)轉(zhuǎn)變類似,穩(wěn)定爆震的建立時間明顯增加;周期性振蕩模態(tài)下,起爆時間與未安裝噴管工況相近,但峰值壓力存在周期性振蕩.

    (4) 在相同供給流量和當(dāng)量比條件下,安裝噴管后的推力明顯提高,典型工況下比未安裝噴管時增加了101.9%;轉(zhuǎn)變模態(tài)下,轉(zhuǎn)變點前的室壓和推力均略大于轉(zhuǎn)變后穩(wěn)定單波模態(tài),推力下降了7.6%;周期性振蕩模態(tài)下,燃燒室工作過程中室壓始終存在小幅波動,增加了噴管入口氣流的非穩(wěn)態(tài)特性,造成推力損失;并且,峰值壓力的周期性振蕩會造成推力輸出的不穩(wěn)定,應(yīng)避免產(chǎn)生此模態(tài).

    猜你喜歡
    爆震周期性燃燒室
    燃燒室形狀對國六柴油機性能的影響
    雷克薩斯車系爆震控制基理介紹
    一種熱電偶在燃燒室出口溫度場的測量應(yīng)用
    電子制作(2019年19期)2019-11-23 08:41:54
    數(shù)列中的周期性和模周期性
    肺爆震傷治療的研究進展
    一類整數(shù)遞推數(shù)列的周期性
    基于擴頻碼周期性的單通道直擴通信半盲分離抗干擾算法
    長距離爆震式點火槍設(shè)計
    焊接(2015年9期)2015-07-18 11:03:52
    高幾何壓縮比活塞的燃燒室形狀探討
    3.0T磁敏感加權(quán)成像在兔顱腦爆震傷早期出血灶檢測及預(yù)后評估中的應(yīng)用
    国产激情偷乱视频一区二区| 亚洲成av人片免费观看| 88av欧美| 我的老师免费观看完整版| 久久久久精品国产欧美久久久| 亚洲欧美清纯卡通| 九九热线精品视视频播放| 国产av在哪里看| 91狼人影院| 国内精品久久久久精免费| 男人的好看免费观看在线视频| 床上黄色一级片| 三级国产精品欧美在线观看| 色综合婷婷激情| 日本三级黄在线观看| 岛国在线免费视频观看| 亚洲,欧美精品.| 亚洲乱码一区二区免费版| 免费在线观看影片大全网站| 99久久99久久久精品蜜桃| 久久精品国产99精品国产亚洲性色| 一个人免费在线观看电影| 国产三级黄色录像| 久久久久久国产a免费观看| 欧美性感艳星| 我的女老师完整版在线观看| av在线蜜桃| 黄色一级大片看看| 九色国产91popny在线| 无人区码免费观看不卡| 最近中文字幕高清免费大全6 | 真人一进一出gif抽搐免费| 国产一区二区三区在线臀色熟女| 99久久精品热视频| 村上凉子中文字幕在线| 成年女人看的毛片在线观看| 黄色日韩在线| 国产成+人综合+亚洲专区| 亚洲美女搞黄在线观看 | 国产免费男女视频| 男女那种视频在线观看| 精品熟女少妇八av免费久了| 色哟哟·www| 男女做爰动态图高潮gif福利片| 三级男女做爰猛烈吃奶摸视频| 成人精品一区二区免费| 成熟少妇高潮喷水视频| 麻豆成人午夜福利视频| 一夜夜www| 变态另类丝袜制服| 亚洲美女黄片视频| 亚洲中文字幕日韩| 亚洲,欧美,日韩| 欧美高清性xxxxhd video| aaaaa片日本免费| 真人一进一出gif抽搐免费| 老鸭窝网址在线观看| 国产伦人伦偷精品视频| 久久久久久大精品| www.色视频.com| 日韩欧美精品v在线| 亚洲专区国产一区二区| 国产精品人妻久久久久久| 99久久成人亚洲精品观看| 久久久国产成人精品二区| 成人亚洲精品av一区二区| 欧美最黄视频在线播放免费| 久久精品国产亚洲av涩爱 | 欧美黄色淫秽网站| 俺也久久电影网| 日韩av在线大香蕉| 欧美一区二区国产精品久久精品| 天堂动漫精品| 亚洲片人在线观看| 97碰自拍视频| 在线看三级毛片| 色精品久久人妻99蜜桃| 高清毛片免费观看视频网站| 我要搜黄色片| 久久精品久久久久久噜噜老黄 | 亚洲五月婷婷丁香| 亚洲性夜色夜夜综合| 露出奶头的视频| 亚洲av日韩精品久久久久久密| 亚洲自偷自拍三级| 久久这里只有精品中国| 久久国产乱子伦精品免费另类| 少妇人妻一区二区三区视频| 亚洲专区国产一区二区| 白带黄色成豆腐渣| 日本熟妇午夜| 国产精品亚洲一级av第二区| 国产精品人妻久久久久久| 真实男女啪啪啪动态图| 精品日产1卡2卡| 美女 人体艺术 gogo| 热99re8久久精品国产| 精品久久久久久久末码| 亚洲最大成人av| 色噜噜av男人的天堂激情| 又粗又爽又猛毛片免费看| 久久精品91蜜桃| 亚洲国产高清在线一区二区三| 亚洲一区二区三区不卡视频| 欧美高清成人免费视频www| 在线播放国产精品三级| 一卡2卡三卡四卡精品乱码亚洲| 99久久精品国产亚洲精品| 久久性视频一级片| 性插视频无遮挡在线免费观看| 亚洲狠狠婷婷综合久久图片| 欧美性猛交╳xxx乱大交人| 亚洲狠狠婷婷综合久久图片| 亚洲国产欧美人成| 美女高潮喷水抽搐中文字幕| 亚洲激情在线av| 一区福利在线观看| 2021天堂中文幕一二区在线观| 国产极品精品免费视频能看的| 一区二区三区激情视频| 亚洲 欧美 日韩 在线 免费| 国产精品1区2区在线观看.| 人妻制服诱惑在线中文字幕| 99精品在免费线老司机午夜| 又黄又爽又刺激的免费视频.| 校园春色视频在线观看| 9191精品国产免费久久| 欧美一区二区国产精品久久精品| 婷婷色综合大香蕉| 亚洲精品成人久久久久久| 亚洲av成人不卡在线观看播放网| 精品一区二区三区视频在线观看免费| 亚洲欧美日韩东京热| 欧美+日韩+精品| 美女被艹到高潮喷水动态| 99精品在免费线老司机午夜| 国产色婷婷99| 亚洲av电影不卡..在线观看| 亚洲男人的天堂狠狠| 久久久精品大字幕| 日韩大尺度精品在线看网址| 观看免费一级毛片| 校园春色视频在线观看| 老司机午夜福利在线观看视频| 免费看日本二区| 中出人妻视频一区二区| 99久久成人亚洲精品观看| 看黄色毛片网站| 国语自产精品视频在线第100页| 亚洲电影在线观看av| 中文亚洲av片在线观看爽| 在线十欧美十亚洲十日本专区| 国产精品乱码一区二三区的特点| 老熟妇乱子伦视频在线观看| 色5月婷婷丁香| 亚洲中文日韩欧美视频| 午夜精品一区二区三区免费看| 蜜桃久久精品国产亚洲av| 在线播放国产精品三级| 直男gayav资源| 亚洲男人的天堂狠狠| 亚洲美女搞黄在线观看 | 久久久久久久久久成人| 9191精品国产免费久久| 美女被艹到高潮喷水动态| 免费观看精品视频网站| 日日摸夜夜添夜夜添小说| 成年女人永久免费观看视频| 久久精品久久久久久噜噜老黄 | 亚洲av一区综合| 国产一区二区三区在线臀色熟女| 国产av不卡久久| 90打野战视频偷拍视频| 亚州av有码| 国产成人福利小说| 长腿黑丝高跟| 亚洲无线在线观看| 听说在线观看完整版免费高清| 床上黄色一级片| 欧美一区二区国产精品久久精品| 日本a在线网址| 毛片女人毛片| 中文字幕免费在线视频6| 亚洲三级黄色毛片| 在线十欧美十亚洲十日本专区| 村上凉子中文字幕在线| 噜噜噜噜噜久久久久久91| 嫩草影院入口| 日韩av在线大香蕉| 性色avwww在线观看| 国产麻豆成人av免费视频| 一级黄片播放器| 夜夜夜夜夜久久久久| 亚洲人成网站在线播| 99久久精品一区二区三区| 激情在线观看视频在线高清| 欧美精品啪啪一区二区三区| 欧美+日韩+精品| 看十八女毛片水多多多| 三级男女做爰猛烈吃奶摸视频| 国内精品美女久久久久久| 久久久久久九九精品二区国产| 一区二区三区高清视频在线| 午夜福利18| 成人三级黄色视频| 在线观看一区二区三区| 五月玫瑰六月丁香| av在线天堂中文字幕| 给我免费播放毛片高清在线观看| 国产精品美女特级片免费视频播放器| 日本撒尿小便嘘嘘汇集6| 99久国产av精品| 欧美国产日韩亚洲一区| 亚洲成av人片在线播放无| 国语自产精品视频在线第100页| 精品人妻视频免费看| 亚洲av.av天堂| 综合色av麻豆| 男女之事视频高清在线观看| 国产精品精品国产色婷婷| 午夜福利在线观看吧| 成人性生交大片免费视频hd| 精品一区二区三区视频在线| 欧美黄色淫秽网站| 欧美潮喷喷水| 搞女人的毛片| 一个人看的www免费观看视频| 在线a可以看的网站| 免费在线观看影片大全网站| av国产免费在线观看| 午夜福利欧美成人| 午夜a级毛片| 免费看美女性在线毛片视频| 国产淫片久久久久久久久 | 午夜免费激情av| 亚洲五月天丁香| 嫩草影院新地址| 嫁个100分男人电影在线观看| 两人在一起打扑克的视频| 亚洲avbb在线观看| 国产精品一区二区免费欧美| 日日夜夜操网爽| 高清毛片免费观看视频网站| 欧美性猛交黑人性爽| 久久精品综合一区二区三区| 51午夜福利影视在线观看| 桃红色精品国产亚洲av| 久久99热6这里只有精品| 特大巨黑吊av在线直播| 精品人妻一区二区三区麻豆 | 波多野结衣巨乳人妻| 少妇高潮的动态图| 亚洲 国产 在线| 男人舔女人下体高潮全视频| 欧美性猛交黑人性爽| 桃红色精品国产亚洲av| 18禁黄网站禁片午夜丰满| 一卡2卡三卡四卡精品乱码亚洲| 久久国产精品人妻蜜桃| 国产成+人综合+亚洲专区| 亚洲精品在线美女| 久久精品国产自在天天线| 欧美日韩中文字幕国产精品一区二区三区| 老鸭窝网址在线观看| 日韩有码中文字幕| 日本在线视频免费播放| 成人精品一区二区免费| 黄色一级大片看看| 嫩草影视91久久| 宅男免费午夜| 最近中文字幕高清免费大全6 | 国产精品一及| 91久久精品电影网| 性插视频无遮挡在线免费观看| 国产成人影院久久av| 热99re8久久精品国产| 欧美绝顶高潮抽搐喷水| 黄色日韩在线| 国产av在哪里看| 亚洲精品粉嫩美女一区| 在线观看美女被高潮喷水网站 | 级片在线观看| 伊人久久精品亚洲午夜| 99久久久亚洲精品蜜臀av| 色综合亚洲欧美另类图片| 搡老妇女老女人老熟妇| 久久久国产成人免费| 日本三级黄在线观看| 亚洲欧美日韩东京热| 好男人电影高清在线观看| 亚洲无线在线观看| 成年人黄色毛片网站| 国产精品久久电影中文字幕| 亚洲无线观看免费| 国产极品精品免费视频能看的| 搡女人真爽免费视频火全软件 | av在线天堂中文字幕| 美女cb高潮喷水在线观看| 美女cb高潮喷水在线观看| 国产久久久一区二区三区| 欧美乱色亚洲激情| 国产三级在线视频| 亚洲av电影在线进入| 久久久久性生活片| www.熟女人妻精品国产| 欧美精品啪啪一区二区三区| 亚洲成人精品中文字幕电影| 黄色视频,在线免费观看| 99久久精品国产亚洲精品| 亚洲午夜理论影院| 日韩精品青青久久久久久| 亚州av有码| 男女之事视频高清在线观看| 夜夜夜夜夜久久久久| 午夜免费男女啪啪视频观看 | 国产精品嫩草影院av在线观看 | 51午夜福利影视在线观看| 国产亚洲av嫩草精品影院| 性欧美人与动物交配| 亚洲男人的天堂狠狠| 又黄又爽又免费观看的视频| 亚洲精品日韩av片在线观看| 在线观看美女被高潮喷水网站 | 人人妻人人看人人澡| 亚洲精品亚洲一区二区| 小说图片视频综合网站| 亚洲自拍偷在线| 丰满乱子伦码专区| 国内精品久久久久久久电影| 91狼人影院| 国产男靠女视频免费网站| 丰满乱子伦码专区| 精品无人区乱码1区二区| 精品国内亚洲2022精品成人| 精品一区二区三区人妻视频| 亚洲内射少妇av| 亚洲人成电影免费在线| 国产精品不卡视频一区二区 | 国产淫片久久久久久久久 | 久久国产精品影院| 色综合亚洲欧美另类图片| 性插视频无遮挡在线免费观看| 最新在线观看一区二区三区| 午夜免费男女啪啪视频观看 | 精品日产1卡2卡| 久久久久国产精品人妻aⅴ院| 日本五十路高清| 综合色av麻豆| 日本黄色片子视频| 欧美精品啪啪一区二区三区| 国产主播在线观看一区二区| 中国美女看黄片| 午夜精品一区二区三区免费看| av在线观看视频网站免费| 免费大片18禁| 日本在线视频免费播放| 亚洲第一区二区三区不卡| 欧美日韩中文字幕国产精品一区二区三区| 少妇人妻精品综合一区二区 | 内地一区二区视频在线| 亚洲经典国产精华液单 | 国产 一区 欧美 日韩| 12—13女人毛片做爰片一| 极品教师在线免费播放| 亚洲成a人片在线一区二区| 国产午夜精品论理片| 久久精品国产自在天天线| 日韩欧美精品免费久久 | 婷婷色综合大香蕉| 欧美xxxx黑人xx丫x性爽| 精品国内亚洲2022精品成人| 国产精品影院久久| 岛国在线免费视频观看| 国产在线精品亚洲第一网站| 91久久精品电影网| 18禁在线播放成人免费| 亚洲人成网站在线播| 欧美三级亚洲精品| 亚洲av日韩精品久久久久久密| 欧美中文日本在线观看视频| 精品福利观看| 有码 亚洲区| 国产成人欧美在线观看| 国产精品一区二区三区四区免费观看 | 老熟妇仑乱视频hdxx| 超碰av人人做人人爽久久| 亚洲国产精品久久男人天堂| 国产老妇女一区| 少妇的逼好多水| 桃色一区二区三区在线观看| 看十八女毛片水多多多| 国产精品,欧美在线| 欧美xxxx性猛交bbbb| 亚洲va日本ⅴa欧美va伊人久久| 在线观看一区二区三区| 久久久久九九精品影院| 赤兔流量卡办理| 精品久久久久久久久亚洲 | 成人午夜高清在线视频| 久久久精品大字幕| 2021天堂中文幕一二区在线观| 久久6这里有精品| 日韩大尺度精品在线看网址| av天堂在线播放| 一边摸一边抽搐一进一小说| 国产在线男女| www日本黄色视频网| 久久久色成人| 中文资源天堂在线| 少妇丰满av| 欧美色视频一区免费| 久久久精品大字幕| 亚洲精品成人久久久久久| 最近在线观看免费完整版| 免费av观看视频| 亚州av有码| 又爽又黄无遮挡网站| 女同久久另类99精品国产91| 国产精品,欧美在线| 亚洲成人久久爱视频| 亚洲成av人片在线播放无| 男女做爰动态图高潮gif福利片| 午夜精品在线福利| 欧美+日韩+精品| 久久99热这里只有精品18| 日本一二三区视频观看| 韩国av一区二区三区四区| 99久久精品热视频| 色视频www国产| 亚洲精品在线美女| 日本a在线网址| 真人做人爱边吃奶动态| av欧美777| 中文字幕精品亚洲无线码一区| 又紧又爽又黄一区二区| 噜噜噜噜噜久久久久久91| 亚洲人成网站在线播| 日韩欧美国产一区二区入口| 日本一二三区视频观看| 内射极品少妇av片p| 长腿黑丝高跟| 久久草成人影院| 欧美一区二区精品小视频在线| 欧美zozozo另类| 可以在线观看毛片的网站| 99久久精品热视频| 久久99热6这里只有精品| 亚洲国产高清在线一区二区三| 中文字幕熟女人妻在线| 亚洲国产精品久久男人天堂| 波多野结衣高清无吗| 男插女下体视频免费在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 亚州av有码| xxxwww97欧美| 在线观看免费视频日本深夜| 观看美女的网站| 夜夜躁狠狠躁天天躁| 一区福利在线观看| 可以在线观看的亚洲视频| 色av中文字幕| 9191精品国产免费久久| 高清日韩中文字幕在线| 国产精品伦人一区二区| 国产高清视频在线播放一区| 毛片一级片免费看久久久久 | 18禁黄网站禁片午夜丰满| 女同久久另类99精品国产91| 日本免费a在线| 小蜜桃在线观看免费完整版高清| 91麻豆精品激情在线观看国产| 美女高潮的动态| 熟女电影av网| 欧美在线一区亚洲| 9191精品国产免费久久| 88av欧美| 日本在线视频免费播放| 成人三级黄色视频| 欧美xxxx性猛交bbbb| 国产精品1区2区在线观看.| 日本免费a在线| 天堂网av新在线| 99久久精品热视频| 国产精品亚洲av一区麻豆| 1000部很黄的大片| 欧美成人a在线观看| 欧美绝顶高潮抽搐喷水| 精品午夜福利在线看| 中文字幕免费在线视频6| 日本五十路高清| 麻豆国产av国片精品| 变态另类成人亚洲欧美熟女| 制服丝袜大香蕉在线| 在线国产一区二区在线| 久久人人爽人人爽人人片va | 国内精品一区二区在线观看| 最近视频中文字幕2019在线8| 亚洲三级黄色毛片| 久久精品国产亚洲av香蕉五月| 一级av片app| av天堂在线播放| 熟女人妻精品中文字幕| 精品一区二区三区视频在线观看免费| 欧美黑人巨大hd| 搞女人的毛片| 亚洲狠狠婷婷综合久久图片| 亚洲三级黄色毛片| 免费电影在线观看免费观看| 欧美另类亚洲清纯唯美| 人人妻人人澡欧美一区二区| 别揉我奶头~嗯~啊~动态视频| 在线观看66精品国产| 男女之事视频高清在线观看| 日韩欧美在线乱码| 精品久久久久久久久av| 高清在线国产一区| 精华霜和精华液先用哪个| 亚洲人成电影免费在线| 国产极品精品免费视频能看的| 久久亚洲真实| 欧美不卡视频在线免费观看| 午夜亚洲福利在线播放| 麻豆一二三区av精品| 国产三级黄色录像| 99热精品在线国产| 亚洲熟妇中文字幕五十中出| 天堂影院成人在线观看| eeuss影院久久| 国产日本99.免费观看| 中文字幕av成人在线电影| 99精品久久久久人妻精品| 亚洲 欧美 日韩 在线 免费| 久久午夜福利片| 一级黄色大片毛片| 色播亚洲综合网| 成人国产综合亚洲| 亚洲欧美日韩卡通动漫| 亚洲,欧美,日韩| 极品教师在线免费播放| 精品久久久久久久人妻蜜臀av| АⅤ资源中文在线天堂| 99久久成人亚洲精品观看| 91av网一区二区| 婷婷丁香在线五月| 18禁黄网站禁片免费观看直播| 中出人妻视频一区二区| 精品人妻视频免费看| 可以在线观看的亚洲视频| 人妻夜夜爽99麻豆av| 97热精品久久久久久| 激情在线观看视频在线高清| 嫩草影院新地址| 老司机深夜福利视频在线观看| 最近视频中文字幕2019在线8| 国产亚洲精品综合一区在线观看| 床上黄色一级片| 国产aⅴ精品一区二区三区波| 精品国产三级普通话版| 精品日产1卡2卡| 亚洲乱码一区二区免费版| 性欧美人与动物交配| 亚洲乱码一区二区免费版| 亚洲成a人片在线一区二区| 99国产极品粉嫩在线观看| 三级毛片av免费| 男人舔奶头视频| 中文字幕精品亚洲无线码一区| 亚洲一区二区三区不卡视频| 宅男免费午夜| 最近最新中文字幕大全电影3| 波多野结衣高清作品| 国产色爽女视频免费观看| 欧美黄色片欧美黄色片| 高清在线国产一区| 精品乱码久久久久久99久播| 99久久成人亚洲精品观看| 国产免费av片在线观看野外av| 国产成人啪精品午夜网站| 久久精品夜夜夜夜夜久久蜜豆| 蜜桃久久精品国产亚洲av| 九九久久精品国产亚洲av麻豆| 最近最新中文字幕大全电影3| 国产成年人精品一区二区| 亚洲三级黄色毛片| 成年女人永久免费观看视频| 免费一级毛片在线播放高清视频| 午夜福利免费观看在线| 韩国av一区二区三区四区| 一个人免费在线观看电影| 亚洲欧美日韩东京热| 一个人观看的视频www高清免费观看| h日本视频在线播放| 成人无遮挡网站| 亚洲 欧美 日韩 在线 免费| 国产精华一区二区三区| 国内精品美女久久久久久| 成人av在线播放网站| 嫩草影视91久久| 亚洲成a人片在线一区二区| 精华霜和精华液先用哪个| 国产成人啪精品午夜网站| 91麻豆精品激情在线观看国产| 国产av一区在线观看免费| 精品乱码久久久久久99久播| 亚洲在线自拍视频| 久久人人爽人人爽人人片va | 亚洲 欧美 日韩 在线 免费| 免费观看的影片在线观看| 午夜福利高清视频| 国产淫片久久久久久久久 | 久久精品综合一区二区三区| 免费av观看视频| 亚洲av熟女| 99视频精品全部免费 在线| 成年人黄色毛片网站|