渦旋流堵漏槳葉位置與堵漏效果關(guān)系機理

韓云東， 高占勝， 楊常青， 謝田華， 杜蓬杉

(海軍大連艦艇學(xué)院 航海系， 遼寧 大連 116018)

渦旋流堵漏槳葉位置與堵漏效果關(guān)系機理

韓云東， 高占勝， 楊常青， 謝田華， 杜蓬杉

(海軍大連艦艇學(xué)院 航海系， 遼寧 大連 116018)

為克服現(xiàn)有堵漏手段的不足，提出一種新技術(shù)方法——船舶渦旋流堵漏方法。根據(jù)物理試驗的具體情況，在利用Gambit軟件完成槳葉實體和計算域模型建模的基礎(chǔ)上，通過Fluent軟件完成剛性槳葉生成渦旋流的數(shù)值仿真，根據(jù)仿真得到的渦旋流場速度矢量和壓強分布，對試驗得出的槳葉工作位置對破口進水量的影響規(guī)律進行理論剖析，進而從運行機理上揭示槳葉位置與渦旋流形成效果的關(guān)系。

船舶;艦船工程；損害管制；計算流體動力學(xué)；堵漏；渦旋流；槳葉

船舶渦旋流堵漏技術(shù)是指通過專用設(shè)備在船舶破口外側(cè)形成局部渦旋流場，使破口處的海水高速旋轉(zhuǎn)，大幅減小海水的壓力，明顯降低(甚至是消除)破口進水速度，同時輔以傳統(tǒng)的堵漏方式，從而使堵漏的實施難度得到明顯降低，進而實現(xiàn)便捷、快速的堵漏，達到提高艦船堵漏效率的目的。[1-5]

金良安等[1]對船舶渦旋流堵漏的思想、基本原理及柔性槳葉和剛性槳葉生成的渦旋流對破口進水量的影響等進行初步研究；韓云東等[5]對剛性槳葉位置與渦旋流形成效果的關(guān)系進行試驗，定量研究槳葉工作位置對破口進水量的影響規(guī)律，揭示槳葉位置對渦旋流堵漏效果的影響。文獻[2]～文獻[4]從試驗數(shù)據(jù)分析的角度研究渦旋流堵漏方法，但缺少渦旋流運行機理方面的分析。本文針對上述研究成果，利用計算流體動力學(xué)(Computation Fluid Dynamics, CFD)技術(shù)對剛性槳葉生成的堵漏渦旋流場進行數(shù)值仿真，根據(jù)仿真結(jié)果進一步從理論層面對渦旋流形成位置和艦船堵漏效果的影響規(guī)律進行機理分析。

1 堵漏渦旋流的物理試驗

1.1 試驗裝置

采用專用的試驗裝置(見圖1)進行物理試驗。將驅(qū)動槳葉旋轉(zhuǎn)的電機的轉(zhuǎn)速及剛性槳葉與破口在水平方向上的距離固定，改變槳葉與破口在垂直方向上的距離，即在破口附近不同高度上裝置相同轉(zhuǎn)速的槳葉，形成高速渦旋流場。在有渦旋流的情況下測量到達模擬的破損艙室內(nèi)指定液位所用的時間。分析比較所用時間的變化，得到剛性槳葉工作位置與渦旋流堵漏效果之間的關(guān)系。

圖1 試驗裝置示意

葉片可調(diào)式渦旋流生成器通過特定的槳葉在水中生成所需的渦旋流，由驅(qū)動電機、軸桿和剛性槳葉組成。剛性槳葉是在直徑為1.5 cm，長度為8 cm的管軸上，沿管軸圓周方向均勻排列4片8 cm×3 cm的鋼制葉片(見圖2)。轉(zhuǎn)動時依靠槳葉與水的相互作用生成渦旋流。

圖2 剛性槳葉

1.2 試驗方法

測定試驗水池水深為53 cm，破口距水面高度17 cm，將槳葉浸深H設(shè)定為9 cm，依次調(diào)整為11 cm，13 cm，15 cm和17 cm。當(dāng)H=13 cm時，槳葉中心高度與破口中心高度平齊。試驗以“破艙進水時間”為指標(biāo)，按照指定的試驗方法，分別測量剛性槳葉入水深度為9～17 cm時在不同破口距離下的“破艙進水時間”值。

1.3 試驗結(jié)果

首先測得靜水下的進水時間T=40 s，然后分別測定剛性槳葉浸深H=9 cm，11 cm，13 cm，15 cm和17 cm時，在不同破口距離L值下的進水時間T值。每組試驗重復(fù)進行3次，取其平均值，測量結(jié)果見表1。

由試驗數(shù)據(jù)可知：當(dāng)剛性槳葉與破口處在同一高度附近時，距離破口越近，生成的渦旋流對延長堵漏時間、降低破口進水量的效果越明顯。

2 堵漏渦旋流的CFD仿真

表1 剛性槳葉不同位置時進水時間的測量值 s

CFD方法是以黏性流理論為基礎(chǔ)，以數(shù)值計算方法為手段，通過數(shù)值仿真試驗獲取相應(yīng)水動力的一種方法。該方法已成為模型試驗的一種輔助手段，并正逐漸成為與模型試驗同等重要的試驗流體力學(xué)的研究方法。[6-8]

2.1 槳葉模型建立

根據(jù)“1.1”節(jié)中的試驗?zāi)Ｐ统叽纾瑒傂詷~是由4片8 cm×3 cm的鋼片均勻地焊接在鋼管上形成的，其結(jié)構(gòu)見圖3，其中：D為剛性槳葉的直徑；di為槳轂內(nèi)徑；do為槳轂外徑；T為槳葉徑向長度；L為槳葉寬度；d為槳葉厚度。

a)側(cè)視圖b)俯視圖

圖3 剛性槳葉結(jié)構(gòu)示意

剛性槳葉的參數(shù)：槳轂內(nèi)徑di=1 cm；槳轂外徑do=2 cm；槳葉徑向長度T=3 cm；槳葉寬度L=8 cm；槳葉厚度d=0.3 cm。對于剛性槳葉需在槳葉固定器上加裝的葉片，也采用相同規(guī)格、厚度d=1.5 mm的白鋼片制作。首先運用軟件Gambit進行槳葉實體建模，包括中心槳軸和4個葉片；然后在槳葉上劃分非結(jié)構(gòu)網(wǎng)格，槳葉網(wǎng)格分布見圖4。

2.2 計算域模型建立

圖4 槳葉網(wǎng)格分布

為實現(xiàn)槳葉旋轉(zhuǎn)過程的動態(tài)數(shù)值模擬，將整個流體計算域分為旋轉(zhuǎn)動域和外部靜域。在數(shù)值模擬過程中，旋轉(zhuǎn)動域隨槳葉一起運動，從而實現(xiàn)CFD計算的動邊界模擬。在旋轉(zhuǎn)動域和外部靜域的交界處設(shè)置Interface邊界，實現(xiàn)內(nèi)外流域的流場數(shù)據(jù)傳遞和交換。同時，為保證計算邊界對槳葉旋轉(zhuǎn)沒有影響，該計算沿x軸和y軸到槳葉的距離為5D(槳葉徑向最大半徑)，沿z軸到槳葉的距離為3h(槳葉軸向高度)，見圖5。

圖5 堵漏渦旋流計算域分區(qū)

分別在旋轉(zhuǎn)動域和外部靜域上劃分三維非結(jié)構(gòu)四面體網(wǎng)格。在槳葉附近采用0.001 m的網(wǎng)格尺寸進行近壁面網(wǎng)格劃分，在此基礎(chǔ)上以一定的比例逐漸向外進行網(wǎng)格尺度遞增，從而在保證近壁面網(wǎng)格數(shù)量足夠的前提下最大限度地減少計算域內(nèi)總體網(wǎng)格的數(shù)目，提高數(shù)值計算效率。最終計算域共包含四面體非結(jié)構(gòu)網(wǎng)格117萬個，最差網(wǎng)格的扭曲度為0.807(滿足三維數(shù)值計算要求)。計算域總體情況見圖6。

2.3 渦旋流仿真

將建好的槳葉旋轉(zhuǎn)模型導(dǎo)入到CFD數(shù)值求解軟件Fluent中，進行數(shù)值計算。計算模型和相關(guān)參數(shù)的具體設(shè)置：

1) 利用1st-order implicit格式的Unsteady非定常數(shù)值計算模型。

2) 基于單元中心的Green-Gauss格式計算流場變量梯度。

圖6 計算域總體情況

3) 建立基于k-ε的兩方程湍流模型，使用Standard Wall Functions進行近壁面處理。

4) 設(shè)置槳葉和旋轉(zhuǎn)動域以400 r/min的轉(zhuǎn)速旋轉(zhuǎn)，根據(jù)Moving Mesh模型處理移動邊界問題，在動域和靜域交界面處設(shè)置Interface邊界條件，實現(xiàn)移動的交錯網(wǎng)格之間的數(shù)值傳遞。

5) 將計算域上部邊界設(shè)置為Pressure-Outlet壓力出口邊界條件，其他外邊界設(shè)置為Symmetry對稱邊界條件。

6) 槳葉設(shè)置為wall壁面邊界條件；流場壓強和速度采用PISO的壓強修正法計算，計算格式為PRESTO；其他計算均采用二階計算格式。

通過仿真就能得到剛性槳葉生成船舶堵漏渦旋流場的速度矢量和壓強分布情況。

3 仿真結(jié)果分析

針對得到的堵漏渦旋流的CFD仿真結(jié)果，從槳葉軸向中心位置剖面和槳葉垂向剖面，分別對仿真生成的渦旋流場的速度矢量和壓強分布進行分析。

3.1 槳葉軸向中心位置剖面的渦旋流場分析

圖7為槳葉軸向中心位置剖面的渦旋流場速度矢量分布。從圖7中可看出，在槳葉旋轉(zhuǎn)過程中，近壁面流體質(zhì)點會隨槳葉一起旋轉(zhuǎn)，形成渦旋流。在葉輪區(qū)內(nèi)部，流體質(zhì)點速度隨到槳葉中心軸的距離的增加而逐漸增大，并在槳葉邊緣處取得最大值(見圖7中靠近槳葉的區(qū)域)。此外，在主渦旋流區(qū)，流體質(zhì)點同樣作旋轉(zhuǎn)運動，并隨到槳葉邊緣距離的增加而呈遞減趨勢。

圖8為槳葉軸向中心位置剖面的渦旋流場壓力分布。從圖8中可看出，在槳葉旋轉(zhuǎn)過程中，高速旋轉(zhuǎn)的渦旋流生成器驅(qū)動水體旋轉(zhuǎn)，形成葉輪區(qū)，在葉輪區(qū)流體的誘導(dǎo)下，主渦旋流區(qū)的流體形成接近于二維流動條件的平面自由渦。葉輪區(qū)的流體在徑向速度的作用下被不斷甩出葉輪區(qū)，從而使葉輪區(qū)軸向中心形成一個低壓的壞境。近壁面流體質(zhì)點會隨槳葉一起旋轉(zhuǎn)，形成渦旋流。在葉輪區(qū)和主渦旋流區(qū)，壓強隨到槳葉中心軸的距離的增加而逐漸增大，但壓強梯度呈逐漸變小的趨勢。

圖7 槳葉軸向中心位置剖面的渦旋流場速度矢量分布

圖8 槳葉軸向中心位置剖面的渦旋流場壓力分布

3.2 槳葉垂向中心位置剖面的渦旋流場分析

圖9為槳葉垂向剖面的渦旋流場速度矢量分布。從圖9中可看出，槳葉產(chǎn)生的渦旋流分布在槳葉周圍，整個水體在槳葉的驅(qū)動下發(fā)生劇烈的湍流運動，宏觀上符合渦旋流的運動形態(tài)。在葉輪區(qū)，流體的速度矢量呈黃綠色，說明該區(qū)域的水體運動矢量最大；在主渦旋流區(qū)，靠近槳葉處流體的速度矢量呈綠色，說明該區(qū)域的水體運動矢量較大，而隨著遠(yuǎn)離槳葉，速度矢量呈淡藍(lán)色，即速度矢量逐漸減小；在上漩內(nèi)渦區(qū)和下漩內(nèi)渦區(qū)(渦旋流的上部和下部)，流體的速度矢量均呈藍(lán)色，該區(qū)域是速度矢量最小的區(qū)域，說明下漩內(nèi)渦區(qū)和上漩內(nèi)渦區(qū)(即旋轉(zhuǎn)槳葉的上下端)是渦旋流速度較小的區(qū)域，即在上下端部呈現(xiàn)三維效應(yīng)，旋轉(zhuǎn)速度較中心位置變低。

圖9 槳葉垂向剖面的渦旋流場速度矢量分布

圖10為槳葉垂向剖面的渦旋流場壓力分布。從圖10中可看出，槳葉產(chǎn)生的渦旋流由于重力的作用，上端渦旋流影響的范圍比較廣，渦旋流影響的范圍自上向下逐漸減小，最下端渦旋流影響的范圍最小。

圖10 槳葉垂向剖面的渦旋流場壓力分布

4 槳葉位置對堵漏效果影響的機理分析

針對物理試驗得出的關(guān)系規(guī)律，結(jié)合上述渦旋流CFD數(shù)值仿真分析結(jié)果，對上述影響規(guī)律的運行機理進行分析。伯努利方程是理想流體穩(wěn)定流動時的動力學(xué)基本規(guī)律[9]，其數(shù)學(xué)表達式為

p+1/2ρv2+ρgH=恒量

(1)

式(1)表明：理想流體在穩(wěn)定流動時，在同一流管不同截面處，單位體積流體的動能、勢能與該處壓強之和都是相等的。[10]根據(jù)伯努利原理，破口處海水沿船舶外板水平切線方向的速度越大，其作用于破口的側(cè)壓力越小。因此，只要設(shè)法增大船舶破口處海水的速度，即可使其側(cè)壓力及破口進水速度減小，從而降低堵漏難度，延長堵漏時間。

由文獻[2]中對艦船渦旋流堵漏效果的影響因素的分析可知，p=p0-ρω2R4/(2r2)。當(dāng)槳葉半徑R一定時，破口處的流場壓強與破口至槳葉中心的距離r(即L)有關(guān)。r越小，p越小，即渦旋流堵漏的效果越好。此外，從CFD數(shù)值模擬的渦旋流速度矢量圖、壓強分布圖及局部放大圖(見圖11)中可看出，破口到槳葉的距離越小，流體旋轉(zhuǎn)速度越大，產(chǎn)生的向心力也就越大。

同時，破口到槳葉的距離越小，流場總壓強越小，到槳葉中心的壓強差越大，渦旋流堵漏的效果也就越明顯。

5 結(jié)束語

針對物理試驗中得出的槳葉位置與渦旋流形成效果之間關(guān)系的影響規(guī)律，利用CFD技術(shù)，在Fluent軟件中按照物理試驗裝置的尺寸完成了渦旋流仿

圖11 軸向中心位置剖面的渦旋流場速度矢量局部放大圖

真，得到了槳葉軸向中心剖面的渦旋流場速度矢量和壓強分布；結(jié)合伯努利方程，從理論層面對渦旋流形成位置與艦船堵漏效果的影響規(guī)律進行了機理分析，揭示了槳葉位置與渦旋流形成效果之間關(guān)系的影響規(guī)律的成因。研究成果可為船舶渦旋流堵漏方法的理論研究和裝備研發(fā)提供必要的理論基礎(chǔ)及數(shù)值依據(jù)。

[1] 金良安, 韓云東, 石侃, 等. 基于人工渦旋流的船舶堵漏新方法[J]. 中國航海, 2011, 34(2):22-26.

[2] 韓云東, 黃若波, 石侃, 等. 艦船堵漏技術(shù)現(xiàn)狀與對策研究[J] . 海軍大連艦艇學(xué)院學(xué)報, 2011, 34(2): 134-138.

[3] 韓云東，石侃，王涌. 艦船渦旋流堵漏方法的實驗研究[J]. 艦船科學(xué)技術(shù), 2011, 33(7):46-49.

[4] 韓云東, 金良安, 侯建軍, 等. 船舶渦旋流堵漏槳葉位置與渦旋流形成效果的關(guān)系研究[J]. 實驗流體力學(xué), 2012,26(1):25-29.

[5] 韓云東，金良安，郭模燦. 剛性槳葉工作位置對渦旋流堵漏效果的影響[J]. 中國造船, 2011,52(2):139-145.

[6] 黃連忠, 林虹兆, 馬冉祺, 等. 基于CFD的圓弧型風(fēng)帆氣動優(yōu)化[J]. 中國航海, 2016, 39(2): 101-105.

[7] MOUSAVI S M, ZAMANKHAN P, JAFARI A. Computer Simulations of Sodium Formate Solution in a Mixing Tank [J]. Communications in Nonlinear Science and Numerical Simulation, 2008(13)：380-399.

[8] JAHODA M, MOSTEK M, KUKUKOVA A, et al. CFD Modeling of Liquid Homogenization in Stirred Tanks with One and Two Impellers Using Large Eddy Simulation [J]. Trans IChemE, Part A, Chemical Engineering Research and Design, 2007, 85(A5): 616-625.

[9] 童秉綱, 尹協(xié)遠(yuǎn), 朱克勤. 渦運動理論[M]. 2版.合肥:中國科學(xué)技術(shù)大學(xué)出版社, 2009.

[10] JHON E F, JOSEPH B F. Fluid Mechanics with Eng-ineering Applications[M]. New York: McGraw-Hill Education, 2002.

MechanismofLeakStoppageBetweenBladePositionandEffectwithSpecialVorticalFlows

HANYundong，GAOZhansheng,YANGChangqing,XIETianhua,DUPengshan

(Department of Navigation, Dalian Naval Academy, Dalian 116018， China)

The leak stoppage method with special vortical flows is a new technological idea to fill the gaps of traditional leak stoppage methods for warships. The model of the rigid blade and the calculation domain based on physical tests with Gambit are constructed and the numerical simulation of the rigid blade produced vortex flow is conducted by Fluent. Based on the calculated vortex flow velocity vector and the pressure distribution, the relationship between the inflow at the break and the blade position are theoretically studied.

ship; naval engineering； damage control； CFD； leak stoppage； vortical flow； blade

2017-06-25

國防科研項目；2015年軍內(nèi)科學(xué)技術(shù)研究項目(2015-JDKYXM-DLJTXY-023)

韓云東(1977—)，男，山東昌邑人，副教授，博士，主要研究方向為船舶安全保障與防護、裝備保障等。E-mail:gyro-han@163.com

1000-4653(2017)04-0024-04

U661.2+3

A