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

    二甲醚摻混對(duì)氨/空氣預(yù)混燃燒特性的影響

    2025-08-15 00:00:00李天喜陳朝陽(yáng)康琦張達(dá)
    重慶大學(xué)學(xué)報(bào) 2025年7期
    關(guān)鍵詞:火焰自由基燃料

    中圖分類號(hào):TK411.2 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1000-582X(2025)07-085-09

    Effect of dimethyl ether addition on NH3. /airpremixed flames

    LI Tianxi,CHEN Zhaoyang,KANG Qi,ZHANG Da (School ofEnergy and Electrical Engineering,Chang'an University,Xi'an 71O064,P.R.China)

    Abstract: To improve the ignition and combustion performance of NH3 in engines,blendingitwith the highreactivity-fuel dimethyl ether (DME) is an effective strategy. This study conducts simulations of NH3/DME/air (20 premixed laminar flames to investigate the efects of DME addition on key combustion characteristics,including laminar flame speed,reaction pathways,and NO formation.The results show that both the adiabatic flame temperature and the laminar flame speed increase significantly with higher DME blending ratios.A strong correlation is observed between laminar flame speed and the concentration of reactive free radicals,suggesting that the increased radical concentration due to DMEisthe main contributor to the enhanced flame speed. Furthermore,the normalized NO concentration in the flame rises significantly with increasing DME content, reaching approximately 50% when the DME blending ratio reaches 80% . Reaction pathway analysis indicates that DME addition inhibits the conversion of nitrogen species to N2 ,thereby leading to increased NO emissions. Sensitivity analysis shows that DME significantly alters the dominant elementary reactions.As the DME ratio increases,thesensitivity of nitrogen-group reactions considerably declines,while hydrogen-and carbon-group reactions increasingly govern the combustion process of the blended fuel.

    Keywords: ammonia; reaction kinetics; dimethyl ether blending; laminar flame speed; reaction pathways

    隨著化石能源的開采和利用,溫室氣體排放引起的全球氣候問(wèn)題成了世界性的重要議題。中國(guó)明確提出 CO2 排放力爭(zhēng)于2030年前達(dá)到峰值,爭(zhēng)取2060年前實(shí)現(xiàn)碳中和。減少交通領(lǐng)域碳排放能有效應(yīng)對(duì)氣候變化、助力實(shí)現(xiàn)碳達(dá)峰碳中和目標(biāo)有重要意義1。氨 ?NH 3 )是一種無(wú)碳儲(chǔ)氫燃料2,含氫質(zhì)量分?jǐn)?shù)為 17.6% ,易于液化,便于儲(chǔ)存和運(yùn)輸,被認(rèn)為是一種有前途的碳中和燃料[-4]。與氫能相比,氨燃料還具有單位儲(chǔ)能成本低、體積能量密度高、安全可靠等優(yōu)點(diǎn)5;但氨反應(yīng)活性低,點(diǎn)火溫度高且燃燒速度慢,直接燃燒利用存在一定困難,通常采用與其他活性物質(zhì)摻燒的方法來(lái)改善其燃燒性能[8--]。二甲醚(DME,分子式 CH3OCH3 是一種低碳含氧燃料,十六烷值高,著火性能好且來(lái)源廣泛,是良好的潔凈替代燃料和燃料活性添加劑[12]。近年來(lái),有學(xué)者將氨/二甲醚混合燃料應(yīng)用在壓燃式發(fā)動(dòng)機(jī)上,探索氨作為發(fā)動(dòng)機(jī)替代燃料的可行性[3-5]。研究發(fā)現(xiàn),二甲醚摻混能夠改善氨燃料的著火燃燒性能,拓寬氨燃料在發(fā)動(dòng)機(jī)上的運(yùn)行范圍,但氮氧化物 (NOx) 的排放較傳統(tǒng)燃料有所升高。

    層流火焰速度是混合氣體的基本特性參數(shù)之一,是構(gòu)建和驗(yàn)證燃料化學(xué)反應(yīng)機(jī)理的重要依據(jù)[。近年來(lái),學(xué)者對(duì)氨及其混合燃料的層流燃燒特性進(jìn)行了研究。Yin等[1在定容燃燒彈中測(cè)量了高溫高壓下 NH3/ DME混合燃料層流燃燒速度,發(fā)現(xiàn)氨基自由基( NH2 )的氧化主要以含碳物質(zhì)為主導(dǎo), 50% 的 NH2 通過(guò) NH2+ CH2O=NH3+CHO 和 NH2+CH3=CH3NH2 反應(yīng)被消耗;Dai等[8]利用快速壓縮機(jī)裝置測(cè)量了 NH3/DME 混合燃料的著火延遲時(shí)間,發(fā)現(xiàn)摻混 2%~5%DME 的混合燃料,在當(dāng)量比 ?=1.0 和 ?=2.0 時(shí),均出現(xiàn)了兩階段著火現(xiàn)象,DME的低溫氧化放熱在 NH3/DME 混合燃料自著火過(guò)程中有著重要的作用。還有研究致力于 NH3/DME 混合燃料化學(xué)反應(yīng)機(jī)理的構(gòu)建與優(yōu)化[19-22]。

    全面提升氨燃料的燃燒與排放性能,了解二甲醚添加對(duì)氨燃料燃燒的影響規(guī)律及作用機(jī)理,對(duì)寬廣初始條件下的 NH3/DME 混合燃料開展系統(tǒng)的基礎(chǔ)燃燒研究非常必要。文中利用化學(xué)反應(yīng)動(dòng)力學(xué)模擬方法,對(duì)初始溫度為 298K ,初始?jí)毫?0.1MPa 時(shí),寬摻混比( 0~100% )和當(dāng)量比(0.7~1.5)條件下的NH/DME/空氣預(yù)混合氣的層流燃燒特性進(jìn)行研究,系統(tǒng)分析了二甲醚摻混對(duì)氨燃料燃燒過(guò)程的影響規(guī)律和動(dòng)力學(xué)機(jī)制,重點(diǎn)考慮二甲醚摻混對(duì)氨燃料燃燒反應(yīng)路徑及NO生成規(guī)律的影響。

    1化學(xué)反應(yīng)機(jī)理的選擇與驗(yàn)證

    近年來(lái),不少學(xué)者對(duì) NH3/DME 混合燃料的化學(xué)反應(yīng)機(jī)理開展了研究。Glarborg等[23]構(gòu)建了考慮碳族與氮族元素交叉反應(yīng)的 NH3/DME 化學(xué)反應(yīng)機(jī)理,優(yōu)化了NO相關(guān)反應(yīng)的速率常數(shù)與熱化學(xué)參數(shù),該機(jī)理經(jīng)過(guò)著火延遲期和物質(zhì)濃度的驗(yàn)證,但在常溫常壓條件下對(duì)層流燃燒速度的預(yù)測(cè)能力欠佳。Meng等[]在Glarborg研究的基礎(chǔ)上,結(jié)合Burke等[24開發(fā)的二甲醚低溫氧化機(jī)理,同時(shí)耦合Dai等[18]提出的機(jī)理中 NH3 和DME的相關(guān)反應(yīng),構(gòu)建了適用于寬溫度范圍的NH3/DME混合燃料機(jī)理(221組分,1597反應(yīng))。Xiao等[2提出了詳細(xì)的 NH3/DME 混合機(jī)理(102組分,594反應(yīng)),并基于層流燃燒速度和著火延遲期實(shí)驗(yàn)結(jié)果對(duì)機(jī)理的可靠性進(jìn)行了驗(yàn)證。Issayev等[21]在Shrestha等[25]所提出機(jī)理的基礎(chǔ)上,對(duì) NH3/DME 混合燃料化學(xué)反應(yīng)機(jī)理進(jìn)行了改進(jìn),并對(duì)改進(jìn)后的機(jī)理開展了廣泛的實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證。

    文中采用化學(xué)反應(yīng)動(dòng)力學(xué)模擬軟件ChemkinI中的PREMIX模塊,對(duì)初始溫度 T0=298K 、初始?jí)毫?p0= 0.1MPa 時(shí)的NH/DME/空氣混合氣的層流燃燒速度進(jìn)行了模擬計(jì)算。將基于Meng機(jī)理[1]、Issayev機(jī)理21和Xiao機(jī)理[2]的模擬預(yù)測(cè)結(jié)果與Xiao等[20]的實(shí)驗(yàn)值進(jìn)行了對(duì)比,如圖1所示。對(duì)比后發(fā)現(xiàn),Issayev機(jī)理的預(yù)測(cè)結(jié)果高于其他2種機(jī)理的預(yù)測(cè)值和實(shí)驗(yàn)值,且在二甲醚摻混比增大時(shí)誤差增大;Xiao機(jī)理的預(yù)測(cè)結(jié)果與實(shí)驗(yàn)值一致性較好,但過(guò)度預(yù)測(cè)了貧燃料側(cè)混合物的層流燃燒速度。Meng機(jī)理的預(yù)測(cè)結(jié)果在多個(gè)摻混比下均能與實(shí)驗(yàn)值高度吻合,表現(xiàn)出了更好的預(yù)測(cè)能力。因此,文中選用Meng機(jī)理對(duì) NH3/DME /空氣預(yù)混層流燃燒特性進(jìn)行動(dòng)力學(xué)分析。為了方便表述,定義DME摻混比 XDME 為二甲醚在混合燃料中的體積分?jǐn)?shù)。

    2 結(jié)果與分析

    2.1層流燃燒速度和絕熱火焰溫度

    層流燃燒速度是混合氣體的重要宏觀燃燒特性參數(shù)。圖1所示為 TL0=298KL,pL0=0.1MPa 時(shí),不同二甲醚摻混比下,NH/DME/空氣混合氣的層流燃燒速度模擬值隨當(dāng)量比 ? 的變化關(guān)系。由圖可知,純 NH3 燃燒時(shí),層流燃燒速度較小,隨著二甲醚摻混比的增大,混合燃料的層流燃燒速度明顯增大,且在 ?=1.1 附近達(dá)到最大。圖2給出了不同當(dāng)量比下,混合燃料層流燃燒速度隨摻混比變化的擬合關(guān)系,可以看出,3個(gè)當(dāng)量比條件下,混合燃料的層流燃燒速度隨DME摻混比的增加呈線性增加的趨勢(shì)。在當(dāng)量比為0.8時(shí)的稀混合氣條件下,擬合直線的斜率較小;隨著當(dāng)量比逐漸增大,擬合線斜率也隨之增大,說(shuō)明二甲醚摻混對(duì)較濃混合氣燃燒過(guò)程的促進(jìn)效果更加明顯。

    圖1不同DME摻混比下層流燃燒速度預(yù)測(cè)結(jié)果與實(shí)驗(yàn)值對(duì)比 (T0=298K,p0=0.1MPa Fig.1Comparison of flame speed by simulations and experiments at different DME mixing ratios 0 T0=298K , p0=0.1 MPa)

    圖2不同當(dāng)量比條件下層流燃燒速度隨摻混比的變化規(guī)律 (T?0=298K,p?0=0.1MPa ) Fig.2Variation of flame speed with DME mixing ratio at different equivalence ratiost (T?0=298K,p?0=0.1 MPa)

    燃燒速度受燃燒溫度的影響很大。圖3所示為不同DME摻混比時(shí), NH3/DME /空氣混合氣絕熱火焰溫度隨當(dāng)量比的變化關(guān)系。可以看出,絕熱火焰溫度隨當(dāng)量比和摻混比的變化規(guī)律與層流燃燒速度一致,純氨燃燒時(shí)絕熱火焰溫度較低,對(duì)應(yīng)于較小的層流燃燒速度;隨著DME摻混比的增大,絕熱火焰溫度逐漸升高,層流燃燒速度也隨之增大,表明在較濃混合氣條件下,DME摻混對(duì)絕熱火焰溫度及層流燃燒速度的影響更為顯著;另一方面,溫度的升高也可能引起NO生成量的增加,火焰面中NO的濃度隨DME摻混比的變化規(guī)律將在后續(xù)內(nèi)容中重點(diǎn)討論。

    圖3不同DME摻混比下絕熱火焰溫度模擬值隨當(dāng)量比的變化規(guī)律( T0=298 K, p0=0.1 MPa) Fig.3Variation of predicted adiabatic flame temperature with equivalence ratioat diferentDME mixing ratios MPa)

    2.2 主要自由基濃度分布

    鏈?zhǔn)椒磻?yīng)理論認(rèn)為,燃料燃燒過(guò)程和層流燃燒速度強(qiáng)烈依賴于混合氣中活性自由基的濃度。 H,O,OH 自由基是最為活躍的活性自由基,對(duì)火焰發(fā)展的影響也最大。圖4為模擬得到的 T?0=298K,p?0=0.1MPa,?= 0.8時(shí), NH3/DME /空氣火焰中 自由基摩爾分?jǐn)?shù)隨火焰距離的變化關(guān)系??梢钥闯觯?種摻混比條件下,火焰中OH自由基的摩爾分?jǐn)?shù)均最大,H自由基次之,O濃度最小,隨著DME摻混比的增加,3種自由基的摩爾分?jǐn)?shù)均呈現(xiàn)增大趨勢(shì),導(dǎo)致氨的層流燃燒速度加快。圖5給出了火焰面上自由基濃度峰值隨DME摻混比的變化關(guān)系,可以看出隨著摻混比的增大,3種自由基濃度峰值均逐漸增大,3種自由基摩爾分?jǐn)?shù)峰值之和隨DME摻混比的變化趨勢(shì)與層流燃燒速度的變化趨勢(shì)相同,二甲醚摻混引起了混合燃料火焰中自由基濃度增大,從而導(dǎo)致層流燃燒速度增大,進(jìn)一步說(shuō)明了火焰中自由基濃度對(duì) NH3/DME /空氣混合氣層流燃燒速度及燃燒過(guò)程的影響顯著。

    圖4不同DME摻混比下H、O、OH自由基摩爾分?jǐn)?shù)模擬值隨火焰距離的變化規(guī)律 Fig.4Predictions of H, O and OH mole fractionat differentDME mixing ratios

    圖5層流燃燒速度與自由基摩爾分?jǐn)?shù)峰值隨DME摻混比的變化規(guī)律 Fig.5Variation of flame speed and peak values of radical mole fraction with DME mixing ratios (2號(hào)

    2.3 氮氧化物濃度分布

    含氮燃料燃燒過(guò)程中,不可避免地會(huì)產(chǎn)生 NOx 排放,探究不同摻混比下 NOx 的濃度分布具有重要意義。圖6給出了模擬得到的 T?0=298K.p?0=0.1MPa.?=0.8 時(shí),不同DME摻混比下 NH3,NO 和 NO2 的摩爾分?jǐn)?shù)隨火焰距離的變化關(guān)系。由圖可知,在 NH3/DME 火焰中,相比于NO, NO2 濃度很小;隨火焰距離的增大, NO2 的摩爾分?jǐn)?shù)呈現(xiàn)先增大后減小的趨勢(shì),NO的摩爾分?jǐn)?shù)從一定火焰位置開始逐漸增大,在 NH3 濃度減小為0的火焰位置時(shí),達(dá)到最大值后保持穩(wěn)定,表明 NO2 為燃燒過(guò)程的中間物種且濃度較小, NH3/DME 燃燒生成的氮氧化物主要是NO。隨著二甲醚摻混比的增大,NO的摩爾分?jǐn)?shù)先增大,但當(dāng)摻混比由 50% 增加到 80% 時(shí),NO的摩爾分?jǐn)?shù)又減小,這主要是由于 NH3 在混合燃料中占比下降,導(dǎo)致氮原子濃度減小。

    注:圖中NO和NO2曲線分別為放大5倍和100的結(jié)果

    圖6不同DME摻混比下 NH3 、NO和 NO2 摩爾分?jǐn)?shù)隨火焰距離的變化規(guī)律 (T?0=298K,p?0=0.1MPa,?=0.8) Fig.6Variation of mole fraction of NHz ,NO and NO2 with flame distanceunder different DME mixing ratios

    為了排除不同摻混比下氨燃料總含量的不同對(duì) NOx 濃度的影響,定義NO、 NO2 的歸一化濃度為NO、 NO2 的摩爾分?jǐn)?shù)與 NH3 的初始摩爾分?jǐn)?shù)之比。不同二甲醚摻混比下,NO和 NO2 的歸一化濃度隨火焰距離的變化如圖7所示,可以看出,隨著二甲醚摻混比的增大,NO和 NO2 的歸一化濃度增加,且摻混比越大,增加幅度越大。當(dāng)DME摻混比為 80% 時(shí),NO的歸一化濃度達(dá)到了 50% ,這是因?yàn)殡S著二甲醚的加入,自由基池中H、OH、O自由基增多,加速了含氮基團(tuán)氧化為NO的進(jìn)程,導(dǎo)致NO濃度升高。

    注:圖中NO2曲線為放大20倍的結(jié)果

    圖7不同DME摻混比下NO、 NO2 歸一化濃度隨火焰距離的變化規(guī)律 (T?0=298K,p?0=0.1 MPa, ?=0.8 門Fig.7VariationofmormalizedconcentrationofNOand NO2 withflamedistanceunderdifferentDMEmixingratios(T?0=298K,p?0=0.1 MPa, ?=0.8 )

    2.4 反應(yīng)路徑分析

    為了明確二甲醚摻混對(duì)氨燃料燃燒過(guò)程及NO生成的影響規(guī)律,文中對(duì)不同摻混比預(yù)混火焰中燃料的化學(xué)反應(yīng)路徑進(jìn)行了對(duì)比分析。圖8給出了 T?0=298K?,p0=0.1MPa,?=0.8 時(shí),不同DME摻混比 (0%.20% 、50%.80%.10% 條件下的 NH3/DME/ 空氣預(yù)混火焰中 NH3 和DME的主要消耗路徑。可以看出,二甲醚主要通過(guò)脫氫與直接裂解反應(yīng)消耗,其中脫氫反應(yīng)路線占主導(dǎo)地位。

    脫氫: CH3OCH3+(H/OH/O)?CH3OCH2+(H2/H2O/OH)

    裂解: CH3OCH3?CH3O+CH3

    二甲醚脫氫產(chǎn)物 CH3OCH2 經(jīng) β 裂解生成甲醛( CH2O )和甲基, CH2O 經(jīng)過(guò)2次脫氫生成CO,CO最終被OH氧化為 CO2

    混合燃料中二甲醚的反應(yīng)路徑與純二甲醚燃燒時(shí)基本相同,摻混比的不同僅會(huì)導(dǎo)致其起始脫氫與裂解反應(yīng)所占比重發(fā)生變化。摻混燃燒時(shí),受氨的脫氫反應(yīng)競(jìng)爭(zhēng),DME脫氫反應(yīng)比重下降,裂解反應(yīng)比重增加;純DME燃燒時(shí),其脫氫和裂解反應(yīng)占比分別為 97.4% 和 2.6% ,添加 20%DME 混合氣時(shí),二甲醚脫氫反應(yīng)占比下降為 85.4% ,裂解反應(yīng)占比上升至 14.6% 。二甲醚脫氫產(chǎn)物的后續(xù)反應(yīng)路徑幾乎不受燃料摻混比變化的影響,但由于氨燃料的引入增加了含氮基團(tuán)參與的反應(yīng)(如HNO和 NO2 與 CH3 的反應(yīng),NO與HCO的反應(yīng)),導(dǎo)致 NO2 與 CH3 反應(yīng)生成NO的占比在DME摻混比為 20%~80% 時(shí)都超過(guò)了 10% 。由此可知,碳族和氮族的交叉反應(yīng)對(duì)整個(gè)燃燒過(guò)程以及NO的排放具有重要的影響。

    氨的消耗起始于燃料分子與活性自由基 (H,OH,O) 的脫氫反應(yīng)。隨后,其脫氫產(chǎn)物 NH2 通過(guò)相繼脫氫(路徑1)、氧化(路徑2)或與NO反應(yīng)(路徑3)消耗,最終生成 N2 。氨的氧化過(guò)程可以為3個(gè)主要通道,其中路徑1和2為主要反應(yīng)通道,2條路徑在不同摻混比下的占比均在 30% 以上。

    路徑 2:NH3?NH2?(H2NO?HNO)/HNO?NO?(N2O?N2)/N2

    路徑 3:NH3?NH2?N2

    由自由基摩爾分?jǐn)?shù)的分析可知,隨著二甲醚摻混比的增大,火焰中自由基濃度顯著增大,這會(huì)促進(jìn)氨的起始脫氫反應(yīng)發(fā)生,提高氨的消耗反應(yīng)速率,使層流火焰速度增大。DME摻混后,路徑1中 NH3 的次級(jí)脫氫反應(yīng)( NH2?NH 占比有所下降,但摻混比進(jìn)一步增加時(shí)變化不大;次級(jí)脫氫產(chǎn)物NH的后續(xù)反應(yīng)占比隨DME摻混變化較大,其中生成NO的反應(yīng)鏈 (NHNO,NHNNO) 占比增大,生成 N2 的反應(yīng)鏈( N2?NH?N2O?N2) 占比在摻混 20%DME 時(shí)有所增大,而當(dāng)摻混比進(jìn)一步增加時(shí),又逐漸減小,表明 20%

    DME摻混,有利于 NH3 經(jīng)路徑1通道向 N2 轉(zhuǎn)化,進(jìn)一步增大摻混比則可能引起NO生成量的增加;路徑2中的關(guān)鍵反應(yīng)步驟 NH2?HNO 和 NH2?H2NO 的占比均隨DME摻混比的增大而增大,分別從純氨時(shí)的 21.5% 和0提升至DME摻混比為 80% 時(shí)的 41.0% 和 19.0% ,該路徑的主要中間物種HNO幾乎全部轉(zhuǎn)化為NO,NO進(jìn)一步與含氮基團(tuán)反應(yīng)最終轉(zhuǎn)化為 N2 或被 HO2 氧化生成 NO2 。當(dāng)DME摻混比為 20% 時(shí),NO的各反應(yīng)通道占比與純氨燃燒時(shí)差別不大,但當(dāng)DME摻混比進(jìn)一步增加時(shí),NO向 N2 轉(zhuǎn)化的反應(yīng)鏈 (NON2,NON2ON2) 占比顯著下降,而NO氧化生成 NO2 的反應(yīng)占比卻顯著增加,原因在于隨著DME摻混比的增加,自由基池中含氮基團(tuán)的比例下降, HO2 自由基濃度上升,表明DME摻混比大于 20% 后,混合燃料的NO排放可能顯著增加,必須考慮一定的凈化措施;路徑3在氨的反應(yīng)過(guò)程中占比較小,隨著DME摻混比由0增大到 80% ,路徑3中 NH2?N2 的反應(yīng)占比由 20.9% 降低至 6.6% 。

    NH3 CH3OCH3 10.2%12.6% +100 21.5% +

    ++ H2NO 04000 NH2 +0 30.2% 30% 0%,3.1%4.2%2.6% CH30 (204 33.1%+0H 99.4%+M,+02 19.2% 9.6%+0H 0% 33.7% 20.3% 7.4%5.4% 8.0% +HO259.6% 29.7% CH20 (2 CH3 - (2號(hào) (204 C2H6 20.9% NH NO 22.7%

    +NO 18.80% 12.1% 52 NH 96.6% 8 +HM 年 +H21.7%+0H42.5%51.6% 39.8% 100% +OH28.3%+0240.6%45.7% 32.2% 100% HCO CH4 N2 +N05806 1399% 553% 248%12% 1 +M2 ↑ 38.6% +NH2 100% +H +H95.5%,98.5%,99.0%,99.2% +OH 100% +N,+NH2 45.4%,46.2%,38.2%,7.4% NO2 CO +OH (204號(hào) CO2 (20

    2.5 敏感性分析

    圖9所示為 T?0=298K.p?0=0.1MPa.?=0.8 時(shí),不同DME摻混比 (0%,20%,50%,80%) 條件下混合氣的質(zhì)量流量敏感性系數(shù)。敏感性系數(shù)為正值表示,當(dāng)該基元反應(yīng)速率增大時(shí),火焰?zhèn)鞑ニ俣仍龃螅砻髟摲磻?yīng)對(duì)混合燃料的燃燒起促進(jìn)作用;相反,敏感性系數(shù)為負(fù)值表明起抑制作用。對(duì)混合燃料的燃燒過(guò)程起促進(jìn)作用的反應(yīng)主要有

    R1: H+O2=OH+O R :27:CO+OH=CO2+H R30:HCO+M=H+CO+M

    其中,R1反應(yīng)在所有摻混比混合氣燃燒過(guò)程中均具有最高的敏感度,且敏感性系數(shù)隨DME摻混比的增大而降低;R27和R30反應(yīng)的敏感性系數(shù)隨DME摻混比的增大而增大,這是因?yàn)殡S著DME摻混比的增大,火焰中含碳組分濃度增加,含碳組分的氧化反應(yīng)對(duì)整個(gè)燃燒過(guò)程的影響越來(lái)越顯著。另外,DME摻混引起火焰中總自由基濃度增加,使得鏈分支反應(yīng)R1的敏感性系數(shù)減小,且摻混 20% DME時(shí),降低幅度最大。對(duì)混合燃料的燃燒過(guò)程起抑制作用的基元反應(yīng)主要有:

    R9: H+O2(+M)=HO2+(M)

    (2號(hào) R31:HCO+O2=CO+HO2

    R9和R31反應(yīng)消耗了 H 和HCO,生成了相對(duì)穩(wěn)定的 HO2 ,屬于鏈終止反應(yīng),抑制了反應(yīng)的進(jìn)行。當(dāng)DME

    摻混比為 20% 時(shí),含氮物質(zhì)的相關(guān)基元反應(yīng)R737和R803的敏感性系數(shù)顯著減小。值得注意的是,在DME摻混比為 80% 時(shí),基元反應(yīng)R803對(duì)燃燒的促進(jìn)作用轉(zhuǎn)變?yōu)橐种谱饔?,說(shuō)明含氮基團(tuán)對(duì)反應(yīng)過(guò)程的影響相對(duì)減弱,表明了碳族基團(tuán)對(duì)活性自由基的爭(zhēng)奪能力強(qiáng)于氮族,隨著DME摻混比的增大,燃燒反應(yīng)越來(lái)越受碳族反應(yīng)主導(dǎo)。這一發(fā)現(xiàn)與Issayev等[2的觀點(diǎn)相一致,即隨著DME摻混比的增大,大多數(shù)氮族基元反應(yīng)的敏感性均存在降低的趨勢(shì)。

    3結(jié)論

    文中使用化學(xué)反應(yīng)動(dòng)力學(xué)模擬軟件ChemkinⅡI,對(duì)初始溫度為 298K 、初始?jí)毫?0.1MPa 、當(dāng)量比為0.7~1.5 條件下的 NH3/DME/ 空氣混合層流火焰進(jìn)行了數(shù)值模擬,分析了二甲醚摻混對(duì)混合物層流燃燒速度、火焰結(jié)構(gòu)、NO生成路徑等的影響規(guī)律,主要結(jié)論如下。

    1)混合物的層流燃燒速度隨著二甲醚摻混比的增大而增大,且在較濃混合氣時(shí)增大的幅度更大。

    2)隨著二甲醚摻混比的增大,火焰中 H,O,OH 自由基摩爾分?jǐn)?shù)增大,混合物層流燃燒速度與摩爾分?jǐn)?shù)濃度密切相關(guān),二甲醚摻混引起的自由基池濃度的增大是層流燃燒速度增大的主要原因。

    3)隨著二甲醚摻混比的增大,火焰中NO的歸一化濃度顯著增大,當(dāng)DME摻混比為 80% 時(shí),NO的歸一化濃度達(dá)到了 50% 。路徑分析也表明,高比例二甲醚摻混不利于氨燃料燃燒中氮元素向 N2 的轉(zhuǎn)化,會(huì)導(dǎo)致NO排放增加。

    4)二甲醚摻混顯著影響了主要基元反應(yīng)的敏感性系數(shù),當(dāng)DME摻混比為 20% 時(shí),氮族基元反應(yīng)的敏感性明顯降低,混合物燃燒過(guò)程中碳族基元主導(dǎo)著燃燒過(guò)程的進(jìn)行。

    參考文獻(xiàn)

    [1]馬冬,竇廣玉,彭頔.我國(guó)機(jī)動(dòng)車碳減排潛力分析及建議[J].世界環(huán)境,2023(3):56-59. Ma D,Dou G Y,Peng D.Analysis ofcarbonemission reduction potentialof motorvehicles in China and therelated recommendations[J].World Environment,2023(3):56-59.(in Chinese)

    [2]AlturaifiSA,MathieuO,PetersenEL.Anexperimentaland modelingstudyofammoniapyrolysis[J].CombustionandFlame, 2022,235:111694.

    [3]Lesmana H,Zhang Z Z,LiXM,et al. NH3 as a transport fuel in internal combustion engines: a technical review[J]. Journal of Energy Resources Technology,2019,141(7): 070703.

    [4]CardosoJS,Silva V,RochaRC,et al.Ammonia asan energyvector:currentand future prospectsforlow-carbon fuel applications in internal combustion engines[J].Journal ofCleaner Production,2021,296:126562.

    [5]KangLWPanWG,ZhangJK,etal.Aeviewonammoniabendscombustionforidustrialappcations[J].Fuel,2,2: 126150.

    [6]高虎,劉凡,李海.碳中和目標(biāo)下氨燃料的機(jī)遇、挑戰(zhàn)及應(yīng)用前景[J].發(fā)電技術(shù),2022,43(3):462-467. Gao H,LiuF,LiH.Opportunities,challengesandapplicationprospectsofammonia fuelunder thetargetofcarbonneutrality [J].Power Generation Technology,2022, 43(3): 462-467. (in Chinese)

    [7]Lhuiller C,BrequignyP,Lamoureux N,etal.Experimentalinvestigationonlaminarburingvelocitisofammonia/hdrogen/ air mixtures at elevated temperatures[J].Fuel, 2020,263:116653.

    [8]許煥煥,葛一,李強(qiáng),等.氨燃料及應(yīng)用技術(shù)研究進(jìn)展[J].東北電力大學(xué)學(xué)報(bào),2022,42(2):1-13. XuHH,GeY,LiQ,etal.Researchprogressofammoniafuelandapplicationtechnology[J].JouralofNortheastElectric PowerUniversity,2022,42(2):1-13. (in Chinese)

    [9]Okafor E C,Naito Y, Colson S,et al. Experimental and numerical study of the laminar burning velocity of CH4 (2 ?NH3 -air premixed flames[J]. Combustion and Flame,2018,187: 185-198.

    [10] Jin B Z,Deng YF,Li G X,etal. Experimental and numerical studyof the laminar burning velocity of NH3/H2/ air premixed flames at elevated pressure and temperature[J]. InternationalJournalofHydrogen Energy,2022,47(85): 36046-36057.

    [11]陳朝陽(yáng),耿莉敏,鞏靜,等.摻氫對(duì)二甲醚預(yù)混層流燃燒特性的影響[J].西安交通大學(xué)學(xué)報(bào),2014,48(6):122-126. Chen ZY,GengLM,Gong J,etal.Influenceof hydrogenaditiononcombustioncharacteristicsof premixedlaminar dimethyl ether flames[J]. Journal of Xi'an Jiaotong University,2014,48(6): 122-126.(in Chinese)

    [12]李東昌,汪映,孟慶斌,等.二甲醚預(yù)混比對(duì)預(yù)混壓燃發(fā)動(dòng)機(jī)性能影響的試驗(yàn)研究[J].西安交通大學(xué)學(xué)報(bào),2013,47(3): 48-52. Li D C, Wang Y,Meng QB,etal.Experimental studyon effects of DME premixed ratioon performance of a PCCI engine[J]. Journal ofXi'an Jiaotong University,2013,47(3): 48-52.(in Chinese)

    [13]GrossCW,Kong SC.Performancecharacteristcsofacompresio-ignitionengineusingdirect-jectionammoia-DME mixtures[J]. Fuel,2013,103: 1069-1079.

    [14]Ryu K H,Zacharakis-Jutz G,Kong SC.Efectsoffuelcompositionson dieselengine performanceusing ammonia-DME mixtures[C]//SAE Technical Paper Series.SAE International,2013:1-6.

    [15]丁夢(mèng)豪,胡以懷,高雨顧,等.氨燃料在發(fā)動(dòng)機(jī)中的應(yīng)用研究[J].能源與節(jié)能,2021(12):163-165. Ding M H,HuYH,GaoYQ,etal.Aplicationof ammonia fuel inengines[J].EnergyandEnergyConservation,2021(12): 163-165.(in Chinese)

    [16]胡二江,黃佐華,姜雪,等.C1-C4烷烴預(yù)混層流燃燒與著火特性研究[J].工程熱物理學(xué)報(bào),2013,34(3):558-562. HuEJ,HuangZH,JiangX,etal.KineticstudyonlaminarburningvelocitiesandignitiondelaytimesofC1-C4alkanes[J]. Journal of Engineering Thermophysics,2013,34(3): 558-562.(in Chinese)

    [17]YinGY,LiJL,ZhouM,etal.Experimentalandkineticstudyonlminarflamespeedsofmmonia/dimethylethrairder high temperature and elevated pressure[J]. Combustion and Flame,2022,238:111915.

    [18]DaiLM,HashemiH,GlarborgP,etal.IgnitiondelaytimesofNH/DMEblendsathigh pressreandlowDMEfraction:RCM experiments and simulations[J]. Combustion and Flame,2021,227: 120-134.

    [19]Meng X Y,Zhang M K,Zhao CH,et al.Studyofcombustion andNO chemicalreaction mechanism inammonia blended with DME[J]. Fuel,2022,319: 123832.

    [20]XiaoHH,LiHZ.Experimentalandkinetic modeling studyofthelaminarburningvelocityofNH/DME/airpremixed flames [J].Combustion and Flame,2022,245:112372.

    [21]IssayevG,GiriBR,ElbazA M,etal.Ignitiondelay time and laminarflame sped measurements ofammonia blended with dimethyl ether: a promising low carbon fuel blend[J].Renewable Energy,2022,181: 1353-1370.

    [22]CaiT,ZhaoD.Enhancingandasesingammonia-aircombustionperformancebyblendingwithdimethylether[J].Renewable and Sustainable Energy Reviews,2022,156:112003.

    [23]GlarborgP,MilerJA,Ruscic B,etal.Modelingnitrogenchemistryincombustion[J].ProgressinEnergyandCombustion Science,2018,67: 31-68.

    [24]Burke U,Somers K P,O'Toole P,etal.Anignition delayand kinetic modeling studyof methane,dimethylether,and their mixtures at high pressures[J].Combustion and Flame,2015,162(2):315-330.

    [25]ShrestaKP,SeidelL,ZeuchT,etal.Detailedkineticmechanismfortheoxidationofammoniaincludingtheformationand reduction of nitrogen oxides[J]. Energy amp; Fuels,2018, 32(10): 10202-10217.

    (編輯 陳移峰)

    猜你喜歡
    火焰自由基燃料
    父親的火把
    金沙江文藝(2025年7期)2025-08-18 00:00:00
    火命的民族 (組詩(shī))
    金沙江文藝(2025年7期)2025-08-18 00:00:00
    蘑菇叢
    金沙江文藝(2025年7期)2025-08-18 00:00:00
    TiO2/BC 光催化活化過(guò)硫酸鹽高效降解鹽酸四環(huán)素
    遼寧化工(2025年7期)2025-08-18 00:00:00
    微波預(yù)處理對(duì)水酶法提取辣木籽油抗氧化活性和品質(zhì)的影響
    森林火災(zāi)燃燒特性對(duì)滅火劑選擇的影響研究
    消防界(2025年4期)2025-08-18 00:00:00
    500bar GDi&Flex Fuel擎醇智行與擎穩(wěn)致遠(yuǎn)
    愛達(dá)荷國(guó)家實(shí)驗(yàn)室測(cè)試高燃耗燃料
    以燃料低碳化為核心的商用車技術(shù)路徑發(fā)展預(yù)期
    專用汽車(2025年7期)2025-08-15 00:00:00
    營(yíng)運(yùn)貨車綜合燃料消耗量不確定度的評(píng)定分析
    專用汽車(2025年7期)2025-08-15 00:00:00
    精品欧美国产一区二区三| 三级毛片av免费| 最近在线观看免费完整版| 日本一二三区视频观看| 精品久久蜜臀av无| 久久午夜综合久久蜜桃| 99久久99久久久精品蜜桃| 曰老女人黄片| 国产男靠女视频免费网站| 18禁美女被吸乳视频| 国产精品免费一区二区三区在线| 亚洲狠狠婷婷综合久久图片| 青草久久国产| 久久这里只有精品19| 成人鲁丝片一二三区免费| 一卡2卡三卡四卡精品乱码亚洲| 级片在线观看| 国产乱人伦免费视频| 日本免费一区二区三区高清不卡| 九色成人免费人妻av| 一个人免费在线观看电影 | 国产成人av激情在线播放| 欧美日本视频| 一二三四在线观看免费中文在| 老汉色∧v一级毛片| 国产爱豆传媒在线观看| 国产私拍福利视频在线观看| 一夜夜www| 男插女下体视频免费在线播放| 99久久精品国产亚洲精品| 免费在线观看日本一区| 一个人看的www免费观看视频| 久久久久久九九精品二区国产| 精华霜和精华液先用哪个| 亚洲av美国av| 观看美女的网站| 免费人成视频x8x8入口观看| 狂野欧美激情性xxxx| 在线播放国产精品三级| 国产亚洲精品av在线| 亚洲av五月六月丁香网| 男人舔女人下体高潮全视频| 18禁美女被吸乳视频| 欧美成狂野欧美在线观看| 97超级碰碰碰精品色视频在线观看| 久久人妻av系列| 天天一区二区日本电影三级| 18禁黄网站禁片免费观看直播| 国产精品日韩av在线免费观看| 免费观看的影片在线观看| 国产单亲对白刺激| 搞女人的毛片| 欧美一区二区精品小视频在线| 国产欧美日韩一区二区三| 少妇裸体淫交视频免费看高清| 别揉我奶头~嗯~啊~动态视频| 成人18禁在线播放| 午夜两性在线视频| 18禁国产床啪视频网站| 精品国产乱码久久久久久男人| 亚洲人成电影免费在线| 丰满人妻熟妇乱又伦精品不卡| 91在线观看av| 精品人妻1区二区| 国产午夜精品久久久久久| 成人无遮挡网站| 亚洲第一欧美日韩一区二区三区| 国产精品乱码一区二三区的特点| 国产一区在线观看成人免费| 老鸭窝网址在线观看| 特大巨黑吊av在线直播| 两性午夜刺激爽爽歪歪视频在线观看| 97碰自拍视频| 哪里可以看免费的av片| 日韩有码中文字幕| 国产精品一区二区免费欧美| 在线观看免费视频日本深夜| 午夜a级毛片| 国产精品精品国产色婷婷| 1024手机看黄色片| 午夜福利成人在线免费观看| 十八禁人妻一区二区| 亚洲欧美日韩高清专用| 人妻夜夜爽99麻豆av| 精品无人区乱码1区二区| АⅤ资源中文在线天堂| 欧美性猛交黑人性爽| 哪里可以看免费的av片| 欧美日本视频| 国产成人精品久久二区二区免费| 中国美女看黄片| 久久99热这里只有精品18| 五月伊人婷婷丁香| 日本五十路高清| 国产极品精品免费视频能看的| 欧美三级亚洲精品| 啦啦啦韩国在线观看视频| 99国产精品一区二区三区| 好看av亚洲va欧美ⅴa在| 亚洲欧美日韩东京热| 九九在线视频观看精品| 熟女电影av网| 午夜福利在线在线| 久久香蕉精品热| 日韩免费av在线播放| 精品人妻1区二区| 欧美又色又爽又黄视频| 亚洲av成人不卡在线观看播放网| 99热只有精品国产| 国产黄a三级三级三级人| 久久久久久久午夜电影| 最近最新中文字幕大全免费视频| 母亲3免费完整高清在线观看| 在线看三级毛片| 国产探花在线观看一区二区| 中文亚洲av片在线观看爽| 日韩欧美精品v在线| 国产成人福利小说| 中亚洲国语对白在线视频| 日本撒尿小便嘘嘘汇集6| 午夜精品在线福利| 在线免费观看不下载黄p国产 | www.熟女人妻精品国产| 久久久久免费精品人妻一区二区| 国产精品1区2区在线观看.| ponron亚洲| 人人妻人人澡欧美一区二区| 他把我摸到了高潮在线观看| 中文字幕人成人乱码亚洲影| 欧美黄色片欧美黄色片| a级毛片在线看网站| 亚洲天堂国产精品一区在线| 国内少妇人妻偷人精品xxx网站 | 麻豆成人av在线观看| 亚洲国产看品久久| 色精品久久人妻99蜜桃| 99久久成人亚洲精品观看| 免费无遮挡裸体视频| 1024香蕉在线观看| 99在线视频只有这里精品首页| 99久久精品一区二区三区| 国产三级在线视频| 母亲3免费完整高清在线观看| 天堂影院成人在线观看| 午夜免费成人在线视频| 天天躁日日操中文字幕| 欧美国产日韩亚洲一区| 午夜福利在线观看免费完整高清在 | 看黄色毛片网站| 午夜影院日韩av| 久久天躁狠狠躁夜夜2o2o| 黑人欧美特级aaaaaa片| 日韩免费av在线播放| 制服丝袜大香蕉在线| 国产成人啪精品午夜网站| 亚洲 欧美 日韩 在线 免费| 午夜亚洲福利在线播放| 级片在线观看| 亚洲乱码一区二区免费版| 亚洲av第一区精品v没综合| 国产成人啪精品午夜网站| 日本黄色视频三级网站网址| 免费看美女性在线毛片视频| 在线播放国产精品三级| 好男人电影高清在线观看| 亚洲九九香蕉| 最新美女视频免费是黄的| 国产爱豆传媒在线观看| 免费看十八禁软件| 国产精品99久久久久久久久| 天堂√8在线中文| 女人被狂操c到高潮| 在线观看66精品国产| 国产成人福利小说| 在线播放国产精品三级| 人人妻,人人澡人人爽秒播| 九九热线精品视视频播放| 日本a在线网址| 少妇人妻一区二区三区视频| 特级一级黄色大片| 99视频精品全部免费 在线 | 亚洲真实伦在线观看| 色综合亚洲欧美另类图片| 午夜精品在线福利| 精品福利观看| 麻豆国产97在线/欧美| 欧美3d第一页| 欧美三级亚洲精品| cao死你这个sao货| 国产成人精品久久二区二区免费| 亚洲,欧美精品.| 非洲黑人性xxxx精品又粗又长| 久久久久国产精品人妻aⅴ院| 久久久久久久久久黄片| 亚洲男人的天堂狠狠| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精品乱码久久久v下载方式 | 亚洲av电影在线进入| 成人三级做爰电影| 大型黄色视频在线免费观看| 99riav亚洲国产免费| 久久精品国产综合久久久| 草草在线视频免费看| 国产真实乱freesex| 不卡av一区二区三区| 欧美黑人欧美精品刺激| 熟女少妇亚洲综合色aaa.| 三级国产精品欧美在线观看 | 曰老女人黄片| 国产亚洲欧美在线一区二区| 欧美大码av| 午夜影院日韩av| 毛片女人毛片| 在线十欧美十亚洲十日本专区| 欧美黑人欧美精品刺激| 国产精品永久免费网站| 人人妻人人澡欧美一区二区| 国产真实乱freesex| 亚洲av五月六月丁香网| 国产精品自产拍在线观看55亚洲| 啦啦啦韩国在线观看视频| 一a级毛片在线观看| 亚洲国产色片| 国产黄片美女视频| 亚洲av第一区精品v没综合| 日韩大尺度精品在线看网址| 久久婷婷人人爽人人干人人爱| 精品一区二区三区四区五区乱码| 身体一侧抽搐| 18美女黄网站色大片免费观看| netflix在线观看网站| 黄片小视频在线播放| e午夜精品久久久久久久| 中文亚洲av片在线观看爽| 欧美日韩乱码在线| 亚洲成av人片在线播放无| 丰满人妻一区二区三区视频av | 黑人欧美特级aaaaaa片| 久久99热这里只有精品18| 很黄的视频免费| 99视频精品全部免费 在线 | 日本免费一区二区三区高清不卡| 给我免费播放毛片高清在线观看| 熟妇人妻久久中文字幕3abv| 90打野战视频偷拍视频| 日韩大尺度精品在线看网址| 国产精品免费一区二区三区在线| 国产精品 欧美亚洲| 男女午夜视频在线观看| 三级毛片av免费| 日韩中文字幕欧美一区二区| 最近最新免费中文字幕在线| 757午夜福利合集在线观看| 欧美在线一区亚洲| 淫妇啪啪啪对白视频| 成人三级黄色视频| 夜夜夜夜夜久久久久| 999精品在线视频| cao死你这个sao货| 国产精品电影一区二区三区| 久久久成人免费电影| 国产亚洲欧美98| 亚洲成人久久性| 最近最新免费中文字幕在线| 国产精品一区二区三区四区免费观看 | 亚洲狠狠婷婷综合久久图片| 一级毛片精品| 免费高清视频大片| 国产aⅴ精品一区二区三区波| 精品久久久久久,| 欧美日韩一级在线毛片| 国产99白浆流出| 亚洲aⅴ乱码一区二区在线播放| 国产成人精品久久二区二区免费| 女生性感内裤真人,穿戴方法视频| 色尼玛亚洲综合影院| 在线视频色国产色| xxxwww97欧美| 国产精品影院久久| 在线免费观看不下载黄p国产 | 又爽又黄无遮挡网站| 国产一区二区三区视频了| 亚洲五月天丁香| 特级一级黄色大片| 日韩欧美免费精品| 国产午夜精品论理片| 国产黄a三级三级三级人| 日韩欧美国产在线观看| 国产蜜桃级精品一区二区三区| 成人鲁丝片一二三区免费| 国产乱人视频| 两性夫妻黄色片| 亚洲av五月六月丁香网| 69av精品久久久久久| 久久久久免费精品人妻一区二区| 舔av片在线| 精品国内亚洲2022精品成人| АⅤ资源中文在线天堂| 日韩大尺度精品在线看网址| 丰满人妻一区二区三区视频av | 国产欧美日韩精品一区二区| 亚洲av美国av| 级片在线观看| 一级a爱片免费观看的视频| 免费在线观看亚洲国产| 男人舔奶头视频| 麻豆国产av国片精品| 国产精品国产高清国产av| 欧美性猛交╳xxx乱大交人| 国产精品野战在线观看| 久久久久久国产a免费观看| 日本一本二区三区精品| 国产成人系列免费观看| 成人鲁丝片一二三区免费| 国产三级中文精品| x7x7x7水蜜桃| 99国产极品粉嫩在线观看| 久99久视频精品免费| 国产一级毛片七仙女欲春2| 亚洲精品久久国产高清桃花| 在线观看66精品国产| 免费人成视频x8x8入口观看| 成年版毛片免费区| a级毛片在线看网站| 国产毛片a区久久久久| 国产精品av久久久久免费| 国产1区2区3区精品| 国产高清视频在线观看网站| 中文字幕久久专区| 亚洲精品国产精品久久久不卡| 久久国产精品人妻蜜桃| 日韩欧美三级三区| 欧美一级a爱片免费观看看| 久久中文看片网| 亚洲精品乱码久久久v下载方式 | 亚洲五月婷婷丁香| 熟女电影av网| 欧美成狂野欧美在线观看| 亚洲熟女毛片儿| 这个男人来自地球电影免费观看| 亚洲国产高清在线一区二区三| 久久久久国内视频| 亚洲电影在线观看av| 国模一区二区三区四区视频 | 在线观看免费午夜福利视频| 91字幕亚洲| 欧美日韩福利视频一区二区| 久久久久久九九精品二区国产| 国产淫片久久久久久久久 | 日日摸夜夜添夜夜添小说| АⅤ资源中文在线天堂| 99热精品在线国产| 19禁男女啪啪无遮挡网站| 看免费av毛片| 欧美一区二区国产精品久久精品| 国产亚洲精品综合一区在线观看| 中文字幕熟女人妻在线| 国产伦在线观看视频一区| 九九热线精品视视频播放| 小蜜桃在线观看免费完整版高清| 国产精品精品国产色婷婷| 欧美日本视频| 亚洲av第一区精品v没综合| 一个人看视频在线观看www免费 | 日本在线视频免费播放| 日韩欧美一区二区三区在线观看| 国产高清videossex| 国产伦精品一区二区三区视频9 | 亚洲av片天天在线观看| 我要搜黄色片| 国产精品国产高清国产av| 免费人成视频x8x8入口观看| 精品不卡国产一区二区三区| 久久九九热精品免费| 亚洲中文字幕一区二区三区有码在线看 | www日本在线高清视频| 91在线精品国自产拍蜜月 | 日本 av在线| 国产精品 欧美亚洲| 国产精品久久久人人做人人爽| 在线播放国产精品三级| 99精品欧美一区二区三区四区| 久久久久久久久免费视频了| 国产精品九九99| 美女被艹到高潮喷水动态| 亚洲精品粉嫩美女一区| 国产亚洲精品综合一区在线观看| 制服丝袜大香蕉在线| 亚洲va日本ⅴa欧美va伊人久久| 狠狠狠狠99中文字幕| 老熟妇乱子伦视频在线观看| 十八禁人妻一区二区| 亚洲精品在线美女| 成人性生交大片免费视频hd| 97碰自拍视频| 久久草成人影院| 日韩 欧美 亚洲 中文字幕| 成人三级做爰电影| 久久精品国产清高在天天线| 精品一区二区三区av网在线观看| 国产高清视频在线观看网站| 一区二区三区国产精品乱码| 日韩欧美免费精品| h日本视频在线播放| 国产精品 国内视频| 99视频精品全部免费 在线 | 啦啦啦免费观看视频1| 色av中文字幕| 久久久成人免费电影| 最近最新中文字幕大全电影3| 亚洲,欧美精品.| 一级毛片女人18水好多| 久久久精品欧美日韩精品| 欧美乱色亚洲激情| 母亲3免费完整高清在线观看| 在线观看免费视频日本深夜| 午夜精品久久久久久毛片777| 国产乱人视频| 久久草成人影院| 精品国产三级普通话版| 亚洲在线自拍视频| 亚洲av电影在线进入| 国产成+人综合+亚洲专区| 国产成人精品久久二区二区91| 国语自产精品视频在线第100页| www.www免费av| 欧美三级亚洲精品| 美女 人体艺术 gogo| 伊人久久大香线蕉亚洲五| 最新在线观看一区二区三区| 国内揄拍国产精品人妻在线| xxx96com| 国产一区二区在线观看日韩 | 看片在线看免费视频| 国产精品爽爽va在线观看网站| 欧美日韩国产亚洲二区| 日本 欧美在线| 动漫黄色视频在线观看| 亚洲aⅴ乱码一区二区在线播放| www日本在线高清视频| 国产精品久久视频播放| 两个人的视频大全免费| xxxwww97欧美| 久久这里只有精品中国| 怎么达到女性高潮| 亚洲九九香蕉| 变态另类成人亚洲欧美熟女| 国产真实乱freesex| 亚洲精品一区av在线观看| 色av中文字幕| 欧美黑人欧美精品刺激| 99久久99久久久精品蜜桃| 欧美成人性av电影在线观看| 精品不卡国产一区二区三区| 亚洲成人久久爱视频| 久久国产乱子伦精品免费另类| 法律面前人人平等表现在哪些方面| 亚洲激情在线av| 啦啦啦观看免费观看视频高清| 嫩草影院精品99| 国产精品一区二区三区四区免费观看 | 国产麻豆成人av免费视频| 午夜福利18| 成人18禁在线播放| 国产一区二区在线av高清观看| 国产精品亚洲美女久久久| 嫩草影院入口| 国产爱豆传媒在线观看| 欧美一级毛片孕妇| 美女高潮喷水抽搐中文字幕| 中文资源天堂在线| 18美女黄网站色大片免费观看| 悠悠久久av| 亚洲成a人片在线一区二区| 午夜福利欧美成人| 国产精品一区二区免费欧美| 最近视频中文字幕2019在线8| 亚洲最大成人中文| 欧美乱妇无乱码| 国产av在哪里看| 欧美极品一区二区三区四区| 岛国视频午夜一区免费看| 亚洲午夜理论影院| 精品久久久久久久末码| 午夜视频精品福利| 狂野欧美激情性xxxx| 91麻豆精品激情在线观看国产| 国产激情偷乱视频一区二区| 久久久久久久精品吃奶| 两个人的视频大全免费| 99久国产av精品| 黄片小视频在线播放| 亚洲av片天天在线观看| 亚洲五月天丁香| 午夜福利高清视频| 特大巨黑吊av在线直播| 一二三四社区在线视频社区8| 香蕉久久夜色| 成人无遮挡网站| 又粗又爽又猛毛片免费看| 欧美另类亚洲清纯唯美| 一二三四社区在线视频社区8| 国产伦一二天堂av在线观看| 一边摸一边抽搐一进一小说| 又粗又爽又猛毛片免费看| 精品99又大又爽又粗少妇毛片 | 免费看光身美女| 午夜福利视频1000在线观看| 国产野战对白在线观看| 最近最新中文字幕大全电影3| 一个人免费在线观看的高清视频| 一二三四在线观看免费中文在| 国产成人aa在线观看| 五月伊人婷婷丁香| 精品国内亚洲2022精品成人| 最近最新中文字幕大全免费视频| 1000部很黄的大片| 亚洲欧美精品综合久久99| 美女高潮喷水抽搐中文字幕| 亚洲专区字幕在线| 亚洲成av人片免费观看| 久久国产乱子伦精品免费另类| 久久人妻av系列| av在线天堂中文字幕| 99视频精品全部免费 在线 | 91麻豆精品激情在线观看国产| 色综合婷婷激情| 麻豆av在线久日| 亚洲 国产 在线| 日本三级黄在线观看| 视频区欧美日本亚洲| 国产黄片美女视频| 非洲黑人性xxxx精品又粗又长| 久久婷婷人人爽人人干人人爱| 波多野结衣巨乳人妻| 欧美黑人巨大hd| 老熟妇仑乱视频hdxx| 人妻久久中文字幕网| 亚洲熟妇熟女久久| 两性午夜刺激爽爽歪歪视频在线观看| 真人一进一出gif抽搐免费| 精华霜和精华液先用哪个| 一级毛片女人18水好多| 亚洲aⅴ乱码一区二区在线播放| 国产精品综合久久久久久久免费| 熟女少妇亚洲综合色aaa.| 国产精品九九99| 国产精品98久久久久久宅男小说| 18禁黄网站禁片午夜丰满| 九九久久精品国产亚洲av麻豆 | 波多野结衣高清作品| 日本熟妇午夜| 操出白浆在线播放| 国产蜜桃级精品一区二区三区| 黑人操中国人逼视频| 男插女下体视频免费在线播放| 在线视频色国产色| 99热精品在线国产| 夜夜爽天天搞| av天堂在线播放| 真实男女啪啪啪动态图| 国产伦人伦偷精品视频| 悠悠久久av| 午夜成年电影在线免费观看| 亚洲狠狠婷婷综合久久图片| 亚洲国产高清在线一区二区三| 精品久久蜜臀av无| 亚洲国产欧美人成| 小说图片视频综合网站| 性欧美人与动物交配| АⅤ资源中文在线天堂| 一个人看的www免费观看视频| 无人区码免费观看不卡| 香蕉av资源在线| 欧美乱码精品一区二区三区| 欧美成狂野欧美在线观看| 国产一区二区在线观看日韩 | 国产乱人伦免费视频| 免费在线观看亚洲国产| 深夜精品福利| 成年版毛片免费区| 国产视频内射| 亚洲国产欧洲综合997久久,| 老司机在亚洲福利影院| 美女 人体艺术 gogo| 精品乱码久久久久久99久播| www日本在线高清视频| 国产麻豆成人av免费视频| 国产私拍福利视频在线观看| 不卡av一区二区三区| 一进一出抽搐动态| 欧美午夜高清在线| 99国产精品一区二区三区| 国产美女午夜福利| 欧美午夜高清在线| 亚洲va日本ⅴa欧美va伊人久久| 操出白浆在线播放| 色吧在线观看| 少妇熟女aⅴ在线视频| 熟女人妻精品中文字幕| 在线永久观看黄色视频| 又黄又粗又硬又大视频| 别揉我奶头~嗯~啊~动态视频| 丰满的人妻完整版| 免费在线观看视频国产中文字幕亚洲| 麻豆成人午夜福利视频| 国产欧美日韩一区二区三| 不卡av一区二区三区| 午夜激情欧美在线| 亚洲18禁久久av| 亚洲欧美精品综合一区二区三区| 美女高潮的动态| 亚洲国产欧美人成| 精品国产亚洲在线| 国产精品 欧美亚洲|