水下管匯阻力系數(shù)和附加質(zhì)量系數(shù)的CFD計算

，， ，

(中國石油大學(xué)(北京) 機械與儲運工程學(xué)院，北京 102249)

·設(shè)計計算·

水下管匯阻力系數(shù)和附加質(zhì)量系數(shù)的CFD計算

田宇，羅曉蘭，戚昱，王麗男

(中國石油大學(xué)(北京) 機械與儲運工程學(xué)院，北京 102249)

深水管匯下放是水下生產(chǎn)系統(tǒng)安全工作的關(guān)鍵技術(shù)之一。針對深水管匯下放時受到的水動力，即阻尼力和附加質(zhì)量力進(jìn)行分析研究。通過CFD軟件Fluent建立管匯水動力計算模型和水下作業(yè)環(huán)境中繞流場的模型，計算管匯下放3個方向阻尼力及其阻力系數(shù)；利用動網(wǎng)格技術(shù)模擬管匯在流場中的加速運動，求得管匯在水中的變速運動帶動周圍的流體一起做不定常運動所形成的附加質(zhì)量系數(shù)。通過與現(xiàn)有DNV經(jīng)驗數(shù)據(jù)對比，計算結(jié)果誤差較小，從而驗證了此計算方法的正確性，為船-纜-體的耦合分析計算及下放安裝的仿真控制提供參考。

水下管匯；CFD；阻尼系數(shù)；附加質(zhì)量系數(shù)

隨著油氣資源的開發(fā)重點向海洋轉(zhuǎn)移，在深水中順利安裝管匯是其關(guān)鍵技術(shù)之一。水動力系數(shù)影響著管匯在水中的動力學(xué)參數(shù)，對管匯所受到的水動力進(jìn)行分析和計算，可以準(zhǔn)確計算下放結(jié)構(gòu)物的載荷，以擴(kuò)大適用船只的范圍，提高作業(yè)安全系數(shù)。管匯的水動力系數(shù)在模擬管匯下放，進(jìn)行船-纜-體的耦合分析時有著極為重要的作用[1]。在管匯的下放過程中，受到的水動力主要包括阻尼力和附加質(zhì)量力。通常對水下物體流體動力系數(shù)的求解主要通過勢流理論等流體力學(xué)知識，進(jìn)行理論經(jīng)驗計算[2]，或者通過水池模型試驗?zāi)M求得結(jié)果，但計算方法在經(jīng)濟(jì)性和可靠性上難以很好的結(jié)合。隨著計算機技術(shù)的發(fā)展，計算流體力學(xué)(CFD)成為新的計算方法，使用計算機數(shù)值模擬求解，在精度和成本上具有極大優(yōu)勢，很適合進(jìn)行相關(guān)的工程計算。本文采用CFD通用軟件Fluent進(jìn)行流場模擬，求解管匯的相關(guān)水動力，計算相關(guān)的阻尼系數(shù)和附加質(zhì)量系數(shù)。

1 管匯計算模型的建立

1.1 受力分析

管匯入水后的運動表現(xiàn)是剛體在六自由度上都可以運動，受力較為復(fù)雜，假設(shè)管匯的水下運動局限于二維垂直平面內(nèi)且是線性的，同時管匯無俯仰運動，其受力分析如圖1所示[3]。

圖1 管匯受力分析

管匯的受力平衡方程為

(1)

作用于水下管匯且對其運動產(chǎn)生影響的水動力主要包括阻尼力和附加質(zhì)量力，根據(jù)力的獨立作用原理，本文中不考慮其他力的影響，單獨對阻尼力和附加質(zhì)量力進(jìn)行分析計算。

1.2 幾何模型

實際管匯模型如圖2所示，在SolidWorks中建立模擬計算的管匯模型。

實際管匯模型的管道密布在中心部分，結(jié)構(gòu)復(fù)雜，為便于計算和建模，對管匯模型進(jìn)行簡化[4]。為了保證水動力特性相似，盡可能減少對主要受力特征的影響，建立的簡化模型各個方向外形相似，投影面積基本保持不變，簡化外部桿件，降低管匯高度，內(nèi)部走線密集的地方簡化為方箱，忽略了內(nèi)部的一些細(xì)微的湍流和旋渦，同時可能會對浮力造成損失，不過在網(wǎng)格劃分和模擬中可以極大地提高計算效率。簡化后的模型尺寸和原模型基本保持一致，長×寬×高為13.26 m×7.3 m×3.0 m，如圖3。

圖2 實際管匯模型

圖3 簡化后的管匯模型

1.3 流體域模型

根據(jù)模擬計算的需要和無界流場的概念，采用便于設(shè)定邊界條件的長方體狀的流體域模型，如圖4所示。管匯的坐標(biāo)原點和流體域模型的幾何中心重合，坐標(biāo)系取x軸為管匯長度方向，y軸為管匯寬度方向，z軸管匯高度方向。

流體域限制了整個模擬過程的邊界，需要選擇合適的尺寸來消除邊界對模擬結(jié)果的影響。一般來說，管匯模型在流體域模型的中心位置時，保證兩者邊界長度比大于5∶1，流體域邊界就不會在模擬試驗中產(chǎn)生影響。本次采用6∶1的邊界比，流體域尺寸為79.56 m×43.8 m×18.0 m。

圖4 流體域模型

1.4 流場網(wǎng)格

在對管匯進(jìn)行水動力分析計算之前，需要在Fluent的前處理軟件ICEM中對管匯及流體域進(jìn)行網(wǎng)格劃分，如圖5所示。

圖5 流場模型網(wǎng)格

由于管匯形狀的復(fù)雜性，整個模型的尺寸較大，網(wǎng)格生成困難，采用結(jié)構(gòu)化網(wǎng)格來劃分，可以減少網(wǎng)格數(shù)量，提高網(wǎng)格質(zhì)量，使水動力的計算快速穩(wěn)定。通過不斷改變參數(shù)，對于遠(yuǎn)離管匯部分的網(wǎng)格劃分較大尺寸，管匯則劃分較小尺寸的網(wǎng)格，這樣處理可以節(jié)約網(wǎng)格數(shù)量，整個流場劃分網(wǎng)格的數(shù)量保持在200萬，網(wǎng)格質(zhì)量保證在0.4以上。

1.5 湍流控制方程

Fluent軟件可以通過模擬管匯的水下繞流場完成虛擬的水池試驗，求解出管匯在繞流場作用下的受力情況?；镜挠嬎阍硎遣捎美字Z平均納維-斯托克斯方程(RANS)，對引入湍流模型后構(gòu)成的封閉方程組求解，得到湍流的時均值。

湍流模型主要有雷諾應(yīng)力模型和渦粘模型2大類，基本方程為

(2)

(3)

雷諾應(yīng)力模型直接構(gòu)建表示雷諾應(yīng)力的方程，與式(2)、(3)聯(lián)立求解，考慮了湍流的各向異性，但數(shù)值穩(wěn)定性一般。渦粘模型不直接處理雷諾應(yīng)力項，而是引入湍動粘度，把湍流應(yīng)力表示成湍流粘度的函數(shù)，整個計算的關(guān)鍵在于確定這種湍動粘度，主要有零方程和k-ε、k-ω二方程湍流模型及其變形。工程中多采用二方程模式，本文采用Realizablek-ε湍流模型[6]，k和ε的方程分別為

(4)

(5)

其中：

(6)

2 阻力系數(shù)

2.1 驗證計算方法的準(zhǔn)確性

管匯固定在定常流域內(nèi)，入水后的總阻力為

(7)

式中：u0為管匯的相對運動速度，m/s；ρ為流體(海水)密度，kg/m3；Ac為物體在來海流流向的投影面積，m2；CDS為穩(wěn)態(tài)海流作用下的總阻力系數(shù)，可利用流體力學(xué)軟件Fluent，模擬管匯的3個方向在水中相對運動得到總阻力Fc進(jìn)行求解。

挪威船級社提供的相關(guān)阻力系數(shù)的理論數(shù)據(jù)(DNV-RP-H103)[7]如圖6。本次研究采用的Fluent模擬加公式計算的方法，選擇其中L/D=2的長方體模型計算阻力系數(shù)，在1 m/s定常流的繞流中，長方體受力為459.09 N，代入公式求得阻力系數(shù)為0.89。DNV提供的理論數(shù)值為0.87，可以看出本次研究所采用的計算方法精確度較好，誤差較小，相對誤差約2.3%。

圖6 阻力系數(shù)的理論數(shù)據(jù)

2.2 邊界條件

在Fluent中的模擬計算，應(yīng)盡可能貼近實際情況，并且要快速準(zhǔn)確。其基本參數(shù)設(shè)置如下：

1) 為了計算穩(wěn)定與節(jié)省計算時間，將求解問題設(shè)為穩(wěn)態(tài)(Steady)。

2) 湍流模型采用Realizablek-ε，壁面函數(shù)采用非平衡(Non-Equilibrium)。

3) 流域材料定義為4 ℃的海水(密度為1 028 kg/m3)。

4) 計算過程中動量方程的離散格式，設(shè)為二階迎風(fēng)格式(Second Order Upwind)，壓力差值格式設(shè)為PRESTO!。

5) 流域邊界入口處選擇速度入口(velocity-inlet)給定速度大小、方向和湍流強度等參數(shù)。

6) 流域邊界出口處選擇自由出流邊界(outflow)。

7) 流域邊界壁面選無滑移的壁面(wall)和對稱面(sym)[8]。

8) 固壁(管匯表面)選無滑移的壁面(wall)。

2.3 阻力系數(shù)計算

假設(shè)管匯固定且始終保持水平姿態(tài)，坐標(biāo)系的x軸、y軸、z軸分別對應(yīng)于管匯的長、寬和高方向，參照一般海流流速范圍，先后在這3個方向設(shè)置流速為0.1、0.2、…、0.9、1.0 m/s的來流，得到管匯受力如圖7。

圖7 管匯3個方向阻力變化曲線

代入數(shù)據(jù)求解管匯的阻力系數(shù)，x、y、z3個方向管匯迎流面積分別為21.90、27.35、50.60 m2。求出x、y、z3個方向管匯的阻力系數(shù)CDSX=1.06、CDSY=1.32、CDSZ=2.41。

可以看出，管匯x、y、z3個方向的阻力系數(shù)由于迎流面積的變大而變大，符合阻力系數(shù)的特點。另外，對比巴西石油公司在管匯模型試驗時所得的阻力值，選擇尺寸相近的管匯，和本次的模擬結(jié)果相比，兩者在同一數(shù)量級上。

3 附加質(zhì)量系數(shù)

3.1 驗證計算方法的準(zhǔn)確性

附加質(zhì)量產(chǎn)生的原因是物體在水中的變速運動帶動周圍的流體一起做不定常運動，這部分流體的運動消耗能量，造成了阻力，使物體的慣性增加。在求解附加質(zhì)量時，需要模擬管匯的非定常運動，有2種思路：一是管匯靜止，流體非定常運動繞流過管匯；二是流體域靜止，管匯做非定常運動[9]?？紤]到流體非定常運動模擬困難，同時Fluent軟件的動網(wǎng)格技術(shù)，便于管匯運動模擬，因此采用第2種方法來求解，附加質(zhì)量系數(shù)為

(8)

式中：FZ為非定常運動中某一時刻管匯的總阻力，N；FD為同一瞬態(tài)對應(yīng)速度下管匯定常運動的總阻力，N；mP為管匯入水后排出水的質(zhì)量，kg，大小為水的密度乘以管匯體積；a為所取瞬態(tài)的加速度，m/s2，通過速度公式求導(dǎo)可得。

DNV提供的附加質(zhì)量系數(shù)的理論數(shù)據(jù)如圖8，選擇和阻力系數(shù)計算中相同的長方體，b/a=2。采用本次研究中使用的計算附加質(zhì)量系數(shù)的方法，求得附加質(zhì)量系數(shù)為0.351。對比DNV提供的理論數(shù)據(jù)0.36，可見研究中采用的附加質(zhì)量系數(shù)的計算方法精確度良好，誤差較小，相對誤差約2.5%。

圖8 附加質(zhì)量系數(shù)的理論數(shù)據(jù)

3.2 附加質(zhì)量系數(shù)計算

動網(wǎng)格技術(shù)可以根據(jù)變化的區(qū)域邊界位置，通過對網(wǎng)格形狀、尺寸的變形來自動更新[10]。將流體域中管匯邊界指定為動網(wǎng)格，使其做周期性的非定常運動，運動形式選擇正弦運動(如圖9)，由用戶自定義函數(shù)UDF給定。其運動方程為

vel[a]=1+sin(6.28*time)

(9)

式中：a取值為0、1、2，代表x、y、z3個方向。

圖9 管匯速度-時間曲線

管匯非定常運動過程中，對特定時間點的數(shù)據(jù)進(jìn)行保存，設(shè)置求解器狀態(tài)為瞬態(tài)，運動總時間取4 s，步長為0.05 s，每20步即1 s進(jìn)行一次autosave，保存這一瞬態(tài)所有計算結(jié)果，計算完成后，讀取管匯非定常運動的總阻力值FZ。

FD的計算和阻力系數(shù)計算中總阻力的求解方法相同。通過模擬管匯的勻速繞流實現(xiàn)定常運動，速度和非定常運動中瞬態(tài)時刻相同，取1 m/s；其余流體域的設(shè)置條件相同。

mP為管匯入水后排出水的質(zhì)量，即海水密度乘以管匯體積，其值為124 167.728 kg。

a為瞬時加速度，通過速度公式求導(dǎo)，在所取時刻1、2、3、4 s時，瞬時加速度均為6.28 m/s2。

經(jīng)計算，管匯x、y、z3個方向的附加質(zhì)量系數(shù)如表1～3。

表1 x方向附加質(zhì)量系數(shù)

表2 y方向附加質(zhì)量系數(shù)

表3 z方向附加質(zhì)量系數(shù)

4 結(jié)語

1) 以水下管匯為研究對象，利用CFD軟件Fluent，計算其阻力系數(shù)和附加質(zhì)量系數(shù)。

2) 對DNV提供的模型進(jìn)行系數(shù)計算，對比其理論數(shù)據(jù)，驗證計算方法的精確性，然后用于管匯的計算。

3) 建立管匯和流體域計算模型，模擬穩(wěn)態(tài)流中的水池試驗，求出了管匯3個平動自由度的阻力系數(shù)和附加質(zhì)量系數(shù)。

4) 本文的計算方法效率高，成本低，精確度也能夠到保證，可為計算管匯的水動力系數(shù)提供參考。

[1] 姚麗琳.深水水下安裝下放運動分析技術(shù)研究[D].上海：上海交通大學(xué)，2011.

[2] 姚保太，康寧，鄭偉奇.變加速運動圓球附加質(zhì)量和阻力仿真分析[J].計算機輔助工程，2014，23(3)：82-87.

[3] 趙景浩.深水水下吊裝過程動力學(xué)建模與仿真[D].青島：中國海洋大學(xué)，2014.

[4] Antonio Carlos Fernandes.Pendulous installation method prospective model testing and numerical analysis[R].OMAE 2006-92655，Hamburg，Germany.

[5] 王兆立，牛江龍，秦再白，等.基于CFD理論的滑行艇阻力數(shù)值計算[C]// 第十四屆中國海洋(岸)工程學(xué)術(shù)討論會論文集(上冊)，2009.

[6] 湯珂，段夢蘭，張新虎，等.基于CFD的管匯非慣性水動力系數(shù)計算[J].石油礦場機械，2014，43(3)：34-37.

[7] DNV-RP-H103，Modelling and Analysis of Marine Operations[S].2011.

[8] 邵峰，金久才，張杰.無人船F(xiàn)LUENT阻力模擬與試驗驗證[J].海洋技術(shù)學(xué)報，2014，33(5)：8-12.

[9] 馮雙雙，沙永祥，魏曉娟，等.物體附加質(zhì)量CFD計算方法[C]//中國計算力學(xué)大會2014暨錢令希計算力學(xué)獎頒獎大會，2014.

[10] 黃昆侖，龐永杰，蘇玉民，等.潛器線性水動力系數(shù)計算方法研究[J].船舶力學(xué)，2008，12(5)：697-703.

CFDBasedDragCoefficientandAdditionalMassCoefficientCalculationofSubseaManifold

TIAN Yu，LUO Xiaolan，QI Yu，WANG Linan

(CollegeofMechanicalEngineering，ChinaUniversityofPetroleum，Beijing102249，China)

The installation of subsea manifold is one of the key technologies for the safety of underwater production systems.In this paper，the drag force and the additional mass force of hydrodynamics are analyzed.The CFD software Fluent is used to establish the model of manifold water and simulate the calculation of flow around the underwater operating environment.The three directions of drag force and its drag force coefficient are calculated.The use of dynamic grid technology to simulate the manifold’s accelerated movement in the flow field，to consider the manifold in the water in the variable speed drive around the fluid together to do with the constant movement of the formation of additional mass coefficient.By comparing with the empirical data of DNV，the error of the calculation result is small，which verifies the correctness of the calculation method and lays the foundation for the coupling analysis of ship-cable-body and the simulation control of the decentralized installation.

subsea manifold；CFD；drag force coefficient；additional mass coefficient

1001-3482(2017)06-0016-05

2017-06-29

國家工信部2013年高校船舶(海洋裝備)科研項目“水下采油樹配套工具研發(fā)”(工信部聯(lián)裝[2013]41號)；國家重點研發(fā)計劃重點專項項目(2016YFC0303700)

田 宇(1995-)，男，甘肅天水人，碩士研究生，主要從事工業(yè)機器人研究，E-mail：ty-cup123@163.com。

TE952

A

10.3969/j.issn.1001-3482.2017.06.004