• <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)展
    欧美高清成人免费视频www| 男女国产视频网站| 免费观看在线日韩| 久久精品国产a三级三级三级| 久久国产亚洲av麻豆专区| 国产极品天堂在线| 91精品一卡2卡3卡4卡| 超碰97精品在线观看| 一区二区三区四区激情视频| 国内揄拍国产精品人妻在线| 色94色欧美一区二区| tube8黄色片| 亚洲av免费高清在线观看| 日韩av免费高清视频| 婷婷色麻豆天堂久久| 不卡视频在线观看欧美| 99热这里只有精品一区| 精品少妇内射三级| 丰满人妻一区二区三区视频av| 日本-黄色视频高清免费观看| 男女无遮挡免费网站观看| 少妇熟女欧美另类| 亚洲欧美日韩卡通动漫| 菩萨蛮人人尽说江南好唐韦庄| 免费观看在线日韩| 国产午夜精品一二区理论片| 久久精品久久久久久久性| 久久久久久久国产电影| 国产 精品1| 97精品久久久久久久久久精品| 日韩一区二区三区影片| 少妇被粗大猛烈的视频| 99热国产这里只有精品6| 日韩三级伦理在线观看| 在线观看国产h片| 国产成人a∨麻豆精品| 亚洲精品日本国产第一区| 少妇 在线观看| 成年av动漫网址| 黄色欧美视频在线观看| www.色视频.com| 在线观看免费日韩欧美大片 | 日本wwww免费看| 少妇人妻 视频| 一边亲一边摸免费视频| 人人妻人人添人人爽欧美一区卜| 一本久久精品| 97精品久久久久久久久久精品| 成年av动漫网址| 亚洲美女黄色视频免费看| 久久免费观看电影| 亚洲怡红院男人天堂| 在线观看三级黄色| 国产黄片视频在线免费观看| 久久国内精品自在自线图片| 人妻系列 视频| 高清在线视频一区二区三区| 亚洲人成网站在线观看播放| 国产成人精品福利久久| 欧美精品一区二区大全| 日本午夜av视频| 国产成人精品婷婷| 亚洲精品久久午夜乱码| 色婷婷久久久亚洲欧美| 美女视频免费永久观看网站| 亚洲第一av免费看| 内射极品少妇av片p| 国产精品.久久久| 国产av精品麻豆| 久久免费观看电影| 日本色播在线视频| 日韩中文字幕视频在线看片| 免费av中文字幕在线| 99九九线精品视频在线观看视频| 精品久久久精品久久久| 性色av一级| 午夜影院在线不卡| 亚洲国产日韩一区二区| 久久久久久久久久久久大奶| 亚洲精品日韩av片在线观看| 一边亲一边摸免费视频| 精品少妇久久久久久888优播| 久久久久精品性色| 18禁动态无遮挡网站| 又爽又黄a免费视频| 美女中出高潮动态图| 国产伦精品一区二区三区四那| 国产欧美日韩精品一区二区| 18+在线观看网站| 亚洲欧美成人精品一区二区| h视频一区二区三区| 国产欧美亚洲国产| 黄色日韩在线| 久久久午夜欧美精品| 菩萨蛮人人尽说江南好唐韦庄| 我要看日韩黄色一级片| 观看美女的网站| 亚洲va在线va天堂va国产| 中文精品一卡2卡3卡4更新| 久久av网站| 伊人久久精品亚洲午夜| 久久亚洲国产成人精品v| 国产精品国产三级国产专区5o| 免费观看a级毛片全部| 寂寞人妻少妇视频99o| 少妇高潮的动态图| av视频免费观看在线观看| 新久久久久国产一级毛片| 国产伦精品一区二区三区四那| 黄色日韩在线| 在线观看免费日韩欧美大片 | 国产精品福利在线免费观看| 免费不卡的大黄色大毛片视频在线观看| 街头女战士在线观看网站| 欧美日韩精品成人综合77777| 日韩伦理黄色片| 国产黄色免费在线视频| 中文字幕人妻丝袜制服| 亚洲中文av在线| 人妻 亚洲 视频| 伦精品一区二区三区| 色5月婷婷丁香| 在线观看免费视频网站a站| 伊人久久国产一区二区| 国产精品久久久久久久久免| 国产精品一区www在线观看| 人妻人人澡人人爽人人| 91在线精品国自产拍蜜月| 国产精品一区二区在线不卡| 欧美精品国产亚洲| 精品久久久久久久久亚洲| 男人舔奶头视频| 国产女主播在线喷水免费视频网站| 看免费成人av毛片| 免费看不卡的av| 如日韩欧美国产精品一区二区三区 | 成人亚洲欧美一区二区av| 国产成人精品福利久久| 国产黄色免费在线视频| av又黄又爽大尺度在线免费看| 少妇裸体淫交视频免费看高清| 18禁在线无遮挡免费观看视频| 韩国av在线不卡| 欧美97在线视频| av福利片在线观看| 国产午夜精品久久久久久一区二区三区| 欧美97在线视频| 波野结衣二区三区在线| 国产 一区精品| 免费av中文字幕在线| av不卡在线播放| 中文字幕制服av| 少妇的逼水好多| 日韩,欧美,国产一区二区三区| 日韩人妻高清精品专区| 亚洲自偷自拍三级| videossex国产| 在线亚洲精品国产二区图片欧美 | 99国产精品免费福利视频| 久久久亚洲精品成人影院| 国产欧美另类精品又又久久亚洲欧美| 毛片一级片免费看久久久久| 极品教师在线视频| 尾随美女入室| 伊人亚洲综合成人网| 国产成人免费观看mmmm| 久久影院123| 国产精品一区www在线观看| 99热这里只有是精品50| 99精国产麻豆久久婷婷| 久久久久精品久久久久真实原创| 天堂8中文在线网| 免费在线观看成人毛片| 国产av一区二区精品久久| 大话2 男鬼变身卡| 99久久中文字幕三级久久日本| 亚洲婷婷狠狠爱综合网| 国产爽快片一区二区三区| 精品人妻熟女毛片av久久网站| 一区二区三区精品91| 狂野欧美激情性xxxx在线观看| av女优亚洲男人天堂| 国产色婷婷99| 国产免费福利视频在线观看| 精品久久久久久电影网| 男女啪啪激烈高潮av片| 国产伦理片在线播放av一区| 岛国毛片在线播放| 在线精品无人区一区二区三| 国产69精品久久久久777片| 久久精品国产亚洲av天美| 色5月婷婷丁香| 久久热精品热| 亚洲内射少妇av| 国产日韩欧美视频二区| 高清视频免费观看一区二区| 国产一区二区在线观看日韩| 六月丁香七月| 亚洲国产成人一精品久久久| av有码第一页| 午夜av观看不卡| 精品国产一区二区久久| 国产淫语在线视频| 成年人免费黄色播放视频 | 国产一区二区三区综合在线观看 | 精品人妻偷拍中文字幕| 亚洲av福利一区| av女优亚洲男人天堂| 日韩成人伦理影院| 日韩伦理黄色片| 夜夜骑夜夜射夜夜干| 午夜免费观看性视频| 自拍偷自拍亚洲精品老妇| 高清av免费在线| 欧美日韩国产mv在线观看视频| 又黄又爽又刺激的免费视频.| 天天操日日干夜夜撸| 97在线视频观看| 国产精品久久久久久久久免| 美女脱内裤让男人舔精品视频| 日本猛色少妇xxxxx猛交久久| 在线观看www视频免费| 日韩熟女老妇一区二区性免费视频| 日韩大片免费观看网站| av女优亚洲男人天堂| 亚洲国产精品999| 精品国产乱码久久久久久小说| 日韩免费高清中文字幕av| 秋霞伦理黄片| 国产精品偷伦视频观看了| 91精品伊人久久大香线蕉| 欧美 日韩 精品 国产| 黄色日韩在线| 国产亚洲欧美精品永久| 久久99热6这里只有精品| 亚洲中文av在线| 日韩三级伦理在线观看| 日韩 亚洲 欧美在线| 国产伦理片在线播放av一区| 日韩不卡一区二区三区视频在线| 国产精品熟女久久久久浪| 国产精品不卡视频一区二区| 精品国产乱码久久久久久小说| 欧美 亚洲 国产 日韩一| 国产精品人妻久久久久久| 又大又黄又爽视频免费| 男人和女人高潮做爰伦理| 日韩视频在线欧美| 男女免费视频国产| 搡老乐熟女国产| 免费黄网站久久成人精品| 国产免费一级a男人的天堂| 久久人人爽av亚洲精品天堂| 老熟女久久久| 天天躁夜夜躁狠狠久久av| 日本爱情动作片www.在线观看| av在线app专区| freevideosex欧美| 久久韩国三级中文字幕| 亚洲第一av免费看| 自拍偷自拍亚洲精品老妇| 欧美日韩国产mv在线观看视频| 国产伦精品一区二区三区四那| 精华霜和精华液先用哪个| 国产精品久久久久成人av| 成年女人在线观看亚洲视频| 国产精品国产三级专区第一集| av卡一久久| 国产欧美日韩精品一区二区| 国产午夜精品久久久久久一区二区三区| 嫩草影院入口| 欧美一级a爱片免费观看看| 日韩av不卡免费在线播放| 丰满迷人的少妇在线观看| 欧美成人精品欧美一级黄| 午夜老司机福利剧场| 中文在线观看免费www的网站| 三级国产精品欧美在线观看| 激情五月婷婷亚洲| 国产无遮挡羞羞视频在线观看| 国产成人精品福利久久| 国产一区二区在线观看av| av女优亚洲男人天堂| 国产午夜精品一二区理论片| 亚洲精品日韩av片在线观看| av国产久精品久网站免费入址| 国产精品麻豆人妻色哟哟久久| 一级,二级,三级黄色视频| 18禁在线播放成人免费| 老熟女久久久| 两个人的视频大全免费| 久久影院123| 自线自在国产av| 高清黄色对白视频在线免费看 | tube8黄色片| 日韩一本色道免费dvd| 国产精品女同一区二区软件| 高清毛片免费看| 亚洲电影在线观看av| 国国产精品蜜臀av免费| 如何舔出高潮| 欧美最新免费一区二区三区| 一本色道久久久久久精品综合| 最近的中文字幕免费完整| 亚洲在久久综合| 亚洲国产av新网站| 男人爽女人下面视频在线观看| 色5月婷婷丁香| 九九久久精品国产亚洲av麻豆| 99久国产av精品国产电影| a级片在线免费高清观看视频| 美女主播在线视频| 大香蕉97超碰在线| 久热这里只有精品99| 五月伊人婷婷丁香| 春色校园在线视频观看| 久久久精品免费免费高清| 国产一区二区在线观看日韩| av不卡在线播放| 日韩成人伦理影院| 日本vs欧美在线观看视频 | 高清黄色对白视频在线免费看 | 国产探花极品一区二区| 十分钟在线观看高清视频www | 日韩熟女老妇一区二区性免费视频| 日韩三级伦理在线观看| 2021少妇久久久久久久久久久| 99久久精品一区二区三区| 黄片无遮挡物在线观看| 啦啦啦视频在线资源免费观看| av线在线观看网站| 久久99一区二区三区| 国产精品人妻久久久影院| 五月天丁香电影| 蜜桃在线观看..| 欧美bdsm另类| 亚洲欧美一区二区三区国产| 少妇高潮的动态图| av又黄又爽大尺度在线免费看| 超碰97精品在线观看| 成年美女黄网站色视频大全免费 | 91精品伊人久久大香线蕉| 成年女人在线观看亚洲视频| 精品国产乱码久久久久久小说| 国产淫片久久久久久久久| 麻豆成人av视频| 最近最新中文字幕免费大全7| 少妇熟女欧美另类| 欧美区成人在线视频| 搡老乐熟女国产| 免费黄网站久久成人精品| 亚洲激情五月婷婷啪啪| 亚洲国产最新在线播放| 精品人妻熟女毛片av久久网站| 国产亚洲欧美精品永久| 午夜激情福利司机影院| 成人免费观看视频高清| 国产亚洲最大av| 亚洲国产毛片av蜜桃av| av免费观看日本| 人妻少妇偷人精品九色| 精品人妻熟女毛片av久久网站| 777米奇影视久久| 亚洲成人一二三区av| 五月开心婷婷网| 黑人猛操日本美女一级片| 亚洲精品国产av蜜桃| 我的老师免费观看完整版| 91在线精品国自产拍蜜月| 免费人成在线观看视频色| 亚洲真实伦在线观看| 18禁在线无遮挡免费观看视频| 熟女av电影| 爱豆传媒免费全集在线观看| 91精品国产九色| 91午夜精品亚洲一区二区三区| 亚洲va在线va天堂va国产| 久久久久久久久久成人| 国产精品免费大片| 少妇熟女欧美另类| 亚洲,欧美,日韩| 国产精品国产三级国产av玫瑰| 91精品国产九色| 久久久亚洲精品成人影院| av国产久精品久网站免费入址| 欧美+日韩+精品| 自拍欧美九色日韩亚洲蝌蚪91 | 一区在线观看完整版| 一区二区三区乱码不卡18| 日韩精品免费视频一区二区三区 | 亚洲四区av| 国产成人freesex在线| 各种免费的搞黄视频| 人妻少妇偷人精品九色| 三级经典国产精品| 男女国产视频网站| 高清毛片免费看| 午夜福利网站1000一区二区三区| 人人妻人人添人人爽欧美一区卜| 一级爰片在线观看| 26uuu在线亚洲综合色| 人妻系列 视频| 人妻夜夜爽99麻豆av| 日日爽夜夜爽网站| 国产精品国产av在线观看| 亚洲av男天堂| 美女视频免费永久观看网站| 天堂中文最新版在线下载| 国产黄片视频在线免费观看| 在线观看免费视频网站a站| 国产精品秋霞免费鲁丝片| 曰老女人黄片| 成人亚洲欧美一区二区av| 一级毛片久久久久久久久女| 夜夜爽夜夜爽视频| 3wmmmm亚洲av在线观看| 简卡轻食公司| 看免费成人av毛片| 97在线视频观看| 国产色爽女视频免费观看| 三级经典国产精品| 国产一区亚洲一区在线观看| 纯流量卡能插随身wifi吗| 亚洲av中文av极速乱| 亚洲av免费高清在线观看| 国产白丝娇喘喷水9色精品| 国产亚洲最大av| 99热这里只有是精品在线观看| 亚洲国产精品999| 日本午夜av视频| 久久韩国三级中文字幕| 精品人妻一区二区三区麻豆| 午夜福利在线观看免费完整高清在| av一本久久久久| 国产一级毛片在线| 国产高清不卡午夜福利| 久久人人爽人人爽人人片va| 丰满饥渴人妻一区二区三| 国产白丝娇喘喷水9色精品| av又黄又爽大尺度在线免费看| 日韩一区二区三区影片| 国产白丝娇喘喷水9色精品| 大片电影免费在线观看免费| 丰满饥渴人妻一区二区三| h日本视频在线播放| 中文精品一卡2卡3卡4更新| 99热这里只有是精品在线观看| www.av在线官网国产| 亚洲一级一片aⅴ在线观看| 国产精品久久久久久精品电影小说| 在线播放无遮挡| 亚洲国产av新网站| 女性被躁到高潮视频| 麻豆成人av视频| 七月丁香在线播放| 少妇熟女欧美另类| 亚洲成人av在线免费| 一级毛片aaaaaa免费看小| 久久99精品国语久久久| 丰满少妇做爰视频| 国产午夜精品一二区理论片| 人人妻人人澡人人看| 欧美精品亚洲一区二区| 人妻制服诱惑在线中文字幕| 亚州av有码| 妹子高潮喷水视频| 全区人妻精品视频| 久久人人爽av亚洲精品天堂| 国产亚洲av片在线观看秒播厂| 国产极品天堂在线| videossex国产| 亚洲精华国产精华液的使用体验| 亚洲av成人精品一二三区| 日韩 亚洲 欧美在线| 成人亚洲精品一区在线观看| 国产精品欧美亚洲77777| 18禁动态无遮挡网站| 永久网站在线| 啦啦啦中文免费视频观看日本| 国产日韩欧美亚洲二区| 国产永久视频网站| 美女福利国产在线| 亚洲丝袜综合中文字幕| 国产一区二区三区av在线| 婷婷色综合www| 亚洲精品一二三| 精品久久久久久电影网| 欧美3d第一页| 国产成人精品无人区| 一区二区三区四区激情视频| 免费人成在线观看视频色| 少妇丰满av| 狂野欧美白嫩少妇大欣赏| 三级经典国产精品| 精品久久久久久久久av| 国产无遮挡羞羞视频在线观看| 国产av码专区亚洲av| 亚洲国产精品一区三区| 99热这里只有是精品50| 精华霜和精华液先用哪个| 老司机影院成人| 少妇高潮的动态图| 亚洲精品国产成人久久av| 国产亚洲av片在线观看秒播厂| 亚洲天堂av无毛| 午夜激情福利司机影院| 在线天堂最新版资源| 久久精品久久精品一区二区三区| 日韩欧美一区视频在线观看 | 精品一区二区三区视频在线| 国产精品.久久久| 在线观看国产h片| 观看免费一级毛片| 最新中文字幕久久久久| 在线免费观看不下载黄p国产| 人体艺术视频欧美日本| 亚洲精品久久午夜乱码| 日韩欧美一区视频在线观看 | 欧美精品亚洲一区二区| 免费大片黄手机在线观看| xxx大片免费视频| 国产精品伦人一区二区| 成年美女黄网站色视频大全免费 | 青春草亚洲视频在线观看| 亚洲精品一区蜜桃| 国产熟女午夜一区二区三区 | 亚洲四区av| .国产精品久久| 午夜免费鲁丝| 国产探花极品一区二区| 9色porny在线观看| 日韩成人伦理影院| 一级毛片久久久久久久久女| 国产精品国产av在线观看| 国产精品国产三级专区第一集| 亚洲第一区二区三区不卡| 一区二区三区乱码不卡18| 一本—道久久a久久精品蜜桃钙片| 麻豆乱淫一区二区| 97精品久久久久久久久久精品| 久久久国产精品麻豆| 日本爱情动作片www.在线观看| 亚洲精华国产精华液的使用体验| 久久久国产一区二区| 亚洲国产毛片av蜜桃av| 国产伦精品一区二区三区视频9| 国产一级毛片在线| 久久久午夜欧美精品| 熟女人妻精品中文字幕| 人妻夜夜爽99麻豆av| 夜夜骑夜夜射夜夜干| av天堂中文字幕网| 91精品国产国语对白视频| 伦精品一区二区三区| 久久午夜综合久久蜜桃| 午夜视频国产福利| 日韩av免费高清视频| 欧美精品高潮呻吟av久久| 99久久中文字幕三级久久日本| av线在线观看网站| 少妇人妻 视频| 精品99又大又爽又粗少妇毛片| 国产精品无大码| 久久精品国产亚洲网站| 午夜福利网站1000一区二区三区| 91在线精品国自产拍蜜月| 国产成人一区二区在线| 成人影院久久| av视频免费观看在线观看| 国产精品三级大全| 成年女人在线观看亚洲视频| 大码成人一级视频| 国产一区有黄有色的免费视频| 女人精品久久久久毛片| 极品少妇高潮喷水抽搐| www.色视频.com| 日日啪夜夜爽| 99久久精品热视频| 国产乱人偷精品视频| 国产精品秋霞免费鲁丝片| 亚洲不卡免费看| 亚洲怡红院男人天堂| 乱人伦中国视频| 九九久久精品国产亚洲av麻豆| 内地一区二区视频在线| 久久精品久久久久久噜噜老黄| 晚上一个人看的免费电影| 五月伊人婷婷丁香| 欧美 日韩 精品 国产| 69精品国产乱码久久久| 香蕉精品网在线| 3wmmmm亚洲av在线观看| 99热这里只有精品一区| 亚洲怡红院男人天堂| 新久久久久国产一级毛片| www.色视频.com| 又爽又黄a免费视频| 色吧在线观看| 少妇精品久久久久久久| 一级毛片 在线播放| 国产亚洲av片在线观看秒播厂| 99久久中文字幕三级久久日本| 丝袜喷水一区| xxx大片免费视频| 男人添女人高潮全过程视频| 日韩av免费高清视频| h视频一区二区三区| a级毛片免费高清观看在线播放| 少妇的逼水好多| 国产精品一区二区在线不卡| 久久亚洲国产成人精品v| 偷拍熟女少妇极品色| 久久精品国产亚洲av涩爱| 国精品久久久久久国模美| 成人亚洲精品一区在线观看|