• <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)化模型
    一级毛片 在线播放| 亚洲怡红院男人天堂| 日本黄色日本黄色录像| 国产精品99久久99久久久不卡 | 免费人妻精品一区二区三区视频| 夫妻性生交免费视频一级片| 国产熟女午夜一区二区三区 | 下体分泌物呈黄色| 国产精品麻豆人妻色哟哟久久| 黄色配什么色好看| 哪个播放器可以免费观看大片| a级毛片免费高清观看在线播放| 国产亚洲精品久久久com| 久久精品国产a三级三级三级| 伊人久久精品亚洲午夜| 街头女战士在线观看网站| 午夜久久久在线观看| 一区二区三区乱码不卡18| 内地一区二区视频在线| 久久久国产精品麻豆| 观看av在线不卡| 18禁在线播放成人免费| 男人狂女人下面高潮的视频| 国产毛片在线视频| 美女主播在线视频| 嫩草影院入口| 国产在线一区二区三区精| 国产精品一区二区三区四区免费观看| 一区二区av电影网| 亚洲国产最新在线播放| 又黄又爽又刺激的免费视频.| 国产精品一区二区在线观看99| 美女主播在线视频| 王馨瑶露胸无遮挡在线观看| 亚洲不卡免费看| 99re6热这里在线精品视频| 永久网站在线| 日本av免费视频播放| 夫妻午夜视频| 亚洲国产精品一区三区| 欧美精品亚洲一区二区| 国产成人午夜福利电影在线观看| 国产精品偷伦视频观看了| 搡老乐熟女国产| 久久精品夜色国产| 国产伦在线观看视频一区| 久久99精品国语久久久| 综合色丁香网| 久久久久久人妻| 黄色配什么色好看| 桃花免费在线播放| 国产老妇伦熟女老妇高清| tube8黄色片| 在线天堂最新版资源| 美女国产视频在线观看| 人妻人人澡人人爽人人| 一个人免费看片子| 男的添女的下面高潮视频| 大片电影免费在线观看免费| 久久精品国产亚洲网站| 国产成人午夜福利电影在线观看| 国产精品国产三级专区第一集| 九九爱精品视频在线观看| 女性生殖器流出的白浆| 亚洲怡红院男人天堂| 国产免费一区二区三区四区乱码| 国产成人精品一,二区| 国产免费福利视频在线观看| 激情五月婷婷亚洲| 国产一级毛片在线| 亚洲国产日韩一区二区| 一区二区av电影网| 亚洲精品日本国产第一区| 日韩精品免费视频一区二区三区 | 午夜免费观看性视频| 22中文网久久字幕| 亚洲av电影在线观看一区二区三区| 国产中年淑女户外野战色| 免费看不卡的av| 日产精品乱码卡一卡2卡三| 9色porny在线观看| 精品99又大又爽又粗少妇毛片| 大香蕉久久网| 国产欧美另类精品又又久久亚洲欧美| 欧美日韩av久久| 午夜影院在线不卡| 91午夜精品亚洲一区二区三区| 男女国产视频网站| 国产精品蜜桃在线观看| 人体艺术视频欧美日本| 国产精品三级大全| 日本黄色片子视频| 九草在线视频观看| 国产成人免费观看mmmm| 熟女人妻精品中文字幕| 国产探花极品一区二区| 色吧在线观看| 狂野欧美激情性bbbbbb| 黄色一级大片看看| 在线观看三级黄色| 国产69精品久久久久777片| 99国产精品免费福利视频| 五月伊人婷婷丁香| 亚洲精品国产av成人精品| 久久国产精品男人的天堂亚洲 | 你懂的网址亚洲精品在线观看| 爱豆传媒免费全集在线观看| 三级国产精品片| 久久久久久久久久成人| 日韩不卡一区二区三区视频在线| 麻豆成人午夜福利视频| 欧美日本中文国产一区发布| 色吧在线观看| 免费在线观看成人毛片| 久久久久久久久久久免费av| 2018国产大陆天天弄谢| 超碰97精品在线观看| 免费大片黄手机在线观看| 又大又黄又爽视频免费| 五月天丁香电影| 秋霞伦理黄片| 日韩伦理黄色片| 午夜福利影视在线免费观看| 久久亚洲国产成人精品v| 日本色播在线视频| 热99国产精品久久久久久7| 两个人免费观看高清视频 | 乱码一卡2卡4卡精品| 高清黄色对白视频在线免费看 | 日韩熟女老妇一区二区性免费视频| 日韩制服骚丝袜av| 最近2019中文字幕mv第一页| 免费看光身美女| 亚洲国产av新网站| 国产精品久久久久久精品电影小说| 少妇被粗大的猛进出69影院 | 日韩,欧美,国产一区二区三区| 亚洲av成人精品一二三区| 国产精品久久久久成人av| 久久婷婷青草| 国产欧美日韩综合在线一区二区 | 黄色一级大片看看| 亚洲精品色激情综合| 最近的中文字幕免费完整| 伊人久久精品亚洲午夜| 女人精品久久久久毛片| 久久99热6这里只有精品| 国产精品一区二区三区四区免费观看| 有码 亚洲区| 天堂俺去俺来也www色官网| 性高湖久久久久久久久免费观看| 国产亚洲91精品色在线| 国产精品久久久久久精品古装| 成人亚洲精品一区在线观看| 成人特级av手机在线观看| 免费看日本二区| 国产av精品麻豆| 黑人猛操日本美女一级片| 美女视频免费永久观看网站| 欧美精品一区二区大全| 啦啦啦视频在线资源免费观看| 国产色爽女视频免费观看| 男人舔奶头视频| 91在线精品国自产拍蜜月| 一个人免费看片子| 亚洲精品日韩av片在线观看| 丝袜脚勾引网站| 99re6热这里在线精品视频| 国产成人免费观看mmmm| 婷婷色av中文字幕| 国产成人精品无人区| 9色porny在线观看| 成人综合一区亚洲| 国产欧美日韩综合在线一区二区 | 久久久久久久久久人人人人人人| 国产亚洲午夜精品一区二区久久| 久久久久久久大尺度免费视频| 国产av码专区亚洲av| 97在线视频观看| 91精品国产九色| 中文字幕人妻丝袜制服| 青青草视频在线视频观看| 亚洲精品一区蜜桃| 欧美精品一区二区免费开放| 在线 av 中文字幕| 人人妻人人澡人人爽人人夜夜| 高清av免费在线| 亚洲国产欧美日韩在线播放 | 妹子高潮喷水视频| 色网站视频免费| 多毛熟女@视频| 晚上一个人看的免费电影| 日本免费在线观看一区| av国产精品久久久久影院| 国产色爽女视频免费观看| 自拍偷自拍亚洲精品老妇| 99九九在线精品视频 | 久久 成人 亚洲| 日韩视频在线欧美| 人人妻人人澡人人看| 丝袜脚勾引网站| 大香蕉久久网| 偷拍熟女少妇极品色| 亚洲,一卡二卡三卡| 久久久久久久亚洲中文字幕| 26uuu在线亚洲综合色| 丝袜脚勾引网站| 三级经典国产精品| 视频中文字幕在线观看| 久久精品国产亚洲网站| av专区在线播放| 黄色视频在线播放观看不卡| 久久精品熟女亚洲av麻豆精品| 亚洲欧美精品自产自拍| 九九久久精品国产亚洲av麻豆| 国产亚洲最大av| 亚洲国产色片| 免费黄频网站在线观看国产| 日韩熟女老妇一区二区性免费视频| 日日啪夜夜撸| 日本猛色少妇xxxxx猛交久久| 亚洲国产最新在线播放| 久久久久久久国产电影| 内射极品少妇av片p| 中文欧美无线码| 久久久久人妻精品一区果冻| 男人爽女人下面视频在线观看| 精品人妻一区二区三区麻豆| 波野结衣二区三区在线| 国产真实伦视频高清在线观看| 多毛熟女@视频| 欧美 亚洲 国产 日韩一| 久久精品国产亚洲av天美| 少妇高潮的动态图| 国内揄拍国产精品人妻在线| 精品少妇内射三级| 亚洲精品乱码久久久v下载方式| 亚洲精品aⅴ在线观看| 久久韩国三级中文字幕| 亚洲,欧美,日韩| 亚洲av福利一区| 亚洲精品国产av蜜桃| 亚洲精品乱码久久久久久按摩| 美女内射精品一级片tv| 国产一区二区三区综合在线观看 | 免费观看在线日韩| 午夜福利视频精品| 伊人久久精品亚洲午夜| 在线观看国产h片| 国产亚洲av片在线观看秒播厂| av在线观看视频网站免费| 国产亚洲精品久久久com| 看免费成人av毛片| 国产高清国产精品国产三级| 久久精品熟女亚洲av麻豆精品| 少妇人妻 视频| 麻豆乱淫一区二区| 男女免费视频国产| 中国国产av一级| 好男人视频免费观看在线| 色哟哟·www| 性色avwww在线观看| 三上悠亚av全集在线观看 | 99久久精品热视频| 六月丁香七月| 热re99久久精品国产66热6| 国产又色又爽无遮挡免| 成人毛片a级毛片在线播放| 少妇的逼好多水| 久久精品国产亚洲网站| 久热这里只有精品99| 亚洲精品自拍成人| 亚洲国产精品一区三区| 麻豆精品久久久久久蜜桃| 丝袜喷水一区| 久久国产精品大桥未久av | 亚洲精品中文字幕在线视频 | 只有这里有精品99| 噜噜噜噜噜久久久久久91| 在线观看www视频免费| 看免费成人av毛片| 七月丁香在线播放| 亚洲人与动物交配视频| 内射极品少妇av片p| 日本vs欧美在线观看视频 | 国产成人精品久久久久久| 国产精品三级大全| 丰满饥渴人妻一区二区三| 欧美激情极品国产一区二区三区 | 国产一区二区三区综合在线观看 | 五月开心婷婷网| 97超视频在线观看视频| 亚洲av综合色区一区| 日韩av不卡免费在线播放| 插阴视频在线观看视频| 色网站视频免费| 国精品久久久久久国模美| 在线播放无遮挡| 亚洲av福利一区| 中文字幕人妻熟人妻熟丝袜美| 欧美bdsm另类| 男人狂女人下面高潮的视频| 午夜av观看不卡| 韩国av在线不卡| 交换朋友夫妻互换小说| 性色avwww在线观看| 五月开心婷婷网| 校园人妻丝袜中文字幕| 麻豆成人av视频| 亚洲精品国产色婷婷电影| 精品视频人人做人人爽| 日韩av免费高清视频| videos熟女内射| 最新中文字幕久久久久| 免费看日本二区| 亚洲精品,欧美精品| 亚洲四区av| 欧美人与善性xxx| 久久久久视频综合| √禁漫天堂资源中文www| 在线观看免费高清a一片| 久久精品久久久久久久性| 日本-黄色视频高清免费观看| 啦啦啦中文免费视频观看日本| 国产色婷婷99| 美女主播在线视频| 乱系列少妇在线播放| 亚洲成人av在线免费| 在线观看三级黄色| 晚上一个人看的免费电影| 亚洲国产精品一区二区三区在线| 久久青草综合色| av专区在线播放| 九九久久精品国产亚洲av麻豆| 在线观看www视频免费| 精品熟女少妇av免费看| 国产高清国产精品国产三级| 一级爰片在线观看| 最新中文字幕久久久久| 国产视频内射| 纵有疾风起免费观看全集完整版| 毛片一级片免费看久久久久| 日韩一区二区视频免费看| 久久久久久久久大av| 亚洲av综合色区一区| 亚洲中文av在线| 国产亚洲午夜精品一区二区久久| 自线自在国产av| 只有这里有精品99| 日本欧美国产在线视频| 看十八女毛片水多多多| 国产亚洲午夜精品一区二区久久| 精品人妻一区二区三区麻豆| 只有这里有精品99| 国产精品熟女久久久久浪| 精品人妻熟女毛片av久久网站| 久久人人爽人人爽人人片va| 久久久久久人妻| av网站免费在线观看视频| 精品少妇内射三级| 精品人妻熟女毛片av久久网站| 国产男女内射视频| 狠狠精品人妻久久久久久综合| 黄色日韩在线| 久久午夜综合久久蜜桃| 新久久久久国产一级毛片| 国产精品国产三级专区第一集| 亚洲四区av| 岛国毛片在线播放| 日韩av免费高清视频| 纯流量卡能插随身wifi吗| 永久免费av网站大全| 国精品久久久久久国模美| 国产欧美日韩精品一区二区| 成人国产麻豆网| 久久精品夜色国产| 国产av码专区亚洲av| 天堂8中文在线网| 欧美一级a爱片免费观看看| 又大又黄又爽视频免费| 亚洲精品第二区| 在线天堂最新版资源| 日本爱情动作片www.在线观看| 国产色婷婷99| 久久久久精品久久久久真实原创| 欧美精品亚洲一区二区| 国产精品一区二区三区四区免费观看| 自拍欧美九色日韩亚洲蝌蚪91 | 五月开心婷婷网| 国内少妇人妻偷人精品xxx网站| 国产成人免费无遮挡视频| 美女中出高潮动态图| 久久99一区二区三区| 久久久久久久精品精品| a级一级毛片免费在线观看| 日本wwww免费看| 久久精品国产亚洲网站| 国产熟女欧美一区二区| 免费大片18禁| 男女无遮挡免费网站观看| 又大又黄又爽视频免费| 免费在线观看成人毛片| 大话2 男鬼变身卡| 毛片一级片免费看久久久久| 久久亚洲国产成人精品v| 亚洲自偷自拍三级| 国产爽快片一区二区三区| 午夜日本视频在线| 国产亚洲精品久久久com| 三级经典国产精品| 成人影院久久| 精品一区二区三卡| 久久久久久伊人网av| videos熟女内射| av又黄又爽大尺度在线免费看| 日韩视频在线欧美| 国产精品一区二区在线不卡| 国产在线视频一区二区| 丰满少妇做爰视频| 高清黄色对白视频在线免费看 | 少妇被粗大猛烈的视频| 亚洲av不卡在线观看| 菩萨蛮人人尽说江南好唐韦庄| 99热这里只有精品一区| 日韩人妻高清精品专区| 国产成人免费观看mmmm| 丰满乱子伦码专区| 老熟女久久久| 久久99蜜桃精品久久| 欧美 亚洲 国产 日韩一| 国产精品麻豆人妻色哟哟久久| 如日韩欧美国产精品一区二区三区 | 午夜影院在线不卡| 九九久久精品国产亚洲av麻豆| 国产精品伦人一区二区| 男人舔奶头视频| 久久久久网色| 欧美区成人在线视频| 91久久精品国产一区二区成人| 欧美+日韩+精品| 成人特级av手机在线观看| 成人综合一区亚洲| 日日撸夜夜添| 亚洲不卡免费看| 成人黄色视频免费在线看| 午夜激情福利司机影院| 中文天堂在线官网| 99久国产av精品国产电影| 欧美日韩综合久久久久久| 日日摸夜夜添夜夜爱| 超碰97精品在线观看| 亚洲熟女精品中文字幕| 国产在线一区二区三区精| 国产男人的电影天堂91| 午夜福利,免费看| 卡戴珊不雅视频在线播放| av.在线天堂| 亚洲欧美一区二区三区国产| 观看av在线不卡| 丁香六月天网| 国产视频首页在线观看| 国产在线免费精品| 老司机亚洲免费影院| 蜜臀久久99精品久久宅男| 在线亚洲精品国产二区图片欧美 | 99久久精品国产国产毛片| 日韩av不卡免费在线播放| 国产一级毛片在线| 18禁裸乳无遮挡动漫免费视频| 国产av一区二区精品久久| 久久久久久久精品精品| 欧美日韩一区二区视频在线观看视频在线| 免费黄色在线免费观看| 国内揄拍国产精品人妻在线| 大又大粗又爽又黄少妇毛片口| 中国国产av一级| 成人亚洲精品一区在线观看| 亚洲欧美中文字幕日韩二区| 日本av免费视频播放| 又爽又黄a免费视频| 国产亚洲午夜精品一区二区久久| 人妻制服诱惑在线中文字幕| 春色校园在线视频观看| 亚洲真实伦在线观看| 黑人巨大精品欧美一区二区蜜桃 | 女人精品久久久久毛片| 国产成人午夜福利电影在线观看| 美女福利国产在线| 国产永久视频网站| 特大巨黑吊av在线直播| 777米奇影视久久| 国产成人一区二区在线| 国语对白做爰xxxⅹ性视频网站| 国内少妇人妻偷人精品xxx网站| 久热久热在线精品观看| 五月玫瑰六月丁香| 国产精品99久久久久久久久| 中文字幕精品免费在线观看视频 | 亚洲国产精品一区二区三区在线| 妹子高潮喷水视频| 97超视频在线观看视频| 国产伦在线观看视频一区| 日韩欧美精品免费久久| 亚洲激情五月婷婷啪啪| 久久精品国产亚洲av天美| 精品久久久久久久久av| 国产成人精品福利久久| 简卡轻食公司| 日韩精品有码人妻一区| 精品一品国产午夜福利视频| 黄色毛片三级朝国网站 | 最近中文字幕2019免费版| 最近2019中文字幕mv第一页| 亚洲av成人精品一二三区| 精品一品国产午夜福利视频| 一二三四中文在线观看免费高清| 两个人的视频大全免费| 中文字幕制服av| 高清av免费在线| 久久鲁丝午夜福利片| 日韩亚洲欧美综合| 午夜av观看不卡| 久久久久久久亚洲中文字幕| 色视频在线一区二区三区| 九九在线视频观看精品| 免费在线观看成人毛片| 成人国产av品久久久| 肉色欧美久久久久久久蜜桃| 国产视频内射| 免费不卡的大黄色大毛片视频在线观看| 熟女av电影| 久久午夜福利片| 亚洲精品久久久久久婷婷小说| 亚洲国产精品国产精品| av天堂久久9| 插逼视频在线观看| 另类精品久久| 久久久久视频综合| 内射极品少妇av片p| 热99国产精品久久久久久7| 国产真实伦视频高清在线观看| 国产日韩欧美亚洲二区| 免费大片18禁| 在线看a的网站| 中文字幕制服av| 夫妻午夜视频| 日韩熟女老妇一区二区性免费视频| 国产成人aa在线观看| 少妇人妻 视频| 人人妻人人看人人澡| 一级毛片我不卡| 街头女战士在线观看网站| 久久久久久久亚洲中文字幕| 免费观看性生交大片5| 国产女主播在线喷水免费视频网站| 在线亚洲精品国产二区图片欧美 | 午夜激情福利司机影院| 国产成人免费观看mmmm| 亚洲久久久国产精品| 在现免费观看毛片| 最新中文字幕久久久久| 国产成人午夜福利电影在线观看| 亚洲综合精品二区| 日本色播在线视频| 99久久精品一区二区三区| 18禁在线无遮挡免费观看视频| 亚洲一区二区三区欧美精品| 国产乱人偷精品视频| 一级二级三级毛片免费看| 制服丝袜香蕉在线| 中文字幕亚洲精品专区| 精品一区在线观看国产| 你懂的网址亚洲精品在线观看| 男人狂女人下面高潮的视频| 国产国拍精品亚洲av在线观看| h日本视频在线播放| 青春草国产在线视频| av又黄又爽大尺度在线免费看| 岛国毛片在线播放| 97精品久久久久久久久久精品| 一级毛片 在线播放| 亚洲在久久综合| 丝袜脚勾引网站| 亚洲av成人精品一二三区| 欧美老熟妇乱子伦牲交| 五月天丁香电影| 国产精品秋霞免费鲁丝片| 色5月婷婷丁香| 精品一区二区三区视频在线| 国产伦精品一区二区三区视频9| 久久av网站| av女优亚洲男人天堂| 菩萨蛮人人尽说江南好唐韦庄| 高清视频免费观看一区二区| 亚洲国产色片| 美女国产视频在线观看| 精品少妇黑人巨大在线播放| 精品人妻熟女av久视频| 我的女老师完整版在线观看| 99国产精品免费福利视频| 熟女人妻精品中文字幕| 成人国产av品久久久| 如日韩欧美国产精品一区二区三区 | 蜜桃久久精品国产亚洲av| 国产高清国产精品国产三级| 国产免费福利视频在线观看| 国产深夜福利视频在线观看| 人人澡人人妻人| 纵有疾风起免费观看全集完整版| 精品亚洲乱码少妇综合久久| 综合色丁香网| 男女啪啪激烈高潮av片| 99热这里只有是精品50| 国产精品久久久久久av不卡| 黄片无遮挡物在线观看|