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

    RP-3航空煤油替代燃料及其化學(xué)反應(yīng)動(dòng)力學(xué)模型

    2015-01-04 12:52:03于維銘鐘北京清華大學(xué)航天航空學(xué)院北京00084北京汽車集團(tuán)產(chǎn)業(yè)投資有限公司北京000
    物理化學(xué)學(xué)報(bào) 2015年4期
    關(guān)鍵詞:傳播速度煤油機(jī)理

    鄭 東 于維銘, 鐘北京,*(清華大學(xué)航天航空學(xué)院,北京00084;北京汽車集團(tuán)產(chǎn)業(yè)投資有限公司,北京000)

    [Article]

    RP-3航空煤油替代燃料及其化學(xué)反應(yīng)動(dòng)力學(xué)模型

    鄭 東1于維銘1,2鐘北京1,*
    (1清華大學(xué)航天航空學(xué)院,北京100084;2北京汽車集團(tuán)產(chǎn)業(yè)投資有限公司,北京100021)

    本文提出了40%(摩爾分?jǐn)?shù),下同)正癸烷、42%正十二烷、13%乙基環(huán)己烷和5%對(duì)二甲苯的四組分RP-3航空煤油替代燃料模型,并通過(guò)實(shí)驗(yàn)充分驗(yàn)證了替代燃料模型與實(shí)際RP-3航空煤油在理化特性上的相似性.采用對(duì)沖火焰實(shí)驗(yàn)臺(tái)架,測(cè)量了RP-3航空煤油以及四組分替代燃料的層流火焰?zhèn)鞑ニ俣?對(duì)比結(jié)果表明本文提出的替代燃料能夠準(zhǔn)確描述實(shí)際RP-3航空煤油的燃燒速率.進(jìn)一步發(fā)展了包含168組分、1089反應(yīng)的半詳細(xì)反應(yīng)動(dòng)力學(xué)模型,驗(yàn)證結(jié)果表明本文機(jī)理能夠準(zhǔn)確預(yù)測(cè)RP-3航空煤油著火延遲時(shí)間和火焰?zhèn)鞑ニ俣?

    RP-3航空煤油;替代燃燒;反應(yīng)機(jī)理;火焰?zhèn)鞑ニ俣?對(duì)沖火焰

    1 引言

    航空發(fā)動(dòng)機(jī)燃燒室中燃油的著火和燃燒過(guò)程是一個(gè)非常復(fù)雜的湍流化學(xué)反應(yīng)過(guò)程.1航空發(fā)動(dòng)機(jī)存在燃燒不穩(wěn)定性,以及高密度、高稀釋度、超稀薄、高馬赫數(shù)等極端條件下點(diǎn)火困難等問(wèn)題,都與燃燒室內(nèi)復(fù)雜的湍流燃燒過(guò)程密切相關(guān).湍流燃燒問(wèn)題的核心是湍流與化學(xué)反應(yīng)的相互作用.要研究湍流與化學(xué)反應(yīng)相互作用,燃料燃燒化學(xué)反應(yīng)機(jī)理是基礎(chǔ).因此發(fā)展航空煤油化學(xué)反應(yīng)動(dòng)力學(xué)模型對(duì)于深入理解航空發(fā)動(dòng)機(jī)湍流燃燒機(jī)制,闡明湍流與化學(xué)反應(yīng)的相互作用機(jī)理具有極其重要的理論和實(shí)際意義.

    ?Editorial office ofActa Physico-Chimica Sinica

    航空煤油是C7-C16大分子碳?xì)鋸?fù)雜烴類的混合物,包含鏈烷烴(直鏈烷烴與支鏈烷烴)、環(huán)烷烴以及芳香烴等成百上千種組分.2由于原油產(chǎn)地、加工工藝以及添加劑的不同,各國(guó)、各型號(hào)、甚至各批次的航空煤油的組分也不盡相同.圖1給出了歐美常用的航空煤油JP-8、Jet-A,以及國(guó)產(chǎn)航空煤油RP-3不同烴類的構(gòu)成比例.2-4由圖可見(jiàn),國(guó)內(nèi)外不同型號(hào)航空煤油的烴類構(gòu)成比例存在較大差異,其具體的組分構(gòu)成差別就更明顯了.鑒于航空煤油成分的復(fù)雜性,無(wú)法針對(duì)某種航空煤油直接構(gòu)建其化學(xué)反應(yīng)動(dòng)力學(xué)模型.行之有效的方法是從復(fù)雜的航空煤油中選取一種或幾種具有代表性的成分作為航空煤油替代燃料的組分,進(jìn)而構(gòu)建其詳細(xì)替代燃料的化學(xué)反應(yīng)機(jī)理.這也是研究汽油、柴油等大分子碳?xì)鋸?fù)雜燃料化學(xué)反應(yīng)機(jī)理的慣用做法.5-7表1列舉了目前部分常見(jiàn)的國(guó)內(nèi)外航空煤油替代燃料及其化學(xué)反應(yīng)動(dòng)力學(xué)模型.8-14

    圖1 JP-8,Jet-A和RP-3烴類組成Fig.1 Molecularclass compositions of JP-8,Jet-A,and RP-3

    表1 航空煤油替代燃料及其化學(xué)反應(yīng)動(dòng)力學(xué)模型Table 1 Aviation kerosene surrogate fuels and chemical reaction kinetic models

    單組分的正癸烷或者正十二烷是最簡(jiǎn)單的航空煤油替代燃料,Dagaut等8在攪拌反應(yīng)器(JSR)中研究了Jet-A和正癸烷的氧化特性,發(fā)展了包含90組分和573反應(yīng)的正癸烷反應(yīng)機(jī)理作為Jet-A替代燃料的反應(yīng)機(jī)理.隨后許多研究者提出了包含鏈烷烴、環(huán)烷烴、芳香烴的多組分航空煤油Jet-A和JP-8替代燃料,并發(fā)展了相應(yīng)的詳細(xì)和簡(jiǎn)化反應(yīng)動(dòng)力學(xué)模型,其規(guī)模從上百組分到幾千組分不等.9-12國(guó)內(nèi)對(duì)于多組分RP-3航空煤油替代燃料及其化學(xué)反應(yīng)機(jī)理的研究還比較少.肖保國(guó)等13提出了一個(gè)79% (摩爾分?jǐn)?shù),下同)正癸烷、13%三甲基環(huán)己烷和8%乙基苯的三組分RP-3航空煤油替代燃料模型(近似分子式C9.71H20.52),并發(fā)展了包含109組分和946反應(yīng)的動(dòng)力學(xué)模型.曾文等14在Dagaut三組分Jet-A替代燃料模型15的基礎(chǔ)上,提出由65%(體積分?jǐn)?shù))正癸烷、10%甲苯、25%丙基環(huán)己烷組成的三組分RP-3航空煤油替代燃料模型,相應(yīng)地發(fā)展了150組分和591反應(yīng)的動(dòng)力學(xué)模型.

    綜合分析國(guó)內(nèi)外研究現(xiàn)狀不難發(fā)現(xiàn),國(guó)產(chǎn)RP-3航空煤油替代燃料及其反應(yīng)機(jī)理的研究與國(guó)外存在較大差距.國(guó)內(nèi)鮮有的幾個(gè)多組分替代模型均是基于物理特性或是參考國(guó)外模型發(fā)展的,其平均碳原子數(shù)與實(shí)際國(guó)產(chǎn)航空煤油相差較大,而且其反應(yīng)機(jī)理的驗(yàn)證也不夠充分.能否描述國(guó)產(chǎn)RP-3航空煤油的化學(xué)反應(yīng)動(dòng)力學(xué)特性尚需進(jìn)一步研究.因此,本文在考慮RP-3航空煤油組成成分、物理特性和化學(xué)特性的基礎(chǔ)上,提出一個(gè)能夠描述實(shí)際RP-3理化特性的多組分替代燃料,通過(guò)測(cè)量火焰?zhèn)鞑ニ俣冗M(jìn)行驗(yàn)證.進(jìn)而發(fā)展相應(yīng)的半詳細(xì)反應(yīng)動(dòng)力學(xué)模型,并驗(yàn)證其著火和燃燒特性.

    2 RP-3航空煤油替代燃料模型

    2.1 RP-3航空煤油成分和物理特性分析

    Edwards和Maurice16建議預(yù)測(cè)實(shí)際燃油物理特性和著火、燃燒、污染物生成等化學(xué)特性的替代燃料應(yīng)采用多組分,而且其組分要與實(shí)際燃油的重要成分相匹配.我們首先對(duì)RP-3航空煤油的組成成分進(jìn)行詳細(xì)分析.RP-3航空煤油樣品來(lái)自于中石化北京燕山石油化工分公司,通過(guò)色譜-質(zhì)譜(GC-MS)聯(lián)用儀確定樣品中含量最大的前100種組分及其相對(duì)體積分?jǐn)?shù).根據(jù)GC-MS分析結(jié)果,國(guó)產(chǎn)航空煤油RP-3的不同烴類體積分?jǐn)?shù)如圖2所示.其中比例最高的是鏈烷烴為53.0%,平均碳原子數(shù)為10.8.其次是環(huán)烷烴為37.7%,平均碳原子數(shù)為10.2.而芳香烴含量只有4.6%,低于國(guó)外JP-8和Jet-A航空煤油.

    圖2 GC-MS分析RP-3的成分組成(體積分?jǐn)?shù))Fig.2 Molecular class compositions(volume fraction)of RP-3 as identified by gas chromatography-mass spectrometer(GC-MS)analysis

    表2 RP-3物理特性分析結(jié)果Table 2 Analysis results for RP-3 physical properties

    然后對(duì)RP-3的一些基本的物理特性進(jìn)行了分析.表2給出了本文測(cè)量的結(jié)果和文獻(xiàn)17中給出的一些物理特性.由表2可見(jiàn),本文實(shí)驗(yàn)測(cè)量結(jié)果與文獻(xiàn)中的結(jié)果較為一致,但也存在一定差異.除了測(cè)量誤差外,主要的原因可能是實(shí)驗(yàn)中所采用的RP-3樣品來(lái)自不同產(chǎn)品批次.

    2.2 RP-3航空煤油替代燃料模型

    發(fā)展一個(gè)多組分RP-3航空煤油替代燃料模型,首先要確定其代表性組分.本文選取代表性組分的兩個(gè)原則是:其一,代表性組分在實(shí)際航空煤油中含量較高,或者能夠綜合反映實(shí)際航空煤油中這種烴類所有組分的平均特性.其二,代表性組分的理化特性與實(shí)際航空煤油接近,或者通過(guò)與其它烴類代表性組分調(diào)配所得到的混合燃料的理化特性與實(shí)際航空煤油接近.首先,由上述成分分析可知, RP-3航空煤油中鏈烷烴比例最高,對(duì)航空煤油理化特性的影響也很大.因此我們選擇正癸烷和正十二烷兩種成分作為RP-3航空煤油鏈烷烴的代表性組分.這樣的選擇一方面是因?yàn)槭樘荚訑?shù)較高,容易與正癸烷調(diào)配出與實(shí)際燃油中比較接近的鏈烷烴平均碳原子數(shù);另一方面是因?yàn)檎榈恼扯群驼舭l(fā)溫度與實(shí)際航空煤油相匹配(見(jiàn)表3).其次,考慮到乙基環(huán)己烷的碳原子數(shù)比甲基環(huán)己烷多,而且乙基環(huán)己烷的密度與RP-3很接近(見(jiàn)表3),因此選擇乙基環(huán)己烷作為RP-3航空煤油環(huán)烷烴的代表性組分.雖然芳香烴類物質(zhì)在RP-3中所占比例不高,但其對(duì)實(shí)際航空煤油的著火與污染物排放有較大影響.相對(duì)于甲苯而言對(duì)二甲苯的碳原子數(shù)較多,且其密度較大(見(jiàn)表3).因此選擇對(duì)二甲苯作為RP-3航空煤油芳香烴的代表性組分.對(duì)二甲苯的加入可以使得替代燃料的密度與RP-3航空煤油密度更接近.由于烯烴和醇類物質(zhì)在RP-3航空煤油的含量較少,現(xiàn)有的航空煤油替代燃料均舍棄這些成分(見(jiàn)表1),因此本文也未考慮這兩種成分.最終本文選擇正癸烷、正十二烷、乙基環(huán)己烷和對(duì)二甲苯四種組分作為RP-3航空煤油的代表性組分,其特性參數(shù)如表3所示.

    表3 RP-3代表性組分特性參數(shù)Table 3 Property parameters of the representative species of RP-3

    確定了RP-3航空煤油的代表性組分后,我們進(jìn)一步需要確定四種代表性組分的比例.方法是:選取實(shí)際燃料重要的理化特性參數(shù),通過(guò)優(yōu)化各代表性組分的比例,使得替代燃料的理化特性參數(shù)值與對(duì)應(yīng)所選取的實(shí)際燃料的特性參數(shù)值相匹配,從而確定替代燃料各組分的比例.由于氣相燃料的擴(kuò)散特性取決于其分子量,18為了能夠再現(xiàn)實(shí)際航空煤油在氣相燃燒中的擴(kuò)散特性,替代燃料的平均分子量應(yīng)該接近實(shí)際燃料的平均分子量;碳?xì)浔确从沉巳紵晌镏蠧O2和H2O的比例,決定了燃料反應(yīng)生成焓以及絕熱火焰溫度.而絕熱火焰溫度與燃燒速率以及其它燃燒現(xiàn)象密切相關(guān).4換而言之,碳?xì)浔仁欠磻?yīng)燃料燃燒速率以及其它燃燒現(xiàn)象的重要特性參數(shù);在兩相流燃燒中,密度和粘度對(duì)于燃料的霧化、蒸發(fā)過(guò)程有較大影響,為了能夠描述實(shí)際航空煤油在兩相流燃燒中燃料的霧化、蒸發(fā)過(guò)程,替代燃料的密度和粘度應(yīng)該與實(shí)際燃料相匹配.因此選擇平均分子量、碳?xì)浔?、密度以及粘度為所要匹配的理化特性參?shù).采用試驗(yàn)和計(jì)算相結(jié)合,經(jīng)過(guò)多次試驗(yàn),最終確定了正癸烷、正十二烷、乙基環(huán)己烷和對(duì)二甲苯四種組分的摩爾分?jǐn)?shù)分別為40%、42%、13%和5%.

    將本文提出的摩爾分?jǐn)?shù)分別為40%正癸烷、42%正十二烷、13%乙基環(huán)己烷和5%對(duì)二甲苯的四組分RP-3航空煤油替代燃料的理化特性進(jìn)行了測(cè)量和計(jì)算,并與實(shí)際RP-3航空煤油的實(shí)驗(yàn)結(jié)果進(jìn)行比較,結(jié)果如表2所示.由表2可見(jiàn),本文提出的四組分RP-3航空煤油替代燃料的平均分子式、平均分子量、碳?xì)浔取⒚芏?、粘度以及蒸發(fā)溫度均與實(shí)際RP-3航空煤油相近.

    3 航空煤油及其替代燃料火焰?zhèn)鞑ニ俣?/h2>

    在上一節(jié)中,我們提出了四組分RP-3航空煤油替代燃料,并通過(guò)對(duì)比驗(yàn)證了部分理化特性參數(shù)與實(shí)際RP-3航空煤油的相似性.本節(jié)我們將進(jìn)一步驗(yàn)證替代燃料和實(shí)際RP-3航空煤油燃燒特性的相似性.眾所周知,火焰?zhèn)鞑ニ俣仁潜碚魅剂先紵匦缘闹匾獏?shù),因此,我們將通過(guò)對(duì)比替代燃料和實(shí)際RP-3航空煤油的火焰?zhèn)鞑ニ俣葋?lái)驗(yàn)證它們?nèi)紵匦缘南嗨菩?

    四組分RP-3航空煤油替代燃料和實(shí)際RP-3航空煤油與空氣預(yù)混氣的層流火焰?zhèn)鞑ニ俣炔捎脤?duì)沖火焰進(jìn)行測(cè)量.圖3所示為液體對(duì)沖火焰實(shí)驗(yàn)臺(tái)系統(tǒng)的示意圖,主要包括燃料輸送系統(tǒng)、空氣輸送系統(tǒng)、空氣加熱裝置、燃料與高溫空氣混合裝置、混合氣輸送系統(tǒng)、示蹤粒子裝置、噴嘴系統(tǒng)、測(cè)量系統(tǒng)等.空氣經(jīng)過(guò)流量計(jì),進(jìn)入預(yù)熱器加熱至預(yù)設(shè)溫度,而后送至混合裝置.與此同時(shí)燃料經(jīng)由注射泵,送入混合裝置.燃料在混合裝置中迅速蒸發(fā),并與空氣充分混合,形成均勻的預(yù)混氣.預(yù)混氣通過(guò)配氣管路,均勻送至上下兩個(gè)對(duì)稱布置的噴嘴.在離噴嘴一定距離處著火燃燒,形成穩(wěn)定的對(duì)沖火焰.通過(guò)三維粒子動(dòng)態(tài)分析儀,測(cè)得對(duì)沖火焰中心軸線的速度分布,進(jìn)而獲得某一拉伸率下的層流火焰?zhèn)鞑ニ俣?然后通過(guò)線性外推,即可得到無(wú)拉伸層流火焰?zhèn)鞑ニ俣?19

    圖3 液體燃料對(duì)沖火焰實(shí)驗(yàn)系統(tǒng)圖Fig.3 System diagram of counterflow twin-flame configuration for liquid fuels

    實(shí)驗(yàn)條件如下:0.1 MPa,燃料/空氣預(yù)混氣初溫設(shè)定為403 K,燃料流量范圍為67.9-238 mL·h-1,空氣流速范圍為0.88-1.44 m3·h-1,化學(xué)計(jì)量比范圍為0.7-1.4.實(shí)驗(yàn)操作、數(shù)據(jù)處理以及誤差分析方法詳情參考文獻(xiàn).19圖4給出了實(shí)驗(yàn)測(cè)量的四組分RP-3航空煤油替代燃料/空氣和實(shí)際RP-3航空煤油/空氣預(yù)混氣在不同化學(xué)計(jì)量比下的層流火焰?zhèn)鞑ニ俣?由圖4可知,在工況相同的情況下,在化學(xué)計(jì)量比0.9-1.3的范圍內(nèi),替代燃料和RP-3煤油的火焰?zhèn)鞑ニ俣炔顒e不大,化學(xué)計(jì)量比小于0.9和大于1.3時(shí),替代燃料的火焰?zhèn)鞑ニ俣嚷孕∮赗P-3煤油.考慮到實(shí)驗(yàn)上的誤差,可以認(rèn)為在預(yù)混氣相同壓力和初始溫度下,四組分替代燃料和RP-3煤油具有相同的火焰?zhèn)鞑ニ俣?說(shuō)明兩者的燃燒特性非常相似.

    圖4 替代燃料和實(shí)際RP-3火焰?zhèn)鞑ニ俣葴y(cè)量結(jié)果Fig.4 Experimental laminar flame speeds of surrogate fuels and real RP-3

    4 替代燃料詳細(xì)化學(xué)反應(yīng)機(jī)理及其驗(yàn)證

    4.1 RP-3替代燃料詳細(xì)反應(yīng)機(jī)理構(gòu)建

    我們依據(jù)上述提出的四組分RP-3航空煤油替代燃料模型構(gòu)建其半詳細(xì)的反應(yīng)動(dòng)力學(xué)模型.由于替代燃料模型中代表性組分均為大分子碳?xì)淙剂?其C3以下反應(yīng)機(jī)理(核心機(jī)理)基本一致,因此構(gòu)建多組分機(jī)理即發(fā)展通用性好的核心機(jī)理和代表性組分的子機(jī)理.這樣所有代表性組分子機(jī)理與核心機(jī)理整合,就得到多組分航空煤油化學(xué)反應(yīng)機(jī)理.

    由于小分子碳?xì)淙剂系幕瘜W(xué)反應(yīng)機(jī)理得到廣泛的研究,其中Wang等20發(fā)展的USC-Mech II是一個(gè)得到普遍應(yīng)用和廣泛認(rèn)可的小分子碳?xì)淙剂戏磻?yīng)機(jī)理.該機(jī)理包含了H2/CO/C1-C4詳細(xì)反應(yīng)機(jī)理,由111組分、784反應(yīng)組成.并經(jīng)過(guò)著火延遲時(shí)間、火焰?zhèn)鞑ニ俣纫约盎鹧娼Y(jié)構(gòu)等大量的實(shí)驗(yàn)驗(yàn)證.因此本文直接采用USC-Mech II作為核心機(jī)理.

    因此,本文的主要工作是在核心機(jī)理的基礎(chǔ)上,進(jìn)一步構(gòu)建替代燃料模型中大分子碳?xì)浣M分的子機(jī)理.首先,針對(duì)鏈烷烴代表性組分正癸烷和正十二烷,從高、低溫兩個(gè)方面進(jìn)行子機(jī)理的構(gòu)建. You等21以USC-Mech II為核心機(jī)理構(gòu)建了包含正癸烷和正十二烷的反應(yīng)機(jī)理,并對(duì)高溫條件下的著火延遲時(shí)間、火焰?zhèn)鞑ニ俣?、組分分布進(jìn)行了詳細(xì)驗(yàn)證.鑒于該機(jī)理采用的核心機(jī)理與本文相同,而且能夠準(zhǔn)確描述高溫氧化過(guò)程.因此,本文通過(guò)反應(yīng)路徑和靈敏度分析(零維均質(zhì)模型,初始溫度1200 K),提取出其C12/C10-C4高溫氧化過(guò)程的重要反應(yīng)路徑及其主要基元反應(yīng),即構(gòu)成了本文正癸烷和正十二烷高溫反應(yīng)的子機(jī)理.低溫子機(jī)理考慮如圖5所示的集總反應(yīng)過(guò)程,對(duì)應(yīng)的反應(yīng)動(dòng)力學(xué)參數(shù)來(lái)源于Jia等22的正辛烷到正十六烷反應(yīng)機(jī)理.

    其次,對(duì)于環(huán)烷烴代表性組分乙基環(huán)己烷,采用上述相同的分析方法,從Sirjean等23發(fā)展的包含352組分、2083反應(yīng)的JetSurF 1.1詳細(xì)反應(yīng)機(jī)理中提取出乙基環(huán)己烷的重要氧化反應(yīng)路徑及其主要基元反應(yīng),即構(gòu)成乙基環(huán)己烷氧化過(guò)程的子機(jī)理.

    最后,對(duì)于芳香烴代表性組分對(duì)二甲苯,則以先前發(fā)展的多組分汽油替代燃料反應(yīng)機(jī)理24,25中甲苯機(jī)理為基礎(chǔ),添加對(duì)二甲苯氧化裂解至甲苯、苯、乙基苯的反應(yīng)機(jī)理.即構(gòu)成對(duì)二甲苯氧化機(jī)理.添加部分的機(jī)理主要來(lái)源于李玉陽(yáng)26發(fā)展的低壓預(yù)混火焰中對(duì)二甲苯氧化、裂解反應(yīng)機(jī)理,以及Gail和Dagaut27提出的攪拌反應(yīng)器中對(duì)二甲苯的氧化機(jī)理.

    至此,本文構(gòu)建了正癸烷、正十二烷、乙基環(huán)己烷和對(duì)二甲苯四種組分的反應(yīng)子機(jī)理(見(jiàn)Supporting Information),通過(guò)與核心機(jī)理進(jìn)行整合,并剔除冗余或重復(fù)的組分和反應(yīng),最終得到一個(gè)包含168組分和1089反應(yīng)的多組分RP-3航空煤油替代燃料的半詳細(xì)反應(yīng)動(dòng)力學(xué)模型,其合理性的驗(yàn)證見(jiàn)下一節(jié).

    圖5 低溫氧化過(guò)程集總反應(yīng)機(jī)理Fig.5 Lumped kinetic scheme of the low-temperature oxidation reactions ket:ketohydroperoxide

    圖6 計(jì)算的RP-3航空煤油和替代燃料火焰?zhèn)鞑ニ俣扰c實(shí)驗(yàn)結(jié)果比較Fig.6 Comparison of experimental and calculated laminar flame speeds of RP-3 aviation kerosene and surrogate fuels

    圖7 不同初始溫度下計(jì)算的RP-3航空煤油著火延遲時(shí)間與實(shí)驗(yàn)結(jié)果比較Fig.7 Comparison of experimental and calculated ignition delay times of RP-3 aviation kerosene at different initial temperatures

    4.2 RP-3替代燃料詳細(xì)反應(yīng)機(jī)理驗(yàn)證

    一個(gè)反應(yīng)機(jī)理有效性的重要標(biāo)志之一就是準(zhǔn)確預(yù)測(cè)燃料著火延遲時(shí)間和火焰?zhèn)鞑ニ俣?本文依據(jù)現(xiàn)有的實(shí)驗(yàn)結(jié)果詳細(xì)比較了RP-3航空煤油的著火延遲時(shí)間、火焰?zhèn)鞑ニ俣?計(jì)算是在化學(xué)反應(yīng)動(dòng)力學(xué)程序CHEMKIN28上完成的.對(duì)于著火延遲時(shí)間和火焰?zhèn)鞑ニ俣鹊挠?jì)算,分別采用零維均質(zhì)模型和一維預(yù)混火焰模型.

    圖6對(duì)比了初始?jí)毫?.1 MPa、溫度為403 K、化學(xué)計(jì)量比范圍為0.7-1.4時(shí)計(jì)算的RP-3航空煤油替代燃料的火焰?zhèn)鞑ニ俣扰c本文的實(shí)驗(yàn)測(cè)量值.由圖6可以看出,本文半詳細(xì)機(jī)理能夠準(zhǔn)確預(yù)測(cè)RP-3替代燃料的火焰?zhèn)鞑ニ俣?并能很好地描述實(shí)際RP-3航空煤油的燃燒速率.

    曾文等14利用激波管實(shí)驗(yàn)臺(tái)架詳細(xì)測(cè)量了不同初始?jí)毫Α⒒瘜W(xué)計(jì)量比下,RP-3著火延遲時(shí)間隨溫度的變化曲線.圖7給出了本文機(jī)理預(yù)測(cè)的RP-3著火延遲時(shí)間,并與曾文等的實(shí)驗(yàn)結(jié)果進(jìn)行了比較.由圖7可見(jiàn),化學(xué)計(jì)量比為0.5時(shí),本文機(jī)理預(yù)測(cè)的RP-3著火延遲時(shí)間與實(shí)驗(yàn)值一致.而化學(xué)計(jì)量比為1.0時(shí),本文機(jī)理的計(jì)算結(jié)果略小于實(shí)驗(yàn)測(cè)量值.但總體而言,本文機(jī)理能夠較為準(zhǔn)確地預(yù)測(cè)實(shí)際RP-3航空煤油的著火延遲時(shí)間,反映了實(shí)際RP-3航空煤油的著火特性.與此同時(shí),也間接反映了四組分替代燃料模型能夠預(yù)測(cè)實(shí)際RP-3航空煤油的著火特性.

    5 結(jié)論

    以國(guó)產(chǎn)航空煤油為對(duì)象,提出了摩爾分?jǐn)?shù)分別為40%正癸烷、42%正十二烷、13%乙基環(huán)己烷和5%對(duì)二甲苯的四組分RP-3航空煤油替代燃料模型,并通過(guò)實(shí)驗(yàn)充分驗(yàn)證了替代燃料模型與實(shí)際RP-3航空煤油在理化特性上的相似性.

    采用對(duì)沖火焰實(shí)驗(yàn)臺(tái)架,測(cè)量了不同化學(xué)計(jì)量比下RP-3航空煤油、以及本文提出的四組分替代燃料的層流火焰?zhèn)鞑ニ俣?對(duì)比結(jié)果表明本文提出的替代燃料能夠準(zhǔn)確描述實(shí)際RP-3航空煤油的燃燒速率.

    在四組分替代燃料模型的基礎(chǔ)上,進(jìn)一步發(fā)展了包含168組分、1089反應(yīng)的半詳細(xì)反應(yīng)動(dòng)力學(xué)模型,并用著火延時(shí)間和火焰?zhèn)魉俣葘?duì)機(jī)理進(jìn)行了充分驗(yàn)證.驗(yàn)證結(jié)果表明本文的反應(yīng)機(jī)理能夠較準(zhǔn)確地預(yù)測(cè)RP-3航空煤油的著火延遲時(shí)間和火焰?zhèn)鞑ニ俣?較好地預(yù)測(cè)RP-3航空煤油著火和燃燒特性.

    Supporting Information:Sub-mechanisms have been included.This information is available free of charge via the internet at http://www.whxb.pku.edu.cn.

    References

    (1) Ma,H.A.;Jie,M.Z.;Zeng,W.;Chen,X.X.Journal of Aerospace Power 2013,28(3),1139.[馬洪安,解茂昭,曾 文,陳瀟瀟.航空動(dòng)力學(xué)報(bào),2013,28(3),1139.]

    (2) Fan,X.J.;Yu,G.Journal of Propulsion Technology 2006,27 (2),187.[范學(xué)軍,俞 剛.推進(jìn)技術(shù),2006,27(2),187.]

    (3) Violi,A.;Yan,S.;Eddings,E.G.;Sarofim,A.F.;Granata,S.; Faravelli,T.;Ranzi,E.Combustion Science and Technology 2002,174(11-12),399.doi:10.1080/00102200215080

    (4) Dooley,S.;Won,S.H.;Heyne,J.;Farouk,T.I.;Ju,Y.;Dryer,F. L.;Kumar,K.;Hui,X.;Sung,C.J.;Wang,H.;Oehlschlaeger, M.A.;Iyer,V.;Iyer,S.;Litzinger,T.A.;Santoro,R.J.; Malewicki,T.;Brezinsky,K.Combustion and Flame 2012,159 (4),1444.doi:10.1016/j.combustflame.2011.11.002

    (5) Zheng,D.;Zhong,B.J.Acta Phys.-Chim.Sin.2012,28(9), 2029. [鄭 東,鐘北京.物理化學(xué)學(xué)報(bào).2012,28(9),2029.] doi:10.3866/PKU.WHXB201207042

    (6) Pitz,W.J.;Cernansky,N.P.;Dryer,F.L.;Egolfopoulos,F.N.; Farrell,J.T.;Friend,D.G.;Pitsch,H.SAE Technical Paper 2007-01-0175,2007.doi:10.4271/2007-01-0175

    (7) Pitz,W.J.;Mueller,C.J.Progress in Energy and Combustion Science 2011,37(3),330.doi:10.1016/j.pecs.2010.06.004

    (8) Dagaut,P.;Reuillon,M.;Boettner,J.C.;Cathonnet,M. Symposium(International)on Combustion 1994,25(1),919. doi:10.1016/S0082-0784(06)80727-7

    (9) Dooley,S.;Won,S.H.;Chaos,M.;Heyne,J.;Ju,Y.;Dryer,F. L.;Kumar,K.;Sung,C.J.;Wang,H.;Oehlschlaeger,M.A.; Santoro,R.J.;Litzinger,T.A.Combustion and Flame 2010,157 (12),2333.doi:10.1016/j.combustflame.2010.07.001

    (10) Guéret,C.;Cathonnet,M.;Boettner,J.C.;Gaillard,F. Symposium(International)on Combustion 1991,23(1),211. doi:10.1016/S0082-0784(06)80261-4

    (11) Montgomery,C.;Cannon,S.;Mawid,M.;Sekar,B.Reduced Chemical Kinetic Mechanisms for JP-8 Combustion.In 40th AIAAAerospace Sciences Meeting&Exhibit,Aerospace Sciences Meetings,Reno,Nevada,USA.Jan 14-17,2002; AIAAMember,Ed.;American Institute ofAeronautics and Astronautics:Reston,VA,2002.doi:10.2514/6.2002-336

    (12) Malewicki,T.;Gudiyella,S.;Brezinsky,K.Combustion and Flame 2013,160(1),17.doi:10.1016/j.combustflame. 2012.09.013

    (13) Xiao,B.G.;Yang,S.H.;Zhao,H.Y.;Qian,W.Q.;Le,J.L. Journal of Aerospace Power 2010,25(9),1948.[肖保國(guó),楊順華,趙慧勇,錢(qián)煒祺,樂(lè)嘉陵.航空動(dòng)力學(xué)報(bào),2010,25(9), 1948.]

    (14) Zeng,W.;Li,H.X.;Ma,H.A.;Liang,S.;Cheng,B.D.Journal of Propulsion Technology 2014,35(8),1139. [曾 文,李海霞,馬洪安,梁 雙,陳保東.推進(jìn)技術(shù),2014,35(8),1139.]

    (15) Dagaut,P.Physical Chemistry Chemical Physics 2002,4(11), 2079.doi:10.1039/b110787a

    (16) Edwards,T.;Maurice,L.Q.Journal of Propulsion and Power 2001,17(2),461.doi:10.2514/2.5765

    (17) Zhu,Y.H.;Yu,C.X.;Li,Z.M.;Mi,Z.T.;Zhang,X.W. Petrochemical Technology 2006,35(12),1152.[朱玉紅,余彩香,李子木,米鎮(zhèn)濤,張香文.石油化工,2006,35(12),1152.]

    (18) Holley,A.T.;You,X.Q.;Dames,E.;Wang,H.;Egolfopoulos, F.N.Proceedings of the Combustion Institute 2009,32(1), 1157.doi:10.1016/j.proci.2008.05.067

    (19) Yu,W.M.;Zhong,B.J.;Yuan,Z.;Wang,G.Z.Journal of Propulsion Technology 2014,35(1),70.[于維銘,鐘北京,袁振,王治國(guó).推進(jìn)技術(shù),2014,35(1),70.]

    (20) Wang,H.;You,X.;Joshi,A.V.;Davis,S.G.;Laskin,A.; Egolfopoulos,F.;Chung K.L.High-Temperature Combustion Reaction Model of H2/CO/C1-C4Compounds.http://ignis.usc. edu/USC_Mech_II.htm(accessed Oct 1,2014).

    (21) You,X.;Egolfopoulos,F.N.;Wang,H.Proceedings of the Combustion Institute 2009,32(1),403.doi:10.1016/j. proci.2008.06.041

    (22) Chang,Y.;Jia,M.;Liu,Y.;Li,Y.;Xie,M.;Yin,H.Energy& Fuels 2013,27(6),3467.doi:10.1021/ef400460d

    (23) Sirjean,B.;Dames,E.;Sheen,D.A.;Egolfopoulos,F.N.; Wang,H.;Davidson,D.F.;Hanson,R.K.;Pitsch,H.;Bowman, C.T.;Law,C.K.;Tsang,W.;Cernansky,N.P.;Miller,D.L.; Violi,A.;Lindstedt,R.P.AHigh-Temperature Chemical Kinetic Model of n-Alkane,Cyclohexane,and Methyl-,Ethyl-, n-Propyl and n-Butyl-cyclohexane Oxidation at High Temperatures,JetSurF version 1.1.http://melchior.usc.edu/ JetSurF/JetSurF1.1(accessed Oct 1,2014).

    (24) Zhong,B.J.;Zheng,D.Fuel 2014,128(15),458.

    (25) Zheng,D.;Zhang,Y.P.;Zhong,B.J.Acta Phys.-Chim.Sin. 2013,29(6),1154.[鄭 東,張?jiān)迄i,鐘北京.物理化學(xué)學(xué)報(bào), 2013,29(6),1154.]doi:10.3866/PKU.WHXB201303201

    (26) Li,Y.Y.Experimental and Kinetic Modeling Study of Premixed Aromatic Hydrocarbon Flames at Low Pressure.Ph.D. Dissertation,University of Science and Technology of China, Hefei,2010.[李玉陽(yáng).芳烴燃料低壓預(yù)混火焰的實(shí)驗(yàn)和動(dòng)力學(xué)模型研究[D].合肥:中國(guó)科學(xué)技術(shù)大學(xué),2010.]

    (27) Ga?l,S.;Dagaut,P.Combustion and Flame 2005,141(3),281. doi:10.1016/j.combustflame.2004.12.020

    (28) Kee,R.J.;Rupley,F.M.;Miller,J.A.CHEMKIN Release 4.1; Reaction Design:San Diego,CA.2006.

    RP-3 Aviation Kerosene Surrogate Fuel and the Chemical Reaction Kinetic Model

    ZHENG Dong1YU Wei-Ming1,2ZHONG Bei-Jing1,*
    (1School of Aerospace,Tsinghua University,Beijing 100084,P.R.China;2BAIC Group Industrial Investment Co.,Ltd.,Beijing 100021,P.R.China)

    A four-component RP-3 aviation kerosene surrogate fuel,comprising 40%n-decane/42%ndodecane/13%ethycyclohexane/5%p-xylene(molar fraction),was presented.Experiments showed the physical and chemical similarity of the surrogate fuel to the real RP-3.Counterflow,twin-flame experiments were used to determine the laminar flame speeds of both the real and the surrogate fuel and showed that the surrogate fuel accurately modeled the burning rate of real RP-3.Asemi-detailed chemical reaction mechanism for ignition and oxidation of the RP-3 surrogate fuel that consists of 168 species and 1089 elementary reactions has been developed.Experimental results validate the model and highlight its ability to accurately predict the ignition delay times and laminar flame speeds of real RP-3.

    RP-3 aviation kerosene;Surrogate fuel;Reaction mechanism;Laminar flame speed; Counter-flow flame

    O643

    10.3866/PKU.WHXB201501231www.whxb.pku.edu.cn

    Rceived:October 13,2014;Revised:January 22,2015;Published on Web:January 23,2015.

    ?Corresponding author.Email:zhongbj@tsinghua.edu.cn;Tel:+86-10-62772928.

    The project was supported by the National Natural Science Foundation of China(91441113).

    國(guó)家自然科學(xué)基金(91441113)資助項(xiàng)目

    猜你喜歡
    傳播速度煤油機(jī)理
    代謝綜合征患者臂踝脈搏波傳播速度與頸動(dòng)脈粥樣硬化的關(guān)系
    隔熱纖維材料的隔熱機(jī)理及其應(yīng)用
    煤層氣吸附-解吸機(jī)理再認(rèn)識(shí)
    新雷
    一類廣義canmassa—Holm方程的無(wú)限傳播速度與漸近行為
    為什么煤油燈需要借助燈芯才能燃燒
    霧霾機(jī)理之問(wèn)
    正丁醇/煤油混合物非預(yù)混燃燒壓力振蕩特性
    高能合成煤油GN-1理化性能及應(yīng)用分析
    載人航天(2016年5期)2016-10-19 09:24:10
    DNTF-CMDB推進(jìn)劑的燃燒機(jī)理
    变态另类丝袜制服| 亚洲成人免费电影在线观看| 国产精品乱码一区二三区的特点| 尾随美女入室| 99久国产av精品| 男女视频在线观看网站免费| 69人妻影院| 69人妻影院| 俄罗斯特黄特色一大片| 国产免费av片在线观看野外av| 深夜精品福利| 色在线成人网| 欧美日本亚洲视频在线播放| 久久精品国产清高在天天线| 国产一区二区激情短视频| 一区福利在线观看| 色吧在线观看| 久久中文看片网| 在线免费观看的www视频| 美女高潮喷水抽搐中文字幕| 3wmmmm亚洲av在线观看| 精品午夜福利在线看| 91在线精品国自产拍蜜月| 亚洲va在线va天堂va国产| 美女高潮的动态| 别揉我奶头 嗯啊视频| 色尼玛亚洲综合影院| 麻豆av噜噜一区二区三区| 男人狂女人下面高潮的视频| 一卡2卡三卡四卡精品乱码亚洲| 很黄的视频免费| 麻豆国产av国片精品| 3wmmmm亚洲av在线观看| 亚洲电影在线观看av| 亚洲精品在线观看二区| eeuss影院久久| 日本精品一区二区三区蜜桃| 露出奶头的视频| 日本 欧美在线| 精品无人区乱码1区二区| 中文字幕免费在线视频6| 国产精品一区二区免费欧美| 在线国产一区二区在线| 国产白丝娇喘喷水9色精品| 日本黄大片高清| 亚洲美女黄片视频| 日本 av在线| 久久精品影院6| 国产精品女同一区二区软件 | 亚洲美女搞黄在线观看 | 在线观看美女被高潮喷水网站| 我要搜黄色片| or卡值多少钱| а√天堂www在线а√下载| 亚洲最大成人手机在线| 干丝袜人妻中文字幕| 亚洲天堂国产精品一区在线| 国产精品野战在线观看| 婷婷亚洲欧美| 身体一侧抽搐| 亚洲成人中文字幕在线播放| 精品乱码久久久久久99久播| 国产精品自产拍在线观看55亚洲| 午夜福利在线观看吧| 韩国av在线不卡| 中文字幕人妻熟人妻熟丝袜美| 免费观看精品视频网站| 欧美精品啪啪一区二区三区| 国产乱人视频| 日本 av在线| 99在线人妻在线中文字幕| 国产精品免费一区二区三区在线| 亚洲国产欧洲综合997久久,| 少妇裸体淫交视频免费看高清| 久久九九热精品免费| 男人狂女人下面高潮的视频| 国产在线男女| 熟女电影av网| 午夜福利欧美成人| 波多野结衣高清作品| 亚洲一级一片aⅴ在线观看| 国产黄a三级三级三级人| 校园人妻丝袜中文字幕| 春色校园在线视频观看| 亚洲欧美精品综合久久99| 黄色丝袜av网址大全| 欧美日本视频| ponron亚洲| 午夜福利欧美成人| 亚洲内射少妇av| 日本a在线网址| 日本欧美国产在线视频| 一本久久中文字幕| av女优亚洲男人天堂| 久久午夜亚洲精品久久| 日韩高清综合在线| 中文字幕精品亚洲无线码一区| 精品午夜福利视频在线观看一区| 亚洲男人的天堂狠狠| 欧美潮喷喷水| 欧美精品啪啪一区二区三区| 国产精品一区二区性色av| 成人国产麻豆网| 亚洲狠狠婷婷综合久久图片| 国产精品久久久久久av不卡| 日日干狠狠操夜夜爽| 国产黄a三级三级三级人| 午夜影院日韩av| 国产精品福利在线免费观看| 成人国产麻豆网| 日本精品一区二区三区蜜桃| 日本一二三区视频观看| 国产伦人伦偷精品视频| 狠狠狠狠99中文字幕| 久久久色成人| 夜夜夜夜夜久久久久| 又爽又黄无遮挡网站| 99久国产av精品| 久久精品国产亚洲av天美| 久久精品国产鲁丝片午夜精品 | 91久久精品电影网| 免费看av在线观看网站| av视频在线观看入口| 久久久精品欧美日韩精品| 俄罗斯特黄特色一大片| 亚洲熟妇熟女久久| 亚洲中文日韩欧美视频| 久久国产精品人妻蜜桃| 淫妇啪啪啪对白视频| 久久精品国产亚洲av涩爱 | 亚洲五月天丁香| 亚洲欧美日韩高清专用| 日韩大尺度精品在线看网址| 两个人的视频大全免费| 免费不卡的大黄色大毛片视频在线观看 | 国产私拍福利视频在线观看| 亚洲精品国产成人久久av| а√天堂www在线а√下载| 一卡2卡三卡四卡精品乱码亚洲| 国产探花极品一区二区| 日韩欧美国产一区二区入口| 欧美性猛交黑人性爽| x7x7x7水蜜桃| netflix在线观看网站| 日韩欧美国产一区二区入口| 欧美日韩中文字幕国产精品一区二区三区| 日韩 亚洲 欧美在线| 大又大粗又爽又黄少妇毛片口| 亚洲中文日韩欧美视频| 内射极品少妇av片p| 欧美最新免费一区二区三区| 亚洲av二区三区四区| 国产精品一区二区三区四区免费观看 | 成人精品一区二区免费| 麻豆久久精品国产亚洲av| 天堂√8在线中文| 久久九九热精品免费| 色尼玛亚洲综合影院| 国产精品1区2区在线观看.| 久久精品国产自在天天线| 99久久无色码亚洲精品果冻| 国产69精品久久久久777片| 色av中文字幕| 成人高潮视频无遮挡免费网站| 亚洲自拍偷在线| 美女xxoo啪啪120秒动态图| 国产真实乱freesex| 成人性生交大片免费视频hd| 国产大屁股一区二区在线视频| 白带黄色成豆腐渣| 亚洲精品日韩av片在线观看| 男女做爰动态图高潮gif福利片| av在线天堂中文字幕| 日本一二三区视频观看| 国产视频内射| 亚洲图色成人| 日日摸夜夜添夜夜添av毛片 | 我要看日韩黄色一级片| 中国美女看黄片| 夜夜夜夜夜久久久久| 九九爱精品视频在线观看| 日韩大尺度精品在线看网址| 色综合婷婷激情| 免费无遮挡裸体视频| 欧美三级亚洲精品| 最新中文字幕久久久久| 九九久久精品国产亚洲av麻豆| 亚洲18禁久久av| 日韩av在线大香蕉| 亚洲av免费高清在线观看| av女优亚洲男人天堂| 国产在线精品亚洲第一网站| 三级男女做爰猛烈吃奶摸视频| 久久99热这里只有精品18| 欧美激情在线99| 成人鲁丝片一二三区免费| 成人一区二区视频在线观看| 免费在线观看影片大全网站| 亚洲国产欧美人成| 一本精品99久久精品77| 婷婷色综合大香蕉| 一进一出抽搐动态| 在线观看免费视频日本深夜| 狠狠狠狠99中文字幕| 国内精品久久久久精免费| 亚洲自拍偷在线| 美女cb高潮喷水在线观看| 亚洲午夜理论影院| 十八禁网站免费在线| 久久草成人影院| 三级国产精品欧美在线观看| 国产高清有码在线观看视频| 久久久久久久久中文| 色综合亚洲欧美另类图片| 三级毛片av免费| 色综合色国产| 精品久久国产蜜桃| 国产 一区 欧美 日韩| 免费看a级黄色片| 美女cb高潮喷水在线观看| 国产精品亚洲美女久久久| 美女高潮喷水抽搐中文字幕| 中出人妻视频一区二区| 午夜免费成人在线视频| 国产人妻一区二区三区在| 国产精品一区二区三区四区免费观看 | 色吧在线观看| 精品99又大又爽又粗少妇毛片 | 亚洲精品456在线播放app | 成人二区视频| 日韩精品青青久久久久久| 嫩草影院入口| 成人永久免费在线观看视频| 国产一区二区在线av高清观看| 欧美极品一区二区三区四区| av中文乱码字幕在线| 亚洲欧美清纯卡通| 色5月婷婷丁香| 精品乱码久久久久久99久播| 日韩欧美在线二视频| 国产精品久久久久久亚洲av鲁大| 日本欧美国产在线视频| 在线观看66精品国产| 亚洲最大成人av| 午夜久久久久精精品| 最好的美女福利视频网| 国产成人av教育| 在线a可以看的网站| 18禁在线播放成人免费| 十八禁国产超污无遮挡网站| avwww免费| 欧美色欧美亚洲另类二区| 欧美激情久久久久久爽电影| 久久精品国产亚洲av香蕉五月| 最近最新免费中文字幕在线| 亚洲在线观看片| 欧美一区二区精品小视频在线| 床上黄色一级片| 在线免费观看不下载黄p国产 | 男女之事视频高清在线观看| 欧美成人免费av一区二区三区| 国语自产精品视频在线第100页| 久久久久久久久大av| 国产精品嫩草影院av在线观看 | 91久久精品电影网| 亚洲狠狠婷婷综合久久图片| 午夜福利欧美成人| 在线观看午夜福利视频| 啦啦啦观看免费观看视频高清| 国内少妇人妻偷人精品xxx网站| 精品久久国产蜜桃| 日本 欧美在线| 级片在线观看| 欧美日韩综合久久久久久 | 日韩欧美一区二区三区在线观看| 最近最新免费中文字幕在线| 国产v大片淫在线免费观看| 我要看日韩黄色一级片| 麻豆国产97在线/欧美| 成年人黄色毛片网站| 免费av毛片视频| 综合色av麻豆| 非洲黑人性xxxx精品又粗又长| a级一级毛片免费在线观看| 欧美最新免费一区二区三区| 亚洲一级一片aⅴ在线观看| 少妇人妻精品综合一区二区 | 中国美白少妇内射xxxbb| 精品久久久久久,| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美日韩无卡精品| 精品欧美国产一区二区三| av福利片在线观看| 少妇猛男粗大的猛烈进出视频 | 久久久久久久久中文| 国产精品人妻久久久影院| 久久精品国产亚洲av天美| 亚洲欧美日韩无卡精品| 欧美3d第一页| 2021天堂中文幕一二区在线观| 成人av在线播放网站| 国产亚洲精品久久久久久毛片| 两个人视频免费观看高清| 国产亚洲精品av在线| 长腿黑丝高跟| 精品久久久噜噜| 亚洲av.av天堂| 欧美另类亚洲清纯唯美| 少妇被粗大猛烈的视频| 国产亚洲欧美98| 黄色配什么色好看| 欧美性猛交黑人性爽| 啦啦啦观看免费观看视频高清| 亚洲国产精品sss在线观看| 春色校园在线视频观看| 最近最新中文字幕大全电影3| 桃红色精品国产亚洲av| a级毛片a级免费在线| 国产淫片久久久久久久久| 大又大粗又爽又黄少妇毛片口| 三级男女做爰猛烈吃奶摸视频| 亚洲av日韩精品久久久久久密| 亚洲久久久久久中文字幕| 国产亚洲精品久久久久久毛片| 国产一区二区激情短视频| 国产精品久久久久久久久免| 天堂网av新在线| 国产爱豆传媒在线观看| 亚洲性夜色夜夜综合| 国产视频内射| 一级黄色大片毛片| 九色成人免费人妻av| 国产单亲对白刺激| 变态另类成人亚洲欧美熟女| 女生性感内裤真人,穿戴方法视频| 午夜激情福利司机影院| 国产爱豆传媒在线观看| bbb黄色大片| av在线蜜桃| 一本久久中文字幕| 日本-黄色视频高清免费观看| 午夜福利在线观看吧| 成年版毛片免费区| 久久午夜福利片| 成人国产综合亚洲| 嫩草影院入口| 女同久久另类99精品国产91| 俺也久久电影网| 欧美日韩综合久久久久久 | 久久久久精品国产欧美久久久| 男人舔女人下体高潮全视频| 国产精品久久电影中文字幕| 成人国产麻豆网| 人妻久久中文字幕网| 日韩欧美精品v在线| 国内精品久久久久久久电影| 日韩人妻高清精品专区| x7x7x7水蜜桃| 村上凉子中文字幕在线| 干丝袜人妻中文字幕| 国产三级在线视频| 干丝袜人妻中文字幕| 午夜福利18| 麻豆av噜噜一区二区三区| 亚洲真实伦在线观看| 午夜激情欧美在线| 2021天堂中文幕一二区在线观| 在线看三级毛片| 人人妻,人人澡人人爽秒播| 亚洲成人久久爱视频| ponron亚洲| 可以在线观看的亚洲视频| 久久久久久九九精品二区国产| 在现免费观看毛片| 国产精品综合久久久久久久免费| 亚洲成人久久爱视频| 国产乱人视频| 一级黄片播放器| 91午夜精品亚洲一区二区三区 | eeuss影院久久| 亚洲av电影不卡..在线观看| 午夜福利欧美成人| 亚洲第一电影网av| 免费电影在线观看免费观看| 成人综合一区亚洲| 国产精品日韩av在线免费观看| 亚洲av成人精品一区久久| 亚洲国产欧洲综合997久久,| 天堂av国产一区二区熟女人妻| 国产精品日韩av在线免费观看| 小蜜桃在线观看免费完整版高清| 乱人视频在线观看| 伊人久久精品亚洲午夜| 亚洲aⅴ乱码一区二区在线播放| 99久久九九国产精品国产免费| 欧美zozozo另类| 色综合色国产| a在线观看视频网站| 熟女电影av网| 成人av在线播放网站| 国产乱人视频| 啦啦啦啦在线视频资源| 久久精品久久久久久噜噜老黄 | 国产一区二区亚洲精品在线观看| 久久久久久久久中文| 最近视频中文字幕2019在线8| 日日摸夜夜添夜夜添av毛片 | 欧美潮喷喷水| 色视频www国产| 亚洲国产精品合色在线| 变态另类丝袜制服| 男女那种视频在线观看| 黄色女人牲交| 国产v大片淫在线免费观看| 国产精品1区2区在线观看.| 99热只有精品国产| 三级毛片av免费| 99视频精品全部免费 在线| 国产一区二区三区av在线 | 色综合站精品国产| 特大巨黑吊av在线直播| 不卡一级毛片| 久久6这里有精品| av专区在线播放| 国产91精品成人一区二区三区| 欧美中文日本在线观看视频| 久久热精品热| 国产亚洲精品久久久com| 亚洲av五月六月丁香网| 听说在线观看完整版免费高清| 亚洲美女视频黄频| 亚洲图色成人| 国产男人的电影天堂91| 久久久久久久久中文| 99九九线精品视频在线观看视频| 欧美+日韩+精品| 亚洲国产精品成人综合色| 精品不卡国产一区二区三区| 国产精品嫩草影院av在线观看 | 中亚洲国语对白在线视频| 精品不卡国产一区二区三区| 午夜爱爱视频在线播放| 亚洲自拍偷在线| 亚洲avbb在线观看| 国产成年人精品一区二区| 99久久精品一区二区三区| 极品教师在线免费播放| 精品福利观看| 国产亚洲精品久久久com| 精品人妻偷拍中文字幕| 国产单亲对白刺激| 久久精品国产清高在天天线| 国产aⅴ精品一区二区三区波| 亚洲熟妇熟女久久| 久久久久久久久久成人| 中国美女看黄片| 全区人妻精品视频| 午夜激情福利司机影院| 成人国产麻豆网| 干丝袜人妻中文字幕| 精品国产三级普通话版| 成人综合一区亚洲| 美女xxoo啪啪120秒动态图| 国产精品一区二区三区四区免费观看 | 亚洲综合色惰| 亚洲欧美日韩高清在线视频| 噜噜噜噜噜久久久久久91| 内射极品少妇av片p| 国产精品久久视频播放| 亚洲专区国产一区二区| 嫩草影院入口| 在线天堂最新版资源| 亚洲成人久久爱视频| 精品一区二区三区人妻视频| 禁无遮挡网站| 三级国产精品欧美在线观看| 国产亚洲欧美98| 成人国产麻豆网| 午夜精品一区二区三区免费看| 国产一区二区三区在线臀色熟女| 国产白丝娇喘喷水9色精品| 桃红色精品国产亚洲av| 人人妻,人人澡人人爽秒播| 一卡2卡三卡四卡精品乱码亚洲| 亚洲一级一片aⅴ在线观看| 国产在线男女| 亚洲国产精品成人综合色| 成年女人看的毛片在线观看| 在线免费观看不下载黄p国产 | 成人美女网站在线观看视频| 俄罗斯特黄特色一大片| 变态另类成人亚洲欧美熟女| 国产成人aa在线观看| 亚洲成人中文字幕在线播放| 久久久久久久亚洲中文字幕| 噜噜噜噜噜久久久久久91| 老熟妇乱子伦视频在线观看| 久9热在线精品视频| 亚洲乱码一区二区免费版| 久久草成人影院| 欧美日韩精品成人综合77777| 婷婷丁香在线五月| av国产免费在线观看| 欧美极品一区二区三区四区| 午夜福利欧美成人| 国产在视频线在精品| 免费不卡的大黄色大毛片视频在线观看 | 男女之事视频高清在线观看| 亚洲性久久影院| 国产精品精品国产色婷婷| 国产欧美日韩精品一区二区| 国产一区二区激情短视频| 免费人成在线观看视频色| АⅤ资源中文在线天堂| 成人美女网站在线观看视频| 日本精品一区二区三区蜜桃| 淫妇啪啪啪对白视频| 国产精品一区二区三区四区久久| av在线老鸭窝| 国产精品一区二区三区四区久久| 国内精品久久久久精免费| 久久精品影院6| 欧美最黄视频在线播放免费| 男女边吃奶边做爰视频| 亚洲,欧美,日韩| 久久天躁狠狠躁夜夜2o2o| 亚洲第一区二区三区不卡| 一区二区三区四区激情视频 | 在线观看66精品国产| 免费看美女性在线毛片视频| 欧美日本视频| 国产伦人伦偷精品视频| 国产精品98久久久久久宅男小说| 毛片女人毛片| 亚洲精品亚洲一区二区| 最新在线观看一区二区三区| 3wmmmm亚洲av在线观看| 尤物成人国产欧美一区二区三区| 亚洲av成人av| 搡老熟女国产l中国老女人| 91在线精品国自产拍蜜月| 国产精品一区二区三区四区免费观看 | 久久精品国产亚洲av涩爱 | 亚洲精品国产成人久久av| 亚洲欧美日韩高清在线视频| 日韩欧美三级三区| 国产又黄又爽又无遮挡在线| 97超视频在线观看视频| 久久久久久久亚洲中文字幕| 十八禁网站免费在线| 国产精品野战在线观看| 成年女人永久免费观看视频| 国产色婷婷99| 一夜夜www| 1000部很黄的大片| 中文字幕免费在线视频6| 国国产精品蜜臀av免费| 国产乱人伦免费视频| 国产成人av教育| 色综合亚洲欧美另类图片| 国产亚洲精品av在线| 桃红色精品国产亚洲av| 亚洲国产精品sss在线观看| 黄色欧美视频在线观看| 91在线精品国自产拍蜜月| 免费在线观看影片大全网站| 欧美另类亚洲清纯唯美| 欧美性猛交黑人性爽| 国产老妇女一区| 啦啦啦啦在线视频资源| 久久人人爽人人爽人人片va| 伦理电影大哥的女人| 亚洲成人久久爱视频| 精品一区二区免费观看| 亚洲国产日韩欧美精品在线观看| 男人的好看免费观看在线视频| 国产 一区精品| 大又大粗又爽又黄少妇毛片口| 日韩av在线大香蕉| 又粗又爽又猛毛片免费看| 久久久久精品国产欧美久久久| 俄罗斯特黄特色一大片| 身体一侧抽搐| 欧美高清性xxxxhd video| 中文字幕熟女人妻在线| 国产熟女欧美一区二区| 国产精品久久久久久久久免| 亚洲国产欧美人成| 欧美区成人在线视频| 欧美中文日本在线观看视频| 亚洲内射少妇av| 国产乱人伦免费视频| 麻豆精品久久久久久蜜桃| 91麻豆精品激情在线观看国产| 老司机深夜福利视频在线观看| 国产淫片久久久久久久久| 男女做爰动态图高潮gif福利片| 在线观看舔阴道视频| 久久久色成人| 久久这里只有精品中国| 欧美日韩国产亚洲二区| 久9热在线精品视频| 久久婷婷人人爽人人干人人爱| 久久久久免费精品人妻一区二区| 亚洲欧美日韩高清专用| 久久婷婷人人爽人人干人人爱| 欧美日韩国产亚洲二区| 成人av在线播放网站| 香蕉av资源在线| 成人高潮视频无遮挡免费网站| 少妇人妻精品综合一区二区 | 精品一区二区三区视频在线观看免费| 成熟少妇高潮喷水视频| 中文字幕熟女人妻在线| 色哟哟哟哟哟哟|