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

    增壓直噴汽油機(jī)中潤(rùn)滑油液滴誘發(fā)早燃機(jī)制的數(shù)值研究

    2016-12-23 01:50:18李耀庭黃勇成孟凡勝張文嘉陳銳
    關(guān)鍵詞:過(guò)氧庚烷汽油機(jī)

    李耀庭,黃勇成,孟凡勝,張文嘉,陳銳

    (1.西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安;2.一汽解放汽車有限公司無(wú)錫柴油機(jī)廠,214026,江蘇無(wú)錫)

    ?

    增壓直噴汽油機(jī)中潤(rùn)滑油液滴誘發(fā)早燃機(jī)制的數(shù)值研究

    李耀庭1,黃勇成1,孟凡勝1,張文嘉1,陳銳2

    (1.西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安;2.一汽解放汽車有限公司無(wú)錫柴油機(jī)廠,214026,江蘇無(wú)錫)

    通過(guò)利用多維計(jì)算流體動(dòng)力學(xué)(CFD)耦合一個(gè)多組分汽油替代物的骨架反應(yīng)機(jī)理,對(duì)增壓直噴汽油機(jī)中潤(rùn)滑油液滴誘發(fā)早燃的機(jī)制進(jìn)行了研究。首先,在模擬增壓高強(qiáng)化汽油機(jī)壓縮上止點(diǎn)熱力學(xué)狀態(tài)的定容圓柱網(wǎng)格內(nèi),研究了潤(rùn)滑油液滴的存在對(duì)汽油/空氣混合氣自燃過(guò)程的影響,結(jié)果表明:以正庚烷和過(guò)氧氫酮分別作為潤(rùn)滑油蒸發(fā)產(chǎn)物時(shí)潤(rùn)滑油液滴的存在均可縮短混合氣的著火延遲時(shí)間,尤其是以過(guò)氧氫酮作為潤(rùn)滑油蒸發(fā)產(chǎn)物時(shí)混合氣的著火延遲時(shí)間縮短更為明顯;隨著潤(rùn)滑油液滴粒徑的增大,混合氣的著火延遲時(shí)間先縮短后延長(zhǎng),潤(rùn)滑油液滴溫度和混合氣的溫度、壓力升高會(huì)進(jìn)一步縮短混合氣的著火延遲時(shí)間。然后,以過(guò)氧氫酮作為潤(rùn)滑油蒸發(fā)產(chǎn)物,在增壓直噴汽油機(jī)動(dòng)網(wǎng)格內(nèi)模擬了潤(rùn)滑油液滴存在時(shí)缸內(nèi)混合氣的自燃過(guò)程,結(jié)果表明:潤(rùn)滑油液滴蒸發(fā)釋放出了著火性能較好的組分,縮短了液滴周圍混合氣的著火延遲時(shí)間,導(dǎo)致液滴周圍混合氣在火花點(diǎn)火之前自燃,從而引發(fā)早燃。最后,根據(jù)以上結(jié)論提出了一個(gè)潤(rùn)滑油液滴誘發(fā)早燃的機(jī)制。

    增壓直噴汽油機(jī);早燃;潤(rùn)滑油液滴;骨架反應(yīng)機(jī)理;數(shù)值研究

    發(fā)動(dòng)機(jī)小型化是車用汽油機(jī)提高燃油經(jīng)濟(jì)性和降低CO2排放的重要措施之一,一般結(jié)合使用缸內(nèi)直噴技術(shù)和廢氣渦輪增壓技術(shù)可以實(shí)現(xiàn)發(fā)動(dòng)機(jī)小型化。隨著增壓比和功率密度的不斷提高,增壓直噴汽油機(jī)在低速、高負(fù)荷工況下出現(xiàn)了一種新的爆震模式——早燃爆震,也稱作超級(jí)爆震,其誘因是缸內(nèi)混合氣在火花點(diǎn)火前已發(fā)生預(yù)燃,從而產(chǎn)生極高的爆壓和壓力波動(dòng)。發(fā)生超級(jí)爆震時(shí),最高缸壓可達(dá)到30 MPa以上,壓力波動(dòng)幅度可達(dá)到20 MPa以上[1],對(duì)發(fā)動(dòng)機(jī)的破壞性極大,甚至造成火花塞電極熔斷、活塞頂部燒蝕、連桿彎曲等嚴(yán)重后果,導(dǎo)致發(fā)動(dòng)機(jī)損毀[2]。超級(jí)爆震具有偶發(fā)性和間歇性特征,無(wú)法通過(guò)加濃混合氣、加強(qiáng)壁面?zhèn)鳠峄虿捎酶咝镣橹灯蛠?lái)避免,所以成為目前發(fā)動(dòng)機(jī)增壓小型化技術(shù)發(fā)展的最大阻礙。

    早燃是火花點(diǎn)火發(fā)動(dòng)機(jī)在常規(guī)點(diǎn)火之前出現(xiàn)的異常著火燃燒現(xiàn)象。王志等人的研究發(fā)現(xiàn),超級(jí)爆震前一定發(fā)生了早燃,為了杜絕超級(jí)爆震的發(fā)生,必須抑制早燃[2]。國(guó)內(nèi)外眾多機(jī)構(gòu)對(duì)早燃進(jìn)行了研究,提出了多種可能的誘發(fā)機(jī)制。通過(guò)總結(jié)歸納這些機(jī)制,可以將早燃點(diǎn)的來(lái)源分為結(jié)構(gòu)熱點(diǎn)(火花塞、氣門、活塞等熾熱表面)和非結(jié)構(gòu)熱點(diǎn)(潤(rùn)滑油液滴、燃燒室顆粒物等)兩大類。到目前為止,雖然對(duì)于早燃的誘發(fā)原因?qū)W術(shù)界尚無(wú)定論,但是近年來(lái)大量的研究結(jié)果表明,結(jié)構(gòu)熱點(diǎn)對(duì)早燃的產(chǎn)生基本沒(méi)有影響,而非結(jié)構(gòu)熱點(diǎn)中的潤(rùn)滑油液滴則被認(rèn)為是最有可能引發(fā)增壓直噴汽油機(jī)缸內(nèi)混合氣早燃的重要物質(zhì)。

    Dahnz等人通過(guò)分析實(shí)驗(yàn)數(shù)據(jù)和一維模擬結(jié)果,最早指出了潤(rùn)滑油液滴是誘發(fā)早燃的可能原因[1]。Kalghatgi等人提出早燃可能是混合氣中混入的潤(rùn)滑油長(zhǎng)鏈成分自燃引起的[3]。Palaveev等人發(fā)現(xiàn),增加噴霧碰壁能夠顯著提升早燃發(fā)生的幾率,說(shuō)明誘發(fā)早燃的物質(zhì)很可能來(lái)自潤(rùn)滑油的蒸發(fā)產(chǎn)物或者由潤(rùn)滑油和燃油相互反應(yīng)生成的[4]。以上研究均表明,誘發(fā)混合氣早燃的物質(zhì)很有可能來(lái)自于潤(rùn)滑油,但是誘發(fā)機(jī)制和具體物質(zhì)的確定還有待進(jìn)一步研究。

    基于以上背景,本文充分發(fā)揮了數(shù)值模擬技術(shù)的獨(dú)特優(yōu)勢(shì),通過(guò)多維計(jì)算流體動(dòng)力學(xué)耦合化學(xué)反應(yīng)動(dòng)力學(xué),從正面對(duì)潤(rùn)滑油液滴誘發(fā)早燃的機(jī)制進(jìn)行了詳細(xì)的數(shù)值研究,開(kāi)創(chuàng)性地探討了不同潤(rùn)滑油蒸發(fā)產(chǎn)物替代物誘發(fā)早燃的可能性,研究了潤(rùn)滑油液滴粒徑、溫度和混合氣溫度、壓力對(duì)汽油/空氣混合氣著火延遲時(shí)間及早燃的影響,最后對(duì)潤(rùn)滑油液滴存在時(shí)增壓直噴汽油機(jī)缸內(nèi)混合氣的自燃過(guò)程進(jìn)行了模擬驗(yàn)證。

    1 計(jì)算模型

    1.1 化學(xué)反應(yīng)動(dòng)力學(xué)骨架機(jī)理模型

    研究早燃主要涉及到缸內(nèi)工質(zhì)的自燃和火焰的傳播過(guò)程,所以需要采用化學(xué)反應(yīng)動(dòng)力學(xué)模型準(zhǔn)確模擬燃料/空氣混合氣的著火延遲時(shí)間和層流火焰?zhèn)鞑ニ俾省?/p>

    本文選用異辛烷和正庚烷的混合物作為汽油替代物,通過(guò)分析和借鑒前人建立的PRF(primary reference fuel)骨架反應(yīng)機(jī)理,利用路徑分析法建立了一個(gè)新的能夠滿足上述要求的化學(xué)反應(yīng)動(dòng)力學(xué)骨架反應(yīng)機(jī)理,構(gòu)建思路如下。

    首先,參考Tsurushima機(jī)理初步建立了正庚烷和異辛烷的低溫反應(yīng)機(jī)理[5]。然后,根據(jù)C3H5是異辛烷低溫反應(yīng)階段的主要烴類生成物[6],將Tsurushima模型中的反應(yīng)

    C6H13CO+O2?C3H6+C3H6+CO+HO2

    修改為

    C6H13CO+O2?C3H7+C3H5+CO+HO2

    再根據(jù)Tanaka的研究[7],添加了正庚烷和異辛烷的耦合反應(yīng)

    C8H18+C7H15?C7H15+C8H17

    高溫反應(yīng)區(qū)借鑒Liu等[8]的思路建立了烷基的高溫裂解反應(yīng),同時(shí)添加了C1-C3以及H2的一些小分子反應(yīng)。

    最終,骨架反應(yīng)機(jī)理包含了41種組分和125個(gè)反應(yīng)[9]。圖1給出了該機(jī)理在混合氣壓力p為4 MPa時(shí)以化學(xué)計(jì)量比混合的不同辛烷值汽油/空氣混合氣著火延遲時(shí)間的計(jì)算值與實(shí)驗(yàn)值的對(duì)比,其中激波管的實(shí)驗(yàn)結(jié)果來(lái)自Fieweger等的結(jié)果[10]。圖2給出了該機(jī)理對(duì)不同辛烷值汽油/空氣混合氣在常規(guī)工況下層流火焰?zhèn)鞑ニ俾实挠?jì)算值與實(shí)驗(yàn)值的對(duì)比,其中實(shí)驗(yàn)數(shù)據(jù)來(lái)自Huang等的結(jié)果[11]。圖中所有計(jì)算值均來(lái)自氣相化學(xué)反應(yīng)動(dòng)力學(xué)程序CHEMKIN PRO的計(jì)算結(jié)果。

    由圖1和圖2可以看出,本文的骨架反應(yīng)機(jī)理能夠較為準(zhǔn)確地模擬不同辛烷值汽油/空氣混合氣的著火延遲時(shí)間和層流火焰?zhèn)鞑ニ俾?可用于預(yù)測(cè)內(nèi)燃機(jī)缸內(nèi)工質(zhì)的自燃及燃燒過(guò)程。所以,下文中建立的多維數(shù)值模型將通過(guò)耦合本文骨架反應(yīng)機(jī)理進(jìn)行早燃誘發(fā)機(jī)制的研究。

    圖1 著火延遲時(shí)間計(jì)算值與實(shí)驗(yàn)值的對(duì)比

    圖2 層流火焰?zhèn)鞑ニ俾视?jì)算值與實(shí)驗(yàn)值的對(duì)比

    1.2 多維數(shù)值模型

    增壓直噴汽油機(jī)普遍存在噴霧碰壁現(xiàn)象,氣缸壁上的潤(rùn)滑油膜被噴霧油束稀釋導(dǎo)致潤(rùn)滑油膜黏度、表面張力下降,進(jìn)而液滴飛濺與蒸發(fā)的可能性增加。潤(rùn)滑油液滴蒸發(fā)后燃油蒸氣與空氣形成可燃混合氣,由于潤(rùn)滑油組分的著火性較好,潤(rùn)滑油液滴周圍混合氣的著火延遲時(shí)間縮短,所以混合氣便有可能在點(diǎn)火之前自燃而引發(fā)早燃[12]。圖3給出了燃燒室內(nèi)潤(rùn)滑油液滴誘發(fā)混合氣早燃的過(guò)程。

    圖3 潤(rùn)滑油液滴誘發(fā)早燃的過(guò)程

    據(jù)此,本文使用通用流體分析軟件包STAR-CD建立了定容圓柱網(wǎng)格,用以研究潤(rùn)滑油液滴的存在及其粒徑、溫度和混合氣溫度、壓力對(duì)汽油/空氣混合氣著火延遲時(shí)間以及早燃的影響,同時(shí)建立了增壓直噴汽油機(jī)的動(dòng)網(wǎng)格,用以驗(yàn)證在汽油機(jī)缸內(nèi)定容圓柱網(wǎng)格中得到的結(jié)論,并對(duì)缸內(nèi)潤(rùn)滑油液滴誘發(fā)早燃的詳細(xì)過(guò)程進(jìn)行研究。

    1.2.1 定容圓柱網(wǎng)格 模型1采用直徑為84 mm(約等于氣缸直徑)、高度為10 mm的定容圓柱網(wǎng)格,該網(wǎng)格用來(lái)模擬增壓直噴汽油機(jī)壓縮上止點(diǎn)的熱力學(xué)狀態(tài),網(wǎng)格總數(shù)為72 000。計(jì)算區(qū)域內(nèi)設(shè)定充滿當(dāng)量比為1的均質(zhì)汽油/空氣混合氣,汽油選用正庚烷和異辛烷的混合物替代,以耦合本文的PRF骨架反應(yīng)機(jī)理。辛烷值為93,是常規(guī)車用汽油辛烷值。通過(guò)用戶子程序droico.f定義了潤(rùn)滑油液滴的位置、粒徑以及溫度等參數(shù),潤(rùn)滑油液滴的位置固定在圓截面右側(cè)5 mm高度處,如圖4所示。

    1.2.2 發(fā)動(dòng)機(jī)動(dòng)網(wǎng)格 模型2為基于某增壓直噴汽油機(jī)(發(fā)動(dòng)機(jī)參數(shù)見(jiàn)表1)建立的發(fā)動(dòng)機(jī)的動(dòng)網(wǎng)格,上止點(diǎn)網(wǎng)格總數(shù)為353 756,下止點(diǎn)網(wǎng)格總數(shù)為602 860,如圖5所示。根據(jù)增壓直噴汽油機(jī)易于產(chǎn)生早燃的低速大負(fù)荷工況,設(shè)定發(fā)動(dòng)機(jī)轉(zhuǎn)速為1 600 r/min,最大扭矩為280 N·m,增壓壓力為0.19 MPa。相應(yīng)地設(shè)置了計(jì)算的初始條件和邊界條件,其中進(jìn)、排氣口溫度分別設(shè)定為320 K和990 K,進(jìn)、排氣口壓力分別設(shè)定為0.19 MPa和0.22 MPa。計(jì)算的湍流模型選擇k-ε模型,算法采用PISO算法。通過(guò)建立多組分噴霧蒸發(fā)模型來(lái)預(yù)測(cè)直噴發(fā)動(dòng)機(jī)噴霧油束的發(fā)展及霧化蒸發(fā)過(guò)程;通過(guò)耦合本文汽油替代物骨架反應(yīng)機(jī)理來(lái)模擬汽油的著火延遲時(shí)間和層流火焰速率;利用特征時(shí)間法建立基于復(fù)雜化學(xué)反應(yīng)動(dòng)力學(xué)的湍流燃燒模型,來(lái)模擬多組分汽油混合氣的湍流燃燒過(guò)程[9]。

    表1 發(fā)動(dòng)機(jī)基本參數(shù)

    為了研究增壓直噴汽油機(jī)工作過(guò)程中潤(rùn)滑油液滴蒸發(fā)對(duì)缸內(nèi)工質(zhì)著火及火焰?zhèn)鞑ミ^(guò)程的影響,假設(shè)在發(fā)動(dòng)機(jī)壓縮行程后半段,由于活塞的減速運(yùn)動(dòng)使得潤(rùn)滑油液滴因慣性力從活塞間隙處的油膜內(nèi)飛濺出來(lái)。潤(rùn)滑油液滴的初始運(yùn)動(dòng)速度定義為活塞運(yùn)動(dòng)的瞬時(shí)最大速度,為7.8 m/s,初始運(yùn)動(dòng)方向定義為與垂直方向成10°的方向,如圖6所示。

    圖5 增壓直噴汽

    圖6 潤(rùn)滑油液滴的初始位油機(jī)動(dòng)網(wǎng)格置及運(yùn)動(dòng)方向示意圖

    2 結(jié)果分析

    本文通過(guò)耦合化學(xué)反應(yīng)動(dòng)力學(xué),依次在定容圓柱網(wǎng)格和增壓直噴汽油機(jī)動(dòng)網(wǎng)格中,針對(duì)潤(rùn)滑油液滴對(duì)汽油/空氣混合氣自燃過(guò)程的影響進(jìn)行了多維數(shù)值模擬研究。值得注意的是,在STAR-CD中,潤(rùn)滑油替代物的選取將直接影響模擬結(jié)果,潤(rùn)滑油液滴的替代物和蒸發(fā)產(chǎn)物需要分別進(jìn)行定義。在本文的骨架反應(yīng)機(jī)理的41種反應(yīng)物中,正庚烷能夠較為恰當(dāng)?shù)啬M潤(rùn)滑油組分良好的著火性能[1],所以選用正庚烷作為潤(rùn)滑油液滴替代物,并在定容圓柱網(wǎng)格內(nèi)研究了潤(rùn)滑油液滴對(duì)汽油/空氣混合氣著火延遲時(shí)間的影響。著火延遲時(shí)間定義為混合氣溫度上升梯度為400 K所經(jīng)歷的時(shí)間。

    圖7給出了潤(rùn)滑油液滴粒徑對(duì)汽油/空氣混合氣著火延遲時(shí)間的影響,其中潤(rùn)滑油蒸發(fā)產(chǎn)物為正庚烷,潤(rùn)滑油液滴溫度為440 K,混合氣溫度為760 K,混合氣壓力為4 MPa(與發(fā)動(dòng)機(jī)增壓壓力為0.2 MPa時(shí)壓縮上止點(diǎn)壓力相當(dāng))。由圖中可以看出:當(dāng)混合氣中存在用正庚烷替代的潤(rùn)滑油液滴時(shí),混合氣的著火延遲時(shí)間比沒(méi)有潤(rùn)滑油液滴時(shí)(3.0 ms)縮短;隨著潤(rùn)滑油液滴粒徑的增大,混合氣的著火延遲時(shí)間先逐漸減小后逐漸增大;當(dāng)液滴的粒徑在45 μm左右時(shí),混合氣的著火延遲時(shí)間(2.21 ms)最短。這是因?yàn)?當(dāng)液滴粒徑較小時(shí),其氧化過(guò)程不足以提供火焰?zhèn)鞑ニ枰臒嶂?火焰不能夠向外傳播;當(dāng)液滴粒徑較大時(shí),液滴蒸發(fā)吸熱現(xiàn)象較明顯,混合氣周圍溫度降低,氧化過(guò)程較緩慢。

    圖7 汽油/空氣混合氣著火延遲時(shí)間與潤(rùn)滑油液滴粒徑的關(guān)系

    值得注意的是,當(dāng)正庚烷直接作為潤(rùn)滑油蒸發(fā)產(chǎn)物時(shí),混合氣的著火延遲時(shí)間最短達(dá)到2.21 ms,比無(wú)潤(rùn)滑油時(shí)最多縮短了0.79 ms。在發(fā)動(dòng)機(jī)轉(zhuǎn)速為1 600 r/min的低速工況下,即使不考慮活塞的下行運(yùn)動(dòng),混合氣著火延遲角度也只能減小7.6°,這說(shuō)明正庚烷作為潤(rùn)滑油蒸發(fā)產(chǎn)物時(shí)汽油/空氣混合氣的著火延遲時(shí)間仍然較長(zhǎng),直接誘發(fā)混合氣自燃的可能性較低,混合氣的早燃可能是由除正庚烷以外的與潤(rùn)滑油相關(guān)的其他物質(zhì)導(dǎo)致的。

    Palaveev等利用著火性能更好的正十六烷作為潤(rùn)滑油蒸發(fā)產(chǎn)物進(jìn)行了模擬,發(fā)現(xiàn)直到壓縮行程末期缸內(nèi)混合氣仍未出現(xiàn)著火現(xiàn)象[4],說(shuō)明誘發(fā)早燃的根本物質(zhì)應(yīng)該不是正構(gòu)直鏈烷烴。烴類低溫燃燒的中間產(chǎn)物過(guò)氧化物的著火延遲時(shí)間比正構(gòu)直鏈烷烴更短,這可能是誘發(fā)缸內(nèi)混合氣早燃的根本物質(zhì)。根據(jù)本文骨架反應(yīng)機(jī)理的路徑分析,在低溫反應(yīng)區(qū),烷烴通過(guò)脫氫、加氧、異構(gòu)化、二次加氧、二次異構(gòu)化和最終的離解反應(yīng)生成了穩(wěn)定的中間基過(guò)氧氫酮,根據(jù)Magar等人的研究,過(guò)氧氫酮在整個(gè)低溫反應(yīng)階段是活性最強(qiáng)的穩(wěn)定中間基[13],所以本文認(rèn)為過(guò)氧氫酮更有可能是誘發(fā)早燃的根本物質(zhì)。

    圖8給出了定容圓柱網(wǎng)格內(nèi)不同辛烷值汽油/空氣混合氣自燃時(shí)過(guò)氧氫酮濃度隨時(shí)間的變化,其中潤(rùn)滑油液滴溫度為760 K,混合氣壓力為4 MPa。由圖中可以看出,隨著辛烷值的降低,即正庚烷比例的增加,混合氣內(nèi)過(guò)氧氫酮濃度峰值的出現(xiàn)時(shí)間提前,這說(shuō)明正庚烷的存在會(huì)促進(jìn)穩(wěn)定中間基過(guò)氧氫酮的生成。因此,在潤(rùn)滑油液滴(正庚烷)的作用下,液滴周圍混合氣會(huì)有較高濃度的過(guò)氧氫酮,當(dāng)其積累達(dá)到一定量之后就可能使混合氣在點(diǎn)火之前發(fā)生自燃。

    圖8 不同辛烷值混合氣自燃時(shí)過(guò)氧氫酮濃度的變化

    基于以上結(jié)論,為了簡(jiǎn)化問(wèn)題,本文選取過(guò)氧氫酮作為潤(rùn)滑油的蒸發(fā)產(chǎn)物,由此研究了潤(rùn)滑油液滴粒徑對(duì)混合氣著火延遲時(shí)間的影響,結(jié)果見(jiàn)圖7。對(duì)比圖中的結(jié)果可以看出:正庚烷、過(guò)氧氫酮隨液滴粒徑的變化趨勢(shì)相同,當(dāng)過(guò)氧氫酮作為潤(rùn)滑油蒸發(fā)產(chǎn)物時(shí),混合氣的著火延遲時(shí)間最短達(dá)到1.43 ms,比無(wú)潤(rùn)滑油液滴時(shí)縮短了1.57 ms,對(duì)應(yīng)的著火延遲角度減小了15°,說(shuō)明過(guò)氧氫酮的存在相較于正庚烷大大縮短了混合氣的著火延遲時(shí)間,過(guò)氧氫酮更有可能是導(dǎo)致缸內(nèi)工質(zhì)早燃的重要物質(zhì)。

    圖9給出了潤(rùn)滑油液滴溫度對(duì)汽油/空氣混合氣著火延遲時(shí)間的影響,其中潤(rùn)滑油蒸發(fā)產(chǎn)物分別為正庚烷和過(guò)氧氫酮,液滴粒徑為45 μm,混合氣溫度為760 K,混合氣壓力為4 MPa。由圖中可以看出,隨著潤(rùn)滑油液滴溫度的升高,混合氣的著火延遲時(shí)間逐漸縮短。這是因?yàn)橐旱蔚谋砻鎻埩宛ざ入S著溫度的升高而不斷降低,由此加快了液滴的破碎速度,縮短了蒸發(fā)所需時(shí)間;液滴溫度升高,蒸發(fā)所需的吸熱量有所降低,由此加快了液滴由液相向氣相的轉(zhuǎn)化,促進(jìn)了周圍混合氣的氧化反應(yīng)。

    圖9 汽油/空氣混合氣著火延遲時(shí)間與潤(rùn)滑油液滴溫度的關(guān)系

    圖10給出了潤(rùn)滑油液滴對(duì)不同溫度混合氣著火延遲時(shí)間的影響,其中潤(rùn)滑油蒸發(fā)產(chǎn)物分別為正庚烷和過(guò)氧氫酮,潤(rùn)滑油液滴溫度為460 K,液滴粒徑為45 μm,混合氣壓力為4 MPa。由圖中可以看出,隨著混合氣溫度的升高,混合氣的著火延遲時(shí)間逐漸縮短。這是因?yàn)榛旌蠚鉁囟鹊纳仙苯蛹铀倭藵?rùn)滑油液滴的蒸發(fā),促進(jìn)了潤(rùn)滑油蒸氣與汽油/空氣混合氣的混合,同時(shí)溫度的升高也直接加速了混合氣的氧化過(guò)程。

    圖10 汽油/空氣混合氣著火延遲時(shí)間與混合氣溫度的關(guān)系

    圖11 汽油/空氣混合氣著火延遲時(shí)間與混合氣壓力的關(guān)系

    圖11給出了潤(rùn)滑油液滴對(duì)不同壓力混合氣著火延遲時(shí)間的影響,其中潤(rùn)滑油蒸發(fā)產(chǎn)物分別為正庚烷和過(guò)氧氫酮,潤(rùn)滑油液滴溫度為460 K,液滴粒徑為45 μm,混合氣溫度為760 K。由圖中可以看出,隨著混合氣壓力的升高,混合氣的著火延遲時(shí)間逐漸縮短。這主要是當(dāng)量比一定時(shí),高壓環(huán)境下混合氣的濃度較高,而高濃度有利于低溫化學(xué)反應(yīng),加速了冷焰火焰的傳播,從而促進(jìn)了潤(rùn)滑油液滴周圍混合氣的低溫氧化過(guò)程和自燃。

    由圖9~圖11仍然可以看出:當(dāng)過(guò)氧氫酮作為潤(rùn)滑油蒸發(fā)產(chǎn)物時(shí),汽油/空氣混合氣的著火延遲時(shí)間相較于正庚烷大大縮短,說(shuō)明過(guò)氧氫酮更有可能是導(dǎo)致缸內(nèi)工質(zhì)早燃的物質(zhì)。

    本文在初次使用正庚烷作為缸內(nèi)飛濺的潤(rùn)滑油液滴的蒸發(fā)產(chǎn)物進(jìn)行缸內(nèi)工質(zhì)的燃燒計(jì)算時(shí)發(fā)現(xiàn),潤(rùn)滑油液滴的存在無(wú)法促進(jìn)缸內(nèi)工質(zhì)在火花點(diǎn)火之前自燃,改用過(guò)氧氫酮作為潤(rùn)滑油蒸發(fā)產(chǎn)物時(shí)則成功地預(yù)測(cè)了潤(rùn)滑油液滴周圍混合氣早燃火核的形成及發(fā)展過(guò)程。

    圖12給出了過(guò)氧氫酮作為潤(rùn)滑油蒸發(fā)產(chǎn)物時(shí)潤(rùn)滑油液滴周圍混合氣在火花點(diǎn)火之前發(fā)生自燃后的缸內(nèi)溫度場(chǎng)分布,計(jì)算時(shí)設(shè)定潤(rùn)滑油液滴的初始粒徑為45 μm,溫度為460 K。由圖中可以看出,火花點(diǎn)火角度在715°時(shí),潤(rùn)滑油附近的混合氣已發(fā)生自燃現(xiàn)象,有明確的火核形成,局部溫度迅速上升,火焰開(kāi)始向四周傳播,到火花點(diǎn)火角度720°時(shí),缸內(nèi)已燃面積已達(dá)到1/3左右?;鸹c(diǎn)火之后,由火花點(diǎn)火形成的火焰面與之前潤(rùn)滑油附近混合氣自燃形成的火焰面相交,使得火焰的傳播速度加快,并迅速充滿了整個(gè)燃燒室。

    圖13給出了715°自燃火核形成時(shí)過(guò)氧氫酮的濃度場(chǎng)分布。由圖中可以看出:活塞減速時(shí)飛濺出來(lái)的潤(rùn)滑油液滴因蒸發(fā)形成的過(guò)氧氫酮逐漸向四周擴(kuò)散,然后與汽油/空氣混合氣混合;火核中心處過(guò)氧氫酮具有較高的濃度,其加速了該區(qū)域的低溫氧化過(guò)程,促進(jìn)了火核的進(jìn)一步形成及火焰的擴(kuò)散。

    圖12 潤(rùn)滑油液滴周圍混合氣在火花點(diǎn)火之前、自燃后缸內(nèi)溫度場(chǎng)分布

    圖13 過(guò)氧氫酮濃度場(chǎng)分布

    圖14給出了715°燃燒流場(chǎng)的OH自由基濃度場(chǎng)分布。OH自由基濃度可以用來(lái)描述火焰面。由圖中可以看出,在自燃的開(kāi)始時(shí)刻,OH自由基在火核中心處有較高的濃度,其濃度場(chǎng)分布與過(guò)氧氫酮的濃度場(chǎng)分布類似。

    圖14 OH濃度場(chǎng)分布

    綜上所述,一個(gè)合理的假設(shè)機(jī)制可以用來(lái)解釋潤(rùn)滑油液滴誘發(fā)增壓直噴汽油機(jī)缸內(nèi)工質(zhì)的早燃:增壓直噴汽油機(jī)的噴霧碰壁現(xiàn)象導(dǎo)致了氣缸壁上的潤(rùn)滑油膜黏度、表面張力下降,使得部分潤(rùn)滑油液滴飛濺、蒸發(fā)后與汽油/空氣混合氣混合,潤(rùn)滑油組分(正庚烷)促進(jìn)了混合氣低溫氧化過(guò)程中穩(wěn)定中間基(過(guò)氧氫酮)的產(chǎn)生,使活塞間隙處殘留的穩(wěn)定中間基(過(guò)氧氫酮)具有較高的濃度,在經(jīng)過(guò)幾個(gè)循環(huán)的積累后進(jìn)入燃燒室,從而縮短了混合氣的著火延遲時(shí)間,導(dǎo)致混合氣發(fā)生早燃。由于整個(gè)過(guò)程需要滿足一定的物理和化學(xué)條件,而且該穩(wěn)定中間基(過(guò)氧氫酮)濃度需要多個(gè)循環(huán)的積累來(lái)提高,所以該假設(shè)也能夠間接解釋增壓直噴汽油機(jī)中超級(jí)爆震具有偶發(fā)性和間歇性的特點(diǎn)。

    3 結(jié) 論

    本文利用化學(xué)反應(yīng)動(dòng)力學(xué)和多維數(shù)值模擬的耦合對(duì)早燃的誘發(fā)機(jī)制進(jìn)行了研究,得出以下結(jié)論。

    (1)本文建立的多組分汽油替代物骨架反應(yīng)機(jī)理能夠較為準(zhǔn)確地模擬不同辛烷值汽油/空氣混合氣的著火延遲時(shí)間和層流火焰?zhèn)鞑ニ俾?可用于預(yù)測(cè)內(nèi)燃機(jī)缸內(nèi)工質(zhì)的自燃及燃燒過(guò)程。

    (2)當(dāng)發(fā)動(dòng)機(jī)缸內(nèi)存在易于蒸發(fā)的潤(rùn)滑油液滴時(shí),缸內(nèi)混合氣的著火延遲時(shí)間明顯縮短,所以潤(rùn)滑油液滴很有可能是增壓直噴汽油機(jī)中缸內(nèi)工質(zhì)發(fā)生早燃的一大誘因。

    (3)正庚烷作為潤(rùn)滑油蒸發(fā)產(chǎn)物時(shí),汽油/空氣混合氣的著火延遲時(shí)間仍較長(zhǎng),說(shuō)明正庚烷直接誘發(fā)早燃的可能性較低。過(guò)氧氫酮作為潤(rùn)滑油蒸發(fā)產(chǎn)物時(shí),混合氣的著火延遲時(shí)間明顯縮短,說(shuō)明過(guò)氧氫酮更有可能是引發(fā)增壓直噴汽油機(jī)缸內(nèi)工質(zhì)早燃的物質(zhì),具體誘發(fā)機(jī)制有待進(jìn)一步研究與驗(yàn)證。

    (4)隨著潤(rùn)滑油液滴粒徑的增大,汽油/空氣混合氣的著火延遲時(shí)間先縮短后延長(zhǎng),說(shuō)明存在一個(gè)最適宜的潤(rùn)滑油液滴粒徑,使得混合氣的著火延遲時(shí)間最短。潤(rùn)滑油液滴溫度和混合氣溫度、壓力的升高也會(huì)進(jìn)一步縮短汽油/空氣混合氣的著火延遲時(shí)間,誘發(fā)混合氣早燃。

    [1] DAHNZ C, HAN KM, SPICHER U. Investigations on pre-ignition in highly supercharged SI engines [J]. SAE International Journal of Engines, 2010, 3(1): 214-224.

    [2] 王志, 龍巖, 王建昕. 增壓汽油機(jī)中早燃和超級(jí)爆震的研究進(jìn)展 [J]. 汽車安全與節(jié)能學(xué)報(bào), 2015(1): 17-29. WANG Zhi, LONG Yan, WANG Jianxin. Research progress of pre-ignition and super-knock in boost gasoline engine [J]. Automotive Safety and Energy, 2015(1): 17-29.

    [3] KALGHATGI G T, BRADLEY D. Pre-ignition and ‘super-knock’ in turbo-charged spark-ignition engines [J]. International Journal of Engine Research, 2012, 13(4): 399-414.

    [4] PALAVEEV S, MAGAR M, KUBACH H, et al. Premature flame initiation in a turbocharged DISI engine: numerical and experimental investigations [J]. SAE International Journal of Engines, 2013, 6(1): 54-66.

    [5] TSURUSHIMA T. A new skeletal PRF kinetic model for HCCI combustion [J]. Proceedings of the Combustion Institute, 2009, 32(2): 2835-2841.

    [6] ZHENG J, YANG W, MILLER D L, et al. A skeletal chemical kinetic model for the HCCI combustion process, SAE 2002-01-0423 [R]. Washington, DC, USA: SAE, 2002.

    [7] TANAKA S, AYALA F, KECK J C. A reduced chemical kinetic model for HCCI combustion of primary reference fuels in a rapid compression machine [J]. Combustion and Flame, 2003, 133(4): 467-481.

    [8] LIU Y D, JIA M, XIE M Z, et al. Improvement on a skeletal chemical kinetic model of iso-octane for internal combustion engine by using a practical methodology [J]. Fuel, 2013, 103: 884-891.

    [9] 陳銳. 潤(rùn)滑油對(duì)增壓高強(qiáng)化汽油機(jī)早燃爆震影響的數(shù)值模擬研究 [D]. 西安: 西安交通大學(xué), 2014.

    [10]FIEWEGER K, BLUMENTHAL R, ADOMEIT G. Self-ignition of SI engine model fuels: a shock tube investigation at high pressure [J]. Combustion and Flame 1997, 109(4): 599-619.

    [11]HUANG Y, SUNG C J, ENG J A. Laminar flame speeds of primary reference fuels and reformer gas mixtures [J]. Combustion and Flame, 2004, 139(3): 239-251.

    [12]STEPINA V, VESELY V. Lubricants and special fluids [M]. Amsterdam, Netherlands: Elsevier, 1992: 12-35.

    [13]LAUER T, HEISS M, BOBICIC N, et al. A comprehensive simulation approach to irregular combustion, SAE 2014-01-1214 [R]. Washington, DC, USA: SAE, 2014.

    [本刊相關(guān)文獻(xiàn)鏈接]

    汪映,黃智勇,柯希春,等.液化石油氣對(duì)二甲醚-柴油雙燃料發(fā)動(dòng)機(jī)燃燒和排放特性的影響.2016,50(5):19-23.[doi:10.7652/xjtuxb201605003]

    鐘兵,洪偉,蘇巖,等.控制參數(shù)對(duì)增壓缸內(nèi)直噴汽油機(jī)部分負(fù)荷下微粒排放特性的影響.2016,50(5):95-100.[doi:10.7652/xjtuxb201605014]

    胡鵬,王偉,鄧康耀,等.排氣再循環(huán)與燃燒邊界條件耦合對(duì)增壓汽油機(jī)有效熱效率的影響.2016,50(1):34-40.[doi:10.7652/xjtuxb201601006]

    麻斌,董偉,高瑩,等.稀燃條件下廢氣再循環(huán)對(duì)缸內(nèi)直噴汽油機(jī)微粒排放粒徑分布的影響.2016,50(5):87-94.[doi:10.7652/xjtuxb201605013]

    高忠權(quán),李春艷,劉兵,等.采用離子電流法的發(fā)動(dòng)機(jī)非正常燃燒診斷.2015,49(5):1-6.[doi:10.7652/xjtuxb201505001]

    (編輯 苗凌)

    Numerical Investigations for Mechanism of Pre-Ignition Induced by Lubricant Oil Droplets in Supercharged DISI Engines

    LI Yaoting1,HUANG Yongcheng1,MENG Fansheng1,ZHANG Wenjia1,CHEN Rui2

    (1. School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China;2. Wuxi Diesel Engine Works, FAW Jiefang Automotive Co. Ltd., Wuxi, Jiangsu 214026, China)

    Multidimensional computational fluid dynamics coupled with a skeletal chemical kinetic model of gasoline substitute is used to investigate the mechanism of pre-ignition induced by lubricant oil droplets in a supercharged direct-injection spark-ignition (DISI) engine under high load and at low speed. The influence of lubricant oil droplet on the auto-ignition of gasoline and air mixture is evaluated in a constant volume combustion bomb with the thermodynamic conditions similar to those in the DISI engine cylinder at the compression top dead center. The results show that the ignition delay time of gasoline and air mixture is shortened by the presence of lubricant oil droplet when n-heptane and hydrogen peroxide ketone (C8KET) are used as lubricant oil evaporation substitutes respectively, and this behavior is more obvious when C8KET is used as lubricant oil evaporation substitute. The ignition delay time of gasoline and air mixture first decreases then increases with the increasing lubricant oil droplet diameter. In addition, the ignition delay time of fuel-air mixture is further shortened with the increasing lubricant oil droplet temperature and the temperature and pressure of the fuel-air mixture. A numerical investigation for auto-ignition of gasoline and air mixture with the presence of lubricant oil droplets in the supercharged DISI engine is carried out by using C8KET as lubricant oil evaporation substitute. The results show that the ignition delay time of the fuel-air mixture around the lubricant oil droplets is shortened by the constituents with good ignition performance evaporating from lubricating oil droplets, which causes auto-ignition of fuel-air mixture in this area before the spark ignition and leads to the occurrence of pre-ignition. A mechanism of pre-ignition induced by lubricant oil droplets is finally proposed following the conclusions above.

    supercharged DISI engine; pre-ignition; lubricant oil droplet; skeletal chemical kinetic model; numerical investigation

    2015-11-29。 作者簡(jiǎn)介:李耀庭(1992—),男,碩士生;黃勇成(通信作者),男,副教授。 基金項(xiàng)目:國(guó)家自然科學(xué)基金資助項(xiàng)目(50776069);陜西省自然科學(xué)基金資助項(xiàng)目(2015JM5236)。

    時(shí)間:2016-05-10

    10.7652/xjtuxb201607009

    TK411.2

    A

    0253-987X(2016)07-0051-07

    網(wǎng)絡(luò)出版地址:http:∥www.cnki.net/kcms/detail/61.1069.T.20160510.1519.006.html

    猜你喜歡
    過(guò)氧庚烷汽油機(jī)
    過(guò)氧自由基氧化對(duì)米糠蛋白體外胃蛋白酶消化性質(zhì)的影響
    引爆生活的一種鹽
    ——過(guò)氧碳酸鈉
    科學(xué)家解析青蒿素類過(guò)氧橋鍵合成機(jī)制,為合成青蒿素帶來(lái)希望
    微型圓管中正庚烷/空氣預(yù)混催化燃燒特性實(shí)驗(yàn)
    1,4-硫氮雜庚烷鹽酸鹽的簡(jiǎn)便合成
    聚丙烯成核劑雙環(huán)[2.2.1]-庚烷-2,3-二羧酸鈉的合成
    萃取精餾分離甲苯-正庚烷混合物的模擬研究
    2013年型Accord轎車用2.4L直列4缸汽油機(jī)的開(kāi)發(fā)
    Ford公司的1.0L Eco-Boost汽油機(jī)
    用于增壓直接噴射汽油機(jī)的噴束引導(dǎo)燃燒過(guò)程
    91午夜精品亚洲一区二区三区| 久久国产乱子免费精品| 91在线观看av| 免费看av在线观看网站| 舔av片在线| 国产伦精品一区二区三区视频9| av中文乱码字幕在线| 精品久久久久久久人妻蜜臀av| 人妻少妇偷人精品九色| 久久久久久久午夜电影| 欧美+亚洲+日韩+国产| 91在线观看av| 一卡2卡三卡四卡精品乱码亚洲| 欧美成人a在线观看| 国产一区亚洲一区在线观看| 国产成人精品久久久久久| 亚洲丝袜综合中文字幕| 久久99热6这里只有精品| 欧美zozozo另类| 日日撸夜夜添| 亚洲成人久久性| 亚洲成av人片在线播放无| 亚洲精品456在线播放app| 少妇猛男粗大的猛烈进出视频 | 啦啦啦啦在线视频资源| 亚洲七黄色美女视频| 长腿黑丝高跟| 久久精品91蜜桃| 国产在线男女| 成年免费大片在线观看| 搡老妇女老女人老熟妇| 久久99热这里只有精品18| 三级经典国产精品| 久久久a久久爽久久v久久| 在线免费观看不下载黄p国产| 最近中文字幕高清免费大全6| 国产蜜桃级精品一区二区三区| 不卡一级毛片| 国产亚洲精品av在线| 欧美极品一区二区三区四区| 久99久视频精品免费| 国产高清有码在线观看视频| 久久久成人免费电影| 一区二区三区高清视频在线| 伦精品一区二区三区| 观看免费一级毛片| 热99re8久久精品国产| 岛国在线免费视频观看| 99热只有精品国产| 在线观看66精品国产| 亚洲在线自拍视频| 亚洲无线观看免费| 白带黄色成豆腐渣| 免费观看精品视频网站| 欧美丝袜亚洲另类| 久久天躁狠狠躁夜夜2o2o| 床上黄色一级片| 亚洲精品成人久久久久久| 中文亚洲av片在线观看爽| 国产精品野战在线观看| 成年av动漫网址| 精品一区二区三区视频在线| 日本a在线网址| 国产一区亚洲一区在线观看| 色在线成人网| 日本黄色视频三级网站网址| 亚洲欧美成人精品一区二区| 亚洲精品乱码久久久v下载方式| 麻豆成人午夜福利视频| 女人十人毛片免费观看3o分钟| 国产成年人精品一区二区| 亚洲自拍偷在线| 色哟哟哟哟哟哟| 免费看光身美女| 精品国内亚洲2022精品成人| 精华霜和精华液先用哪个| 婷婷六月久久综合丁香| 免费看a级黄色片| 菩萨蛮人人尽说江南好唐韦庄 | 欧美最黄视频在线播放免费| 一个人看的www免费观看视频| 自拍偷自拍亚洲精品老妇| 99久久中文字幕三级久久日本| 嫩草影院精品99| 一本一本综合久久| 亚洲av第一区精品v没综合| 男人的好看免费观看在线视频| 在线免费观看不下载黄p国产| 国产aⅴ精品一区二区三区波| 老司机午夜福利在线观看视频| 久久久国产成人免费| 久久韩国三级中文字幕| 男插女下体视频免费在线播放| 欧美中文日本在线观看视频| 国产精品久久电影中文字幕| 亚洲高清免费不卡视频| 精品人妻熟女av久视频| 亚洲国产精品久久男人天堂| 亚洲成av人片在线播放无| 99视频精品全部免费 在线| 18禁黄网站禁片免费观看直播| 变态另类丝袜制服| 夜夜爽天天搞| 亚洲成a人片在线一区二区| 丝袜喷水一区| 特级一级黄色大片| 日本免费一区二区三区高清不卡| 婷婷六月久久综合丁香| 国产一区二区亚洲精品在线观看| 日本精品一区二区三区蜜桃| 男女之事视频高清在线观看| 秋霞在线观看毛片| 午夜福利成人在线免费观看| 免费人成在线观看视频色| 欧美人与善性xxx| 亚洲中文字幕日韩| 3wmmmm亚洲av在线观看| 国产精品国产三级国产av玫瑰| 久久精品影院6| 两个人视频免费观看高清| 18禁在线播放成人免费| 国产人妻一区二区三区在| 晚上一个人看的免费电影| 亚洲不卡免费看| 伊人久久精品亚洲午夜| av在线蜜桃| 少妇丰满av| 舔av片在线| 亚洲国产精品成人久久小说 | 成人鲁丝片一二三区免费| 国产色爽女视频免费观看| 精品少妇黑人巨大在线播放 | 欧美国产日韩亚洲一区| 亚洲人成网站高清观看| 真实男女啪啪啪动态图| 日日摸夜夜添夜夜添av毛片| 国产精品久久久久久久电影| 草草在线视频免费看| 一夜夜www| 别揉我奶头~嗯~啊~动态视频| 久久久久久久久久久丰满| 国产精品乱码一区二三区的特点| 波野结衣二区三区在线| 午夜福利18| 偷拍熟女少妇极品色| 欧美色视频一区免费| 国产单亲对白刺激| 国产伦精品一区二区三区四那| 永久网站在线| 色综合亚洲欧美另类图片| 欧美一区二区亚洲| 亚洲欧美日韩高清专用| 尤物成人国产欧美一区二区三区| 免费一级毛片在线播放高清视频| 熟女电影av网| 91av网一区二区| 在线观看美女被高潮喷水网站| 午夜福利视频1000在线观看| 精品久久久久久成人av| 亚洲国产日韩欧美精品在线观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲久久久久久中文字幕| 日本黄色片子视频| 成人鲁丝片一二三区免费| 搡老熟女国产l中国老女人| 内射极品少妇av片p| 国产伦一二天堂av在线观看| 熟妇人妻久久中文字幕3abv| 色综合色国产| 免费观看在线日韩| 精品一区二区免费观看| 亚洲精品粉嫩美女一区| 免费不卡的大黄色大毛片视频在线观看 | av在线蜜桃| 久久亚洲精品不卡| 免费看美女性在线毛片视频| 国产精品美女特级片免费视频播放器| 日韩成人av中文字幕在线观看 | 12—13女人毛片做爰片一| 久久精品影院6| 国产欧美日韩一区二区精品| 国内精品久久久久精免费| 身体一侧抽搐| 天堂√8在线中文| 国产片特级美女逼逼视频| 国产av麻豆久久久久久久| av在线观看视频网站免费| 欧美日韩乱码在线| 一级黄片播放器| 欧美xxxx性猛交bbbb| 亚洲一区高清亚洲精品| 亚洲精品乱码久久久v下载方式| 丝袜喷水一区| 亚洲中文日韩欧美视频| 内射极品少妇av片p| 变态另类成人亚洲欧美熟女| 国产女主播在线喷水免费视频网站 | 亚洲av电影不卡..在线观看| 极品教师在线视频| 国产精品亚洲美女久久久| 久久精品91蜜桃| 热99在线观看视频| 黄色视频,在线免费观看| 校园人妻丝袜中文字幕| а√天堂www在线а√下载| 成人午夜高清在线视频| 精品一区二区三区视频在线| 国内久久婷婷六月综合欲色啪| 婷婷亚洲欧美| 日韩精品有码人妻一区| 丰满乱子伦码专区| av视频在线观看入口| 色视频www国产| 青春草视频在线免费观看| 在线观看免费视频日本深夜| 国产成人91sexporn| 一本久久中文字幕| 在线播放国产精品三级| 精品福利观看| 国产不卡一卡二| 久久久精品欧美日韩精品| 国模一区二区三区四区视频| 少妇人妻精品综合一区二区 | 黑人高潮一二区| 久久久精品欧美日韩精品| 韩国av在线不卡| 亚洲天堂国产精品一区在线| 亚洲经典国产精华液单| 欧美一区二区精品小视频在线| 97人妻精品一区二区三区麻豆| 亚洲精品粉嫩美女一区| 人妻久久中文字幕网| 国产成人精品久久久久久| 国产精品av视频在线免费观看| 淫秽高清视频在线观看| 免费搜索国产男女视频| 午夜精品在线福利| 美女被艹到高潮喷水动态| 小说图片视频综合网站| 精品午夜福利在线看| 国产成人福利小说| av女优亚洲男人天堂| 欧美另类亚洲清纯唯美| 又黄又爽又免费观看的视频| 国产成人精品久久久久久| 大型黄色视频在线免费观看| 97人妻精品一区二区三区麻豆| 99热网站在线观看| 国产一区二区亚洲精品在线观看| 亚洲av第一区精品v没综合| 18禁裸乳无遮挡免费网站照片| 国产精品美女特级片免费视频播放器| 91精品国产九色| 老师上课跳d突然被开到最大视频| 午夜视频国产福利| 亚洲色图av天堂| 久久久久久九九精品二区国产| 91av网一区二区| 国产精品,欧美在线| 欧美性猛交黑人性爽| 少妇猛男粗大的猛烈进出视频 | 成人特级黄色片久久久久久久| 男女啪啪激烈高潮av片| 最新中文字幕久久久久| 日韩av在线大香蕉| 久久久国产成人免费| 精品午夜福利视频在线观看一区| 亚洲在线自拍视频| 中国美白少妇内射xxxbb| 免费观看在线日韩| 久久久久久久久大av| 国产成人福利小说| 69av精品久久久久久| 国产人妻一区二区三区在| 亚洲av不卡在线观看| 免费观看精品视频网站| 人妻制服诱惑在线中文字幕| 色哟哟哟哟哟哟| 五月伊人婷婷丁香| 有码 亚洲区| 美女内射精品一级片tv| 日韩 亚洲 欧美在线| 国产精品永久免费网站| 精品久久久久久久人妻蜜臀av| 网址你懂的国产日韩在线| 91午夜精品亚洲一区二区三区| 九九久久精品国产亚洲av麻豆| 久久精品国产清高在天天线| 男插女下体视频免费在线播放| 又爽又黄a免费视频| 18+在线观看网站| 国产老妇女一区| 国产精品无大码| 日韩欧美三级三区| 超碰av人人做人人爽久久| 亚洲最大成人手机在线| 特大巨黑吊av在线直播| 亚洲av中文字字幕乱码综合| 一级黄色大片毛片| 色在线成人网| 99在线人妻在线中文字幕| 99国产精品一区二区蜜桃av| 国产人妻一区二区三区在| 日本在线视频免费播放| 97超碰精品成人国产| 亚洲中文字幕日韩| 欧美日韩精品成人综合77777| 国产精品三级大全| 夜夜爽天天搞| 色视频www国产| 亚洲丝袜综合中文字幕| 婷婷精品国产亚洲av| 亚洲成人精品中文字幕电影| 国产黄a三级三级三级人| videossex国产| 老司机午夜福利在线观看视频| 免费观看精品视频网站| 国产高清不卡午夜福利| 成人三级黄色视频| 国产一区二区亚洲精品在线观看| 性色avwww在线观看| 老司机福利观看| 看黄色毛片网站| 日本-黄色视频高清免费观看| 最新中文字幕久久久久| 亚洲av一区综合| 成人性生交大片免费视频hd| 中文在线观看免费www的网站| 两性午夜刺激爽爽歪歪视频在线观看| 一区二区三区高清视频在线| av免费在线看不卡| 91久久精品国产一区二区成人| 色视频www国产| 中文字幕av在线有码专区| 精品熟女少妇av免费看| 一级a爱片免费观看的视频| 亚洲内射少妇av| 俺也久久电影网| 中文字幕免费在线视频6| 亚洲性久久影院| 日韩人妻高清精品专区| 亚洲综合色惰| 一区二区三区四区激情视频 | 乱系列少妇在线播放| 亚洲av成人精品一区久久| 小说图片视频综合网站| 成人国产麻豆网| 一区二区三区四区激情视频 | 99热6这里只有精品| 一级毛片aaaaaa免费看小| 欧美日韩一区二区视频在线观看视频在线 | 精品人妻视频免费看| 最近视频中文字幕2019在线8| 人妻少妇偷人精品九色| 淫妇啪啪啪对白视频| 熟女电影av网| 亚洲av.av天堂| 国模一区二区三区四区视频| 免费黄网站久久成人精品| 国内精品一区二区在线观看| 少妇人妻一区二区三区视频| 欧美性猛交╳xxx乱大交人| 秋霞在线观看毛片| 午夜视频国产福利| 22中文网久久字幕| 国产av不卡久久| 1024手机看黄色片| 亚洲欧美日韩高清在线视频| 亚洲中文字幕日韩| 99久久中文字幕三级久久日本| 久久午夜福利片| 尤物成人国产欧美一区二区三区| 三级经典国产精品| 国产高清激情床上av| 色播亚洲综合网| 免费黄网站久久成人精品| 少妇高潮的动态图| 一a级毛片在线观看| 69人妻影院| 久久草成人影院| 亚洲欧美日韩东京热| 精品人妻视频免费看| 一夜夜www| 人妻久久中文字幕网| 国产精品一区二区三区四区免费观看 | 变态另类丝袜制服| 别揉我奶头~嗯~啊~动态视频| 大又大粗又爽又黄少妇毛片口| 亚洲欧美精品自产自拍| 国产亚洲91精品色在线| 国产淫片久久久久久久久| 99精品在免费线老司机午夜| 亚洲第一区二区三区不卡| 国产成年人精品一区二区| 欧美色视频一区免费| 日韩一区二区视频免费看| 久久久久久九九精品二区国产| 男女边吃奶边做爰视频| 少妇人妻精品综合一区二区 | 乱系列少妇在线播放| 看黄色毛片网站| 午夜影院日韩av| 日日干狠狠操夜夜爽| 精品久久久久久成人av| 亚洲专区国产一区二区| 久久草成人影院| 又黄又爽又免费观看的视频| 国产不卡一卡二| 欧美xxxx黑人xx丫x性爽| 国产aⅴ精品一区二区三区波| 精品一区二区三区视频在线观看免费| 国产成年人精品一区二区| 成人av在线播放网站| 日韩在线高清观看一区二区三区| 日产精品乱码卡一卡2卡三| 欧美性猛交╳xxx乱大交人| 亚洲精品456在线播放app| 国语自产精品视频在线第100页| 九九久久精品国产亚洲av麻豆| 最近手机中文字幕大全| 国产精品电影一区二区三区| 国产精品久久久久久亚洲av鲁大| 我要看日韩黄色一级片| 如何舔出高潮| 成人亚洲欧美一区二区av| 99热只有精品国产| 亚洲成人久久爱视频| 色5月婷婷丁香| 日韩精品有码人妻一区| 俺也久久电影网| 免费看日本二区| 特级一级黄色大片| 日产精品乱码卡一卡2卡三| 99热只有精品国产| 国产麻豆成人av免费视频| 亚洲精华国产精华液的使用体验 | 日韩欧美一区二区三区在线观看| 夜夜夜夜夜久久久久| 老熟妇乱子伦视频在线观看| 女生性感内裤真人,穿戴方法视频| 直男gayav资源| 日产精品乱码卡一卡2卡三| 国产一区二区三区av在线 | 久久99热这里只有精品18| 免费观看的影片在线观看| 中文字幕av在线有码专区| 亚洲精品粉嫩美女一区| 久久久久久久亚洲中文字幕| 看非洲黑人一级黄片| 男人和女人高潮做爰伦理| 成人精品一区二区免费| 精品久久久久久久久亚洲| 久久午夜福利片| 国产精品野战在线观看| av国产免费在线观看| 亚洲精品国产av成人精品 | 亚洲自偷自拍三级| 少妇人妻精品综合一区二区 | 欧美丝袜亚洲另类| 少妇高潮的动态图| 国产精品一区二区三区四区久久| 三级毛片av免费| 成年版毛片免费区| 国产av不卡久久| 国产精品久久久久久久久免| 国产高清不卡午夜福利| 久久精品综合一区二区三区| 女的被弄到高潮叫床怎么办| 九九在线视频观看精品| 亚洲国产欧美人成| 亚洲久久久久久中文字幕| 国产精品综合久久久久久久免费| 亚洲av第一区精品v没综合| 丝袜喷水一区| 国产精品福利在线免费观看| 国模一区二区三区四区视频| 国产在线男女| 亚洲欧美成人综合另类久久久 | 你懂的网址亚洲精品在线观看 | 成人精品一区二区免费| 久久精品久久久久久噜噜老黄 | 国产精品99久久久久久久久| av在线观看视频网站免费| 欧美成人a在线观看| 欧美极品一区二区三区四区| 在线a可以看的网站| 中文字幕熟女人妻在线| 亚洲在线观看片| 久久精品夜夜夜夜夜久久蜜豆| a级毛色黄片| 国产高潮美女av| 亚洲va在线va天堂va国产| av.在线天堂| 欧美性感艳星| 在线播放无遮挡| 有码 亚洲区| 国产成人a区在线观看| 日韩欧美三级三区| 99热这里只有是精品50| 久久精品夜夜夜夜夜久久蜜豆| 最近的中文字幕免费完整| 精品一区二区三区视频在线观看免费| videossex国产| 一级黄色大片毛片| 超碰av人人做人人爽久久| 久久人妻av系列| 美女大奶头视频| 久99久视频精品免费| 亚洲经典国产精华液单| 最近视频中文字幕2019在线8| 久久久a久久爽久久v久久| 亚洲精品日韩av片在线观看| 国产淫片久久久久久久久| 国产白丝娇喘喷水9色精品| 亚洲熟妇中文字幕五十中出| 国产91av在线免费观看| 波多野结衣高清无吗| 男人的好看免费观看在线视频| 国产探花极品一区二区| 国产一区二区在线av高清观看| 别揉我奶头 嗯啊视频| 日本与韩国留学比较| 男人和女人高潮做爰伦理| 色综合色国产| 欧美一区二区亚洲| 99精品在免费线老司机午夜| 久久婷婷人人爽人人干人人爱| 美女 人体艺术 gogo| 国产极品精品免费视频能看的| 久久6这里有精品| 91在线观看av| 国产亚洲91精品色在线| 亚洲第一电影网av| 一级a爱片免费观看的视频| 热99re8久久精品国产| 精品熟女少妇av免费看| 色播亚洲综合网| 欧美+日韩+精品| a级毛片免费高清观看在线播放| 国产 一区 欧美 日韩| a级毛色黄片| 久久久久久九九精品二区国产| 国产成人影院久久av| 国产成人freesex在线 | 乱码一卡2卡4卡精品| 欧美日韩精品成人综合77777| 精品日产1卡2卡| 精华霜和精华液先用哪个| 成年av动漫网址| 亚洲精品色激情综合| 看黄色毛片网站| 久久综合国产亚洲精品| 欧美最黄视频在线播放免费| 国产精品一及| 精品99又大又爽又粗少妇毛片| 久久精品国产自在天天线| 成年免费大片在线观看| 欧美一区二区精品小视频在线| 乱系列少妇在线播放| 91久久精品电影网| 秋霞在线观看毛片| 3wmmmm亚洲av在线观看| 精品久久久久久久久久久久久| 亚洲精品粉嫩美女一区| 亚洲精华国产精华液的使用体验 | 亚洲av免费高清在线观看| 老熟妇仑乱视频hdxx| 秋霞在线观看毛片| 亚洲精品日韩av片在线观看| 在线播放国产精品三级| 国产成人福利小说| 国产精品亚洲美女久久久| 亚洲高清免费不卡视频| 国产精品一区二区性色av| aaaaa片日本免费| 久久精品夜夜夜夜夜久久蜜豆| 国产在线精品亚洲第一网站| 国产高清视频在线观看网站| 十八禁国产超污无遮挡网站| 我的女老师完整版在线观看| 亚洲国产精品国产精品| 婷婷六月久久综合丁香| 精品久久久久久成人av| 男女下面进入的视频免费午夜| 欧美成人精品欧美一级黄| 亚洲精品日韩av片在线观看| 亚洲一区高清亚洲精品| 午夜影院日韩av| 成年免费大片在线观看| 国产精品福利在线免费观看| 色视频www国产| 日韩欧美精品v在线| 色视频www国产| 精品人妻视频免费看| 免费看日本二区| 婷婷色综合大香蕉| 99久国产av精品| 97人妻精品一区二区三区麻豆| 性欧美人与动物交配| 亚洲美女黄片视频| 男人狂女人下面高潮的视频| 国产精品无大码| 如何舔出高潮| 热99在线观看视频| 国产乱人偷精品视频| 久久精品国产清高在天天线| 久久久a久久爽久久v久久| 丰满的人妻完整版| 国产乱人偷精品视频| 国产精品免费一区二区三区在线| 观看免费一级毛片| 波多野结衣高清作品| 亚洲丝袜综合中文字幕| 丰满乱子伦码专区| 亚洲成av人片在线播放无|