霍夢(mèng)佳，鞏明鑫，王國(guó)亮，吳宏雷

（中國(guó)電建集團(tuán)山東電力建設(shè)有限公司，山東濟(jì)南250001）

隨著船舶行業(yè)的持續(xù)高位運(yùn)行，世界各國(guó)大型干船塢不斷興建，而大船塢排水配套設(shè)施成為影響船塢性能的主要因素之一。在船廠項(xiàng)目中，船塢是最重要的生產(chǎn)設(shè)施，排水系統(tǒng)的合理設(shè)計(jì)及優(yōu)化是對(duì)船塢的生產(chǎn)過(guò)程至關(guān)重要。

沙特國(guó)王港項(xiàng)目是在建的世界上規(guī)模最大的“超級(jí)船廠”，項(xiàng)目新建了1座10萬(wàn)t級(jí)造船塢，其基本尺寸為：塢長(zhǎng)374 m，塢寬90 m，塢底高程-10.8 m，屬于大型船塢。船塢排水系統(tǒng)配備3臺(tái)立式混流泵，本文借助CFD技術(shù)對(duì)船塢的排水系統(tǒng)進(jìn)行模擬計(jì)算，預(yù)測(cè)排水管路及泵入口處的流型，檢查是否有渦流存在，并對(duì)排水系統(tǒng)性能進(jìn)行分析和評(píng)價(jià)。

1 模擬方法

1.1 計(jì)算軟件及方法

本模擬研究在商業(yè)模擬軟件ANSYS CFX-19.0中進(jìn)行，采用雷諾平均Navier-Stokes方程求解三維定常、不可壓縮、黏性流體模型[1]。利用有限體積法（FVM）對(duì)控制方程進(jìn)行離散，并用簡(jiǎn)單算法進(jìn)行求解。為了得到高質(zhì)量的計(jì)算結(jié)果，此模型共劃分了2 300萬(wàn)個(gè)網(wǎng)格單元。

1.2 建立模型

建立排水系統(tǒng)的模型如圖1所示，包括干船塢（Dry Dock）、地下流道（Culvert）、導(dǎo)流片（Vanes）、防渦流裝置（Anti-vortex Devices）和泵室（Pump bay）等。

圖1 干船塢泵房計(jì)算模型

在泵座上安裝通風(fēng)口以排除空氣，并將圖2所示的灰色部分采用混凝土填充，因?yàn)槿绻肟诤聿可戏娇臻g過(guò)大，在較大的流速差異下，在泵入口處容易產(chǎn)生旋流，將對(duì)排水性能產(chǎn)生不利影響。

圖2 泵倉(cāng)尺寸

模型的入口邊界設(shè)置在干船塢的自由水面，假定流量均勻；出口邊界位于泵上方的排出管處，距泵的吸入口有足夠的距離。

根據(jù)水泵流量，將入口邊界條件設(shè)置為統(tǒng)一的指定速度，出口邊界條件之一設(shè)置為壓力指定條件，其他邊界分別根據(jù)泵的流量設(shè)置為速度指定條件[2]。

在模型中，干船塢的空氣—水自由面采用“自由滑移墻邊界（Free slip wall boundary）”，這是在水面變形很小且可以忽略不計(jì)時(shí)的常規(guī)簡(jiǎn)化方法[3，4]。

1.3 計(jì)算實(shí)例

操作條件及計(jì)算實(shí)例，詳見(jiàn)表1。

表1 操作條件及計(jì)算實(shí)例

1.4 計(jì)算結(jié)果的檢查方法

根據(jù)工程經(jīng)驗(yàn)和美國(guó)國(guó)家泵入口設(shè)計(jì)標(biāo)準(zhǔn)（ANSI/HI 9.8-2012）進(jìn)行模型測(cè)試，通過(guò)以下方法檢查CFD計(jì)算結(jié)果。

（1）來(lái)自船塢流道的水流量方向應(yīng)指向水泵入口，并在到達(dá)泵入口前實(shí)現(xiàn)最小渦流。

（2）在泵入口附近未觀察到明顯的導(dǎo)致旋渦發(fā)生的回流。

（3）在吸入管中未發(fā)現(xiàn)明顯的旋流，面積平均旋流角表示的預(yù)旋強(qiáng)度應(yīng)小于5°。其中，旋流角θ定義如下：

式中：Vθ和Vα分別指的是泵軸線的切線速度和軸向速度（m/s）；A指的是吸入管的橫截面積（m2）。根據(jù)圖3所示的泵房剖面，檢查泵室的上述流量標(biāo)準(zhǔn)；在喉部的吸入口處檢查泵吸入管中的流體旋轉(zhuǎn)強(qiáng)度，即旋流角。圖3中各參數(shù)代表的含義如下：ZS為吸水口喉部橫截面，以此截面檢查預(yù)旋強(qiáng)度；Z1為喇叭口和泵室底部之間的中段平面，以此截面檢查是否有渦流在泵室底部發(fā)生；Y為管道中心的垂直平面，用于檢查流動(dòng)均勻性和速度分布。

圖3 泵房剖面檢查

2 結(jié)果分析

2.1 實(shí)例1

在此例中，所有的泵都處于運(yùn)行狀態(tài)，計(jì)算結(jié)果如圖4—12所示。

圖4—5顯示了喇叭口和泵室底部之間的中段面Z1處的速度云圖和矢量分布。由圖4—5可以看出，每個(gè)泵室的流量分布幾乎均勻，在吸入口處未觀察到明顯的回流。由于流動(dòng)的均勻性和抗渦流裝置的有效性，沒(méi)有渦流進(jìn)入泵中。

圖4 截面Z1的速度云圖

圖5 截面Z1的速度矢量

圖6—8顯示了截面YA、YB和YC處的速度矢量分布。由圖6—8可以看出，每個(gè)泵室中的流量分布均勻，在吸入口附近沒(méi)有觀察到明顯的回流。

圖6 截面YA的速度矢量

圖7 截面YB的速度矢量

圖8 截面YC的速度矢量

圖9顯示了從泵室底部和側(cè)壁吸入口附近的流線分布，可以看出泵961-G-4421A、泵961-G-4421B和泵961-G-4421C的喇叭口處無(wú)流線回流。

圖9 泵倉(cāng)底部的流線分布

圖10—12顯示了進(jìn)水口喉部的軸向速度、二次流和旋流角分布。從圖10—12可以看出，進(jìn)場(chǎng)流量被均勻地吸入，面積平均旋流角小于先前的接受標(biāo)準(zhǔn)5°。

圖10 泵961-G-4421A截面ZS處的軸速度云圖、二次流及旋流角分布

2.2 實(shí)例2

在此例中，所有的泵都處于運(yùn)行狀態(tài)，計(jì)算結(jié)果如圖13—21所示。

圖13—14顯示了喇叭口和泵室底部之間的中剖面Z1處的速度云圖和矢量分布。由圖13—14可以看出，每個(gè)泵室的流量分布幾乎均勻，在吸入口處未觀察到明顯的回流。由于流動(dòng)的均勻性和抗渦流裝置的有效性，沒(méi)有渦流進(jìn)入泵中。

圖13 截面Z1的速度云圖

圖11 泵961-G-4421B截面ZS處的軸速度云圖、二次流及旋流角分布

圖12 泵961-G-4421C截面ZS處的軸速度云圖、二次流及旋流角分布

圖14 截面Z1的速度矢量

圖16 截面YB的速度矢量

圖17 截面YC的速度矢量

圖15—17顯示了截面YA、YB和YC處的速度矢量分布。由圖15—17可以看出，每個(gè)泵室中的流量分布均勻，在吸入口附近沒(méi)有觀察到明顯的回流。

圖15 截面YA的速度矢量

圖18顯示了從泵室底部和側(cè)壁吸入口附近的流線分布，可以看出泵961-G-4421A、泵961-G-4421B和泵961-G-4421C的喇叭口處無(wú)流線回流。

圖18 泵倉(cāng)底部的流線分布

圖19—21顯示了進(jìn)水口喉部的軸向速度、二次流和旋流角分布。從圖19—21可以看出，進(jìn)場(chǎng)流量被均勻地吸入，面積平均旋流角小于先前的接受標(biāo)準(zhǔn)5°。

2.3 實(shí)例3

在此例中，泵961-G-4421A處于運(yùn)行狀態(tài)，泵961-G-4421B和泵961-G-4421C處于待機(jī)狀態(tài)，計(jì)算結(jié)果如圖22—26所示。

圖20 泵961-G-4421B截面ZS處的軸速度云圖、二次流及旋流角分布

圖21 泵961-G-4421C截面ZS處的軸速度云圖、二次流及旋流角分布

圖22—23顯示了喇叭口和泵室底部之間的中剖面Z1處的速度云圖和矢量分布。由圖22—23可以看出，每個(gè)泵室的流量分布幾乎均勻，在吸入口處未觀察到明顯的回流。由于流動(dòng)的均勻性和抗渦流裝置的有效性，沒(méi)有渦流進(jìn)入泵961-G-4421A中。

圖22 截面Z1的速度云圖

圖23 截面Z1的速度矢量

圖24顯示了在截面YA上的速度矢量分布，可以看出每個(gè)泵室中的流量均勻，在吸入口附近沒(méi)有觀察到明顯的回流。圖25顯示了泵室底部和側(cè)壁的吸入口附近的流線，泵961-G-4421A的喇叭口處沒(méi)有流線回流。圖26顯示了吸入口喉部的軸向速度、二次流和旋流角分布，進(jìn)場(chǎng)流量被均勻地吸入，面積平均旋流角小于先前的接受程度標(biāo)準(zhǔn)5°。

圖24 截面YA的速度矢量

圖25 泵倉(cāng)底部的流線分布

3 結(jié)論

通過(guò)對(duì)船塢排水系統(tǒng)進(jìn)行建模和計(jì)算，分析了3種工況下的管路和泵吸入口處的流型和速度分布，得到了以下結(jié)論。

（1）在所有工況下，從泵艙側(cè)壁和底部都沒(méi)有觀察到明顯的回流現(xiàn)象，這表明導(dǎo)流片和防渦流裝置對(duì)避免渦流的產(chǎn)生是有效的。

（2）各泵在所有計(jì)算工況下的旋流角均小于5°且在可接受標(biāo)準(zhǔn)之內(nèi)，詳見(jiàn)表2。

表2 各計(jì)算實(shí)例下的吸入管口的旋流角

圖26 泵961-G-4421A截面ZS處的軸速度云圖、二次流及旋流角分布

A處船塢排水系統(tǒng)的性能對(duì)船塢的生產(chǎn)作業(yè)具有十分重要的作用，本文通過(guò)對(duì)沙特國(guó)王港項(xiàng)目的干船塢排水系統(tǒng)進(jìn)行CFD分析，得到了各操作條件下的排水系統(tǒng)水流流型及速度分布等，驗(yàn)證了防渦流裝置對(duì)避免渦流產(chǎn)生的必要性，可有效提高排水系統(tǒng)性能。

本文采用的模型設(shè)置、計(jì)算方法及研究思路將對(duì)國(guó)內(nèi)外的工程有一定的指導(dǎo)意義。