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

    橫向聲波擾動下的乙醇燃燒火焰結(jié)構(gòu)和振蕩特性

    2022-07-12 01:13:18張玉濤林國鋮史學(xué)強張玉杰王興明
    工程科學(xué)學(xué)報 2022年8期
    關(guān)鍵詞:油池渦旋頻域

    張玉濤,林國鋮,史學(xué)強,張玉杰,王興明

    西安科技大學(xué)安全科學(xué)與工程學(xué)院,西安 710054

    如何實現(xiàn)高效清潔的滅火一直是火災(zāi)科學(xué)領(lǐng)域的熱點問題[1].相對于傳統(tǒng)常見的滅火技術(shù),無污染的聲波滅火技術(shù)近年來逐漸進(jìn)入學(xué)者們的研究視野[2-3].2002年Baillot和Demare[4]研究了2 ~3000 Hz范圍內(nèi)的聲波對于擴散射流火焰的影響,結(jié)果表明小于50 Hz的聲波會顯著降低火焰的穩(wěn)定性;此外,聲波對火焰的抑制行為已被證明在微重力狀態(tài)下有顯著的效果,未來可將其推廣應(yīng)用于滅火易造成二次損失的封閉或半封閉場所[5].

    關(guān)于聲波對火焰行為影響的研究最早源自于內(nèi)燃機和鍋爐燃燒引發(fā)的熱聲不穩(wěn)定現(xiàn)象[6-9],學(xué)者們從改善燃燒效率的角度出發(fā),研究封閉狀態(tài)下的聲波對火焰行為的控制效果:Hakim等[6]和Han等[7]發(fā)現(xiàn)燃燒噪聲的不同聲學(xué)模式會導(dǎo)致火焰不穩(wěn)定振蕩,造成局部熄火的現(xiàn)象;Vignat等[8]針對甲烷燃燒器內(nèi)火焰對聲擾動的響應(yīng)進(jìn)行了研究,發(fā)現(xiàn)當(dāng)波長與火焰厚度相當(dāng)時,火焰對聲學(xué)擾動的敏感性最大,火焰形態(tài)和燃燒的延續(xù)都會受到影響;Kypraiou等[9]利用平面激光誘導(dǎo)熒光技術(shù)(Planer laser induced fluorescence,PLIF)對非預(yù)混火焰在振蕩聲強迫環(huán)境下的熄滅條件進(jìn)行了研究.

    在此基礎(chǔ)上,近年來已有部分研究從針對封閉環(huán)境的聲振蕩影響衍伸至開放環(huán)境聲波對火焰的熄滅和抑制,進(jìn)而研究聲波在開放環(huán)境下的滅火應(yīng)用:Xiong等[10]研究了聲波對于下落液滴火焰的熄滅效果,并將Damkohler(Da數(shù))作為判定熄滅的標(biāo)準(zhǔn)分析下落液滴火焰的熄滅機制;Zong等[11]實驗測試了橫向聲強迫對射流擴散火焰的抑制效果,發(fā)現(xiàn)橫向低頻聲波能明顯降低擴散火焰的火焰高度.Friedman和Stoliarov[12]對線形火焰在聲波擾動作用下的燃燒行為進(jìn)行了研究;此外,魏珠萍[13]對蠟燭以及香棒在聲波作用下的燃燒特性的研究表明,不同聲壓級范圍的聲波對固體擴散火焰的燃燒速率具有促進(jìn)和抑制的雙重作用;Fujisawa等[14]通過紋影技術(shù)分析了聲波作用下火焰形態(tài)的轉(zhuǎn)變過程;Beisner等[5]的研究表明,在微重力狀態(tài)下聲波對火焰的有良好的抑制效果.這些研究為聲波滅火的實際運用和機理研究提供了一定的借鑒.

    就現(xiàn)有研究而言,聲波對于實驗室級別的小尺度火焰的熄滅、控制效果已經(jīng)得到文獻(xiàn)證實,但對滅火和抑制火焰的具體尺度范圍還未有準(zhǔn)確的定義,火焰的具體抑制原理也尚不明晰.同時,常見的固體、液體和氣體擴散火焰由于其不同的燃燒特點,在聲波強迫作用下的發(fā)展特點也不盡相同[15].因此,本文針對油池火焰,研究了3、4和5 cm直徑的油池火焰在30~90 Hz頻率聲波作用下的幾何形態(tài)、頻域分布和火焰特征參數(shù).為聲強迫抑制火焰行為的具體條件及聲波滅火的應(yīng)用提供了參考.

    1 實驗裝置

    實驗設(shè)備包括信號發(fā)生器、功放、揚聲器、測聲裝置、油池及高速攝像機,實驗臺布置如圖1(a)所示.油池為不銹鋼杯,高度均為8 cm,直徑分別為3、4和5 cm,燃料為無水乙醇,點燃時液面高度距離油池杯口0.5 cm,環(huán)境溫度為28 ℃.揚聲器直徑為30 cm,放置位置距離火焰15 cm,中心正對不銹鋼杯口.聲波信號為正弦信號,實驗頻率范圍的選取參考已有研究的頻率范圍[10-14]并結(jié)合預(yù)實驗,結(jié)果發(fā)現(xiàn)火焰在100 ~500 Hz頻率的聲波作用下呈現(xiàn)與自然燃燒近似一致的火焰行為,因此選取了 30 ~90 Hz(間隔為 10 Hz)的聲波作為實驗的頻率范圍.由于相同的揚聲器聲功率在播放不同聲頻率時聲壓級會存在較大差異,實驗通過檢測并調(diào)節(jié)聲壓級大小的方式使不同聲頻率下的聲壓相近.測聲裝置采用多通道動態(tài)信號分析系統(tǒng)(北京聲望聲電技術(shù)公司),聲探頭布置在油池上方10 cm處,聲信號測試在每次點燃火焰前進(jìn)行,各頻率下的聲壓級如圖1(b)所示.高速攝像機型號為臺灣衡欣9501,采樣頻率設(shè)定為每秒200幀,分辨率為800像素×600像素,設(shè)置在距油池1.5 m處且垂直于油池與揚聲器連線,攝像頭正對火焰.根據(jù)奈奎斯特定理[16],高速攝像機的幀數(shù)及測聲裝置的采樣頻率均大于聲信號最大頻率的2倍,以保證能捕捉火焰與聲波耦合作用下的整體信號.

    圖1 實驗裝置示意圖及聲波參數(shù)圖.(a)裝置示意圖;(b)聲參數(shù)圖Fig.1 Diagram of experimental system and sound parameters: (a) experimental system; (b) sound parameters

    2 結(jié)果與分析

    2.1 火焰形態(tài)

    油池火由于重力作用,豎直方向上存在壓力梯度,而在水平方向上由于燃燒產(chǎn)生的熱氣體和環(huán)境氣體的密度差異而產(chǎn)生密度梯度,兩者梯度的不重合導(dǎo)致斜壓渦量產(chǎn)生,進(jìn)而形成渦旋[17].聲波擾動會造成火焰在徑向方向的密度梯度和壓力梯度的不穩(wěn)定變化,引起火焰面的彎曲變形和火焰渦旋的動態(tài)演變,而不同尺度的油池火焰也會呈現(xiàn)不同的形態(tài)特征.通過高速攝像機拍攝聲波作用下的火焰圖像,結(jié)果如圖2所示.由于火焰形狀呈現(xiàn)周期性交替改變,選取每種頻率下出現(xiàn)頻次最高的火焰形態(tài)發(fā)展圖像序列進(jìn)行說明.高速攝像機對每組實驗拍攝時長為7 ~ 10 s,每張圖像的時間間隔為0.005 s.

    圖2 無聲波強迫下與聲波強迫下不同直徑油池火焰部分火焰圖像(左、中、右分別為3,4和5 cm直徑油池火焰;a,b,c,d,e序列分別為無聲波,30,50,70 和 90 Hz聲波)Fig.2 Partial flame images of pool flame under an acoustic force (left, middle, and right pool flames are 3, 4, and 5 cm in diameter, respectively; a, b, c,d, and e sequences are nonacoustic with 30, 50, 70, and 90 Hz acoustic force, respectively)

    3 cm直徑的乙醇油池火焰在部分頻率下的火焰連續(xù)圖像如圖2(左)所示,無聲波激勵下(a序列)火焰呈現(xiàn)規(guī)律脈動的狀態(tài),火焰兩側(cè)根部出現(xiàn)軸對稱側(cè)渦旋結(jié)構(gòu),逐漸上升至火焰頂部,在火焰末端處脫落,此時渦旋的上升周期時間相對長,這種行為也被稱為“膨脹”(puffing)[18];在30 Hz 和50 Hz聲波作用下,火焰向水平方向發(fā)生偏轉(zhuǎn)(b序列和c序列).這是由于側(cè)渦旋的運動使得火焰的垂直動量隨火焰高度的增加而增加,因此在火焰根部具有最明顯的水平偏轉(zhuǎn).同時由于火焰在水平方向的偏移,渦旋對的不對稱運動抽運了更多的空氣和燃燒產(chǎn)物,引起火焰在中部的斷裂,造成局部熄火并在水平氣流的作用下形成整體截斷狀態(tài)(b8),這與Pretrel等[19]的實驗現(xiàn)象一致;火焰在70 Hz與90 Hz(d序列和e序列)聲強迫下較為穩(wěn)定,圖像與無聲波狀態(tài)下的火焰相似,但火焰外圍輪廓亮度較大,火焰呈現(xiàn)尖、長的特點.由于聲波脈動頻率的增加影響了水平方向動量的積累,限制了渦旋的尺度,因而在火焰表面形成更多更小直徑的渦[20],促進(jìn)了空氣和可燃物的摻混使得火焰面燃燒更為完全,火焰也更為明亮.

    由圖2(中)可以看出,與3 cm直徑相比,4 cm直徑油池火焰由于熱釋放量的增加,在豎直方向的壓力梯度增加,無聲波狀態(tài)下(a序列)的膨脹行為(puffing)更為明顯;在30 Hz聲頻率作用下,與3 cm直徑油池火焰類似,火焰在不對稱渦旋的卷吸作用下出現(xiàn)橫向截斷(b序列);在50 Hz聲頻率作用下時,火焰在頂部呈現(xiàn)“T”型(c序列),可能是由于聲壓的正弦交替轉(zhuǎn)變至反方向,引起背聲側(cè)渦量的增大,原本在頂端處脫落的渦旋進(jìn)一步卷吸火焰,使原本依附于火焰輪廓線的渦旋能夠引起火焰面的褶皺變形;火焰在70 Hz和90 Hz聲波環(huán)境下(d序列和e序列),同樣處于燃燒較為穩(wěn)定的狀態(tài),亮度也略大于無聲波時.

    5 cm直徑油池的火焰(圖2右)在所有實驗頻率下都不存在明顯的橫向截斷狀態(tài),在30 Hz聲波作用下,火焰雖然發(fā)生偏移和斷裂,但與3 cm和4 cm直徑油池火不同,燃燒尺度的增加能更快補充被渦旋卷吸而破碎的火焰(b序列);此外,火焰面積的增加促使沿火焰面移動的渦旋從單一方向的運動向多方向的運動轉(zhuǎn)變(c序列);火焰在70 Hz和90 Hz聲波強迫下不發(fā)生偏移,處于穩(wěn)定燃燒狀態(tài),這與3 cm、4 cm直徑油池火相似.

    由于聲波引起的渦旋擾動以及火焰本身脈動不穩(wěn)定性會引起火焰形態(tài)的動態(tài)變化,因此當(dāng)火焰在特定聲頻率下出現(xiàn)間歇截斷的火焰形態(tài)時,定義火焰此時為間歇截斷狀態(tài),該聲強迫的頻率為間歇截斷頻率;出現(xiàn)火焰相對于中心軸線的偏轉(zhuǎn)時,則定義為偏轉(zhuǎn)狀態(tài);火焰在水平方向上不發(fā)生擾動時則為穩(wěn)定燃燒狀態(tài).總體而言,隨著聲波頻率的增加,渦旋逐漸由影響火焰結(jié)構(gòu)的內(nèi)部區(qū)域轉(zhuǎn)移至火焰表面的輪廓線,渦旋直徑也逐漸變小,引起的燃燒不穩(wěn)定性逐漸減弱.此外,3 cm和4 cm直徑油池火焰的火焰形態(tài)呈現(xiàn)由間歇截斷到偏轉(zhuǎn),再到穩(wěn)定燃燒的狀態(tài),5 cm直徑火焰不存在間歇截斷頻率,火焰形態(tài)直接由偏轉(zhuǎn)狀態(tài)向穩(wěn)定狀態(tài)轉(zhuǎn)變;而隨著油池直徑增加,渦旋移動的空間方向性顯著增加,而火焰對聲波的敏感性逐漸降低,聲波作用下的火焰形態(tài)變化示意圖如圖3和表1所示.

    圖3 油池火焰在聲波作用下火焰形態(tài)變化示意圖(I、II、III、和IV分別表示無聲波作用、間歇截斷、偏轉(zhuǎn)和穩(wěn)定燃燒的狀態(tài))Fig.3 Schematic diagram of flame shape change with acoustic force (I,II, III, and IV represent the nonacoustic, intermittent, deflective, and stable states, respectively)

    表1 不同直徑油池火焰在不同頻率聲波作用下的燃燒狀態(tài)Table 1 State of pool flames with different diameters and different frequencies of acoustic force

    2.2 火焰幾何尺寸

    由于原始火焰圖像的數(shù)量較大,為了進(jìn)一步定量分析火焰幾何尺寸,利用圖像二值化技術(shù)批量處理火焰圖像.其原理為:將圖像轉(zhuǎn)化為像素點的二維矩陣,再設(shè)置一個或多個灰度閾值,比較圖像中各像素的灰度值與閾值的大小關(guān)系,根據(jù)比較結(jié)果重新對矩陣賦值,從而將圖像的對應(yīng)像素分成兩類或多類[21].本文選取40作為亮度閾值將圖像的像素分為火焰部分(灰度值為1)和非火焰部分(灰度值為0),如圖4(a)所示,確定豎直方向與水平方向上灰度值為1的邊界點,根據(jù)像素點數(shù)量計算火焰尺寸,火焰面積定義火焰部分像素點數(shù)量之和,火焰寬度定義為邊界點在水平方向上所間隔的最大的像素點數(shù)量,火焰高度定義為邊界點在豎直方向上所間隔的最大的像素點數(shù)量.60 Hz頻率聲波作用下4 cm直徑油池的火焰面積如圖4(b)所示.由于聲強迫作用下火焰往復(fù)振蕩,為更好地反映火焰在一段時間內(nèi)的整體特征,計算聲波作用時間段內(nèi)火焰尺寸的平均數(shù)值,圖5為不同聲頻率下8 s內(nèi)的平均火焰面積、平均火焰寬度和平均火焰高度及標(biāo)準(zhǔn)差.

    圖4 火焰圖像二值化處理及火焰面積像素點數(shù).(a)圖像二值化處理;(b)火焰面積像素點數(shù)Fig.4 Binarization of flame image and number of flame area pixels: (a) binarization of flame image; (b) number of flame area pixels

    圖5 不同直徑油池火焰在不同聲頻率下的火焰面積、火焰寬度、火焰高度以及分形維數(shù)的平均數(shù)值.(a)火焰面積; (b)火焰寬度; (c)火焰高度;(d)分形維數(shù)Fig.5 Average flame area, flame width, flame height and fractal of pool flames with different diameters under acoustic wave of different frequencies: (a)flame area; (b) flame width; (c) flame height; (d) fractal dimension

    由圖5可以看出,3 cm和4 cm直徑油池火焰具有相同的幾何形態(tài)變化趨勢:在間歇截斷頻率下,平均火焰面積小幅度增加,在進(jìn)入偏轉(zhuǎn)頻率時火焰面積迅速下降,在進(jìn)入穩(wěn)定燃燒頻率時保持穩(wěn)定.5 cm直徑油池的平均火焰面積在偏轉(zhuǎn)頻率下持續(xù)增加,進(jìn)入穩(wěn)定燃燒頻率時面積下降并保持穩(wěn)定.這是由于5 cm直徑油池火焰具有更大火焰尺度,聲波擾動引起的流場擾動卷吸預(yù)混的助燃作用大于熱量損失帶來的影響.平均火焰寬度和平均火焰高度變化趨勢也印證了這一點,三種直徑油池火焰的平均寬度在增加的情況下,5 cm直徑油池的火焰高度相對穩(wěn)定,而3 cm和4 cm直徑油池火焰的火焰高度顯著降低,表明前者火焰體積增加而后者則無.當(dāng)頻率增加至70 Hz以上時,三種直徑油池的平均火焰寬度與自由燃燒時相近,但火焰高度和火焰面積都有所增加.由2.1的分析可知,隨著聲波頻率的增加,火焰表面的渦旋數(shù)量增多尺度減小,摻混效率和燃燒效率提升,羽流的上升速度在表面渦旋的作用下增加,因此火焰高度也有了提升.

    火焰在聲波作用下加劇幾何表面的皺褶程度,通過分形維數(shù)的計算能夠分析火焰面的皺褶情況[22],分形理論的數(shù)學(xué)表示為:

    式中,F(xiàn)為n維空間下的有界集合,F(xiàn)有Nr個互不重合的子集,r為所有坐標(biāo)方向上的尺度因子,DF為F的分形維數(shù).本文中分形維數(shù)計算采用傳統(tǒng)的數(shù)盒子法[23],利用MATLAB計算三種直徑的油池火焰在不同聲頻率下的平均分形維數(shù),計算結(jié)果如圖5(d)所示.

    由圖5(d)可知,火焰分形維數(shù)隨油池直徑的增大而增大,火焰羽流范圍的增加引起火焰面拓展,出現(xiàn)更多皺褶空間.三種直徑油池火焰均在30 Hz聲強迫下的平均分形維數(shù)最大,由前文分析可知,由于水平方向的壓力梯度變化引起渦量加強和渦旋運動顯著影響火焰區(qū)域,火焰處于高度不穩(wěn)定狀態(tài),火焰面彎曲皺褶,幾何特征復(fù)雜;與5 cm直徑油池不同的是,3 cm和4 cm直徑油池火焰均在聲強迫頻率為60 Hz時分形維數(shù)最小,這是由于火焰面積的減小限制了表面皺褶的能力.

    2.3 火焰信號頻域特性

    火焰的不穩(wěn)定燃燒可通過火焰的熱釋放速率、火焰高度、火焰面積等參數(shù)的頻域特性來反映[24].為分析火焰信號在不同狀態(tài)的頻域特性,對不同頻率的火焰寬度、火焰面積及火焰高度的時域信號進(jìn)行快速傅里葉變換(Fast fourier transform,F(xiàn)FT)分析.圖6為三種直徑油池火焰在30、50和70 Hz聲波條件下的快速傅里葉變化結(jié)果.

    圖6 不同直徑油池火焰在聲波強迫下的頻域信號(a、b和c前綴分別表示在30、50和70 Hz聲波作用下;1、2和3后綴分別表示火焰面積、火焰寬度和火焰高度信號)Fig.6 Frequency signals of pool flames with different diameters under an acoustic force (a, b, and c prefixes denote acoustic wave at 30, 50, and 70 Hz,respectively; 1, 2, and 3 suffixes indicate the flame area, flame width, and flame height signal, respectively)

    從幅值上看,火焰面積、寬度和高度數(shù)值的FFT結(jié)果中的幅值隨油池直徑的增加而增加;此外,輸入的聲波信號頻率越低,在低于火焰振蕩頻率的頻域內(nèi)(0 ~ 10 Hz)峰值越集中,表明低頻聲波對火焰結(jié)構(gòu)的破壞更為明顯,導(dǎo)致火焰信號極不穩(wěn)定;當(dāng)輸入的聲信號為70 Hz時,由于火焰重新趨于穩(wěn)定,頻域分布趨于自由燃燒下的頻域分布,但此時主頻和次頻的幅值降低,次頻不顯著,表明此時聲波降低了正常的火焰脈動幅度.

    從峰值頻率上看,間歇截斷狀態(tài)和穩(wěn)定燃燒狀態(tài)下的聲信號在火焰面積的頻域分布中不易識別,而偏轉(zhuǎn)狀態(tài)下峰值較為明顯,這是由于偏轉(zhuǎn)狀態(tài)下火焰結(jié)構(gòu)更為完整,周期規(guī)律性更強.火焰高度無法識別出70 Hz的聲信號,火焰寬度則不受限制,聲信號在火焰寬度頻域分布中始終突出.這是由于聲波方向與火焰羽流方向垂直,聲波擾動引起橫向壓力梯度相對于受控于重力和浮力作用引起的垂直壓力梯度更大,火焰在橫向方向的擾動更為顯著;70 Hz聲頻率下火焰的頻域分布除火焰自身的脈動頻率外還存在噪點,表明火焰在穩(wěn)定燃燒頻率下依然存在小振幅脈動.

    2.4 火焰理查德森數(shù)

    針對橫風(fēng)作用下的油池火焰,Majeski等[25]和Shang 等[26]通過 Richardson(Ri)數(shù),即火焰受到浮力與慣性力之比分析火焰不穩(wěn)定性.

    式中:Δρ為側(cè)風(fēng)密度與燃燒產(chǎn)物的密度差,kg·m-3;ρf為火焰密度,kg·m-3;L為油池特征尺寸,一般取油池直徑,m;g為重力加速度,m·s-2;U為側(cè)風(fēng)風(fēng)速,m·s-1.Ri數(shù)值越大,表示火焰處于越穩(wěn)定的狀態(tài),反之火焰越不穩(wěn)定.

    劉鵬翔[27]通過推導(dǎo)和研究發(fā)現(xiàn)小尺度油池火焰在環(huán)境風(fēng)作用下的Ri-1和火焰傾角的正切值存在一次函數(shù)關(guān)系:

    式中,tanθ為火焰傾角的正切值,定義為水平位移與和火焰高度的比值;A和B為常數(shù)項.

    與側(cè)風(fēng)類似,橫向聲波同樣會引起火焰水平動量擾動和火焰傾角的特征性變化.此外,聲波作用下的火焰不穩(wěn)定性還對頻率有明顯依賴性.因此借助火焰傾角與Ri數(shù)的關(guān)系分析聲波作用下的Ri數(shù)的形式,將式(3)的U帶入當(dāng)?shù)乜諝庹駝铀俣?,通過當(dāng)?shù)財_動聲壓的有效聲壓表示[21]:

    式中,P為當(dāng)?shù)財_動有效聲壓,Pa;ρ為介質(zhì)靜態(tài)密度,kg·m-3,通過下式得到:

    式中,P0為靜態(tài)氣壓,Pa;T為介質(zhì)溫度,K.經(jīng)由實驗測得3、4和5 cm直徑油池在自然狀態(tài)下的火焰中心溫度分別為384.2,453.6和564.7 ℃,取火焰中心溫度與28 ℃空氣溫度的平均溫度作為介質(zhì)溫度即火焰平均溫度,進(jìn)而得到介質(zhì)密度和以及與外側(cè)密度之差.

    式(4)中火焰傾角正切值在本文采用平均火焰寬度的偏移與平均火焰高度之比表示.帶入數(shù)值得到火焰傾角與聲頻率的關(guān)系如圖7(a)所示.由此認(rèn)為聲波頻率與火焰傾角的正切值存在近似反正切函數(shù)關(guān)系,在此基礎(chǔ)上增加頻率修正系數(shù)Z,使火焰傾角的正切值與Ria數(shù)的倒數(shù)滿足一次函數(shù)關(guān)系,即:

    圖7 聲波頻率、火焰傾角和Ria-1值之間的關(guān)系.(a)不同頻率下的火焰傾角正切值; (b) Ria-1值與火焰傾角正切值的擬合關(guān)系; (c)火焰Ria-1值Fig.7 Relationship between acoustic frequency, flame angle, and Ria-1: (a) tangent value of flame inclination angle at different frequencies; (b) fitting relationship between Ria-1 and the tangent value of flame inclination angle; (c) flame Ria-1 value

    式中,Ria表示聲波作用下的火焰Richardson數(shù);Z為頻率修正系數(shù),根據(jù)實驗結(jié)果擬合得到;ρ0為環(huán)境氣體密度,kg·m-3;f為聲頻率,Hz.火焰傾角正切值與Ria數(shù)倒數(shù)擬合得到的一次函數(shù)關(guān)系如圖7(b).為進(jìn)一步分析聲頻率對火焰不穩(wěn)定性的影響,得到聲頻率與Ria數(shù)的關(guān)系如圖7(c),為方便與無聲波作用下的火焰Ri值進(jìn)行對比(無聲波作用時空氣振動速度u≈0,Ri趨于無窮),Ria值均采用倒數(shù)的形式進(jìn)行討論.

    由圖7(c)可知,Ria-1值隨著油池的直徑增加而降低,對應(yīng)了油池尺度對聲波強迫的敏感性隨著尺度的增加而降低;當(dāng)外加聲壓大小相同時,隨著聲信號頻率的增加,Ria-1值呈降低的趨勢,對應(yīng)了火焰不穩(wěn)定性隨聲頻率的增加逐漸減弱的演變關(guān)系.此外,三種直徑油池火焰Ria-1值收斂于接近0的數(shù)值對應(yīng)當(dāng)聲波頻率超過一定數(shù)值時火焰近乎不再發(fā)生擾動,即與自然燃燒近乎無異(Ri=0).因此,增加頻率修正系數(shù)后的Ria數(shù)的結(jié)果與實驗現(xiàn)象有較好的符合,在描述聲波作用下的油池火焰穩(wěn)定性方面具有一定可靠性.在此基礎(chǔ)上可以發(fā)現(xiàn),在50 ~70 Hz之間,Ria-1值隨頻率增加下降的速率最大,表明此時火焰對聲波頻率的響應(yīng)尤為顯著,當(dāng)聲波頻率高于或低于該段頻率時,對火焰的影響存在邊際效應(yīng).此外,結(jié)合表1可以發(fā)現(xiàn),間歇截斷狀態(tài),偏轉(zhuǎn)狀態(tài)和穩(wěn)定燃燒狀態(tài)存在臨界Ria-1值,分別為10.32和2.92.

    3 結(jié)論

    (1)隨著聲波頻率的增加,三種直徑油池火的燃燒不穩(wěn)定性逐漸減弱,渦旋逐漸由影響火焰結(jié)構(gòu)的內(nèi)部區(qū)域轉(zhuǎn)移至火焰表面的輪廓線,渦旋直徑也逐漸變?。? cm和4 cm直徑油池火焰的間歇截斷和偏轉(zhuǎn)狀態(tài)的臨界頻率分別頻率為50 Hz和40 Hz,5 cm直徑油池火焰由于熱釋放速率的增加,對聲波擾動的敏感性降低,不存在間歇截斷頻率;三種油池火焰均在70 Hz時由偏轉(zhuǎn)狀態(tài)轉(zhuǎn)變?yōu)榉€(wěn)定燃燒狀態(tài).

    (2)3 cm直徑和4 cm直徑油池火焰在聲波作用下有相似的幾何形態(tài)變化趨勢,當(dāng)聲波頻率增加至70 Hz以上時,三種直徑油池的平均火焰寬度與自由燃燒時相近,火焰高度和火焰面積都有所增加;分形維數(shù)分析表明火焰在間歇截斷狀態(tài)下和不穩(wěn)定偏轉(zhuǎn)狀態(tài)下的皺褶程度最高.

    (3)火焰信號的頻域分析表明輸入的聲波信號頻率越低,在低于火焰振蕩頻率的頻域范圍內(nèi)峰值越集中;相較于間歇截斷狀態(tài)和穩(wěn)定燃燒狀態(tài),聲波頻率在偏轉(zhuǎn)狀態(tài)下的頻域信號中更明顯;聲頻率在火焰寬度的頻域分布中始終突出.

    (4)火焰傾角正切值與聲波頻率近似存在反正切函數(shù)關(guān)系,通過火焰傾角與Ri數(shù)的關(guān)系得到聲波作用下的火焰Ria數(shù)形式,Ria數(shù)值結(jié)果與實驗現(xiàn)象具有較好的符合;對于同一水平聲壓大小下的聲波,火焰存在對于頻率的顯著響應(yīng)區(qū)間;間歇截斷狀態(tài),偏轉(zhuǎn)狀態(tài)和穩(wěn)定燃燒狀態(tài)存在臨界Ria-1值,分別為10.32和2.92.

    猜你喜歡
    油池渦旋頻域
    油浸變壓器事故油池的優(yōu)化設(shè)計
    球盤點接觸區(qū)外潤滑油分布的試驗研究
    基于PM算法的渦旋電磁波引信超分辨測向方法
    “L”形油池火旋風(fēng)仿真和中尺度試驗研究
    頻域稀疏毫米波人體安檢成像處理和快速成像稀疏陣列設(shè)計
    光渦旋方程解的存在性研究
    變電站主變事故油池設(shè)計施工新思考
    基于改進(jìn)Radon-Wigner變換的目標(biāo)和拖曳式誘餌頻域分離
    一種基于頻域的QPSK窄帶干擾抑制算法
    基于頻域伸縮的改進(jìn)DFT算法
    電測與儀表(2015年3期)2015-04-09 11:37:24
    国产97色在线日韩免费| 久久人人爽av亚洲精品天堂| 黄色片一级片一级黄色片| 亚洲七黄色美女视频| 亚洲男人天堂网一区| 99国产极品粉嫩在线观看| 九色亚洲精品在线播放| 精品高清国产在线一区| 多毛熟女@视频| 美女 人体艺术 gogo| av超薄肉色丝袜交足视频| 精品久久久久久久毛片微露脸| 久久久精品免费免费高清| 丝瓜视频免费看黄片| av免费在线观看网站| 午夜日韩欧美国产| 男人的好看免费观看在线视频 | 看片在线看免费视频| 无人区码免费观看不卡| 最近最新中文字幕大全免费视频| 欧美日本中文国产一区发布| 欧美在线黄色| 日本黄色视频三级网站网址 | 免费在线观看完整版高清| 久久国产乱子伦精品免费另类| 手机成人av网站| 一级毛片高清免费大全| 欧美精品av麻豆av| 亚洲一区高清亚洲精品| 国产精品偷伦视频观看了| 男女之事视频高清在线观看| 亚洲熟妇中文字幕五十中出 | 日韩欧美一区视频在线观看| 午夜视频精品福利| 一二三四在线观看免费中文在| 丁香六月欧美| 国产片内射在线| 狂野欧美激情性xxxx| 日韩有码中文字幕| 日韩大码丰满熟妇| 亚洲国产精品一区二区三区在线| 日韩熟女老妇一区二区性免费视频| 国产精品久久久人人做人人爽| 久久国产乱子伦精品免费另类| 少妇粗大呻吟视频| 麻豆成人av在线观看| 高清欧美精品videossex| 国产1区2区3区精品| 国产日韩一区二区三区精品不卡| 99久久综合精品五月天人人| 日本一区二区免费在线视频| 国产精华一区二区三区| 啪啪无遮挡十八禁网站| 99国产综合亚洲精品| 老汉色av国产亚洲站长工具| 国产又色又爽无遮挡免费看| 亚洲伊人色综图| ponron亚洲| 精品国产亚洲在线| 欧美日韩av久久| 高清欧美精品videossex| 一级毛片高清免费大全| 夜夜躁狠狠躁天天躁| 欧美国产精品一级二级三级| 18禁裸乳无遮挡动漫免费视频| 99热国产这里只有精品6| 91国产中文字幕| 亚洲一码二码三码区别大吗| 亚洲人成77777在线视频| 亚洲av成人av| a级片在线免费高清观看视频| 国产又色又爽无遮挡免费看| 亚洲成人免费av在线播放| 好男人电影高清在线观看| 国产精品久久久久久精品古装| 麻豆成人av在线观看| 777久久人妻少妇嫩草av网站| 欧美乱色亚洲激情| 午夜日韩欧美国产| 国产1区2区3区精品| 国产99久久九九免费精品| 久久国产精品人妻蜜桃| 国产高清videossex| 国产精品免费大片| 国产亚洲欧美精品永久| 女人被狂操c到高潮| 亚洲人成伊人成综合网2020| 国产极品粉嫩免费观看在线| 国产av又大| 青草久久国产| 十八禁人妻一区二区| 久久精品91无色码中文字幕| 国产成人精品久久二区二区91| 久久人妻熟女aⅴ| 精品国产乱码久久久久久男人| 午夜福利在线免费观看网站| 人人妻人人澡人人爽人人夜夜| 亚洲精品粉嫩美女一区| 99re6热这里在线精品视频| 国产91精品成人一区二区三区| 黄色女人牲交| 国产免费av片在线观看野外av| 久热这里只有精品99| 久久久精品免费免费高清| 国产主播在线观看一区二区| 一区福利在线观看| av在线播放免费不卡| 国产成人免费无遮挡视频| 久久国产精品人妻蜜桃| 成人18禁在线播放| 在线播放国产精品三级| 国产精品免费大片| svipshipincom国产片| 精品欧美一区二区三区在线| 狂野欧美激情性xxxx| 777米奇影视久久| av网站在线播放免费| 法律面前人人平等表现在哪些方面| 亚洲精品国产一区二区精华液| 色综合欧美亚洲国产小说| 亚洲avbb在线观看| www.精华液| 成人18禁在线播放| 成人影院久久| 最新美女视频免费是黄的| 成人精品一区二区免费| 久久久久国产精品人妻aⅴ院 | 国内毛片毛片毛片毛片毛片| 免费在线观看完整版高清| 1024香蕉在线观看| 在线观看免费视频网站a站| 涩涩av久久男人的天堂| 久久久久久久国产电影| 国产又色又爽无遮挡免费看| 午夜免费成人在线视频| 国产无遮挡羞羞视频在线观看| 又黄又爽又免费观看的视频| 国产激情久久老熟女| 美女 人体艺术 gogo| 在线观看免费午夜福利视频| 视频区图区小说| 成人国语在线视频| 18禁裸乳无遮挡动漫免费视频| 女警被强在线播放| av超薄肉色丝袜交足视频| 露出奶头的视频| 日韩有码中文字幕| 国产成人影院久久av| 国内毛片毛片毛片毛片毛片| 黄色视频不卡| 少妇被粗大的猛进出69影院| 国产野战对白在线观看| 免费在线观看影片大全网站| 777久久人妻少妇嫩草av网站| 老熟妇乱子伦视频在线观看| 又黄又爽又免费观看的视频| 国产精品免费视频内射| 热re99久久精品国产66热6| 女人被躁到高潮嗷嗷叫费观| 亚洲全国av大片| 老熟女久久久| 日韩成人在线观看一区二区三区| av在线播放免费不卡| 国产不卡av网站在线观看| www.999成人在线观看| 国产精品一区二区在线观看99| 精品国产一区二区三区久久久樱花| 男女免费视频国产| 在线观看午夜福利视频| 久9热在线精品视频| 性少妇av在线| 高清av免费在线| 精品第一国产精品| 一级毛片精品| 好看av亚洲va欧美ⅴa在| 岛国在线观看网站| 女同久久另类99精品国产91| 久久久久久人人人人人| 91成人精品电影| av视频免费观看在线观看| 99re6热这里在线精品视频| 午夜福利免费观看在线| 国产成人精品久久二区二区91| 超色免费av| 美女 人体艺术 gogo| 高清av免费在线| 精品少妇久久久久久888优播| 99久久99久久久精品蜜桃| 9热在线视频观看99| 999久久久精品免费观看国产| 成人国产一区最新在线观看| 亚洲国产看品久久| 高清视频免费观看一区二区| 国产男女超爽视频在线观看| 夜夜爽天天搞| 久久香蕉激情| 成人国语在线视频| 国产aⅴ精品一区二区三区波| 国产1区2区3区精品| www日本在线高清视频| 妹子高潮喷水视频| 嫩草影视91久久| 啦啦啦免费观看视频1| 深夜精品福利| 黑丝袜美女国产一区| 亚洲一区二区三区欧美精品| 99香蕉大伊视频| 极品人妻少妇av视频| 黑人猛操日本美女一级片| 69精品国产乱码久久久| 老司机午夜十八禁免费视频| 国产精品欧美亚洲77777| 亚洲 欧美一区二区三区| 在线十欧美十亚洲十日本专区| 国产亚洲一区二区精品| av中文乱码字幕在线| 91麻豆精品激情在线观看国产 | 国产欧美亚洲国产| 午夜亚洲福利在线播放| 十八禁人妻一区二区| 亚洲成a人片在线一区二区| 日韩欧美在线二视频 | 99精品欧美一区二区三区四区| 亚洲av电影在线进入| 久久久久久久精品吃奶| 精品熟女少妇八av免费久了| 久久久国产成人精品二区 | 黄色a级毛片大全视频| 亚洲午夜理论影院| 国产视频一区二区在线看| 999精品在线视频| 一a级毛片在线观看| 精品人妻熟女毛片av久久网站| 日韩欧美免费精品| 大香蕉久久成人网| 精品久久蜜臀av无| 亚洲 欧美一区二区三区| 国产成人免费无遮挡视频| 老司机深夜福利视频在线观看| 亚洲精品国产区一区二| 日韩 欧美 亚洲 中文字幕| 亚洲中文av在线| 中国美女看黄片| 99久久精品国产亚洲精品| 久久人人爽av亚洲精品天堂| 亚洲va日本ⅴa欧美va伊人久久| 黄色 视频免费看| 在线观看免费视频网站a站| 国产免费现黄频在线看| 十八禁人妻一区二区| 狂野欧美激情性xxxx| 丁香六月欧美| 黑人巨大精品欧美一区二区mp4| 亚洲少妇的诱惑av| 免费在线观看完整版高清| 亚洲av第一区精品v没综合| 亚洲人成电影免费在线| 伦理电影免费视频| 超碰成人久久| 女人被躁到高潮嗷嗷叫费观| 丰满的人妻完整版| 少妇裸体淫交视频免费看高清 | 777久久人妻少妇嫩草av网站| 婷婷精品国产亚洲av在线 | 精品第一国产精品| 成年人午夜在线观看视频| 日韩制服丝袜自拍偷拍| 午夜免费观看网址| 精品免费久久久久久久清纯 | 老汉色∧v一级毛片| 国产蜜桃级精品一区二区三区 | 美女福利国产在线| 久久久精品区二区三区| 午夜福利乱码中文字幕| 免费高清在线观看日韩| 成人黄色视频免费在线看| av一本久久久久| 欧美国产精品va在线观看不卡| 久久人人97超碰香蕉20202| 亚洲成人国产一区在线观看| 久久精品国产亚洲av高清一级| 老司机午夜十八禁免费视频| 人妻丰满熟妇av一区二区三区 | 欧美黑人精品巨大| 黄频高清免费视频| 精品国产一区二区久久| 亚洲在线自拍视频| 水蜜桃什么品种好| 免费黄频网站在线观看国产| 久久久国产精品麻豆| 香蕉国产在线看| av天堂久久9| 国产精品98久久久久久宅男小说| 精品亚洲成a人片在线观看| 美女高潮到喷水免费观看| 久久亚洲真实| 久久久久久亚洲精品国产蜜桃av| 国产精品电影一区二区三区 | 亚洲精品美女久久久久99蜜臀| 中国美女看黄片| 欧美中文综合在线视频| 老司机亚洲免费影院| videos熟女内射| 日韩视频一区二区在线观看| 在线观看免费午夜福利视频| 国产人伦9x9x在线观看| 亚洲午夜精品一区,二区,三区| 天天躁狠狠躁夜夜躁狠狠躁| 天天添夜夜摸| 无人区码免费观看不卡| 国产精品一区二区免费欧美| 国产乱人伦免费视频| 天天影视国产精品| 99精国产麻豆久久婷婷| 女警被强在线播放| 亚洲精品久久午夜乱码| 纯流量卡能插随身wifi吗| 亚洲性夜色夜夜综合| 一进一出抽搐动态| 国产男女超爽视频在线观看| 亚洲伊人色综图| 久久精品人人爽人人爽视色| 美女扒开内裤让男人捅视频| 真人做人爱边吃奶动态| 人人妻,人人澡人人爽秒播| 欧美老熟妇乱子伦牲交| 精品国产超薄肉色丝袜足j| 一a级毛片在线观看| 俄罗斯特黄特色一大片| 性少妇av在线| 亚洲午夜理论影院| 国产日韩欧美亚洲二区| 精品免费久久久久久久清纯 | 日日摸夜夜添夜夜添小说| 欧美亚洲日本最大视频资源| 狂野欧美激情性xxxx| 免费久久久久久久精品成人欧美视频| 黄网站色视频无遮挡免费观看| 美女高潮到喷水免费观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲色图av天堂| 亚洲精品粉嫩美女一区| 一级片'在线观看视频| 老司机靠b影院| 久久精品国产亚洲av高清一级| 黄色毛片三级朝国网站| 丁香欧美五月| 在线观看免费午夜福利视频| 久久久久国产精品人妻aⅴ院 | 啦啦啦视频在线资源免费观看| 精品视频人人做人人爽| 热99re8久久精品国产| 亚洲精品一卡2卡三卡4卡5卡| 久久国产精品大桥未久av| av天堂久久9| 亚洲国产毛片av蜜桃av| 男人操女人黄网站| 99国产精品一区二区三区| 亚洲精品国产精品久久久不卡| 成人18禁在线播放| 亚洲精品在线观看二区| 99riav亚洲国产免费| 在线十欧美十亚洲十日本专区| 色94色欧美一区二区| 最新美女视频免费是黄的| 亚洲全国av大片| 国产伦人伦偷精品视频| 成熟少妇高潮喷水视频| 老司机亚洲免费影院| 午夜福利欧美成人| 夜夜爽天天搞| 操美女的视频在线观看| 啦啦啦在线免费观看视频4| 亚洲三区欧美一区| 日韩欧美一区视频在线观看| 精品国产乱码久久久久久男人| 午夜老司机福利片| 人妻丰满熟妇av一区二区三区 | 老司机在亚洲福利影院| 精品视频人人做人人爽| 中国美女看黄片| 免费av中文字幕在线| 岛国在线观看网站| 精品国产一区二区三区四区第35| 亚洲人成77777在线视频| 亚洲欧美激情综合另类| 欧美乱码精品一区二区三区| av视频免费观看在线观看| 如日韩欧美国产精品一区二区三区| 免费一级毛片在线播放高清视频 | 亚洲中文av在线| 欧美黄色淫秽网站| 色老头精品视频在线观看| 久久天堂一区二区三区四区| 午夜精品在线福利| 色综合婷婷激情| 99re在线观看精品视频| 最新的欧美精品一区二区| 精品国内亚洲2022精品成人 | 俄罗斯特黄特色一大片| 欧美中文综合在线视频| 999精品在线视频| 久久中文字幕人妻熟女| 日韩欧美在线二视频 | 青草久久国产| 在线国产一区二区在线| 巨乳人妻的诱惑在线观看| 超色免费av| 午夜日韩欧美国产| 久久精品熟女亚洲av麻豆精品| 桃红色精品国产亚洲av| 丁香六月欧美| 精品国产一区二区三区四区第35| 一区二区三区精品91| 亚洲国产欧美网| 亚洲精品久久成人aⅴ小说| 久久精品亚洲av国产电影网| 亚洲国产欧美网| 人妻 亚洲 视频| 日韩有码中文字幕| 757午夜福利合集在线观看| 亚洲av美国av| 亚洲免费av在线视频| a级毛片黄视频| www.999成人在线观看| 亚洲一区高清亚洲精品| 国产精品98久久久久久宅男小说| 成年人免费黄色播放视频| 精品久久蜜臀av无| 999久久久精品免费观看国产| 国产精品自产拍在线观看55亚洲 | 久久精品熟女亚洲av麻豆精品| 国产av又大| 超碰成人久久| 亚洲熟妇中文字幕五十中出 | 黄色女人牲交| 久久久久久久精品吃奶| 精品国产乱子伦一区二区三区| 涩涩av久久男人的天堂| 又大又爽又粗| 亚洲av片天天在线观看| 日韩精品免费视频一区二区三区| 欧美人与性动交α欧美精品济南到| 高清欧美精品videossex| 午夜福利,免费看| 欧美黄色淫秽网站| 一区二区三区激情视频| 叶爱在线成人免费视频播放| 天天操日日干夜夜撸| 搡老乐熟女国产| 国产野战对白在线观看| 久久性视频一级片| 99国产精品一区二区三区| www.999成人在线观看| 99久久国产精品久久久| 亚洲专区国产一区二区| 欧美老熟妇乱子伦牲交| 一边摸一边抽搐一进一小说 | 一区福利在线观看| 夜夜爽天天搞| 王馨瑶露胸无遮挡在线观看| 1024香蕉在线观看| 首页视频小说图片口味搜索| 国产无遮挡羞羞视频在线观看| 亚洲欧美一区二区三区黑人| 国产99久久九九免费精品| av欧美777| 三级毛片av免费| 亚洲欧美激情在线| 国产伦人伦偷精品视频| 久久青草综合色| 国产高清激情床上av| 精品国产一区二区久久| 一级黄色大片毛片| 中文字幕人妻丝袜制服| 精品国内亚洲2022精品成人 | 久久久久久久久久久久大奶| 色综合婷婷激情| 亚洲av美国av| 国产男女内射视频| 黑人巨大精品欧美一区二区蜜桃| a在线观看视频网站| 欧美精品啪啪一区二区三区| 国产精品.久久久| 久久人妻av系列| 精品一区二区三卡| 日韩一卡2卡3卡4卡2021年| 不卡一级毛片| 亚洲五月天丁香| 欧美 亚洲 国产 日韩一| 亚洲精品国产精品久久久不卡| 国产人伦9x9x在线观看| 欧美在线黄色| 日本撒尿小便嘘嘘汇集6| 久久婷婷成人综合色麻豆| 精品福利永久在线观看| 国产成+人综合+亚洲专区| 99在线人妻在线中文字幕 | 久久久精品免费免费高清| 成人国语在线视频| 久久人人97超碰香蕉20202| 91国产中文字幕| 97人妻天天添夜夜摸| 久久人妻熟女aⅴ| 在线国产一区二区在线| 亚洲精品美女久久久久99蜜臀| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产毛片av蜜桃av| 国产成人一区二区三区免费视频网站| 亚洲精品在线美女| 国产精品久久视频播放| 老司机亚洲免费影院| 俄罗斯特黄特色一大片| 看片在线看免费视频| www.999成人在线观看| 久久午夜综合久久蜜桃| 在线观看免费视频网站a站| 成人国语在线视频| 久久精品国产a三级三级三级| 免费观看a级毛片全部| 亚洲熟女毛片儿| 中国美女看黄片| 欧美日韩av久久| 满18在线观看网站| 9191精品国产免费久久| 国产三级黄色录像| 脱女人内裤的视频| 亚洲av美国av| 国产日韩一区二区三区精品不卡| 精品人妻在线不人妻| 欧美日韩亚洲国产一区二区在线观看 | 亚洲熟女毛片儿| 国产精品成人在线| 亚洲三区欧美一区| 久久天躁狠狠躁夜夜2o2o| 亚洲精品一卡2卡三卡4卡5卡| 国产成人影院久久av| 国产精品亚洲一级av第二区| 高清黄色对白视频在线免费看| 成人永久免费在线观看视频| 90打野战视频偷拍视频| 亚洲精品av麻豆狂野| 欧美大码av| 成年人免费黄色播放视频| 人人澡人人妻人| 丰满的人妻完整版| 亚洲avbb在线观看| 欧美av亚洲av综合av国产av| 美女 人体艺术 gogo| 12—13女人毛片做爰片一| 一进一出抽搐动态| 午夜影院日韩av| 精品国产一区二区久久| 丰满饥渴人妻一区二区三| 高清黄色对白视频在线免费看| 热re99久久国产66热| 欧美黑人欧美精品刺激| 高清av免费在线| 国产日韩欧美亚洲二区| 在线观看日韩欧美| 日韩中文字幕欧美一区二区| 亚洲欧美一区二区三区黑人| 午夜激情av网站| 欧美日韩成人在线一区二区| 亚洲欧美精品综合一区二区三区| 18在线观看网站| 老司机靠b影院| 亚洲人成电影免费在线| ponron亚洲| 18禁黄网站禁片午夜丰满| 伊人久久大香线蕉亚洲五| 人人妻人人澡人人看| 老熟女久久久| 十分钟在线观看高清视频www| 午夜福利视频在线观看免费| 最新在线观看一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 狂野欧美激情性xxxx| 国产不卡一卡二| 欧美国产精品一级二级三级| 他把我摸到了高潮在线观看| 国产精品一区二区在线不卡| 午夜福利,免费看| 黄色a级毛片大全视频| 老汉色av国产亚洲站长工具| 窝窝影院91人妻| 一个人免费在线观看的高清视频| 一本综合久久免费| 国产成人精品无人区| 成年人午夜在线观看视频| 国产精品 国内视频| 久久国产精品影院| 男人的好看免费观看在线视频 | 国产男女内射视频| 两个人免费观看高清视频| 欧美精品av麻豆av| 深夜精品福利| 一级毛片精品| 99精国产麻豆久久婷婷| 久久精品亚洲熟妇少妇任你| 老鸭窝网址在线观看| 亚洲精品成人av观看孕妇| 欧美成狂野欧美在线观看| 老熟妇仑乱视频hdxx| 99久久精品国产亚洲精品| 国产精品久久久av美女十八| 欧美另类亚洲清纯唯美| 午夜免费观看网址| xxxhd国产人妻xxx| 国产精品一区二区在线观看99| 国产在视频线精品| 交换朋友夫妻互换小说| 亚洲熟女毛片儿| 国产极品粉嫩免费观看在线| 免费黄频网站在线观看国产|