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

    鋁粉/氫氣/空氣混合爆轟現(xiàn)象試驗研究1)

    2023-12-16 11:49:04張曉源盧子寅李進(jìn)平張仕忠陸星宇
    力學(xué)學(xué)報 2023年11期
    關(guān)鍵詞:波面鋁粉聲速

    張曉源 盧子寅 李進(jìn)平 張仕忠 陸星宇 陳 宏

    (中國科學(xué)院力學(xué)研究所高溫氣體動力學(xué)國家重點實驗室,北京 100190)

    引言

    有關(guān)粉塵爆炸的試驗研究,始于對工業(yè)安全的考慮.特定條件下粉塵與空氣形成可燃混合物,一旦遇到點火源,就會引發(fā)燃燒并產(chǎn)生爆炸.在19 世紀(jì)中葉就開始有研究煤粉、鋁粉和農(nóng)作物等粉塵爆炸特性的相關(guān)試驗[1-4].為了更好地理解和預(yù)防粉塵爆炸,先后開展了大量的試驗研究.根據(jù)參與化學(xué)反應(yīng)的物質(zhì)形態(tài)的區(qū)別,可將粉末與氣體混合物的爆轟現(xiàn)象分為3 類[5]:(1)多相爆轟,在爆轟過程中發(fā)生的化學(xué)反應(yīng)是固態(tài)顆粒與氣體之間的多相反應(yīng),通常可燃物僅為固態(tài)顆粒;(2)混合爆轟,可燃物既有固態(tài)顆粒又有氣體,在爆轟過程中,既有氣體間的化學(xué)反應(yīng),又有固態(tài)顆粒與氣體的反應(yīng),區(qū)別于多相爆轟,混合爆轟情況更加復(fù)雜;(3)含塵氣體爆轟,固態(tài)顆粒是惰性的,不參加化學(xué)反應(yīng),爆轟過程中反應(yīng)只發(fā)生于氣體之間.

    金屬粉末作為一種新型燃料,具有密度大和燃燒焓高等優(yōu)勢,并且綠色環(huán)保,易于運輸和保存,在航天推進(jìn)中扮演著重要的角色,有效增大了固體推進(jìn)劑的比沖[6-7].而金屬粉末在應(yīng)用中最主要的問題是其往往難以點燃,對點火和起爆的要求高.金屬粉末與燃料混合使用,將有希望克服這樣的缺點.作為新型空天技術(shù)的發(fā)展方向之一,爆轟推進(jìn)技術(shù)通過爆轟波點燃燃料,與傳統(tǒng)的基于燃燒的推進(jìn)系統(tǒng)相比,可以實現(xiàn)更快、更完全的燃燒,提高發(fā)動機(jī)效率和性能[8-9].目前,在脈沖爆轟、旋轉(zhuǎn)爆轟和斜爆轟等方面的研究與發(fā)展中,爆轟推進(jìn)技術(shù)已經(jīng)取得了實質(zhì)性的進(jìn)展[10-13],但其中固體燃料的應(yīng)用研究相對較少[14-16].有研究表明,將某些金屬粉末與其他燃料混合使用,通過提高火焰速度、溫度和壓力上升率,可以提高爆轟發(fā)動機(jī)的燃燒效率和推力水平[17].當(dāng)金屬粉末和可燃?xì)怏w同時存在,并與空氣混合發(fā)生爆轟時,金屬粉末和可燃?xì)怏w都將反應(yīng)釋熱,其中存在復(fù)雜的物理和化學(xué)過程.金屬粉末在爆轟推進(jìn)技術(shù)上的應(yīng)用基于對爆轟波的理論認(rèn)識,需要充分的試驗研究成果來支撐.

    金屬粉末/空氣爆轟研究領(lǐng)域擁有豐富的試驗研究成果[18-22],這些成果揭示了金屬粉末爆轟的特性、波陣面結(jié)構(gòu)以及爆燃轉(zhuǎn)爆轟等重要現(xiàn)象,加深了對金屬粉末爆轟機(jī)理的理解.針對金屬粉末/氣體混合爆轟現(xiàn)象的研究,Zhang 等[22-23]和Veyssiere等[24-26]兩組研究團(tuán)隊分別采用水平和直立式結(jié)構(gòu)的爆轟管,選用了μm 級的鋁粉,混合氣體多為富氧的乙炔/空氣或者氫氣/空氣進(jìn)行實驗,試驗中發(fā)現(xiàn)了穩(wěn)定的雙波面爆轟波結(jié)構(gòu).在對混合爆轟現(xiàn)象的復(fù)雜性與多樣性進(jìn)行分析之后,分別給出了對混合爆轟波結(jié)構(gòu)的分類,使我們對單波面、雙波面等混合爆轟現(xiàn)象有了基本的認(rèn)識與了解.

    本文關(guān)注高溫水蒸氣為主體的燃?xì)庵袖X粉點火燃燒特性對混合爆轟現(xiàn)象的影響,開展了鋁粉/氫氣/空氣混合爆轟試驗研究.實驗在當(dāng)量比的氫氣/空氣混合氣體中進(jìn)行,氣相爆轟后高溫產(chǎn)物中的主要氧化劑成分是水蒸氣.高溫燃?xì)鈼l件下鋁粉的點火燃燒特性直接對混合爆轟現(xiàn)象產(chǎn)生影響.通過對混合爆轟現(xiàn)象的深入分析,歸納混合爆轟波的結(jié)構(gòu)的影響規(guī)律,了解金屬燃料粉末在復(fù)雜環(huán)境中的點火與燃燒特性.對混合爆轟現(xiàn)象的深入認(rèn)知,將利于爆轟技術(shù)的發(fā)展與應(yīng)用.

    1 混合爆轟現(xiàn)象

    在混合爆轟現(xiàn)象中,燃料既來源于混合物中的可燃?xì)怏w,也來源于散布在氣體中的可燃顆粒,爆轟波依靠這兩部分的化學(xué)釋熱來維持.顆粒相的點火延遲和燃燒時間與氣相的相比,相差一個甚至數(shù)個量級,因而混合爆轟的釋熱過程復(fù)雜.在頭激波之后,化學(xué)反應(yīng)區(qū)域更長,釋熱順序的不同將引起混合爆轟多種不同的波系結(jié)構(gòu).

    根據(jù)ZND 理論,激波后是一段反應(yīng)區(qū),反應(yīng)區(qū)末端是化學(xué)反應(yīng)完成后的爆轟產(chǎn)物[27].混合爆轟中氣相間的化學(xué)反應(yīng)在波陣面后迅速完成,而金屬粉末由于傳熱和相變等物理過程,反應(yīng)慢一步發(fā)生.氣相間爆轟首先產(chǎn)生高溫燃?xì)?之后鋁粉與氣相爆轟產(chǎn)物中的氧化劑反應(yīng)釋熱.富氧的氣體爆轟所得的燃?xì)庵?存在O2,H2O 及CO2等氧化劑.相對于氣相爆轟波,鋁粉點火延遲時間的不同,導(dǎo)致了混合爆轟波結(jié)構(gòu)存在不同的類型.

    根據(jù)實驗結(jié)果與理論分析,Zhang 等[22-23]與Veyssiere 等[24-26]在各自的文獻(xiàn)中,對混合爆轟波的結(jié)構(gòu)進(jìn)行了分類.(1)單波面結(jié)構(gòu):在單波結(jié)構(gòu)的情況下,爆轟波由氣相間的反應(yīng)釋熱和顆粒相的反應(yīng)釋熱同時支持,僅有一個波面.該情況下顆粒相的化學(xué)反應(yīng)延遲時間較短,能夠在氣相爆轟波音速面之前開始反應(yīng).(2)雙波面結(jié)構(gòu):特點是存在兩個波面.第一道波面由氣相間的化學(xué)反應(yīng)釋熱支持,第二道波面由顆粒相與氣相產(chǎn)物的化學(xué)反應(yīng)釋熱支持.顆粒相的反應(yīng)點火延遲在氣相爆轟波聲速面之后,因而不能影響第一道波;但是顆粒相的反應(yīng)釋熱卻足夠劇烈,因此產(chǎn)生了第二道波.(3)偽氣相爆轟:當(dāng)顆粒相的點火延遲時間較長,在爆轟波音速面之后遠(yuǎn)處才發(fā)生顆粒相的化學(xué)反應(yīng),爆轟波就只能得到氣相間化學(xué)反應(yīng)釋熱的支持.對于爆轟波而言,顆粒相的化學(xué)反應(yīng)毫無作用,這種情形稱為偽氣相爆轟.

    在混合爆轟中既含有氣相間的反應(yīng),又含有兩相反應(yīng),而金屬顆粒由于傳熱和相變等物理過程,其反應(yīng)慢一步發(fā)生.對于氫/氧/氮氣體與鋁金屬粉末混合物,在激波后氣相的氫/氧反應(yīng)較快完成,而固體鋁粉顆粒的氧化反應(yīng)則相對慢很多,因而鋁粉反應(yīng)區(qū)較長,爆轟波結(jié)構(gòu)模型如圖1 所示.

    圖1 混合爆轟結(jié)構(gòu)示意圖Fig.1 Schematic diagram of hybrid detonation

    2 試驗方案

    2.1 試驗系統(tǒng)

    混合爆轟試驗系統(tǒng)結(jié)構(gòu)見圖2,采用臥式爆轟管,由爆轟試驗段、點火段、真空泄爆段和揚(yáng)塵系統(tǒng)等構(gòu)成,并配套有真空與供氣系統(tǒng)、測量與控制系統(tǒng).點火段、爆轟試驗段與真空泄爆段三者之間由膜片隔開.爆轟試驗段內(nèi)徑224 mm,長度13 m.在當(dāng)量比的氫氣/空氣混合氣體中氮氣占比較大,由小口徑的射流發(fā)火管直接起始爆轟很困難[28],如果點火能量不足,容易產(chǎn)生爆燃轉(zhuǎn)爆轟的現(xiàn)象,不利于試驗觀測,因此采用等直徑的點火段將試驗段起爆.點火段長3 m,管中充入當(dāng)量比的氫氣與氧氣,易于起爆.由射流發(fā)火管首先將點火段起爆,爆轟波擊破點火段與試驗段之間的膜片后,再起爆試驗段.真空泄爆段用于收集和清理粉塵,并且可以降低反射波壓力峰,提高試驗的安全性.

    圖2 混合爆轟管結(jié)構(gòu)示意圖Fig.2 Schematic diagram of the hybrid detonation tube

    2.2 揚(yáng)塵系統(tǒng)

    使用揚(yáng)塵系統(tǒng)可以獲得分布較均勻的粉塵氣體混合物,圖3 給出了揚(yáng)塵結(jié)構(gòu)剖面圖,U 型粉槽固定在炮管的底部,用于盛放鋁粉,揚(yáng)塵管架于粉槽的正上方.在揚(yáng)塵管正下方開有一列直徑1 mm 的小孔正對粉槽,小孔間隔為5 cm.如圖3(b)所示,小孔噴出氮氣,將粉槽內(nèi)的鋁粉揚(yáng)起.揚(yáng)塵管內(nèi)通入的氮氣由圖2 中的電磁閥進(jìn)行控制.選用2 m 長的有機(jī)玻璃管代替爆轟管進(jìn)行揚(yáng)塵效果試驗驗證,觀測鋁粉的揚(yáng)起和塵降情況.利用遮光法,采用光強(qiáng)傳感器測量粉塵的相對均勻性.根據(jù)實際的揚(yáng)塵效果,確定了揚(yáng)塵氮氣的壓力為2.0 MPa,電磁閥開閥時間600 ms.圖4 給出了試驗過程視頻的截圖,鋁粉塵的濃度約100 g/m3.揚(yáng)塵氮氣噴入前,鋁粉均勻鋪在粉槽內(nèi);電磁閥開啟后氮氣噴入,粉塵沿管底同時躍起;之后粉塵充分散布于管內(nèi),管一端的光源發(fā)出的光線,幾乎被粉塵完全遮擋而變暗.

    圖3 揚(yáng)塵管結(jié)構(gòu)示意圖Fig.3 Schematic diagram of dispersion tube

    圖4 揚(yáng)塵效果試驗視頻截圖Fig.4 Screenshot of dust effect test video

    試驗前管內(nèi)預(yù)先加入鋁粉,根據(jù)鋁粉總質(zhì)量和管內(nèi)體積可估算揚(yáng)塵鋁粉濃度.試驗時首先利用揚(yáng)塵結(jié)構(gòu)將鋁粉揚(yáng)起,電磁閥開閥信號持續(xù)600 ms,閥關(guān)閉后點火.發(fā)火管先將點火段起爆,點火段內(nèi)爆轟波向下游傳播,破膜后引爆試驗段.

    2.3 鋁粉

    考慮到μm 級、nm 級的球形鋁粉在含能材料領(lǐng)域有廣泛的應(yīng)用和發(fā)展前景.采用μm 級和nm 級的球形鋁粉開展試驗,4 種μm 級鋁粉的名義粒徑分別為40 μm,20 μm,10 μm 和1 μm,nm 級鋁粉的名義粒徑100 nm.利用掃描電子顯微鏡觀測鋁粉的粒徑分布情況,同時進(jìn)行能譜分析獲取鋁粉的實際鋁含量.圖5 是通過掃描電子顯微鏡觀測5 種規(guī)格鋁粉獲得的照片.可以看到5 種鋁粉顆粒形貌基本接近球形,粒徑小的鋁粉球形度更好,并且無明顯的團(tuán)簇現(xiàn)象.

    圖5 鋁粉掃描電鏡照與鋁粉顆粒尺寸分布Fig.5 Aluminum powder scanning electron microscope photo and particle size distribution

    圖5 鋁粉掃描電鏡照與鋁粉顆粒尺寸分布 (續(xù))Fig.5 Aluminum powder scanning electron microscope photo and particle size distribution (continued)

    由鋁粉的掃描電鏡照片可以統(tǒng)計得出各規(guī)格鋁粉的粒徑尺寸分布情況,如圖5 所示.名義粒徑40 μm鋁粉,質(zhì)量占比89.0%粒徑在20~50 μm 內(nèi);名義粒徑20 μm 鋁粉有89.6%在13~31 μm 內(nèi);名義粒徑10 μm 鋁粉有88.6% 在6~14 μm 內(nèi);名義粒徑1 μm 鋁粉有88.0% 在0.5~2.1 μm 內(nèi);名義粒徑100 nm 鋁粉有85.6%在70~150 nm 內(nèi).5 種鋁粉的粒徑分布區(qū)別明顯,粒徑的不同直接決定了鋁粉的點火延遲與燃燒特性,從而將對混合爆轟波產(chǎn)生影響.

    微米級鋁粉表面為自然氧化,粒徑100 nm 表面采用烴類包覆處理.通過掃描電鏡能譜分析,給出鋁粉樣本的元素組成.名義粒徑10 μm,20 μm 和40 μm的鋁粉中的鋁含量均達(dá)到98%以上;名義粒徑1 μm鋁粉的鋁含量為91.4%;名義粒徑100 nm 鋁粉活性鋁含量在85%~90%之間.

    3 試驗結(jié)果與分析

    3.1 試驗狀態(tài)

    混合爆轟試驗預(yù)充氣體為當(dāng)量比的氫氣和空氣,因揚(yáng)塵需要充入的氮氣量可依據(jù)試驗段的壓力增量測算得出,試驗中實際H2/O2/N2混合氣體的摩爾比為2:1:4.3.在沒有鋁粉的氣相爆轟對比試驗中,保證混合氣體比例相同.混合氣體初始壓力為1 bar(1 bar=100 kPa),根據(jù)爆轟波理論,在CJ 條件下氣相爆轟波計算結(jié)果見表1,CJ 壓力14.87 bar,理論爆速1923 m/s.點火段的氣體是當(dāng)量比的氫氣和氧氣,初始壓力0.8 bar,略低于試驗段混合氣體,以防止試驗段形成過爆轟.微米級鋁粉濃度在300 g/m3左右,100 nm 鋁粉濃度約150 g/m3,試驗狀態(tài)參數(shù)見表2.

    表1 氣相爆轟理論計算結(jié)果Table 1 Theoretical results of gas phase detonation

    表2 混合爆轟實驗工況Table 2 Hybrid detonation test condition

    3.2 雙波面爆轟實驗結(jié)果

    名義粒徑40 μm 和20 μm 的鋁粉/氫氣/空氣混合爆轟試驗中,混合爆轟波呈現(xiàn)出雙波峰的結(jié)構(gòu).圖6 所示的壓力曲線是沿激波管布置的壓力傳感器所測得的,該壓力曲線給出了爆轟波面經(jīng)過各測點位置而產(chǎn)生的壓力變化.S為測點位置與點火段膜片間的距離,第一個測點緊靠膜片S=0 m,各測點之間間隔2 m,共7 個測點.實驗段起爆后,爆轟波經(jīng)歷了兩個階段:分別是發(fā)展階段和穩(wěn)定階段.試驗開始時,爆轟波為普通的單峰結(jié)構(gòu),見圖6(a) 和圖6(b)中S=0 m 處的壓力曲線.發(fā)展階段爆轟波向下游的傳播過程中,從單峰中逐漸剝離出第二道波峰,并最終形成了穩(wěn)定的雙波面結(jié)構(gòu).

    圖6 壓力曲線測量結(jié)果Fig.6 Pressure records of a hybrid detonation

    名義粒徑40 μm 時,在S=6 m 位置形成了雙波面混合爆轟波結(jié)構(gòu),見圖6(a),圖中同時給出了沒有鋁粉的純氣相爆轟波壓力曲線作為對比.雙波結(jié)構(gòu)第一道波的波速為1.87 km/s,與純氣相爆轟波的波速1.88 km/s 接近,且在混合爆轟波的發(fā)展和穩(wěn)定階段保持勻速.從S=0~6 m 混合爆轟波的發(fā)展中,可以觀察到,第二道波由第一道波中逐漸剝離而出.在形成穩(wěn)定的雙波結(jié)構(gòu)后,第一道波壓力峰值略低于純氣相的爆轟波.由S=6~12 m 的壓力跳躍,兩波之間的距離在逐漸拉大,第二道波波速為1.67 km/s,比第一道波慢.在重復(fù)實驗中,該值為1.69 km/s,實驗結(jié)果具有重復(fù)性.

    名義粒徑20 μm 時,混合爆轟波發(fā)展更快,由圖6(b)可見,在S=4 m 位置第二道波峰就已經(jīng)形成了.第一道波波速為1.86 km/s,同樣與純氣相爆轟的波速接近,第一道波峰壓力略低于純氣相爆轟波.雙波峰間的距離基本不變,第二道波的波速為1.89 km/s,與第一道波接近.20 μm 鋁粉混合爆轟第二道波峰壓力增幅較大,明顯高于40 μm 鋁粉.

    由名義粒徑40 μm 和20 μm 鋁粉的混合爆轟試驗結(jié)果可知,雙波面混合爆轟波結(jié)構(gòu)中第一道波是由氣相反應(yīng)支持的爆轟波,與純氣相爆轟對比試驗獲得的爆速與爆壓接近.對于第一道波,鋁粉并沒有參與化學(xué)反應(yīng)釋熱,相反地還從中吸收能量,增加了自己的動能與熱能,所以第一道波爆速與爆壓會略低于純氣相爆轟.粉塵被燃?xì)饬骷訜峒铀?將從氣流中吸收熱能和動能,對爆轟波而言是能量的損耗,從而造成了爆速降低.在粉塵的質(zhì)量濃度相同時,粒徑越小動能與熱量的傳遞越快,粉塵吸收的能量越多,從而爆轟波的能量損失越大.經(jīng)過第一道波后,氣相的爆轟產(chǎn)物近似達(dá)到化學(xué)平衡狀態(tài),由于氫和空氣是當(dāng)量比的,所以主要的氧化劑是H2O,以及微量的O2,NO 等.因此,第一道波后的混合物是鋁粉與氣相爆轟產(chǎn)物的混合.隨著爆轟產(chǎn)物高溫燃?xì)鈱︿X粉的加熱,當(dāng)鋁粉達(dá)到了點火條件便發(fā)生了點火釋熱,第二道波就是由鋁粉化學(xué)反應(yīng)釋熱而產(chǎn)生的.鋁粉粒徑的減小,增大了接觸面,加速了鋁粉與燃?xì)忾g的熱量傳遞,因此鋁粉的點火延遲時間縮短.正如圖1 所示,鋁粉的點火延遲逐漸接近第一道波,但只要鋁粉的釋熱發(fā)生于聲速面之后,就不會影響到第一道波.

    3.3 單波面爆轟實驗結(jié)果

    在名義粒徑10 μm,1 μm 和100 nm 的實驗中,鋁粉/氫氣/空氣混合爆轟波不再是雙波面結(jié)構(gòu),壓力曲線表現(xiàn)為單波峰,如圖7 所示.由圖7(a)中試驗段管壁上的壓力測量數(shù)據(jù)可得,10 μm 鋁粉混合爆轟波波速1.84 km/s,略低于純氣相爆速,但是爆轟波壓力峰值有明顯的增加,增幅達(dá)到約30%.爆速沒有增加,說明激波的支持依然只受到氣相化學(xué)反應(yīng)釋熱,并沒有其他新的能量來源對激波產(chǎn)生影響.而壓力峰值卻大幅地增加,又說明了鋁粉點火燃燒確實發(fā)生了.出現(xiàn)這種情形只能是鋁粉在氣相爆轟波聲速面之后發(fā)生了點火燃燒,但是由于雙波結(jié)構(gòu)間隔很近幾乎重合,在壓力曲線上無法辨識.由圖6 雙波結(jié)構(gòu)可以看到,隨著粒徑的減小,第二道波的速度是增加的,在名義粒徑20 μm 時,雙波面間隔就穩(wěn)定在了0.3 ms 左右.氫氧氣相爆轟CJ 聲速面距離激波面僅為毫米量級[27,29],當(dāng)前激波速度1.9 km/s 的情況下,時間間隔約0.5 μs.如果雙波面的第二道波緊貼第一道聲速面,即使有理想的數(shù)據(jù)采集條件,壓力傳感器的端面尺寸至少得小于1 mm 才有可能捕捉的到兩次壓力階躍,目前還沒有這樣的實驗測試能力.由圖7(a)壓力曲線可見,10 μm 鋁粉的點燃進(jìn)一步提前,但依然在聲速面之后,不能支持第一道激波.雖然試驗測量中無法辨識,但是可以判定該情況下雙波結(jié)構(gòu)緊貼,在試驗壓力曲線上發(fā)生了重合,呈現(xiàn)出單波結(jié)構(gòu).

    圖7 壓力曲線測量結(jié)果Fig.7 Pressure records of a hybrid detonation

    1 μm 鋁粉也表現(xiàn)為單波面爆轟波,但仔細(xì)比較可以發(fā)現(xiàn)存在明顯的區(qū)別.由圖7(b)所示,比較有無鋁粉時的實驗結(jié)果可以發(fā)現(xiàn),加入鋁粉后,爆速有了明顯增加,爆轟波從S=0 m 到達(dá)S=12 m 所需要的時間為6.23 ms,比無鋁粉時減少了0.16 ms,爆轟波速度達(dá)到了1.93 km/s.爆速的增加是由于激波面除了受到氣相化學(xué)反應(yīng)釋熱的支持以外,有了新的能量來源,唯一可能的能量來源就是鋁粉的燃燒釋熱.在聲速面之前就已經(jīng)發(fā)生了鋁粉的點火,并且開始釋放出化學(xué)能.但爆速的增加也只是小幅的增加,鋁粉的燃燒釋熱在聲速面前只是部分的釋出.從壓力峰值可以看出,1 μm 鋁粉混合爆轟波壓力峰值增量比10 μm 鋁粉小,僅10%左右.

    名義粒徑100 nm 鋁粉的表面由烴類包覆,作為納米級鋁粉且無氧化層,其點火延遲更短.但試驗結(jié)果中并沒有得到爆速與爆壓都明顯增加的混合爆轟波.爆速1.91 km/s,比沒有鋁粉的純氣相爆轟略快,但沒有快于1 μm 鋁粉;壓力峰值比1 μm 鋁粉高,但低于10 μm 鋁粉的壓力峰值.該現(xiàn)象表明100 nm 鋁粉也只有少部分釋熱對激波產(chǎn)生支持,由于激波音速面距離很短,即使是達(dá)到了納米級的鋁粉,也不能在1 mm 的音速面之前將完整的化學(xué)反應(yīng)釋熱過程完成.名義粒徑10 μm,1 μm 和100 nm 的混合爆轟波雖然試驗得到的都是單波面結(jié)構(gòu),但爆速和爆壓上區(qū)別明顯.鋁粉點燃時刻與聲速面的位置關(guān)系以及鋁粉燃燒速率決定了壓力峰值大小.壓力峰值是任何爆轟裝置都需要關(guān)注的參數(shù),試驗結(jié)果表明相同條件下混合爆轟波并不是金屬粉塵粒徑越小,壓力峰值就越高.在本文的研究條件范圍內(nèi),當(dāng)鋁粉點火延遲位于第一道聲速面時,將產(chǎn)生最大的壓力峰值.

    3.4 混合爆轟波結(jié)構(gòu)分析

    將鋁粉/氫氣/空氣混合爆轟波的測試結(jié)果歸納于表3 之中.鋁粉粒徑較大時,混合爆轟波呈現(xiàn)雙波面結(jié)構(gòu);粒徑較小時為單波面的結(jié)構(gòu).基于試驗結(jié)果的分析表明,鋁粉的點火燃燒特性,影響了混合爆轟波結(jié)構(gòu).特別是鋁粉的點火延遲相對于氣相爆轟聲速面的位置關(guān)系,對爆轟波結(jié)構(gòu)起了決定性作用.當(dāng)鋁粉的點火延遲遠(yuǎn)離氣相爆轟聲速面,可以得到雙波面混合爆轟波結(jié)構(gòu);當(dāng)其點火延遲接近氣相爆轟聲速面但稍有落后時,由于傳感器尺度與響應(yīng)有限,測量結(jié)果為單波峰,波爆速沒有增加但壓力峰值增加很大;當(dāng)鋁粉點火延遲在氣相爆轟聲速面之前,鋁粉的燃燒對第一道激波產(chǎn)生了支持,混合爆轟為單波面結(jié)果,爆速與壓力峰值都較純氣相狀態(tài)有明顯的提升.

    如圖1 所示的鋁粉/氫氣/空氣混合爆轟波結(jié)構(gòu)中,混合均勻的氫氣與空氣在頭激波的反應(yīng)迅速,而鋁粉則慢一步與高溫的氣相燃燒產(chǎn)物H2O 繼續(xù)反應(yīng).為了說明鋁粉在高溫燃燒產(chǎn)物H2O 中的點火延遲對混合爆轟波結(jié)構(gòu)起到的關(guān)鍵作用.基于一維的ZND 模型,針對如圖1 所示的定常的混合爆轟波雙波結(jié)構(gòu),計算混合爆轟波中鋁粉的點火延遲.建立氣相和顆粒相各自的守恒方程,引入交互項以表征兩項之間的動量交換和熱傳導(dǎo)[30-31].計算模型中假設(shè):鋁粉顆粒為球形,尺寸相同,溫度單一且均勻;顆粒體積忽略,氣相為完全氣體;不考慮管壁的摩擦與傳熱.顆粒相守恒方程

    其中,σp和vp是質(zhì)量濃度與顆粒速度,ep是顆粒內(nèi)能,f和q分別表示顆粒和氣體之間的動能交換和熱傳導(dǎo).氣相守恒方程

    氣體的狀態(tài)方程為

    式中,ρg和vg是氣體的濃度與速度,γ和R分別為氣體比熱比和氣體常數(shù),Qgas為氣相反應(yīng)釋熱量.采用二階段化學(xué)反應(yīng)模型來描述氣相反應(yīng)[32]

    其中Qg為氣相化學(xué)反應(yīng)放熱,β表征了放熱反應(yīng)階段的進(jìn)程.

    動量交換項和熱傳導(dǎo)項分別為

    其中,dp為顆粒直徑,Cd是粒子的拖拽力系數(shù),Nu是努塞爾數(shù),λg是氣體導(dǎo)熱系數(shù),ε是鋁粉顆粒的輻射率,σ是玻爾茲曼常數(shù).

    采用迭代求解,計算過程中首先假定爆速,代入由式(1)~式(3)組成的常微分方程組中進(jìn)行求解.通過非理想爆轟的廣義CJ 條件[30],反復(fù)校正爆速,最終得到確定且唯一的爆速,從而得出鋁粉的動量與能量交換過程.圖8 給出了10 μm 鋁粉顆粒溫度TAl的計算結(jié)果,爆速取值為1990.1 m/s 時,滿足非理想爆轟的廣義CJ 條件,馬赫數(shù)為1 時有效釋熱率dQeff/dt=0,即化學(xué)反應(yīng)放熱量與熱損失達(dá)到平衡.聲速面在0.96 mm 處,聲速面顆粒溫度上升至約340 K.

    圖8 10 μm 鋁粉顆粒溫度的計算結(jié)果Fig.8 Calculation result of TAl (10 μm)

    選取兩個代表溫度來分析鋁粉顆粒的點火,一是鋁的熔點933 K,另一個是綜合考慮相關(guān)研究成果,給出的μm 級鋁粉的點火溫度1350 K[33].顆粒受氣體加熱到達(dá)點火溫度所需要的波后距離即為點火延遲距離.圖9 給出了不同粒徑鋁粉,在混合爆轟中達(dá)到熔點933 K 和點火溫度1350 K,所需要的點火延遲距離.很明顯隨粒徑的增大熱量傳導(dǎo)速度減慢.在混合爆轟波中,2 μm 鋁粉顆粒僅需幾毫米就達(dá)到了點火溫度,已經(jīng)接近氣相氫氧爆轟CJ 面1 mm左右的距離;而50 μm 鋁粉顆粒則需要0.34 m 才能到達(dá)熔點933 K,需要0.78 m 到達(dá)點火溫度1350 K.根據(jù)圖9 中的鋁粉點火延遲曲線,可以更深入地分析混合爆轟試驗中得到的爆轟波結(jié)構(gòu).

    圖9 不同粒徑鋁粉顆粒的點火延遲Fig.9 Ignition delay of aluminum powder with different particle sizes

    鋁粉顆粒2 μm 的點火延遲距離已經(jīng)接近氣相氫氧爆轟CJ 面,隨著粒徑的繼續(xù)減小,必然鋁粉顆??梢栽诼曀倜嬷斑_(dá)到點火溫度并被點燃.結(jié)合試驗結(jié)果,鋁粉粒徑較小的情況下(100 nm 和1 μm),對比同樣狀態(tài)下的純氣相爆轟,鋁粉的加入使得波后壓力增高,爆速有所增加.可見,鋁粉點火釋熱開始于聲速面之前,因而爆轟波由氣相間和顆粒相兩部分的反應(yīng)釋熱支持,燃?xì)?鋁粉混合爆轟呈現(xiàn)單波面結(jié)構(gòu).

    在鋁粉粒徑較大的情況下(20 μm 和40 μm),由圖9 的計算結(jié)果,顆粒在第一聲速面之后點火.實驗結(jié)果出現(xiàn)雙波面結(jié)構(gòu),氣相間化學(xué)反應(yīng)釋熱支持第一道波,而鋁粉反應(yīng)釋熱支持第二道波.

    在鋁粉粒徑10 μm 情況下,實驗所得曲線雖然也是單波面結(jié)構(gòu),但根據(jù)鋁粉點火延遲的計算結(jié)果,此時鋁粉顆粒的燃燒反應(yīng)發(fā)生于第一聲速面之后,因而其釋熱對爆轟波不產(chǎn)生支持,雖然峰值壓力有大幅提高,但爆速并不增加.因為10 μm 鋁粉的點火延遲距離在10 mm 量級,尺寸上接近了壓力傳感器,所以測量結(jié)果無法分辨,但其實也是雙波面結(jié)構(gòu),是兩波面距離很近的混合爆轟波.結(jié)合鋁粉點火延遲計算結(jié)果來分析混合爆轟現(xiàn)象,可更加清楚地知道,爆轟波壓力曲線的變化可以明顯反映鋁粉在燃?xì)庵械狞c火與燃燒;另一方面,鋁粉在燃?xì)庵械狞c火燃燒特性決定了混合爆轟波的結(jié)構(gòu).

    4 結(jié)論

    鋁粉/氫氣/空氣混合爆轟現(xiàn)象具有復(fù)雜性和多樣性,通過鋁粉/氫氣/空氣混合爆轟波試驗,對于不同粒徑的球形鋁粉顆粒,獲得了單波面和雙波面的混合爆轟波結(jié)構(gòu).當(dāng)量比的混合均勻的氫氣與空氣,氣相反應(yīng)在激波面后迅速完成,因而鋁粉的點火燃燒是氫氧爆轟燃?xì)鈼l件下發(fā)生的.混合爆轟波結(jié)構(gòu)受到鋁粉在燃?xì)庵悬c火燃燒特性的影響,而存在不同的結(jié)構(gòu).只有顆粒相反應(yīng)釋熱發(fā)生在爆轟波聲速面之前才能對爆轟波的首道激波產(chǎn)生影響,使得爆速增加;而當(dāng)顆粒相的反應(yīng)釋熱落后于聲速面,且在爆轟波的觀測時間內(nèi)顆粒相的點火釋熱能夠發(fā)生,就可以產(chǎn)生第二道激波,得到雙波面的混合爆轟波結(jié)構(gòu).

    試驗結(jié)果表明,鋁粉粒徑的大小對鋁粉點火延遲與燃燒時間影響很大,從而影響混合爆轟波的結(jié)構(gòu).鋁粉粒徑較小時(100 nm 和1 μm),鋁粉化學(xué)反應(yīng)快于聲速面之前發(fā)生,爆轟波由氣相和顆粒燃燒兩部分的反應(yīng)釋熱支持,對比同樣狀態(tài)下的純氣相爆轟,鋁粉的加入使得爆速和波峰壓力增高,此時呈單波爆轟結(jié)構(gòu).鋁粉粒徑較大時則出現(xiàn)了雙波面爆轟的情況.粒徑20 μm 時,第二道波與第一道波速度接近,兩波間距在向下游傳遞過程中基本不變;40 μm時,第二道波波速較慢,因而兩波間距逐漸拉開.粒徑10 μm 的混合爆轟波本質(zhì)上是兩波面距離很近的雙波面結(jié)構(gòu),峰值壓力有大幅提高,但爆速并沒有增加.在鋁粉/氫氣/空氣混合爆轟現(xiàn)象中,反映了鋁粉在復(fù)雜條件下的燃燒特性,通過對混合爆轟試驗結(jié)果的分析,明確了鋁粉的點火燃燒特性對混合爆轟現(xiàn)象的影響機(jī)理.

    致謝

    本文工作的研究思路和試驗方法得到了俞鴻儒院士的建議和幫助,在此表示感謝.

    猜你喜歡
    波面鋁粉聲速
    納米鋁粉的反應(yīng)性研究進(jìn)展及趨勢
    基于恒定陡度聚焦波模型的分析與討論
    水道港口(2020年6期)2020-02-22 11:33:50
    納米鋁粉對RDX基炸藥爆速的影響
    火工品(2019年1期)2019-04-29 03:03:44
    多普勒效應(yīng)中觀察者接收頻率的計算
    淺談光的干涉和衍射的區(qū)別和聯(lián)系
    中文信息(2018年2期)2018-05-30 11:45:10
    聲速是如何測定的
    基于20 L球形爆炸裝置的微米級鋁粉爆炸特性實驗
    防止球形鋁粉導(dǎo)致的安全事故
    山西化工(2015年1期)2015-08-15 00:50:11
    跨聲速風(fēng)洞全模顫振試驗技術(shù)
    機(jī)翼跨聲速抖振研究進(jìn)展
    人体艺术视频欧美日本| 99久久精品国产国产毛片| 99热6这里只有精品| 99热这里只有是精品在线观看| 国产成人精品一,二区| 成人高潮视频无遮挡免费网站| 熟女电影av网| 国产成人91sexporn| 一二三四中文在线观看免费高清| 国产一级毛片在线| 精品一品国产午夜福利视频| 身体一侧抽搐| 中文欧美无线码| 一级a做视频免费观看| 国语对白做爰xxxⅹ性视频网站| 国产精品伦人一区二区| 在线播放无遮挡| www.av在线官网国产| 国产在线免费精品| 欧美精品一区二区免费开放| 亚洲av男天堂| 成人无遮挡网站| 国产爽快片一区二区三区| 不卡视频在线观看欧美| 少妇丰满av| 欧美性感艳星| 91精品国产九色| 色5月婷婷丁香| 亚洲av中文av极速乱| 欧美另类一区| 日本欧美国产在线视频| 欧美3d第一页| 欧美日韩精品成人综合77777| 高清黄色对白视频在线免费看 | 国产日韩欧美在线精品| 欧美 日韩 精品 国产| 久久99热6这里只有精品| 久久精品国产亚洲av天美| 少妇的逼水好多| 精品少妇黑人巨大在线播放| 大话2 男鬼变身卡| av免费观看日本| 成人国产av品久久久| 国产免费一级a男人的天堂| 一级毛片黄色毛片免费观看视频| 亚洲欧美精品自产自拍| 大陆偷拍与自拍| 欧美三级亚洲精品| 免费高清在线观看视频在线观看| 国产精品无大码| 日韩制服骚丝袜av| 国产成人免费观看mmmm| 亚洲欧美精品专区久久| 久久久成人免费电影| 少妇丰满av| 色哟哟·www| 美女福利国产在线 | 青春草视频在线免费观看| 内射极品少妇av片p| 两个人的视频大全免费| 精品国产三级普通话版| 97在线视频观看| 少妇的逼好多水| 久久久久久久亚洲中文字幕| 99热这里只有精品一区| 内地一区二区视频在线| 美女中出高潮动态图| 欧美变态另类bdsm刘玥| 亚洲精华国产精华液的使用体验| 波野结衣二区三区在线| 男女边摸边吃奶| 国产精品一区二区三区四区免费观看| 欧美日韩亚洲高清精品| a级毛色黄片| 国产无遮挡羞羞视频在线观看| 日日啪夜夜爽| 国产大屁股一区二区在线视频| 秋霞伦理黄片| 欧美日本视频| 国产伦在线观看视频一区| 在线播放无遮挡| 国产精品一区二区三区四区免费观看| 国产精品不卡视频一区二区| 免费观看的影片在线观看| 99久久精品国产国产毛片| 人妻系列 视频| 新久久久久国产一级毛片| 久久亚洲国产成人精品v| 精品国产露脸久久av麻豆| 久久精品人妻少妇| 最新中文字幕久久久久| 插阴视频在线观看视频| 一个人免费看片子| 免费大片黄手机在线观看| 永久网站在线| 婷婷色麻豆天堂久久| 最黄视频免费看| 日韩伦理黄色片| 免费观看a级毛片全部| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产av新网站| 亚洲精品,欧美精品| 少妇人妻一区二区三区视频| 国产精品国产三级国产av玫瑰| 亚洲成人手机| 欧美xxxx性猛交bbbb| 成年av动漫网址| 国产白丝娇喘喷水9色精品| 免费播放大片免费观看视频在线观看| videossex国产| 国产亚洲av片在线观看秒播厂| 永久网站在线| 日韩成人av中文字幕在线观看| 特大巨黑吊av在线直播| 日本黄大片高清| 日本wwww免费看| 夜夜看夜夜爽夜夜摸| 一个人免费看片子| 高清在线视频一区二区三区| 精品一品国产午夜福利视频| 国产在线视频一区二区| 久久女婷五月综合色啪小说| 国产av码专区亚洲av| 亚洲av.av天堂| 在线观看一区二区三区| 男女下面进入的视频免费午夜| 狂野欧美激情性xxxx在线观看| 久久久久网色| 亚洲色图综合在线观看| 亚洲精品国产av成人精品| 亚洲人与动物交配视频| 男人添女人高潮全过程视频| 亚洲电影在线观看av| 亚洲国产欧美在线一区| 国产淫片久久久久久久久| 国产日韩欧美亚洲二区| 99热这里只有是精品50| 一个人看视频在线观看www免费| 91午夜精品亚洲一区二区三区| 一级毛片 在线播放| 久久6这里有精品| 性色av一级| 国产av一区二区精品久久 | 激情五月婷婷亚洲| 久久精品国产亚洲av天美| 欧美xxⅹ黑人| 综合色丁香网| av黄色大香蕉| 亚洲成色77777| 乱码一卡2卡4卡精品| 老司机影院成人| 国产成人aa在线观看| 激情 狠狠 欧美| 久久久久精品久久久久真实原创| 人人妻人人看人人澡| 最近最新中文字幕大全电影3| 午夜视频国产福利| 欧美高清成人免费视频www| 国产精品一区二区性色av| 汤姆久久久久久久影院中文字幕| 男人舔奶头视频| 观看av在线不卡| 亚洲内射少妇av| 国产精品嫩草影院av在线观看| 涩涩av久久男人的天堂| 观看美女的网站| 一本一本综合久久| 在线观看免费高清a一片| 最黄视频免费看| 亚洲av福利一区| 色视频在线一区二区三区| 免费观看的影片在线观看| 丝瓜视频免费看黄片| 欧美高清成人免费视频www| 九九在线视频观看精品| 国产真实伦视频高清在线观看| 七月丁香在线播放| 一级二级三级毛片免费看| 亚洲婷婷狠狠爱综合网| 免费在线观看成人毛片| 特大巨黑吊av在线直播| 亚洲人成网站高清观看| 成人午夜精彩视频在线观看| 中文资源天堂在线| 精品国产露脸久久av麻豆| 日本色播在线视频| 激情五月婷婷亚洲| 日韩 亚洲 欧美在线| 亚洲第一av免费看| 久久人妻熟女aⅴ| 国产成人免费观看mmmm| 欧美日韩国产mv在线观看视频 | 黑丝袜美女国产一区| 我的女老师完整版在线观看| 久久ye,这里只有精品| 在线免费观看不下载黄p国产| 亚洲三级黄色毛片| 看免费成人av毛片| 国产免费福利视频在线观看| 少妇人妻一区二区三区视频| 精品久久久久久电影网| 中国国产av一级| 色婷婷久久久亚洲欧美| 国产无遮挡羞羞视频在线观看| 国产日韩欧美亚洲二区| 人人妻人人看人人澡| kizo精华| 夫妻性生交免费视频一级片| 久久国产乱子免费精品| 久久综合国产亚洲精品| 丰满少妇做爰视频| 午夜激情久久久久久久| 大片免费播放器 马上看| 久久热精品热| 最黄视频免费看| 一本—道久久a久久精品蜜桃钙片| 欧美少妇被猛烈插入视频| 嫩草影院入口| 超碰97精品在线观看| 成人午夜精彩视频在线观看| 22中文网久久字幕| 亚洲欧洲日产国产| 狂野欧美激情性bbbbbb| 18+在线观看网站| 国产黄频视频在线观看| 久久国产亚洲av麻豆专区| 国产淫语在线视频| www.色视频.com| 99热全是精品| 两个人的视频大全免费| 插阴视频在线观看视频| 性高湖久久久久久久久免费观看| 一级二级三级毛片免费看| 国产精品麻豆人妻色哟哟久久| 久久国产乱子免费精品| 久久精品国产亚洲网站| 中文天堂在线官网| 亚洲四区av| 黄色配什么色好看| 我要看日韩黄色一级片| 久久婷婷青草| 搡女人真爽免费视频火全软件| 国产中年淑女户外野战色| 国产有黄有色有爽视频| 91精品伊人久久大香线蕉| 亚洲精华国产精华液的使用体验| 亚洲精品色激情综合| 一级av片app| 男女免费视频国产| 久久午夜福利片| 小蜜桃在线观看免费完整版高清| av又黄又爽大尺度在线免费看| 国产欧美日韩一区二区三区在线 | 国产欧美日韩一区二区三区在线 | 久久久欧美国产精品| 色综合色国产| 这个男人来自地球电影免费观看 | 激情 狠狠 欧美| 国产免费一区二区三区四区乱码| 久久国内精品自在自线图片| 97精品久久久久久久久久精品| 又爽又黄a免费视频| 人妻制服诱惑在线中文字幕| 老司机影院毛片| 性高湖久久久久久久久免费观看| 久久久久久伊人网av| 亚洲国产精品成人久久小说| 国产精品不卡视频一区二区| 自拍偷自拍亚洲精品老妇| 高清不卡的av网站| 少妇精品久久久久久久| 日韩人妻高清精品专区| 观看美女的网站| 亚洲精品一区蜜桃| 精品熟女少妇av免费看| 高清av免费在线| 黄色视频在线播放观看不卡| 国产精品99久久久久久久久| av播播在线观看一区| 免费观看无遮挡的男女| 国产精品久久久久久精品古装| 99re6热这里在线精品视频| 观看av在线不卡| 91久久精品国产一区二区成人| 国产v大片淫在线免费观看| 久久久欧美国产精品| 国产精品无大码| 国产精品久久久久久av不卡| av黄色大香蕉| 欧美区成人在线视频| 亚洲精品色激情综合| 下体分泌物呈黄色| 亚洲人成网站高清观看| 国产精品一区二区在线观看99| 中文字幕av成人在线电影| 国产高清不卡午夜福利| 内地一区二区视频在线| 久久久久网色| 欧美变态另类bdsm刘玥| 国产真实伦视频高清在线观看| 亚洲欧美日韩卡通动漫| 国产久久久一区二区三区| 最近中文字幕2019免费版| 啦啦啦在线观看免费高清www| 在线免费观看不下载黄p国产| 成人高潮视频无遮挡免费网站| 色吧在线观看| 午夜激情福利司机影院| 国产永久视频网站| 简卡轻食公司| 欧美精品一区二区大全| 亚洲国产色片| 午夜福利网站1000一区二区三区| 大香蕉97超碰在线| av.在线天堂| 99久久综合免费| 伊人久久精品亚洲午夜| 观看美女的网站| 日韩成人伦理影院| 偷拍熟女少妇极品色| av专区在线播放| 啦啦啦在线观看免费高清www| 免费黄频网站在线观看国产| 日日撸夜夜添| 老女人水多毛片| 久久99热这里只有精品18| 欧美3d第一页| av国产免费在线观看| 干丝袜人妻中文字幕| 免费高清在线观看视频在线观看| 99热这里只有是精品在线观看| 午夜福利高清视频| 国产淫片久久久久久久久| 国产成人精品久久久久久| 插阴视频在线观看视频| 天堂中文最新版在线下载| 久久午夜福利片| 成人免费观看视频高清| 国产成人a∨麻豆精品| 黄色一级大片看看| 六月丁香七月| 熟女电影av网| 男人和女人高潮做爰伦理| 欧美日韩亚洲高清精品| 免费播放大片免费观看视频在线观看| 中国美白少妇内射xxxbb| 一级毛片 在线播放| 51国产日韩欧美| 日本-黄色视频高清免费观看| 噜噜噜噜噜久久久久久91| 乱码一卡2卡4卡精品| 欧美日韩视频精品一区| 国产精品嫩草影院av在线观看| 免费大片黄手机在线观看| 国产精品三级大全| 久久精品国产亚洲av天美| 91aial.com中文字幕在线观看| 在线 av 中文字幕| h日本视频在线播放| 丰满迷人的少妇在线观看| 国产色婷婷99| 婷婷色综合大香蕉| 亚洲人与动物交配视频| 亚洲国产高清在线一区二区三| 高清不卡的av网站| h视频一区二区三区| 亚洲怡红院男人天堂| 久久久久性生活片| 大陆偷拍与自拍| 搡老乐熟女国产| 男人狂女人下面高潮的视频| 国产成人精品福利久久| 亚洲美女视频黄频| 99久久综合免费| 免费黄频网站在线观看国产| 1000部很黄的大片| 最近中文字幕高清免费大全6| 国产91av在线免费观看| 日本-黄色视频高清免费观看| 国产亚洲最大av| 国产一级毛片在线| 夜夜看夜夜爽夜夜摸| 噜噜噜噜噜久久久久久91| 亚洲美女搞黄在线观看| 亚洲av国产av综合av卡| 一级毛片黄色毛片免费观看视频| 国产女主播在线喷水免费视频网站| 日韩伦理黄色片| av女优亚洲男人天堂| 国产成人精品福利久久| 亚洲精品国产成人久久av| 欧美成人午夜免费资源| 欧美精品一区二区大全| 少妇猛男粗大的猛烈进出视频| 夜夜爽夜夜爽视频| 午夜福利在线在线| 国模一区二区三区四区视频| 久久久精品免费免费高清| 精品久久久久久电影网| 免费在线观看成人毛片| 久久久久久久久久人人人人人人| 丰满乱子伦码专区| 精品亚洲乱码少妇综合久久| 男女边摸边吃奶| 久久午夜福利片| 日韩欧美一区视频在线观看 | .国产精品久久| 色哟哟·www| 亚洲图色成人| 日韩欧美 国产精品| 久久综合国产亚洲精品| 99久久人妻综合| 日韩av不卡免费在线播放| 国产黄片视频在线免费观看| 免费黄频网站在线观看国产| freevideosex欧美| 最后的刺客免费高清国语| 人人妻人人爽人人添夜夜欢视频 | 纵有疾风起免费观看全集完整版| 美女内射精品一级片tv| 久久av网站| 日本猛色少妇xxxxx猛交久久| 在线观看人妻少妇| 视频区图区小说| 欧美日韩亚洲高清精品| 久久久久久久久大av| 蜜桃久久精品国产亚洲av| 只有这里有精品99| 丝袜脚勾引网站| 一级av片app| 五月开心婷婷网| 国产精品欧美亚洲77777| 久久精品国产亚洲av天美| 日本与韩国留学比较| 欧美成人一区二区免费高清观看| 午夜福利影视在线免费观看| 国产伦理片在线播放av一区| 久久久国产一区二区| 欧美成人午夜免费资源| 亚洲欧美日韩无卡精品| 97精品久久久久久久久久精品| 国产成人精品久久久久久| 欧美成人a在线观看| 欧美激情极品国产一区二区三区 | 国产免费一区二区三区四区乱码| 国产男女超爽视频在线观看| 久久久久精品性色| 国产 一区精品| 精品亚洲乱码少妇综合久久| 国产精品免费大片| 十八禁网站网址无遮挡 | 乱码一卡2卡4卡精品| 一边亲一边摸免费视频| 熟女av电影| 欧美zozozo另类| 欧美日韩一区二区视频在线观看视频在线| 最近的中文字幕免费完整| 亚洲熟女精品中文字幕| 老女人水多毛片| 黑人高潮一二区| 久久精品国产鲁丝片午夜精品| 国产精品一区二区性色av| 日本色播在线视频| 乱系列少妇在线播放| 最近中文字幕2019免费版| 免费观看性生交大片5| 久热这里只有精品99| 黑人高潮一二区| 国产久久久一区二区三区| 黄片无遮挡物在线观看| 国产成人aa在线观看| 在线看a的网站| 丰满人妻一区二区三区视频av| 欧美精品一区二区免费开放| 成人二区视频| 美女国产视频在线观看| 美女中出高潮动态图| 免费大片18禁| 精品酒店卫生间| 国精品久久久久久国模美| 蜜桃久久精品国产亚洲av| 亚洲,欧美,日韩| 国产男人的电影天堂91| av黄色大香蕉| 日韩制服骚丝袜av| 欧美 日韩 精品 国产| 亚洲欧美一区二区三区黑人 | 赤兔流量卡办理| 国产 精品1| 日本-黄色视频高清免费观看| 99九九线精品视频在线观看视频| 男人添女人高潮全过程视频| 国内揄拍国产精品人妻在线| 国产精品一二三区在线看| 全区人妻精品视频| 少妇猛男粗大的猛烈进出视频| 国产成人午夜福利电影在线观看| 国产精品国产三级国产av玫瑰| 亚洲精品456在线播放app| 只有这里有精品99| 一本一本综合久久| 在线天堂最新版资源| 午夜福利视频精品| 性高湖久久久久久久久免费观看| 国精品久久久久久国模美| 伊人久久精品亚洲午夜| 日韩成人av中文字幕在线观看| 熟女av电影| 久久综合国产亚洲精品| 国产精品99久久久久久久久| 一本久久精品| 久久久午夜欧美精品| 直男gayav资源| 久久久精品94久久精品| 久久女婷五月综合色啪小说| 久久久久国产精品人妻一区二区| 国产精品av视频在线免费观看| 午夜免费观看性视频| 免费观看无遮挡的男女| 国产永久视频网站| 亚洲不卡免费看| 国产精品免费大片| 亚洲色图av天堂| 国产美女午夜福利| 日韩制服骚丝袜av| 亚洲精品一二三| av卡一久久| 亚洲精品乱码久久久久久按摩| 成年av动漫网址| 97在线视频观看| 人人妻人人澡人人爽人人夜夜| 免费观看的影片在线观看| 91aial.com中文字幕在线观看| 嘟嘟电影网在线观看| 亚洲欧洲国产日韩| 中文字幕人妻熟人妻熟丝袜美| 久久久久久伊人网av| 18禁动态无遮挡网站| 国产精品偷伦视频观看了| 熟女人妻精品中文字幕| 身体一侧抽搐| 亚洲欧美清纯卡通| 人妻一区二区av| 国产精品女同一区二区软件| 久久久欧美国产精品| 高清av免费在线| 亚洲国产毛片av蜜桃av| 免费不卡的大黄色大毛片视频在线观看| 欧美高清性xxxxhd video| 亚洲精品日韩av片在线观看| 成年免费大片在线观看| 只有这里有精品99| 婷婷色av中文字幕| 欧美高清成人免费视频www| 国内揄拍国产精品人妻在线| 新久久久久国产一级毛片| 777米奇影视久久| 另类亚洲欧美激情| 视频中文字幕在线观看| 色5月婷婷丁香| 午夜免费鲁丝| 国产成人freesex在线| 最近中文字幕高清免费大全6| 亚洲成色77777| 我要看黄色一级片免费的| 男的添女的下面高潮视频| 青春草国产在线视频| 只有这里有精品99| 如何舔出高潮| 久久久成人免费电影| 在线观看人妻少妇| 国产亚洲欧美精品永久| 国产在线男女| 观看美女的网站| 免费大片黄手机在线观看| av国产久精品久网站免费入址| 久久人妻熟女aⅴ| 国产黄色视频一区二区在线观看| 自拍偷自拍亚洲精品老妇| 中文字幕制服av| 最近最新中文字幕免费大全7| 青春草视频在线免费观看| 国产亚洲最大av| 女人久久www免费人成看片| 亚洲精品成人av观看孕妇| 高清日韩中文字幕在线| 国产成人freesex在线| 亚洲人与动物交配视频| 亚洲一级一片aⅴ在线观看| 亚洲国产日韩一区二区| 黄色视频在线播放观看不卡| 国产精品一及| 国产欧美亚洲国产| 一级黄片播放器| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美日韩东京热| 少妇 在线观看| 亚洲精品色激情综合| 亚洲色图av天堂| 久久精品久久久久久久性| 九色成人免费人妻av| 欧美日韩综合久久久久久| 婷婷色av中文字幕| 免费av中文字幕在线| 国产有黄有色有爽视频| 99热6这里只有精品| 国产精品无大码| 亚洲国产精品国产精品| 日韩av免费高清视频| 高清视频免费观看一区二区| 一区二区三区四区激情视频| 久久久久久久久久成人| 国产成人免费无遮挡视频| 99精国产麻豆久久婷婷| 黄片wwwwww|