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

    爆震波通過環(huán)形通道傳播模式試驗(yàn)研究

    2021-02-15 10:33:34袁雪強(qiáng)蔣露欣劉世杰
    火箭推進(jìn) 2021年6期

    袁雪強(qiáng),蔣露欣,張 多,劉世杰

    (國防科技大學(xué) 高超聲速沖壓發(fā)動機(jī)技術(shù)重點(diǎn)試驗(yàn)室,湖南 長沙 410073)

    0 引言

    爆震燃燒作為一種近似于等容燃燒的燃燒形式,具有能量釋放速率快、熱力循環(huán)效率高等優(yōu)點(diǎn),因此以爆震燃燒為動力的爆震發(fā)動機(jī)在高超聲速推進(jìn)領(lǐng)域中有很好的應(yīng)用前景。其中,旋轉(zhuǎn)爆震發(fā)動機(jī)(rotating detonation engine,RDE)是爆震波以旋轉(zhuǎn)形式傳播為特點(diǎn)的動力裝置,這種發(fā)動機(jī)的燃燒室通常采用環(huán)形通道結(jié)構(gòu),爆震波在環(huán)形燃燒室內(nèi)沿周向旋轉(zhuǎn)傳播,形成的高溫氣體產(chǎn)物通過尾噴管排出從而產(chǎn)生推力。旋轉(zhuǎn)爆震發(fā)動機(jī)具有結(jié)構(gòu)簡單、質(zhì)量輕等優(yōu)勢,并利用爆震波傳播速度快的特點(diǎn),可通過一次點(diǎn)火使發(fā)動機(jī)的工作頻率達(dá)到數(shù)千赫茲,從而形成近乎連續(xù)的推力,因此是目前最有可能實(shí)現(xiàn)工程化應(yīng)用的爆震發(fā)動機(jī)形式。對于旋轉(zhuǎn)爆震發(fā)動機(jī)而言,環(huán)形燃燒室徑向尺寸的改變會顯著影響爆震波的傳播模態(tài)、傳播頻率以及波頭高度等參數(shù),是影響其推力性能的關(guān)鍵因素之一。為了深入探究徑向尺寸對爆震波傳播的影響機(jī)制,可將爆震波沿周向的旋轉(zhuǎn)傳播近似看作靜止氣體爆震波在環(huán)形通道內(nèi)的傳播過程,通過研究環(huán)形通道內(nèi)爆震波的傳播過程,可為旋轉(zhuǎn)爆震發(fā)動機(jī)環(huán)形燃燒室的設(shè)計提供有力指導(dǎo)。

    目前有關(guān)爆震波通過環(huán)形通道傳播過程的研究表明,受到內(nèi)外壁面的共同影響,爆震波在環(huán)形通道內(nèi)表現(xiàn)出多種不同的傳播模式。Thomas等最早通過試驗(yàn)發(fā)現(xiàn)了爆震波能夠在環(huán)形通道中穩(wěn)定地傳播,但同時也可能形成爆震波解耦和再起爆交替產(chǎn)生的不穩(wěn)定傳播。Lee等的仿真研究發(fā)現(xiàn),當(dāng)通道的曲率半徑小于某個臨界值時,爆震波表現(xiàn)為類似于單頭螺旋爆震的不穩(wěn)定傳播模式,而當(dāng)曲率半徑大于該臨界值時,爆震波能夠保持原有的橫波結(jié)構(gòu)穩(wěn)定地傳播。Kudo等通過試驗(yàn)觀測將爆震波的傳播模式按照波面是否發(fā)生解耦劃分為不穩(wěn)定傳播和穩(wěn)定傳播模式,并在對試驗(yàn)結(jié)果進(jìn)行統(tǒng)計后得到了兩種傳播模式發(fā)生轉(zhuǎn)換的內(nèi)壁面臨界曲率半徑。Nakayama等則基于Kudo的研究成果,研究將爆震波的傳播模式細(xì)分為不穩(wěn)定傳播模式、過渡傳播模式和穩(wěn)定傳播模式,并統(tǒng)計得到了3種模式之間發(fā)生相互轉(zhuǎn)換的內(nèi)壁面臨界曲率半徑分別為14

    λ

    和23

    λ

    ,其中

    λ

    為爆震胞格寬度。Pan等使用螺旋扁平方管對環(huán)形通道中爆震波的傳播實(shí)現(xiàn)了更為全面的觀測,所得結(jié)果與Nakayama的研究相似,同時還得到了形成爆震的極限曲率半徑在2.6

    λ

    到4.8

    λ

    之間。Sugiyama等和Yuan等分別使用兩步反應(yīng)和詳細(xì)反應(yīng)模型對不同曲率半徑與通道寬度比的環(huán)形通道下爆震波的傳播過程進(jìn)行了仿真,并運(yùn)用曲面爆震理論探討了傳播模式之間轉(zhuǎn)換的臨界條件。Short等和Xia等則在研究中發(fā)現(xiàn)了多種新的傳播模式,包括爆震波部分解耦后又恢復(fù)穩(wěn)定的不穩(wěn)定傳播模式以及曲面爆震波與馬赫干共同存在的穩(wěn)定傳播模式。

    從目前的相關(guān)研究來看,現(xiàn)有對環(huán)形通道內(nèi)爆震波的傳播模式的劃分方式較為簡單,其本質(zhì)上只考慮了內(nèi)壁面的爆震波衍射是否會使爆震波面發(fā)生解耦,無法全面體現(xiàn)爆震波受環(huán)形通道內(nèi)外壁面影響下的傳播特性,因此需建立更為詳細(xì)的劃分方法。基于上述問題,本文對爆震波通過環(huán)形通道的傳播現(xiàn)象進(jìn)行了試驗(yàn)觀測及分析,對環(huán)形通道內(nèi)爆震波的傳播模式進(jìn)行了更為系統(tǒng)詳細(xì)的劃分,并依次分析了各傳播模式的形成機(jī)制、波系結(jié)構(gòu)及速度特性,獲得了各參數(shù)對傳播模式分布的影響規(guī)律,所得結(jié)果可為環(huán)形燃燒室設(shè)計提供參考。

    1 試驗(yàn)系統(tǒng)及觀測方法

    本文的試驗(yàn)系統(tǒng)如圖1所示。該系統(tǒng)主要由點(diǎn)火裝置、爆震管和供氣系統(tǒng)3部分組成,其基本的工作過程為通過供氣系統(tǒng)對爆震管內(nèi)的氣體進(jìn)行抽吸,并充入試驗(yàn)用的預(yù)混氣,而后通過點(diǎn)火裝置點(diǎn)燃預(yù)混氣,點(diǎn)燃后的氣體在爆震管中發(fā)展為爆震波,并在試驗(yàn)段中通過觀測設(shè)備進(jìn)行觀測,最后從爆震管出口排出。

    圖1 彎曲爆震管試驗(yàn)系統(tǒng)組成示意圖Fig.1 Schematic of curved detonation tube experiment system

    彎曲爆震管結(jié)構(gòu)如圖2所示。整個爆震管可分為起爆段、過渡段、穩(wěn)定段、試驗(yàn)段和出口段,起爆段采用內(nèi)徑為20 mm的圓形直管,其前端安裝點(diǎn)火用的火花塞,后端連接過渡段,兩側(cè)分別開孔并與進(jìn)氣和排氣的管道相連。為了增加流場中的擾動,加速實(shí)現(xiàn)爆燃向爆震轉(zhuǎn)變的過程,在起爆段內(nèi)安裝了一段長度為400 mm的螺旋狀金屬擾流絲。過渡段為擴(kuò)張角較小的平滑擴(kuò)張構(gòu)型,以確保爆震波在傳入穩(wěn)定段的過程中不發(fā)生熄滅。穩(wěn)定段由截面尺寸為156 mm×20 mm,長度為800 mm的兩段方形直管連接而成,用于使爆震波在進(jìn)入試驗(yàn)段之前得到充分發(fā)展并實(shí)現(xiàn)穩(wěn)定傳播,其中靠近試驗(yàn)段的部分內(nèi)部可安裝長條插塊,以改變爆震波進(jìn)入試驗(yàn)段時的通道高度。試驗(yàn)段為內(nèi)部可以安裝插塊的方形箱體結(jié)構(gòu),通過更換內(nèi)部的插塊來改變通道構(gòu)型,其框架尺寸為280 mm×280 mm×82 mm,左側(cè)和頂部分別與穩(wěn)定段和出口段相連。插塊通過螺絲固定,再通過觀測玻璃蓋板進(jìn)行密封。玻璃蓋板內(nèi)部安裝石英玻璃,從而實(shí)現(xiàn)對流場的觀測和拍攝。由于觀測區(qū)域的限制,試驗(yàn)段中的環(huán)形通道只有90°。試驗(yàn)采用的氣體為298 K溫度下體積比為H∶O∶Ar=2∶1∶2的預(yù)混氣。試驗(yàn)初始壓力

    p

    分別設(shè)置為30 kPa、40 kPa和50 kPa,通道寬度

    d

    分別設(shè)置為20 mm、40 mm和60 mm,內(nèi)壁面的曲率半徑

    R

    設(shè)置為20~100 mm,間隔為20 mm。出口段尾部使用PVC膜進(jìn)行密封,爆震波形成后可通過破壞PVC膜排出,以減少反射波系對觀測的影響。

    圖2 彎曲爆震管和試驗(yàn)段結(jié)構(gòu)剖視圖Fig.2 Sectional view of curved detonation tube and experimental section

    激光陰影技術(shù)是本文試驗(yàn)的主要觀測手段,其測量光路如圖3所示。由于陰影反映的是流場密度的二階梯度,因此使用陰影技術(shù)能夠捕捉到流場中變化最為劇烈的爆震波波面信息。本試驗(yàn)中光源是由最大輸出功率為100 mW的半導(dǎo)體激光器發(fā)出的波長為532 nm的連續(xù)激光,拍攝爆震波傳播過程所使用的高速相機(jī)為Photron Fastcam SA-X2高速數(shù)字?jǐn)z影儀,并通過電腦端的操作軟件實(shí)現(xiàn)控制、圖像處理與數(shù)據(jù)存儲。本文試驗(yàn)采用外部觸發(fā)的拍攝模式,拍攝的圖像分辨率為256×256,拍攝速度為21 600 fps。另外,試驗(yàn)中還在相機(jī)鏡頭前安裝了只能通過532 nm激光的濾光片,以過濾燃燒產(chǎn)生的強(qiáng)烈火光。

    圖3 激光陰影原理示意圖Fig.3 Schematic of laser shadowgraph technique

    2 結(jié)果與分析

    過去的研究中,爆震波的傳播模式一般按照波面是否發(fā)生解耦劃分為不穩(wěn)定傳播模式和穩(wěn)定傳播模式。不穩(wěn)定傳播模式表現(xiàn)為傳播過程中爆震波面發(fā)生解耦,對應(yīng)的胞格圖上出現(xiàn)爆震胞格消失的區(qū)域,如圖4(a)、(b)所示,而穩(wěn)定模式則表現(xiàn)為爆震波在傳播過程中能始終保持波面耦合,對應(yīng)胞格圖中無胞格消失的區(qū)域,如圖4(c)、(d)所示。事實(shí)上,從環(huán)形通道內(nèi)外壁面對爆震波的影響機(jī)制來看,爆震波在環(huán)形通道內(nèi)能否穩(wěn)定傳播主要取決于內(nèi)壁面的衍射是否會使爆震波發(fā)生解耦,因此這樣的劃分方法忽略了外壁面附近馬赫干的發(fā)展對爆震波傳播的影響。為了更為深入地研究爆震波在環(huán)形通道內(nèi)的傳播模式及其特性,根據(jù)馬赫干的變化情況將現(xiàn)有的兩種傳播模式又細(xì)分為馬赫干增長型、馬赫干平穩(wěn)型和馬赫干衰減型,由此可得共計6種傳播模式。下面對這6種傳播模式的形成機(jī)制、波面結(jié)構(gòu)以及速度分布情況進(jìn)行詳細(xì)的分析說明。

    圖4 現(xiàn)有研究中的兩種傳播模式下的密度紋影和數(shù)值胞格圖對比Fig.4 Comparison of density schlieren and numerical soot foil for the two propagation modes in previous study

    2.1 不穩(wěn)定傳播模式

    2.1.1 馬赫干增長型不穩(wěn)定傳播模式

    馬赫干增長型不穩(wěn)定傳播模式(IMUP)不同時刻的密度陰影如圖5所示,對應(yīng)工況為

    p

    =30 kPa,

    d

    =40 mm,

    R

    =40 mm。從圖中可以看出,這種傳播模式的形成機(jī)制為內(nèi)壁面附近的爆震波面因衍射而發(fā)生解耦,而外壁面處形成的馬赫反射三波點(diǎn)與解耦波面發(fā)生接觸,由于馬赫干后方區(qū)域經(jīng)外壁面的壓縮溫度和壓力均較高,因此能夠點(diǎn)燃解耦區(qū)域內(nèi)的預(yù)混氣并形成局部起爆,起爆產(chǎn)生的能量會加速反射三波點(diǎn)向內(nèi)壁面的傳播,最終導(dǎo)致了馬赫干高度的迅速增長。該傳播模式的結(jié)構(gòu)特點(diǎn)是爆震波面上橫波結(jié)構(gòu)基本消失,但存在一個主導(dǎo)的馬赫反射三波點(diǎn)結(jié)構(gòu)在內(nèi)外壁面之間來回碰撞反射,且內(nèi)壁面附近的爆震波面始終存在解耦區(qū)域。

    圖5 馬赫干增長型不穩(wěn)定傳播模式(IMUP) 密度陰影圖Fig.5 Density shadowgraph of increasing Mach-stem unstable propagation (IMUP) mode

    2.1.2 馬赫干平穩(wěn)型不穩(wěn)定傳播模式

    馬赫干平穩(wěn)型不穩(wěn)定傳播模式(SMUP)的密度陰影如圖6所示,對應(yīng)工況為

    p

    =40 kPa,

    d

    =40 mm,

    R

    =40 mm。相較于IMUP模式,這種傳播模式的初始壓力

    p

    更高,所以內(nèi)壁面附近的衍射爆震波能夠自發(fā)形成再起爆,但還無法實(shí)現(xiàn)穩(wěn)定的傳播,因此波面上存在著有限的解耦區(qū)域。另外較高的初始壓力還導(dǎo)致爆震波誘導(dǎo)區(qū)和反應(yīng)區(qū)長度的縮短,從而使馬赫干高度增長緩慢,這就延長了馬赫反射三波點(diǎn)到達(dá)內(nèi)壁面的運(yùn)動距離,因此解耦區(qū)域與外壁面附近的馬赫反射三波點(diǎn)并未發(fā)生接觸,無法形成局部起爆以促使馬赫干的增長,而當(dāng)外壁面相對于爆震波面的夾角處于一定范圍內(nèi)時就可以形成高度基本不變的馬赫干。這種傳播模式的結(jié)構(gòu)特點(diǎn)為內(nèi)壁面附近衍射爆震波的再起爆與外壁面附近馬赫干的傳播相互獨(dú)立,互不影響。

    圖6 馬赫干平穩(wěn)型不穩(wěn)定傳播模式(SMUP) 密度陰影圖Fig.6 Density shadowgraph of steady Mach-stem unstable propagation (SMUP) mode

    2.1.3 馬赫干衰減型不穩(wěn)定傳播模式

    馬赫干衰減型不穩(wěn)定傳播模式(DMUP)的試驗(yàn)結(jié)果如圖7所示,對應(yīng)工況為

    p

    =40 kPa,

    d

    =60 mm,

    R

    =20 mm。該傳播模式的形成機(jī)理與SMUP模式相似,較高的

    p

    使內(nèi)壁面附近衍射爆震波能自發(fā)形成再起爆,但無法實(shí)現(xiàn)穩(wěn)定傳播,導(dǎo)致波面存在解耦區(qū),而由于

    R

    /d

    相比SMUP模式有所減小,因此不僅導(dǎo)致解耦區(qū)域與馬赫反射三波點(diǎn)無接觸,還相應(yīng)增大了外壁面相對于衍射波面的夾角,從而使馬赫干高度隨著傳播而逐漸發(fā)生衰減。最終馬赫干會完全消失,意味著馬赫反射向規(guī)則反射發(fā)生了轉(zhuǎn)化。這種傳播模式的結(jié)構(gòu)特點(diǎn)為發(fā)展穩(wěn)定后的爆震波面不存在馬赫干,外壁面處的反射為規(guī)則反射。

    圖7 馬赫干衰減型不穩(wěn)定傳播模式(DMUP) 密度陰影圖Fig.7 Density shadowgraph of decreasing Mach-stem unstable propagation (DMUP) mode

    2.2 穩(wěn)定傳播模式

    2.2.1 馬赫干增長型穩(wěn)定傳播模式

    圖8所示為馬赫干增長型穩(wěn)定傳播模式(IMSP)的密度陰影圖,對應(yīng)工況為

    p

    =30 kPa,

    d

    =20 mm,

    R

    =80 mm。該傳播模式的形成機(jī)制分為兩種情況:①如果內(nèi)壁面無法達(dá)到衍射爆震波穩(wěn)定傳播的條件,則衍射爆震波仍能夠維持一段距離的耦合狀態(tài)之后才發(fā)生解耦,而由于

    R

    /d

    較大,馬赫反射三波點(diǎn)到達(dá)內(nèi)壁面的運(yùn)動距離較短。因此,反射三波點(diǎn)可在衍射爆震波解耦之前通過在波面上的運(yùn)動到達(dá)內(nèi)壁面,從而保持爆震波不發(fā)生解耦地穩(wěn)定傳播。②如果衍射爆震波本身就能夠無解耦地穩(wěn)定傳播,那么只需保證反射三波點(diǎn)能夠通過運(yùn)動到達(dá)內(nèi)壁面即可實(shí)現(xiàn)馬赫干增長型的穩(wěn)定傳播。在這種傳播模式中,爆震波面實(shí)質(zhì)上是在內(nèi)外壁面來回碰撞的馬赫干,由于馬赫干的彎曲程度較小,且這種傳播模式多出現(xiàn)在通道寬度較小而曲率半徑較大的工況中,因此波面形態(tài)基本為垂直于內(nèi)外壁面的平直面。

    圖8 馬赫干增長型穩(wěn)定傳播模式(IMSP) 密度陰影圖Fig.8 Density shadowgraph of increasing Mach-stem stable propagation (IMSP) mode

    2.2.2 馬赫干平穩(wěn)型穩(wěn)定傳播模式

    圖9所示為馬赫干平穩(wěn)型穩(wěn)定傳播模式(SMSP)的密度陰影,對應(yīng)工況為

    p

    =40 kPa,

    d

    =40 mm,

    R

    =80 mm。相較于SMUP模式,該傳播模式下初始壓力

    p

    更高,因此衍射爆震波能夠在內(nèi)壁面附近實(shí)現(xiàn)無解耦的穩(wěn)定傳播,而較小的

    R

    /d

    能使外壁面附近形成高度基本不變的馬赫干。在這種傳播模式下,穩(wěn)定后的爆震波面表現(xiàn)為曲面爆震波與馬赫干組合的結(jié)構(gòu),且這種結(jié)構(gòu)在之后的傳播過程中基本保持恒定。

    圖9 馬赫干平穩(wěn)型穩(wěn)定傳播模式(SMSP) 密度陰影圖Fig.9 Density shadowgraph of steady Mach-stem stable propagation (SMSP) mode

    2.2.3 馬赫干衰減型穩(wěn)定傳播模式

    圖10所示為馬赫干衰減型穩(wěn)定傳播模式(DMSP)的試驗(yàn)拍攝結(jié)果,對應(yīng)工況為

    p

    =50 kPa,

    d

    =40 mm,

    R

    =40 mm。該傳播模式的形成機(jī)制與SMSP模式基本相同,足夠高的初始壓力

    p

    使內(nèi)壁面附近的衍射爆震波能夠?qū)崿F(xiàn)穩(wěn)定傳播,而

    R

    /d

    比SMSP模式進(jìn)一步減小,從而使波面與外壁面的夾角進(jìn)一步增大,導(dǎo)致馬赫干高度逐漸衰減,最終馬赫反射會向規(guī)則反射轉(zhuǎn)化,從而使馬赫干完全消失。這種模式的結(jié)構(gòu)特點(diǎn)為穩(wěn)定后的爆震波僅由平滑過渡的彎曲波面構(gòu)成,且該曲面結(jié)構(gòu)在之后傳播過程中基本保持恒定,外壁面處可見明顯的規(guī)則反射結(jié)構(gòu)。

    圖10 馬赫干衰減型穩(wěn)定傳播模式(DMSP) 密度陰影圖Fig.10 Density shadowgraph of decreasing Mach-stem stable propagation (DMSP) mode

    2.3 爆震波傳播速度分布規(guī)律

    通過對密度陰影中爆震波位置的提取,得到了各種傳播模式下爆震波在內(nèi)外壁面的傳播速度隨壁面中心夾角

    θ

    的分布情況,如圖11所示。圖中各傳播模式分別對應(yīng)于2.1節(jié)和2.2節(jié)中的各試驗(yàn)工況,其中

    D

    為爆震波面在壁面位置的傳播速度,

    D

    為該工況下的CJ爆震理論傳播速度。從內(nèi)壁面的傳播速度來看,在不穩(wěn)定傳播模式下內(nèi)壁面的傳播速度發(fā)生了明顯的衰減,各工況下的

    D

    /D

    均下降到了0.5左右,這表明爆震波解耦所造成的速度虧損非常嚴(yán)重。此外,在DMUP模式下觀測到傳播速度下降后上升的情況,推測可能為再起爆形成的局部爆震在內(nèi)壁面的碰撞使爆震波面重新耦合在一起,由此引發(fā)了短暫的速度躍升,這一現(xiàn)象在過去的研究中也有報道。相比之下,穩(wěn)定傳播模式由于內(nèi)壁面附近波面無解耦,因此其速度分布要穩(wěn)定許多,

    D

    /D

    基本保持在0.8~1之間,仍有一定的速度虧損,這是由于內(nèi)壁面附近爆震波面發(fā)生彎曲所導(dǎo)致的。

    圖11 各傳播模式下的爆震波傳播速度分布Fig.11 Propagation velocity evolution of detonation wave for different propagation modes

    對于外壁面的傳播速度,從圖中可以很明顯地看到,不論在何種傳播模式下,外壁面附近爆震波的傳播速度均大于內(nèi)壁面,且大于

    D

    ,處于過驅(qū)狀態(tài)。這是由于外壁面在垂直于放射三波點(diǎn)方向上的速度與

    D

    相當(dāng),由速度分解可知垂直于波面方向上的速度僅為爆震波傳播速度的一個分量,如圖12所示,假設(shè)馬赫干為平直波面且垂直于外壁面,則外壁面出的波面?zhèn)鞑ニ俣扰c反射三波點(diǎn)處速度相等,因此在外壁面上測得的傳播速度會大于

    D

    。另外,不穩(wěn)定模式和穩(wěn)定模式在外壁面的速度分布上沒有明顯的差異,其過驅(qū)度隨馬赫干的衰減而呈上升的趨勢,尤其是在馬赫干衰減型的傳播模式下,由于馬赫反射向規(guī)則反射發(fā)生了轉(zhuǎn)換,其過驅(qū)度顯著高于其他模式,達(dá)到2.0左右。此外還能發(fā)現(xiàn),同一波面在內(nèi)外壁面的速度相位差也與馬赫干的變化有關(guān),且也隨著馬赫干的衰減而增大,速度相位差主要由衍射而發(fā)生彎曲的爆震波面在內(nèi)外壁面的位置不同引起,而馬赫干的衰減會使彎曲波面的長度增加,從而導(dǎo)致了速度相位差的增大。

    圖12 環(huán)形通道內(nèi)爆震波面速度分布示意圖Fig.12 Schematic of detonation wave front velocity distribution in the annular channel

    以上的分析表明,所有的不穩(wěn)定傳播模式均會在內(nèi)壁面出現(xiàn)速度驟減和躍升的劇烈振蕩,且還伴隨著爆震波解耦和再起爆等不穩(wěn)定燃燒現(xiàn)象,如果旋轉(zhuǎn)爆震發(fā)動機(jī)環(huán)形燃燒室內(nèi)的爆震燃燒出現(xiàn)這樣的不穩(wěn)定傳播模式,不僅會嚴(yán)重影響爆震發(fā)動機(jī)的推力性能,還有可能導(dǎo)致燃燒室結(jié)構(gòu)的燒蝕與損壞。相比之下,穩(wěn)定傳播模式中內(nèi)外壁面的速度分布均較為平穩(wěn),能夠?qū)崿F(xiàn)較為固定的旋轉(zhuǎn)傳播頻率,且燃料的燃燒也更為充分,因此在對環(huán)形燃燒室的尺寸參數(shù)進(jìn)行設(shè)計時,應(yīng)當(dāng)盡量避免燃燒室內(nèi)出現(xiàn)不穩(wěn)定傳播,確保爆震波在穩(wěn)定傳播模式下平穩(wěn)地旋轉(zhuǎn)傳播。此外,爆震波面較大的速度相位差可能導(dǎo)致瞬時推力沿徑向發(fā)生偏移,而較大過驅(qū)度所形成的高溫高壓則可能導(dǎo)致燃燒室壁面發(fā)生燒蝕。

    綜合來看,IMSP模式下內(nèi)外壁面的傳播速度均較為平穩(wěn),且均接近

    D

    ,速度相位差也很小,是環(huán)形燃燒室內(nèi)爆震波較為理想的傳播模式。

    2.4 爆震波傳播模式變化規(guī)律

    通過對所有試驗(yàn)數(shù)據(jù)的統(tǒng)計,可得爆震波在不同

    p

    、

    d

    R

    條件下的傳播模式分布如圖13所示,從其分布中可以得到爆震波傳播模式的變化規(guī)律。

    圖13 所有試驗(yàn)工況下爆震波傳播模式分布Fig.13 Distribution of detonation propagation modes for all experiment cases

    首先對于不穩(wěn)定傳播模式與穩(wěn)定傳播模式之間的轉(zhuǎn)換,這主要由內(nèi)壁面附近爆震波面是否發(fā)生解耦決定,因此初始壓力

    p

    的增大提高了爆震波的強(qiáng)度,從而防止其在衍射過程中發(fā)生解耦,因此穩(wěn)定傳播模式的分布范圍擴(kuò)大,這與過去的研究是相符的。而隨著曲率半徑

    R

    的增大,傳播模式也會從不穩(wěn)定向穩(wěn)定轉(zhuǎn)變,這是由于內(nèi)壁面附近衍射波的稀疏效應(yīng)隨

    R

    的增大而減弱,從而降低了爆震解耦的發(fā)生。但通道寬度

    d

    的改變似乎對不穩(wěn)定模式與穩(wěn)定模式之間的轉(zhuǎn)換沒有顯著的影響,這一現(xiàn)象可通過曲面爆震理論進(jìn)行有效解釋。對于馬赫干的變化類型,則主要由外壁面附近的爆震波馬赫反射過程所決定。當(dāng)

    p

    增大時,馬赫干會呈現(xiàn)增長—平穩(wěn)—衰減的轉(zhuǎn)換趨勢,這是由于

    p

    的增大會縮短爆震波誘導(dǎo)區(qū)和反應(yīng)區(qū)的長度,從而加劇馬赫干的衰減。

    R

    的增大能夠增大反射時馬赫干的絕對高度,而通道寬度

    d

    的減小則能夠縮短反射三波點(diǎn)到達(dá)內(nèi)壁面的距離,因此

    R

    /d

    的增大會使馬赫干呈現(xiàn)衰減—平穩(wěn)—增長的轉(zhuǎn)換趨勢。

    最后,基于2.1節(jié)和2.2節(jié)中對爆震波各種傳播模式形成機(jī)制的分析,可形成圖14所示的各傳播模式工況范圍的判別流程,通過該流程就能夠?qū)Σ煌r下環(huán)形通道內(nèi)爆震波的傳播模式進(jìn)行判別,從而獲得每種傳播模式的工況分布范圍,如果可以從理論上確定流程中每個物理過程發(fā)生的臨界條件,則可以對環(huán)形通道中爆震波的傳播模式進(jìn)行有效的預(yù)測,這將在今后開展進(jìn)一步的深入研究。

    圖14 環(huán)形通道內(nèi)爆震波傳播模式判別流程圖Fig.14 Flowchart for distinguishing different propagation modes of detonation wave in annular channel

    3 結(jié)論

    本文對爆震波通過環(huán)形通道的傳播現(xiàn)象進(jìn)行了試驗(yàn)研究。根據(jù)爆震波受內(nèi)外壁面的影響機(jī)制對其傳播模式進(jìn)行了詳細(xì)的劃分,并依次分析了各傳播模式的形成機(jī)制、結(jié)構(gòu)及速度特性,獲得了傳播模式分布隨參數(shù)的變化規(guī)律,主要結(jié)論如下:

    1)爆震波傳播模式按照波面是否發(fā)生解耦可劃分為不穩(wěn)定傳播模式和穩(wěn)定傳播模式,其中不穩(wěn)定傳播模式下內(nèi)壁面附近的爆震波會發(fā)生解耦,從而導(dǎo)致內(nèi)壁面附近爆震波速度出現(xiàn)虧損和振蕩,相比之下,穩(wěn)定傳播模式中爆震波不會發(fā)生解耦,因此在內(nèi)外壁面附近的速度分布均更加平穩(wěn)。

    2)根據(jù)馬赫干的變化情況可將傳播模式細(xì)分為馬赫干增長型、平穩(wěn)型和衰減型。外壁面附近波面的過驅(qū)度隨馬赫干的衰減呈現(xiàn)上升的趨勢,尤其在馬赫干衰減型的傳播模式下,爆震波會在外壁面形成規(guī)則反射,導(dǎo)致外壁面附近波面過驅(qū)度遠(yuǎn)高于其他模式,同時內(nèi)外壁面速度相位差也隨馬赫干的衰減而增大。

    3)隨著初始壓力和壁面曲率半徑的增大,爆震波實(shí)現(xiàn)穩(wěn)定傳播的能力得到增強(qiáng),但其穩(wěn)定傳播能力與通道寬度無關(guān)。此外,初始壓力的增大使馬赫干呈現(xiàn)增長—平穩(wěn)—衰減的轉(zhuǎn)換趨勢,而壁面曲率半徑與通道寬度比的增大會縮短馬赫反射三波點(diǎn)到達(dá)內(nèi)壁面的運(yùn)動距離,使馬赫干呈現(xiàn)衰減—平穩(wěn)—增長的變化規(guī)律。

    高潮久久久久久久久久久不卡| 国产精品一及| 国产黄a三级三级三级人| 亚洲av熟女| 国产精品美女特级片免费视频播放器 | aaaaa片日本免费| 亚洲精品中文字幕一二三四区| 日韩高清综合在线| 成年人黄色毛片网站| 丝袜人妻中文字幕| 国产97色在线日韩免费| 黄色a级毛片大全视频| 在线十欧美十亚洲十日本专区| 一区二区三区高清视频在线| 999久久久国产精品视频| 欧美av亚洲av综合av国产av| 日韩三级视频一区二区三区| 长腿黑丝高跟| 香蕉国产在线看| 国模一区二区三区四区视频 | 成年人黄色毛片网站| 又爽又黄无遮挡网站| 亚洲av电影不卡..在线观看| 国产1区2区3区精品| 国产精品亚洲一级av第二区| 大型黄色视频在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品影院6| 亚洲黑人精品在线| 黄色片一级片一级黄色片| 法律面前人人平等表现在哪些方面| 国产精品久久电影中文字幕| 午夜精品一区二区三区免费看| 人人妻人人看人人澡| 9191精品国产免费久久| 无限看片的www在线观看| 韩国av一区二区三区四区| 亚洲午夜理论影院| 国产一级毛片七仙女欲春2| 亚洲国产精品sss在线观看| 亚洲18禁久久av| 午夜福利免费观看在线| 啪啪无遮挡十八禁网站| 国产成人aa在线观看| 亚洲国产精品久久男人天堂| 国产av一区在线观看免费| 亚洲午夜理论影院| 亚洲中文字幕一区二区三区有码在线看 | 搡老熟女国产l中国老女人| 亚洲欧美精品综合一区二区三区| 欧美+亚洲+日韩+国产| 宅男免费午夜| 男女之事视频高清在线观看| 中国美女看黄片| 成年人黄色毛片网站| 日韩欧美在线二视频| 两个人看的免费小视频| 久久精品国产综合久久久| 黄色成人免费大全| 一二三四在线观看免费中文在| 亚洲精品在线观看二区| 久久精品国产清高在天天线| 日日干狠狠操夜夜爽| 国产精品美女特级片免费视频播放器 | 国产成年人精品一区二区| 久久人妻福利社区极品人妻图片| 精品免费久久久久久久清纯| 精品久久久久久久末码| 在线观看免费午夜福利视频| 欧美极品一区二区三区四区| 日日摸夜夜添夜夜添小说| 婷婷亚洲欧美| 18禁观看日本| 国产蜜桃级精品一区二区三区| 一个人免费在线观看的高清视频| 亚洲精品久久国产高清桃花| 美女 人体艺术 gogo| 日韩欧美在线二视频| 女人爽到高潮嗷嗷叫在线视频| 精品国产超薄肉色丝袜足j| 村上凉子中文字幕在线| 一二三四社区在线视频社区8| 精品国产乱码久久久久久男人| 亚洲aⅴ乱码一区二区在线播放 | 曰老女人黄片| 亚洲成人免费电影在线观看| 久久精品人妻少妇| 一区福利在线观看| 亚洲精品美女久久久久99蜜臀| 国产一区二区三区在线臀色熟女| 日韩大码丰满熟妇| 欧美3d第一页| 久久久久久久久中文| 国产一级毛片七仙女欲春2| 91字幕亚洲| 精品国产美女av久久久久小说| 亚洲熟妇熟女久久| 岛国在线免费视频观看| www.www免费av| 成人国语在线视频| 精品久久久久久久久久久久久| 女人爽到高潮嗷嗷叫在线视频| 亚洲一区二区三区色噜噜| 国产一区二区在线观看日韩 | 国产人伦9x9x在线观看| 香蕉久久夜色| 999久久久国产精品视频| 一进一出好大好爽视频| 日日干狠狠操夜夜爽| 可以在线观看的亚洲视频| 色尼玛亚洲综合影院| 特级一级黄色大片| 999精品在线视频| 黄色女人牲交| 亚洲精品粉嫩美女一区| 麻豆国产97在线/欧美 | 欧美日韩福利视频一区二区| 免费在线观看成人毛片| 女人被狂操c到高潮| 国产高清videossex| 久久午夜综合久久蜜桃| bbb黄色大片| 在线观看美女被高潮喷水网站 | 哪里可以看免费的av片| 亚洲成人免费电影在线观看| 国产精品亚洲美女久久久| 特大巨黑吊av在线直播| 日韩大码丰满熟妇| 少妇裸体淫交视频免费看高清 | 黄片小视频在线播放| 国产精品一区二区精品视频观看| 久久久久国产一级毛片高清牌| 大型av网站在线播放| 亚洲成人国产一区在线观看| 欧美日韩黄片免| 国产1区2区3区精品| 亚洲成a人片在线一区二区| 啪啪无遮挡十八禁网站| АⅤ资源中文在线天堂| 99国产极品粉嫩在线观看| 国产伦在线观看视频一区| 欧美色欧美亚洲另类二区| 亚洲精品中文字幕一二三四区| 99久久精品热视频| 国产精品久久久久久久电影 | 好看av亚洲va欧美ⅴa在| 日本一本二区三区精品| 99热6这里只有精品| 亚洲黑人精品在线| 精华霜和精华液先用哪个| 亚洲电影在线观看av| 亚洲精品av麻豆狂野| 国产激情偷乱视频一区二区| 亚洲一区二区三区不卡视频| 亚洲精品一卡2卡三卡4卡5卡| 两个人的视频大全免费| 亚洲乱码一区二区免费版| 精品久久久久久久久亚洲| 免费看a级黄色片| 国产精品久久久久久精品电影小说 | 欧美bdsm另类| 国产高清激情床上av| 久久热精品热| 日韩中字成人| 黄色视频,在线免费观看| 欧美xxxx黑人xx丫x性爽| 美女 人体艺术 gogo| 亚洲av成人精品一区久久| 亚洲在久久综合| 免费看a级黄色片| 高清日韩中文字幕在线| 国产成人精品婷婷| 中文字幕熟女人妻在线| 国产成人影院久久av| 婷婷色av中文字幕| 中文字幕久久专区| 内射极品少妇av片p| 一区二区三区高清视频在线| 亚洲精品久久久久久婷婷小说 | 日韩欧美精品免费久久| 中国美女看黄片| 国产视频内射| 天堂av国产一区二区熟女人妻| 人妻夜夜爽99麻豆av| 搡女人真爽免费视频火全软件| 长腿黑丝高跟| 久久午夜福利片| 免费看av在线观看网站| 成人综合一区亚洲| 97人妻精品一区二区三区麻豆| 亚洲国产精品成人综合色| 91精品国产九色| 99热全是精品| 国模一区二区三区四区视频| 一级黄片播放器| 99久久成人亚洲精品观看| 麻豆成人av视频| 久久99精品国语久久久| 舔av片在线| 韩国av在线不卡| 欧美高清成人免费视频www| 日韩一本色道免费dvd| 日本撒尿小便嘘嘘汇集6| 亚洲在线观看片| 精品一区二区三区人妻视频| 亚洲一级一片aⅴ在线观看| 亚洲最大成人手机在线| 啦啦啦观看免费观看视频高清| 身体一侧抽搐| www.色视频.com| 亚洲精品影视一区二区三区av| 中文字幕人妻熟人妻熟丝袜美| 国产爱豆传媒在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲欧美成人精品一区二区| 人妻少妇偷人精品九色| 在线a可以看的网站| 国产一区二区三区在线臀色熟女| 精品久久久久久久久久免费视频| 少妇被粗大猛烈的视频| 国产中年淑女户外野战色| 欧美最新免费一区二区三区| 一个人免费在线观看电影| 国产日韩欧美在线精品| 美女大奶头视频| 日本欧美国产在线视频| 一个人看的www免费观看视频| 国产精品三级大全| 又黄又爽又刺激的免费视频.| 亚洲成人精品中文字幕电影| 1024手机看黄色片| 人体艺术视频欧美日本| 色播亚洲综合网| 国产激情偷乱视频一区二区| 日韩三级伦理在线观看| 亚洲国产精品久久男人天堂| 男女啪啪激烈高潮av片| 老司机影院成人| 看免费成人av毛片| 偷拍熟女少妇极品色| 国产日韩欧美在线精品| 国产蜜桃级精品一区二区三区| 国产视频内射| 日韩av不卡免费在线播放| videossex国产| av在线亚洲专区| 黑人高潮一二区| 国产v大片淫在线免费观看| 国产极品精品免费视频能看的| a级毛片免费高清观看在线播放| 久久精品夜色国产| 午夜视频国产福利| 一区二区三区高清视频在线| 深夜精品福利| 三级毛片av免费| 国产伦理片在线播放av一区 | 成年女人永久免费观看视频| 99热全是精品| 国产精品乱码一区二三区的特点| 国内揄拍国产精品人妻在线| 好男人视频免费观看在线| 久久人妻av系列| 一区二区三区高清视频在线| 淫秽高清视频在线观看| 成人二区视频| 少妇高潮的动态图| 91麻豆精品激情在线观看国产| 老司机影院成人| 尾随美女入室| 国产精品av视频在线免费观看| 午夜福利视频1000在线观看| 国产欧美日韩精品一区二区| 亚洲av免费在线观看| 特大巨黑吊av在线直播| 亚洲欧美日韩无卡精品| 麻豆精品久久久久久蜜桃| 亚洲国产高清在线一区二区三| а√天堂www在线а√下载| 波多野结衣巨乳人妻| 97超视频在线观看视频| 乱码一卡2卡4卡精品| 亚洲成人精品中文字幕电影| 日韩一区二区三区影片| 亚洲精品乱码久久久久久按摩| 亚洲精品国产av成人精品| 久久精品国产亚洲av涩爱 | 黄色配什么色好看| 精品免费久久久久久久清纯| 亚洲,欧美,日韩| 久久久精品94久久精品| 欧美最黄视频在线播放免费| 中文精品一卡2卡3卡4更新| 国产精品永久免费网站| 男人舔女人下体高潮全视频| 国产精品蜜桃在线观看 | 精品欧美国产一区二区三| 九草在线视频观看| 一区二区三区四区激情视频 | 夜夜爽天天搞| 老熟妇乱子伦视频在线观看| 成人午夜高清在线视频| 久久精品国产亚洲av涩爱 | 网址你懂的国产日韩在线| 超碰av人人做人人爽久久| 欧美变态另类bdsm刘玥| 丰满乱子伦码专区| 麻豆一二三区av精品| 深夜精品福利| 国产日韩欧美在线精品| 五月玫瑰六月丁香| 天堂网av新在线| 久久精品国产清高在天天线| 国产精品久久电影中文字幕| 一区福利在线观看| 中文字幕免费在线视频6| 久久99蜜桃精品久久| 久久精品夜色国产| 亚洲精品亚洲一区二区| 在线免费观看不下载黄p国产| 成人毛片60女人毛片免费| 性色avwww在线观看| 亚洲经典国产精华液单| 两个人视频免费观看高清| 欧美日韩国产亚洲二区| 国产一区二区三区在线臀色熟女| 久久久欧美国产精品| 国产黄a三级三级三级人| 国产又黄又爽又无遮挡在线| 99热只有精品国产| 国产日韩欧美在线精品| 长腿黑丝高跟| 国产精品电影一区二区三区| 男女边吃奶边做爰视频| 可以在线观看毛片的网站| 亚洲熟妇中文字幕五十中出| 日韩大尺度精品在线看网址| АⅤ资源中文在线天堂| 99久久精品一区二区三区| 波多野结衣高清作品| 老女人水多毛片| 菩萨蛮人人尽说江南好唐韦庄 | 日韩视频在线欧美| 亚洲av成人av| 国产激情偷乱视频一区二区| 国产一区二区激情短视频| 国产色婷婷99| 男女边吃奶边做爰视频| 国产精品人妻久久久久久| 色综合站精品国产| 97超碰精品成人国产| 国产精品三级大全| 久久久欧美国产精品| 国国产精品蜜臀av免费| 波多野结衣巨乳人妻| 亚洲国产欧美人成| 国产精品人妻久久久久久| 九九热线精品视视频播放| av专区在线播放| 哪里可以看免费的av片| 欧美变态另类bdsm刘玥| 黄色视频,在线免费观看| 国产亚洲精品av在线| 最近视频中文字幕2019在线8| 国产成人午夜福利电影在线观看| 亚洲最大成人av| 亚洲精华国产精华液的使用体验 | 网址你懂的国产日韩在线| 欧美日本亚洲视频在线播放| 欧美丝袜亚洲另类| 蜜臀久久99精品久久宅男| 人妻系列 视频| 2022亚洲国产成人精品| 又爽又黄无遮挡网站| 成人三级黄色视频| 亚洲av不卡在线观看| eeuss影院久久| 老司机福利观看| avwww免费| 国产精品久久电影中文字幕| 美女国产视频在线观看| 国产亚洲av片在线观看秒播厂 | 永久网站在线| 亚洲欧美中文字幕日韩二区| 99国产极品粉嫩在线观看| a级毛色黄片| 亚洲乱码一区二区免费版| 女人被狂操c到高潮| 国产 一区精品| 国产精品久久久久久久久免| 久久精品国产亚洲av香蕉五月| 99久久人妻综合| 国产成人91sexporn| 中文在线观看免费www的网站| 国产精品,欧美在线| 中文在线观看免费www的网站| 日韩欧美三级三区| 国产欧美日韩精品一区二区| 99九九线精品视频在线观看视频| 久久久久国产网址| 性插视频无遮挡在线免费观看| 三级男女做爰猛烈吃奶摸视频| 国模一区二区三区四区视频| 国产一级毛片七仙女欲春2| 精品一区二区三区视频在线| 久久欧美精品欧美久久欧美| 嘟嘟电影网在线观看| 久久精品国产自在天天线| 日本三级黄在线观看| 91麻豆精品激情在线观看国产| 国产色爽女视频免费观看| 亚洲国产欧美在线一区| 欧美变态另类bdsm刘玥| 日本av手机在线免费观看| 国产精品99久久久久久久久| 毛片一级片免费看久久久久| 欧美日韩一区二区视频在线观看视频在线 | 亚洲图色成人| 在线免费观看的www视频| 中文字幕久久专区| 舔av片在线| 一个人看视频在线观看www免费| 成人午夜高清在线视频| 插逼视频在线观看| 亚洲欧美中文字幕日韩二区| 欧美日本亚洲视频在线播放| 中文字幕免费在线视频6| 国产精品久久久久久精品电影| 国产精品蜜桃在线观看 | a级毛片a级免费在线| 最近2019中文字幕mv第一页| h日本视频在线播放| 男女那种视频在线观看| 久久鲁丝午夜福利片| 舔av片在线| 99热精品在线国产| 亚洲av男天堂| 国产91av在线免费观看| 欧美日本亚洲视频在线播放| 国产精品一区www在线观看| 天堂中文最新版在线下载 | 午夜老司机福利剧场| 久久久久国产网址| 免费大片18禁| 亚洲精华国产精华液的使用体验 | 欧美最黄视频在线播放免费| 亚洲一级一片aⅴ在线观看| 高清毛片免费看| 日韩中字成人| 男人狂女人下面高潮的视频| 日本五十路高清| 国产人妻一区二区三区在| 永久网站在线| 国产精华一区二区三区| 熟女电影av网| 在线观看午夜福利视频| 免费无遮挡裸体视频| 丝袜美腿在线中文| www.av在线官网国产| 91在线精品国自产拍蜜月| 欧美最新免费一区二区三区| 亚洲国产精品成人久久小说 | 欧美高清性xxxxhd video| 国产精品免费一区二区三区在线| 天堂av国产一区二区熟女人妻| 免费大片18禁| 在线观看免费视频日本深夜| 能在线免费观看的黄片| 亚洲精品成人久久久久久| 国产精品永久免费网站| 精品久久国产蜜桃| 亚洲欧美日韩东京热| 一边摸一边抽搐一进一小说| 午夜福利高清视频| 一级av片app| 久久久久久久久大av| 成人欧美大片| 美女xxoo啪啪120秒动态图| 久久久久久久久中文| 成年女人永久免费观看视频| 国产免费男女视频| 97热精品久久久久久| 人人妻人人澡人人爽人人夜夜 | 看非洲黑人一级黄片| 一本精品99久久精品77| 淫秽高清视频在线观看| 亚洲色图av天堂| 欧美色视频一区免费| 久久99蜜桃精品久久| 久久久久久久久久久丰满| 国产黄a三级三级三级人| 久久精品国产亚洲av涩爱 | 亚洲精品久久久久久婷婷小说 | 99精品在免费线老司机午夜| 国产黄片视频在线免费观看| 国产午夜福利久久久久久| 一级黄片播放器| 久久亚洲精品不卡| 99精品在免费线老司机午夜| 国产高清有码在线观看视频| 精品一区二区三区人妻视频| 男人狂女人下面高潮的视频| 国产精品免费一区二区三区在线| 国产精品不卡视频一区二区| 中国美女看黄片| 91午夜精品亚洲一区二区三区| 中文字幕制服av| 欧美精品国产亚洲| 精品少妇黑人巨大在线播放 | 男女那种视频在线观看| 亚洲最大成人av| 夫妻性生交免费视频一级片| 热99re8久久精品国产| 久久精品久久久久久噜噜老黄 | 日本免费一区二区三区高清不卡| 老司机福利观看| 长腿黑丝高跟| 岛国毛片在线播放| 亚洲欧美中文字幕日韩二区| 少妇丰满av| 99精品在免费线老司机午夜| 女人十人毛片免费观看3o分钟| 色哟哟哟哟哟哟| 嘟嘟电影网在线观看| 亚洲18禁久久av| 国产免费男女视频| 亚洲精品国产av成人精品| av天堂在线播放| 久久久久久久久久久丰满| 久久久精品94久久精品| 人人妻人人澡人人爽人人夜夜 | 中文字幕熟女人妻在线| 欧美日韩综合久久久久久| 免费在线观看成人毛片| 国产精品一及| 亚洲性久久影院| 国产精华一区二区三区| 99国产极品粉嫩在线观看| 别揉我奶头 嗯啊视频| 亚洲无线在线观看| 在线观看66精品国产| 秋霞在线观看毛片| 美女国产视频在线观看| 尾随美女入室| 久久国内精品自在自线图片| 黄片无遮挡物在线观看| 国产成人精品一,二区 | 美女大奶头视频| 国产色婷婷99| 日韩av在线大香蕉| 亚洲精品乱码久久久久久按摩| 99久久人妻综合| 久久久精品欧美日韩精品| 美女脱内裤让男人舔精品视频 | 全区人妻精品视频| 少妇人妻精品综合一区二区 | 精品日产1卡2卡| 欧美精品国产亚洲| 日本熟妇午夜| 精品人妻熟女av久视频| 少妇被粗大猛烈的视频| 麻豆一二三区av精品| 国模一区二区三区四区视频| 欧美日本亚洲视频在线播放| 国产色爽女视频免费观看| 校园人妻丝袜中文字幕| 精品熟女少妇av免费看| 麻豆av噜噜一区二区三区| 亚洲一级一片aⅴ在线观看| 国产午夜精品论理片| 99久国产av精品国产电影| 亚洲乱码一区二区免费版| 成年av动漫网址| 精品久久国产蜜桃| 在线观看66精品国产| 三级毛片av免费| 我的女老师完整版在线观看| 一个人看视频在线观看www免费| 麻豆成人午夜福利视频| 一区福利在线观看| 97超碰精品成人国产| 99riav亚洲国产免费| 国产 一区精品| 中文字幕久久专区| 男女视频在线观看网站免费| 男人舔奶头视频| 国产高清三级在线| av在线播放精品| 亚洲激情五月婷婷啪啪| 老女人水多毛片| 91在线精品国自产拍蜜月| 天天躁夜夜躁狠狠久久av| 国内久久婷婷六月综合欲色啪| 一级二级三级毛片免费看| 国产精品一区二区性色av| 日本三级黄在线观看| 深爱激情五月婷婷| 国产女主播在线喷水免费视频网站 | 午夜亚洲福利在线播放| 我的老师免费观看完整版| 别揉我奶头 嗯啊视频| 一级毛片久久久久久久久女| 亚洲一区二区三区色噜噜| 九九爱精品视频在线观看| 成人漫画全彩无遮挡| 亚洲欧美精品自产自拍| 亚洲av不卡在线观看| 狂野欧美激情性xxxx在线观看| 亚洲欧美精品自产自拍| 内地一区二区视频在线| 深爱激情五月婷婷| 日韩欧美 国产精品| 日韩欧美国产在线观看| av在线播放精品| 国产极品天堂在线| 国产成人a区在线观看|