李琪琪，梁彬烽，郭泳盈，曾 皓，蘇掌旭，彭 力

(華南師范大學(xué) a.物理與電信工程學(xué)院;b.信息光電子科技學(xué)院，廣東 廣州 510006)

射流是指從排泄口射出或靠機(jī)械推動(dòng)射入另一流體域內(nèi)的運(yùn)動(dòng)流體. 工業(yè)、農(nóng)業(yè)和生活中的廢棄流體，通常以射流的形式排放到環(huán)境流體中，例如：煙囪排入大氣的廢氣，河流和海洋水域的污水排放，等等. 實(shí)驗(yàn)研究方面，Kamotani和Grebe[1]對(duì)流速比R(R為射流流速與橫流流速之比)較大范圍內(nèi)的紊動(dòng)射流宏觀特性進(jìn)行測(cè)量，主要結(jié)論是射流流速場(chǎng)和溫度場(chǎng)不同導(dǎo)致射流軸線不同，并給出了射流基于速度和溫度的軸線方程；Subramanya和Porey[2]，Pratte和Baines[3]系統(tǒng)地研究了射流軸線方程，側(cè)重于射流的時(shí)均量研究；Andreopoules和Rodi[4-5]研究了射流的時(shí)均結(jié)構(gòu)，重點(diǎn)測(cè)量高階紊動(dòng)量. 在流場(chǎng)顯示方面，黃真理等[6]利用平面激光誘導(dǎo)熒光，對(duì)明槽水流中底孔排放的單孔和多孔射流濃度場(chǎng)進(jìn)行了研究，觀測(cè)到低速比時(shí)射流出現(xiàn)分叉現(xiàn)象和尾渦，并給出了馬蹄渦結(jié)構(gòu)沿水深的變化過(guò)程. 在數(shù)值研究方面,Sykes等[7]運(yùn)用方程湍流模式，探討在高流速比下橫流中射流的渦動(dòng)力學(xué)特性；Demuren[8]運(yùn)用二階動(dòng)量矩湍流模式結(jié)合多重網(wǎng)格技術(shù)計(jì)算低流速比下射流的流場(chǎng)和溫度場(chǎng)；彭文啟[9]采用k-ε湍流模式，結(jié)合多重網(wǎng)格的混合有限分析法對(duì)流速比R=2的流場(chǎng)進(jìn)行了研究，檢驗(yàn)多重網(wǎng)格下混合有限分析法適用復(fù)雜流場(chǎng). 目前國(guó)內(nèi)外對(duì)射流沖擊橫流產(chǎn)生渦旋的實(shí)驗(yàn)方法，渦旋三維特征的測(cè)量研究還很少，并且多數(shù)實(shí)驗(yàn)為單點(diǎn)測(cè)量，較難捕捉到流場(chǎng)中大尺度渦旋結(jié)構(gòu)的瞬時(shí)特征. 本文基于射流沖擊橫流原理，在水槽里產(chǎn)生了穩(wěn)定的渦旋，并采用雙攝像頭和3D-PIV技術(shù)對(duì)渦旋的三維流場(chǎng)特征和結(jié)構(gòu)進(jìn)行研究分析，計(jì)算得到渦旋的速度環(huán)量、渦量場(chǎng)、線速度場(chǎng)以及黏度.

1 實(shí)驗(yàn)原理

1.1 射流沖擊橫流產(chǎn)生渦旋原理

Gaster和Crighton[10]應(yīng)用了穩(wěn)定性理論對(duì)射流進(jìn)行分析，在流動(dòng)顯示中得出射流通過(guò)渦卷吸入環(huán)境水體. Pettersen(1978)[11]研究發(fā)現(xiàn)射流沖擊橫流會(huì)產(chǎn)生沿射流軸線運(yùn)動(dòng)的螺旋渦，并與渦環(huán)疊加. 所以當(dāng)射流沖擊橫流時(shí)，在向壁面發(fā)展過(guò)程中射流邊界存在渦旋的剪切層，在此區(qū)域射流與環(huán)境流體產(chǎn)生卷吸和混合. 射流到達(dá)壁面形成對(duì)底壁的沖擊，在沖擊區(qū)內(nèi)流線彎曲并折轉(zhuǎn)過(guò)渡到平行于壁面的壁射流流動(dòng)形態(tài)，在這個(gè)過(guò)程射流與橫流、固體壁面之間產(chǎn)生復(fù)雜的卷吸作用，從而在流動(dòng)中產(chǎn)生了射流剪切層、橫流繞流等復(fù)雜流動(dòng)特征[12]. 環(huán)境橫流遇到射流的阻礙形成繞流，射流由于前后邊界存在壓力差發(fā)生偏轉(zhuǎn)，對(duì)環(huán)境水體為半無(wú)限水深情況，在水域的上下會(huì)存在壓力差，導(dǎo)致產(chǎn)生渦旋.

下面利用RNGk-湍流模型[9,13]對(duì)射流沖擊橫流的結(jié)構(gòu)進(jìn)行分析. 圖1為半無(wú)限橫流紊動(dòng)射流示意圖.

圖1 半無(wú)限橫流紊動(dòng)射流示意圖

1)在出口邊界x=x2滿(mǎn)足

2)在底部邊界z=0滿(mǎn)足

a.當(dāng)x2+y2≤D2/4時(shí)，u=v=0,ω=uj,φ=c0;

3)在上邊界z=z2滿(mǎn)足

4)在對(duì)稱(chēng)面y=0，v=0滿(mǎn)足

5)在側(cè)平面y=y2滿(mǎn)足v=0,ω=0,u=u0,且

由RNGk-ε湍流模型的公式分析可知，在不同邊界面上(例如出口邊界、底部邊界、上邊界等)渦旋滿(mǎn)足的條件不相同，說(shuō)明不同的邊界面會(huì)產(chǎn)生不同的渦旋結(jié)構(gòu). 而模型中的邊界面是相對(duì)射流與橫流的沖擊起點(diǎn)而言的，如果將起點(diǎn)看做是射流的射水口，由此可知離射水口不同的位置會(huì)有不同的渦旋結(jié)構(gòu)，即射水口的位置可以影響渦旋參量. 除此之外公式中涉及到了橫流速度u0，射流速度uj，因此可知橫流、射流速度可以影響渦旋參量.

1.2 3D-PIV技術(shù)原理

采用粒子圖像測(cè)速(PIV)技術(shù)[14]獲得速度場(chǎng)信息，進(jìn)而根據(jù)速度場(chǎng)分析渦旋的速度環(huán)量、渦量場(chǎng)和黏度，可以避免測(cè)試裝置對(duì)裝置內(nèi)原始流場(chǎng)的影響，還可以對(duì)空間瞬態(tài)流場(chǎng)進(jìn)行實(shí)時(shí)測(cè)量，深入研究裝置內(nèi)液體流動(dòng)特性. 為了實(shí)現(xiàn)三維粒子圖像測(cè)速(3D-PIV)[15]，對(duì)左右相機(jī)各自連續(xù)2幀粒子圖像進(jìn)行運(yùn)動(dòng)分析. 首先對(duì)左右相機(jī)同一時(shí)刻的粒子圖像作粒子匹配和融合后找到對(duì)應(yīng)區(qū)域，然后根據(jù)同一相機(jī)前后不同時(shí)刻對(duì)應(yīng)的粒子圖像作相關(guān)性分析獲得粒子在平面上的位移矢量信息，最后對(duì)矢量結(jié)果進(jìn)行三維重構(gòu)后得到運(yùn)動(dòng)粒子的空間位移.

在實(shí)驗(yàn)中，用CMOS相機(jī)替代PIV技術(shù)中常用的CCD相機(jī)，用細(xì)小泡沫替代PIV測(cè)速技術(shù)中的玻璃微珠作為示蹤粒子. 用穩(wěn)定光源照亮流場(chǎng)，并照像記錄粒子圖像. 由已知的時(shí)間間隔分析得到粒子位移，即可得到示蹤粒子在某點(diǎn)的速度矢量. 如果示蹤粒子在某個(gè)已知的時(shí)間間隔出內(nèi)，由位置(x1,y1,z1)運(yùn)動(dòng)至(x2,y2,z2)，則由示蹤粒子所在處的某流體質(zhì)點(diǎn)三維速度分量u,v,w為

當(dāng)Δt→0時(shí)，即

PIV記錄粒子的Lagrange速度，當(dāng)粒子對(duì)流體的跟隨性非常好時(shí)，可認(rèn)為此速度代表流體的Euler速度. 即PIV通過(guò)測(cè)得示蹤粒子的速度矢量來(lái)定量顯示整個(gè)流場(chǎng)的流動(dòng)特征.

2 理論模擬

2.1 Autodesk CFD模擬方法

CFD計(jì)算流體力學(xué)(Computational fluid dynamics,CFD)從計(jì)算方法出發(fā)，利用計(jì)算機(jī)得到流體控制方程的近似解.

模擬過(guò)程為：

a.用solid works軟件建立三維模型，該模型材料為PVC，厚度為8 mm，尺寸為800 mm×300 mm×400 mm，與實(shí)際的裝置相同.

b.將三維模型導(dǎo)入Autodesk CFD，設(shè)置好模型實(shí)驗(yàn)裝置的初始溫度、壓強(qiáng)、邊界條件(即橫流與射流的速度)以及運(yùn)動(dòng)的方向等，得出水流模擬的軌跡運(yùn)動(dòng)圖像. 最后分析圖像，找出流速比、水箱深度與射水口位置最佳參量.

CFD建立理想化的實(shí)驗(yàn)?zāi)Ｐ褪峭ㄟ^(guò)改變初始條件和邊界條件等參量即可得出不同情況下的模擬效果，這可降低探究實(shí)驗(yàn)成本，提高實(shí)際實(shí)驗(yàn)的可行性.

2.2 模擬實(shí)驗(yàn)參量設(shè)計(jì)

隨著流速比的增大，渦旋近區(qū)范圍內(nèi)射流向下游的偏轉(zhuǎn)程度有所減小，射流入水后射流主體兩側(cè)剪切渦間的相互作用對(duì)環(huán)境流體的卷吸作用明顯，更加容易形成渦旋[16]. 在模擬的過(guò)程中，結(jié)合實(shí)際情況可達(dá)到的最大射流速度為8.00 m/s. 首先在模擬裝置中裝入20 cm高的水，再?gòu)目趶綖榫匦蔚纳渌谠O(shè)置射流，在模擬過(guò)程中，保持射流方向與橫流方向一致，從左邊射向右邊. 模擬實(shí)驗(yàn)參量設(shè)置如表1所示.

表1 參量設(shè)置

部分結(jié)果模擬如圖2～3所示.三維坐標(biāo)圖中，不同顏色代表不同的速度，紅色框內(nèi)為形成的渦旋.

由圖2的紅色框內(nèi)渦旋的形態(tài)、速度大小和明顯程度可以看出：當(dāng)射水口接近水面、射流速度為8.00 m/s時(shí)，橫流速度越小，渦旋產(chǎn)生越明顯，當(dāng)橫流速度為0.05 m/s時(shí)可以產(chǎn)生巨大而且穩(wěn)定渦旋.

(a)橫流0.50 m/s

由圖3的紅色框內(nèi)渦旋的形態(tài)、速度大小和明顯程度可以看出：當(dāng)射水口靠近水底、射流速度為8.00 m/s時(shí)，橫流速度越小，渦旋產(chǎn)生越明顯，當(dāng)橫流速度為0.05 m/s時(shí)產(chǎn)生的渦旋更加明顯.

(a)橫流0.50 m/s

2.3 模擬結(jié)果分析

當(dāng)射流速度固定在8.00 m/s時(shí)，橫流速度越小，即流速比更大，更易產(chǎn)生巨大穩(wěn)定的渦旋，因此得到符合實(shí)驗(yàn)的流速為：橫流流速0.05 m/s，射流流速8.00 m/s. 同時(shí)比較射水口位置接近水面和靠近水底的三維圖像，可知射水口接近水面產(chǎn)生的渦旋在射水口兩側(cè)呈正立狀態(tài)，射水口靠近水底時(shí)，產(chǎn)生的渦旋在水底下呈臥鋪狀態(tài)，考慮到渦旋的可視化和測(cè)量，最終實(shí)驗(yàn)選擇射水口接近水面.

3 實(shí)驗(yàn)裝置設(shè)計(jì)

3.1 渦旋產(chǎn)生裝置

渦旋產(chǎn)生裝置如圖4所示，整體由水箱、隔板、水泵、水管、射水口和2個(gè)攝像頭組成. 水箱和隔板使用PVC材料，厚度為8 mm. 水箱尺寸為800 mm×300 mm×400 mm. 隔板尺寸為400 mm×300 mm. 隔板均分水箱，使整個(gè)水箱呈回型結(jié)構(gòu)，在盡量小的空間模擬環(huán)境橫流. 水管連接射水口和水泵，水泵置于水箱角落，射水口置于對(duì)角的另角落，避免水泵吸水造成干擾，射水口距離水箱底部5～10 cm.

圖4 渦旋產(chǎn)生裝置整體效果圖

3.2 渦旋測(cè)量裝置

3.2.1 測(cè)量裝置設(shè)計(jì)

使用紅色泡沫為標(biāo)記物，實(shí)驗(yàn)裝置周?chē)藗?cè)壁和頂部，其余部分均用黑色卡紙圍住，來(lái)避免外部光源干擾. 攝像頭使用普通USB攝像頭，攝像頭分辨率為1 280 pixel×720 pixel，幀數(shù)為20 s-1. 2個(gè)攝像頭在同一平面上并保持拍攝角度垂直，分別拍攝裝置側(cè)視圖和俯視圖. 拍攝俯視圖的攝像頭距水面高度10 cm，距離側(cè)壁15 cm. 拍攝側(cè)視圖的攝像頭水平位置為側(cè)壁與隔板的中心，即距離隔板7.5 cm. 2個(gè)攝像頭通過(guò)2個(gè)USB接口與同一電腦連接.

3.2.2 測(cè)量算法設(shè)計(jì)

使用opencv顏色識(shí)別算法提取泡沫顏色(紅色)區(qū)域，并且轉(zhuǎn)為灰度圖. 使用Sobel算子，在x方向求梯度并且進(jìn)行threshold二值化處理，接著dilate膨脹1次、再erode腐蝕1次，最后再dilate膨脹1次，獲得特征明顯的細(xì)小泡沫區(qū)域圖像.

通過(guò)findContours標(biāo)記白色區(qū)域在圖像中的像素點(diǎn)位置，并根據(jù)比例換算成實(shí)際空間坐標(biāo)點(diǎn). 根據(jù)2個(gè)攝像頭分別獲得的2組二維坐標(biāo)點(diǎn)合成三維坐標(biāo)點(diǎn)，并由三維坐標(biāo)點(diǎn)判定渦旋中心坐標(biāo). 可由連續(xù)2幀圖像所獲得的三維坐標(biāo)點(diǎn)計(jì)算速度向量場(chǎng)v. 計(jì)算某點(diǎn)到渦旋中心距離，并近似認(rèn)為半徑a. 根據(jù)速度場(chǎng)和以a為半徑的圓計(jì)算速度環(huán)量Γ. 根據(jù)a、速度場(chǎng)v計(jì)算渦量場(chǎng)Ω. 再根據(jù)渦量場(chǎng)Ω計(jì)算時(shí)間t的偏導(dǎo)以及拉普拉斯算子▽2計(jì)算渦旋的擴(kuò)散系數(shù)γ.

測(cè)量算法流程圖如圖5所示.

圖5 測(cè)量算法流程圖

4 實(shí)驗(yàn)操作流程

實(shí)驗(yàn)裝置如圖6所示.

圖6 實(shí)際裝置圖

實(shí)驗(yàn)操作流程分為2部分：1)放入細(xì)小紅色泡沫，用流速較大的射流沖擊橫流來(lái)產(chǎn)生可視化的巨大穩(wěn)定渦旋；2)對(duì)巨大穩(wěn)定渦旋進(jìn)行相關(guān)測(cè)量.

4.1 射流沖擊橫流產(chǎn)生渦旋實(shí)驗(yàn)步驟

向水箱中注入20 cm高的水，加入用于可視化的紅色泡沫.

打開(kāi)水泵，射水口開(kāi)始噴水，射水口噴射的水流速約為8.00 m/s. 水泵抽水和射水口所射出的水流會(huì)帶動(dòng)裝置中的靜水，最終使水箱中的水流速穩(wěn)定在0.10 m/s，射水口噴水之后，再微調(diào)節(jié)射水口的深度和方向，直至等待一段時(shí)間后生成明顯且穩(wěn)定的大型渦旋，效果圖如圖7所示.

(a)

4.2 渦旋結(jié)構(gòu)測(cè)量操作流程

打開(kāi)2個(gè)USB攝像頭，并連接到筆記本電腦上，根據(jù)攝像頭拍攝情況來(lái)調(diào)整攝像頭的焦距和位置，使2個(gè)攝像頭的拍攝的圖像清晰、互相垂直且寬度相等.

在生成渦旋且紅色泡沫分布在渦旋附近時(shí)，點(diǎn)擊運(yùn)行渦旋測(cè)量程序，系統(tǒng)將自動(dòng)處理該時(shí)刻的圖像并獲取速度環(huán)量場(chǎng)、線速度場(chǎng)、渦量場(chǎng)和黏度分布的數(shù)據(jù)和三維可視化圖像，電腦的測(cè)量效果圖如圖8～9所示.

圖8 帶有紅色泡沫的渦旋原圖

圖9 坐標(biāo)點(diǎn)標(biāo)記圖

5 實(shí)驗(yàn)數(shù)據(jù)結(jié)果與分析

經(jīng)過(guò)設(shè)計(jì)的算法，系統(tǒng)自動(dòng)對(duì)渦旋的坐標(biāo)點(diǎn)進(jìn)行處理，得到渦旋的參量.

圖10是速度環(huán)量場(chǎng). 速度環(huán)量是標(biāo)量，有正負(fù)號(hào)，規(guī)定沿曲線順時(shí)針繞行的方向?yàn)檎较?，沿曲線逆時(shí)針繞行的方向?yàn)樨?fù)方向. 圖10中可看出渦旋各個(gè)標(biāo)記點(diǎn)的速度環(huán)量，其中箭頭表示方向，大部分的標(biāo)記點(diǎn)沿著順時(shí)針?lè)较蜣D(zhuǎn)動(dòng)，方向向右，有小部分的標(biāo)記點(diǎn)逆時(shí)針旋轉(zhuǎn)，方向向左. 線段的長(zhǎng)度表示速度環(huán)量的大小，由圖10可知當(dāng)標(biāo)記點(diǎn)的半徑越大時(shí)速度越大. 因?yàn)閷?shí)驗(yàn)過(guò)程中用紅色泡沫為標(biāo)記點(diǎn)，標(biāo)記點(diǎn)的運(yùn)動(dòng)是隨機(jī)的，不一定總是沿著渦旋結(jié)構(gòu)運(yùn)動(dòng)，所以導(dǎo)致標(biāo)記點(diǎn)的運(yùn)動(dòng)方向并不統(tǒng)一，出現(xiàn)圖10中有部分標(biāo)記點(diǎn)運(yùn)動(dòng)方向相反的情況，說(shuō)明測(cè)量結(jié)果的方向符合實(shí)際實(shí)驗(yàn). 由于流體的黏性、斜壓性以及無(wú)勢(shì)外力將使速度環(huán)量隨著渦旋的深入有數(shù)值上的減小，即速度環(huán)量隨著曲線半徑的增大而增大，因此圖10的速度大小符合理論和實(shí)驗(yàn)結(jié)果，測(cè)量準(zhǔn)確.

圖10 速度環(huán)量場(chǎng)

圖11是線速度場(chǎng). 其方向沿運(yùn)動(dòng)軌道的切線方向，線段長(zhǎng)度表示速度大小，可以表示渦旋標(biāo)記點(diǎn)運(yùn)動(dòng)的快慢和方向. 由圖11可以看出箭頭方向沿著切線方向，渦旋外側(cè)速度大于內(nèi)部速度，這與圖10中測(cè)量的速度結(jié)果相一致.

圖11 線速度場(chǎng)

圖12是渦量場(chǎng). 圖中線段長(zhǎng)度表示渦量的大小，等于渦旋各點(diǎn)繞渦旋中心旋轉(zhuǎn)的平均角速度的2倍，箭頭方向表示環(huán)量方向，與渦旋的瞬時(shí)轉(zhuǎn)動(dòng)軸線重合[17]. 由圖12可以看出渦旋各點(diǎn)的渦量方向都是豎直向上的，在實(shí)驗(yàn)中如果紅色泡沫繞著渦旋做圓周運(yùn)動(dòng)，那么標(biāo)記點(diǎn)的角速度方向用右手螺旋定則可知與轉(zhuǎn)動(dòng)軸的方向相同，即豎直向上. 由圖12中還可以觀察到渦旋中心軸的渦量最大，遠(yuǎn)離中心軸渦量變小，在實(shí)驗(yàn)觀察中可知渦旋中心軸上的角速度最大，遠(yuǎn)離中心軸時(shí)角速度變小，因此參量結(jié)果符合實(shí)際實(shí)驗(yàn).

圖13是渦旋的黏度的測(cè)量結(jié)果. 實(shí)驗(yàn)中采用的液體是水，查表可知水的黏度在常溫(取24.9 ℃)的標(biāo)準(zhǔn)值為[18]1 mm2·s-1，如圖13紅色線所畫(huà)的數(shù)值. 實(shí)驗(yàn)中每個(gè)標(biāo)記點(diǎn)的黏度的值由圖13可知大部分集中在紅色線附近，計(jì)算誤差允許范圍內(nèi)的標(biāo)記點(diǎn)的平均黏度為0.898 2 mm2·s-1，與水的標(biāo)準(zhǔn)值大致相同，實(shí)驗(yàn)測(cè)量結(jié)果接近標(biāo)準(zhǔn)值.

圖12 渦量場(chǎng)

圖13 黏度

6 結(jié) 論

利用Autodesk CFD軟件對(duì)實(shí)驗(yàn)參量進(jìn)行模擬，探究了橫流和射流的流速比以及射水口的位置對(duì)實(shí)驗(yàn)的影響，并且最終確定當(dāng)射水口的位置接近水面，橫流速度為0.05 m/s，射流速度為8.00 m/s時(shí)的能產(chǎn)生巨大穩(wěn)定的渦旋. 使用PVC材料搭建了尺寸為800 mm×300 mm×400 mm實(shí)驗(yàn)裝置，利用射流沖擊橫流的原理產(chǎn)生渦旋，并在3D-PIV技術(shù)的基礎(chǔ)上設(shè)計(jì)算法測(cè)量渦旋的三維結(jié)構(gòu)，測(cè)量出渦旋的速度環(huán)量場(chǎng)、線速度場(chǎng)、渦量場(chǎng)以及黏度，分析發(fā)現(xiàn)速度環(huán)量場(chǎng)、渦量場(chǎng)符合理論和實(shí)際實(shí)驗(yàn)的情況，線速度與速度環(huán)量場(chǎng)結(jié)果一致，測(cè)得粒子的平均黏度為0.898 2 mm2·s-1，與水的標(biāo)準(zhǔn)值1 mm2·s-1大致相同.