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

    稀釋氣體對二甲醚著火延遲特性的影響

    2017-06-05 15:00:30張紅光李佳政石智成
    關(guān)鍵詞:二甲醚熱效應(yīng)混合氣體

    張紅光 李佳政 石智成 高 翔

    (1.北京工業(yè)大學(xué)環(huán)境與能源工程學(xué)院, 北京 100124; 2.北京電動車輛協(xié)同創(chuàng)新中心, 北京 100124)

    稀釋氣體對二甲醚著火延遲特性的影響

    張紅光1,2李佳政1,2石智成1,2高 翔1,2

    (1.北京工業(yè)大學(xué)環(huán)境與能源工程學(xué)院, 北京 100124; 2.北京電動車輛協(xié)同創(chuàng)新中心, 北京 100124)

    在上止點(diǎn)壓力1.5 MPa,上止點(diǎn)溫度670~795 K,當(dāng)量比為1的實(shí)驗(yàn)條件下,利用快速壓縮機(jī)(RCM)研究了N2、N2/Ar(50%/50%)和Ar/CO2(61.2%/38.8%)稀釋氣體對二甲醚(DME)著火延遲期的影響。利用CHMKIN-PRO軟件在較大溫度范圍內(nèi)進(jìn)行了模擬研究。結(jié)果表明:不同稀釋氣體組分對二甲醚混合氣體第一階段著火延遲期影響很小,但對總著火延遲期具有較大的影響,尤其在負(fù)溫度系數(shù)(NTC)區(qū)間的起點(diǎn)溫度附近。與N2相比,N2/Ar(50%/50%)使得實(shí)驗(yàn)中總著火延遲期縮短30%。稀釋氣體的化學(xué)效應(yīng)對第一階段著火延遲期和總著火延遲期的影響均很小。在低溫區(qū)間和NTC區(qū)間,稀釋氣體的熱效應(yīng)起主導(dǎo)作用,而當(dāng)上止點(diǎn)溫度超過NTC區(qū)間時,CO2的化學(xué)效應(yīng)影響增強(qiáng)且超過熱效應(yīng)而起到主導(dǎo)作用。隨稀釋率增加,DME混合氣體第一階段著火延遲期略有延長,而總著火延遲期出現(xiàn)明顯延長,且總著火延遲期的NTC現(xiàn)象變得更加明顯。

    稀釋氣體; 二甲醚; 著火延遲期; 快速壓縮機(jī)

    引言

    日益嚴(yán)重的環(huán)境問題使得人們更加重視內(nèi)燃機(jī)的燃燒與污染排放。為了減少內(nèi)燃機(jī)的有害排放,理解燃料燃燒時發(fā)生的化學(xué)反應(yīng)是十分必要的??焖賶嚎s機(jī)(RCM)能夠模擬內(nèi)燃機(jī)中活塞的單次壓縮過程,且便于精準(zhǔn)控制壓縮過程中的各項(xiàng)熱力學(xué)參數(shù),從而成為研究燃料燃燒特性的理想實(shí)驗(yàn)平臺[1]。在快速壓縮機(jī)實(shí)驗(yàn)過程中,為達(dá)到不同的測試工況,常常通過調(diào)節(jié)稀釋氣體混合比例來實(shí)現(xiàn)。但是,由于不同稀釋氣體比熱容比、熱擴(kuò)散率、導(dǎo)熱性等的不同,相同燃料采用不同稀釋氣體所得實(shí)驗(yàn)結(jié)果有所差異。

    針對稀釋氣體對快速壓縮機(jī)中燃料燃燒過程的影響,已有學(xué)者展開相關(guān)研究[2-3]。目前已有研究主要針對單階段著火或高溫燃燒,而稀釋氣體對于具有兩階段著火及負(fù)溫度系數(shù)(NTC)現(xiàn)象燃料燃燒影響的研究很少,近來只有DI等[4]通過實(shí)驗(yàn)和數(shù)值模擬研究了稀釋氣體熱效應(yīng)和化學(xué)效應(yīng)對異辛烷和正庚烷兩種燃料著火過程的影響。研究發(fā)現(xiàn):不同稀釋氣體組分對正庚烷和異辛烷著火影響很大,尤其是在兩級著火區(qū)間。低到中溫范圍內(nèi),二甲醚(DME)燃料通常具有明顯的兩階段燃燒和NTC現(xiàn)象,而且,不同燃料對采用稀釋氣體種類的敏感性不同[5],稀釋氣體組分對二甲醚燃料燃燒過程的影響尚不清楚。基于此,本文利用RCM開展稀釋氣體組分對二甲醚著火延遲期影響的實(shí)驗(yàn)研究,并結(jié)合CHEMKIN-PRO軟件進(jìn)行模擬與化學(xué)動力學(xué)分析。

    1 實(shí)驗(yàn)裝置與參數(shù)定義

    1.1 快速壓縮機(jī)

    使用的RCM參見文獻(xiàn)[6],主要由氣壓驅(qū)動系統(tǒng)、液壓控制系統(tǒng)、燃燒系統(tǒng)、配氣系統(tǒng)、數(shù)據(jù)采集系統(tǒng)組成。燃燒缸直徑50 mm,驅(qū)動氣壓0.2~0.8 MPa,壓縮比8.42~16.9。本文RCM的燃燒活塞采用文獻(xiàn)[7-13]中的“creviced piston”,這種活塞設(shè)計(jì)可有效減小活塞運(yùn)動過程中引起的渦流對燃燒室的影響,保證“絕熱核心假設(shè)”的成立,從而精確地計(jì)算上止點(diǎn)溫度Tc。Tc計(jì)算公式為

    (1)

    式中p0——混合氣體的初始壓力,MPaT0——初始溫度,Kpc——上止點(diǎn)壓力,MPaγ——混合氣體的比熱容比

    圖1為實(shí)驗(yàn)臺架原理圖。

    圖1 RCM實(shí)驗(yàn)臺架原理圖Fig.1 Schematic diagram of rapid compression machine (RCM) test bench1.同步信號發(fā)生器 2.高速攝影機(jī) 3.數(shù)據(jù)采集 4.加熱保溫系統(tǒng) 5.磁力預(yù)混罐 6.氮?dú)?7.氧氣 8.二甲醚 9.真空泵10.燃燒缸 11.液壓缸 12.空氣壓縮機(jī) 13.儲氣缸 14.驅(qū)動氣缸 15.位移傳感器 16.限位墊片 17.液壓控制站 18.余隙墊片 19.示波器 20.電荷放大器

    1.2 著火延遲期定義

    圖2 著火延遲期的定義Fig.2 Definition of ignition delay

    二甲醚混合氣著火延遲期的定義如圖2所示,定義時間0點(diǎn)為活塞到達(dá)上止點(diǎn)的時刻,從時間0點(diǎn)到燃燒壓力變化率的第一個峰值定義為第一階段著火延遲期,從燃燒壓力變化率的第一個峰值到第二個峰值定義為第二階段著火延遲期。從圖中可以明顯看出,二甲醚混合氣體燃燒呈現(xiàn)出明顯的兩階段放熱和兩階段著火延遲,總著火延遲期τ=τ1+τ2。

    1.3 實(shí)驗(yàn)條件及實(shí)驗(yàn)可重復(fù)性

    稀釋氣體對著火過程的影響分為熱效應(yīng)(比熱容)、化學(xué)效應(yīng)(第三體碰撞)和稀釋效應(yīng)(燃料和O2濃度)。本文利用快速壓縮機(jī)研究N2、N2/Ar(50%/50%)和Ar/CO2(61.2%/38.8%)稀釋氣體對二甲醚著火延遲期的影響。為了將稀釋氣體熱效應(yīng)和化學(xué)效應(yīng)區(qū)分開研究,將Ar和CO2以61.2%/38.8%的摩爾比例混合,其比熱容與N2基本相同。實(shí)驗(yàn)中氣體純度依次為:DME 99.5%、O299.999%、N299.999%、Ar 99.999%、CO299.999%。實(shí)驗(yàn)中上止點(diǎn)壓力通過調(diào)節(jié)混合氣體初始壓力來控制,上止點(diǎn)溫度通過調(diào)節(jié)壓縮比和混合氣體初始溫度來控制,實(shí)驗(yàn)條件如表1所示。

    實(shí)驗(yàn)過程中,缸內(nèi)燃燒壓力數(shù)據(jù)由kistler6125C型缸壓傳感器與kistler5011型電荷放大器進(jìn)行采集,并用TektronixMS04000型示波器記錄數(shù)據(jù)。每組實(shí)驗(yàn)條件根據(jù)情況重復(fù)3~5次,以保證實(shí)驗(yàn)數(shù)據(jù)的可靠性。本文中同一工況實(shí)驗(yàn)可重復(fù)性按著火延遲期相對誤差Δα來評價[4],經(jīng)計(jì)算圖3中著火延遲期相對誤差Δα為2.5%,誤差在可接受范圍內(nèi),且從圖中可以看出,DME混合氣體壓縮及燃燒曲線吻合較好,即基于本RCM的燃燒實(shí)驗(yàn)可重復(fù)性良好。

    表1 實(shí)驗(yàn)條件

    圖3 實(shí)驗(yàn)可重復(fù)性Fig.3 Repeatability of experiment

    2 數(shù)值模擬

    利用CHEMKIN-PRO軟件模擬二甲醚混合氣在快速壓縮機(jī)內(nèi)燃燒時的物理及化學(xué)反應(yīng)條件,模擬過程中采用ZHAO機(jī)理[14],該機(jī)理包含55種組分和290步可逆基元反應(yīng),是國際上公認(rèn)的、可靠的二甲醚燃燒詳細(xì)機(jī)理,并在快速壓縮機(jī)[15]、激波管[16]和流動反應(yīng)器[17]等燃燒裝置中得到驗(yàn)證。由于快速壓縮機(jī)壓縮過程非常快,若忽略壓縮過程中的傳熱損失,則活塞到達(dá)上止點(diǎn)后近似為絕熱定容反應(yīng),因此,本文反應(yīng)模型采用0-D均質(zhì)定容絕熱反應(yīng)器。

    反應(yīng)器初始參數(shù)設(shè)為實(shí)驗(yàn)中RCM壓縮上止點(diǎn)的參數(shù),利用MITTAL等[18]采用的“體積膨脹法”來模擬傳熱損失。在計(jì)算中將壓力的變化以等熵假設(shè)轉(zhuǎn)換為體積的變化。瞬時體積計(jì)算公式為

    Veff(t)=Veff(0)vp(t) (t>0)

    (2)

    式中Veff(0)——活塞達(dá)到上止點(diǎn)時的燃燒室體積Veff(t)——活塞達(dá)到上止點(diǎn)后的瞬時體積vp(t)——體積膨脹系數(shù)

    具體方法參見文獻(xiàn)[19]。

    3 結(jié)果與討論

    3.1 稀釋氣體的熱效應(yīng)與化學(xué)效應(yīng)

    圖4為稀釋氣體組分對二甲醚第一階段著火延遲期和總著火延遲期影響的實(shí)驗(yàn)與模擬結(jié)果。如圖4所示,ZHAO機(jī)理能夠較好地預(yù)測稀釋氣體組分對DME著火延遲期的影響,但是在較高上止點(diǎn)溫度下,尤其對于第一階段著火延遲期,計(jì)算值較實(shí)驗(yàn)值偏高。這主要是由于實(shí)驗(yàn)過程中采用加熱的方法提高混合氣體初始溫度,而第一階段著火延遲期較短,是因?yàn)樵趬嚎s過程中存在潛在的化學(xué)反應(yīng)。

    圖4 稀釋氣體對二甲醚著火延遲期的影響Fig.4 Effects of diluent gas composition on ignition delays of DME

    如圖4所示,不同稀釋氣體組分下,二甲醚混合氣體第一階段著火延遲期基本相同,即稀釋氣體組分對二甲醚混合氣體第一階段著火延遲期影響很小。從圖中還可以看出,二甲醚混合氣體總著火延遲期呈現(xiàn)明顯的NTC現(xiàn)象。在低溫區(qū)間,稀釋氣體組分對二甲醚總著火延遲期的影響同樣很小,但在NTC區(qū)間,稀釋氣體組分對總著火延遲期具有較大的影響,尤其在NTC區(qū)間的起點(diǎn)溫度附近。與N2相比,N2/Ar(50%/50%)使得實(shí)驗(yàn)結(jié)果與模擬結(jié)果總著火延遲期分別縮短30%和15%。

    圖5為N2、N2/Ar(50%/50%)和Ar/CO2(61.2%/38.8%)作為稀釋氣體,上止點(diǎn)壓力為1.5 MPa、上止點(diǎn)溫度為700 K時壓縮及燃燒過程的實(shí)驗(yàn)壓力曲線。不同稀釋氣體組分下,二甲醚混合氣體第一階段著火延遲期基本相同,但是第一階段著火后的壓力和溫度升高差異很大。第一階段著火后的壓力和溫度升高由混合氣體的累計(jì)放熱量和比熱容決定,而比熱容起到主導(dǎo)作用。相同實(shí)驗(yàn)條件下,N2累計(jì)放熱量更高[4],但是Ar具有更小的比熱容,從而第一階段著火后壓力和溫度升高幅度更大。如圖5所示,N2作為稀釋氣體時,第一階段著火后壓力升高0.32 MPa,對應(yīng)溫度升高223 K,而N2/Ar(50%/50%)為稀釋氣體時壓力升高0.42 MPa, 對應(yīng)溫度升高 293 K。由于第二階段著火延遲與第一階段著火相耦合,第一階段著火后壓力和溫度的升高直接影響第二階段著火延遲,進(jìn)而影響總的著火延遲期。與N2相比,N2/Ar(50%/50%)作為稀釋氣體時第一階段著火后的壓力、溫度升高幅度更大,化學(xué)反應(yīng)速率加快,第二階段著火延遲期縮短,進(jìn)而總著火延遲期縮短。在低溫區(qū)間,第一階段著火延遲期較長,而第二階段著火延遲期較短,第一階段著火后溫度和壓力升高對第二階段著火延遲期影響很小,因此總著火延遲期基本相同。而在NTC區(qū)間,由于第二階段著火延遲期較長,第一階段著火后溫度和壓力升高對第二階段著火延遲期影響增加,因此,與N2相比,N2/Ar(50%/50%)使得總著火延遲期出現(xiàn)明顯縮短。

    圖5 稀釋氣體組分對二甲醚著火過程影響的實(shí)驗(yàn)結(jié)果Fig.5 Experiment results of effects of diluent gas composition on ignition process of DME

    已有研究表明[15],對二甲醚著火延遲期和放熱影響較大的2個三體碰撞反應(yīng)為

    表2為ZHAO機(jī)理中幾種稀釋氣體在這2個反應(yīng)中的三體碰撞系數(shù),從表中可以看出,Ar/CO2(61.2%/38.8%)與N2作為稀釋氣體時三體碰撞系數(shù)差異很大,Ar/CO2(61.2%/38.8%)的三體碰撞系數(shù)在2個反應(yīng)中分別為N2的2.09倍和1.87倍。而從圖4中可以看出,Ar/CO2(61.2%/38.8%)作為稀釋氣體時,其第一階段著火延遲期和總著火延遲期與N2幾乎相同,說明在本文實(shí)驗(yàn)溫度范圍內(nèi)由第三體碰撞反應(yīng)引起的化學(xué)效應(yīng)對二甲醚著火延遲期的影響很小。值得注意的是,模擬上止點(diǎn)溫度高于830 K時總著火延遲期有所差異,即稀釋氣體的化學(xué)效應(yīng)開始變得明顯。

    表2 三體反應(yīng)碰撞系數(shù)

    為進(jìn)一步研究稀釋氣體對DME第一階段著火延遲期和總著火延遲期的影響,將N2的比熱容Cp設(shè)置為分別與Ar和CO2相同,在上止點(diǎn)溫度800 K時進(jìn)行數(shù)值模擬。熱效應(yīng)指數(shù)定義為Δτt,Ar(τmodifiedN2-τN2)和Δτt,CO2(τmodifiedN2-τN2),化學(xué)效應(yīng)指數(shù)定義為Δτc,Ar(τAr-τmodifiedN2)和Δτc,CO2(τCO2-τmodifiedN2)。如圖6所示,稀釋氣體的熱效應(yīng)和化學(xué)效應(yīng)對DME第一階段著火延遲期影響均很小。Ar的熱效應(yīng)使得總著火延遲期縮短0.79 ms,而CO2的熱效應(yīng)使得總著火延遲期延長2.02 ms。這主要是由于不同稀釋氣體比熱容的差異引起的,如與N2相比,CO2具有較高的比熱容,第一階段著火后壓力和溫度升高幅度較小,化學(xué)反應(yīng)速率降低,第二階段著火延遲期延長,從而總著火延遲期延長。從圖中還可以看出,Ar的化學(xué)效應(yīng)使得總著火延遲期延長0.14 ms,而CO2的化學(xué)效應(yīng)使得總著火延遲期縮短0.28 ms。這主要是由于不同稀釋氣體三體碰撞系數(shù)的差異引起的,如與N2相比,CO2具有更高的三體碰撞系數(shù)(該系數(shù)是N2的3.8倍),從而化學(xué)反應(yīng)速率加快,著火延遲期縮短。

    圖6 熱效應(yīng)和化學(xué)效應(yīng)對二甲醚著火延遲期影響Fig.6 Thermal and chemical effects on ignition delays of DME

    為進(jìn)一步研究不同溫度工況下稀釋氣體對DME著火延遲期的影響,在650~1 000 K溫度范圍內(nèi)進(jìn)行模擬研究,并將結(jié)果進(jìn)行歸一化處理。歸一化處理公式為

    η=Δτ/τtotal

    (3)

    式中 Δτ——前文中提到的Δτt,Ar、Δτc,Ar、Δτt,CO2和Δτc,CO2

    τtotal——Ar和CO2作為稀釋氣體時的總著火延遲期

    如圖7和圖8所示,稀釋氣體的熱效應(yīng)和化學(xué)效應(yīng)對第一階段著火延遲期影響很小。從圖7可以看出,隨上止點(diǎn)溫度增加,Ar的熱效應(yīng)和化學(xué)效應(yīng)均呈現(xiàn)先增強(qiáng)后減小趨勢,但是在整個溫度區(qū)間,Ar的熱效應(yīng)對DME總著火延遲期的影響起主導(dǎo)作用。從圖8可以看出,與Ar類似,在低溫區(qū)間和NTC區(qū)間CO2的熱效應(yīng)起主導(dǎo)作用。但是,當(dāng)上止點(diǎn)溫度超過NTC區(qū)間時(大于930 K),CO2的化學(xué)效應(yīng)影響增強(qiáng)且超過熱效應(yīng)而起到主導(dǎo)作用,CO2作為稀釋氣體時的總著火延遲期小于N2。

    圖7 Ar的熱效應(yīng)和化學(xué)效應(yīng)對二甲醚著火延遲期的影響Fig.7 Thermal and chemical effects of Ar on ignition delays of DME

    圖8 CO2的熱效應(yīng)和化學(xué)效應(yīng)對二甲醚著火延遲期的影響Fig.8 Thermal and chemical effects of CO2 on ignition delays of DME

    3.2 稀釋氣體的稀釋效應(yīng)

    圖9為N2稀釋率對二甲醚第一階段著火延遲期和總著火延遲期影響的實(shí)驗(yàn)與模擬結(jié)果。由圖9可得,ZHAO機(jī)理能夠較好地預(yù)測N2稀釋率對DME著火延遲的影響,但是在較高上止點(diǎn)溫度下,第一階段著火延遲期計(jì)算值明顯高于實(shí)驗(yàn)值。如圖9所示,隨著N2稀釋率增加,第一階段著火延遲期略有延長,但是總著火延遲期出現(xiàn)明顯延長。

    圖9 N2稀釋率對二甲醚著火延遲期的影響Fig.9 Effects of N2 dilution ratio on ignition delays of DME

    圖10所示為不同N2稀釋率下DME第一階段著火放熱(用壓力升高代表)的實(shí)驗(yàn)與模擬結(jié)果,模擬與實(shí)驗(yàn)結(jié)果趨勢一致,實(shí)驗(yàn)值小于模擬值。如前文所述,第二階段著火延遲與第一階段著火相耦合,第一階段著火產(chǎn)生的活性分子數(shù)量主導(dǎo)第二階段著火過程。隨N2稀釋率增加,DME和O2濃度降低,第一階段著火延遲期略有延長,但如圖10所示,DME第一階段著火放熱明顯降低,累計(jì)活性分子數(shù)量降低,從而第二階段著火延遲期延長。因此,隨N2稀釋率增加總著火延遲期明顯延長。從圖10還可以看出,隨上止點(diǎn)溫度增加,DME第一階段著火放熱近似呈線性變化降低,且不同N2稀釋率下第一階段著火放熱差異縮小。

    圖10 N2稀釋率對第一階段著火壓力升高的影響Fig.10 Effects of N2 dilution ratio on pressure rise during the first-stage ignition

    從圖9還可以看出,DME總著火延遲期出現(xiàn)明顯的NTC現(xiàn)象,且隨N2稀釋率增加,NTC現(xiàn)象變得更加明顯。已有研究結(jié)果表明[20]形成過氧烷基(RO2)和烷基(R)反應(yīng)的比例決定NTC現(xiàn)象。隨N2稀釋率增加,O2濃度降低,更易形成穩(wěn)定分子烷基。因此,隨N2稀釋率增加,NTC區(qū)間向低溫方向移動,NTC現(xiàn)象變得更加明顯。

    N2稀釋率的實(shí)驗(yàn)與模擬結(jié)果對比證明ZHAO機(jī)理能夠較好地預(yù)測稀釋率對DME著火延遲期的影響。因此,通過數(shù)值模擬進(jìn)一步研究了Ar稀釋率和CO2稀釋率對DME著火延遲期的影響。如圖11和圖12所示,Ar稀釋率和CO2稀釋率對DME著火延遲期的影響與N2呈相似的變化趨勢。此外,與N2和Ar相比,CO2作為稀釋氣體時總著火延遲期的NTC現(xiàn)象不明顯。

    圖11 Ar稀釋率對二甲醚著火延遲期的影響Fig.11 Effects of Ar dilution ratio on ignition delays of DME

    圖12 CO2稀釋率對二甲醚著火延遲期的影響Fig.12 Effects of CO2 dilution ratio on ignition delays of DME

    4 結(jié)論

    (1) 實(shí)驗(yàn)與模擬結(jié)果表明:稀釋氣體組分對DME第一階段著火延遲期影響很小,在低溫區(qū)間,稀釋氣體組分對總著火延遲期的影響同樣很小。而在NTC區(qū)間,稀釋氣體組分對總著火延遲期具有較大的影響,尤其在NTC區(qū)間的起點(diǎn)溫度附近。與N2相比,N2/Ar(50%/50%)使得實(shí)驗(yàn)結(jié)果與模擬結(jié)果總的著火延遲期分別縮短30%和15%。

    (2) 對比N2和Ar/CO2(61.2%/38.8%)作為稀釋氣體時的DME著火延遲期,實(shí)驗(yàn)結(jié)果表明稀釋氣體的化學(xué)效應(yīng)對第一階段著火延遲期和總著火延遲期的影響很小,而模擬結(jié)果表明上止點(diǎn)溫度超過830 K時,化學(xué)效應(yīng)開始變得明顯。

    (3) 分別用N2、Ar和CO2作稀釋氣體的模擬結(jié)果表明:在整個溫度區(qū)間范圍內(nèi),Ar的熱效應(yīng)對DME總著火延遲期的影響起主導(dǎo)作用;在低溫區(qū)間和NTC區(qū)間,CO2的熱效應(yīng)同樣起主導(dǎo)作用,但是當(dāng)上止點(diǎn)溫度超過NTC區(qū)間時,CO2的化學(xué)效應(yīng)影響增強(qiáng)且超過熱效應(yīng)而起到主導(dǎo)作用。

    (4) 隨稀釋率增加,第一階段著火延遲期略有延長,而總著火延遲期出現(xiàn)明顯延長,且隨稀釋率增加,總著火延遲期的NTC現(xiàn)象變得更加明顯。

    1 SUNG C J, CURRAN H J. Using rapid compression machines for chemical kinetic studies [J]. Progress in Energy and Combustion Science, 2014, 44: 1-18.

    2 WURMEL J, SILKE E J, CERRAN H J, et al. The effect of diluent gases on ignition delay times in the shock tube and in the rapid compression machine [J]. Combustion and Flame, 2007,151(1-2):289-302.

    3 SHEN H-P S, VANDEROVER J, OEHLSCHLAEGER M A. A shock tube study of iso-octane ignition at elevated pressures: the influence of diluent gases [J].Combustion and Flame, 2008,155(4):739-755.

    4 DI H S, HE X, ZHANG P, et al. Effects of buffer gas composition on low temperature ignition of iso-octane and n-heptane [J]. Combustion and Flame, 2014,161(10):2531-2538.

    5 PETERSEN E L, DAVIDSON D F, HANSON R K. Kinetics modeling of shock-induced ignition in low-dilution CH4/O2mixtures at high pressures and intermediate temperatures [J].Combustion and Flame, 1999,117(1-2):272-290.

    6 LIU H, ZHANG H G, SHI Z C, et al. Performance characterization and auto-ignition performance of a rapid compression machine [J]. Energies, 2014,7(9):6083-6104.

    7 MITTAL G, SUNG C J. Aerodynamics inside a rapid compression machine [J]. Combustion and Flame, 2006,145(1-2):160-180.

    8 SANTE R D. Measurements of the auto-ignition of n-heptane/toluene mixture using a rapid compression machine [J]. Combustion and Flame, 2012,159(1):55-63.

    9 LEE D,HOCHGREB S. Rapid compression machines: heat transfer and suppression of corner vortex [J].Combustion and Flame, 1998,114(3): 531-545.

    10 WERLER M, CANCINO L R, SCHIESSL R, et al. Ignition delay times of diethyl ether measured in a high-pressure shock tube and a rapid compression machine [J]. Proceedings of the Combustion Institute, 2015, 35(1):259-266.

    11 ALLEN C, MITTAL G, SUNG C J, et al. An aerosol rapid compression machine for studying energetic-nanoparticle-enhanced combustion of liquid fuels [J]. Proceedings of the Combustion Institute, 2011, 33:3367-3374.

    12 YU Y, VANHOVE G, GRIFFITHS J F, et al. Influence of EGR and syngas components on the autoignition of natural gas in a rapid compression machine: a detailed experimental study[J]. Energy Fuels, 2013,27(7): 3988-3996.

    13 AN H, CHUNG J, LEE S, et al. The effects of hydrogen addition on the auto-ignition delay of homogeneous primary reference fuel/air mixtures in a rapid compression machine[J]. International Journal of Hydrogen Energy,2015,40(40): 13994-14005.

    14 ZHAO Z, CHAOS M, KAZAKOV A, et al. Thermal decomposition reaction and a comprehensive kinetic model of dimethyl ether [J]. International Journal of Chemical Kinetics, 2008,40(1):1-18.

    15 ULTAN B, KIERAN P S, PETER O’TOOLE, et al. An ignition delay and kinetic modeling study of methane, dimethyl ether, and their mixtures at high pressures [J]. Combustion and Flame,2015,162(2):315-330.

    16 PAN L, HU E J, TIAN Z M, et al. Experimental and kinetic study on ignition delay times of dimethyl ether at high temperatures [J]. Energy Fuels, 2015, 29(5): 3495-3506.

    17 CURRAN H J, FISCHER S L, DRYER F L. The reaction kinetics of dimethyl ether. II: low-temperature oxidation in flow reactors [J]. International Journal of Chemical Kinetic, 2000,32(12):741-759.

    18 MITTAL G, CHAOS M, SUNG C J, et al. Dimethyl ether autoignition in a rapid compression machine: experiments and chemical kinetic modeling [J]. Fuel Processing Technology, 2008,89(12):1244-1254.

    19 張紅光,石智成,盧海濤,等.基于快速壓縮機(jī)的甲烷著火延遲期的研究[J]. 北京工業(yè)大學(xué)學(xué)報(bào),2016,42(4):577-584. ZHANG Hongguang, SHI Zhicheng, LU Haitao, et al. Study of ignition delay time of methane based on a rapid compression machine [J]. Journal of Beijing University of Technology, 2016,42(4):577-584. (in Chinese)

    20 WAGNON S W, WOOLDRIDGE M S. Effects of buffer gas composition on autoignition [J]. Combustion and Flame, 2014, 161(4): 898-907.

    Effects of Diluent Gases on Ignition Delays Characteristics of Dimethyl Ether

    ZHANG Hongguang1,2LI Jiazheng1,2SHI Zhicheng1,2GAO Xiang1,2

    (1.CollegeofEnvironmentalandEnergyEngineering,BeijingUniversityofTechnology,Beijing100124,China2.CollaborativeInnovationCenterofElectricVehiclesinBeijing,Beijing100124,China)

    Effects of diluent gases composition on ignition delays of dimethyl ether (DME) were investigated at compressed pressure of 1.5 MPa, compressed temperature from 670 K to 795 K and equivalence ratio of 1 by using a rapid compression machine (RCM). The diluent gases considered were nitrogen (N2), mixture of N2and argon (Ar) at mole ratio of 50% to 50% and mixture of Ar and carbon dioxide (CO2) at mole ratio of 61.2% to 38.8%. The simulation study was performed over a wider temperature range based on CHMKIN-PRO software. The results showed that diluent gases composition had little impact on the first-stage ignition delay. However, significant differences in the total ignition delay were observed, especially in the negative temperature coefficient (NTC) region. Compared with N2, mixture of N2and Ar decreased the total ignition delay by 30%. The chemical effect of diluent gases composition on the first-stage and total ignition delays was negligible. The thermal effect was dominant factor in the low temperature and NTC regions, however, the chemical effect of CO2enhanced and exceeded the thermal effect in the temperature region beyond the NTC region. With the increase of N2dilution ratio, the first-stage ignition delay was slightly increased, while a significant increase in the total ignition delay was observed. Moreover, the NTC behavior of total ignition delay was noted to become more pronounced at high N2dilution ratio.

    diluent gases; dimethyl ether; ignition delay; rapid compression machine

    10.6041/j.issn.1000-1298.2017.04.044

    2016-08-01

    2016-08-28

    國家自然科學(xué)基金項(xiàng)目(51376011)、北京市自然科學(xué)基金項(xiàng)目(3152005)和北京市教育委員會科技計(jì)劃重點(diǎn)項(xiàng)目(KZ201410005003)

    張紅光(1970—),男,教授,博士生導(dǎo)師,主要從事內(nèi)燃機(jī)燃燒與測試研究,E-mail: zhg5912@263.net

    TK401

    A

    1000-1298(2017)04-0335-07

    猜你喜歡
    二甲醚熱效應(yīng)混合氣體
    SF6/N2混合氣體負(fù)流注放電特性的模擬研究
    吉林電力(2022年2期)2022-11-10 09:24:46
    摻氫對二甲醚層流燃燒特性的影響
    化學(xué)反應(yīng)熱效應(yīng)類試題解析
    NO、O2及NO2混合氣體考點(diǎn)歸納
    混合氣體在聚合物注射成型保壓階段中擴(kuò)散的分子動力學(xué)模擬
    中國塑料(2016年8期)2016-06-27 06:34:58
    加載速率對合成纖維力學(xué)性能影響的力熱效應(yīng)
    基于分子弛豫模型的混合氣體多物理場二維重建算法
    氣相法二甲醚生產(chǎn)腐蝕問題探討
    二甲醚汽車研制:領(lǐng)先者的創(chuàng)新與困境
    二甲醚燃料的發(fā)展前景
    亚洲欧美成人精品一区二区| 日韩三级伦理在线观看| 伦理电影免费视频| 午夜福利视频精品| 成人18禁高潮啪啪吃奶动态图 | 国产一级毛片在线| 嘟嘟电影网在线观看| 少妇的逼水好多| 肉色欧美久久久久久久蜜桃| 亚洲国产精品专区欧美| 亚洲va在线va天堂va国产| 成年av动漫网址| 日韩在线高清观看一区二区三区| 久久99蜜桃精品久久| 国产精品欧美亚洲77777| 国产亚洲精品久久久com| 亚洲精品乱码久久久v下载方式| 99re6热这里在线精品视频| 亚洲av国产av综合av卡| 一边亲一边摸免费视频| 简卡轻食公司| 大又大粗又爽又黄少妇毛片口| 黄色日韩在线| 女人精品久久久久毛片| 好男人视频免费观看在线| 国产精品蜜桃在线观看| 久久女婷五月综合色啪小说| 最新中文字幕久久久久| 插阴视频在线观看视频| 国产日韩欧美在线精品| 久久热精品热| av天堂久久9| 精品人妻熟女av久视频| 熟女av电影| 午夜日本视频在线| 偷拍熟女少妇极品色| 99热全是精品| 99re6热这里在线精品视频| 国产精品偷伦视频观看了| 丰满乱子伦码专区| 精品99又大又爽又粗少妇毛片| 熟女av电影| 国产极品天堂在线| 久久久久久久精品精品| 乱人伦中国视频| 色94色欧美一区二区| 看非洲黑人一级黄片| 嫩草影院新地址| 国产乱来视频区| 日韩欧美一区视频在线观看 | 韩国av在线不卡| 成人午夜精彩视频在线观看| 免费人妻精品一区二区三区视频| 日本爱情动作片www.在线观看| 啦啦啦啦在线视频资源| 国产黄片美女视频| 日韩中文字幕视频在线看片| 高清欧美精品videossex| 又爽又黄a免费视频| 国产黄片美女视频| 国产男人的电影天堂91| 久久人人爽人人片av| 少妇被粗大的猛进出69影院 | 国产在线男女| 街头女战士在线观看网站| 日韩强制内射视频| 欧美丝袜亚洲另类| 欧美+日韩+精品| 欧美日韩在线观看h| 9色porny在线观看| 国产爽快片一区二区三区| 日本黄色日本黄色录像| 亚洲欧洲国产日韩| 2021少妇久久久久久久久久久| 又黄又爽又刺激的免费视频.| 插逼视频在线观看| xxx大片免费视频| 亚洲四区av| 午夜影院在线不卡| 亚洲精品视频女| 我要看日韩黄色一级片| 亚洲国产成人一精品久久久| 精品国产露脸久久av麻豆| 爱豆传媒免费全集在线观看| 日韩 亚洲 欧美在线| 国产在线男女| 看非洲黑人一级黄片| 99久久综合免费| 亚洲av在线观看美女高潮| 成人亚洲欧美一区二区av| 久久久久久伊人网av| 亚洲综合精品二区| 人人妻人人看人人澡| 汤姆久久久久久久影院中文字幕| 一级毛片久久久久久久久女| 99久久人妻综合| 美女脱内裤让男人舔精品视频| 国产欧美另类精品又又久久亚洲欧美| 国产 精品1| 亚洲精品一二三| 中文欧美无线码| 18禁在线播放成人免费| 国产爽快片一区二区三区| 99精国产麻豆久久婷婷| 在线观看三级黄色| 国产精品国产三级国产专区5o| 熟女av电影| 日本午夜av视频| 国产成人a∨麻豆精品| 精品亚洲乱码少妇综合久久| 久久久亚洲精品成人影院| 久久久精品94久久精品| 国产欧美日韩一区二区三区在线 | 国产成人91sexporn| 国语对白做爰xxxⅹ性视频网站| 亚洲第一av免费看| 自拍偷自拍亚洲精品老妇| 欧美日韩视频高清一区二区三区二| 香蕉精品网在线| 精品酒店卫生间| 亚洲综合色惰| 最近2019中文字幕mv第一页| 日韩在线高清观看一区二区三区| 日本vs欧美在线观看视频 | av播播在线观看一区| 日本爱情动作片www.在线观看| 久久精品国产亚洲av天美| 国产av码专区亚洲av| 五月玫瑰六月丁香| 久久精品久久久久久久性| 黄色配什么色好看| xxx大片免费视频| 97精品久久久久久久久久精品| 在线播放无遮挡| 婷婷色综合大香蕉| 热99国产精品久久久久久7| 精品一品国产午夜福利视频| 日本午夜av视频| 久久午夜福利片| 亚洲三级黄色毛片| 欧美日韩精品成人综合77777| 国产免费一区二区三区四区乱码| 久久久久久久亚洲中文字幕| 成人亚洲欧美一区二区av| 这个男人来自地球电影免费观看 | 久热久热在线精品观看| 久久久a久久爽久久v久久| 国产在线视频一区二区| 亚洲av中文av极速乱| 午夜激情久久久久久久| 最后的刺客免费高清国语| 丝袜在线中文字幕| 赤兔流量卡办理| 三级国产精品欧美在线观看| 99久久精品热视频| 在线观看国产h片| 亚洲国产色片| 97在线人人人人妻| 国产一区二区三区综合在线观看 | 欧美激情极品国产一区二区三区 | 搡女人真爽免费视频火全软件| 最近手机中文字幕大全| 一个人看视频在线观看www免费| 少妇丰满av| 大香蕉97超碰在线| 久久精品国产亚洲av涩爱| 九九久久精品国产亚洲av麻豆| 国产淫片久久久久久久久| 男人和女人高潮做爰伦理| 国产成人精品一,二区| 国产精品国产三级专区第一集| av国产精品久久久久影院| 精品人妻一区二区三区麻豆| 国产精品国产三级专区第一集| 91精品国产国语对白视频| 99热这里只有是精品50| 欧美bdsm另类| 免费在线观看成人毛片| 国产有黄有色有爽视频| 一级二级三级毛片免费看| 久久精品国产鲁丝片午夜精品| 午夜福利网站1000一区二区三区| 少妇裸体淫交视频免费看高清| 十八禁高潮呻吟视频 | 韩国高清视频一区二区三区| 日本黄大片高清| 亚洲情色 制服丝袜| 内地一区二区视频在线| 最新的欧美精品一区二区| 男男h啪啪无遮挡| 大又大粗又爽又黄少妇毛片口| 国产精品无大码| 91久久精品国产一区二区成人| 国产高清三级在线| 观看av在线不卡| 男人爽女人下面视频在线观看| 久久久久久久久久成人| 亚州av有码| 亚洲国产色片| 欧美日韩综合久久久久久| 精品久久国产蜜桃| 一个人免费看片子| 免费人妻精品一区二区三区视频| 日本欧美国产在线视频| 成人免费观看视频高清| 激情五月婷婷亚洲| 欧美精品高潮呻吟av久久| av在线观看视频网站免费| 亚洲婷婷狠狠爱综合网| 国产欧美日韩精品一区二区| 亚洲国产精品专区欧美| 国产精品久久久久久精品古装| 最近中文字幕2019免费版| 亚洲欧洲日产国产| 纯流量卡能插随身wifi吗| 日本免费在线观看一区| 噜噜噜噜噜久久久久久91| 成年av动漫网址| 精品熟女少妇av免费看| 亚洲欧美日韩卡通动漫| 蜜桃在线观看..| 精华霜和精华液先用哪个| 男人爽女人下面视频在线观看| 欧美激情极品国产一区二区三区 | 色婷婷av一区二区三区视频| 亚洲精品中文字幕在线视频 | 六月丁香七月| 国产成人freesex在线| 亚洲自偷自拍三级| 精品人妻一区二区三区麻豆| 久久精品国产自在天天线| 水蜜桃什么品种好| 国产免费福利视频在线观看| 久久午夜福利片| 中国美白少妇内射xxxbb| 最近最新中文字幕免费大全7| 日本黄色片子视频| 男女啪啪激烈高潮av片| 国产精品熟女久久久久浪| www.av在线官网国产| 国产亚洲午夜精品一区二区久久| 亚洲国产精品国产精品| 99热国产这里只有精品6| 欧美日韩av久久| 亚洲经典国产精华液单| 少妇裸体淫交视频免费看高清| 亚洲三级黄色毛片| 一本色道久久久久久精品综合| 美女视频免费永久观看网站| 自拍偷自拍亚洲精品老妇| 国内精品宾馆在线| 国产熟女欧美一区二区| 91精品国产国语对白视频| 国产精品人妻久久久久久| 日韩大片免费观看网站| 亚洲精品乱码久久久v下载方式| 亚洲精品久久久久久婷婷小说| 午夜福利影视在线免费观看| 精品人妻熟女毛片av久久网站| 久久久国产一区二区| av线在线观看网站| 看非洲黑人一级黄片| 亚洲精品第二区| 中文字幕人妻丝袜制服| 日韩欧美精品免费久久| 久久久久国产网址| 亚洲国产精品一区三区| 久久久欧美国产精品| xxx大片免费视频| 欧美成人午夜免费资源| 日韩伦理黄色片| 亚洲人与动物交配视频| 在线看a的网站| 亚洲色图综合在线观看| 在线观看一区二区三区激情| 国产伦精品一区二区三区视频9| 亚洲精品国产av成人精品| 亚洲性久久影院| 亚洲欧洲国产日韩| 日本黄色片子视频| 国产永久视频网站| 国语对白做爰xxxⅹ性视频网站| 亚洲人成网站在线观看播放| 亚洲伊人久久精品综合| 超碰97精品在线观看| 久久99一区二区三区| 搡老乐熟女国产| 精品亚洲乱码少妇综合久久| 最新中文字幕久久久久| 亚洲欧美日韩另类电影网站| 一级毛片黄色毛片免费观看视频| 免费观看性生交大片5| 赤兔流量卡办理| 女人精品久久久久毛片| 狠狠精品人妻久久久久久综合| 这个男人来自地球电影免费观看 | 一级av片app| 高清午夜精品一区二区三区| 最近最新中文字幕免费大全7| 亚洲高清免费不卡视频| 麻豆乱淫一区二区| videos熟女内射| 97超碰精品成人国产| 国产亚洲5aaaaa淫片| 久久久久视频综合| 男男h啪啪无遮挡| 下体分泌物呈黄色| 色94色欧美一区二区| 亚洲情色 制服丝袜| 成年女人在线观看亚洲视频| 一本色道久久久久久精品综合| 国产亚洲一区二区精品| 国产黄色免费在线视频| 丝袜脚勾引网站| av天堂中文字幕网| 亚洲成人手机| av福利片在线观看| 成人无遮挡网站| 亚洲,欧美,日韩| av天堂久久9| 成人午夜精彩视频在线观看| 国产男人的电影天堂91| 欧美成人午夜免费资源| 国产真实伦视频高清在线观看| 丝袜脚勾引网站| 久久久久久久久久人人人人人人| 免费观看无遮挡的男女| 伦精品一区二区三区| 国产成人精品婷婷| 黄色欧美视频在线观看| 欧美精品国产亚洲| 国产色婷婷99| 亚洲美女视频黄频| 麻豆乱淫一区二区| 国内少妇人妻偷人精品xxx网站| 99久久综合免费| 国产伦在线观看视频一区| 丰满乱子伦码专区| 国产欧美日韩精品一区二区| 777米奇影视久久| 自线自在国产av| 亚洲真实伦在线观看| 性色avwww在线观看| 国产精品国产av在线观看| av又黄又爽大尺度在线免费看| 亚洲精品国产av蜜桃| 免费人妻精品一区二区三区视频| 免费黄色在线免费观看| 久久久久久久久久久久大奶| 国产片特级美女逼逼视频| 在线天堂最新版资源| 亚洲婷婷狠狠爱综合网| av福利片在线观看| 日韩免费高清中文字幕av| 欧美区成人在线视频| 99久久综合免费| 女人久久www免费人成看片| 国精品久久久久久国模美| 免费看光身美女| 日韩成人av中文字幕在线观看| 欧美人与善性xxx| 人人澡人人妻人| 91精品一卡2卡3卡4卡| 国产亚洲最大av| 人人妻人人看人人澡| 亚洲精品国产成人久久av| av视频免费观看在线观看| 亚洲内射少妇av| 在线观看免费视频网站a站| 菩萨蛮人人尽说江南好唐韦庄| 精品99又大又爽又粗少妇毛片| 3wmmmm亚洲av在线观看| 国产无遮挡羞羞视频在线观看| 在线观看免费日韩欧美大片 | 国产成人午夜福利电影在线观看| 国产精品女同一区二区软件| 99国产精品免费福利视频| 免费大片18禁| 黄色一级大片看看| av国产久精品久网站免费入址| 大陆偷拍与自拍| 国语对白做爰xxxⅹ性视频网站| 人体艺术视频欧美日本| 国产精品熟女久久久久浪| 久久99一区二区三区| 男人爽女人下面视频在线观看| 纵有疾风起免费观看全集完整版| 免费黄色在线免费观看| 国产一区二区三区av在线| 国产 精品1| 日日摸夜夜添夜夜添av毛片| av在线老鸭窝| 国产女主播在线喷水免费视频网站| 久久精品久久精品一区二区三区| 欧美一级a爱片免费观看看| 成人国产麻豆网| 国产日韩欧美视频二区| 99精国产麻豆久久婷婷| 狂野欧美激情性xxxx在线观看| 老熟女久久久| 自拍偷自拍亚洲精品老妇| 麻豆成人午夜福利视频| 国产探花极品一区二区| 午夜激情久久久久久久| 亚洲精品国产av蜜桃| 亚洲精品视频女| 国产精品久久久久成人av| 国产在线免费精品| 一级毛片电影观看| 老女人水多毛片| 婷婷色综合大香蕉| 国产乱来视频区| 永久免费av网站大全| 最后的刺客免费高清国语| 亚洲在久久综合| 欧美国产精品一级二级三级 | 国产黄色免费在线视频| 欧美 亚洲 国产 日韩一| 久久久久久久久久久免费av| 美女xxoo啪啪120秒动态图| 在线观看免费视频网站a站| 亚洲综合色惰| 丁香六月天网| 国产亚洲午夜精品一区二区久久| 高清视频免费观看一区二区| 国产成人a∨麻豆精品| 人体艺术视频欧美日本| 亚洲国产精品成人久久小说| 国产极品天堂在线| 亚洲精品国产av蜜桃| 久久久久久伊人网av| 国产色婷婷99| 欧美性感艳星| 色94色欧美一区二区| 美女cb高潮喷水在线观看| 全区人妻精品视频| 乱系列少妇在线播放| 在线观看免费日韩欧美大片 | 女人久久www免费人成看片| av国产久精品久网站免费入址| 午夜福利在线观看免费完整高清在| 男人爽女人下面视频在线观看| 亚洲av中文av极速乱| 美女视频免费永久观看网站| 日日撸夜夜添| a级一级毛片免费在线观看| 亚洲情色 制服丝袜| 国产精品.久久久| 一本久久精品| 天天躁夜夜躁狠狠久久av| 欧美丝袜亚洲另类| www.色视频.com| 男人添女人高潮全过程视频| 亚洲国产精品一区三区| .国产精品久久| 色婷婷久久久亚洲欧美| 亚洲精品自拍成人| 亚洲国产成人一精品久久久| 免费久久久久久久精品成人欧美视频 | 99re6热这里在线精品视频| 午夜激情久久久久久久| 国产美女午夜福利| 欧美3d第一页| 国产精品99久久久久久久久| 国产黄色免费在线视频| 人人澡人人妻人| 午夜免费观看性视频| 在线观看三级黄色| 久久久国产欧美日韩av| 久久99蜜桃精品久久| 国产中年淑女户外野战色| 日韩一本色道免费dvd| 国产片特级美女逼逼视频| 久久人妻熟女aⅴ| 午夜老司机福利剧场| av福利片在线观看| 精品卡一卡二卡四卡免费| 黄色配什么色好看| 啦啦啦中文免费视频观看日本| 蜜桃在线观看..| 丰满乱子伦码专区| 久久久精品免费免费高清| 高清av免费在线| 亚洲三级黄色毛片| 午夜激情久久久久久久| 青春草视频在线免费观看| 久久人人爽人人爽人人片va| 国产成人精品一,二区| 日韩三级伦理在线观看| av在线老鸭窝| 18禁裸乳无遮挡动漫免费视频| 老女人水多毛片| 99久久中文字幕三级久久日本| 日日啪夜夜爽| 国产伦精品一区二区三区视频9| 国产精品久久久久久精品古装| freevideosex欧美| 夜夜骑夜夜射夜夜干| 日日啪夜夜撸| 在线 av 中文字幕| 精品人妻熟女毛片av久久网站| 国产精品三级大全| 欧美一级a爱片免费观看看| xxx大片免费视频| 精品亚洲成a人片在线观看| 精品国产乱码久久久久久小说| 97精品久久久久久久久久精品| 纵有疾风起免费观看全集完整版| 五月天丁香电影| 纵有疾风起免费观看全集完整版| 中文乱码字字幕精品一区二区三区| 欧美xxⅹ黑人| 亚洲精品日本国产第一区| 亚洲精品aⅴ在线观看| 中文乱码字字幕精品一区二区三区| 欧美精品人与动牲交sv欧美| 国产中年淑女户外野战色| 国产 一区精品| 蜜桃久久精品国产亚洲av| 又大又黄又爽视频免费| 久热久热在线精品观看| 精品久久久久久久久亚洲| 亚洲av不卡在线观看| 3wmmmm亚洲av在线观看| 三上悠亚av全集在线观看 | 美女cb高潮喷水在线观看| 日本-黄色视频高清免费观看| 国内精品宾馆在线| 高清黄色对白视频在线免费看 | 国产有黄有色有爽视频| 色婷婷久久久亚洲欧美| 精华霜和精华液先用哪个| 夜夜爽夜夜爽视频| 91在线精品国自产拍蜜月| 制服丝袜香蕉在线| 91精品国产九色| 日韩中文字幕视频在线看片| 国产探花极品一区二区| 在线 av 中文字幕| 少妇丰满av| 久久久精品免费免费高清| 男男h啪啪无遮挡| 91午夜精品亚洲一区二区三区| 国产深夜福利视频在线观看| 伦理电影免费视频| 九九在线视频观看精品| 妹子高潮喷水视频| 欧美精品国产亚洲| 欧美激情极品国产一区二区三区 | 最近手机中文字幕大全| 欧美日韩在线观看h| 精品久久久噜噜| 最近中文字幕高清免费大全6| 国产成人精品无人区| 中文字幕免费在线视频6| 高清av免费在线| 成人免费观看视频高清| 亚洲欧美日韩另类电影网站| av免费在线看不卡| 精品少妇黑人巨大在线播放| 国产老妇伦熟女老妇高清| 日本午夜av视频| 秋霞伦理黄片| 久久韩国三级中文字幕| a级毛片在线看网站| 久久精品国产鲁丝片午夜精品| 青春草亚洲视频在线观看| 免费观看在线日韩| 18禁动态无遮挡网站| 国国产精品蜜臀av免费| av天堂中文字幕网| 国产69精品久久久久777片| 亚洲精品色激情综合| 免费大片黄手机在线观看| 国产美女午夜福利| 国产成人免费观看mmmm| 一个人看视频在线观看www免费| 男人爽女人下面视频在线观看| 人人妻人人澡人人看| 赤兔流量卡办理| 亚洲va在线va天堂va国产| 如日韩欧美国产精品一区二区三区 | 性色av一级| 亚洲av福利一区| 久久久国产精品麻豆| 少妇熟女欧美另类| av又黄又爽大尺度在线免费看| 久久精品国产亚洲av涩爱| 99九九在线精品视频 | 一本大道久久a久久精品| 大话2 男鬼变身卡| 国产精品一区www在线观看| 国内少妇人妻偷人精品xxx网站| 最近的中文字幕免费完整| 日日啪夜夜撸| 啦啦啦啦在线视频资源| 草草在线视频免费看| 成人美女网站在线观看视频| 国产视频首页在线观看| 亚洲国产日韩一区二区| av卡一久久| 日本-黄色视频高清免费观看| 在线免费观看不下载黄p国产| 观看av在线不卡| 一个人免费看片子| 一级a做视频免费观看| 十八禁网站网址无遮挡 | 精品国产露脸久久av麻豆| 熟女人妻精品中文字幕| 欧美日韩在线观看h| 97在线人人人人妻| 国产成人精品婷婷| 99国产精品免费福利视频| 最后的刺客免费高清国语| 久久免费观看电影|