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

    氣體射流不穩(wěn)定性的建模和實(shí)驗(yàn)研究

    2016-11-12 06:41:14劉泛函王華徐建新
    化工進(jìn)展 2016年11期
    關(guān)鍵詞:李雅普諾夫氣液

    劉泛函,王華,徐建新

    (1省部共建復(fù)雜有色金屬資源清潔利用國家重點(diǎn)實(shí)驗(yàn)室,云南 昆明650093;2昆明理工大學(xué)冶金與能源工程學(xué)院,云南 昆明650093;3昆明理工大學(xué)質(zhì)量發(fā)展研究院,云南 昆明 650093)

    氣體射流不穩(wěn)定性的建模和實(shí)驗(yàn)研究

    劉泛函1,2,王華1,2,徐建新1,3

    (1省部共建復(fù)雜有色金屬資源清潔利用國家重點(diǎn)實(shí)驗(yàn)室,云南 昆明650093;2昆明理工大學(xué)冶金與能源工程學(xué)院,云南 昆明650093;3昆明理工大學(xué)質(zhì)量發(fā)展研究院,云南 昆明 650093)

    基于圖像處理技術(shù)和混沌理論,提出了一種氣體射流圖像的量化方法,用以反映氣液混合性能。本文使用大津法來獲得自適應(yīng)閥值,并得到最佳的圖像分割效果?;旌闲阅埽∕)用以表征氣體射流的空間分布,其定義為氣體射流像素占比(W)與射流區(qū)域內(nèi)的灰度標(biāo)準(zhǔn)差(N)之比。通過混合性能的非線性時間序列來表征射流噴射過程的不穩(wěn)定變化,并得到了以下結(jié)果:在射流區(qū)制下,射流的穩(wěn)定性隨著修正弗魯?shù)聰?shù)的增加而增加,證明了所提出方法的可靠性;其次,在同一工況中射流的穩(wěn)定性與最大李雅普諾夫指數(shù)之間具有明顯的線性關(guān)系,其線性相關(guān)系數(shù)為0.954,并構(gòu)建了射流的穩(wěn)定性與最大李雅普諾夫指數(shù)的線性關(guān)系模型。證明了最大李雅普諾夫指數(shù)不僅可以判斷混合過程是否呈現(xiàn)混沌狀態(tài),而且其數(shù)值的大小反應(yīng)了氣體射流穩(wěn)定性的強(qiáng)弱。

    多相流;混合;穩(wěn)定性;氣泡;混沌

    在鋼鐵和有色冶金領(lǐng)域內(nèi),為了滿足工業(yè)生產(chǎn)中的各種要求,高速氣體噴射到熔融金屬熔池內(nèi)的技術(shù)被廣泛地應(yīng)用。如轉(zhuǎn)爐煉鋼和冰銅的連續(xù)吹煉等過程都是利用高速氧氣射流噴入金屬熔池內(nèi),使得反應(yīng)動力學(xué)條件改善,可以在較短的時間內(nèi)完成各種冶煉反應(yīng),同時熔池中的氣體流動行為對于冶金反應(yīng)、氣體和夾雜物的排出以及熔池?cái)嚢瓒加泻艽蟮挠绊?。人們發(fā)現(xiàn)在熔池熔煉過程中,無論采用向上底吹還是水平側(cè)吹,在噴射冶金的過程中都存在著一個難以解決的問題,即風(fēng)口或噴槍周圍的耐火材料蝕損十分嚴(yán)重,極大地降低了爐襯的使用壽命,氣體射流的回?fù)衄F(xiàn)象被認(rèn)為是造成耐火材料腐蝕的一個重要原因。

    AOKI等[1]發(fā)現(xiàn)了氣體射流的回?fù)衄F(xiàn)象并對將其定義為“back-attack”,回?fù)羰撬職怏w射流間歇性的出現(xiàn)氣體攜帶液體反向運(yùn)動撞擊噴嘴的現(xiàn)象,這也解釋了耐火材料的腐蝕。YANG等[2-3]對氣體射流的不穩(wěn)定性進(jìn)行拍攝研究,并提出了一個“空穴模型”的理論用來解釋風(fēng)口腐蝕的原因。MORI等[4]研究了水下氣體射流的動力學(xué)特性,他們發(fā)現(xiàn)隨著噴嘴處氣體噴射速度的增加射流的回?fù)纛l率逐漸減小。WEILAND等[5]通過高速攝像技術(shù)對水下氣體射流氣液界面的穩(wěn)定性進(jìn)行了研究。戚隆溪等[6]使用探針測量了射流軸線的靜壓分布,結(jié)果表明在欠膨脹工況下,流場內(nèi)存在著復(fù)雜的膨脹壓縮波系區(qū)域。王柏懿、DAI等[7-8]拍攝了回?fù)衄F(xiàn)象的演化過程,完整的回?fù)暨^程包括了頸縮、脹鼓、分離和回?fù)?,并指出回?fù)衄F(xiàn)象產(chǎn)生的原因是水下氣體射流剪切層的不穩(wěn)定性。王曉剛等[9]研究了高速氣體射流的形貌,發(fā)現(xiàn)了脹鼓現(xiàn)象,并指出了射流回?fù)糁靶枰?jīng)歷多次脹鼓,且回?fù)舻膹?qiáng)度隨著脹鼓次數(shù)和幅度的增加而增大。郭強(qiáng)等[10-11]發(fā)現(xiàn)了一種特殊的回?fù)衄F(xiàn)象,它包含有頸縮、脹鼓和分離而不發(fā)生回?fù)?,認(rèn)為回?fù)舻念l率隨著馬赫數(shù)的增加而降低。施紅輝等[12-13]對脹鼓與回?fù)糁g的關(guān)聯(lián)性進(jìn)行了研究,認(rèn)為脹鼓和回?fù)舳际怯蓺怏w射流內(nèi)部的壓力震蕩所引發(fā)的,脹鼓是能量積聚的表現(xiàn),積聚到一定程度便引發(fā)回?fù)簟_@些研究成果對減弱或消除耐火材料的腐蝕具有重要的意義。

    高速氣體射流是一個復(fù)雜的非線性運(yùn)動,在噴射過程中其內(nèi)部的激波和膨脹波將會引起流形振蕩。目前,對于這類非線性運(yùn)動人們一般采用混沌理論對其進(jìn)行研究。寧偉征[14]利用偏心射流攪拌和機(jī)械攪拌的共同作用,使得攪拌槽內(nèi)的流場形成混沌混合,并通過宏觀不穩(wěn)定頻率、分形維數(shù)和最大李雅普諾夫指數(shù)作為混沌混合的判據(jù)對混合過程進(jìn)行研究,提高了攪拌槽的混合效率。徐建新[15]提出了一種多相流混合效果貝蒂數(shù)評價(jià)方法,對多相混合的均勻性以及非均勻性進(jìn)行了有效量化,并使用最大李雅普諾夫指數(shù)在內(nèi)的多個指標(biāo)對多相混合效果的混沌特征進(jìn)行了分析。劉政等[16]基于流動軌跡的Kolmogorov熵和分形維數(shù)對鋁合金熔體中混沌對流進(jìn)行表征,研究了電磁場作用下混沌對流對半固態(tài)A356鋁合金熔體形貌演變的影響。GOTODA等[17-18]基于非線性動力學(xué)研究了在旋流作用下的火焰振蕩,他發(fā)現(xiàn)隨著旋流速度的增加,火焰從周期震蕩經(jīng)歷準(zhǔn)周期震蕩向混沌轉(zhuǎn)變,也表明所采用的關(guān)聯(lián)維數(shù)法可以有效的量化復(fù)雜火焰的不穩(wěn)定運(yùn)動。綜上所述,混沌是這類非線性動力學(xué)系統(tǒng)所特有的一種運(yùn)動形式,混沌特征對系統(tǒng)性能有較大的影響,李雅普諾夫指數(shù)法作為其中最經(jīng)典的一種判別方法,在混沌特征分析中有著廣泛的應(yīng)用。

    目前,對高速氣體射流的研究中,還沒有一個科學(xué)可靠的量化方法。前人所采用統(tǒng)計(jì)脹鼓和回?fù)纛l率的研究方法容易出現(xiàn)誤判,同時這種方法也不能有效的區(qū)分出每次脹鼓和回?fù)衄F(xiàn)象的所造成不穩(wěn)定程度的不同。本文基于圖像灰度統(tǒng)計(jì)提出了一種氣液混合性能的量化方法,通過分析混合性能的時間序列,來量化氣體射流不穩(wěn)定性。結(jié)果表明,該方法對氣體射流的穩(wěn)定性提出了量化標(biāo)準(zhǔn)并能夠客觀、準(zhǔn)確地進(jìn)行評價(jià),可以為氣液混合過程的分析及優(yōu)化提供一定的指導(dǎo),具有良好的適用性。

    1 實(shí)驗(yàn)裝置

    高溫熔渣具有很強(qiáng)的腐蝕性,使得對熔池內(nèi)部的研究存在著很大的困難,因此本文通過水力學(xué)模型對側(cè)吹氣體射流的行為進(jìn)行研究。實(shí)驗(yàn)裝置包括一個瓦紐科夫爐水力學(xué)模型、供氣系統(tǒng)以及數(shù)據(jù)/圖像采集系統(tǒng),如圖1所示。瓦紐科夫爐水力學(xué)模型按照實(shí)際瓦紐科夫爐1∶10的比例制作,具體結(jié)構(gòu)尺寸如圖2所示。在實(shí)驗(yàn)過程中氣體通過固定在爐體兩側(cè)的風(fēng)口處的L形噴槍對爐內(nèi)的熔池進(jìn)行噴吹攪拌,噴槍由帶孔膠塞固定在風(fēng)口處。水力學(xué)模型內(nèi)部液體為水,液面深度為20cm,噴嘴距液面處的垂直距為5cm。小型空氣壓縮機(jī)通過PVC透明管與玻璃轉(zhuǎn)子流量計(jì)相連接,玻璃轉(zhuǎn)子流量計(jì)再通過PVC透明管連接到噴嘴上,其中玻璃轉(zhuǎn)子流量計(jì)的型號為LZB-15,其測量范圍是0.25~2.5m3/h,小型空氣壓縮機(jī)的型號為羅威Z-0.12/8,其額定功率為2.5HP(1HP=0.7457kW),額定壓力為0.8MPa,容積流量為120L/min,并帶有一個24L、具有穩(wěn)定壓力作用的儲氣罐。高速攝像儀放置在水力模型的后側(cè)進(jìn)行拍攝,水力模型的另一側(cè)(渣室側(cè))貼有黑紙可以最大限度地減少外部光線的影響,有助于提高所拍攝照片的對比度,高速攝像儀的型號為德國PCO.dimax HD,在全分辨率1920×1080像素下的最高攝像速度為2128幀/s下可捕獲到超高清晰的圖像,1008×1000像素時為3822幀/s,最短快門時間是1.5μs,其自帶內(nèi)存為36Gb??紤]到多組噴槍共同工作時,射流之間不僅會相互影響同時也會相互遮擋,使得后續(xù)的量化處理存在較大的困難,因此實(shí)驗(yàn)過程中只研究一組相對噴槍的混合攪拌性能,其他不使用的風(fēng)口通過膠塞進(jìn)行封堵,實(shí)驗(yàn)具體工況如表1所示。

    圖1 實(shí)驗(yàn)裝置

    圖2 瓦紐科夫爐水力學(xué)模型結(jié)構(gòu)尺寸(單位:m)

    2 圖像分析

    2.1圖像分割

    在一幅完整的彩色圖像中,每種顏色都是由紅色(R)、綠色(G)和藍(lán)色(B)3種原色構(gòu)成。每個原色的亮度范圍都是0~255,其數(shù)值越大,色彩越明亮,每個像素點(diǎn)所顯示出來的顏色都是3種原色亮度的疊加,RGB顏色模型的空間示意圖如圖3所示。而在實(shí)驗(yàn)過程中,高速攝像儀拍攝所得到的是灰度數(shù)字圖,圖像中每個像素都顯示為從最暗黑色到最亮的白色的灰度,其亮度范圍是0~255。從理論上來講,任意一種顏色的不同深淺,以及處于不同亮度上的不同顏色都可以轉(zhuǎn)換成灰度。根據(jù)人眼對三種原色敏感程度的不同,以不同的權(quán)值進(jìn)行加權(quán)平均就能得到合理的灰度圖像,見式(1)。

    圖3 RGB顏色模型的空間示意圖

    圖像分割是從圖像處理到圖像分析過程中的一個重要步驟,由于圖片中氣體比背景要亮得多,那么通過閥值就可以簡單、快速、準(zhǔn)確地對二者進(jìn)行分割。不同的分割算法得到不同的閾值使得分割圖像的結(jié)果具有很大的差別,本文是通過大津法獲得的最佳閾值。大津法是基于聚類的思想,假設(shè)圖象包含L個灰度級(0,1,…,L-1),從L個灰度級遍歷閥值T,當(dāng)T為某個值的時候,使得前景和背景的類間方差最大,這也就意味著錯分的可能性最小,此時這個T值便是所要求得的最佳閾值,方差的計(jì)算公式如式(2)、式(3)。

    式中,T為設(shè)定的閾值;w0為前景像素點(diǎn)的占比;w1為背景像素點(diǎn)的占比;u0為前景像素點(diǎn)的平均灰度;u1為背景像素點(diǎn)的平均灰度;u為圖像總的平均灰度。

    2.2混合性能指數(shù)計(jì)算

    氣體射流形狀的改變被高速攝像儀以數(shù)字的形式進(jìn)行記錄,所有圖像都實(shí)時記錄并保存在計(jì)算機(jī)中,這給后續(xù)的分析工作提供了豐富的數(shù)據(jù)??梢岳靡曨l圖像中豐富的信息,對其進(jìn)行編碼并提取特征量,這些特征量的簡單直觀,便于進(jìn)行相關(guān)的處理,圖像分析的關(guān)鍵步驟如下所述。

    步驟1:定義圖像的工作區(qū)域,只截取被研究的區(qū)域,對結(jié)果有干擾的區(qū)域?qū)粍h除,值得注意的是,一組連續(xù)工況中所有圖像所定義的工作區(qū)域,其位置和大小是相同的,使得所有圖像可以在同一體系中進(jìn)行分析和比較。

    步驟2:最佳閾值T是借助于Matlab軟件通過大津法進(jìn)行計(jì)算所獲得的,并對圖像進(jìn)行進(jìn)一步的分割處理(Pjet>T),射流區(qū)域的像素占比用于表征射流混合反應(yīng)區(qū)域的大小,它的表達(dá)式為式(4)。

    步驟3:氣液混合是一個復(fù)雜的隨機(jī)過程,人們逐漸的意識到在氣體射流與熔液的混合反應(yīng)過程中,氣液接觸面積對于促進(jìn)冶煉反應(yīng)的重要性。如圖4所示氣體射流內(nèi)的氣體將會出現(xiàn)以下3種情況,即大氣泡、小氣泡和氣液摻混,氣液摻混可以極大地增加氣液間接觸面積,從而加快氣液間反應(yīng)。

    對分割區(qū)域內(nèi)像素點(diǎn)的亮度特征進(jìn)行分析,發(fā)現(xiàn)氣泡大小的不同對分割區(qū)域內(nèi)像素灰度值的均勻性具有很大的影響,大氣泡的灰度標(biāo)準(zhǔn)差(N)最大,小氣泡次之,而氣液摻混最小,那么就通過分割區(qū)域內(nèi)像素灰度值的標(biāo)準(zhǔn)差來表征氣體分布的均勻性,即標(biāo)準(zhǔn)差越小,氣液混合越均勻,那么氣液間的反應(yīng)就越快。

    步驟4:混合性能指數(shù)可以通過計(jì)算在分割區(qū)域中像素灰度值的標(biāo)準(zhǔn)差得到,其表達(dá)式如式(5)。

    基于整個混合過程的采樣圖像就可以獲得混合性能指數(shù)的時間序列M(t)(t=0,1,2,…,k),通過混合性能指數(shù)的平均值來衡量整個混合過程的混合性能,混合性能指數(shù)的平均值表達(dá)式如式(6)。

    圖4 氣體射流中不同大小的氣泡

    3 結(jié)果與討論

    3.1氣體射流的穩(wěn)定性

    在高速氣體射流的拍攝過程中,應(yīng)盡量減小分辨率以增加拍攝的持續(xù)時間,同時曝光時間不應(yīng)過長,否則照片容易模糊。通過高速攝像儀自帶的軟件對其參數(shù)設(shè)置如下,拍攝速度為1000幀/s,分辨率為720×400,曝光時間為800μs。圖5展示了亞音速工況下高速氣體射流流態(tài)的演化過程。

    在0.001~0.006s可清晰地觀察到射流的脹鼓形態(tài),噴嘴出口附近產(chǎn)生某種小幅度膨脹,這種膨脹的瞬間發(fā)生,然后會被調(diào)整射流迅速帶離并消失,整個過程中射流不發(fā)生破碎。其發(fā)生機(jī)制是在雷諾數(shù)較高且密度差較大的條件下,混合剪切層的不穩(wěn)定性使得噴嘴出口處的射流發(fā)生間歇性頸縮,射流通道不暢從而造成氣體聚集,在這個過程中射流內(nèi)部的壓力急劇增大,當(dāng)壓力達(dá)到一定程度后就會導(dǎo)致射流脹鼓,脹鼓發(fā)生的演化過程通過示意圖進(jìn)行展示。

    在0.015~0.023s可觀測到典型的水下氣體射流回?fù)衄F(xiàn)象的演化過程,射流經(jīng)過多次的脹鼓后迅速膨脹,由于液體慣性的約束使得使膨脹主要發(fā)生在軸向方向上。當(dāng)膨脹達(dá)到一定程度時,氣體射流斷裂成為兩部分,其中一部分繼續(xù)向前運(yùn)動,而另一部分則向后運(yùn)動并撞擊噴嘴。其運(yùn)動過程通過示意圖進(jìn)行展示,后者在撞擊噴嘴的過程中破碎成大量的氣泡并上浮。

    在瓦紐科夫爐水力模型雙側(cè)吹試驗(yàn)過程中,由于氣體射流發(fā)生間歇性的脹鼓與回?fù)簦淞鞯膭恿l件不斷發(fā)生變化,使得射流的混合區(qū)域發(fā)生大幅度振蕩,射流內(nèi)部的氣體相互碰撞并合并形成大氣泡,大氣泡的形成減小了氣液間的接觸面積,造成混合過程中混合性能的惡化,如圖6所示。大氣泡迅速上浮又對整個氣體射流的流場及形態(tài)變化造成一定的影響。因此氣體射流的不穩(wěn)定性與射流的形態(tài)變化具有關(guān)聯(lián)性,量化射流的形態(tài)特征用以反映射流的混合性能對射流穩(wěn)定性的研究具有重要的意義。

    3.2不同修正弗魯?shù)聰?shù)下的混合過程

    在氣液混合實(shí)驗(yàn)的過程中,采樣圖像之間的間隔時間是影響量化的另一個重要因素,如果間隔時間過長就會丟失混合過程的連貫性信息,如果間隔時間過短則會花費(fèi)大量的計(jì)算時間。綜合考慮結(jié)果的準(zhǔn)確度和計(jì)算量之間的關(guān)系,將采樣時間間隔設(shè)置為0.1s,圖7顯示不同的修正弗魯?shù)聰?shù)下混合性能的時間序列,其計(jì)算結(jié)果如表1所示。其中,D0為噴嘴直徑,mm;Q為氣體體積流量,m3/h;Fr′為修正的弗魯?shù)聰?shù);σ為時間序列的標(biāo)準(zhǔn)差;V為比σ,用來評價(jià)氣體射流的不穩(wěn)定性,1/V用來評價(jià)氣體射流的穩(wěn)定性。

    圖5 亞音速工況下高速氣體射流流態(tài)的演化過程

    圖6 大氣泡對射流形態(tài)的影響

    圖7 混合性能指數(shù)的時間序列

    施紅輝等[19]用“氣袋”模型形象的對氣體射流的不穩(wěn)定性進(jìn)行了解釋,他認(rèn)為氣體射流是由不可壓縮的流體包裹著一個可壓縮的“氣袋”,由此產(chǎn)生非線性的震蕩壓力,使得射流的氣液邊界不穩(wěn)定。在亞音速工況下,“氣袋”的長度隨修正弗魯?shù)聰?shù)的增加有著較為明顯的增長,“氣袋”越長其內(nèi)部的能量越難以積聚,脹鼓和回?fù)衄F(xiàn)象也就越難以發(fā)生,因此隨著修正弗魯?shù)聰?shù)的增加,氣體射流越穩(wěn)定,這與本文所得到的實(shí)驗(yàn)結(jié)果是相吻合的。

    表1 不同修正弗魯?shù)聰?shù)下氣體射流的混合性能和穩(wěn)定性

    如圖8和圖9所示,在修正弗魯?shù)聰?shù)較大的射流區(qū)制下,射流的穩(wěn)定性隨著修正弗魯?shù)聰?shù)的增加而增加。但在修正弗魯?shù)聰?shù)較小的泡流區(qū)制下,修正弗魯?shù)聰?shù)并不能完全解釋混合性能以及穩(wěn)定性的變化,此時氣體體積流量對混合性能以及穩(wěn)定性所造成的影響同樣十分明顯。從總體來看,隨著修正弗魯?shù)聰?shù)的增加,射流穩(wěn)定性的總趨勢是增加的,這說明修正弗魯?shù)聰?shù)是影響射流穩(wěn)定性變化的主要因素。L3工況的修正弗魯?shù)聰?shù)最大,同時也是實(shí)驗(yàn)中混合反應(yīng)性能和穩(wěn)定性的最好的工況。

    相關(guān)系數(shù)是用于測量在多元統(tǒng)計(jì)中兩組變量之間線性度的一個量綱為1指數(shù),它包括正相關(guān)、不相關(guān)和負(fù)相關(guān)。通過一定的計(jì)算,得到相關(guān)系數(shù)取值范圍為[0,1]時,說明兩組數(shù)據(jù)呈正相關(guān),且取值越大相關(guān)性越強(qiáng),相關(guān)系數(shù)的計(jì)算公式為式(7)~式(9)。

    其中

    圖8 混合性能變化趨勢

    圖9 穩(wěn)定性變化趨勢

    通過對圖10的觀察,可以明顯地發(fā)現(xiàn)混合性能指數(shù)和修正弗魯?shù)聰?shù)具有相似的趨勢,認(rèn)為混合性能指數(shù)和修正弗魯?shù)聰?shù)之間具有線性關(guān)系,通過最小二乘法擬合出參數(shù),其擬合方程為:=5× 10-6Fr′+0.003,線性相關(guān)系數(shù)為0.979。混合性能指數(shù)與修正弗魯?shù)聰?shù)的變化趨勢具有良好的線性關(guān)系,說明混合性能指數(shù)可以有效的量化側(cè)吹射流的噴吹混合過程,證明了該方法的可靠性。

    3.3最大李雅普諾夫指數(shù)

    為了進(jìn)一步定量研究不同工況下射流不穩(wěn)定性的差異,通過非線性時間序列分析手段對混合性能進(jìn)行了混沌特征分析,本文采用李雅普諾夫指數(shù)法來進(jìn)行混沌特征的表征。李雅普諾夫指數(shù)反映了相鄰軌跡的分離率,它的重要功能之一是估算系統(tǒng)的混沌行為。李雅普諾夫指數(shù)始終是一個實(shí)數(shù),可以是正數(shù)、負(fù)數(shù)或零。最大李雅普諾夫指數(shù)(LLE)為正則顯示了系統(tǒng)是混沌的,混沌系統(tǒng)對初始條件十分敏感,在初期略微的改變都會使得后期發(fā)生重大的變化。最大李雅普諾夫指數(shù)為負(fù)則顯示了系統(tǒng)是規(guī)則的,說明系統(tǒng)較為穩(wěn)定,對初始條件不敏感。而最大李雅普諾夫指數(shù)為零時說明系統(tǒng)處于臨界狀態(tài)。

    通過Visual Recurrence Analysis軟件使用虛假鄰點(diǎn)法和平均互信息法計(jì)算時間序列的嵌入維數(shù)(d)和延遲時間(t),并對相空間進(jìn)行重新構(gòu)建。圖11所示為利用最大李雅普諾夫指數(shù)法[20]判斷混合性能指數(shù)的序列是否達(dá)到混沌狀態(tài),如果最大李雅普諾夫指數(shù)等于零,則系統(tǒng)處于臨界狀態(tài),此時最大李雅普諾夫指數(shù)無法判斷系統(tǒng)是否處于混沌狀態(tài),就需要利用0-1test方法[21-23]來進(jìn)一步驗(yàn)證。不同工況下最大李雅普諾夫指數(shù)的計(jì)算結(jié)果如表2所示,表明在所有的工況下的混合過程都是混沌的,隨著修正弗魯?shù)聰?shù)的增加,系統(tǒng)變得越來越復(fù)雜,系統(tǒng)的混沌程度增大。

    圖11 時間序列曲線中獲取的李雅普諾夫指數(shù)(LLE)

    表2 不同修正弗魯?shù)聰?shù)下氣體射流的最大李雅普諾夫指數(shù)和穩(wěn)定性

    通過觀察圖12,可以發(fā)現(xiàn)穩(wěn)定性與最大李雅普諾夫指數(shù)之間同樣具有線性關(guān)系。通過最小二乘法進(jìn)行擬合并構(gòu)建了穩(wěn)定性和最大李雅普諾夫指數(shù)之間的線性模型,穩(wěn)定性和最大李雅普諾夫指數(shù)之間的關(guān)系為LLE=62.726/V-0.113,線性相關(guān)系數(shù)為0.954。證明了最大李雅普諾夫指數(shù)不僅可以判斷混合過程是否呈現(xiàn)混沌狀態(tài),而且其數(shù)值的大小反應(yīng)了混合過程穩(wěn)定性的強(qiáng)弱。

    圖12 最大李雅普諾夫指數(shù)(LLE)與穩(wěn)定性之間的線性關(guān)系

    4 結(jié) 論

    本文對瓦紐科夫爐水力學(xué)模型內(nèi)側(cè)吹氣體的行為進(jìn)行了實(shí)驗(yàn)研究,利用高速攝像儀對不同的修正弗魯?shù)聰?shù)下的氣體射流形態(tài)變化進(jìn)行了記錄?;趫D像處理技術(shù)和混沌理論,提出了一種新的方法來量化所采集的圖像中包含的氣體射流非線性運(yùn)動。該方法的主要步驟包括定義工作區(qū)、分割圖像、提取特征參數(shù)和計(jì)算混合性能指數(shù)。其中,混合性能指數(shù)被定義為氣體像素占比與氣體像素灰度標(biāo)準(zhǔn)差之比,用以表征氣體射流在空間上的分布,該方法的可靠性通過比較混合性能指數(shù)與修正弗魯?shù)聰?shù)的變化趨勢進(jìn)行驗(yàn)證。分析了混合性能指數(shù)的時間序列,結(jié)果表明氣體射流的穩(wěn)定性隨著修正弗魯?shù)聰?shù)的增大而增加,證明了最大李雅普諾夫指數(shù)數(shù)值的大小同樣也可以反映混合過程穩(wěn)定性的強(qiáng)弱。該方法在研究氣體射流過程中可以為工業(yè)設(shè)計(jì)冶金熔煉爐以及其他多相混合性能評價(jià)體系提供一種新的借鑒。

    [1] AOKI T,MASUDA S,HATANO A,et al. Characteristics of submerged gas jets and a new type bottom blowing tuyere[C]// Injection Phenomena in Extraction and Refining,Newcastle upon Tyne,England,1982.

    [2] YANG Q X,GUSTAVSSON H,BURSTROM E. Erosion of refractory during gas injection. A cavitation based model[J]. Scandinavian Journal of Metallurgy,1990,19:127-136.

    [3] YANG Q X,GUSTAVSSON H. Effects of gas jet instability on refractory wear[J]. Scandinavian Journal of Metallurgy,1992,21:15-26.

    [4] MORI K,OZAWA Y,SANO M. Characterization of gas jet behavior at a submerged orifice in liquid metal.[J]. Isij International,1982,22:377-384.

    [5] WEILAND C,YAGLA J,VLACHOS P. Experimental study of the stability of a high-speed gas jet under the influence of liquid cross-flow[C]// ASME/JSME 2007 5th Joint Fluids Engineering Conference,2007:599-610.

    [6] 戚隆溪,曹勇,王柏懿. 水下欠膨脹高速氣體射流的實(shí)驗(yàn)研究[J].力學(xué)學(xué)報(bào),2000,32(6):667-675.

    [7] 王柏懿,戴振卿,戚隆溪,等. 水下超聲速氣體射流回?fù)衄F(xiàn)象的實(shí)驗(yàn)研究[J]. 力學(xué)學(xué)報(bào),2007,39(2):267-272.

    [8] DAI Z Q,WANG B Y,QI L X,et al. Experimental study on hydrodynamic behaviors of high-speed gas jets in still water[J]. Acta Mechanica Sinica,2006,22(5):443-448.

    [9] 王曉剛,王超,郭強(qiáng),等. 二維水槽中高速氣體射流的振蕩流流型研究[J]. 浙江理工大學(xué)學(xué)報(bào),2009,26(4):613-618.

    [10] 郭強(qiáng),施紅輝,王超,等. 準(zhǔn)二維水下超聲速垂直過膨脹射流研究[J]. 實(shí)驗(yàn)流體力學(xué),2010,24(3):6-18.

    [11] 郭強(qiáng),施紅輝,王超,等. 水下超聲速氣體射流氣液兩相復(fù)雜流動研究[J]. 工程熱物理學(xué)報(bào),2012,33(5):809-812.

    [12] SHI H H,GUO Q,WANG C,et al. Oscillation flow induced by underwater supersonic gas jets[J]. Shock Waves,2010,20:47-352.

    [13] 施紅輝,郭強(qiáng),王超,等. 水下超音速氣體射流脹鼓和回?fù)舻年P(guān)聯(lián)性研究[J]. 力學(xué)學(xué)報(bào),2010,42(6):1206-1210.

    [14] 寧偉征. 偏心射流強(qiáng)化流體混沌混合的研究[D]. 重慶:重慶大學(xué),2012.

    [15] 徐建新. 多相體系攪拌混合效果評價(jià)方法及其應(yīng)用研究[D]. 昆明:昆明理工大學(xué),2012.

    [16] 劉政,張嘉藝,羅浩林,等. 混沌對流下的半固態(tài)A356鋁合金初生相形貌演變研究[J]. 金屬學(xué)報(bào),2015,52(2):177-183.

    [17] GOTODA H,ASANO Y,CHUAH K H,et al. Nonlinear analysis on dynamic behavior of buoyancy-induced flame oscillation under swirling flow[J]. International Journal of Heat and Mass Transfer,2009,52(23):5423-5432.

    [18] GOTOD H,UEDA T. Transition from periodic to non-periodic motion of a bunsen-type premixed flame tip with burner rotation[J]. Proceedings of the Combustion Institute,2002,29(2):1503-1509.

    [19] 施紅輝,王柏懿,戴振卿. 水下超聲速氣體射流的力學(xué)機(jī)制研究[J].中國科學(xué)(物理學(xué)力學(xué)天文學(xué)),2010,40(1):92-100.

    [20] WOLF A,SWIFT J B,SWINNEY H L,et al. Determining Lyapunov exponents from a time series[J]. Physica D Nonlinear Phenomena,1985,16(3):285-317.

    [21] GOTTWALD G A,MELBOURNE I. Comment on “Reliability of the 0-1 test for chaos”[J]. Physical Review E,2008,77(2):028201.

    [22] FALCONER I,GOTTWALD G A,MELBOURNE I,et al. Application of the 0-1 test for chaos to experimental data[J]. SIAM Journal on Applied Dynamical Systems,2007,6(2):395-402.

    [23] ASCANI F,KEELER P,KUBIAK R,et al. Detection of low-dimensional chaos in quasi-periodic time series:The 0-1 test[R]. Technical Report,Santa Fe Institute Complex Systems Summer School,2008.

    Experiments and modelling of gas jet instability using image processing and chaos theory

    LIU Fanhan1,2,WANG Hua1,2,XU Jianxin1,3
    (1State Key Laboratory of Complex Nonferrous Metal Resources Clean Utilization,Kunming 650093,Yunnan,China;2Faculty of Metallurgical and Energy Engineering,Kunming University of Science and Technology,Kunming 650093,Yunnan,China;3Quality Development Institure,Kunming University of Science and Technology,Kunming 650093,Yunnan,China)

    Based on image processing technology and chaos theory, a novel method for quantitative characterization of gas jet images was proposed,which reflects the mixing performance. The Otsu method was used to find the best threshold value and hence the best segmented images. The mixing performance(M) is defined as the ratio of proportion of gas jet pixels(W) to standard deviation of the grayscale values in gas jet region(N) for characterizing spatial distribution of gas jet. Unsteady motions of the gas jet were characterized by the nonlinear time series of the mixing performances. As have shown,the stability of gas jet(1/V) decrease with the increasing modified Froude number under jetting regime,which confirms the validity of our approach. On top of that,the linear relation between stability and largest Lyapunov exponent(LLE) was obvious with a correlation coefficient of 0.954. For further research a linear model of stability and LLE were constructed. The results showed that the LLE can not only judge the mixing process is nonchaotic or chaotic,also the numerical size of LLE reflected the stability of gas jet.

    multiphase flow;mixing;stability;bubble;chaos

    TF 02

    A

    1000-6613(2016)11-3433-08

    10.16085/j.issn.1000-6613.2016.11.007

    2016-04-11;修改稿日期:2016-05-24。

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2014CB460605)及國家自然科學(xué)基金(51174105,51206071,51366005,51406071,51206071)項(xiàng)目。

    劉泛函(1986—),男,博士研究生。聯(lián)系人:徐建新,副教授。E-mail xujianxina@163.com。

    猜你喜歡
    李雅普諾夫氣液
    李雅普諾夫:彼得堡數(shù)學(xué)學(xué)派的健將
    基于增廣Lyapunov 泛函的時變時滯T-S模糊系統(tǒng)穩(wěn)定性分析
    微重力下兩相控溫型儲液器內(nèi)氣液界面仿真分析
    氣液分離罐液位計(jì)接管泄漏分析
    系統(tǒng)H∞范數(shù)計(jì)算:Lyapunov函數(shù)的直接優(yōu)化方法
    CO2 驅(qū)低液量高氣液比井下氣錨模擬與優(yōu)化
    Novoselova瑪莎·諾夫斯洛娃【英國】
    博客天下(2014年24期)2014-09-25 02:47:29
    采用李雅普諾夫函數(shù)的電液伺服系統(tǒng)反饋線性化控制
    氣液接觸法制備球形氫氧化鎂
    幽默與漫畫
    海外文摘(2009年11期)2009-12-31 08:33:26
    婷婷成人精品国产| 窝窝影院91人妻| 亚洲成国产人片在线观看| 亚洲激情五月婷婷啪啪| av网站免费在线观看视频| 俄罗斯特黄特色一大片| 亚洲精品中文字幕一二三四区 | 宅男免费午夜| 亚洲视频免费观看视频| 精品人妻熟女毛片av久久网站| 成人免费观看视频高清| 免费女性裸体啪啪无遮挡网站| 免费观看人在逋| 亚洲七黄色美女视频| 免费人妻精品一区二区三区视频| 色精品久久人妻99蜜桃| 日本撒尿小便嘘嘘汇集6| 日本vs欧美在线观看视频| 各种免费的搞黄视频| 美女扒开内裤让男人捅视频| 男人爽女人下面视频在线观看| 日本91视频免费播放| 久久女婷五月综合色啪小说| 性少妇av在线| 亚洲中文日韩欧美视频| 男人舔女人的私密视频| 视频区欧美日本亚洲| 久久久国产成人免费| 精品一区二区三区av网在线观看 | 在线观看www视频免费| 80岁老熟妇乱子伦牲交| 18在线观看网站| 法律面前人人平等表现在哪些方面 | 国产淫语在线视频| 国内毛片毛片毛片毛片毛片| 色94色欧美一区二区| 9热在线视频观看99| 欧美日韩国产mv在线观看视频| 亚洲欧美清纯卡通| 黄片播放在线免费| 自拍欧美九色日韩亚洲蝌蚪91| 99精品久久久久人妻精品| 欧美人与性动交α欧美精品济南到| 久久这里只有精品19| 美女扒开内裤让男人捅视频| 精品少妇内射三级| 女性被躁到高潮视频| 一边摸一边做爽爽视频免费| 午夜精品国产一区二区电影| 国产一区有黄有色的免费视频| 久久天堂一区二区三区四区| 国产在线免费精品| 男女下面插进去视频免费观看| av超薄肉色丝袜交足视频| 国产极品粉嫩免费观看在线| 2018国产大陆天天弄谢| 精品一区二区三区四区五区乱码| 国产精品一二三区在线看| 精品卡一卡二卡四卡免费| 国产成人精品久久二区二区免费| 黄色毛片三级朝国网站| 又黄又粗又硬又大视频| 狂野欧美激情性bbbbbb| 国产成人欧美| 精品视频人人做人人爽| 大香蕉久久网| 亚洲av欧美aⅴ国产| 久久久国产成人免费| 亚洲综合色网址| 亚洲五月色婷婷综合| 在线 av 中文字幕| 丰满饥渴人妻一区二区三| 午夜福利乱码中文字幕| 久久性视频一级片| 啦啦啦在线免费观看视频4| 18禁观看日本| 一个人免费在线观看的高清视频 | 国产免费视频播放在线视频| 日本91视频免费播放| 成年美女黄网站色视频大全免费| 欧美大码av| 亚洲三区欧美一区| 亚洲国产日韩一区二区| e午夜精品久久久久久久| 亚洲欧美一区二区三区久久| 欧美在线黄色| a在线观看视频网站| 久久久久久久久免费视频了| 99国产极品粉嫩在线观看| 黑丝袜美女国产一区| 一本色道久久久久久精品综合| 无限看片的www在线观看| 亚洲国产毛片av蜜桃av| 久久免费观看电影| 日日摸夜夜添夜夜添小说| 亚洲全国av大片| 久久精品人人爽人人爽视色| 久久久久视频综合| 侵犯人妻中文字幕一二三四区| 久久av网站| 如日韩欧美国产精品一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 免费在线观看影片大全网站| 大香蕉久久网| 久久精品亚洲av国产电影网| 99精品欧美一区二区三区四区| 欧美 日韩 精品 国产| 两性午夜刺激爽爽歪歪视频在线观看 | 在线永久观看黄色视频| 19禁男女啪啪无遮挡网站| 永久免费av网站大全| 亚洲人成电影免费在线| 熟女少妇亚洲综合色aaa.| 91精品国产国语对白视频| 一级黄色大片毛片| 涩涩av久久男人的天堂| 欧美黑人欧美精品刺激| 男女下面插进去视频免费观看| 欧美中文综合在线视频| 啦啦啦啦在线视频资源| 久久国产精品大桥未久av| 久久久久国产一级毛片高清牌| 一边摸一边抽搐一进一出视频| av视频免费观看在线观看| 最新的欧美精品一区二区| 91老司机精品| 精品亚洲成国产av| 亚洲五月婷婷丁香| 国产亚洲精品一区二区www | 国产又爽黄色视频| 两性夫妻黄色片| 啦啦啦在线免费观看视频4| 极品少妇高潮喷水抽搐| 十八禁人妻一区二区| xxxhd国产人妻xxx| 新久久久久国产一级毛片| 男女之事视频高清在线观看| 9色porny在线观看| videos熟女内射| 高清黄色对白视频在线免费看| 欧美大码av| 天天躁夜夜躁狠狠躁躁| 久久免费观看电影| 9191精品国产免费久久| 啦啦啦免费观看视频1| 男女边摸边吃奶| 淫妇啪啪啪对白视频 | 成年人黄色毛片网站| 久久国产亚洲av麻豆专区| 丁香六月天网| 成人手机av| 精品高清国产在线一区| 亚洲精品久久成人aⅴ小说| a在线观看视频网站| 91精品国产国语对白视频| 亚洲第一av免费看| 视频在线观看一区二区三区| 久久青草综合色| 国产日韩欧美亚洲二区| 免费高清在线观看视频在线观看| 欧美日韩国产mv在线观看视频| 精品视频人人做人人爽| 热99re8久久精品国产| 中文字幕精品免费在线观看视频| 午夜福利,免费看| 老司机午夜十八禁免费视频| 国产亚洲av高清不卡| 成年人免费黄色播放视频| 99久久精品国产亚洲精品| a级毛片黄视频| 男女床上黄色一级片免费看| 狂野欧美激情性xxxx| 少妇 在线观看| 性高湖久久久久久久久免费观看| 婷婷丁香在线五月| 免费在线观看影片大全网站| 在线看a的网站| 9色porny在线观看| 巨乳人妻的诱惑在线观看| 免费看十八禁软件| 咕卡用的链子| 国产精品二区激情视频| 亚洲成国产人片在线观看| 午夜两性在线视频| 满18在线观看网站| 少妇猛男粗大的猛烈进出视频| 午夜影院在线不卡| 国产极品粉嫩免费观看在线| 黄色视频在线播放观看不卡| 国产精品一区二区在线观看99| 多毛熟女@视频| 欧美 日韩 精品 国产| 51午夜福利影视在线观看| 制服人妻中文乱码| 亚洲欧美色中文字幕在线| 国产色视频综合| 咕卡用的链子| 99re6热这里在线精品视频| 成人国语在线视频| 啦啦啦啦在线视频资源| 国产免费福利视频在线观看| 久久综合国产亚洲精品| 一个人免费看片子| 最黄视频免费看| 国产精品免费视频内射| 岛国毛片在线播放| av网站在线播放免费| 午夜两性在线视频| 精品国产乱码久久久久久男人| 两个人看的免费小视频| 久久狼人影院| 性少妇av在线| 欧美日韩黄片免| 精品免费久久久久久久清纯 | 国产激情久久老熟女| 欧美日韩亚洲高清精品| 最近最新免费中文字幕在线| 男男h啪啪无遮挡| 亚洲精品久久成人aⅴ小说| 如日韩欧美国产精品一区二区三区| 国产精品久久久久成人av| 久久av网站| 日韩中文字幕视频在线看片| 亚洲全国av大片| 一级,二级,三级黄色视频| 精品少妇内射三级| 精品人妻熟女毛片av久久网站| 成人三级做爰电影| videos熟女内射| 超色免费av| 午夜福利视频精品| 亚洲国产精品999| 精品欧美一区二区三区在线| 黑人巨大精品欧美一区二区蜜桃| 一区在线观看完整版| 免费一级毛片在线播放高清视频 | 人人妻人人澡人人爽人人夜夜| 国产精品 欧美亚洲| www.自偷自拍.com| 香蕉国产在线看| 岛国在线观看网站| 国产亚洲欧美在线一区二区| 高潮久久久久久久久久久不卡| 午夜福利视频在线观看免费| 国产色视频综合| 在线 av 中文字幕| 肉色欧美久久久久久久蜜桃| 午夜福利视频精品| 蜜桃在线观看..| 亚洲第一欧美日韩一区二区三区 | 亚洲第一av免费看| 两人在一起打扑克的视频| 久久人人爽人人片av| 女性生殖器流出的白浆| 亚洲成国产人片在线观看| 一本久久精品| 中文字幕制服av| 日韩熟女老妇一区二区性免费视频| 国产精品 国内视频| svipshipincom国产片| 高清欧美精品videossex| 亚洲伊人色综图| 国产av精品麻豆| 美女脱内裤让男人舔精品视频| 女人爽到高潮嗷嗷叫在线视频| 十八禁高潮呻吟视频| 午夜福利乱码中文字幕| 欧美变态另类bdsm刘玥| 午夜福利在线观看吧| 高清在线国产一区| 飞空精品影院首页| 亚洲人成77777在线视频| 最新在线观看一区二区三区| 午夜激情av网站| 日韩欧美一区二区三区在线观看 | 91字幕亚洲| 人妻 亚洲 视频| 亚洲av电影在线进入| 午夜久久久在线观看| 亚洲中文字幕日韩| 国产精品.久久久| 中文精品一卡2卡3卡4更新| 亚洲熟女毛片儿| www.999成人在线观看| 香蕉国产在线看| 看免费av毛片| 亚洲专区中文字幕在线| 中文字幕另类日韩欧美亚洲嫩草| 国产精品免费大片| 一边摸一边做爽爽视频免费| 19禁男女啪啪无遮挡网站| e午夜精品久久久久久久| 9热在线视频观看99| 美女高潮喷水抽搐中文字幕| 女警被强在线播放| 国产在视频线精品| 一区福利在线观看| 国产一区二区三区av在线| 美女福利国产在线| 国产日韩欧美在线精品| 18在线观看网站| 中文字幕高清在线视频| h视频一区二区三区| 91国产中文字幕| 亚洲自偷自拍图片 自拍| av有码第一页| 日日夜夜操网爽| 欧美黄色片欧美黄色片| 中文字幕最新亚洲高清| 91麻豆精品激情在线观看国产 | 国产精品二区激情视频| 欧美久久黑人一区二区| 人妻 亚洲 视频| 亚洲综合色网址| 免费黄频网站在线观看国产| 国产麻豆69| 黄色 视频免费看| 免费人妻精品一区二区三区视频| 午夜两性在线视频| 一区二区三区精品91| 下体分泌物呈黄色| 亚洲国产欧美一区二区综合| 亚洲精品一卡2卡三卡4卡5卡 | netflix在线观看网站| 精品人妻熟女毛片av久久网站| 久久国产精品人妻蜜桃| 国产91精品成人一区二区三区 | 国产亚洲av高清不卡| 亚洲黑人精品在线| 一边摸一边抽搐一进一出视频| 久久人妻熟女aⅴ| 美国免费a级毛片| 免费高清在线观看日韩| 精品少妇内射三级| 黄色视频,在线免费观看| 中文字幕av电影在线播放| 人人妻人人爽人人添夜夜欢视频| 亚洲国产精品一区三区| 亚洲久久久国产精品| 国产日韩一区二区三区精品不卡| 宅男免费午夜| 淫妇啪啪啪对白视频 | 99国产精品一区二区蜜桃av | 免费在线观看日本一区| 精品少妇黑人巨大在线播放| 国产av国产精品国产| 高清视频免费观看一区二区| 丰满迷人的少妇在线观看| 久久精品国产亚洲av香蕉五月 | 久久久久视频综合| 国产xxxxx性猛交| 国产欧美日韩一区二区精品| 国产精品熟女久久久久浪| 伦理电影免费视频| 欧美精品亚洲一区二区| 18禁国产床啪视频网站| 麻豆国产av国片精品| av电影中文网址| 久久狼人影院| 99热国产这里只有精品6| 岛国在线观看网站| 日韩,欧美,国产一区二区三区| 一本—道久久a久久精品蜜桃钙片| 又紧又爽又黄一区二区| 亚洲性夜色夜夜综合| 国产高清videossex| 韩国精品一区二区三区| 母亲3免费完整高清在线观看| 免费一级毛片在线播放高清视频 | 欧美国产精品va在线观看不卡| 亚洲国产精品999| 午夜激情久久久久久久| 亚洲国产欧美日韩在线播放| 99国产精品一区二区蜜桃av | 亚洲精品国产一区二区精华液| 精品国产一区二区三区四区第35| 免费久久久久久久精品成人欧美视频| 少妇被粗大的猛进出69影院| 久久热在线av| 91精品三级在线观看| 亚洲精品久久午夜乱码| 岛国毛片在线播放| 午夜91福利影院| 免费高清在线观看日韩| 欧美 亚洲 国产 日韩一| 少妇裸体淫交视频免费看高清 | 久久人妻熟女aⅴ| 国产精品熟女久久久久浪| 亚洲一码二码三码区别大吗| 又紧又爽又黄一区二区| 国产av又大| 999精品在线视频| 亚洲人成电影免费在线| 女人高潮潮喷娇喘18禁视频| 各种免费的搞黄视频| 日本欧美视频一区| av一本久久久久| 国产精品国产av在线观看| 亚洲欧洲精品一区二区精品久久久| 国产精品久久久人人做人人爽| 亚洲伊人久久精品综合| 欧美国产精品va在线观看不卡| 人妻人人澡人人爽人人| 两个人免费观看高清视频| 成人国语在线视频| 国产精品1区2区在线观看. | 免费不卡黄色视频| 欧美日韩亚洲高清精品| 亚洲av国产av综合av卡| 中文字幕高清在线视频| 天天影视国产精品| 在线观看免费午夜福利视频| 日韩一区二区三区影片| 亚洲色图 男人天堂 中文字幕| 久久久久久人人人人人| 国产成人免费观看mmmm| 免费少妇av软件| 午夜激情久久久久久久| 欧美老熟妇乱子伦牲交| 久久99热这里只频精品6学生| 国产三级黄色录像| 91九色精品人成在线观看| 人人妻人人澡人人爽人人夜夜| 一本综合久久免费| a级毛片在线看网站| 黄网站色视频无遮挡免费观看| 在线观看免费日韩欧美大片| 天天躁夜夜躁狠狠躁躁| 成年女人毛片免费观看观看9 | 丝袜脚勾引网站| 色精品久久人妻99蜜桃| 国产一级毛片在线| av片东京热男人的天堂| 精品国产一区二区三区四区第35| 久久影院123| 99国产精品免费福利视频| 亚洲色图综合在线观看| 国产av精品麻豆| 成年人午夜在线观看视频| 这个男人来自地球电影免费观看| 午夜视频精品福利| 欧美黄色淫秽网站| 日韩欧美一区二区三区在线观看 | 日本一区二区免费在线视频| 精品亚洲成国产av| 欧美少妇被猛烈插入视频| 黄片小视频在线播放| xxxhd国产人妻xxx| 国产有黄有色有爽视频| 黄色视频在线播放观看不卡| 欧美日韩视频精品一区| 乱人伦中国视频| 18禁裸乳无遮挡动漫免费视频| 国产主播在线观看一区二区| 人人澡人人妻人| 人人妻人人澡人人爽人人夜夜| 不卡av一区二区三区| 母亲3免费完整高清在线观看| 动漫黄色视频在线观看| 叶爱在线成人免费视频播放| 精品国产乱子伦一区二区三区 | 黄色 视频免费看| 男女免费视频国产| 窝窝影院91人妻| 亚洲精品久久成人aⅴ小说| av片东京热男人的天堂| 女人被躁到高潮嗷嗷叫费观| 久久精品成人免费网站| 亚洲va日本ⅴa欧美va伊人久久 | 国产欧美日韩综合在线一区二区| 妹子高潮喷水视频| 亚洲激情五月婷婷啪啪| 欧美精品啪啪一区二区三区 | 亚洲伊人久久精品综合| 欧美 日韩 精品 国产| 麻豆av在线久日| 五月天丁香电影| 国产一区二区三区在线臀色熟女 | 亚洲精品国产区一区二| 欧美精品人与动牲交sv欧美| 老司机深夜福利视频在线观看 | 精品卡一卡二卡四卡免费| 国产精品九九99| 看免费av毛片| 天堂中文最新版在线下载| 99热全是精品| 久久这里只有精品19| 国产淫语在线视频| 亚洲精品久久午夜乱码| 免费女性裸体啪啪无遮挡网站| 国产亚洲av高清不卡| 女警被强在线播放| 亚洲av国产av综合av卡| 自线自在国产av| 2018国产大陆天天弄谢| 亚洲国产日韩一区二区| 色婷婷久久久亚洲欧美| 亚洲五月婷婷丁香| av在线app专区| 男女国产视频网站| www.av在线官网国产| 两性午夜刺激爽爽歪歪视频在线观看 | 高清欧美精品videossex| 丝袜人妻中文字幕| 人人妻,人人澡人人爽秒播| 亚洲免费av在线视频| 乱人伦中国视频| 青草久久国产| 午夜老司机福利片| 12—13女人毛片做爰片一| 久久九九热精品免费| 丰满饥渴人妻一区二区三| 免费高清在线观看视频在线观看| 天堂俺去俺来也www色官网| 青草久久国产| 久久久久久久大尺度免费视频| 精品国产乱码久久久久久小说| 亚洲全国av大片| 99久久精品国产亚洲精品| 国产精品麻豆人妻色哟哟久久| 大码成人一级视频| 国产xxxxx性猛交| 亚洲第一欧美日韩一区二区三区 | 我的亚洲天堂| 亚洲精品美女久久av网站| 中文字幕高清在线视频| 99国产综合亚洲精品| 夜夜骑夜夜射夜夜干| 国产精品一区二区在线观看99| 男男h啪啪无遮挡| 肉色欧美久久久久久久蜜桃| 超碰97精品在线观看| 狂野欧美激情性xxxx| 亚洲成av片中文字幕在线观看| 无限看片的www在线观看| 91国产中文字幕| tube8黄色片| 国产成人av激情在线播放| 天堂中文最新版在线下载| 国产精品久久久久久精品古装| 国产97色在线日韩免费| 99精品欧美一区二区三区四区| 欧美人与性动交α欧美软件| 亚洲熟女精品中文字幕| 亚洲第一av免费看| 国精品久久久久久国模美| 啪啪无遮挡十八禁网站| 老熟妇仑乱视频hdxx| www.精华液| 90打野战视频偷拍视频| 一级a爱视频在线免费观看| 亚洲国产av新网站| 国产免费av片在线观看野外av| 亚洲av成人一区二区三| 视频在线观看一区二区三区| 男女午夜视频在线观看| 我的亚洲天堂| 永久免费av网站大全| 男人舔女人的私密视频| 国产片内射在线| 国产一区二区三区在线臀色熟女 | 日本vs欧美在线观看视频| a级片在线免费高清观看视频| 久久精品国产综合久久久| 又大又爽又粗| 男女午夜视频在线观看| 我的亚洲天堂| 精品福利永久在线观看| 亚洲精品一区蜜桃| 久久久久久久久免费视频了| 高清黄色对白视频在线免费看| 18禁观看日本| 搡老乐熟女国产| 久久中文字幕一级| 黄色 视频免费看| 亚洲欧洲日产国产| 大陆偷拍与自拍| 天堂俺去俺来也www色官网| 亚洲国产欧美在线一区| 国产欧美亚洲国产| 亚洲五月色婷婷综合| 亚洲成国产人片在线观看| 成年人免费黄色播放视频| 日韩 亚洲 欧美在线| e午夜精品久久久久久久| 亚洲久久久国产精品| 9热在线视频观看99| 日韩一卡2卡3卡4卡2021年| 免费少妇av软件| av天堂久久9| 久久精品国产a三级三级三级| 国产精品一区二区在线观看99| 欧美人与性动交α欧美软件| 桃红色精品国产亚洲av| 嫩草影视91久久| 国产亚洲午夜精品一区二区久久| 高清视频免费观看一区二区| 高潮久久久久久久久久久不卡| 欧美日本中文国产一区发布| 久久国产亚洲av麻豆专区| 久久女婷五月综合色啪小说| 欧美日韩国产mv在线观看视频| 在线看a的网站| 秋霞在线观看毛片| 18禁黄网站禁片午夜丰满| a级片在线免费高清观看视频| 老司机影院成人| 国产成人一区二区三区免费视频网站| 久久天堂一区二区三区四区| 久久精品亚洲熟妇少妇任你| 国产成人一区二区三区免费视频网站| 欧美黄色片欧美黄色片| 曰老女人黄片| 少妇人妻久久综合中文| 男女床上黄色一级片免费看|