• <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
    亚洲成人国产一区在线观看| 人人妻人人澡欧美一区二区| 国产精品影院久久| 日本 av在线| 欧美性猛交黑人性爽| 免费看十八禁软件| 99久久精品国产亚洲精品| av片东京热男人的天堂| 久久久久久久久中文| 日韩精品中文字幕看吧| 韩国精品一区二区三区| 久久久久久久午夜电影| 免费看十八禁软件| 中文在线观看免费www的网站 | 亚洲成国产人片在线观看| 黄色a级毛片大全视频| 少妇 在线观看| 91字幕亚洲| 亚洲国产精品999在线| 国产精品香港三级国产av潘金莲| 久久狼人影院| 亚洲国产看品久久| 在线永久观看黄色视频| 亚洲第一青青草原| 久久久久久久午夜电影| 少妇裸体淫交视频免费看高清 | 在线视频色国产色| 听说在线观看完整版免费高清| 亚洲精华国产精华精| 中文字幕人妻熟女乱码| 久久久精品国产亚洲av高清涩受| 99久久精品国产亚洲精品| 午夜福利成人在线免费观看| 老司机午夜福利在线观看视频| 亚洲第一电影网av| 日本在线视频免费播放| 九色国产91popny在线| 啦啦啦韩国在线观看视频| 国产又黄又爽又无遮挡在线| 黄片播放在线免费| 中文字幕精品免费在线观看视频| 99久久99久久久精品蜜桃| 黄色a级毛片大全视频| 亚洲av成人一区二区三| 欧美精品啪啪一区二区三区| 看黄色毛片网站| 99久久99久久久精品蜜桃| 国产激情欧美一区二区| 亚洲国产欧美一区二区综合| a级毛片在线看网站| 欧美在线一区亚洲| 两性夫妻黄色片| 非洲黑人性xxxx精品又粗又长| 亚洲精品国产一区二区精华液| 成人18禁高潮啪啪吃奶动态图| 婷婷精品国产亚洲av| 国产精品1区2区在线观看.| 国产黄片美女视频| 久久久久久免费高清国产稀缺| 欧美+亚洲+日韩+国产| 久久久久亚洲av毛片大全| 国产蜜桃级精品一区二区三区| 精品乱码久久久久久99久播| 免费在线观看黄色视频的| 欧美国产精品va在线观看不卡| 国产精品亚洲美女久久久| 757午夜福利合集在线观看| 欧美大码av| 中文字幕av电影在线播放| 免费在线观看视频国产中文字幕亚洲| 一级黄色大片毛片| 此物有八面人人有两片| 美女午夜性视频免费| 黄色视频不卡| 黄频高清免费视频| 午夜精品在线福利| 欧美黄色淫秽网站| 免费观看精品视频网站| 少妇熟女aⅴ在线视频| 动漫黄色视频在线观看| 欧美又色又爽又黄视频| 亚洲狠狠婷婷综合久久图片| 国产精品亚洲av一区麻豆| 免费在线观看视频国产中文字幕亚洲| 婷婷六月久久综合丁香| 18禁裸乳无遮挡免费网站照片 | 欧美不卡视频在线免费观看 | 亚洲精品中文字幕一二三四区| 桃色一区二区三区在线观看| 亚洲av片天天在线观看| 欧美在线黄色| 俺也久久电影网| 九色国产91popny在线| 色哟哟哟哟哟哟| 亚洲熟妇中文字幕五十中出| 久久久久九九精品影院| 中文字幕最新亚洲高清| 999久久久精品免费观看国产| 哪里可以看免费的av片| 亚洲成a人片在线一区二区| av电影中文网址| 欧美在线黄色| tocl精华| 又黄又粗又硬又大视频| 欧美黑人精品巨大| av超薄肉色丝袜交足视频| 美女扒开内裤让男人捅视频| 精品少妇一区二区三区视频日本电影| 国产欧美日韩精品亚洲av| 在线观看日韩欧美| 香蕉久久夜色| 国产免费av片在线观看野外av| 免费看十八禁软件| 亚洲中文日韩欧美视频| 一区二区日韩欧美中文字幕| 后天国语完整版免费观看| 一区二区三区国产精品乱码| 99re在线观看精品视频| 欧美人与性动交α欧美精品济南到| 亚洲在线自拍视频| 人人妻,人人澡人人爽秒播| a级毛片a级免费在线| 日本撒尿小便嘘嘘汇集6| 韩国av一区二区三区四区| 欧美在线黄色| 麻豆国产av国片精品| 黄色毛片三级朝国网站| 国产亚洲精品第一综合不卡| 日韩av在线大香蕉| 50天的宝宝边吃奶边哭怎么回事| 最好的美女福利视频网| 夜夜爽天天搞| 最近最新中文字幕大全电影3 | 欧美黑人巨大hd| 精品第一国产精品| 免费看a级黄色片| 18禁美女被吸乳视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品中文字幕在线视频| 亚洲欧美一区二区三区黑人| 午夜激情福利司机影院| 日韩欧美三级三区| 国产精华一区二区三区| 亚洲av成人一区二区三| 久久天堂一区二区三区四区| 国产极品粉嫩免费观看在线| 欧美日韩精品网址| 免费av毛片视频| 欧美乱妇无乱码| 黄片播放在线免费| 亚洲第一av免费看| 欧美日韩一级在线毛片| 欧美激情久久久久久爽电影| 国产欧美日韩精品亚洲av| 国产日本99.免费观看| 精品卡一卡二卡四卡免费| 免费看a级黄色片| 久久人妻福利社区极品人妻图片| 一进一出抽搐动态| 欧美中文日本在线观看视频| avwww免费| 国产黄a三级三级三级人| 成在线人永久免费视频| 在线观看午夜福利视频| 亚洲成人久久爱视频| 夜夜躁狠狠躁天天躁| 亚洲色图av天堂| 国产野战对白在线观看| 亚洲精品美女久久av网站| 午夜福利视频1000在线观看| 伦理电影免费视频| 色婷婷久久久亚洲欧美| 波多野结衣高清无吗| 国产精品一区二区三区四区久久 | 一进一出抽搐动态| 成人一区二区视频在线观看| 日本黄色视频三级网站网址| avwww免费| 欧美在线一区亚洲| 成人亚洲精品av一区二区| 国产精品亚洲美女久久久| 日韩欧美一区视频在线观看| 波多野结衣高清无吗| 精品日产1卡2卡| 午夜福利成人在线免费观看| 久久久精品国产亚洲av高清涩受| 国产精品一区二区免费欧美| 国产精品98久久久久久宅男小说| 久久国产乱子伦精品免费另类| 国产精品爽爽va在线观看网站 | 午夜福利成人在线免费观看| 女性被躁到高潮视频| 日韩精品免费视频一区二区三区| 亚洲七黄色美女视频| 欧美色欧美亚洲另类二区| 国产单亲对白刺激| 亚洲专区字幕在线| 巨乳人妻的诱惑在线观看| 成人午夜高清在线视频 | 激情在线观看视频在线高清| 国产真人三级小视频在线观看| 99国产综合亚洲精品| 女警被强在线播放| 成人永久免费在线观看视频| 黄色 视频免费看| 久热这里只有精品99| 国产真实乱freesex| 男男h啪啪无遮挡| 巨乳人妻的诱惑在线观看| 香蕉国产在线看| 亚洲一区二区三区不卡视频| 欧美激情极品国产一区二区三区| 99热这里只有精品一区 | 在线av久久热| 精品国产超薄肉色丝袜足j| 天天添夜夜摸| 观看免费一级毛片| 啦啦啦免费观看视频1| 99久久无色码亚洲精品果冻| 国产成人精品久久二区二区免费| 国产精品二区激情视频| 国产精品亚洲av一区麻豆| 日本黄色视频三级网站网址| 午夜亚洲福利在线播放| 日本一本二区三区精品| 少妇的丰满在线观看| 亚洲午夜精品一区,二区,三区| 亚洲国产看品久久| 国产亚洲av嫩草精品影院| 日日爽夜夜爽网站| 欧美+亚洲+日韩+国产| 男女之事视频高清在线观看| www.精华液| 日韩中文字幕欧美一区二区| 欧美黑人欧美精品刺激| 亚洲成人久久爱视频| 熟妇人妻久久中文字幕3abv| 99re在线观看精品视频| 免费搜索国产男女视频| 亚洲人成网站高清观看| 999久久久国产精品视频| 欧美性猛交黑人性爽| 成年女人毛片免费观看观看9| 国产免费男女视频| 免费av毛片视频| 好男人在线观看高清免费视频 | 人人妻人人澡人人看| 国产亚洲精品av在线| 好男人电影高清在线观看| 国产av不卡久久| 在线观看免费午夜福利视频| 亚洲欧美日韩高清在线视频| 中文字幕高清在线视频| 午夜福利18| 亚洲一码二码三码区别大吗| 99re在线观看精品视频| 99在线人妻在线中文字幕| 成人手机av| 午夜免费鲁丝| 国产精品综合久久久久久久免费| www日本黄色视频网| 老汉色∧v一级毛片| 两个人免费观看高清视频| 欧美av亚洲av综合av国产av| av中文乱码字幕在线| 久久亚洲精品不卡| 伊人久久大香线蕉亚洲五| 亚洲av熟女| 国产人伦9x9x在线观看| 欧美日韩黄片免| 久久精品人妻少妇| 两个人视频免费观看高清| 丰满人妻熟妇乱又伦精品不卡| 少妇熟女aⅴ在线视频| 婷婷六月久久综合丁香| 国产色视频综合| 麻豆一二三区av精品| 我的亚洲天堂| 一进一出抽搐gif免费好疼| 日韩欧美一区二区三区在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲在线自拍视频| 男人舔奶头视频| 岛国视频午夜一区免费看| 一进一出抽搐动态| 亚洲自偷自拍图片 自拍| 中文字幕人成人乱码亚洲影| 日韩av在线大香蕉| 黄片播放在线免费| 日韩欧美国产一区二区入口| а√天堂www在线а√下载| 久久精品国产99精品国产亚洲性色| 亚洲色图av天堂| www.999成人在线观看| 18禁裸乳无遮挡免费网站照片 | 国产亚洲精品一区二区www| 亚洲全国av大片| 日本a在线网址| 久久伊人香网站| 午夜两性在线视频| 亚洲人成77777在线视频| 国产单亲对白刺激| 午夜两性在线视频| 90打野战视频偷拍视频| 国产精品亚洲av一区麻豆| 一级a爱视频在线免费观看| 久久久水蜜桃国产精品网| 制服诱惑二区| 国产激情欧美一区二区| 国产精品久久久av美女十八| 午夜福利在线观看吧| 日本一本二区三区精品| 婷婷精品国产亚洲av在线| 欧美日韩乱码在线| 久久国产精品人妻蜜桃| 黄色 视频免费看| 免费观看人在逋| 欧美黄色淫秽网站| 国产亚洲精品久久久久久毛片| 久久久久久大精品| 国产激情久久老熟女| 动漫黄色视频在线观看| 黄色视频,在线免费观看| 国产av一区在线观看免费| 在线观看www视频免费| 国产午夜精品久久久久久| 久久草成人影院| 女性被躁到高潮视频| www.999成人在线观看| 女人被狂操c到高潮| 精品一区二区三区视频在线观看免费| 亚洲国产看品久久| 在线观看日韩欧美| 手机成人av网站| 免费电影在线观看免费观看| 一边摸一边做爽爽视频免费| 搡老岳熟女国产| 免费搜索国产男女视频| 一级毛片女人18水好多| 亚洲精品一区av在线观看| 一个人观看的视频www高清免费观看 | 可以在线观看的亚洲视频| 国产亚洲av嫩草精品影院| 久久久久久大精品| 后天国语完整版免费观看| 亚洲av片天天在线观看| 无遮挡黄片免费观看| а√天堂www在线а√下载| 一卡2卡三卡四卡精品乱码亚洲| 一二三四在线观看免费中文在| 动漫黄色视频在线观看| 亚洲精品在线美女| 哪里可以看免费的av片| 最新在线观看一区二区三区| 曰老女人黄片| 亚洲av成人av| 国产区一区二久久| 又紧又爽又黄一区二区| 亚洲自拍偷在线| 日韩 欧美 亚洲 中文字幕| 欧美色视频一区免费| 亚洲av成人一区二区三| 在线观看一区二区三区| 99久久国产精品久久久| 国产乱人伦免费视频| 亚洲久久久国产精品| 亚洲成a人片在线一区二区| 亚洲三区欧美一区| videosex国产| 久久久久精品国产欧美久久久| 999久久久国产精品视频| 欧美黑人巨大hd| 亚洲专区字幕在线| 日韩精品免费视频一区二区三区| 国产精品美女特级片免费视频播放器 | 嫩草影视91久久| 亚洲av熟女| 真人做人爱边吃奶动态| 日本五十路高清| 午夜福利视频1000在线观看| 欧美乱色亚洲激情| 91麻豆精品激情在线观看国产| 国产成人一区二区三区免费视频网站| 少妇 在线观看| 听说在线观看完整版免费高清| 观看免费一级毛片| 变态另类成人亚洲欧美熟女| www日本黄色视频网| 在线观看日韩欧美| 亚洲午夜精品一区,二区,三区| 琪琪午夜伦伦电影理论片6080| svipshipincom国产片| 国产av在哪里看| 麻豆成人午夜福利视频| 日本熟妇午夜| 免费高清视频大片| 国产97色在线日韩免费| 亚洲成人久久爱视频| 亚洲国产精品合色在线| 亚洲,欧美精品.| 老司机福利观看| 天天躁夜夜躁狠狠躁躁| 99久久无色码亚洲精品果冻| av电影中文网址| 久久久水蜜桃国产精品网| 老司机福利观看| 日本黄色视频三级网站网址| 欧美性猛交╳xxx乱大交人| 黄色丝袜av网址大全| 大型黄色视频在线免费观看| 在线观看免费日韩欧美大片| 欧美乱码精品一区二区三区| 国产精品一区二区三区四区久久 | 精品日产1卡2卡| 不卡一级毛片| 国产亚洲精品久久久久5区| xxx96com| 丝袜在线中文字幕| 大型av网站在线播放| 窝窝影院91人妻| 91麻豆精品激情在线观看国产| 国产日本99.免费观看| 一进一出抽搐gif免费好疼| 麻豆av在线久日| 国产免费男女视频| 深夜精品福利| 精品久久久久久成人av| 亚洲成人免费电影在线观看| 国产伦人伦偷精品视频| 黄色毛片三级朝国网站| 国产主播在线观看一区二区| 人人妻人人澡人人看| 黄色成人免费大全| a级毛片a级免费在线| 久久久久精品国产欧美久久久| 亚洲欧美日韩无卡精品| 欧美成狂野欧美在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲人成网站高清观看| 成人三级黄色视频| 真人一进一出gif抽搐免费| 又紧又爽又黄一区二区| 男女之事视频高清在线观看| 午夜久久久在线观看| 午夜免费鲁丝| 精品久久蜜臀av无| 亚洲专区中文字幕在线| 十八禁人妻一区二区| netflix在线观看网站| 一区二区日韩欧美中文字幕| 国产亚洲欧美98| 亚洲成人久久性| 男女视频在线观看网站免费 | 两人在一起打扑克的视频| 中文在线观看免费www的网站 | 丝袜在线中文字幕| 免费搜索国产男女视频| 悠悠久久av| 高清毛片免费观看视频网站| 一本大道久久a久久精品| 亚洲av电影在线进入| 久久伊人香网站| 免费在线观看亚洲国产| 精品久久久久久久久久免费视频| 成人国产一区最新在线观看| 国产av又大| 一边摸一边做爽爽视频免费| a级毛片在线看网站| 国产aⅴ精品一区二区三区波| 热re99久久国产66热| 亚洲精品美女久久久久99蜜臀| 男女那种视频在线观看| 宅男免费午夜| 91麻豆av在线| 久久人人精品亚洲av| 好男人在线观看高清免费视频 | 国产午夜精品久久久久久| 高潮久久久久久久久久久不卡| 亚洲成av片中文字幕在线观看| 久久热在线av| 国产高清激情床上av| 中出人妻视频一区二区| 在线视频色国产色| 亚洲成国产人片在线观看| 国产三级在线视频| 日韩一卡2卡3卡4卡2021年| 此物有八面人人有两片| 国产精品美女特级片免费视频播放器 | 香蕉丝袜av| 国产伦一二天堂av在线观看| 视频区欧美日本亚洲| 人妻丰满熟妇av一区二区三区| 亚洲,欧美精品.| 黑人欧美特级aaaaaa片| 亚洲成人国产一区在线观看| √禁漫天堂资源中文www| 国产高清激情床上av| 校园春色视频在线观看| 亚洲国产精品999在线| 亚洲第一欧美日韩一区二区三区| 亚洲第一av免费看| 精品日产1卡2卡| 美女午夜性视频免费| 色av中文字幕| 国产精品香港三级国产av潘金莲| 久久久久国产一级毛片高清牌| 欧美日本亚洲视频在线播放| e午夜精品久久久久久久| 最好的美女福利视频网| 桃色一区二区三区在线观看| 午夜免费成人在线视频| e午夜精品久久久久久久| 亚洲真实伦在线观看| 国产视频内射| 日韩欧美国产一区二区入口| 日韩一卡2卡3卡4卡2021年| 国产爱豆传媒在线观看 | 国产亚洲欧美98| 亚洲国产毛片av蜜桃av| 99久久国产精品久久久| 此物有八面人人有两片| 国产国语露脸激情在线看| 亚洲 国产 在线| 久久九九热精品免费| 欧美三级亚洲精品| 99在线视频只有这里精品首页| 亚洲国产欧美日韩在线播放| 亚洲色图 男人天堂 中文字幕| 亚洲成人久久性| 99精品久久久久人妻精品| 男人舔女人的私密视频| 99热6这里只有精品| 亚洲成人久久性| 无人区码免费观看不卡| 97碰自拍视频| 欧美丝袜亚洲另类 | 日本熟妇午夜| 母亲3免费完整高清在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 一进一出好大好爽视频| 制服人妻中文乱码| 亚洲激情在线av| 精品一区二区三区av网在线观看| 他把我摸到了高潮在线观看| 黄片播放在线免费| 亚洲熟妇熟女久久| 99久久99久久久精品蜜桃| 免费观看人在逋| 婷婷亚洲欧美| 日韩av在线大香蕉| 女同久久另类99精品国产91| av视频在线观看入口| 国产精品1区2区在线观看.| 国产亚洲av嫩草精品影院| 一区二区三区高清视频在线| 中文字幕精品亚洲无线码一区 | 国产成人精品无人区| av免费在线观看网站| 成人手机av| 亚洲专区中文字幕在线| 十分钟在线观看高清视频www| 精华霜和精华液先用哪个| 欧美日韩亚洲国产一区二区在线观看| 精品国产美女av久久久久小说| 一本精品99久久精品77| 免费电影在线观看免费观看| 精品第一国产精品| 19禁男女啪啪无遮挡网站| 精品国产国语对白av| 精品国产乱码久久久久久男人| 国产精品国产高清国产av| 男女视频在线观看网站免费 | 日韩欧美 国产精品| 国产国语露脸激情在线看| 亚洲自拍偷在线| 欧洲精品卡2卡3卡4卡5卡区| 在线观看一区二区三区| 最近最新中文字幕大全电影3 | 香蕉丝袜av| 国产精品99久久99久久久不卡| 高清毛片免费观看视频网站| 欧美一区二区精品小视频在线| 夜夜夜夜夜久久久久| 一区二区日韩欧美中文字幕| 男女那种视频在线观看| 露出奶头的视频| 欧美绝顶高潮抽搐喷水| 视频在线观看一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 国产亚洲欧美98| 日韩欧美一区视频在线观看| 精品一区二区三区av网在线观看| 亚洲av电影在线进入| 久久久久国产精品人妻aⅴ院| 免费看十八禁软件| 午夜福利高清视频| 国产精品亚洲av一区麻豆| 一区福利在线观看| 欧美一级毛片孕妇| 日韩欧美三级三区| 真人做人爱边吃奶动态| 欧美日韩中文字幕国产精品一区二区三区| 国产精品久久电影中文字幕| 欧美av亚洲av综合av国产av| 伦理电影免费视频| 国产午夜精品久久久久久| 国产精品亚洲美女久久久| 国产成+人综合+亚洲专区| 久久久国产精品麻豆| 午夜日韩欧美国产| 搡老妇女老女人老熟妇| 一边摸一边做爽爽视频免费| 操出白浆在线播放| 亚洲第一av免费看| 国产精品一区二区免费欧美|