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

    乙烯燃料超燃沖壓發(fā)動機燃燒過程研究

    2021-03-20 03:37:26鐘富宇樂嘉陵岳茂雄
    實驗流體力學(xué) 2021年1期
    關(guān)鍵詞:凹腔激波燃燒室

    鐘富宇, 樂嘉陵, 田 野, 岳茂雄

    中國空氣動力研究與發(fā)展中心 高超聲速沖壓發(fā)動機技術(shù)重點實驗室, 四川 綿陽 621000

    0 引 言

    吸氣式高超聲速飛行器是實現(xiàn)全球范圍內(nèi)快速遠(yuǎn)程民用運輸、可突防遠(yuǎn)程快速打擊以及近地軌道空間運輸?shù)鹊闹匾d體[1]。超燃沖壓發(fā)動機是這類飛行器最具推廣意義的推進裝置,尤其是在馬赫數(shù)大于5.0時具有非常廣闊的應(yīng)用前景。超燃沖壓發(fā)動機燃燒室的可靠點火和穩(wěn)定燃燒過程一直是人們關(guān)注的焦點。在超聲速燃燒過程中,對流場結(jié)構(gòu)和火焰發(fā)展變化過程的研究能夠加深對點火和燃燒機理的認(rèn)識,需要采用多種非接觸式燃燒流場診斷技術(shù)對其進行觀測,以達(dá)到多種測量結(jié)果相互補充和印證的目的。

    田野等[2]通過試驗研究了不同當(dāng)量比的氫燃料對燃燒模態(tài)和燃燒流場結(jié)構(gòu)的影響,采用了多種非接觸光學(xué)測量方法進行同步測量,結(jié)果表明:當(dāng)量比大于0.17時,燃燒流場結(jié)構(gòu)不穩(wěn)定,火焰分布呈現(xiàn)破碎狀,火焰在燃燒室上下壁面之間來回傳播;當(dāng)量比小于或等于0.17時,燃燒流場結(jié)構(gòu)穩(wěn)定,火焰呈現(xiàn)連續(xù)分布,火焰穩(wěn)定分布于凹槽下部剪切層內(nèi)。吳戈等[3]綜合利用OH-PLIF及CH-PLIF研究了乙烯燃料在凹腔流動的幾個截面火焰結(jié)構(gòu),結(jié)果表明:當(dāng)量比較低時,燃燒的火焰反應(yīng)區(qū)主要在凹腔中,OH沿中心軸對稱分布,當(dāng)量比較高時OH主要沿燃燒室兩側(cè)壁分布;CH-PLIF的結(jié)果顯示放熱區(qū)呈現(xiàn)高度褶皺的破碎狀,放熱區(qū)分布范圍比反應(yīng)區(qū)更窄。Ruan等[4]在基于凹腔的超燃沖壓發(fā)動機燃燒室構(gòu)型中進行了反應(yīng)流和無反應(yīng)流的試驗,結(jié)果表明:燃燒首先在擴散火焰控制下進行,然后是充分混合后的高效率燃燒,在凹腔內(nèi)也出現(xiàn)了明顯的擴散火焰控制燃燒現(xiàn)象。Gordon和Mastorakos[5-7]研究了不同來流條件下的點火和火焰穩(wěn)定,發(fā)現(xiàn)火焰穩(wěn)定性對來流溫度和速度很敏感,高總焓的來流條件有利于燃料的自點火和火焰穩(wěn)定。張灣洲等[8]采用試驗方法研究了飛行馬赫數(shù)Ma=4.0、總溫935 K來流參數(shù)下的超燃發(fā)動機乙烯點火試驗,結(jié)果表明:回流區(qū)有利于點火,剪切層和凹槽后部是穩(wěn)焰的主要區(qū)域;點火成功后,剪切層內(nèi)和凹槽后部持續(xù)卷吸氧化劑,能夠維持穩(wěn)定的燃燒。Brieschenk等[9-10]對超燃沖壓發(fā)動機激光誘導(dǎo)等離子體(LIP)點火行為進行了試驗研究,采用OH-PLIF表征等離子點火區(qū)域的演變過程,結(jié)果表明高超聲速流動中等離子點火可以促進羥基的形成。Kumaran等[11]用數(shù)值方法研究了化學(xué)模型對超聲速燃燒的影響,計算結(jié)果表明,多步化學(xué)反應(yīng)預(yù)測比單步化學(xué)反應(yīng)預(yù)測的釋熱量更高、釋熱區(qū)域更大,多步化學(xué)預(yù)測的燃燒過程更詳細(xì)、更復(fù)雜。Nakaya等[12]采用試驗方法研究了正十二烷中各組分對點火和燃燒的影響,發(fā)現(xiàn)在所有的超聲速燃燒情況下,點火均是在擴張截面處的邊界層中發(fā)生,火焰向上游傳播最終穩(wěn)定在凹腔內(nèi),基于壓力和OH基的熒光發(fā)現(xiàn)了2種燃燒模式:射流尾跡穩(wěn)定模式和凹腔火焰穩(wěn)定模式。

    何粲等[13]對一種乙烯燃料矩形截面超燃沖壓發(fā)動機進行了數(shù)值模擬,研究了該模型在不同燃燒模態(tài)下的流動特性,與試驗結(jié)果吻合良好,其結(jié)果表明:發(fā)動機處于雙模態(tài)超燃或雙模態(tài)亞燃模態(tài)時,隨著激波串結(jié)構(gòu)的形成與前移,部分燃燒可能在隔離段內(nèi)完成;而對于純超燃模態(tài),燃燒僅發(fā)生在凹槽與擴張段內(nèi),化學(xué)反應(yīng)與高溫區(qū)的分布相對更集中。Qin等[14]研究了沖壓發(fā)動機燃燒室內(nèi)的火焰穩(wěn)定性,結(jié)果表明:新型支板尾部產(chǎn)生的激波會形成高溫高壓區(qū)域有助于燃燒,但也不可避免地帶來了總壓損失,激波只影響到支板的附近區(qū)域。F?rster等[15]對超燃沖壓發(fā)動機燃燒室內(nèi)超聲速和雙模態(tài)燃燒的現(xiàn)象進行了試驗研究,他們發(fā)現(xiàn),在反應(yīng)區(qū),火焰結(jié)構(gòu)和局部激波形態(tài)不同,不能僅通過評估靜壁壓力分布來區(qū)分燃燒模態(tài)。An等[16]對超聲速燃燒室激光點火和火花塞放電點火進行了對比試驗研究,結(jié)果表明,火花塞放電點火的點火過程明顯慢于激光點火。

    綜上所述,此前研究者已經(jīng)采用不同的方法對超燃沖壓發(fā)動機燃燒室的點火和火焰穩(wěn)定進行了研究,本文則采用多種非接觸式光學(xué)測量手段瞬時同步測量并詳細(xì)研究流場和火焰結(jié)構(gòu)整個動態(tài)發(fā)展過程,獲得了詳細(xì)試驗數(shù)據(jù),可為LES的大渦模擬計算提供驗證。

    1 試驗及測量方法介紹

    1.1 試驗設(shè)備和發(fā)動機構(gòu)型

    本文的試驗研究工作在中國空氣動力研究與發(fā)展中心(CARDC)的直連式超聲速脈沖燃燒風(fēng)洞(如圖1所示)上開展。該風(fēng)洞在之前的研究工作中[17-20]已經(jīng)有所介紹。試驗中,富氧空氣與氫氣在加熱器中混合燃燒,以產(chǎn)生高焓來流,流量為2.68 kg/s,總壓pt=0.82 MPa,總溫Tt=950 K。O2、H2O和N2的摩爾分?jǐn)?shù)分別為21%、12%和67%。來流氣體通過二維噴管加速到Ma=2.0。

    圖1 中國空氣動力研究與發(fā)展中心直連式脈沖燃燒風(fēng)洞Fig. 1 Direct-connected pulse combustion wind tunnel of CARDC

    如圖2所示,超燃沖壓發(fā)動機模型由長300 mm的隔離段和長773 mm的燃燒室組成。隔離段由長220 mm的等直段和長80 mm的擴張段組成,該擴張段的擴張角為1.4°,入口截面為30 mm×150 mm的矩形。凹腔深度11 mm,長121 mm(后斜坡角為21.1°,腔內(nèi)帶有1.4°的發(fā)散角),在凹腔內(nèi),距臺階下游71 mm處設(shè)有壓力監(jiān)測點,用于監(jiān)測燃燒狀態(tài)(圖2中Pressure monitor所示位置)。凹腔下游為擴張段,擴張段分為4小段:第一段長102.5 mm,帶有1.4°擴張角;第二段長145 mm,帶有2.0°擴張角;第三段長141.5 mm,帶有8.0°的擴張角;最后一段長263.0 mm,帶有15.0°的擴張角。2個火花塞位于凹腔臺階下游75 mm處,用于點燃先鋒氫。先鋒氫通過15個直徑為1.0 mm的圓孔噴注進入燃燒室,噴孔展向均布于凹腔臺階上游10 mm處。常溫乙烯通過直徑1.0 mm的10個圓孔從凹腔臺階下游25 mm處噴注進入燃燒室,進入后被先鋒氫火焰引燃。先鋒氫的當(dāng)量比約為0.33,乙烯的當(dāng)量比約為0.10。

    圖2 超燃沖壓發(fā)動機內(nèi)流道型面及燃燒室示意圖Fig. 2 Geometry of internal flow passage and the schematic diagram of a scramjet combustor

    設(shè)備運行的有效試驗時間約為440 ms,試驗操作時序如表1所示。為方便描述時序,將先鋒氫噴注進入燃燒室的時刻記為t=0,在t=-32.0 ms時冷流流場建立穩(wěn)定?;鸹ㄈぷ髁?次,第一次是在t=27.0 ms時,隨后先鋒氫被成功點燃。乙烯在t=100.0 ms時噴注進入燃燒室,隨后被先鋒氫火焰引燃。第二次是在t=207.0 ms時,先鋒氫停止噴注,乙烯繼續(xù)燃燒一段時間,在t=308.0 ms停止噴注乙烯。

    表1 試驗操作時序Table 1 Experimental operation sequence

    1.2 測量方法

    試驗采用壁面壓力測量、紋影、CH自發(fā)光、火焰自發(fā)光和OH-PLIF等非接觸式光學(xué)測量手段來顯示流場結(jié)構(gòu)和火焰發(fā)展過程。紋影拍攝的曝光時間設(shè)置為4.62 μs,拍攝幀率為10 000 幀/s。CH自發(fā)光和火焰自發(fā)光拍攝幀率都設(shè)置為5000 幀/s,曝光時間均為0.2 ms。CH自發(fā)光可用于標(biāo)記乙烯反應(yīng)區(qū)。燃燒室上壁面布置了一系列壓力傳感器,用于記錄燃燒室的沿程壁面壓力。壓力傳感器(量程0~700 kPa)的靈敏度為1%。采用平面激光誘導(dǎo)熒光法(PLIF)得到火焰點火后的OH*分布情況,它可以用來表征燃燒存在的區(qū)域和火焰的位置信息,并為其提供一個定性指標(biāo)。OH-PLIF系統(tǒng)的激光器采用哈爾濱工業(yè)大學(xué)研制的高重頻大能量Nd:YAG固體激光器,激光重復(fù)頻率為500 Hz,單脈沖激光能量30 mJ@532 nm,532 nm激光經(jīng)泵浦染料激光器倍頻輸出283.553 nm激光,單脈沖能量為1.6 mJ,此波段激發(fā)光具有較低的溫度敏感性。光路系統(tǒng)及相關(guān)測試方法在我們之前的工作中[2]已經(jīng)進行過詳細(xì)的說明。

    2 結(jié)果討論

    2.1 整個燃燒過程討論

    在整個燃燒過程中,監(jiān)測點(位于凹腔臺階下游71 mm處,圖2中Preesure monitor所示位置)隨時間變化的壓力曲線如圖3所示。圖中還給出了先鋒氫和乙烯的噴前壓力信號,用于標(biāo)記氫氣和乙烯的噴注時間,并結(jié)合監(jiān)測點的壓力變化標(biāo)示出各階段的時間節(jié)點。

    圖3 燃料(氫氣、乙烯)、加熱器和監(jiān)測點壓力隨時間的變化Fig. 3 Pressure of fuel (hydrogen, ethylene), heater and monitor point variation over time

    火花塞共工作了2次,分別是在t=27.0 ms和127.0 ms時。工作2次是為了保證點火的效果,當(dāng)出現(xiàn)第一次點火失敗的情況,2種燃料混合后再進行一次點火,可以確保燃料被引燃。若第一次點火成功,則第二次點火對燃燒過程和流場結(jié)構(gòu)幾乎沒有影響。由于火花塞放電的干擾,先鋒氫和乙烯的噴前壓力曲線上產(chǎn)生了2個較大峰值的干擾信號,可以用來標(biāo)記點火時刻。當(dāng)先鋒氫開始注入燃燒室時,監(jiān)測點的壓力在t=0時開始緩慢上升,然后保持穩(wěn)定,直到t=27.0 ms時先鋒氫被火花塞成功點燃,監(jiān)測點的壓力才突然上升。此時監(jiān)測點壓力值p達(dá)到了較高水平,約為0.31 MPa,說明先鋒氫點火成功。與此同時,可以看到乙烯的噴前壓力也略有上升,這是由于氫氣燃燒產(chǎn)生的壓升影響到了乙烯的噴前壓力所致。在t=100.0 ms時,監(jiān)測點壓力繼續(xù)升高,這是因為乙烯噴注進入燃燒室,立刻被先鋒氫成功引燃并和先鋒氫共同燃燒,此時監(jiān)測點壓力p在整個燃燒過程中達(dá)到最高水平,約為0.35 MPa?;旌先紵^程中監(jiān)測點壓力略有下降,原因是乙烯剛開始噴注產(chǎn)生的過沖現(xiàn)象在穩(wěn)定后被緩和。在t=207.0 ms時,先鋒氫停止噴注,監(jiān)測點壓力大幅度下降,降至略高于無反應(yīng)流的水平。此時雖然先鋒氫已停止噴注,但乙烯依舊能夠穩(wěn)定燃燒。t=308.0 ms時,乙烯停止噴注,火焰才逐漸熄滅。雖然控制噴注乙烯的電磁閥關(guān)閉了,但在噴注孔和電磁閥之間存在一段長度約400 mm、通徑為6 mm的金屬軟管,其中尚有余氣,從圖中可以看到,軟管中的乙烯完全排放干凈需要大約40.0 ms。

    綜上所述,整個燃燒過程可分為6個階段,如表2和圖3所示。第一階段為先鋒氫噴注進入燃燒室之前的無反應(yīng)流,時間為t=-32.0~0 ms(Part 1);第二階段為先鋒氫開始噴注后直到火花塞工作以前,用于顯示氫氣噴注時的流動特性,時間為t=0~27.0 ms(Part 2);第三階段用于描述氫氣的燃燒過程,包括點火和火焰穩(wěn)定,時間為t=27.0~100.0 ms(Part 3);第四階段為氫氣和乙烯共同燃燒的過程,此時燃燒最劇烈,時間為t=100.0~207.0 ms(Part 4);第五階段闡述了乙烯的燃燒特性,當(dāng)先鋒氫停止噴注后,乙烯能夠繼續(xù)穩(wěn)定燃燒,時間為t=207.0~308.0 ms(Part 5);第六階段是乙烯火焰熄滅過程直到試驗結(jié)束,時間為t=308.0~398.0 ms(Part 6)。

    表2 試驗過程的描述Table 2 Description of the test process

    第一至第五階段的時均沿程壁面壓力分布和發(fā)動機流道擴張面積曲線如圖4所示。因第六階段的沿程壁面壓力和第一階段差別不大,故沒有在圖中展示。雖然第三和第四階段的壁面壓力比其他幾個階段大得多,但這2個階段之間的差異并不大,反壓擾動的距離均約為0.2 m (圖4中x=0.1 m處),這是因為先鋒氫的當(dāng)量比(0.33)比乙烯的當(dāng)量比(0.10)大很多,在混合燃燒過程中,相對較少的乙烯在噴注進入燃燒室后造成的反壓增加效果相對不明顯。

    圖4 各階段上壁面中心線時均沿程壓力分布Fig. 4 The average pressure distribution on the centre line of the stages

    第一、第二和第五階段的主要區(qū)別在凹腔部分,這是因為3個階段的流場結(jié)構(gòu)不同,下文會作詳細(xì)分析。通過對比沿程壁面壓力,給出了不同階段壓力大小的排序。壓力從小到大依次是:第一階段、第二階段、第五階段、第三階段和第四階段。由于第六階段乙烯熄滅后的流場結(jié)構(gòu)和第一階段相同,故不再敘述。

    2.2 第一和第二階段的燃燒過程討論

    圖5為第一階段的試驗紋影照片和數(shù)值模擬流線圖。第一階段為無反應(yīng)流動,超聲速氣流經(jīng)隔離段進入凹腔后,在前臺階的剪切層(A)以及凹腔內(nèi)的低速區(qū)兩者共同作用下,產(chǎn)生了一道斜激波(B),這道斜激波在下壁面反射,并與邊界層相互作用,產(chǎn)生了一個分離區(qū)(E)。高速氣流通過分離區(qū)產(chǎn)生一道分離激波(C),并與斜激波(B)相互作用。同時,在凹腔前部和后部會產(chǎn)生2個回流區(qū)(D、F)。

    流場結(jié)構(gòu)周期性振蕩,本文條件下振蕩頻率約為400 Hz,流動振蕩主要是凹腔造成的。在一個振蕩周期內(nèi)(2.5 ms),分離激波(C)后高壓使得激波向上游運動,此時激波上行。當(dāng)激波繼續(xù)上行至凹腔后部斜坡時,波后的高壓與凹腔內(nèi)的低速區(qū)域混合,壓力降低導(dǎo)致激波減弱,分離區(qū)(E)也同時減弱,凹腔內(nèi)低速區(qū)域擴展,形成通道,入口速度恢復(fù)正常,此時激波下行。

    圖5 第一階段凹腔內(nèi)的紋影圖和計算流線圖Fig. 5 The schlieren image and the streamline of cavity in the first stage

    第二階段的試驗紋影圖像如圖6所示,它描述了先鋒氫噴注的全過程。在t=0時,由于氫氣的噴入,產(chǎn)生了一道激波(F),在t=0.3 ms時,可以看到它先是通過下壁面反射進而與斜激波(B)相互作用,這是第二階段的監(jiān)測點壓力略高于第一階段的原因。與后3幅圖像相比,除了在t=20.0 ms時圖像中出現(xiàn)了一些水(H)以外,其余圖像的流場結(jié)構(gòu)基本沒有發(fā)生變化。這是因為來流中含有高溫水蒸氣,在經(jīng)過溫度較低的玻璃窗時凝結(jié)成了小水滴,并順著玻璃向下游流動,此現(xiàn)象對流場結(jié)構(gòu)和燃燒過程幾乎沒有影響。

    圖6 第二階段凹腔內(nèi)的紋影照片F(xiàn)ig. 6 Schlieren images of cavity in the second stage

    2.3 第三階段的燃燒過程討論

    本節(jié)討論先鋒氫的點火和火焰發(fā)展過程,并使用紋影圖像、火焰自發(fā)光圖像以及OH-PLIF圖像來更好地了解流場結(jié)構(gòu)和火焰特性。如圖7所示,將點火過程定義為從火花塞開始工作一直到壁面壓力達(dá)到恒定水平的過程,點火過程歷時約為26.0 ms(t=27.0~53.0 ms)。第三階段整個過程的紋影圖像如圖8所示,燃燒穩(wěn)定后流場結(jié)構(gòu)變化不大,而點火過程中的流場結(jié)構(gòu)有很大變化?;鸹ㄈぷ鲿r(t=27.0 ms),可以看到一束強光(H)。先鋒氫立刻被點燃,燃燒釋熱非常迅速,僅0.1 ms后,即t=27.1 ms時,流場結(jié)構(gòu)便發(fā)生了很大變化。反射的激波(J)被壓縮到接近下壁面,在燃燒區(qū)之下(紅線以下)與下壁面之間來回反射。在t=27.2 ms時,因噴入氫氣而產(chǎn)生的激波(F)開始被推入隔離段中,激波的形狀發(fā)生了較大的變化,燃燒變得更加強烈,在t=29.0 ms時產(chǎn)生了一道“X”型的激波串(K)。激波串頭部被推到隔離段出口上游0.04 m處僅用了2.0 ms(t=29.0 ~31.0 ms),因此激波串的移動速度約為20 m/s。由于燃燒產(chǎn)生的反壓不斷增大,激波串在隔離段中不斷向上游移動,最終在t=39.0 ms時超出觀測范圍。對比后3幅圖像,在t=45.0 ms后,流場結(jié)構(gòu)趨于穩(wěn)定,紋影圖像幾乎沒有變化。在t=85.0 ms時,可以看到紋影圖像中產(chǎn)生了一道非常接近于正激波的激波(U),其后是一道斜激波(V),表明在此階段的燃燒非常強烈。

    圖7 第三階段氫氣點火過程Fig. 7 Hydrogen ignition process of the third stage

    圖8 第三階段的紋影圖像Fig. 8 Schlieren images of the third stage

    圖9為火焰的自發(fā)光圖像,用于顯示第三階段過程中的火焰分布及其發(fā)展過程。在t=27.0 ms時,可以看到一束非常亮的光,此時火花塞點燃先鋒氫。在t=29.0 ms時,火焰充滿整個凹腔,之后火焰區(qū)域繼續(xù)擴大,直到t=31.0 ms時,火焰接觸到下壁面,而此時激波串也被推入隔離段內(nèi)。火焰沿著下壁面向上游傳播,來回振蕩并且不穩(wěn)定,當(dāng)t=45.0 ms時又被推到了下游。對比后2張圖可以發(fā)現(xiàn),主燃燒區(qū)域還是在凹腔內(nèi),火焰非常明亮。和紋影圖像對比來看,火焰形狀的變化遠(yuǎn)遠(yuǎn)大于流場結(jié)構(gòu)的變化。

    圖9 第三階段火焰自發(fā)光圖像Fig. 9 Flame luminosity images of the third stage

    圖10為OH-PLIF圖像,用來描述化學(xué)反應(yīng)的主要區(qū)域。在t=27.0 ms時,火花塞首先點燃了火花塞周圍的氫氣,在火花塞附近可以看到OH自由基信號。在2.0 ms后,當(dāng)t=29.0 ms時,凹腔內(nèi)充滿了OH自由基信號。隨著時間的推移,OH自由基信號變化不大,主要分布在凹腔內(nèi),穩(wěn)焰模式為凹腔回流區(qū)穩(wěn)定燃燒。

    圖10 第三階段OH-PLIF圖像Fig. 10 OH-PLIF images of the third stage

    2.4 第四階段的燃燒過程討論

    本節(jié)介紹先鋒氫和乙烯在燃燒室中混合燃燒的特性。測量上增加了CH自發(fā)光圖像,用于標(biāo)記乙烯火焰的位置。第四階段的紋影圖像如圖11所示,在t=100.0 ms時,可以清楚地看到凹腔內(nèi)噴注乙烯所產(chǎn)生的噴射柱(L)。對比圖9中t=85.0 ms時的紋影圖像,在t=100.0 ms時凹腔內(nèi)的流場結(jié)構(gòu)變化不大。在此之后,燃燒變得更加強烈,激波串被推入隔離段內(nèi),在t=200.0 ms時超出觀測范圍。

    圖11 第四階段的紋影圖像Fig. 11 Schlieren images of the fourth stage

    第四階段的火焰自發(fā)光圖像如圖12所示。在t=100.0 ms時,凹腔內(nèi)可見明顯的強光(S),乙烯火焰為白色,氫氣火焰為黃色。對比接下來的3張圖片可以發(fā)現(xiàn),乙烯在注入燃燒室后立刻被點燃,這是因為先鋒氫燃燒產(chǎn)生了較高的溫度和壓力,并且產(chǎn)生了一個低速區(qū)域,有利于乙烯燃燒。乙烯火焰存在于凹腔臺階下游附近的剪切層內(nèi),穩(wěn)焰模式為凹腔剪切層穩(wěn)定燃燒。由于乙烯的當(dāng)量比遠(yuǎn)低于先鋒氫,因此穩(wěn)焰模式與先鋒氫(凹腔回流區(qū)穩(wěn)焰模式)有所不同。乙烯注入燃燒室后,2種氣體混合燃燒更加劇烈,火焰擾入隔離段。第四階段燃燒過程中CH自由基的自發(fā)光圖像如圖13所示。CH自由基是乙烯燃燒的中間產(chǎn)物,可以用來監(jiān)測乙烯燃料的燃燒效率。在t=100.0 ms時,凹腔內(nèi)可見少量CH自由基(R),并沿凹腔剪切層向下游擴散,最終火焰位于凹腔剪切層內(nèi)(S),與t=100.0~200.0 ms時的火焰自發(fā)光圖像分析得出的結(jié)果一致。

    圖12 第四階段的火焰自發(fā)光圖像Fig.12 Flame luminosity images of the fourth stage

    圖13 第四階段的CH自發(fā)光圖像Fig. 13 CH luminosity images of the fourth stage

    2.5 第五階段的燃燒過程討論

    當(dāng)先鋒氫在t=207.0 ms停止噴注時,由于乙烯比先鋒氫的當(dāng)量比小很多,因此監(jiān)測點壓力突然下降??梢钥闯?,在圖14中t=210.0 ms時,此前還在上游觀測范圍以外的激波串(M)回到了隔離段內(nèi)。隨著高速氣流繼續(xù)推動激波串向下游運動,在t=211.7 ms時激波串從隔離段中被推回到凹腔內(nèi)。最后,在t=250.0 ms時,激波串消失,替換成了2道斜激波。其中一道斜激波(N)是乙烯注入高速氣流而產(chǎn)生的;另一道斜激波(P)是乙烯燃燒產(chǎn)生的,因為燃燒釋熱將剪切層“抬進了”核心主流中,導(dǎo)致超聲速氣流流動通道變窄,從而產(chǎn)生了這道斜激波。

    圖14 第五階段的紋影圖像Fig. 14 Schlieren images of the fifth stage

    第五階段的火焰自發(fā)光圖像如圖15所示,由于先鋒氫停止噴注,燃燒變?nèi)?。與隔離段中的激波串一樣,火焰也在t=202.6 ms時從隔離段中被推回凹腔內(nèi)。在t=203.0 ms時,靠近下壁面底部的火焰消失,主要分布在凹腔內(nèi)和靠近上壁面處。最后,在t=250.0 ms時,火焰只存在于凹腔內(nèi)。在此階段,火焰區(qū)域變得越來越小,火焰的主色由黃色變?yōu)樗{(lán)色,乙烯火焰(S)在腔內(nèi)的位置從上游往下游移動。

    圖15 第五階段的火焰自發(fā)光圖像Fig. 15 Flame luminosity images of the fifth stage

    第五階段燃燒過程中的CH自發(fā)光圖像如圖16所示。無論先鋒氫是否停止噴注,乙烯的當(dāng)量比均保持恒定,所以CH自由基的發(fā)光面積也不會發(fā)生變化。與圖13相比,唯一的區(qū)別是CH自由基從靠近凹腔臺階的位置移動到靠近凹腔后部斜坡的位置,這與圖15中火焰自發(fā)光圖像的分析是一致的。圖17中,OH自由基信號變化較大,變化率與CH自由基信號相比要大很多。先鋒氫停止噴注前,OH自由基面積較大,在t=201.0和203.0 ms時主要分布在凹腔中后部和凹腔下游的上壁面。當(dāng)t=241.0 ms時,乙烯單獨燃燒,此時OH自由基僅由乙烯燃燒產(chǎn)生,因此,最終在t=251.0 ms時,OH自由基的面積變得很小,僅存在于凹腔的后部。

    圖16 第五階段的CH自發(fā)光圖像Fig. 16 CH luminosity images of the fifth stage

    圖17 第五階段的OH-PLIF圖像Fig. 17 OH-PLIF images of the fifth stage

    3 結(jié) 論

    采用試驗的方法在來流馬赫數(shù)Ma=2.0的條件下研究了以乙烯為燃料的直連式凹腔燃燒室的燃燒過程,結(jié)合壁面壓力曲線、紋影圖、火焰自發(fā)光照相和CH自發(fā)光圖像以及OH-PLIF圖像瞬態(tài)同步詳細(xì)分析了氫燃料點火、氫與乙烯共同燃燒以及乙烯單獨燃燒的流場結(jié)構(gòu)和火焰發(fā)展過程。燃燒過程分6個階段,第一階段描述了先鋒氫注入之前的無反應(yīng)流動,試驗測得流場結(jié)構(gòu)周期性振蕩的頻率約為400 Hz。第二階段描述了先鋒氫被點燃之前的流動特性。第三階段描述了氫氣的燃燒過程,包括點火和火焰發(fā)展過程。第四階段描述了氫氣和乙烯的混合燃燒過程,此時燃燒最劇烈,當(dāng)先鋒氫停止噴注后,乙烯可以繼續(xù)穩(wěn)定燃燒。第五階段描述了乙烯的燃燒特性。

    本文研究重點關(guān)注以下3個方面:1) 無反應(yīng)流動(第一階段)的流場結(jié)構(gòu)呈周期性振蕩,本文條件下振蕩頻率約為400 Hz,這是超聲速來流下凹腔燃燒室構(gòu)型固有特性。2) 噴氫點火過程(第三階段),火花塞點火(t=27.0 ms)僅0.1 ms后,即t=27.1 ms時,流場與火焰結(jié)構(gòu)發(fā)生劇變。燃燒非常強烈,2.0 ms內(nèi)(t=29.0~ 31.0 ms),激波串推動距離為0.04 m,移動速度約為20 m/s,顯示了氫點火反應(yīng)動力學(xué)的重大影響。3) 氫和乙烯混合燃燒及乙烯單獨燃燒中,CH自由基是乙烯燃燒的中間產(chǎn)物,初步分析表明,乙烯的CH自發(fā)光圖像可以用來判斷其燃燒效率的高低。

    上述整個燃燒流場的詳細(xì)試驗測量結(jié)果可為動態(tài)CFD計算與驗證提供依據(jù)。

    猜你喜歡
    凹腔激波燃燒室
    凹腔對高超聲速邊界層穩(wěn)定性的影響
    燃燒室形狀對國六柴油機性能的影響
    一種基于聚類分析的二維激波模式識別算法
    基于HIFiRE-2超燃發(fā)動機內(nèi)流道的激波邊界層干擾分析
    一種熱電偶在燃燒室出口溫度場的測量應(yīng)用
    電子制作(2019年19期)2019-11-23 08:41:54
    縫翼凹腔擋板氣動性能和降噪效果數(shù)值研究
    斜激波入射V形鈍前緣溢流口激波干擾研究
    適于可壓縮多尺度流動的緊致型激波捕捉格式
    壁面噴射當(dāng)量比對支板凹腔耦合燃燒的影響
    高幾何壓縮比活塞的燃燒室形狀探討
    18禁国产床啪视频网站| 熟女电影av网| 国产精品综合久久久久久久免费| 搡老岳熟女国产| 美女黄网站色视频| 亚洲精品中文字幕一二三四区| 日韩高清综合在线| 国产视频内射| 两个人视频免费观看高清| 日韩大码丰满熟妇| 成人三级做爰电影| 岛国在线观看网站| 久久亚洲精品不卡| 日韩欧美免费精品| 亚洲成a人片在线一区二区| 精品日产1卡2卡| 欧美久久黑人一区二区| 亚洲色图av天堂| 午夜两性在线视频| 亚洲av成人av| 在线观看免费午夜福利视频| 国产片内射在线| 国产av麻豆久久久久久久| 日韩大尺度精品在线看网址| 久久精品国产亚洲av高清一级| 麻豆久久精品国产亚洲av| 婷婷亚洲欧美| 亚洲av成人精品一区久久| 国产亚洲欧美98| bbb黄色大片| 1024手机看黄色片| 在线观看免费日韩欧美大片| 一级片免费观看大全| 18禁国产床啪视频网站| 在线视频色国产色| 国产野战对白在线观看| 久久久久久人人人人人| 天天躁夜夜躁狠狠躁躁| 男人舔女人下体高潮全视频| 国内毛片毛片毛片毛片毛片| 免费在线观看亚洲国产| 精品国内亚洲2022精品成人| 精品不卡国产一区二区三区| 亚洲全国av大片| 99久久精品热视频| 亚洲av熟女| 亚洲成人久久爱视频| 国产精品 欧美亚洲| 亚洲av成人一区二区三| 国内少妇人妻偷人精品xxx网站 | 欧美日韩亚洲国产一区二区在线观看| 韩国av一区二区三区四区| 日本黄大片高清| 成人三级黄色视频| 九九热线精品视视频播放| 免费看十八禁软件| 亚洲第一电影网av| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲精品第一综合不卡| 九色成人免费人妻av| 国产精品久久久久久人妻精品电影| 18禁黄网站禁片午夜丰满| 欧美大码av| 757午夜福利合集在线观看| 麻豆久久精品国产亚洲av| 国产av一区二区精品久久| 欧美性猛交╳xxx乱大交人| 亚洲五月婷婷丁香| 欧洲精品卡2卡3卡4卡5卡区| 久久精品亚洲精品国产色婷小说| 少妇粗大呻吟视频| 天堂√8在线中文| 一个人观看的视频www高清免费观看 | 午夜视频精品福利| 国产三级黄色录像| 丝袜人妻中文字幕| 美女扒开内裤让男人捅视频| 又黄又爽又免费观看的视频| 啦啦啦免费观看视频1| 久久精品综合一区二区三区| 一本久久中文字幕| 午夜精品在线福利| 久久99热这里只有精品18| 黄色毛片三级朝国网站| 精品久久久久久久久久久久久| 免费在线观看亚洲国产| 日本撒尿小便嘘嘘汇集6| 无人区码免费观看不卡| 十八禁人妻一区二区| 午夜福利在线观看吧| 一个人观看的视频www高清免费观看 | 床上黄色一级片| 欧美+亚洲+日韩+国产| 999精品在线视频| 国产三级中文精品| av中文乱码字幕在线| 日韩欧美国产一区二区入口| 国产高清视频在线播放一区| 又粗又爽又猛毛片免费看| 人人妻,人人澡人人爽秒播| 精品久久久久久,| 欧美在线黄色| 国产精品久久久av美女十八| 老鸭窝网址在线观看| 国产成人欧美在线观看| 国产精品 欧美亚洲| 成年免费大片在线观看| 国产一区二区三区视频了| 久久久久免费精品人妻一区二区| 超碰成人久久| 变态另类成人亚洲欧美熟女| 老司机深夜福利视频在线观看| 在线播放国产精品三级| 美女扒开内裤让男人捅视频| 久久午夜综合久久蜜桃| 欧美大码av| 看片在线看免费视频| av片东京热男人的天堂| 国产成人欧美在线观看| 国产av又大| 国产97色在线日韩免费| 成人三级做爰电影| 精品一区二区三区四区五区乱码| 国产熟女午夜一区二区三区| 精品无人区乱码1区二区| 久久久久久亚洲精品国产蜜桃av| 日本在线视频免费播放| 国产精品一区二区免费欧美| 怎么达到女性高潮| 久久久久久久久中文| 国产精品一区二区三区四区免费观看 | 日韩欧美三级三区| 老汉色av国产亚洲站长工具| 日韩高清综合在线| 免费在线观看成人毛片| 成人精品一区二区免费| 激情在线观看视频在线高清| 99久久久亚洲精品蜜臀av| 在线永久观看黄色视频| 99精品欧美一区二区三区四区| 99久久无色码亚洲精品果冻| 国内精品久久久久久久电影| 啦啦啦韩国在线观看视频| 不卡av一区二区三区| 两个人的视频大全免费| 村上凉子中文字幕在线| 丰满的人妻完整版| 国产在线观看jvid| 亚洲国产看品久久| 国产午夜精品久久久久久| 亚洲人成电影免费在线| www.精华液| 欧美日韩国产亚洲二区| 悠悠久久av| 三级男女做爰猛烈吃奶摸视频| 麻豆成人午夜福利视频| 三级男女做爰猛烈吃奶摸视频| 桃色一区二区三区在线观看| 久久这里只有精品19| 精品熟女少妇八av免费久了| 他把我摸到了高潮在线观看| 久久久久久久久久黄片| 中文字幕人成人乱码亚洲影| 国产主播在线观看一区二区| 国产三级在线视频| 亚洲片人在线观看| 18禁国产床啪视频网站| 欧美午夜高清在线| 黄色视频,在线免费观看| 18禁裸乳无遮挡免费网站照片| 免费在线观看亚洲国产| 国产精品亚洲av一区麻豆| 精品欧美一区二区三区在线| 久久香蕉精品热| 老汉色av国产亚洲站长工具| 狂野欧美激情性xxxx| 岛国在线观看网站| 国产精品久久久人人做人人爽| 久久午夜亚洲精品久久| 超碰成人久久| 哪里可以看免费的av片| 老司机福利观看| 午夜两性在线视频| 香蕉丝袜av| 久久伊人香网站| 亚洲一区二区三区色噜噜| 国产精品爽爽va在线观看网站| 久久久国产成人精品二区| 欧美日韩精品网址| 亚洲国产欧美人成| 国产久久久一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 成人av在线播放网站| 国产亚洲精品av在线| 人人妻,人人澡人人爽秒播| 精品久久久久久,| 男人舔奶头视频| 丝袜人妻中文字幕| 色综合欧美亚洲国产小说| 国产亚洲精品一区二区www| 色哟哟哟哟哟哟| 午夜免费观看网址| 精品久久久久久久毛片微露脸| 欧美性猛交╳xxx乱大交人| 美女黄网站色视频| 国产真人三级小视频在线观看| 99在线人妻在线中文字幕| av超薄肉色丝袜交足视频| 亚洲一区高清亚洲精品| netflix在线观看网站| 免费在线观看影片大全网站| 国产精品香港三级国产av潘金莲| 九九热线精品视视频播放| 丰满的人妻完整版| 国产一区在线观看成人免费| 国产精品一及| 黄色视频,在线免费观看| 久久久久国内视频| 女人被狂操c到高潮| 国产在线精品亚洲第一网站| 亚洲av片天天在线观看| 特级一级黄色大片| 大型黄色视频在线免费观看| 夜夜看夜夜爽夜夜摸| 91字幕亚洲| 亚洲精品美女久久av网站| a级毛片在线看网站| 欧美成狂野欧美在线观看| 成年免费大片在线观看| 人人妻人人看人人澡| 久久久久久人人人人人| 亚洲激情在线av| 俺也久久电影网| 欧美极品一区二区三区四区| 国产精品亚洲美女久久久| 国产男靠女视频免费网站| 国产成年人精品一区二区| 一本久久中文字幕| 91老司机精品| 成人特级黄色片久久久久久久| 一进一出好大好爽视频| 日本a在线网址| 中出人妻视频一区二区| 久久精品国产综合久久久| 午夜免费成人在线视频| 亚洲人成77777在线视频| 别揉我奶头~嗯~啊~动态视频| 成人一区二区视频在线观看| 国产亚洲精品久久久久久毛片| 色综合亚洲欧美另类图片| 久久精品成人免费网站| 国产aⅴ精品一区二区三区波| 日本成人三级电影网站| 久久国产乱子伦精品免费另类| 欧美黄色片欧美黄色片| 久久天躁狠狠躁夜夜2o2o| 国产成人啪精品午夜网站| 亚洲欧洲精品一区二区精品久久久| 首页视频小说图片口味搜索| 久久草成人影院| 啦啦啦观看免费观看视频高清| 一二三四在线观看免费中文在| 亚洲欧美日韩无卡精品| videosex国产| 久久久精品国产亚洲av高清涩受| 人人妻人人看人人澡| 亚洲美女黄片视频| 18禁美女被吸乳视频| 国产激情久久老熟女| 国产亚洲精品久久久久久毛片| 亚洲熟妇中文字幕五十中出| 黄色片一级片一级黄色片| 校园春色视频在线观看| 国产精品久久久久久人妻精品电影| 99久久精品国产亚洲精品| 国产激情欧美一区二区| 久久久久久免费高清国产稀缺| 亚洲激情在线av| 俄罗斯特黄特色一大片| 欧美丝袜亚洲另类 | 嫁个100分男人电影在线观看| 亚洲人成网站高清观看| 老司机靠b影院| 国产精品爽爽va在线观看网站| 最新在线观看一区二区三区| 日韩国内少妇激情av| avwww免费| 久久久久国内视频| 脱女人内裤的视频| 午夜福利在线在线| 制服诱惑二区| 动漫黄色视频在线观看| 国模一区二区三区四区视频 | 天天躁狠狠躁夜夜躁狠狠躁| 97超级碰碰碰精品色视频在线观看| 不卡一级毛片| 午夜福利视频1000在线观看| 中文亚洲av片在线观看爽| 狠狠狠狠99中文字幕| 久久久久久久久中文| 国产精品一区二区三区四区久久| 久久午夜综合久久蜜桃| 99久久精品热视频| 亚洲国产中文字幕在线视频| 国产又色又爽无遮挡免费看| 少妇的丰满在线观看| 伊人久久大香线蕉亚洲五| 亚洲色图 男人天堂 中文字幕| 精品不卡国产一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 91国产中文字幕| 免费无遮挡裸体视频| 制服丝袜大香蕉在线| 中文字幕高清在线视频| 久久久久久久午夜电影| 国产av麻豆久久久久久久| 在线免费观看的www视频| 亚洲精品国产一区二区精华液| 人人妻人人澡欧美一区二区| 夜夜夜夜夜久久久久| 18禁观看日本| 亚洲人与动物交配视频| 国产精品综合久久久久久久免费| 1024手机看黄色片| 国产一区二区三区视频了| 国产精品久久视频播放| 亚洲专区国产一区二区| 亚洲五月天丁香| 国产精品乱码一区二三区的特点| 国产高清有码在线观看视频 | 琪琪午夜伦伦电影理论片6080| xxx96com| 1024香蕉在线观看| 欧美中文综合在线视频| 国产在线观看jvid| 97超级碰碰碰精品色视频在线观看| 不卡一级毛片| 禁无遮挡网站| 麻豆成人av在线观看| 一夜夜www| 99国产精品99久久久久| 久久久久久久午夜电影| 国产精品电影一区二区三区| www.999成人在线观看| 成人三级做爰电影| 亚洲欧美日韩无卡精品| aaaaa片日本免费| 亚洲精品在线美女| 一本久久中文字幕| 麻豆成人av在线观看| 免费人成视频x8x8入口观看| 欧美日韩瑟瑟在线播放| 91成年电影在线观看| 精品第一国产精品| 免费电影在线观看免费观看| 日本三级黄在线观看| 村上凉子中文字幕在线| 免费在线观看亚洲国产| 天天躁夜夜躁狠狠躁躁| 在线观看日韩欧美| 在线观看www视频免费| a级毛片在线看网站| 亚洲自偷自拍图片 自拍| 欧美一区二区国产精品久久精品 | 999久久久精品免费观看国产| 久久久水蜜桃国产精品网| 色尼玛亚洲综合影院| 99re在线观看精品视频| 国产v大片淫在线免费观看| 久久精品人妻少妇| 午夜老司机福利片| 999精品在线视频| 久久久久九九精品影院| 激情在线观看视频在线高清| 18禁美女被吸乳视频| 99国产精品99久久久久| 国产一区二区三区在线臀色熟女| 人人妻人人看人人澡| 19禁男女啪啪无遮挡网站| 国产人伦9x9x在线观看| 国产主播在线观看一区二区| 亚洲国产中文字幕在线视频| 久久亚洲精品不卡| 日韩欧美在线二视频| 波多野结衣高清无吗| 91国产中文字幕| 亚洲国产高清在线一区二区三| 18禁裸乳无遮挡免费网站照片| 岛国在线观看网站| 午夜亚洲福利在线播放| 精品高清国产在线一区| 欧美性猛交黑人性爽| 国产高清videossex| 色在线成人网| 欧美黄色淫秽网站| 亚洲乱码一区二区免费版| 在线观看免费午夜福利视频| 亚洲av成人不卡在线观看播放网| 亚洲人与动物交配视频| 国产精品久久久久久亚洲av鲁大| 国内精品久久久久久久电影| 男人舔奶头视频| 女警被强在线播放| 999久久久国产精品视频| 麻豆成人午夜福利视频| 亚洲 国产 在线| 久久九九热精品免费| 亚洲色图av天堂| 国内毛片毛片毛片毛片毛片| 国产免费男女视频| av视频在线观看入口| 麻豆一二三区av精品| 国产一区二区激情短视频| 91老司机精品| 变态另类丝袜制服| 国产成人精品久久二区二区免费| 天天添夜夜摸| 亚洲一区二区三区不卡视频| 午夜亚洲福利在线播放| 欧美激情久久久久久爽电影| 亚洲av第一区精品v没综合| 国产熟女xx| 国产一区二区三区在线臀色熟女| 女人爽到高潮嗷嗷叫在线视频| 久久久久久大精品| 天堂av国产一区二区熟女人妻 | 18禁黄网站禁片午夜丰满| 亚洲中文字幕日韩| 色播亚洲综合网| 欧美性猛交╳xxx乱大交人| 精品午夜福利视频在线观看一区| 国产黄片美女视频| 制服丝袜大香蕉在线| 91在线观看av| www日本黄色视频网| 免费在线观看亚洲国产| 在线观看www视频免费| 久久久久国内视频| 久久久久免费精品人妻一区二区| 亚洲精品粉嫩美女一区| 制服人妻中文乱码| 亚洲精品久久成人aⅴ小说| 久久久国产成人精品二区| 欧美又色又爽又黄视频| 久久香蕉激情| 亚洲国产高清在线一区二区三| 久久精品91无色码中文字幕| 757午夜福利合集在线观看| 国产亚洲精品久久久久5区| 久久精品国产综合久久久| 国产一区二区激情短视频| 一级毛片女人18水好多| 18禁国产床啪视频网站| 高清在线国产一区| 久久精品aⅴ一区二区三区四区| 久久久精品国产亚洲av高清涩受| 日本一本二区三区精品| 国产不卡一卡二| 在线看三级毛片| 18禁观看日本| 后天国语完整版免费观看| 三级男女做爰猛烈吃奶摸视频| 亚洲aⅴ乱码一区二区在线播放 | 午夜免费观看网址| 日本在线视频免费播放| 巨乳人妻的诱惑在线观看| 免费在线观看成人毛片| 老熟妇仑乱视频hdxx| 成年人黄色毛片网站| 国产精品亚洲一级av第二区| 精品久久久久久久末码| 老熟妇乱子伦视频在线观看| 午夜亚洲福利在线播放| 女人爽到高潮嗷嗷叫在线视频| 久久久精品欧美日韩精品| 成人一区二区视频在线观看| www.999成人在线观看| 久久国产精品人妻蜜桃| 黄色女人牲交| 精品久久久久久成人av| 99re在线观看精品视频| 久久久精品大字幕| 免费看美女性在线毛片视频| 欧美成人午夜精品| 久久久精品国产亚洲av高清涩受| 久久亚洲精品不卡| 国产av麻豆久久久久久久| 亚洲av第一区精品v没综合| 久久香蕉精品热| 中文字幕精品亚洲无线码一区| 五月伊人婷婷丁香| 亚洲人成网站高清观看| 精品国产乱码久久久久久男人| 国产乱人伦免费视频| 好男人在线观看高清免费视频| 婷婷精品国产亚洲av在线| 午夜福利成人在线免费观看| АⅤ资源中文在线天堂| 熟妇人妻久久中文字幕3abv| 18禁裸乳无遮挡免费网站照片| 18禁观看日本| www.精华液| 天堂√8在线中文| 香蕉久久夜色| 男插女下体视频免费在线播放| 久久久久亚洲av毛片大全| 成人亚洲精品av一区二区| 亚洲精品美女久久久久99蜜臀| 欧美日韩国产亚洲二区| 亚洲精品久久国产高清桃花| 亚洲av成人不卡在线观看播放网| 欧美国产日韩亚洲一区| 淫秽高清视频在线观看| 国产乱人伦免费视频| 精品熟女少妇八av免费久了| 国产精品香港三级国产av潘金莲| 窝窝影院91人妻| 成人午夜高清在线视频| 午夜视频精品福利| 精品久久久久久成人av| 国产亚洲欧美在线一区二区| 欧美不卡视频在线免费观看 | 搡老岳熟女国产| 亚洲av成人精品一区久久| 91av网站免费观看| 亚洲熟女毛片儿| 中文字幕精品亚洲无线码一区| 男女做爰动态图高潮gif福利片| 免费一级毛片在线播放高清视频| 亚洲av熟女| 999精品在线视频| 国产v大片淫在线免费观看| 天堂动漫精品| 99国产精品一区二区三区| 亚洲五月婷婷丁香| 亚洲国产精品sss在线观看| 男女那种视频在线观看| 久久人妻福利社区极品人妻图片| 日本 av在线| 国内精品久久久久精免费| 一夜夜www| 不卡一级毛片| 人妻久久中文字幕网| 麻豆国产97在线/欧美 | 少妇裸体淫交视频免费看高清 | 夜夜爽天天搞| 天天躁夜夜躁狠狠躁躁| 成人18禁在线播放| 欧美性猛交╳xxx乱大交人| 全区人妻精品视频| 国产97色在线日韩免费| 在线免费观看的www视频| 一本大道久久a久久精品| 制服诱惑二区| 男女下面进入的视频免费午夜| 国产激情偷乱视频一区二区| 久久亚洲真实| 成人国产一区最新在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲性夜色夜夜综合| 中文亚洲av片在线观看爽| 久久久国产成人精品二区| 亚洲熟妇中文字幕五十中出| 亚洲午夜精品一区,二区,三区| 国产真实乱freesex| 国产精品av视频在线免费观看| 日韩欧美在线乱码| 国产亚洲av嫩草精品影院| 巨乳人妻的诱惑在线观看| 日日摸夜夜添夜夜添小说| 91av网站免费观看| 成人av一区二区三区在线看| 亚洲精品国产一区二区精华液| 欧美一区二区国产精品久久精品 | 日日干狠狠操夜夜爽| 人妻夜夜爽99麻豆av| 精品久久久久久成人av| 亚洲人成77777在线视频| 久久精品91蜜桃| 狂野欧美激情性xxxx| 国产精品九九99| 亚洲一区中文字幕在线| 高清毛片免费观看视频网站| 亚洲成人久久性| 午夜福利在线观看吧| 亚洲第一欧美日韩一区二区三区| 午夜两性在线视频| а√天堂www在线а√下载| 亚洲av五月六月丁香网| 桃红色精品国产亚洲av| 每晚都被弄得嗷嗷叫到高潮| 在线观看免费午夜福利视频| 波多野结衣高清无吗| 非洲黑人性xxxx精品又粗又长| 久久久久亚洲av毛片大全| 久久久久国产一级毛片高清牌| 99久久精品热视频| 国产99久久九九免费精品| 国产精品av久久久久免费| 亚洲色图av天堂| 成人手机av| 亚洲av成人不卡在线观看播放网| 精品久久蜜臀av无| 99久久精品热视频| 一级毛片女人18水好多| 日韩欧美免费精品| av免费在线观看网站| 午夜精品久久久久久毛片777| 99国产极品粉嫩在线观看| 亚洲精品国产一区二区精华液| 亚洲精品美女久久久久99蜜臀| 9191精品国产免费久久| 老鸭窝网址在线观看| 人妻夜夜爽99麻豆av| 国产成人欧美在线观看|