劉泛函，王華，徐建新

（1省部共建復(fù)雜有色金屬資源清潔利用國家重點實驗室，云南 昆明650093；2昆明理工大學冶金與能源工程學院，云南 昆明650093；3昆明理工大學質(zhì)量發(fā)展研究院，云南 昆明 650093）

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

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

（1省部共建復(fù)雜有色金屬資源清潔利用國家重點實驗室，云南 昆明650093；2昆明理工大學冶金與能源工程學院，云南 昆明650093；3昆明理工大學質(zhì)量發(fā)展研究院，云南 昆明 650093）

基于圖像處理技術(shù)和混沌理論，提出了一種氣體射流圖像的量化方法，用以反映氣液混合性能。本文使用大津法來獲得自適應(yīng)閥值，并得到最佳的圖像分割效果。混合性能（M）用以表征氣體射流的空間分布，其定義為氣體射流像素占比（W）與射流區(qū)域內(nèi)的灰度標準差（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)定性的強弱。

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

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

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

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

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

1　實驗裝置

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

圖1　實驗裝置

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

2　圖像分析

2.1圖像分割

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

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

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

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

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

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

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

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

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

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

步驟4：混合性能指數(shù)可以通過計算在分割區(qū)域中像素灰度值的標準差得到，其表達式如式（5）。

基于整個混合過程的采樣圖像就可以獲得混合性能指數(shù)的時間序列M（t）（t=0，1，2，…，k），通過混合性能指數(shù)的平均值來衡量整個混合過程的混合性能，混合性能指數(shù)的平均值表達式如式（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ā)生機制是在雷諾數(shù)較高且密度差較大的條件下，混合剪切層的不穩(wěn)定性使得噴嘴出口處的射流發(fā)生間歇性頸縮，射流通道不暢從而造成氣體聚集，在這個過程中射流內(nèi)部的壓力急劇增大，當壓力達到一定程度后就會導致射流脹鼓，脹鼓發(fā)生的演化過程通過示意圖進行展示。

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

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

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

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

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

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

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

施紅輝等［19］用“氣袋”模型形象的對氣體射流的不穩(wěn)定性進行了解釋，他認為氣體射流是由不可壓縮的流體包裹著一個可壓縮的“氣袋”，由此產(chǎn)生非線性的震蕩壓力，使得射流的氣液邊界不穩(wěn)定。在亞音速工況下，“氣袋”的長度隨修正弗魯?shù)聰?shù)的增加有著較為明顯的增長，“氣袋”越長其內(nèi)部的能量越難以積聚，脹鼓和回擊現(xiàn)象也就越難以發(fā)生，因此隨著修正弗魯?shù)聰?shù)的增加，氣體射流越穩(wě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ù)最大，同時也是實驗中混合反應(yīng)性能和穩(wěn)定性的最好的工況。

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

其中

圖8　混合性能變化趨勢

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

通過對圖10的觀察，可以明顯地發(fā)現(xiàn)混合性能指數(shù)和修正弗魯?shù)聰?shù)具有相似的趨勢，認為混合性能指數(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ù)

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

通過Visual Recurrence Analysis軟件使用虛假鄰點法和平均互信息法計算時間序列的嵌入維數(shù)（d）和延遲時間（t），并對相空間進行重新構(gòu)建。圖11所示為利用最大李雅普諾夫指數(shù)法［20］判斷混合性能指數(shù)的序列是否達到混沌狀態(tài)，如果最大李雅普諾夫指數(shù)等于零，則系統(tǒng)處于臨界狀態(tài)，此時最大李雅普諾夫指數(shù)無法判斷系統(tǒng)是否處于混沌狀態(tài)，就需要利用0-1test方法［21-23］來進一步驗證。不同工況下最大李雅普諾夫指數(shù)的計算結(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)系。通過最小二乘法進行擬合并構(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)定性的強弱。

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

4　結(jié) 論

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

［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］ 戚隆溪，曹勇，王柏懿. 水下欠膨脹高速氣體射流的實驗研究［J］.力學學報，2000，32（6）：667-675.

［7］ 王柏懿，戴振卿，戚隆溪，等. 水下超聲速氣體射流回擊現(xiàn)象的實驗研究［J］. 力學學報，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］ 王曉剛，王超，郭強，等. 二維水槽中高速氣體射流的振蕩流流型研究［J］. 浙江理工大學學報，2009，26（4）：613-618.

［10］ 郭強，施紅輝，王超，等. 準二維水下超聲速垂直過膨脹射流研究［J］. 實驗流體力學，2010，24（3）：6-18.

［11］ 郭強，施紅輝，王超，等. 水下超聲速氣體射流氣液兩相復(fù)雜流動研究［J］. 工程熱物理學報，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］ 施紅輝，郭強，王超，等. 水下超音速氣體射流脹鼓和回擊的關(guān)聯(lián)性研究［J］. 力學學報，2010，42（6）：1206-1210.

［14］ 寧偉征. 偏心射流強化流體混沌混合的研究［D］. 重慶：重慶大學，2012.

［15］ 徐建新. 多相體系攪拌混合效果評價方法及其應(yīng)用研究［D］. 昆明：昆明理工大學，2012.

［16］ 劉政，張嘉藝，羅浩林，等. 混沌對流下的半固態(tài)A356鋁合金初生相形貌演變研究［J］. 金屬學報，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］ 施紅輝，王柏懿，戴振卿. 水下超聲速氣體射流的力學機制研究［J］.中國科學（物理學力學天文學），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。

國家重點基礎(chǔ)研究發(fā)展計劃（2014CB460605）及國家自然科學基金（51174105，51206071，51366005，51406071，51206071）項目。

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