• <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
    国产成人精品久久久久久| 久久精品国产a三级三级三级| 别揉我奶头 嗯啊视频| 久久久久国产精品人妻一区二区| 色吧在线观看| 看非洲黑人一级黄片| 亚洲精品久久久久久婷婷小说| 亚洲一区二区三区欧美精品 | 美女主播在线视频| 美女主播在线视频| 久久精品国产亚洲av天美| 天堂俺去俺来也www色官网| 天美传媒精品一区二区| 美女cb高潮喷水在线观看| 97精品久久久久久久久久精品| 国产精品一区www在线观看| 美女高潮的动态| 一级av片app| 久久久久久久国产电影| 国产91av在线免费观看| 日日啪夜夜爽| 成人午夜精彩视频在线观看| kizo精华| 久久久久国产网址| 一个人看的www免费观看视频| 欧美xxⅹ黑人| 日韩av在线免费看完整版不卡| 国产午夜精品久久久久久一区二区三区| 亚洲av成人精品一二三区| 精品久久久久久久久av| 国产精品精品国产色婷婷| kizo精华| 日韩人妻高清精品专区| 夫妻午夜视频| 性色av一级| 熟女av电影| 午夜日本视频在线| 视频中文字幕在线观看| 色视频www国产| 日韩在线高清观看一区二区三区| 日产精品乱码卡一卡2卡三| 国产精品爽爽va在线观看网站| 亚洲国产色片| 国产精品人妻久久久久久| 麻豆久久精品国产亚洲av| 国产欧美另类精品又又久久亚洲欧美| 国产片特级美女逼逼视频| 亚洲三级黄色毛片| 国内精品美女久久久久久| 国产一区亚洲一区在线观看| 日韩亚洲欧美综合| 青春草视频在线免费观看| 高清午夜精品一区二区三区| 国产精品人妻久久久影院| 国产亚洲精品久久久com| 久久久a久久爽久久v久久| 国产熟女欧美一区二区| 国产成人精品婷婷| 久久精品夜色国产| 另类亚洲欧美激情| 哪个播放器可以免费观看大片| 校园人妻丝袜中文字幕| 干丝袜人妻中文字幕| 一级爰片在线观看| 黄色配什么色好看| 国产免费一级a男人的天堂| 久久99精品国语久久久| 亚洲怡红院男人天堂| 97在线视频观看| 久久韩国三级中文字幕| 欧美xxxx性猛交bbbb| 欧美成人午夜免费资源| 中文字幕久久专区| 亚洲最大成人中文| 在线a可以看的网站| 国产免费福利视频在线观看| 26uuu在线亚洲综合色| 嫩草影院新地址| 午夜精品国产一区二区电影 | 九九久久精品国产亚洲av麻豆| 免费观看性生交大片5| 春色校园在线视频观看| 97热精品久久久久久| 七月丁香在线播放| 亚洲伊人久久精品综合| 免费观看a级毛片全部| 国产成年人精品一区二区| 国产av国产精品国产| 中文字幕制服av| 成人免费观看视频高清| 一级毛片 在线播放| 久久国内精品自在自线图片| 人人妻人人看人人澡| av在线观看视频网站免费| 色视频在线一区二区三区| 亚洲精品中文字幕在线视频 | 晚上一个人看的免费电影| 国内精品宾馆在线| 又黄又爽又刺激的免费视频.| 天堂中文最新版在线下载 | 亚洲人成网站在线播| 国产精品一区二区在线观看99| 一级毛片 在线播放| 亚洲四区av| 日韩一区二区视频免费看| 国产黄色免费在线视频| 欧美激情在线99| 亚洲av在线观看美女高潮| 熟女电影av网| 亚洲精品国产成人久久av| 免费看av在线观看网站| 插逼视频在线观看| 综合色av麻豆| 久久精品国产亚洲av天美| 亚洲电影在线观看av| 又粗又硬又长又爽又黄的视频| 熟女人妻精品中文字幕| 看免费成人av毛片| 久久99热这里只有精品18| 禁无遮挡网站| 尾随美女入室| 精品国产乱码久久久久久小说| 日日摸夜夜添夜夜爱| 又爽又黄a免费视频| 国产免费福利视频在线观看| 天美传媒精品一区二区| 国产日韩欧美亚洲二区| 神马国产精品三级电影在线观看| 日日摸夜夜添夜夜爱| 国产综合懂色| 我的老师免费观看完整版| 91久久精品国产一区二区三区| 成年免费大片在线观看| 色网站视频免费| 爱豆传媒免费全集在线观看| 超碰97精品在线观看| 白带黄色成豆腐渣| 国产成人福利小说| 国产v大片淫在线免费观看| 99久久九九国产精品国产免费| 亚洲国产精品专区欧美| 五月天丁香电影| 在线播放无遮挡| 婷婷色综合大香蕉| 亚洲四区av| 久热这里只有精品99| 久久久午夜欧美精品| 亚洲精品久久午夜乱码| 尤物成人国产欧美一区二区三区| 日韩av不卡免费在线播放| 在线观看美女被高潮喷水网站| 亚洲不卡免费看| 日本免费在线观看一区| 男女国产视频网站| 久久久久久久大尺度免费视频| 黄色日韩在线| 五月玫瑰六月丁香| 干丝袜人妻中文字幕| 国产老妇女一区| 欧美国产精品一级二级三级 | 久久久久国产精品人妻一区二区| 亚洲av二区三区四区| 中文精品一卡2卡3卡4更新| 麻豆乱淫一区二区| 色综合色国产| 亚洲最大成人av| 亚洲人成网站在线播| 日本黄大片高清| 色哟哟·www| 亚洲精品中文字幕在线视频 | 观看免费一级毛片| 亚洲精品日韩在线中文字幕| 国语对白做爰xxxⅹ性视频网站| 99久久精品一区二区三区| 亚洲熟女精品中文字幕| 亚洲欧美中文字幕日韩二区| 亚洲最大成人中文| 一级毛片电影观看| 在线观看人妻少妇| 亚洲精品亚洲一区二区| 国语对白做爰xxxⅹ性视频网站| 国产午夜精品一二区理论片| 99热网站在线观看| 亚洲精品456在线播放app| 99精国产麻豆久久婷婷| 精品少妇久久久久久888优播| 最近最新中文字幕大全电影3| 午夜福利在线在线| 国产亚洲精品久久久com| 中国美白少妇内射xxxbb| 免费av毛片视频| 狂野欧美激情性bbbbbb| 人人妻人人澡人人爽人人夜夜| 少妇的逼好多水| 欧美3d第一页| 男的添女的下面高潮视频| 欧美激情在线99| 男女无遮挡免费网站观看| 日本色播在线视频| 51国产日韩欧美| 成人免费观看视频高清| 精品一区在线观看国产| 中文在线观看免费www的网站| 欧美亚洲 丝袜 人妻 在线| 美女内射精品一级片tv| 亚洲在久久综合| 丰满人妻一区二区三区视频av| 99热6这里只有精品| 三级男女做爰猛烈吃奶摸视频| 成人特级av手机在线观看| 九九爱精品视频在线观看| 亚洲av欧美aⅴ国产| 欧美激情在线99| av专区在线播放| 成人国产麻豆网| 狂野欧美白嫩少妇大欣赏| 国产白丝娇喘喷水9色精品| 亚洲内射少妇av| 精品99又大又爽又粗少妇毛片| 欧美精品一区二区大全| av国产久精品久网站免费入址| 国产午夜精品一二区理论片| 欧美精品国产亚洲| 直男gayav资源| 国产乱人视频| 超碰av人人做人人爽久久| 毛片一级片免费看久久久久| 国产成人91sexporn| 一边亲一边摸免费视频| 亚洲婷婷狠狠爱综合网| 亚洲人成网站在线播| 麻豆成人av视频| 亚洲精品,欧美精品| 激情 狠狠 欧美| 好男人在线观看高清免费视频| 亚洲三级黄色毛片| 身体一侧抽搐| 欧美精品国产亚洲| 国产淫片久久久久久久久| 欧美一级a爱片免费观看看| 中文天堂在线官网| 国国产精品蜜臀av免费| 久久久久精品性色| kizo精华| 日本午夜av视频| 性插视频无遮挡在线免费观看| 丝瓜视频免费看黄片| 99久久九九国产精品国产免费| 22中文网久久字幕| 亚洲国产精品专区欧美| 插阴视频在线观看视频| videos熟女内射| 九九在线视频观看精品| 亚洲精品乱久久久久久| 春色校园在线视频观看| 国产高清三级在线| 精品熟女少妇av免费看| 有码 亚洲区| 51国产日韩欧美| av在线亚洲专区| 亚洲欧美精品专区久久| 啦啦啦中文免费视频观看日本| 国产爽快片一区二区三区| 一级毛片黄色毛片免费观看视频| 大片电影免费在线观看免费| 97精品久久久久久久久久精品| 精品国产一区二区三区久久久樱花 | 男女无遮挡免费网站观看| 一区二区av电影网| 亚洲成人一二三区av| 欧美精品人与动牲交sv欧美| 热99国产精品久久久久久7| 国产毛片a区久久久久| 嫩草影院精品99| 天堂网av新在线| 亚洲欧洲国产日韩| 欧美极品一区二区三区四区| 国产午夜精品久久久久久一区二区三区| av女优亚洲男人天堂| 免费av不卡在线播放| 久久综合国产亚洲精品| 午夜免费男女啪啪视频观看| 大码成人一级视频| 丝袜喷水一区| 日韩制服骚丝袜av| 丝瓜视频免费看黄片| 亚洲成人精品中文字幕电影| 建设人人有责人人尽责人人享有的 | 内射极品少妇av片p| 国产亚洲午夜精品一区二区久久 | 欧美少妇被猛烈插入视频| av在线老鸭窝| 成人国产av品久久久| 国产成人a∨麻豆精品| 三级国产精品片| 欧美最新免费一区二区三区| a级毛片免费高清观看在线播放| 午夜精品国产一区二区电影 | 日韩欧美精品免费久久| 国产69精品久久久久777片| 一级毛片黄色毛片免费观看视频| 91久久精品国产一区二区成人| 亚洲欧美成人精品一区二区| 91久久精品国产一区二区三区| 精品久久久噜噜| 国产精品久久久久久精品电影| 国产欧美日韩一区二区三区在线 | 亚洲不卡免费看| 欧美日韩视频精品一区| 欧美丝袜亚洲另类| 一区二区三区精品91| 水蜜桃什么品种好| 中文字幕av成人在线电影| 久久久色成人| 午夜福利在线观看免费完整高清在| 王馨瑶露胸无遮挡在线观看| 老司机影院毛片| 国产成人免费无遮挡视频| 国产伦在线观看视频一区| 亚洲内射少妇av| 国产成人a区在线观看| 亚洲av男天堂| 久久久a久久爽久久v久久| 日韩制服骚丝袜av| 新久久久久国产一级毛片| 精品国产三级普通话版| av线在线观看网站| 26uuu在线亚洲综合色| 免费看日本二区| 日产精品乱码卡一卡2卡三| 水蜜桃什么品种好| 亚洲精品自拍成人| 国产乱来视频区| 亚洲一区二区三区欧美精品 | 五月伊人婷婷丁香| 国产欧美日韩一区二区三区在线 | 真实男女啪啪啪动态图| 亚洲精品乱久久久久久| 久久久久久久亚洲中文字幕| 亚洲国产欧美在线一区| 精品一区二区三卡| 亚洲欧洲日产国产| 18+在线观看网站| 夫妻午夜视频| 免费看不卡的av| a级毛片免费高清观看在线播放| 成人一区二区视频在线观看| 亚洲国产欧美人成| 国产精品人妻久久久久久| 欧美另类一区| 亚洲精品一二三| 男女国产视频网站| av在线亚洲专区| 99九九线精品视频在线观看视频| 特级一级黄色大片| 69人妻影院| 一个人看视频在线观看www免费| 免费大片18禁| 人妻制服诱惑在线中文字幕| 一个人观看的视频www高清免费观看| 国产成人精品福利久久| 六月丁香七月| 你懂的网址亚洲精品在线观看| 日韩欧美精品免费久久| 亚洲第一区二区三区不卡| 高清欧美精品videossex| 日本免费在线观看一区| 国产片特级美女逼逼视频| 我的老师免费观看完整版| 日韩亚洲欧美综合| 国产乱人视频| 少妇人妻精品综合一区二区| 免费高清在线观看视频在线观看| 美女高潮的动态| 国产精品无大码| 国产一级毛片在线| 日韩不卡一区二区三区视频在线| 日日啪夜夜爽| 在线观看三级黄色| 成人午夜精彩视频在线观看| av在线天堂中文字幕| 人人妻人人爽人人添夜夜欢视频 | 成人黄色视频免费在线看| 亚洲精品中文字幕在线视频 | 另类亚洲欧美激情| 色吧在线观看| 欧美成人精品欧美一级黄| 亚洲一区二区三区欧美精品 | 在线精品无人区一区二区三 | 国产精品女同一区二区软件| 国产大屁股一区二区在线视频| 久久久亚洲精品成人影院| 六月丁香七月| 人妻一区二区av| 国产精品人妻久久久影院| 国产欧美日韩一区二区三区在线 | 亚洲欧洲日产国产| 日日啪夜夜爽| 日韩欧美精品v在线| 亚洲精品日本国产第一区| .国产精品久久| 精品久久久精品久久久| 成人毛片60女人毛片免费| 简卡轻食公司| 2021天堂中文幕一二区在线观| 又爽又黄无遮挡网站| 五月伊人婷婷丁香| 新久久久久国产一级毛片| 亚洲av免费在线观看| 男人添女人高潮全过程视频| av一本久久久久| av在线亚洲专区| 嫩草影院新地址| 女人久久www免费人成看片| 王馨瑶露胸无遮挡在线观看| 大话2 男鬼变身卡| 久久人人爽人人片av| 亚洲成人精品中文字幕电影| 高清毛片免费看| 亚洲伊人久久精品综合| 国产色婷婷99| 日韩强制内射视频| 天天一区二区日本电影三级| 日韩欧美精品免费久久| 中文天堂在线官网| 日韩一区二区三区影片| 欧美bdsm另类| 老女人水多毛片| 亚洲国产欧美在线一区| 亚洲精华国产精华液的使用体验| 免费观看性生交大片5| 国产白丝娇喘喷水9色精品| 性色avwww在线观看| 国产精品一区二区性色av| 久久久久国产精品人妻一区二区| 亚洲天堂av无毛| 99久久人妻综合| 丰满人妻一区二区三区视频av| 亚洲三级黄色毛片| 99九九线精品视频在线观看视频| 成年版毛片免费区| 日韩免费高清中文字幕av| 欧美精品一区二区大全| 成人鲁丝片一二三区免费| 久久精品久久久久久噜噜老黄| 新久久久久国产一级毛片| 一级片'在线观看视频| 人人妻人人看人人澡| 亚洲精品久久午夜乱码| 插逼视频在线观看| 国产 精品1| 国产黄色免费在线视频| 在线亚洲精品国产二区图片欧美 | 免费少妇av软件| 丰满乱子伦码专区| 国产一区二区三区综合在线观看 | 国产成人aa在线观看| 黄片无遮挡物在线观看| 尤物成人国产欧美一区二区三区| 日韩av免费高清视频| 亚洲精品一二三| 国产高清有码在线观看视频| 欧美日韩一区二区视频在线观看视频在线 | 在线观看人妻少妇| 男女边吃奶边做爰视频| 免费看光身美女| 97热精品久久久久久| 日韩欧美 国产精品| 麻豆精品久久久久久蜜桃| 男女下面进入的视频免费午夜| tube8黄色片| 色哟哟·www| 国产精品99久久久久久久久| videossex国产| 亚洲av欧美aⅴ国产| 超碰97精品在线观看| 精品一区二区免费观看| 久久精品熟女亚洲av麻豆精品| 日韩av在线免费看完整版不卡| 亚洲欧美日韩另类电影网站 | 国产成年人精品一区二区| 18+在线观看网站| 久久6这里有精品| 成年女人看的毛片在线观看| 久久久久性生活片| 又大又黄又爽视频免费| 黄色欧美视频在线观看| 午夜免费鲁丝| 国产精品久久久久久久电影| 久久精品国产亚洲网站| 99久久精品热视频| 大又大粗又爽又黄少妇毛片口| 亚洲av电影在线观看一区二区三区 | 亚洲精品国产成人久久av| 九草在线视频观看| 国产精品嫩草影院av在线观看| 自拍偷自拍亚洲精品老妇| 久久精品国产亚洲av天美| 七月丁香在线播放| 一级爰片在线观看| 丰满人妻一区二区三区视频av| 97在线视频观看| 国产乱人偷精品视频| 亚洲国产欧美在线一区| 中文字幕久久专区| 国产av不卡久久| 一级毛片aaaaaa免费看小| 男女边吃奶边做爰视频| 亚洲人成网站在线观看播放| 99九九线精品视频在线观看视频| 午夜免费男女啪啪视频观看| 亚洲国产欧美人成| 在线观看一区二区三区激情| 各种免费的搞黄视频| 亚洲在线观看片| 大码成人一级视频| 亚洲三级黄色毛片| 久久久精品94久久精品| 日韩在线高清观看一区二区三区| 听说在线观看完整版免费高清| 日韩在线高清观看一区二区三区| 2018国产大陆天天弄谢| 精品亚洲乱码少妇综合久久| 亚洲精品日韩av片在线观看| 亚洲精品456在线播放app| 精华霜和精华液先用哪个| 如何舔出高潮| 国产一区二区在线观看日韩| 欧美极品一区二区三区四区| 国内揄拍国产精品人妻在线| 丝瓜视频免费看黄片| 亚洲av.av天堂| 高清在线视频一区二区三区| 亚洲精品日韩av片在线观看| 性色av一级| 亚洲熟女精品中文字幕| 一级黄片播放器| 日产精品乱码卡一卡2卡三| 国产精品嫩草影院av在线观看| 久久女婷五月综合色啪小说 | 禁无遮挡网站| .国产精品久久| 卡戴珊不雅视频在线播放| 99视频精品全部免费 在线| 精品久久久精品久久久| 男的添女的下面高潮视频| 内地一区二区视频在线| 欧美少妇被猛烈插入视频| 波野结衣二区三区在线| 三级国产精品片| 国产精品国产三级国产av玫瑰| 伦精品一区二区三区| 一本一本综合久久| 精品国产乱码久久久久久小说| 欧美潮喷喷水| 美女内射精品一级片tv| 只有这里有精品99| 久久久久久久亚洲中文字幕| 91在线精品国自产拍蜜月| 在线看a的网站| 女的被弄到高潮叫床怎么办| 国产一区二区亚洲精品在线观看| 女人久久www免费人成看片| 久久久久国产精品人妻一区二区| 久久久久久伊人网av| 国产伦在线观看视频一区| 国产精品熟女久久久久浪| eeuss影院久久| 波多野结衣巨乳人妻| 亚洲国产最新在线播放| 国产精品成人在线| 大又大粗又爽又黄少妇毛片口| 亚洲国产精品成人综合色| 22中文网久久字幕| av黄色大香蕉| 亚洲av一区综合| 1000部很黄的大片| 内地一区二区视频在线| 我要看日韩黄色一级片| 国产午夜福利久久久久久| 晚上一个人看的免费电影| 大码成人一级视频| 白带黄色成豆腐渣| 国产精品国产三级国产专区5o| 亚州av有码| 大香蕉久久网| 亚洲成人久久爱视频| 少妇人妻久久综合中文| 特级一级黄色大片| 黄色视频在线播放观看不卡| 久久久成人免费电影| 国产伦理片在线播放av一区| 男人和女人高潮做爰伦理| 成人国产麻豆网| 亚洲国产精品999| 精品一区二区三区视频在线| 国产伦精品一区二区三区视频9| 欧美高清成人免费视频www| 美女cb高潮喷水在线观看| 欧美bdsm另类| 亚洲欧美成人精品一区二区| 在现免费观看毛片| 久久久久精品久久久久真实原创| 亚洲欧美成人精品一区二区| 五月天丁香电影| 精品一区二区免费观看| 亚洲在久久综合| 久久精品夜色国产| 久久久久久久久大av| 成人高潮视频无遮挡免费网站| 亚洲精品第二区| 18禁裸乳无遮挡动漫免费视频 | 国产在视频线精品| 少妇被粗大猛烈的视频| 波多野结衣巨乳人妻|