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

    燃燒室形狀對(duì)雙燃料發(fā)動(dòng)機(jī)性能影響的模擬分析

    2016-12-12 02:05:19王忠恕李慧楊舟李偉峰
    關(guān)鍵詞:圓柱形雙燃料啞鈴

    王忠恕, 李慧, 楊舟, 李偉峰

    (吉林大學(xué)汽車(chē)仿真與控制國(guó)家重點(diǎn)實(shí)驗(yàn)室, 吉林 長(zhǎng)春 130025)

    ?

    ·設(shè)計(jì)計(jì)算·

    燃燒室形狀對(duì)雙燃料發(fā)動(dòng)機(jī)性能影響的模擬分析

    王忠恕, 李慧, 楊舟, 李偉峰

    (吉林大學(xué)汽車(chē)仿真與控制國(guó)家重點(diǎn)實(shí)驗(yàn)室, 吉林 長(zhǎng)春 130025)

    利用STAR-CD軟件模擬研究了3種燃燒室形狀對(duì)柴油-天然氣雙燃料發(fā)動(dòng)機(jī)性能的影響,3種燃燒室分別為ω形燃燒室、八邊啞鈴形燃燒室和圓柱形燃燒室。研究發(fā)現(xiàn),八邊啞鈴形燃燒室因?yàn)闇p小了喉口直徑,增加了擠流強(qiáng)度,使得氣缸內(nèi)的湍動(dòng)能增強(qiáng),火焰?zhèn)鞑ニ俣燃涌?,燃料利用率提高,同時(shí),在燃燒室的底部設(shè)計(jì)凸臺(tái),能引導(dǎo)燃燒室內(nèi)的氣流運(yùn)動(dòng),并引導(dǎo)柴油向燃燒室的底部擴(kuò)散,促進(jìn)著火點(diǎn)的廣泛分布。因此,八邊啞鈴形燃燒室的缸內(nèi)平均壓力、平均溫度和指示熱效率最高,天然氣剩余比例最小。

    雙燃料發(fā)動(dòng)機(jī); 燃燒室; 燃燒; 排放; 仿真

    能源危機(jī)和環(huán)境污染不斷加劇,排放法規(guī)日趨嚴(yán)格,開(kāi)發(fā)利用新型清潔代用燃料迫在眉睫。天然氣因資源豐富、價(jià)格低廉、環(huán)保性好而備受關(guān)注,成為石油的主要替代品。天然氣的辛烷值較高,具有較好的抗爆性,與汽油機(jī)相比,天然氣發(fā)動(dòng)機(jī)可以采用更高的壓縮比以提高熱效率;與柴油機(jī)相比,天然氣發(fā)動(dòng)機(jī)中天然氣更容易與空氣混合,炭煙排放較少[1-2]。根據(jù)點(diǎn)火方式,天然氣發(fā)動(dòng)機(jī)可以分為火花點(diǎn)燃式發(fā)動(dòng)機(jī)和壓縮燃燒式發(fā)動(dòng)機(jī)。壓燃式天然氣發(fā)動(dòng)機(jī)大多由柴油機(jī)改裝而成,柴油-天然氣雙燃料發(fā)動(dòng)機(jī)是壓燃式天然氣發(fā)動(dòng)機(jī)的一種重要形式,在壓縮上止點(diǎn)附近,噴入氣缸內(nèi)的柴油自燃著火,并點(diǎn)燃?xì)飧變?nèi)的天然氣和空氣混合氣,使天然氣燃燒,因此柴油-天然氣雙燃料發(fā)動(dòng)機(jī)具有壓燃式和點(diǎn)燃式發(fā)動(dòng)機(jī)的雙重特點(diǎn)[3-5]。柴油-天然氣雙燃料發(fā)動(dòng)機(jī)大多采用原柴油機(jī)的壓縮比,可以充分發(fā)揮天然氣辛烷值高的優(yōu)點(diǎn),提高熱效率,此外柴油-天然氣雙燃料發(fā)動(dòng)機(jī)可以降低NOx和炭煙排放量,環(huán)保性好[6-7]。但由于柴油-天然氣雙燃料發(fā)動(dòng)機(jī)采用改裝前原柴油機(jī)的燃燒室,使得發(fā)動(dòng)機(jī)在雙燃料模式時(shí),存在一系列問(wèn)題,其中較為突出的問(wèn)題之一是柴油-天然氣雙燃料發(fā)動(dòng)機(jī)在低負(fù)荷時(shí),CH4排放量高,熱效率低,經(jīng)濟(jì)性差。

    當(dāng)前,國(guó)內(nèi)外學(xué)者對(duì)柴油-天然氣雙燃料發(fā)動(dòng)機(jī)進(jìn)行了廣泛的研究,包括天然氣的替代率、噴油時(shí)刻和轉(zhuǎn)速等對(duì)發(fā)動(dòng)機(jī)的影響。然而,關(guān)于柴油-天然氣雙燃料發(fā)動(dòng)機(jī)燃燒室的研究還相對(duì)較少。國(guó)內(nèi)外研究表明:燃燒室的形狀對(duì)燃燒過(guò)程有重要影響[8-10]。在已有文獻(xiàn)中,對(duì)于天然氣發(fā)動(dòng)機(jī)燃燒室的研究,更多的是針對(duì)天然氣單一燃料發(fā)動(dòng)機(jī)。如BijanYadollahi等利用STAR-CD仿真軟件模擬研究了5種燃燒室形狀對(duì)天然氣和空氣混合氣形成和分布的影響[11]。天津大學(xué)的焦運(yùn)景等針對(duì)火花點(diǎn)燃式天然氣發(fā)動(dòng)機(jī),設(shè)計(jì)了3種不同形狀的燃燒室,對(duì)不同燃燒室內(nèi)的氣流運(yùn)動(dòng)和燃燒過(guò)程進(jìn)行了數(shù)值模擬,并在試驗(yàn)臺(tái)架上進(jìn)行了試驗(yàn)驗(yàn)證[12]。也有部分學(xué)者對(duì)雙燃料發(fā)動(dòng)機(jī)燃燒室進(jìn)行了研究,但主要是對(duì)柴油-汽油或生物柴油-二甲醚等雙燃料發(fā)動(dòng)機(jī)燃燒室的研究,目前針對(duì)柴油-天然氣雙燃料發(fā)動(dòng)機(jī)燃燒室的詳細(xì)研究相對(duì)較少[13]。

    為解決低負(fù)荷時(shí)柴油-天然氣雙燃料發(fā)動(dòng)機(jī)熱效率低、CH4排放量高等問(wèn)題,本研究從改變?nèi)紵业男螤钪郑讦匦稳紵?改裝前原柴油機(jī)燃燒室),重新設(shè)計(jì)了圓柱形燃燒室和八邊啞鈴形燃燒室,利用商用軟件STAR-CD對(duì)ω形燃燒室、八邊啞鈴形燃燒室和圓柱形燃燒室進(jìn)行了數(shù)值模擬。為尋求適用于柴油-天然氣雙燃料發(fā)動(dòng)機(jī)最佳燃燒室提供數(shù)據(jù)支持。

    1 模型的建立和計(jì)算條件的選擇

    1.1 模型的建立

    試驗(yàn)所用柴油-天然氣雙燃料發(fā)動(dòng)機(jī)是由1臺(tái)電控高壓共軌、增壓中冷柴油機(jī)開(kāi)發(fā)而成,發(fā)動(dòng)機(jī)技術(shù)參數(shù)見(jiàn)表1。天然氣是在進(jìn)氣道噴射,從噴嘴噴出后,在氣道內(nèi)與空氣混合,天然氣與空氣混合得較為充分,進(jìn)入氣缸內(nèi)的混合氣趨于均勻。為了節(jié)省計(jì)算成本,本研究不考慮發(fā)動(dòng)機(jī)的進(jìn)氣階段和排氣階段,只考慮進(jìn)氣門(mén)關(guān)閉后和排氣門(mén)開(kāi)啟前的工作階段,即計(jì)算模型不含有進(jìn)排氣道。

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

    計(jì)算起始角度選為進(jìn)氣門(mén)關(guān)閉時(shí)刻,即壓縮上止點(diǎn)前120°曲軸轉(zhuǎn)角,結(jié)束角度選為壓縮上止點(diǎn)后70°,360°為壓縮上止點(diǎn)。活塞運(yùn)行至下止點(diǎn)時(shí),計(jì)算網(wǎng)格總數(shù)約為3×105個(gè),網(wǎng)格尺寸為0.8 mm。在研究過(guò)程中對(duì)網(wǎng)格進(jìn)行了相關(guān)性分析,研究結(jié)果表明,網(wǎng)格數(shù)對(duì)計(jì)算結(jié)果影響不大。圖1示出活塞在下止點(diǎn)時(shí)的計(jì)算網(wǎng)格。

    圖1 計(jì)算網(wǎng)格

    1.2 模型及算法的選擇

    由于κ-ε/RNG模型很好地考慮了近壁區(qū)域的流動(dòng),具有較高的模擬精度,因此本研究湍流模型選為κ-ε/RNG模型。Bai模型對(duì)干壁面和濕壁面采用了不同的處理方法,可以精確反映燃油粒子的撞壁運(yùn)動(dòng),因此本研究液滴撞壁模型選為Bai模型。霧化模型選為Reitz-Diwakar模型。STAR-CD仿真軟件提供了多種燃燒模型,包括CFM模型、ECFM模型、ECFM-3Z模型、PVM-MF模型、DARS模型和G-equation模型。其中ECFM-3Z燃燒模型主要由混合模型、火焰?zhèn)鞑ツP?、后火焰和排放模型、火花點(diǎn)火和自燃/爆震模型四部分組成,并且同時(shí)考慮了自燃、預(yù)混合燃燒和擴(kuò)散燃燒,可以模擬多種燃料同時(shí)燃燒,因此,本研究燃燒模型選為ECFM-3Z模型。由于ECFM-3Z模型內(nèi)置有著火模型和排放模型,因此著火模型選為Delay模型,NOx排放模型選為擴(kuò)展的澤爾多維奇機(jī)理,SOOT的排放模型為ERC模型。STAR-CD軟件提供了兩種求解算法,即SIMPLE算法和PISO算法,SIMPLE多用于穩(wěn)態(tài)計(jì)算,而PISO常用于瞬態(tài)計(jì)算,本研究選取PISO算法。

    1.3 模擬計(jì)算條件

    初始?jí)毫τ稍囼?yàn)得到的示功圖獲得,初始溫度由理想氣體狀態(tài)方程求得。初始流場(chǎng)以渦流比的形式給出,渦流比為1.2,壁面溫度邊界采用絕熱邊界。

    2 數(shù)值模擬結(jié)果

    2.1 模型驗(yàn)證

    為了驗(yàn)證模型的有效性,進(jìn)行了試驗(yàn)驗(yàn)證,所選工況轉(zhuǎn)速為1 335 r/min,扭矩為218 N·m,天然氣能量替代率為90%。圖2示出試驗(yàn)測(cè)得的缸內(nèi)壓力和模擬計(jì)算得出的缸內(nèi)壓力對(duì)比,圖3示出試驗(yàn)測(cè)得的CH4排放量和模擬計(jì)算得出的CH4排放量的對(duì)比。從缸內(nèi)壓力曲線和CH4排放量柱狀圖可以看出,試驗(yàn)結(jié)果和模擬結(jié)果吻合較好,說(shuō)明所選模型和計(jì)算方法是合理的,可以進(jìn)行下一步研究。

    圖2 試驗(yàn)和模擬缸內(nèi)壓力的對(duì)比

    圖3 試驗(yàn)和模擬CH4排放量的對(duì)比

    天然氣能量替代率計(jì)算公式為

    (1)

    2.2 燃燒室形狀影響規(guī)律

    對(duì)ω形燃燒室、八邊啞鈴形燃燒室和圓柱形燃燒室進(jìn)行了數(shù)值模擬,以探究燃燒室的形狀對(duì)柴油-天然氣雙燃料發(fā)動(dòng)機(jī)燃燒過(guò)程的影響規(guī)律。模擬計(jì)算除燃燒室形狀不同外,其他邊界條件均一致,壓縮比保持為17.2∶1。燃燒室形狀見(jiàn)圖4。

    圖4 3種燃燒室形狀示意

    圖5示出3種燃燒室的缸內(nèi)平均壓力、平均溫度和放熱率隨曲軸轉(zhuǎn)角變化的模擬計(jì)算結(jié)果。從圖5a和圖5b可以看出,在上止點(diǎn)附近,著火初期,ω形燃燒室缸內(nèi)平均壓力和平均溫度最高,八邊啞鈴形燃燒室次之,圓柱形燃燒室最低。隨著活塞逐漸下行,八邊啞鈴形燃燒室的缸內(nèi)平均壓力和平均溫度逐漸變?yōu)樽罡撸瑘A柱形燃燒室次之,ω形燃燒室最低。從圖5c可以看出,在上止點(diǎn)附近,ω形燃燒室的第一放熱率峰值最高,其次是八邊啞鈴形燃燒室,圓柱形燃燒室第一放熱率峰值最低。同時(shí),八邊啞鈴形燃燒室第二放熱率峰值高于圓柱形燃燒室,ω形燃燒室第二放熱率峰值不明顯。

    圖5 燃燒室形狀對(duì)缸內(nèi)參數(shù)的影響

    ω形燃燒室由于氣流速度及湍動(dòng)能較小,有利于火核的穩(wěn)定形成,使得滯燃期較短。同時(shí)柴油分布區(qū)域相對(duì)較小,區(qū)域內(nèi)柴油混合氣相對(duì)較濃,一經(jīng)著火便迅速燃燒,并快速引燃柴油周?chē)烊粴猓變?nèi)壓力和溫度快速上升。而八邊啞鈴形燃燒室由于較強(qiáng)的擠流作用,氣流速度較大,形成的火核易被吹熄,使得柴油滯燃期相對(duì)較長(zhǎng),因此,當(dāng)活塞在上止點(diǎn)附近時(shí),八邊啞鈴形燃燒室內(nèi)溫度及壓力較ω形燃燒室低。隨著活塞遠(yuǎn)離上止點(diǎn),八邊啞鈴形燃燒室內(nèi)由于著火點(diǎn)分布廣泛,湍動(dòng)能較大,燃燒速度加快,更多的天然氣參與燃燒,放熱量增多,缸內(nèi)壓力和溫度較高。而ω形燃燒室由于燃燒初期便大量燃燒放熱,著火點(diǎn)分布相對(duì)較少,不利于更多的天然氣參與燃燒,氣缸周?chē)腿紵业撞坑写罅刻烊粴馐S?,缸?nèi)壓力和溫度相對(duì)較低,第二放熱率峰值不明顯。圓柱形燃燒室燃燒放熱情況介于八邊啞鈴形和ω形之間。

    圖6示出燃料剩余百分比隨曲軸轉(zhuǎn)角的變化,圖7示出3種燃燒室指示熱效率對(duì)比。在燃料剩余百分比曲線圖中,柴油剩余百分比定義為當(dāng)前時(shí)刻氣缸內(nèi)剩余的柴油質(zhì)量與到當(dāng)前時(shí)刻為止累計(jì)噴入氣缸內(nèi)柴油總質(zhì)量的比值,天然氣剩余百分比定義為當(dāng)前時(shí)刻氣缸內(nèi)剩余的天然氣質(zhì)量與到當(dāng)前時(shí)刻為止累計(jì)進(jìn)入氣缸內(nèi)天然氣總質(zhì)量的比值。從圖6可以看出,ω形燃燒室內(nèi)天然氣剩余量最多,圓柱形燃燒室次之,八邊啞鈴形燃燒室內(nèi)天然氣剩余量最少;同時(shí)可以看出,ω形燃燒室內(nèi)柴油剩余量最低,圓柱形燃燒室內(nèi)柴油剩余量最高。這驗(yàn)證了八邊形燃燒室更有利于天然氣的燃燒。從圖7可以看出,八邊啞鈴形燃燒指示熱效率最高,約為45.3%,圓柱形燃燒室次之,ω形燃燒室最低。

    圖6 燃料剩余百分比隨曲軸轉(zhuǎn)角的變化

    圖7 3種燃燒室指示熱效率對(duì)比

    2.3 微觀場(chǎng)分析

    2.3.1 速度場(chǎng)分析

    圖8示出3種燃燒室內(nèi)氣流速度在不同曲軸轉(zhuǎn)角時(shí)的微觀場(chǎng),為了讀圖方便,對(duì)每一微觀場(chǎng)采用了固定的標(biāo)尺,徑向取圖位置為過(guò)活塞頂下部2 mm處,軸向取圖位置恰好經(jīng)過(guò)兩個(gè)噴油嘴中心連線,如圖中實(shí)線。從圖中可以看出,當(dāng)活塞運(yùn)行至上止點(diǎn)附近時(shí),八邊啞鈴形燃燒室內(nèi)氣流速度最大,并在燃燒室底部仍具有較強(qiáng)的氣流速度,這是由于八邊啞鈴形燃燒室喉口直徑最小,產(chǎn)生的擠流最強(qiáng)??梢钥闯?,燃燒室喉口直徑對(duì)擠流強(qiáng)度影響很大。較強(qiáng)的擠流有利于增加湍動(dòng)能(這從對(duì)應(yīng)的湍動(dòng)能場(chǎng)可以看出),加快燃燒室內(nèi)已燃?xì)怏w與未燃?xì)怏w混合。天津大學(xué)的焦運(yùn)景和江蘇大學(xué)的潘劍鋒等在燃燒室研究中也得出了類(lèi)似觀點(diǎn)[12-14]。從圓柱形燃燒室底部速度場(chǎng)分布可以看出,燃燒室底部速度尤其是燃燒室底部的中心部位氣流速度較小。八邊啞鈴形燃燒室和ω形燃燒室由于燃燒室底部設(shè)計(jì)有凸臺(tái),對(duì)燃燒室內(nèi)氣流運(yùn)動(dòng)有引導(dǎo)作用,并且八邊啞鈴形燃燒室內(nèi)凸臺(tái)導(dǎo)流作用更明顯。燃燒室底部剩余的天然氣量與氣流運(yùn)動(dòng)和底部所存在的天然氣多少有關(guān),在燃燒室底部設(shè)計(jì)凸臺(tái),一方面減小了氣流速度較小區(qū)域,另一方面可以引導(dǎo)燃燒室內(nèi)氣流運(yùn)動(dòng),進(jìn)而影響燃料與空氣的混合、柴油的分布和火焰的傳播,減少燃燒室底部燃料剩余。日本德島大學(xué)的Y.Kinoguchi等也在研究中得出,較小的縮口并結(jié)合燃燒室底部凸臺(tái)能延長(zhǎng)燃燒室內(nèi)較高的湍動(dòng)能持續(xù)期,有效促進(jìn)燃料的燃燒[15]。

    2.3.2 湍動(dòng)能場(chǎng)分析

    圖9示出3種燃燒室內(nèi)湍動(dòng)能在不同曲軸轉(zhuǎn)角時(shí)的微觀場(chǎng)。從圖中可以看出,在壓縮過(guò)程后期,活塞接近上止點(diǎn)時(shí),八邊啞鈴形燃燒室湍動(dòng)能最大且分布最廣,圓柱形燃燒室次之,ω形燃燒室內(nèi)湍動(dòng)能最小。在燃燒初期,較高的湍動(dòng)能有利于天然氣和空氣混合得更加均勻,以及柴油在燃燒室內(nèi)分布得更加廣泛,有利于增加著火面積,但較高的湍動(dòng)能不利于形成穩(wěn)定的火核。因此,在燃燒初期,柴油的滯燃期相對(duì)較長(zhǎng)。隨著活塞下行,逐漸遠(yuǎn)離上止點(diǎn),八邊啞鈴形燃燒室內(nèi)高湍動(dòng)能區(qū)域進(jìn)一步增大,并在活塞頂和氣缸蓋之間保持較高的湍動(dòng)能,這不僅有利于天然氣與空氣的快速混合,而且有利于火焰的傳播,加快天然氣的燃燒,提高天然氣的利用率??梢钥闯靓匦稳紵覛飧變?nèi)部湍動(dòng)能相對(duì)較小,而圓柱形燃燒室氣缸內(nèi)湍動(dòng)能介于兩者之間。

    2.3.3 柴油濃度場(chǎng)分析

    圖10示出3種燃燒室內(nèi)柴油質(zhì)量分?jǐn)?shù)在不同曲軸轉(zhuǎn)角時(shí)的微觀場(chǎng)。柴油噴入燃燒室后,經(jīng)過(guò)著火前的物理和化學(xué)反應(yīng),在高溫高壓條件下自燃著火。從圖中可以看出,隨著活塞上移,ω形燃燒室內(nèi)柴油逐漸向燃燒室壁面方向擴(kuò)散,而八邊啞鈴形燃燒室內(nèi)柴油向凸臺(tái)方向以及沿著凸臺(tái)壁面向下擴(kuò)散,而圓柱形燃燒室內(nèi)柴油則在噴射位置處向四周擴(kuò)散,3種燃燒室內(nèi)柴油的擴(kuò)散方向主要是由燃燒室的形狀和氣流運(yùn)動(dòng)方向所決定。與八邊啞鈴形燃燒室和圓柱形燃燒室相比,ω形燃燒室內(nèi)柴油剩余較少,這是因?yàn)棣匦稳紵覂?nèi)柴油滯燃期較短,燃燒速度較快,八邊啞鈴形燃燒室和圓柱形燃燒室內(nèi)柴油滯燃期相對(duì)較長(zhǎng),燃燒速度較慢,在活塞到達(dá)上止點(diǎn)時(shí),仍有較多的柴油剩余。柴油燃燒后放出大量的熱量,使得缸內(nèi)溫度升高,有利于缸內(nèi)燃料的燃燒。隨著活塞下行,柴油繼續(xù)燃燒,可以看出,當(dāng)活塞運(yùn)行到365°時(shí),ω形燃燒室內(nèi)柴油消耗殆盡,而圓柱形燃燒室內(nèi)仍有少量的柴油剩余。同時(shí),從圖10中柴油的擴(kuò)散方向可以得出,在燃燒室的底部設(shè)計(jì)凸臺(tái),通過(guò)引導(dǎo)氣流運(yùn)動(dòng)可以引導(dǎo)柴油擴(kuò)散方向,改善柴油在燃燒室內(nèi)的分布。

    圖8 速度場(chǎng)

    圖9 湍動(dòng)能場(chǎng)

    2.3.4 天然氣濃度場(chǎng)分析

    圖11示出3種燃燒室內(nèi)天然氣質(zhì)量分?jǐn)?shù)在不同曲軸轉(zhuǎn)角時(shí)的微觀場(chǎng)。從圖中可以看出,在柴油著火位置,天然氣濃度開(kāi)始下降,這說(shuō)明天然氣是由柴油引燃的。當(dāng)活塞上行到達(dá)壓縮上止點(diǎn)時(shí),在同一曲軸轉(zhuǎn)角下,ω形燃燒室內(nèi)天然氣濃度下降最多,八邊啞鈴形燃燒室次之,圓柱形燃燒室內(nèi)天然氣濃度下降最少,這是由于ω形燃燒室內(nèi)柴油滯燃期較短,燃燒速度較快,并不斷引燃周?chē)烊粴猓固烊粴飧绺斓厝紵?。隨著活塞離開(kāi)壓縮上止點(diǎn)下行,由于八邊啞鈴形燃燒室和圓柱形燃燒室內(nèi)柴油滯燃期較長(zhǎng),氣缸內(nèi)湍動(dòng)能較大,使得柴油在缸內(nèi)分布更加廣泛,有利于引燃更多的天然氣??梢钥闯?,當(dāng)活塞運(yùn)行至400°時(shí),八邊啞鈴形燃燒室內(nèi)天然氣剩余最少,僅在燃燒室底部靠近凸臺(tái)部位有少量剩余,這是由于八邊啞鈴形燃燒室內(nèi)柴油分布廣泛,并且湍流較大,促進(jìn)了已燃?xì)怏w和未燃?xì)怏w的交換,提高了火焰?zhèn)鞑ニ俾剩龠M(jìn)了燃燒室內(nèi)燃料的充分燃燒。ω形燃燒室和圓柱形燃燒室底部和靠近氣缸軸線的區(qū)域有大量的天然氣剩余,由于圓柱形燃燒室氣缸內(nèi)接近氣缸壁處湍動(dòng)能相對(duì)ω形燃燒室要大,燃燒速率較ω形燃燒室要快,天然氣燃燒更為充分,天然氣剩余量較ω形燃燒室要少。

    2.3.5 溫度場(chǎng)分析

    圖12示出3種燃燒室內(nèi)溫度在不同曲軸轉(zhuǎn)角時(shí)的微觀場(chǎng)。從圖中可以看出,柴油開(kāi)始著火前,缸內(nèi)溫度分布均勻,僅在柴油噴入和貫穿的區(qū)域溫度下降。燃燒之后,在同一曲軸轉(zhuǎn)角下,ω形燃燒室內(nèi)已燃區(qū)域內(nèi)的溫度最高,八邊啞鈴形燃燒室次之,圓柱形燃燒室最低,這是由于ω形燃燒室內(nèi)柴油滯燃期最短。缸內(nèi)柴油和天然氣燃燒后,釋放出大量的熱量,使得缸內(nèi)高溫區(qū)域不斷擴(kuò)大,逐漸向壁面?zhèn)鞑ァ?梢钥闯觯?80°時(shí),3種燃燒室底部的溫度均較低。400°時(shí),對(duì)于ω形燃燒室和圓柱形燃燒室,燃燒室底部和靠近氣缸軸線的區(qū)域溫度均較低,而八邊啞鈴形燃燒室內(nèi)溫度分布則較為均勻,這是由于八邊啞鈴形燃燒室內(nèi)柴油分布廣泛,同時(shí)氣缸內(nèi)氣流速度快,湍動(dòng)能大,有利于天然氣的充分燃燒,放出更多的熱量,因此八邊啞鈴形燃燒室內(nèi)溫度相對(duì)較高。

    圖11 天然氣濃度場(chǎng)

    圖12 溫度場(chǎng)

    2.4 燃燒室形狀對(duì)CH4排放的影響分析

    圖13示出3種燃燒室所對(duì)應(yīng)的CH4排放量(CH4排放量是指計(jì)算結(jié)束時(shí)刻430°時(shí)CH4的剩余量)。從圖中可以看出,ω形燃燒室CH4排放量最高,八邊啞鈴形燃燒室CH4排放量最低。這是由于八邊啞鈴形燃燒室內(nèi)天然氣著火點(diǎn)分布較為廣泛,湍動(dòng)能較大,天然氣燃燒充分,因此排放出的CH4最少。這從上文中燃料剩余比例曲線以及缸內(nèi)天然氣微觀場(chǎng)分析也能看出。

    圖13 3種燃燒室內(nèi)CH4排放量

    3 結(jié)論

    a) 八邊啞鈴形燃燒室缸內(nèi)平均壓力和溫度最高,圓柱形燃燒室次之,ω形燃燒室最低;同時(shí),ω形燃燒室內(nèi)柴油滯燃期較短,天然氣燃燒速度較慢,八邊啞鈴形燃燒室和圓柱形燃燒室內(nèi)柴油滯燃期相對(duì)較長(zhǎng),天然氣燃燒速度較快;八邊啞鈴形燃燒室內(nèi),天然氣剩余比例最小,天然氣利用率高,燃料的經(jīng)濟(jì)性較好;此外,八邊啞鈴形燃燒室內(nèi)指示熱效率最高,為45.3%;

    b) 從模擬計(jì)算微觀場(chǎng)中可以看出,八邊啞鈴形燃燒室內(nèi)氣流運(yùn)動(dòng)速度最大,湍動(dòng)能最強(qiáng),圓柱形燃燒室次之,ω形燃燒室最低;ω形燃燒室的CH4排放量最多,其中CH4主要是在燃燒室的底部和氣缸壁附近剩余較多;八邊啞鈴形燃燒室CH4排放量最低,僅在燃燒室底部,貼近凸臺(tái)處剩余少量CH4;對(duì)于柴油-天然氣雙燃料發(fā)動(dòng)機(jī),為保證天然氣的可靠著火,在天然氣著火初期,湍動(dòng)能不應(yīng)過(guò)大;天然氣快速燃燒時(shí),可增加湍動(dòng)能強(qiáng)度和延長(zhǎng)湍動(dòng)能持續(xù)期,以保證天然氣的快速、充分燃燒,提高柴油-天然氣的熱效率,降低CH4排放;

    c) 通過(guò)減小燃燒室喉口直徑,可以增強(qiáng)燃燒室內(nèi)的擠流強(qiáng)度,增加氣缸內(nèi)的湍動(dòng)能,提高缸內(nèi)燃料的燃燒速度;在燃燒室的底部中心部位設(shè)計(jì)凸臺(tái),能引導(dǎo)燃燒室內(nèi)的氣流運(yùn)動(dòng),擴(kuò)大著火點(diǎn)分布范圍,改善燃燒室底部燃料的燃燒,并在較大擠流作用下凸臺(tái)的導(dǎo)流效果更加明顯。

    [1] 宋均,張武高,黃震.車(chē)用天然氣發(fā)動(dòng)機(jī)技術(shù)與性能研究[J].車(chē)用發(fā)動(dòng)機(jī),2001(6):8-12.

    [2] 徐富水,秦江濤.CNG -柴油雙燃料電控發(fā)動(dòng)機(jī)試驗(yàn)研究[J].車(chē)用發(fā)動(dòng)機(jī),2012(1):60-62.

    [3] 周龍保,劉忠長(zhǎng),高宗英.內(nèi)燃機(jī)學(xué)[M].北京:機(jī)械工業(yè)出版社,2010.

    [4] Papagiannakis R G,Hountalas D T.Experimental and Theoretical Analysis of the Combustion and Pollutants Formation Mechanisms in Dual Fuel DI Diesel Engines [C].SAE Paper 2005-01-1726.

    [5] 林志強(qiáng),蘇萬(wàn)華.柴油引燃式天然氣發(fā)動(dòng)機(jī)最佳引燃柴油量及過(guò)量空氣系數(shù)濃限、稀限的研究[J].內(nèi)燃機(jī)學(xué)報(bào),2002,20(6):505-510.

    [6] Kusaka J, Okamoto T,Daisho Y,et al.Combustion and exhaust gas emission characteristics of a diesel engine dual-fuelled with natural gas[C].SAE Review October, 2000, 21(4):489-496.

    [7] 金華玉,劉忠長(zhǎng),王忠恕,等.柴油機(jī)燃燒過(guò)程模擬分析[J].吉林大學(xué)學(xué)報(bào)(工學(xué)版), 2007,37(5):1028-1033.

    [8] Stefano Cordiner,Michele Gambino,SabatoIannaccone.Numerical and Experimental Analysis of Combustion and Exhaust[J].Energy & Fuels, 2008,22:1418-1424.

    [9] Antonio Paolo Carlucci,Domenico Laforgia, Roberto Saracino,et al.Study of Combustion Development in Methane-Diesel Dual Fuel Engines, Based on the Analysis of In-Cylinder Luminance[C].SAE Paper 2010-01-1297.

    [10] Splitter D A,Wissink M, Kokjohn S L, et al. Effect of compression ratio and piston geometry on RCCI load limits and efficiency[C].SAE Paper 2012-01-0383.

    [11] BijanYadollahi,Masound Boroomand.The Effect of Piston Head Geometry on Natural Gas Direct Injection and Mixture Formation in a SIEngine with Centrally Mounted Single-Hole Injector[C].SAE Paper 2011-01-2448.

    [12] 焦運(yùn)景,董宏,張惠明,等.燃燒室形狀對(duì)天然氣發(fā)動(dòng)機(jī)燃燒過(guò)程影響的研究[J].內(nèi)燃機(jī)工程,2008,30(4):28-33.

    [13] 杜潛,袁霞.燃燒室結(jié)構(gòu)對(duì)XN2100雙燃料發(fā)動(dòng)機(jī)性能的影響[J].河南機(jī)電高等專(zhuān)科學(xué)校學(xué)報(bào),2007,15(2):83-85.[14] 潘劍鋒,盧清波.柴油機(jī)燃燒過(guò)程的數(shù)值模擬及燃燒室改進(jìn)[J].江蘇大學(xué)學(xué)報(bào), 2012,33(4):390-395.

    [15] Kidoguchi Y,Sanda M,Miwa K.Experimental and Theoretical Optimization of Combustion Chamber and Fuel Distribution for the Low Emission Direct-Injection Diesel Engine[J]. Journal of Engineering for Gas Turbines and Power,2003,125:351-357.

    [編輯: 李建新]

    Simulation Analysis of Influence of Combustion Chamber Shape on Dual-fuel Engine Performance

    WANG Zhongshu, LI Hui, YANG Zhou, LI Weifeng

    (State Key Laboratory of Automobile Simulation and Control,Jilin University, Changchun 130025, China)

    The influences of three kinds of combustion chamber shape on the performance of diesel and natural gas dual-fuel engine were studied with STAR-CD software. The three kinds of combustion chamber were ω-shaped combustion chamber, octagonal dumbbell-shaped combustion chamber and cylinder-shaped combustion chamber. The results show that the octagonal dumbbell-shaped combustion chamber has better performance. The in-cylinder turbulent kinetic energy strengthens, the flame propagation accelerates and the fuel utilization ratio increases due to the squish flow caused by reducing dumbbell throat diameter. Meanwhile, the designed boss at the bottom of the combustion chamber can guide the mixture movement, and promote diesel diffusion to the chamber bottom, resulting a wide distribution of ignition point. Therefore, the in-cylinder average pressure,average temperature and indicated thermal efficiency of octagonal dumbbell-shaped combustion chamber is the highest and the proportion of remained natural gas is the least.

    dual-fuel engine; combustion chamber; combustion; emission; simulation

    2015-06-23;

    2015-09-16

    國(guó)家重大基礎(chǔ)研究發(fā)展規(guī)劃(“973”)資助項(xiàng)目(2013CB228402);國(guó)家自然科學(xué)基金資助項(xiàng)目(50906033);吉林省科技發(fā)展計(jì)劃資助項(xiàng)目(20126009);吉林大學(xué)研究生創(chuàng)新基金資助項(xiàng)目(2015092)

    王忠恕(1977—),男,教授,博士生導(dǎo)師,主要研究方向?yàn)閮?nèi)燃機(jī)燃燒與排放控制;wangzhongshu@jlu.edu.cn。

    李慧(1988—),女,碩士,主要研究方向?yàn)閮?nèi)燃機(jī)燃燒與排放控制;lhui157@163.com。

    10.3969/j.issn.1001-2222.2016.01.001

    TK422

    B

    1001-2222(2016)01-0001-08

    猜你喜歡
    圓柱形雙燃料啞鈴
    基于修正函數(shù)法的雙燃料散貨船船型優(yōu)化
    樹(shù)干為什么要長(zhǎng)成圓柱形
    雙燃料主機(jī)缸套水冷卻系統(tǒng)優(yōu)化設(shè)計(jì)
    我給爸爸當(dāng)“啞鈴”
    樹(shù)干為什么是圓柱形的
    橫臥啞鈴形Rathke囊腫1例
    79800 DWT雙燃料散貨船LNG儲(chǔ)氣罐位置確定
    船海工程(2015年4期)2016-01-05 15:53:50
    去贅肉又強(qiáng)身的啞鈴操(上)
    去贅肉又強(qiáng)身的啞鈴操(上)
    重載鐵路圓柱形橋墩橫向振動(dòng)試驗(yàn)研究
    成人毛片a级毛片在线播放| 成年版毛片免费区| av女优亚洲男人天堂| 99国产精品一区二区蜜桃av| 亚洲av中文字字幕乱码综合| 九九久久精品国产亚洲av麻豆| 亚洲av不卡在线观看| 亚洲av成人av| 麻豆久久精品国产亚洲av| 99热网站在线观看| 丝袜美腿在线中文| 精品一区二区免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 级片在线观看| 少妇人妻精品综合一区二区 | 美女被艹到高潮喷水动态| 综合色丁香网| 亚洲成人久久性| 欧美高清成人免费视频www| 久久午夜福利片| 亚洲精品日韩av片在线观看| 亚洲av五月六月丁香网| 在线播放国产精品三级| 高清毛片免费观看视频网站| 久久精品国产自在天天线| 国产亚洲精品久久久久久毛片| 又黄又爽又刺激的免费视频.| 国产亚洲91精品色在线| 亚洲性久久影院| 啦啦啦韩国在线观看视频| 黄色日韩在线| 成人av在线播放网站| 久久午夜亚洲精品久久| 99久国产av精品国产电影| 色哟哟·www| 午夜福利高清视频| 欧美高清成人免费视频www| 午夜日韩欧美国产| 亚洲成人久久爱视频| 在线观看免费视频日本深夜| 国产真实乱freesex| 精品午夜福利在线看| 男人舔女人下体高潮全视频| 成人精品一区二区免费| 最后的刺客免费高清国语| a级毛片免费高清观看在线播放| 国产视频内射| 寂寞人妻少妇视频99o| 日本一本二区三区精品| 亚洲天堂国产精品一区在线| 麻豆乱淫一区二区| 九九热线精品视视频播放| 变态另类成人亚洲欧美熟女| 看片在线看免费视频| 国内久久婷婷六月综合欲色啪| 直男gayav资源| 变态另类丝袜制服| 一区二区三区高清视频在线| 国产亚洲精品av在线| 久久久久性生活片| 国产亚洲av嫩草精品影院| 欧美性猛交黑人性爽| 深夜精品福利| 色5月婷婷丁香| 日韩一区二区视频免费看| 狂野欧美白嫩少妇大欣赏| 久久久久久久午夜电影| 18禁黄网站禁片免费观看直播| 欧美三级亚洲精品| 国产成人一区二区在线| 欧美一级a爱片免费观看看| 特大巨黑吊av在线直播| 99热全是精品| 欧美日韩乱码在线| 女人十人毛片免费观看3o分钟| 亚洲美女视频黄频| 在线观看免费视频日本深夜| 久久久久国产网址| 成年版毛片免费区| 欧美3d第一页| 成人无遮挡网站| 最新中文字幕久久久久| 日韩欧美精品v在线| 18禁黄网站禁片免费观看直播| 成人三级黄色视频| 精品不卡国产一区二区三区| 亚洲欧美日韩东京热| 国产亚洲精品久久久久久毛片| 永久网站在线| 久久欧美精品欧美久久欧美| 九九在线视频观看精品| 麻豆av噜噜一区二区三区| av在线天堂中文字幕| 久久久久久大精品| 亚洲av一区综合| 99国产精品一区二区蜜桃av| 国产一区二区三区av在线 | 亚洲va在线va天堂va国产| 亚洲最大成人手机在线| 亚洲av成人精品一区久久| 日韩欧美免费精品| 国产精品av视频在线免费观看| 亚洲自偷自拍三级| 精品久久国产蜜桃| ponron亚洲| 亚洲经典国产精华液单| 日韩欧美在线乱码| 啦啦啦观看免费观看视频高清| 乱码一卡2卡4卡精品| 精品国产三级普通话版| 色视频www国产| 亚洲成人中文字幕在线播放| 色哟哟哟哟哟哟| 永久网站在线| 欧美潮喷喷水| 99久久中文字幕三级久久日本| 日日啪夜夜撸| 久久精品人妻少妇| 麻豆久久精品国产亚洲av| 国产免费一级a男人的天堂| 国产精品久久久久久久电影| 亚洲激情五月婷婷啪啪| 国产成人一区二区在线| 啦啦啦观看免费观看视频高清| 别揉我奶头~嗯~啊~动态视频| 亚洲精品影视一区二区三区av| 91在线精品国自产拍蜜月| 日韩强制内射视频| 一a级毛片在线观看| 99在线人妻在线中文字幕| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲欧美成人精品一区二区| 女人被狂操c到高潮| 久久这里只有精品中国| av中文乱码字幕在线| 91午夜精品亚洲一区二区三区| 插逼视频在线观看| 高清午夜精品一区二区三区 | 天堂av国产一区二区熟女人妻| 免费看日本二区| 在线免费观看不下载黄p国产| 中文字幕精品亚洲无线码一区| 精品一区二区免费观看| 国产午夜精品久久久久久一区二区三区 | av在线播放精品| 亚洲在线观看片| 国产欧美日韩精品亚洲av| 97超视频在线观看视频| 国产精品,欧美在线| 久久精品国产清高在天天线| 精品一区二区免费观看| 一级a爱片免费观看的视频| 人人妻人人澡欧美一区二区| 日本五十路高清| av中文乱码字幕在线| 波多野结衣高清无吗| 此物有八面人人有两片| 久久热精品热| 此物有八面人人有两片| 午夜福利在线在线| 午夜影院日韩av| 亚洲人成网站在线观看播放| 国产高潮美女av| 国产精品国产高清国产av| 99在线视频只有这里精品首页| 我要搜黄色片| 精品欧美国产一区二区三| 国产单亲对白刺激| 女的被弄到高潮叫床怎么办| 麻豆av噜噜一区二区三区| 晚上一个人看的免费电影| 精品一区二区三区av网在线观看| 国产伦在线观看视频一区| 99热这里只有是精品50| 日韩成人av中文字幕在线观看 | 国产精品无大码| 99久久中文字幕三级久久日本| 国产极品精品免费视频能看的| 亚洲av二区三区四区| av在线蜜桃| 国产一区二区在线观看日韩| 成年女人看的毛片在线观看| 精品久久久久久久久亚洲| 亚洲精品粉嫩美女一区| 我要搜黄色片| 18禁裸乳无遮挡免费网站照片| 欧美xxxx性猛交bbbb| 免费看av在线观看网站| 极品教师在线视频| 黄色欧美视频在线观看| 搡老熟女国产l中国老女人| 国产久久久一区二区三区| 国产久久久一区二区三区| 国内精品久久久久精免费| 国产精品久久电影中文字幕| 男人和女人高潮做爰伦理| 三级男女做爰猛烈吃奶摸视频| 蜜桃亚洲精品一区二区三区| 国产淫片久久久久久久久| 夜夜夜夜夜久久久久| 偷拍熟女少妇极品色| 一夜夜www| 嫩草影院入口| 禁无遮挡网站| 精品无人区乱码1区二区| 人妻久久中文字幕网| 中文在线观看免费www的网站| 在线播放无遮挡| 成人综合一区亚洲| 欧美日韩国产亚洲二区| 午夜精品一区二区三区免费看| 精品人妻一区二区三区麻豆 | 亚洲av成人av| 精品人妻偷拍中文字幕| 九九在线视频观看精品| 99久久久亚洲精品蜜臀av| 亚洲精品一区av在线观看| 97人妻精品一区二区三区麻豆| 日本与韩国留学比较| 亚洲精品在线观看二区| 国产 一区精品| 亚洲内射少妇av| 国产一区二区激情短视频| 精品久久久久久久久亚洲| 日本成人三级电影网站| av在线亚洲专区| 此物有八面人人有两片| 精品乱码久久久久久99久播| 99国产精品一区二区蜜桃av| 99riav亚洲国产免费| 午夜免费激情av| 精品99又大又爽又粗少妇毛片| 日韩在线高清观看一区二区三区| 国产极品精品免费视频能看的| 美女黄网站色视频| 欧美一区二区国产精品久久精品| 变态另类成人亚洲欧美熟女| 别揉我奶头 嗯啊视频| 一本精品99久久精品77| 黄色配什么色好看| 亚洲乱码一区二区免费版| 女生性感内裤真人,穿戴方法视频| 一边摸一边抽搐一进一小说| 欧美色欧美亚洲另类二区| 禁无遮挡网站| 女的被弄到高潮叫床怎么办| 99热精品在线国产| 人人妻人人澡人人爽人人夜夜 | 午夜激情福利司机影院| 久久久久久久久久久丰满| 亚洲在线自拍视频| 日韩国内少妇激情av| 日本黄大片高清| 亚洲自拍偷在线| 久久久成人免费电影| 69人妻影院| 国产视频内射| 国产精品人妻久久久久久| 悠悠久久av| 又粗又爽又猛毛片免费看| 伦精品一区二区三区| 国产黄片美女视频| 两个人的视频大全免费| 99热只有精品国产| 亚洲最大成人手机在线| 久久久久精品国产欧美久久久| a级一级毛片免费在线观看| 又黄又爽又刺激的免费视频.| 91在线观看av| 中文亚洲av片在线观看爽| 自拍偷自拍亚洲精品老妇| 精品乱码久久久久久99久播| 一个人免费在线观看电影| 能在线免费观看的黄片| 日本黄色片子视频| 国产乱人偷精品视频| 两性午夜刺激爽爽歪歪视频在线观看| 久久久午夜欧美精品| 人妻夜夜爽99麻豆av| 成人特级黄色片久久久久久久| 欧美又色又爽又黄视频| av免费在线看不卡| 国产黄色小视频在线观看| 男人和女人高潮做爰伦理| 特大巨黑吊av在线直播| 亚洲丝袜综合中文字幕| 亚洲经典国产精华液单| 波多野结衣高清无吗| 免费av毛片视频| 乱人视频在线观看| 久久国内精品自在自线图片| 国产精品久久久久久久久免| 日韩大尺度精品在线看网址| 伦理电影大哥的女人| 国产精品久久久久久亚洲av鲁大| 丰满的人妻完整版| 超碰av人人做人人爽久久| 国产淫片久久久久久久久| 亚洲精品乱码久久久v下载方式| 午夜视频国产福利| 小说图片视频综合网站| 成人特级黄色片久久久久久久| 一级a爱片免费观看的视频| 国产一区二区在线av高清观看| 欧美在线一区亚洲| 亚洲精品久久国产高清桃花| 久久精品夜色国产| 99久久精品一区二区三区| 国产高潮美女av| 国产一区二区激情短视频| 日韩 亚洲 欧美在线| 麻豆av噜噜一区二区三区| 又黄又爽又刺激的免费视频.| 99热精品在线国产| 午夜久久久久精精品| 亚洲三级黄色毛片| 又黄又爽又免费观看的视频| 成人国产麻豆网| av视频在线观看入口| 99热这里只有精品一区| 国语自产精品视频在线第100页| 日日啪夜夜撸| 成年女人永久免费观看视频| 最近的中文字幕免费完整| 国产在线男女| 午夜影院日韩av| 成人午夜高清在线视频| 精品日产1卡2卡| 免费看光身美女| 久久久久久久久中文| 成人高潮视频无遮挡免费网站| 日韩欧美精品v在线| 亚洲美女视频黄频| 久久久久九九精品影院| 免费人成在线观看视频色| 2021天堂中文幕一二区在线观| 99九九线精品视频在线观看视频| 国产成人aa在线观看| 女的被弄到高潮叫床怎么办| 两性午夜刺激爽爽歪歪视频在线观看| av中文乱码字幕在线| 我要看日韩黄色一级片| 国产日本99.免费观看| 亚洲精华国产精华液的使用体验 | 久久久久国产精品人妻aⅴ院| 午夜a级毛片| 黄色视频,在线免费观看| 秋霞在线观看毛片| 免费一级毛片在线播放高清视频| 国产男靠女视频免费网站| 激情 狠狠 欧美| 亚洲中文日韩欧美视频| 国产乱人视频| 综合色丁香网| 直男gayav资源| 寂寞人妻少妇视频99o| 国产亚洲91精品色在线| 一级毛片我不卡| 如何舔出高潮| 日本熟妇午夜| 亚洲av中文av极速乱| 免费看美女性在线毛片视频| 欧美另类亚洲清纯唯美| 精品熟女少妇av免费看| 精华霜和精华液先用哪个| 一进一出好大好爽视频| 国产成年人精品一区二区| 99久国产av精品国产电影| 日韩大尺度精品在线看网址| 12—13女人毛片做爰片一| 深夜精品福利| 麻豆成人午夜福利视频| 99久久精品一区二区三区| 日韩精品青青久久久久久| 欧美性猛交╳xxx乱大交人| 给我免费播放毛片高清在线观看| 国产在线精品亚洲第一网站| 国产极品精品免费视频能看的| 淫秽高清视频在线观看| 非洲黑人性xxxx精品又粗又长| 在线播放国产精品三级| 国产高清不卡午夜福利| 国产精品一及| 日韩欧美精品免费久久| 国产aⅴ精品一区二区三区波| 亚洲国产高清在线一区二区三| 亚洲av免费在线观看| 在现免费观看毛片| 精品国内亚洲2022精品成人| 久久热精品热| 最近视频中文字幕2019在线8| 最新在线观看一区二区三区| 久久久久久久久久成人| 久久国产乱子免费精品| 一区二区三区四区激情视频 | 舔av片在线| 成人鲁丝片一二三区免费| 日韩,欧美,国产一区二区三区 | 成人二区视频| 99久久九九国产精品国产免费| 国产高清激情床上av| 插阴视频在线观看视频| 黄色欧美视频在线观看| 日本免费一区二区三区高清不卡| 日韩强制内射视频| 免费在线观看成人毛片| 联通29元200g的流量卡| 亚洲成人中文字幕在线播放| а√天堂www在线а√下载| 国产 一区 欧美 日韩| 国产精品福利在线免费观看| 18禁裸乳无遮挡免费网站照片| 赤兔流量卡办理| 国产爱豆传媒在线观看| 听说在线观看完整版免费高清| 男人的好看免费观看在线视频| 国产综合懂色| av在线观看视频网站免费| av在线播放精品| 欧美最黄视频在线播放免费| 午夜影院日韩av| 免费搜索国产男女视频| 级片在线观看| 成人毛片a级毛片在线播放| 深爱激情五月婷婷| 国产在视频线在精品| 亚洲精品亚洲一区二区| 亚洲精品影视一区二区三区av| 高清毛片免费看| 国产伦精品一区二区三区四那| 欧美绝顶高潮抽搐喷水| 色在线成人网| 天美传媒精品一区二区| 丝袜喷水一区| 日本黄色视频三级网站网址| 日韩欧美三级三区| 简卡轻食公司| 亚洲欧美日韩高清在线视频| 亚洲经典国产精华液单| 亚洲精品色激情综合| 一级av片app| 99久久久亚洲精品蜜臀av| 久久韩国三级中文字幕| 免费看光身美女| 国产成人福利小说| 亚洲一区高清亚洲精品| 欧美日韩综合久久久久久| 狂野欧美激情性xxxx在线观看| 蜜桃久久精品国产亚洲av| 亚洲精品一区av在线观看| 乱系列少妇在线播放| 两个人的视频大全免费| 最后的刺客免费高清国语| 国产伦精品一区二区三区视频9| 男女下面进入的视频免费午夜| 欧美一区二区国产精品久久精品| 久久久久精品国产欧美久久久| 亚洲真实伦在线观看| 精品久久国产蜜桃| 欧美成人免费av一区二区三区| 国产视频一区二区在线看| 欧美一区二区国产精品久久精品| 精品日产1卡2卡| 日韩欧美 国产精品| 女人十人毛片免费观看3o分钟| 中文字幕久久专区| 性欧美人与动物交配| 精品熟女少妇av免费看| 俄罗斯特黄特色一大片| 国产精品一区二区免费欧美| 18禁黄网站禁片免费观看直播| 少妇人妻精品综合一区二区 | 日本黄色片子视频| 熟女电影av网| 精品免费久久久久久久清纯| 最近最新中文字幕大全电影3| 国产精品久久久久久亚洲av鲁大| 欧美在线一区亚洲| 国产av麻豆久久久久久久| 国产一区二区在线观看日韩| 欧美色视频一区免费| 国产精品99久久久久久久久| 亚洲综合色惰| videossex国产| 一级av片app| av福利片在线观看| 我要看日韩黄色一级片| 午夜久久久久精精品| 久久久午夜欧美精品| 在线观看午夜福利视频| 日本熟妇午夜| 国产亚洲91精品色在线| 亚洲美女搞黄在线观看 | 我要看日韩黄色一级片| 99热精品在线国产| 国产精品免费一区二区三区在线| 如何舔出高潮| 99国产极品粉嫩在线观看| 久99久视频精品免费| 在线a可以看的网站| 黄色视频,在线免费观看| 真实男女啪啪啪动态图| 亚洲欧美日韩无卡精品| 简卡轻食公司| 国产精品久久久久久久电影| 九九久久精品国产亚洲av麻豆| 国产亚洲欧美98| 波多野结衣高清无吗| 久久久国产成人免费| 一级黄色大片毛片| 亚洲av第一区精品v没综合| 午夜视频国产福利| 久久久久久久久久黄片| 国产乱人偷精品视频| www日本黄色视频网| 国产精品一二三区在线看| 国产女主播在线喷水免费视频网站 | 亚洲欧美成人精品一区二区| 国产成人91sexporn| 国产69精品久久久久777片| 国产高潮美女av| 禁无遮挡网站| 日本成人三级电影网站| 日韩av在线大香蕉| 搞女人的毛片| 偷拍熟女少妇极品色| 日韩在线高清观看一区二区三区| 欧美精品国产亚洲| 色尼玛亚洲综合影院| 久久久午夜欧美精品| av女优亚洲男人天堂| 夜夜看夜夜爽夜夜摸| 香蕉av资源在线| 国产单亲对白刺激| 蜜臀久久99精品久久宅男| 国产麻豆成人av免费视频| 亚洲一区二区三区色噜噜| 午夜福利视频1000在线观看| 欧美不卡视频在线免费观看| 在线观看美女被高潮喷水网站| 亚洲国产日韩欧美精品在线观看| 午夜a级毛片| 免费看av在线观看网站| 91午夜精品亚洲一区二区三区| 亚洲av免费高清在线观看| 波野结衣二区三区在线| 久久韩国三级中文字幕| av在线亚洲专区| 淫秽高清视频在线观看| 亚洲精品乱码久久久v下载方式| 精品久久久久久久久亚洲| 久久草成人影院| 日本色播在线视频| 啦啦啦啦在线视频资源| 色在线成人网| 麻豆av噜噜一区二区三区| 欧美日本亚洲视频在线播放| 可以在线观看的亚洲视频| 91久久精品电影网| 一级毛片我不卡| 特级一级黄色大片| 一级毛片电影观看 | 日韩国内少妇激情av| 日本免费a在线| 伊人久久精品亚洲午夜| 国产一区亚洲一区在线观看| 老师上课跳d突然被开到最大视频| 日本成人三级电影网站| 精品免费久久久久久久清纯| 搞女人的毛片| 久久综合国产亚洲精品| 亚洲欧美成人精品一区二区| 国产高清视频在线观看网站| 国产在线精品亚洲第一网站| 91麻豆精品激情在线观看国产| 久久久久久久久久久丰满| 热99在线观看视频| 国产久久久一区二区三区| 色综合亚洲欧美另类图片| 国产久久久一区二区三区| 热99在线观看视频| 九九爱精品视频在线观看| 搡女人真爽免费视频火全软件 | 午夜日韩欧美国产| 一级毛片电影观看 | 成人毛片a级毛片在线播放| 国产淫片久久久久久久久| 三级国产精品欧美在线观看| 国产亚洲精品久久久com| 精品熟女少妇av免费看| 高清毛片免费观看视频网站| 人人妻人人澡人人爽人人夜夜 | 中文字幕av成人在线电影| 乱系列少妇在线播放| 国产爱豆传媒在线观看| 国产真实乱freesex| 国产高清视频在线播放一区| 最近在线观看免费完整版| 寂寞人妻少妇视频99o| 午夜精品在线福利| 三级国产精品欧美在线观看| 久久国内精品自在自线图片| 亚洲电影在线观看av| 男人和女人高潮做爰伦理| 色哟哟哟哟哟哟| 国产精品一区www在线观看| 国产精品综合久久久久久久免费| 在线免费十八禁| 又粗又爽又猛毛片免费看| 狂野欧美白嫩少妇大欣赏| 亚洲av中文av极速乱| 日本撒尿小便嘘嘘汇集6| 村上凉子中文字幕在线| 99热6这里只有精品| 亚洲激情五月婷婷啪啪|