• <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| 亚洲精品中文字幕在线视频| 亚洲人成网站高清观看| 男人的好看免费观看在线视频 | av国产免费在线观看| 高清在线国产一区| 国产欧美日韩精品亚洲av| 黑人欧美特级aaaaaa片| 免费在线观看成人毛片| 国产精品日韩av在线免费观看| 两个人的视频大全免费| 国产激情欧美一区二区| 成年人黄色毛片网站| 国产精品久久久人人做人人爽| 999久久久精品免费观看国产| 国产精品亚洲av一区麻豆| 久久久精品欧美日韩精品| 身体一侧抽搐| 男女那种视频在线观看| 在线观看一区二区三区| 免费在线观看日本一区| 母亲3免费完整高清在线观看| 欧美日本亚洲视频在线播放| 窝窝影院91人妻| 日韩av在线大香蕉| 久久国产乱子伦精品免费另类| 丁香六月欧美| 欧美极品一区二区三区四区| 我要搜黄色片| 国内揄拍国产精品人妻在线| 久久久久久久久中文| 国产成人一区二区三区免费视频网站| 久久久久久久午夜电影| 成人高潮视频无遮挡免费网站| 国产成人啪精品午夜网站| 欧美日韩乱码在线| 身体一侧抽搐| 亚洲av中文字字幕乱码综合| 一卡2卡三卡四卡精品乱码亚洲| 国产精品爽爽va在线观看网站| 国产成人欧美在线观看| 精品人妻1区二区| 看片在线看免费视频| 后天国语完整版免费观看| 老汉色av国产亚洲站长工具| 国产精品一区二区三区四区免费观看 | 人人妻人人看人人澡| x7x7x7水蜜桃| 日韩av在线大香蕉| 亚洲国产日韩欧美精品在线观看 | 哪里可以看免费的av片| 看片在线看免费视频| 亚洲人成电影免费在线| 中文字幕熟女人妻在线| 嫩草影视91久久| 久久午夜综合久久蜜桃| 99在线视频只有这里精品首页| 少妇裸体淫交视频免费看高清 | 美女黄网站色视频| 老司机靠b影院| 亚洲av成人一区二区三| 日本三级黄在线观看| 亚洲精品久久国产高清桃花| 国产三级黄色录像| ponron亚洲| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩有码中文字幕| 少妇裸体淫交视频免费看高清 | 日本 av在线| 麻豆成人av在线观看| 99精品在免费线老司机午夜| 黄片大片在线免费观看| 90打野战视频偷拍视频| 国产野战对白在线观看| 曰老女人黄片| 麻豆一二三区av精品| 国产区一区二久久| 国产麻豆成人av免费视频| a级毛片在线看网站| 叶爱在线成人免费视频播放| 国产精品九九99| 亚洲国产精品999在线| 中文亚洲av片在线观看爽| 成人欧美大片| 可以在线观看的亚洲视频| 狠狠狠狠99中文字幕| 国产av一区二区精品久久| 12—13女人毛片做爰片一| 成人精品一区二区免费| 精品欧美一区二区三区在线| 香蕉丝袜av| 久久久久久九九精品二区国产 | 久久久久久大精品| 色噜噜av男人的天堂激情| 亚洲av成人不卡在线观看播放网| 亚洲av片天天在线观看| 精品久久久久久久久久久久久| 最近在线观看免费完整版| 国产精品一区二区三区四区久久| 亚洲av成人不卡在线观看播放网| 国产黄片美女视频| 看黄色毛片网站| 亚洲第一电影网av| 一本久久中文字幕| 中国美女看黄片| 99国产综合亚洲精品| 人人妻人人澡欧美一区二区| 热99re8久久精品国产| 欧美又色又爽又黄视频| 在线十欧美十亚洲十日本专区| 亚洲男人天堂网一区| 久久中文看片网| 一二三四社区在线视频社区8| 亚洲国产欧美一区二区综合| 日本一区二区免费在线视频| 久久久久精品国产欧美久久久| 久久精品国产综合久久久| 不卡一级毛片| 国产精品影院久久| 久久香蕉国产精品| 香蕉国产在线看| 国产视频内射| 免费观看精品视频网站| 国产精品久久久久久久电影 | 免费高清视频大片| a在线观看视频网站| 亚洲av电影不卡..在线观看| 欧美绝顶高潮抽搐喷水| 久久久久久九九精品二区国产 | 精品午夜福利视频在线观看一区| 国产亚洲欧美98| 亚洲av片天天在线观看| 伊人久久大香线蕉亚洲五| 波多野结衣高清作品| 啪啪无遮挡十八禁网站| 正在播放国产对白刺激| 国产视频一区二区在线看| 天天躁狠狠躁夜夜躁狠狠躁| 性色av乱码一区二区三区2| 狂野欧美白嫩少妇大欣赏| 欧洲精品卡2卡3卡4卡5卡区| 久久国产精品影院| 2021天堂中文幕一二区在线观| 国产精品一区二区三区四区久久| www.www免费av| 99久久无色码亚洲精品果冻| 给我免费播放毛片高清在线观看| 亚洲国产中文字幕在线视频| 一级黄色大片毛片| av在线播放免费不卡| 免费看美女性在线毛片视频| 首页视频小说图片口味搜索| 午夜日韩欧美国产| 麻豆久久精品国产亚洲av| 国产精品久久久久久精品电影| 久久中文字幕一级| 欧美丝袜亚洲另类 | 日韩免费av在线播放| 在线永久观看黄色视频| 19禁男女啪啪无遮挡网站| 叶爱在线成人免费视频播放| 这个男人来自地球电影免费观看| 一本精品99久久精品77| 欧美 亚洲 国产 日韩一| 日韩成人在线观看一区二区三区| 免费看美女性在线毛片视频| 非洲黑人性xxxx精品又粗又长| 又黄又爽又免费观看的视频| 国产又色又爽无遮挡免费看| 亚洲精品色激情综合| av超薄肉色丝袜交足视频| 欧美另类亚洲清纯唯美| 男男h啪啪无遮挡| 每晚都被弄得嗷嗷叫到高潮| 日韩免费av在线播放| 五月玫瑰六月丁香| 丝袜人妻中文字幕| 黄色片一级片一级黄色片| 亚洲av电影不卡..在线观看| 在线观看66精品国产| 男女做爰动态图高潮gif福利片| 午夜a级毛片| 国产成人精品久久二区二区91| 又爽又黄无遮挡网站| 日本三级黄在线观看| 亚洲成人久久性| 国产av不卡久久| 后天国语完整版免费观看| 国产成人欧美在线观看| 成人永久免费在线观看视频| 日韩欧美 国产精品| 99久久久亚洲精品蜜臀av| 亚洲男人天堂网一区| 1024视频免费在线观看| 最新美女视频免费是黄的| 色哟哟哟哟哟哟| av有码第一页| 国产精品av视频在线免费观看| 一个人免费在线观看电影 | 免费高清视频大片| 国产精品久久视频播放| 丰满人妻熟妇乱又伦精品不卡| 在线a可以看的网站| 丰满的人妻完整版| www.www免费av| 亚洲avbb在线观看| 国产精品国产高清国产av| 99国产极品粉嫩在线观看| 1024香蕉在线观看| 丝袜人妻中文字幕| 黄色a级毛片大全视频| 精品一区二区三区视频在线观看免费| 在线免费观看的www视频| 久久精品91无色码中文字幕| 怎么达到女性高潮| 中文亚洲av片在线观看爽| 欧美午夜高清在线| 欧美av亚洲av综合av国产av| 欧美三级亚洲精品| 香蕉av资源在线| 黄片大片在线免费观看| 男女那种视频在线观看| 亚洲专区中文字幕在线| 国内毛片毛片毛片毛片毛片| 精品欧美国产一区二区三| 国产av一区在线观看免费| 99久久99久久久精品蜜桃| 日韩欧美精品v在线| а√天堂www在线а√下载| 国产区一区二久久| 搡老熟女国产l中国老女人| 亚洲自拍偷在线| 国内毛片毛片毛片毛片毛片| 97超级碰碰碰精品色视频在线观看| 日日摸夜夜添夜夜添小说| 国产精品综合久久久久久久免费| 听说在线观看完整版免费高清| 国产麻豆成人av免费视频| 亚洲中文日韩欧美视频| 国产野战对白在线观看| 12—13女人毛片做爰片一| 宅男免费午夜| 日本一区二区免费在线视频| 国产亚洲精品一区二区www| ponron亚洲| 国产精品自产拍在线观看55亚洲| 精品日产1卡2卡| 精品少妇一区二区三区视频日本电影| 宅男免费午夜| 最新在线观看一区二区三区| 欧美性长视频在线观看| 欧美日本亚洲视频在线播放| 亚洲国产欧美一区二区综合| 精品第一国产精品| ponron亚洲| 俄罗斯特黄特色一大片| 亚洲av中文字字幕乱码综合| 妹子高潮喷水视频| 国产黄a三级三级三级人| 国产成人aa在线观看| 国产免费男女视频| 国产av一区二区精品久久| 99在线视频只有这里精品首页| 亚洲av片天天在线观看| av欧美777| x7x7x7水蜜桃| 国产乱人伦免费视频| 日韩欧美一区二区三区在线观看| 麻豆国产av国片精品| 亚洲欧美日韩东京热| 亚洲中文av在线| 天堂影院成人在线观看| 日韩 欧美 亚洲 中文字幕| 日本一二三区视频观看| 免费无遮挡裸体视频| 在线免费观看的www视频| 国产成人一区二区三区免费视频网站| 在线观看66精品国产| 国产欧美日韩精品亚洲av| 国产成人aa在线观看| 国产精品一区二区三区四区免费观看 | avwww免费| 国产亚洲av嫩草精品影院| 欧美高清成人免费视频www| 精华霜和精华液先用哪个| 男女视频在线观看网站免费 | 91九色精品人成在线观看| 成人av一区二区三区在线看| 中文字幕熟女人妻在线| 国内久久婷婷六月综合欲色啪| 全区人妻精品视频| 99在线人妻在线中文字幕| 中出人妻视频一区二区| 亚洲自偷自拍图片 自拍| 少妇被粗大的猛进出69影院| 国产精品美女特级片免费视频播放器 | 亚洲国产精品sss在线观看| 90打野战视频偷拍视频| 成年女人毛片免费观看观看9| 麻豆成人av在线观看| 久久久久久国产a免费观看| 亚洲avbb在线观看| 在线观看www视频免费| 琪琪午夜伦伦电影理论片6080| 中出人妻视频一区二区| 免费在线观看影片大全网站| 欧美黑人欧美精品刺激| 97超级碰碰碰精品色视频在线观看| 国语自产精品视频在线第100页| 久久中文看片网| xxxwww97欧美| 国产亚洲精品久久久久久毛片| 黄色视频不卡| 免费在线观看黄色视频的| 97超级碰碰碰精品色视频在线观看| 欧美精品啪啪一区二区三区| 视频区欧美日本亚洲| 又爽又黄无遮挡网站| 人妻夜夜爽99麻豆av| 亚洲欧美一区二区三区黑人| 成人欧美大片| 91麻豆精品激情在线观看国产| 免费在线观看成人毛片| 欧美激情久久久久久爽电影| 长腿黑丝高跟| 99热这里只有精品一区 | 久久久精品欧美日韩精品| 少妇的丰满在线观看| 好看av亚洲va欧美ⅴa在| 99re在线观看精品视频| 毛片女人毛片| 婷婷亚洲欧美| 久久精品影院6| 人成视频在线观看免费观看| 亚洲专区字幕在线| 波多野结衣高清作品| 国产69精品久久久久777片 | 国产熟女xx| 欧美日本亚洲视频在线播放| 亚洲色图 男人天堂 中文字幕| 久久久久免费精品人妻一区二区| av天堂在线播放| 欧美性猛交黑人性爽| 校园春色视频在线观看| 欧美+亚洲+日韩+国产| cao死你这个sao货| 丁香六月欧美| 国产精品香港三级国产av潘金莲| 国产欧美日韩一区二区精品| 在线观看66精品国产| 亚洲av日韩精品久久久久久密| 床上黄色一级片| svipshipincom国产片| 日韩精品青青久久久久久| 久久中文字幕人妻熟女| 色综合欧美亚洲国产小说| 小说图片视频综合网站| 精品免费久久久久久久清纯| 国产av一区在线观看免费| 欧美成人性av电影在线观看| 日日夜夜操网爽| 2021天堂中文幕一二区在线观| 黄片大片在线免费观看| 欧美丝袜亚洲另类 | 草草在线视频免费看| 三级毛片av免费| 欧美日韩精品网址| 色综合站精品国产| 久久人妻av系列| 日韩成人在线观看一区二区三区| 午夜精品久久久久久毛片777| www日本在线高清视频| 美女午夜性视频免费| 中文字幕精品亚洲无线码一区| 国产视频一区二区在线看| 欧美av亚洲av综合av国产av| 国产成人av教育| 成人欧美大片| 午夜影院日韩av| 欧美人与性动交α欧美精品济南到| 国产精品一区二区三区四区免费观看 | 99久久国产精品久久久| 国产亚洲av高清不卡| 亚洲精品中文字幕一二三四区| 国产激情欧美一区二区| 欧美激情久久久久久爽电影| 亚洲最大成人中文| 久久久水蜜桃国产精品网| 黄色视频,在线免费观看| 精品久久久久久久久久久久久| 国产精品久久久久久亚洲av鲁大| 亚洲全国av大片| 久久久国产欧美日韩av| 免费搜索国产男女视频| 国产免费男女视频| 亚洲精品国产精品久久久不卡| 一本精品99久久精品77| 久久精品国产亚洲av香蕉五月| 国产69精品久久久久777片 | 色综合亚洲欧美另类图片| 女警被强在线播放| 精品一区二区三区视频在线观看免费| 老鸭窝网址在线观看| 香蕉久久夜色| 午夜视频精品福利| 91成年电影在线观看| 久久人妻av系列| 一区福利在线观看| 麻豆国产av国片精品| 少妇的丰满在线观看| 一区二区三区高清视频在线| www.www免费av| 日韩国内少妇激情av| 国产精品一区二区精品视频观看| 日韩成人在线观看一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 久久香蕉精品热| 亚洲熟妇熟女久久| 麻豆av在线久日| 国产午夜精品论理片| 亚洲精品色激情综合| 亚洲人成伊人成综合网2020| 怎么达到女性高潮| 一本精品99久久精品77| 性色av乱码一区二区三区2| 亚洲一区高清亚洲精品| 欧美 亚洲 国产 日韩一| 亚洲性夜色夜夜综合| 亚洲av片天天在线观看| 成人国产一区最新在线观看| 国产成人精品久久二区二区91| 亚洲av第一区精品v没综合| 夜夜爽天天搞| 观看免费一级毛片| 一区二区三区高清视频在线| 免费看美女性在线毛片视频| 色播亚洲综合网| 在线观看免费日韩欧美大片| xxx96com| 欧美又色又爽又黄视频| 色播亚洲综合网| 淫妇啪啪啪对白视频| 午夜两性在线视频| av国产免费在线观看| 亚洲人与动物交配视频| 老司机靠b影院| 丰满的人妻完整版| 老司机午夜福利在线观看视频| 日本免费一区二区三区高清不卡| 麻豆成人午夜福利视频| 亚洲色图av天堂| а√天堂www在线а√下载| 免费在线观看日本一区| 天堂动漫精品| 怎么达到女性高潮| 国产免费av片在线观看野外av| 国产激情久久老熟女| 美女大奶头视频| 99久久国产精品久久久| 国产爱豆传媒在线观看 | 欧美av亚洲av综合av国产av| 国内毛片毛片毛片毛片毛片| 久久精品国产综合久久久| 国产精品爽爽va在线观看网站| 亚洲在线自拍视频| 国产精品免费一区二区三区在线| 50天的宝宝边吃奶边哭怎么回事| 男女床上黄色一级片免费看| 桃红色精品国产亚洲av| 最新在线观看一区二区三区| 两个人的视频大全免费| 岛国在线免费视频观看| 中文字幕最新亚洲高清| 免费电影在线观看免费观看| 色尼玛亚洲综合影院| 免费在线观看影片大全网站| 男人的好看免费观看在线视频 | 国模一区二区三区四区视频 | 日本 欧美在线| 亚洲色图av天堂| 波多野结衣高清作品| 成人亚洲精品av一区二区| 天天一区二区日本电影三级| 国内毛片毛片毛片毛片毛片| 国产伦一二天堂av在线观看| 90打野战视频偷拍视频| 亚洲av成人精品一区久久| 国产精品久久久av美女十八| 青草久久国产| 搡老熟女国产l中国老女人| 国产精品免费一区二区三区在线| 好男人在线观看高清免费视频| 成人国产综合亚洲| 69av精品久久久久久| 欧美绝顶高潮抽搐喷水| 性色av乱码一区二区三区2| 欧美乱码精品一区二区三区| 国产一区二区三区在线臀色熟女| 中文字幕久久专区| 久久久水蜜桃国产精品网| 久久香蕉激情| 免费在线观看视频国产中文字幕亚洲| 亚洲欧美一区二区三区黑人| 99久久综合精品五月天人人| 黑人巨大精品欧美一区二区mp4| 亚洲精品久久国产高清桃花| 成人精品一区二区免费| а√天堂www在线а√下载| 曰老女人黄片| 国产伦人伦偷精品视频| 老熟妇乱子伦视频在线观看| 免费人成视频x8x8入口观看| 白带黄色成豆腐渣| 丰满的人妻完整版| 高清在线国产一区| 99热只有精品国产| 高潮久久久久久久久久久不卡| 亚洲18禁久久av| 精品久久久久久久久久免费视频| 欧美黄色片欧美黄色片| 中文字幕久久专区| 91成年电影在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲成人免费电影在线观看| 色播亚洲综合网| 成年版毛片免费区| 亚洲av熟女| 一级黄色大片毛片| 给我免费播放毛片高清在线观看| 岛国在线观看网站| 18禁国产床啪视频网站| 50天的宝宝边吃奶边哭怎么回事| 国内少妇人妻偷人精品xxx网站 | 波多野结衣高清无吗| 成年免费大片在线观看| 国产在线观看jvid| 国内久久婷婷六月综合欲色啪| 色综合亚洲欧美另类图片| 正在播放国产对白刺激| www日本黄色视频网| 亚洲人成伊人成综合网2020| 欧美久久黑人一区二区| 一进一出好大好爽视频| 欧美久久黑人一区二区| 丰满人妻一区二区三区视频av | 精品久久久久久,| 最近在线观看免费完整版| 中文亚洲av片在线观看爽| 久久亚洲真实| 欧美色欧美亚洲另类二区| 国产又色又爽无遮挡免费看| 99热这里只有精品一区 | 欧美日本亚洲视频在线播放| 最近视频中文字幕2019在线8| 亚洲最大成人中文| 亚洲aⅴ乱码一区二区在线播放 | 亚洲国产欧洲综合997久久,| 真人做人爱边吃奶动态| 99国产综合亚洲精品| 1024视频免费在线观看| 每晚都被弄得嗷嗷叫到高潮| 国内揄拍国产精品人妻在线| 亚洲欧美精品综合久久99| 真人做人爱边吃奶动态| 久久精品91无色码中文字幕| 精品一区二区三区视频在线观看免费| 别揉我奶头~嗯~啊~动态视频| 黄色a级毛片大全视频| 色综合站精品国产| 怎么达到女性高潮| 亚洲美女视频黄频| 久久人妻av系列| 亚洲中文av在线| 免费在线观看视频国产中文字幕亚洲| 色噜噜av男人的天堂激情| 99在线视频只有这里精品首页| 午夜免费成人在线视频| 久久久久久久午夜电影| 啦啦啦观看免费观看视频高清| 日韩 欧美 亚洲 中文字幕| 精品熟女少妇八av免费久了| 岛国在线观看网站| 国产私拍福利视频在线观看| 久久久久久久久免费视频了| 精品乱码久久久久久99久播| 日日爽夜夜爽网站| 岛国视频午夜一区免费看| 一级毛片精品| 男女之事视频高清在线观看| av福利片在线观看| 岛国在线免费视频观看| 在线观看www视频免费| 无人区码免费观看不卡| 亚洲专区中文字幕在线| 国产精品精品国产色婷婷| netflix在线观看网站| 亚洲一码二码三码区别大吗| 欧美乱码精品一区二区三区| 日韩欧美三级三区| 久久精品国产综合久久久| 一级片免费观看大全| 色哟哟哟哟哟哟| 99久久综合精品五月天人人| 亚洲国产精品sss在线观看| 国产av麻豆久久久久久久| 51午夜福利影视在线观看| 亚洲国产高清在线一区二区三| 日本一区二区免费在线视频| 中文字幕人成人乱码亚洲影| 久久人妻av系列| 亚洲男人天堂网一区|