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

    乙烯/空氣在激波管中自點(diǎn)火流場(chǎng)顯示研究

    2019-01-31 00:36:44劉二偉徐勝利
    實(shí)驗(yàn)流體力學(xué) 2019年1期
    關(guān)鍵詞:觀察窗陣面激波

    劉二偉, 徐勝利

    (1. 中國(guó)科學(xué)技術(shù)大學(xué) 工程科學(xué)學(xué)院, 合肥 230026; 2. 清華大學(xué) 航天航空學(xué)院, 北京 100084)

    0 引 言

    和傳統(tǒng)加壓加熱方式相比,激波管能提供實(shí)驗(yàn)時(shí)間1~10ms、溫度和壓力均勻的近似一維流場(chǎng),避免了加熱的復(fù)雜性及其產(chǎn)生的傳熱效應(yīng),同時(shí)壓力和溫度變化范圍較大[1-4],因此,激波管常用于研究燃料自點(diǎn)火特性,包括燃料點(diǎn)火延時(shí)和自點(diǎn)火火焰?zhèn)鞑ヌ匦?。其中,自點(diǎn)火流場(chǎng)圖像具有三維形貌特征,顯示了反應(yīng)區(qū)位置和尺度以及表面不同尺度漩渦結(jié)構(gòu)等唯象特征,不僅可以和點(diǎn)火延時(shí)、火焰?zhèn)鞑ニ俣鹊葴y(cè)量數(shù)據(jù)相互驗(yàn)證,從而優(yōu)化測(cè)量方法,而且能定性考察壓力、溫度和當(dāng)量比等參數(shù)對(duì)火焰形成和發(fā)展過程的影響,也為燃燒數(shù)值模擬研究提供了更多對(duì)比信息。自點(diǎn)火流場(chǎng)顯示研究通常選擇矩形截面激波管(或在圓形截面激波管上添加過渡段),以便在激波管壁安裝玻璃觀察窗,從而利用高速攝影相機(jī)進(jìn)行拍攝。限于相機(jī)性能,已有自點(diǎn)火流場(chǎng)顯示研究不能充分反映流場(chǎng)特征。如李蘭等[5]采用高速攝影相機(jī)拍攝了RP-3煤油點(diǎn)火圖像,曝光時(shí)間高達(dá)43μs,圖像存在強(qiáng)烈的時(shí)間積分效應(yīng)。廖欽等[6]和王高峰等[7]采用ICCD相機(jī)拍攝了大量煤油點(diǎn)火圖像,但單次激波管實(shí)驗(yàn)只能獲得一張圖像,需要多次重復(fù)實(shí)驗(yàn)才能得到不同延時(shí)的流場(chǎng)圖像,這就依賴于激波管流場(chǎng)的重復(fù)性。Herzler等[8]采用ICCD相機(jī)拍攝了激波管截面丙烷點(diǎn)火圖像,但工況較少。因此,系統(tǒng)拍攝不同工況下的自點(diǎn)火流場(chǎng)圖像是有必要的。

    作為煤油催化裂解產(chǎn)物之一,乙烯由于結(jié)構(gòu)簡(jiǎn)單且含有不飽和鍵(C=C鍵),基元反應(yīng)特征明顯,是研究碳?xì)淙剂宵c(diǎn)火延時(shí)和化學(xué)動(dòng)力學(xué)機(jī)理的常用氣態(tài)燃料[9]。國(guó)內(nèi)外已有大量乙烯點(diǎn)火研究,相關(guān)工況范圍較為完善。如Barker等[10]測(cè)量了當(dāng)量比0.13~2.00、C2H4體積百分比0.25%~2.00%、壓力0.3~1.2MPa、溫度1050~1550K的C2H4/O2/Ar點(diǎn)火延時(shí)。Suzuki等[11]測(cè)量了溫度800~1400K的C2H4/O2/Ar點(diǎn)火延時(shí),發(fā)現(xiàn)高溫和低溫工況具有不同的活化能。Brown等[12]測(cè)量了當(dāng)量比1.0、稀釋氣體濃度75%~96%、壓力0.13~0.50MPa、溫度1073~2211K的C2H4/O2/Ar和C2H4/O2/N2點(diǎn)火延時(shí),分析了不同稀釋氣體對(duì)乙烯點(diǎn)火延時(shí)的影響,還研究了乙烯點(diǎn)火過程中的爆燃轉(zhuǎn)爆轟(DDT)現(xiàn)象。Cadman等[13]測(cè)量了當(dāng)量比1和1.5、壓力0.2~0.6MPa、溫度800~1620K的C2H4/O2/N2點(diǎn)火延時(shí)。Saxena等[14]測(cè)量了當(dāng)量比1和3、稀釋氣體濃度93%~98%、壓力0.2~1.8MPa、溫度1000~1650K的C2H4/O2/Ar點(diǎn)火延時(shí),并擬合得到了適用范圍較廣的乙烯點(diǎn)火延時(shí)關(guān)系式。胡宏浩等[9]、梁金虎等[15]和Deng等[16]在點(diǎn)火延時(shí)實(shí)驗(yàn)測(cè)量的基礎(chǔ)上,分別考察了污染物H2O和CO2、N2O和O2對(duì)乙烯點(diǎn)火延時(shí)的影響,并和多種化學(xué)動(dòng)力學(xué)機(jī)理預(yù)測(cè)結(jié)果進(jìn)行了對(duì)比。因此,選擇乙烯作為拍攝對(duì)象,不僅可以研究自點(diǎn)火現(xiàn)象,也有利于對(duì)比已有點(diǎn)火延時(shí)數(shù)據(jù),從而驗(yàn)證實(shí)驗(yàn)測(cè)量精度。

    綜上所述,為考察溫度和壓力等參數(shù)對(duì)自點(diǎn)火流場(chǎng)的影響,在矩形激波管上拍攝不同工況下的乙烯/空氣自點(diǎn)火流場(chǎng)序列圖像并進(jìn)行分析研究。

    1 實(shí)驗(yàn)裝置和測(cè)量方法

    1.1 矩形激波管和氣路系統(tǒng)

    圖1給出了矩形激波管和氣路系統(tǒng)示意圖,詳見文獻(xiàn)[6]。激波管截面尺寸為130mm(高)×80mm (寬),壁厚10mm,相鄰內(nèi)壁面圓倒角為R10mm,高、低壓段長(zhǎng)度分別為2和2.6m。實(shí)驗(yàn)段長(zhǎng)0.4m,觀察窗尺寸為120mm(長(zhǎng))×80mm(高),其材質(zhì)為熔融石英,右側(cè)和激波管反射端面重合。

    圖1 激波管和氣路系統(tǒng)示意圖Fig.1 Sketch of shock tube and gas distribution

    乙烯/空氣預(yù)混氣在混氣罐(約0.2945m3)中配制并靜置12h后使用。高壓段充入氦氣和氮?dú)饣旌衔?,壓力采用精密?shù)字壓力表(BSK811,量程0~1MPa,精度0.05%)測(cè)量。低壓段充入乙烯/空氣預(yù)混氣,采用數(shù)字真空計(jì)(CPCA-130Z,量程0~10kPa,精度1Pa)顯示充氣壓力。破膜方式采用針刺破膜,膜片材質(zhì)為BOPP(厚度為50和100μm)。實(shí)驗(yàn)前,通過真空泵組(Leybold TRIVAC旋片式真空泵和RUVAC羅茨真空泵)將高、低壓段抽真空至余壓約60Pa,漏氣率約10Pa/min。

    1.2 測(cè)量系統(tǒng)

    圖2給出了實(shí)驗(yàn)段和測(cè)量系統(tǒng)示意圖。PCB1、PCB2和PCB3為壓電傳感器(PCB 112A),和反射端面距離分別為657、237和20mm,壓電信號(hào)經(jīng)電荷放大器(YE5853)放大后分別送入示波器(Tektronix 2014C)。光纖與PCB3和反射端面距離相同。光信號(hào)經(jīng)光纖(Ocean optics FIB-600-DUV,波段200~1100nm,芯徑600μm)送入光電倍增管(Hamamatsu R928,波段185~900nm,側(cè)窗型PMT)。改變PMT偏置電壓,可改變PMT輸出電壓幅值,偏置電壓通常設(shè)置為-800V。光信號(hào)經(jīng)PMT轉(zhuǎn)換為電壓信號(hào)送入示波器,由三通將該信號(hào)同時(shí)送入八通道延時(shí)器DG645,從而輸出指定延時(shí)TTL電平到高速ICCD相機(jī)(Stanford Computer Optics,XXRapidFrame with 3 Nos. 4 Picos)外觸發(fā)端口。XXRapidFrame相機(jī)包含4個(gè)獨(dú)立ICCD相機(jī)和同步控制系統(tǒng),各ICCD相機(jī)可單獨(dú)設(shè)置拍攝延時(shí),因此,單次激波管實(shí)驗(yàn)可得到4張不同延時(shí)的自點(diǎn)火流場(chǎng)圖像。ICCD相機(jī)增益范圍為600至1000,為提高圖像亮度,統(tǒng)一設(shè)置為1000。

    圖2 自點(diǎn)火流場(chǎng)信號(hào)測(cè)量系統(tǒng)示意圖Fig.2 Sketch of signal detection for auto-ignition flow field

    需要說明的是,XXRapidFrame相機(jī)要求外觸發(fā)信號(hào)上升沿在短時(shí)間內(nèi)(<10μs)達(dá)到1.6V,考慮到低溫點(diǎn)火條件下的PMT信號(hào)上升較為緩慢,難以保證每次可靠地觸發(fā)XXRapidFrame相機(jī),而DG645觸發(fā)電平則可以任意設(shè)置,因此,統(tǒng)一采用PMT信號(hào)來觸發(fā)DG645,再由DG645觸發(fā)XXRapidFrame相機(jī)。

    圖3給出了PMT信號(hào)、DG645輸出信號(hào)、XXRapidFrame相機(jī)拍攝延時(shí)和曝光時(shí)間時(shí)序圖。其中,Vt為DG645觸發(fā)電平,通常設(shè)為Vt=-0.2V。Δt1為DG645輸出延時(shí),Δt2為XXRapidFrame拍攝延時(shí),通常設(shè)為Δt1=Δt2=0;ti為幅間隔,考慮到典型反射激波速度約為400m/s,通常設(shè)為ti=25μs,對(duì)部分低溫點(diǎn)火條件,ti=50μs;te為曝光時(shí)間,通常設(shè)te=1μs,以凍結(jié)流場(chǎng)。由上述延時(shí)設(shè)置,得到序列圖像拍攝延時(shí)為Δt=Δt1+Δt2+(n-1)·ti。

    圖3 測(cè)量系統(tǒng)信號(hào)時(shí)序圖Fig.3 Sequence of used signals

    2 結(jié)果分析與討論

    2.1 自點(diǎn)火流場(chǎng)圖像重復(fù)性

    表1給出了2個(gè)相同初始工況的實(shí)驗(yàn)參數(shù)列表,其中,φ為乙烯/空氣預(yù)混氣當(dāng)量比,Mas為入射激波馬赫數(shù),τig為點(diǎn)火延時(shí),定義為反射激波到達(dá)測(cè)點(diǎn)(PCB3)時(shí)刻和自發(fā)光信號(hào)(PMT)最大斜率起始點(diǎn)時(shí)刻的時(shí)間差。從表1可以看出:case 1和case 2的Mas、p5和T5偏差小于2%,可近似認(rèn)為2個(gè)工況參數(shù)相同。

    表1 相同工況參數(shù)列表(case 1和2)Table 1 Parameters for the same cases (case 1 and 2)

    需要說明的是:5區(qū)混合物定壓比熱和溫度、各組元質(zhì)量百分?jǐn)?shù)有關(guān)。理想激波關(guān)系假設(shè)化學(xué)反應(yīng)是凍結(jié)的,且定壓比熱和溫度無關(guān)。因此,由理想激波關(guān)系得到的p5和T5,與考慮溫度、組元質(zhì)量百分?jǐn)?shù)時(shí)的數(shù)值是不同的。表1給出了后者(由Gaseq程序計(jì)算)得到的p5和T5,而由理想激波關(guān)系得到的p5和T5分別為104kPa和1373K(case 1)、 104kPa和1397K(case 2)。以T5為例,分別相差163K(case 1)和168K(case 2)。胡宏浩等[9]、梁金虎等[15]的研究指出,采用溫度和組元質(zhì)量百分?jǐn)?shù)修正的p5和T5更合理。下文計(jì)算均采用Gaseq程序計(jì)算給出的p5和T5。

    圖4給出了case 1對(duì)應(yīng)的自點(diǎn)火流場(chǎng)序列圖像。為便于圖像分析,取x和y軸為水平和豎直方向。激波管反射端面位于x=0。光纖和PCB3位于x=20mm,其位置用虛線表示?;鹧嫖恢萌=40mm的火焰右邊沿。反射激波自左向右傳播,用實(shí)線表示。由入射激波實(shí)測(cè)速度,結(jié)合激波關(guān)系和圖像延時(shí),可計(jì)算出反射激波位置。具體地,由PMT信號(hào)得到Vt=-0.2V對(duì)應(yīng)時(shí)刻t1,由PCB3信號(hào)得到反射激波到達(dá)時(shí)刻t2,由入射激波速度得到馬赫數(shù),再根據(jù)激波關(guān)系,計(jì)算得到反射激波速度ur,則圖4的反射激波位置xr=20+ur(t1-t2+Δt)。

    需要說明的是:(1) 對(duì)比觀察窗和激波管截面尺寸可知,觀察窗上下邊沿距管壁25mm,因此,自點(diǎn)火流場(chǎng)圖像不包含邊界層;(2) 受場(chǎng)地限制,實(shí)驗(yàn)過程中激波管高壓段在右側(cè),低壓段在左側(cè),因此,自點(diǎn)火流場(chǎng)圖像和圖1水平方向相反。

    圖4(a)表明:當(dāng)Δt=0,在反射端面附近,即4mm40mm)厚度較均勻(約6mm),下部(y<40mm)厚度逐漸變小(約3mm)。注意到火焰尚未到達(dá)光纖測(cè)點(diǎn),但相機(jī)已開始拍攝,原因是壁面附近存在“λ”形分叉激波,波后火焰已到達(dá)光纖測(cè)點(diǎn),而DG645觸發(fā)電平設(shè)置較低(-0.2V),火焰亮度較弱但足夠觸發(fā)DG645,從而啟動(dòng)相機(jī)拍攝。圖4(b)表明:當(dāng)Δt=25μs,火焰?zhèn)鞑ブ?5mm63mm)略微向下游(左側(cè))傾斜,下部(y<17mm)出現(xiàn)彎曲形狀?;鹧嫦掠稳源嬖诹炼容^低且沿y軸分布較均勻的區(qū)域,這表明速率較低的化學(xué)反應(yīng)仍在進(jìn)行中。距火焰越遠(yuǎn)的下游區(qū)域亮度越低,光強(qiáng)趨近于0,這表明當(dāng)?shù)鼗瘜W(xué)反應(yīng)趨于結(jié)束。圖4(c)表明:當(dāng)Δt=50μs,火焰?zhèn)鞑ブ?0mm

    圖5給出了case 2對(duì)應(yīng)的自點(diǎn)火流場(chǎng)圖像。由圖可知:火焰產(chǎn)生和傳播特征與case 1類似。不同的是:當(dāng)Δt=50μs,火焰陣面和反射激波位置近似重合(圖5(c))。當(dāng)Δt=75μs,火焰已越過反射激波,出現(xiàn)這種偏差的原因是反射激波位置計(jì)算未考慮化學(xué)反應(yīng)產(chǎn)生的當(dāng)?shù)貕毫Σ▽?duì)反射激波的推動(dòng)作用,即忽略了反射激波和火焰之間的耦合作用。同時(shí),t2實(shí)際為分叉激波波頭到達(dá)測(cè)點(diǎn)時(shí)刻,不能準(zhǔn)確表示反射激波中心部位的位置。實(shí)際上,由PCB2和PCB3測(cè)得的p-t曲線得到實(shí)測(cè)值ur為916m/s,遠(yuǎn)大于由激波關(guān)系得到的計(jì)算值ur(416m/s),因此,圖4和5的反射激波位置存在較大誤差,下一步實(shí)驗(yàn)研究應(yīng)采用紋影法顯示反射激波的準(zhǔn)確位置。實(shí)際現(xiàn)象應(yīng)該是:當(dāng)火焰和反射激波位置接近或重合時(shí),兩者耦合向前傳播,形成波前有來流的右行爆轟波。

    圖4 自點(diǎn)火流場(chǎng)圖像(case 1)Fig.4 Images of auto-ignition flow field (case 1)

    圖5 自點(diǎn)火流場(chǎng)圖像(case 2)Fig.5 Images of auto-ignition flow field (case 2)

    根據(jù)圖4和5火焰位置和圖像幅間隔時(shí)間,表2給出了計(jì)算得到的火焰?zhèn)鞑ソ^對(duì)速度u和相對(duì)2區(qū)氣流速度Δu。其中,下標(biāo)a、b、c和d對(duì)應(yīng)圖像編號(hào)。表2表明:和case 1相比,case 2的ubc和ucd較大,這是由于case 2的T5較大,導(dǎo)致化學(xué)反應(yīng)速率較高,使得火焰?zhèn)鞑ニ俣仍龃?。另外,u隨時(shí)間推進(jìn)并不是恒定值。造成誤差的主要原因是:圖4和5火焰都不是嚴(yán)格的平面火焰,火焰陣面存在渦結(jié)構(gòu),在不同位置取“測(cè)速點(diǎn)”,就會(huì)對(duì)應(yīng)著不同的火焰?zhèn)鞑ニ俣取H簟皽y(cè)速點(diǎn)”左右偏差1mm,又ti=25μs,則u計(jì)算值相差40m/s。

    表2 火焰平均傳播速度(case 1和2)Table 2 Mean velocities of flame propagation (case 1 and 2)

    綜上所述,當(dāng)p5=106kPa、T5=1210K和p5=108kPa、T5=1229K,火焰產(chǎn)生和傳播特征近似相同。自點(diǎn)火首先發(fā)生在反射端面附近,形成近似平面火焰并右行傳播,火焰表面存在小尺度渦結(jié)構(gòu),即火焰在傳播過程中表面出現(xiàn)失穩(wěn)現(xiàn)象,火焰陣面逐漸垂直于激波管軸線。

    2.2 溫度影響

    表3給出了壓力近似不變、不同溫度的工況參數(shù)列表。從表3看出,φ值不變(1.0),p5為104±6kPa,T5為1294~1077K,對(duì)應(yīng)τig為12~728μs。

    圖6~9給出了和表3對(duì)應(yīng)的自點(diǎn)火流場(chǎng)圖像。

    表3 不同溫度工況參數(shù)列表(case 3~6)Table 3 Parameters for cases at different temperatures (case 3~6)

    圖6表明:當(dāng)T5=1294K,火焰產(chǎn)生和傳播特征和圖4、5類似。

    圖7表明:當(dāng)T5=1138K,點(diǎn)火首先發(fā)生在激波管反射端面附近,形成右行傳播火焰(圖7(a))。隨時(shí)間推進(jìn),火焰厚度增大(圖7(a)~(c))。當(dāng)Δt=0,反射端面附近(x<10mm)出現(xiàn)亮度較低火焰。當(dāng)Δt=50μs,火焰分布在x<20mm區(qū)域。觀察窗上側(cè)(x=30mm,70mm

    圖8表明:當(dāng)T5=1109K,點(diǎn)火同時(shí)出現(xiàn)在激波管反射端面附近和遠(yuǎn)離端面區(qū)域,形成多處離散火焰(圖8(a))。隨時(shí)間推進(jìn),離散火焰逐漸融合后繼續(xù)右行傳播(圖8(b)~(d))。和圖7相比,火焰范圍明顯增大。當(dāng)Δt=0,觀察窗內(nèi)出現(xiàn)多個(gè)不規(guī)則離散火焰,空間分布范圍較廣(x=10、30和70mm)。當(dāng)Δt=50μs,火焰向周圍擴(kuò)散,仍呈離散分布且尚未融合。當(dāng)Δt=100μs,不同區(qū)域火焰已融合(0

    圖9表明:當(dāng)T5=1077K,點(diǎn)火首先發(fā)生在遠(yuǎn)離激波管反射端面處(圖9(a)),火焰左行傳播(圖9(b)~(d))。當(dāng)Δt=0,右側(cè)流場(chǎng)(x>20mm)出現(xiàn)火焰區(qū),左側(cè)流場(chǎng)(x<20mm)無明顯火焰(亮度接近于0),表明左側(cè)流場(chǎng)未出現(xiàn)自點(diǎn)火。當(dāng)Δt=50μs,左側(cè)流場(chǎng)出現(xiàn)火焰區(qū),且存在明顯的孔洞結(jié)構(gòu)。當(dāng)Δt=100μs,火焰區(qū)擴(kuò)散至整個(gè)觀察窗區(qū)域,且右側(cè)流場(chǎng)亮度較高。當(dāng)Δt=150μs,左側(cè)流場(chǎng)的火焰亮度增大。這和文獻(xiàn)給出的自點(diǎn)火流場(chǎng)圖像是類似的,如廖欽[6]得到的煤油及其裂解產(chǎn)物的自點(diǎn)火圖像,Herzler等[8]得到的丙烷自點(diǎn)火圖像,Pfahl等[17]得到的異辛烷自點(diǎn)火陰影圖像。這些圖像均表明:低溫自點(diǎn)火不是從激波反射端面開始的。原因可能是:(1) 當(dāng)點(diǎn)火延時(shí)較長(zhǎng),反射端面附近流場(chǎng)和壁面存在傳熱效應(yīng)。(2) 考慮到氣體粘性效應(yīng),反射端面附近入射激波后的氣流不是瞬時(shí)滯止,存在反射激波和入射激波誘導(dǎo)邊界層的復(fù)雜相互作用,導(dǎo)致邊界層的低溫反應(yīng)物和反射端面附近存在輸運(yùn)關(guān)系。兩者同時(shí)作用,降低了反射端面附近的化學(xué)反應(yīng)速率,但其具體機(jī)理尚待驗(yàn)證。

    圖6 自點(diǎn)火流場(chǎng)圖像(case 3)Fig.6 Images of auto-ignition flow field (case 3)

    圖7 自點(diǎn)火流場(chǎng)圖像(case 4)Fig.7 Images of auto-ignition flow field (case 4)

    綜上所述,當(dāng)p5和φ相同,T5降低導(dǎo)致火焰區(qū)厚度增大、光強(qiáng)減弱,初始火焰形狀由近似平面火焰演化為不規(guī)則的多個(gè)離散火焰區(qū),且自點(diǎn)火區(qū)域也遠(yuǎn)離反射端面。在低溫條件下,火焰首先出現(xiàn)在遠(yuǎn)離反射端面處并左行傳播。

    2.3 壓力影響

    表4給出了不同壓力對(duì)應(yīng)的工況參數(shù)列表(case 7~14),其中,case 7和8、case 9和10、case 11和12、case 13和14分別為對(duì)比工況。從表4看出,p5分別為265±11kPa和419±44kPa,對(duì)應(yīng)T5分別為1252~1079K和1384~1055K。這表明:當(dāng)p5增大,p5和T5偏差范圍也增大。原因是:隨著p5增大,膜片厚度也增大,實(shí)驗(yàn)重復(fù)性變差,導(dǎo)致p5和對(duì)比工況的T5偏差增大,例如,case 7和8的T5相差132K。

    圖10~17給出了對(duì)應(yīng)表4的自點(diǎn)火流場(chǎng)圖像。

    圖10和11表明:隨p5增大,對(duì)應(yīng)的反應(yīng)物濃度增大,火焰光強(qiáng)明顯增大,火焰陣面為明顯的平面形狀,渦結(jié)構(gòu)近似消失,同時(shí)火焰下游流場(chǎng)光強(qiáng)也變大,表明燃燒仍在劇烈進(jìn)行。

    圖12表明:當(dāng)Δt=0,反射端面附近(x<8mm)出現(xiàn)亮度較低的近似平面火焰。當(dāng)Δt=25μs,火焰分布在x<20mm,上部(y=50mm)和下部(y=20mm)均存在渦結(jié)構(gòu)。當(dāng)Δt=50μs,上部和下部漩渦區(qū)分別形成向外傳播的火焰面,并相互穿透(20mm

    表4 不同壓力工況參數(shù)列表(case 7~14)Table 4 Parameters for cases at different pressures (case 7~14)

    圖10 自點(diǎn)火流場(chǎng)圖像(case 7)Fig.10 Images of auto-ignition flow field (case 7)

    圖11 自點(diǎn)火流場(chǎng)圖像(case 8)Fig.11 Images of auto-ignition flow field (case 8)

    圖12 自點(diǎn)火流場(chǎng)圖像(case 9)Fig.12 Images of auto-ignition flow field (case 9)

    圖13表明:當(dāng)Δt=0,反射端面附近(x<10mm)出現(xiàn)亮度較低的近似平面火焰。當(dāng)Δt=25μs,火焰分布在x<20mm?;鹧骊嚸娣譃樯喜亢拖虏?個(gè)半圓形狀,在y=45mm處相連。當(dāng)Δt=50μs,火焰分布在x<80mm,表現(xiàn)為2個(gè)圓弧交叉形狀。當(dāng)Δt=75μs,火焰分布在x<117mm,火焰陣面近似垂直于激波管軸線。和圖12相比,火焰光強(qiáng)明顯增大,且上部和下部出現(xiàn)更明顯的局部爆炸現(xiàn)象,并近似沿激波管軸線對(duì)稱,最終形成統(tǒng)一的火焰陣面。

    在圖14和15中,火焰的產(chǎn)生和傳播特征與圖8類似。不同的是:當(dāng)p5=375kPa,觀察窗下部(x=20mm,y=20mm)出現(xiàn)局部爆炸(圖15(b)),火焰向外傳播并擴(kuò)散到反應(yīng)較慢區(qū)域(圖15(c)),火焰陣面趨于垂直激波管軸線(圖15(d))。

    圖16和17表明:火焰的產(chǎn)生和傳播特征與圖9類似。不同的是:當(dāng)p5=444kPa,觀察窗上部(x=40mm,y=70mm)出現(xiàn)局部爆炸并左行傳播(圖17(c)),使反射端面附近出現(xiàn)較強(qiáng)的自發(fā)光(圖17(d))。

    圖13 自點(diǎn)火流場(chǎng)圖像(case 10)Fig.13 Images of auto-ignition flow field (case 10)

    圖14 自點(diǎn)火流場(chǎng)圖像(case 11)Fig.14 Images of auto-ignition flowfield (case 11)

    圖15 自點(diǎn)火流場(chǎng)圖像(case 12)Fig.15 Images of auto-ignition flow field (case 12)

    圖16 自點(diǎn)火流場(chǎng)圖像(case 13)Fig.16 Images of auto-ignition flow field (case 13)

    圖17 自點(diǎn)火流場(chǎng)圖像(case 14)Fig.17 Images of auto-ignition flow field (case 14)

    綜上所述,當(dāng)給定p5,T5對(duì)火焰的影響近似相同;當(dāng)T5不變,隨p5增大,火焰光強(qiáng)增大,化學(xué)反應(yīng)更劇烈;當(dāng)p5=265和419kPa,易發(fā)生局部爆炸并向周圍傳播,多個(gè)局部爆炸區(qū)在傳播過程中逐漸融合,經(jīng)壁面反射后形成近似平面火焰。

    3 結(jié) 論

    (1) 相同工況自點(diǎn)火流場(chǎng)圖像表明:當(dāng)p5=106kPa、T5=1210K和p5=108kPa、T5=1229K,火焰的產(chǎn)生和傳播特征近似相同。點(diǎn)火首先發(fā)生在反射端面附近,形成近似平面火焰并右行傳播,火焰表面存在小尺度渦結(jié)構(gòu)。

    (2) 當(dāng)p5和φ相同,不同T5自點(diǎn)火流場(chǎng)圖像表明:當(dāng)T5降低,火焰厚度增大,光強(qiáng)減弱,初始火焰由平面形狀演變?yōu)椴灰?guī)則的離散區(qū)域,點(diǎn)火位置逐漸遠(yuǎn)離激波管端面。當(dāng)T5較低,火焰首先出現(xiàn)在觀察窗右側(cè)并左行傳播。

    (3) 當(dāng)T5和φ相同,不同p5自點(diǎn)火流場(chǎng)圖像表明:隨p5增大,火焰光強(qiáng)增大,化學(xué)反應(yīng)更劇烈。當(dāng)p5=265和419kPa,火焰出現(xiàn)局部爆炸并向周圍傳播,多個(gè)局部爆炸區(qū)在傳播過程中逐漸融合,經(jīng)壁面反射后形成近似平面火焰。

    猜你喜歡
    觀察窗陣面激波
    錐臺(tái)形觀察窗有機(jī)玻璃結(jié)構(gòu)特性分析與優(yōu)化設(shè)計(jì)
    大型柔性陣面陣架動(dòng)力學(xué)分析*
    直升機(jī)觀察窗屈曲分析
    基于相鄰一維線陣干涉儀陣面的測(cè)向補(bǔ)償算法研究
    一種基于聚類分析的二維激波模式識(shí)別算法
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    斜激波入射V形鈍前緣溢流口激波干擾研究
    適于可壓縮多尺度流動(dòng)的緊致型激波捕捉格式
    耐壓結(jié)構(gòu)觀察窗應(yīng)力分析
    陣面分布不均勻的相控陣天線維修優(yōu)化模型
    亚洲中文字幕一区二区三区有码在线看| 69av精品久久久久久| 免费av毛片视频| 精品日产1卡2卡| 亚洲成人中文字幕在线播放| 久久久久性生活片| 成人毛片a级毛片在线播放| 久久久国产成人免费| 日韩欧美 国产精品| 午夜爱爱视频在线播放| 日本黄大片高清| 麻豆国产av国片精品| 一进一出抽搐gif免费好疼| 成人av在线播放网站| 亚洲激情五月婷婷啪啪| 少妇熟女欧美另类| 99国产极品粉嫩在线观看| 免费看av在线观看网站| 国产 一区精品| 午夜精品国产一区二区电影 | 女同久久另类99精品国产91| 久久午夜福利片| 午夜亚洲福利在线播放| 国产伦精品一区二区三区视频9| 人人妻,人人澡人人爽秒播| 亚洲精品久久国产高清桃花| 又爽又黄a免费视频| 在线免费十八禁| 国产亚洲精品久久久com| 久久久精品94久久精品| 国产精品人妻久久久影院| 亚洲色图av天堂| 亚洲精品影视一区二区三区av| 欧美区成人在线视频| 给我免费播放毛片高清在线观看| 国产成人精品久久久久久| 熟女人妻精品中文字幕| 亚洲av中文av极速乱| 日本免费一区二区三区高清不卡| 精品熟女少妇av免费看| 在线天堂最新版资源| 99riav亚洲国产免费| 亚洲电影在线观看av| 欧美色视频一区免费| 22中文网久久字幕| 中文字幕av在线有码专区| 国产午夜精品论理片| 99久久无色码亚洲精品果冻| 久久国内精品自在自线图片| 成年女人看的毛片在线观看| 能在线免费观看的黄片| 亚洲精品在线观看二区| 老司机福利观看| 亚洲精品色激情综合| 中文字幕av在线有码专区| 两个人视频免费观看高清| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美日韩高清专用| 日日啪夜夜撸| 别揉我奶头 嗯啊视频| 亚洲熟妇中文字幕五十中出| 日韩欧美在线乱码| 欧美zozozo另类| 寂寞人妻少妇视频99o| aaaaa片日本免费| 午夜a级毛片| 天堂网av新在线| 久久综合国产亚洲精品| 亚洲人成网站在线播放欧美日韩| 成人性生交大片免费视频hd| 深夜a级毛片| 亚洲在线观看片| 天美传媒精品一区二区| 亚洲人成网站在线播放欧美日韩| 久久精品夜色国产| 搡女人真爽免费视频火全软件 | 久久久久久久亚洲中文字幕| 亚洲精品国产成人久久av| 少妇高潮的动态图| 99在线视频只有这里精品首页| 国产精品无大码| 日本一本二区三区精品| 美女cb高潮喷水在线观看| 国产成人aa在线观看| avwww免费| 成人精品一区二区免费| 欧美激情在线99| 国产真实伦视频高清在线观看| 精品一区二区免费观看| 亚洲国产精品sss在线观看| 不卡一级毛片| 亚洲最大成人中文| 成人精品一区二区免费| 天堂av国产一区二区熟女人妻| 可以在线观看的亚洲视频| 性色avwww在线观看| 亚洲精品影视一区二区三区av| 最好的美女福利视频网| 大又大粗又爽又黄少妇毛片口| 69av精品久久久久久| 午夜福利在线在线| 国产高清三级在线| 日韩一区二区视频免费看| 日本免费a在线| 悠悠久久av| av在线亚洲专区| 熟女电影av网| 亚洲成人av在线免费| 在线观看av片永久免费下载| 大型黄色视频在线免费观看| 人人妻,人人澡人人爽秒播| 欧美日韩一区二区视频在线观看视频在线 | 亚洲av不卡在线观看| 免费大片18禁| 久久久午夜欧美精品| 99精品在免费线老司机午夜| 亚洲va在线va天堂va国产| 日韩中字成人| 久久久久久久久久成人| 最近的中文字幕免费完整| 国产男人的电影天堂91| av天堂中文字幕网| 如何舔出高潮| 久久综合国产亚洲精品| 男女啪啪激烈高潮av片| 精品一区二区三区av网在线观看| 久久韩国三级中文字幕| 国产精品福利在线免费观看| 精品久久久久久久人妻蜜臀av| 一本久久中文字幕| 久久这里只有精品中国| 国产免费男女视频| 麻豆乱淫一区二区| 成年女人永久免费观看视频| 日本-黄色视频高清免费观看| 九九爱精品视频在线观看| 观看美女的网站| 日本欧美国产在线视频| 色视频www国产| 国产精品久久久久久久久免| 成年av动漫网址| 国产中年淑女户外野战色| 97超视频在线观看视频| 国产视频一区二区在线看| 男人狂女人下面高潮的视频| 亚洲精品国产成人久久av| 日本免费一区二区三区高清不卡| 精品久久久噜噜| 美女cb高潮喷水在线观看| 日日撸夜夜添| 日本成人三级电影网站| 内地一区二区视频在线| 欧美zozozo另类| 久久精品国产自在天天线| 日韩,欧美,国产一区二区三区 | 婷婷精品国产亚洲av在线| 色吧在线观看| 国产成人福利小说| 中国美女看黄片| 少妇猛男粗大的猛烈进出视频 | 国产 一区 欧美 日韩| 欧美在线一区亚洲| eeuss影院久久| 亚洲国产精品合色在线| 老司机午夜福利在线观看视频| a级一级毛片免费在线观看| 成人永久免费在线观看视频| 国产视频内射| 国产精品精品国产色婷婷| 亚洲精品在线观看二区| 亚洲精华国产精华液的使用体验 | av视频在线观看入口| avwww免费| 欧美日韩综合久久久久久| 看黄色毛片网站| 欧美区成人在线视频| 又黄又爽又刺激的免费视频.| 亚洲精品国产av成人精品 | 久久99热6这里只有精品| 秋霞在线观看毛片| 亚洲av.av天堂| 国内精品美女久久久久久| 啦啦啦观看免费观看视频高清| 又黄又爽又免费观看的视频| 亚洲久久久久久中文字幕| 亚洲高清免费不卡视频| 99久久中文字幕三级久久日本| 在线天堂最新版资源| 免费观看精品视频网站| 最近在线观看免费完整版| 淫秽高清视频在线观看| 亚洲成人精品中文字幕电影| 少妇高潮的动态图| 老女人水多毛片| 欧美色视频一区免费| 久久久精品欧美日韩精品| 无遮挡黄片免费观看| 国产毛片a区久久久久| 中文在线观看免费www的网站| 精品欧美国产一区二区三| 欧美成人精品欧美一级黄| 91狼人影院| 高清毛片免费看| 毛片女人毛片| 国产一区二区亚洲精品在线观看| 色吧在线观看| 亚洲av成人av| 少妇猛男粗大的猛烈进出视频 | 国产v大片淫在线免费观看| 久久天躁狠狠躁夜夜2o2o| 欧美精品国产亚洲| 婷婷六月久久综合丁香| 国产精品女同一区二区软件| 欧美成人一区二区免费高清观看| 亚洲va在线va天堂va国产| 在现免费观看毛片| 一级黄片播放器| 中出人妻视频一区二区| 欧美日本视频| 日本五十路高清| 国产精品久久久久久久电影| 亚洲精华国产精华液的使用体验 | 大香蕉久久网| 又粗又爽又猛毛片免费看| 婷婷精品国产亚洲av在线| 亚洲专区国产一区二区| 欧美bdsm另类| 综合色丁香网| 国产精品久久久久久久久免| 如何舔出高潮| 亚洲在线观看片| 免费看av在线观看网站| 俺也久久电影网| 国产精品女同一区二区软件| 色吧在线观看| 一级a爱片免费观看的视频| 亚洲国产精品成人综合色| 中出人妻视频一区二区| 亚洲欧美日韩东京热| 久久久久久久久中文| 内地一区二区视频在线| 午夜激情福利司机影院| av专区在线播放| 国产老妇女一区| 亚洲精品一区av在线观看| 热99在线观看视频| or卡值多少钱| 韩国av在线不卡| 美女被艹到高潮喷水动态| 亚洲中文日韩欧美视频| 婷婷亚洲欧美| 久久欧美精品欧美久久欧美| 亚洲欧美精品自产自拍| 看十八女毛片水多多多| 亚洲无线观看免费| 女同久久另类99精品国产91| 黑人高潮一二区| 少妇人妻精品综合一区二区 | .国产精品久久| 色5月婷婷丁香| 99精品在免费线老司机午夜| 欧美国产日韩亚洲一区| 麻豆av噜噜一区二区三区| 国产精品嫩草影院av在线观看| 一个人观看的视频www高清免费观看| 国产精品三级大全| 乱系列少妇在线播放| 性欧美人与动物交配| 能在线免费观看的黄片| 男人的好看免费观看在线视频| 亚洲av免费高清在线观看| 亚洲av成人精品一区久久| 久久久久国产精品人妻aⅴ院| 精品人妻偷拍中文字幕| 免费在线观看成人毛片| 99riav亚洲国产免费| 欧美激情在线99| 18禁在线播放成人免费| 哪里可以看免费的av片| 一进一出好大好爽视频| 亚洲国产精品sss在线观看| 日韩精品青青久久久久久| 国产午夜福利久久久久久| 老熟妇乱子伦视频在线观看| 中国美白少妇内射xxxbb| 一级毛片aaaaaa免费看小| 国产乱人偷精品视频| 人妻丰满熟妇av一区二区三区| 久久久久久九九精品二区国产| 久久久久久国产a免费观看| 中文字幕免费在线视频6| 又粗又爽又猛毛片免费看| 嫩草影院入口| 久久久久久久久久久丰满| 女同久久另类99精品国产91| 国产av不卡久久| 亚洲欧美日韩高清专用| 97超碰精品成人国产| 成人一区二区视频在线观看| 国产 一区 欧美 日韩| 99久久成人亚洲精品观看| 精品99又大又爽又粗少妇毛片| 国产视频内射| 国产中年淑女户外野战色| 99久久精品一区二区三区| 久久亚洲精品不卡| 国产淫片久久久久久久久| 国产精品亚洲美女久久久| 不卡视频在线观看欧美| 91在线观看av| 国产精品不卡视频一区二区| 日韩欧美在线乱码| 午夜免费男女啪啪视频观看 | 欧美日本亚洲视频在线播放| 久久欧美精品欧美久久欧美| 欧美+亚洲+日韩+国产| 一区二区三区高清视频在线| 人妻夜夜爽99麻豆av| 亚洲成人久久性| 国产女主播在线喷水免费视频网站 | 校园人妻丝袜中文字幕| 国产成人aa在线观看| 国产免费男女视频| 国产高清视频在线观看网站| 亚洲va在线va天堂va国产| 成人欧美大片| 亚洲真实伦在线观看| 国产精品福利在线免费观看| 日韩欧美精品v在线| 一级a爱片免费观看的视频| 亚洲精品成人久久久久久| www日本黄色视频网| 久久精品久久久久久噜噜老黄 | 一级毛片我不卡| 久久九九热精品免费| 老司机影院成人| 亚洲国产精品国产精品| 精品一区二区三区视频在线观看免费| 日本黄色视频三级网站网址| 成年女人看的毛片在线观看| 国内精品美女久久久久久| 22中文网久久字幕| 免费看a级黄色片| 日韩中字成人| 亚洲无线在线观看| 麻豆乱淫一区二区| 亚洲中文日韩欧美视频| 欧美日韩在线观看h| 午夜福利在线观看吧| 中国美白少妇内射xxxbb| 精品国产三级普通话版| 97超视频在线观看视频| 亚洲精品粉嫩美女一区| 欧美性猛交黑人性爽| 午夜亚洲福利在线播放| av卡一久久| 99久国产av精品国产电影| av女优亚洲男人天堂| 精品无人区乱码1区二区| 蜜桃亚洲精品一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 亚洲美女搞黄在线观看 | 日韩强制内射视频| 天堂影院成人在线观看| 日日啪夜夜撸| 一区福利在线观看| 免费不卡的大黄色大毛片视频在线观看 | 别揉我奶头 嗯啊视频| 99热这里只有是精品50| a级一级毛片免费在线观看| 国产精品亚洲美女久久久| 亚洲经典国产精华液单| 1000部很黄的大片| 少妇人妻一区二区三区视频| 一进一出好大好爽视频| 男女边吃奶边做爰视频| 日本爱情动作片www.在线观看 | 不卡视频在线观看欧美| .国产精品久久| 国产乱人视频| 国产一区二区在线观看日韩| 免费不卡的大黄色大毛片视频在线观看 | 国内久久婷婷六月综合欲色啪| 亚洲三级黄色毛片| 亚洲中文字幕一区二区三区有码在线看| 亚洲熟妇中文字幕五十中出| 久久久国产成人免费| 日韩一区二区视频免费看| 男女啪啪激烈高潮av片| a级毛片a级免费在线| 国产亚洲欧美98| 不卡视频在线观看欧美| 午夜福利在线在线| 国产精品,欧美在线| 一级毛片电影观看 | 一级毛片我不卡| 国产日本99.免费观看| 国内久久婷婷六月综合欲色啪| 午夜视频国产福利| 一级黄色大片毛片| 少妇的逼好多水| 欧洲精品卡2卡3卡4卡5卡区| 全区人妻精品视频| 亚洲性夜色夜夜综合| 午夜精品国产一区二区电影 | 亚洲av中文av极速乱| 一区二区三区免费毛片| 午夜福利成人在线免费观看| 亚洲国产日韩欧美精品在线观看| 国产精品免费一区二区三区在线| 一本一本综合久久| 国产大屁股一区二区在线视频| 国产精品亚洲美女久久久| 婷婷色综合大香蕉| 午夜爱爱视频在线播放| 欧美激情国产日韩精品一区| 久久久成人免费电影| 男插女下体视频免费在线播放| 亚洲av成人精品一区久久| 12—13女人毛片做爰片一| 中文字幕熟女人妻在线| 成人鲁丝片一二三区免费| 国产精品国产高清国产av| 亚洲乱码一区二区免费版| 天堂√8在线中文| 亚洲va在线va天堂va国产| 人妻久久中文字幕网| 亚洲中文日韩欧美视频| 国产极品精品免费视频能看的| 欧美bdsm另类| 久久综合国产亚洲精品| 亚洲成人久久爱视频| 久久久久久久亚洲中文字幕| 网址你懂的国产日韩在线| 日本与韩国留学比较| 精品久久久久久久久av| 国产亚洲av嫩草精品影院| 高清毛片免费看| 国产成人freesex在线 | av视频在线观看入口| 久久久久性生活片| 看片在线看免费视频| 高清日韩中文字幕在线| 国产久久久一区二区三区| 精品人妻视频免费看| 午夜激情福利司机影院| 老熟妇仑乱视频hdxx| 午夜福利视频1000在线观看| 久久精品国产亚洲av天美| 国产探花极品一区二区| 亚洲欧美日韩无卡精品| 深夜精品福利| 性色avwww在线观看| 亚洲国产精品久久男人天堂| 天天一区二区日本电影三级| 神马国产精品三级电影在线观看| 亚洲国产欧洲综合997久久,| 午夜老司机福利剧场| 免费观看精品视频网站| 99热网站在线观看| 中文字幕免费在线视频6| 深夜a级毛片| 久久99热6这里只有精品| av在线亚洲专区| 国产男靠女视频免费网站| 色综合站精品国产| 国产真实乱freesex| 国产成人a区在线观看| 日本黄色视频三级网站网址| 我的女老师完整版在线观看| 一区二区三区免费毛片| 亚洲第一电影网av| 亚洲av二区三区四区| 久久精品国产自在天天线| 日本爱情动作片www.在线观看 | 欧美日韩一区二区视频在线观看视频在线 | 成人综合一区亚洲| 日本黄大片高清| 22中文网久久字幕| 久久久久久国产a免费观看| 国产精品美女特级片免费视频播放器| 亚洲三级黄色毛片| 成人毛片a级毛片在线播放| 天堂√8在线中文| 久99久视频精品免费| 成人欧美大片| 精品人妻偷拍中文字幕| 国产精品爽爽va在线观看网站| 美女 人体艺术 gogo| 亚洲人成网站在线观看播放| 女生性感内裤真人,穿戴方法视频| 啦啦啦观看免费观看视频高清| 搡老岳熟女国产| 99国产精品一区二区蜜桃av| 蜜桃久久精品国产亚洲av| 亚洲av二区三区四区| 国产精品免费一区二区三区在线| 成人特级av手机在线观看| 国内精品宾馆在线| 免费无遮挡裸体视频| 哪里可以看免费的av片| 国产91av在线免费观看| 18禁黄网站禁片免费观看直播| 国产v大片淫在线免费观看| 国产午夜精品久久久久久一区二区三区 | 免费人成视频x8x8入口观看| 国产白丝娇喘喷水9色精品| 亚洲av.av天堂| 国产淫片久久久久久久久| 在线国产一区二区在线| 少妇人妻一区二区三区视频| 免费人成视频x8x8入口观看| 亚洲天堂国产精品一区在线| 色哟哟哟哟哟哟| 日韩强制内射视频| 欧美成人精品欧美一级黄| 欧美日韩综合久久久久久| 国产欧美日韩精品亚洲av| 99热全是精品| 亚洲国产精品合色在线| 久久中文看片网| 国产美女午夜福利| 91狼人影院| 国产伦精品一区二区三区视频9| 久久国内精品自在自线图片| 欧美潮喷喷水| 精品久久久久久久久亚洲| 高清毛片免费观看视频网站| 亚洲三级黄色毛片| 十八禁网站免费在线| 日韩欧美国产在线观看| 国产精品免费一区二区三区在线| 真实男女啪啪啪动态图| 欧美又色又爽又黄视频| 免费人成视频x8x8入口观看| 久久久精品大字幕| 最好的美女福利视频网| 18禁在线无遮挡免费观看视频 | 激情 狠狠 欧美| 欧美激情国产日韩精品一区| 深夜精品福利| 免费看av在线观看网站| 国产熟女欧美一区二区| 嫩草影院精品99| 99久久成人亚洲精品观看| 91午夜精品亚洲一区二区三区| 一本久久中文字幕| 国产极品精品免费视频能看的| 可以在线观看的亚洲视频| 婷婷亚洲欧美| 女生性感内裤真人,穿戴方法视频| 成人一区二区视频在线观看| 久久久精品94久久精品| 内射极品少妇av片p| 亚洲自偷自拍三级| 级片在线观看| 欧美精品国产亚洲| 亚洲av中文字字幕乱码综合| 91久久精品国产一区二区三区| 久久久久久久久久久丰满| 欧美又色又爽又黄视频| 午夜老司机福利剧场| 日韩成人伦理影院| 超碰av人人做人人爽久久| 伊人久久精品亚洲午夜| 国国产精品蜜臀av免费| 色尼玛亚洲综合影院| 人人妻,人人澡人人爽秒播| 97热精品久久久久久| 97超碰精品成人国产| 97超视频在线观看视频| 91精品国产九色| 国产乱人视频| 亚洲va在线va天堂va国产| 特级一级黄色大片| 我要看日韩黄色一级片| 国产免费男女视频| 卡戴珊不雅视频在线播放| 一级黄色大片毛片| 日本免费a在线| a级一级毛片免费在线观看| 丰满的人妻完整版| 国产单亲对白刺激| 在线观看av片永久免费下载| 成人无遮挡网站| 日韩国内少妇激情av| 久久久久久国产a免费观看| 午夜久久久久精精品| av国产免费在线观看| 国产黄色视频一区二区在线观看 | 国产黄a三级三级三级人| 国产精品野战在线观看| av在线天堂中文字幕| 少妇熟女aⅴ在线视频| 性插视频无遮挡在线免费观看| 99视频精品全部免费 在线| 男人的好看免费观看在线视频| 夜夜看夜夜爽夜夜摸| 给我免费播放毛片高清在线观看| 亚洲电影在线观看av| 国内精品宾馆在线| 我的老师免费观看完整版| 亚洲国产精品成人久久小说 | 亚洲精品亚洲一区二区| 男人狂女人下面高潮的视频| 国产精品久久久久久精品电影| 午夜久久久久精精品| 男女做爰动态图高潮gif福利片| 晚上一个人看的免费电影| 日本与韩国留学比较| 国产中年淑女户外野战色| 国产女主播在线喷水免费视频网站 | 两性午夜刺激爽爽歪歪视频在线观看| 淫秽高清视频在线观看| 99久久精品一区二区三区|