水深對潮流能電站載體水動力系數的影響

陳鷥鷺，謝永和，許頌捷

(1.江蘇科技大學船舶與海洋工程學院，江蘇鎮江 212000;2.浙江海洋學院船舶與海洋工程學院，浙江舟山 316000）

水深對潮流能電站載體水動力系數的影響

陳鷥鷺1，謝永和2，許頌捷2

(1.江蘇科技大學船舶與海洋工程學院，江蘇鎮江 212000;2.浙江海洋學院船舶與海洋工程學院，浙江舟山 316000）

利用三維線性勢流理論和源匯分布法研究潮流能電站載體的輻射問題，對不同水深下載體的附加質量和附加阻尼進行數值計算。計算結果表明，3個水平面運動(橫蕩、縱蕩和首搖）附加質量隨水深增大變化較為平緩，而3個垂向運動(垂蕩、橫搖和縱搖）附加質量隨水深變化顯著。載體的6個運動附加阻尼都隨水深增加而增加，且波浪周期在4～9 s時增幅較大。

潮流能電站;水動力系數;附加質量;附加阻尼

0 引言

潮流能是一種清潔無污染、蘊藏量豐富的可再生新能源，隨著化石能源不斷枯竭和環保要求的不斷提高，世界各國對潮流能的開發利用愈加重視。我國潮流能資源極為豐富，據《中國沿海農村海洋能資源區劃》統計，中國沿岸130個水道的理論平均功率約為1.4×104MW，其中浙江、福建、山東和臺灣海峽是世界上潮流能資源最豐富的地區之一。

漂浮式潮流能電站是一種新型的潮流能利用轉換裝置，其載體承受水輪機的工作載荷以及整個機組的重量，載體的水動力性能是確定潮流能電站在波浪中運動性能優良的重要指標。本文利用三維線性勢流理論，在頻域內研究潮流能電站載體的輻射問題，應用三維源匯分布法對不同水深下潮流能電站載體運動的水動力系數，包括附加質量和附加阻尼進行數值計算與分析。

1 理論基礎

假定結構物周圍的流場是無旋、無粘、不可壓縮的理想流體，流動的基本方程為關于速度勢的線性Laplace方程［1］。在微幅運動下，假定物體是在平衡位置附近作簡諧搖蕩運動，圓頻率為ω，則采用分離變量法，并將定常空間速度勢按式

進行分解。φj為載體單位振幅運動的輻射勢。定解條件如下:

采用源匯分步法，將速度勢表示成物面上的分布源形式［2］:

面元法即Hess-Smith方法［4］是求解分布源積分方程的一種常用數值方法。它是通過對船體平均濕表面的離散，將上述源積分方程轉換成線性代數方程組來求解不同運動模態的源強分布，從而根據式(2）確定場內各點的輻射速度勢。然后按照定義計算出載體的水動力系數:

2 載體外形設計

考慮到雙體船具有耐波性好，拖航阻力小等特點，本文選用雙體船作為潮流能電站的載體［5］。設計時由水輪機固定工作直徑的大小確定片體間的距離及船長，為避免載體有產生的縱傾從而引起水輪機葉片上緣出水，及水輪機工作時產生較大的側向力導致載體船首或船尾出水［6］，設計載體的最終主尺度參數見表1。

表1 潮流能電站載體的主要參數Tab.1 Main parameters of the carrier

3 計算模型

本次數值模擬采用SESAM軟件的建模模塊GeniE及水動力計算模塊HydroD，在HydroD中調用WADAM程序計算，POSTRESP進行結果分析。

首先在GeniE中建模，由于載體可視為大尺寸構件，采用板模型進行建模，將模型外表面設為濕表面，創建工況加載水壓，設置網格密度為0.5，生成網格，導出有限元模型文件 (T*.FEM）如圖1所示。

圖1 載體有限元模型Fig.1 Finite elementmodel

在一般情況下，常把海浪視為無限多個頻率不等、方向不同、振幅變化而相位雜亂的微幅簡諧波疊加而成的不規則波系，這些特征可以用海浪譜表示出來［7］。本文波浪載荷的直接計算采用Bretschneider譜，波浪周期為4～38 s，步長為2 s。考慮水深對潮流能電站載體水動力系數的影響，分別選取水深為10 m，20 m，30 m及100 m進行計算。在HydroD中定義上述環境參數，將有限元模型文件 (T*.FEM）導入 (見圖2），即得到潮流能電站載體水動力模型。

圖2 載體水動力模型Fig.2 Hydrodynamic model

4 數值計算結果

附加質量系數和附加阻尼是船體形狀、振動頻率和前進速度的函數，附加質量系數和附加阻尼均由一個6×6的矩陣表示，本文主要分析與潮流能電站載體水動力性能有關系的6個對角線上的系數。附加質量以μjk表示，μ11，μ22和μ33分別為縱蕩、橫蕩、和垂蕩3個運動模態的附加質量，μ44，μ55和μ66為橫搖、縱搖和首搖3個角運動的附加慣性矩，同上附加阻尼以λjk表示。計算完成后從POSTRESP中顯示計算結果如圖3～圖5所示。

圖3 載體的附加質量Fig.3 Added mass of the carrier

6個運動模態附加質量如圖3所示。波浪周期在4～9 s時，水深越大附加質量越大;波浪周期在9～38 s時，水深越大附加質量越小。在同樣深度下，橫搖附加質量最大，而橫蕩時附加質量最小。這與載體雙體船船型特征相符合［8］。

從圖3中可以看出，水平面運動附加質量隨水深變化規律并不明顯，而垂向運動附加質量隨水深變化顯著。若以μ/μ∞代表淺水和深水運動下附加質量的比值［9］，其隨水深的變化曲線如圖4所示，水深由10 m變化到20 m時，垂蕩、橫搖、縱搖附加質量的變化，分別減少了19.4%，13.5%，10.6%;縱蕩、橫蕩、首搖的附加質量變化相對平緩，分別減少了2%，8%，0.1%，水深由20 m增加到30 m時，垂蕩、橫搖、縱搖附加質量分別減少了7.1%，2.6%，1.7%，而水平面運動的附加質量變化很小。

圖4 μ/μ∞隨水深變化規律Fig.4 μ/μ∞ variation with depth

圖5 載體的附加阻尼Fig.5 Added damping of the carrier

從圖5可知，所有6個運動模態附加阻尼隨水深變化規律大致相同且在波浪周期在4～9 s時變化較為明顯，附加阻尼隨著水深增大而增大。所有自由度的附加阻尼中從數值上比較，載體橫搖、縱搖附加阻尼最大，縱蕩、橫蕩附加阻尼最小。橫搖與橫蕩的附加阻尼相差3個數量級。

5 結語

本文通過對潮流能電站載體在深水及幾種淺水水深下的數值計算，得到附加質量和附加阻尼等水動力系數的變化規律，得出如下結論:

1）在所有水深下，載體的3個水平面運動附加質量變化較為平緩，而3個垂向運動附加質量變化顯著。

2）相對于深水條件下，所有6個運動附加質量都在水深小于20m時變化劇烈，水深大于20m時則趨于平緩。

3）數值上，載體6個運動附加阻尼都隨水深增加而增加，且波浪周期在4～9 s時增幅較大。

4）總體上，由于載體在淺水中附加質量和附加阻尼的較大增加，垂蕩、橫搖、縱搖運動幅度將相對趨于緩和，從而給潮流能電站載體提供更安全穩定的航行環境，避免水輪機出水或觸底現象的發生。

［1］LIDong-li，YANG Liang，ZHANG Hong-yu，et al.Prediction of ship maneuverability based on CFD method［J］.Journal of Wuhan University of Technology，2009，31(24）:120 －123.

［2］倪崇本，朱仁傳，繆國平，等.一種基于CFD的船舶總阻力預報方法［J］.水動力學研究與進展，A輯，2010，25(5）:580－586.

NIChong-ben，ZHU Ren-chuan，MIAO Guo-ping，et al.A method for ship resistance prediction based on CFD computation［J］.Chinese Journal of Hydrodynamics，2010，25(5）:580－586.

［3］HUANG Hao，GUO Hai-qiang，ZHU Ren-chuan，et al.Computations of ship roll damping in viscous flow［J］.Journal of Ship Mechanics，2008，12(4）:568 －573.

［4］詹成勝，劉祖源，程細德.潛艇水動力系數數值計算［J］.船海工程，2008，6(3）:1 －4.

ZHAN Cheng-sheng，LIU Zu-yuan，CHENG Xi-de.Numerical calculation of the submarine's hydrodynamic coefficients［J］.Ship and Ocean Engineering，2008，6(3）:1 －4.

［5］荊豐梅，盛其虎，張亮，等.漂浮式潮流電站總體設計［C］.中國可再生能源學會海洋能專業委員會第三屆學術討論會論文集，溫嶺，2010:116－124.

［6］龐永杰，楊路春，李宏偉，等.潛體水動力導數的CFD計算方法研究［J］.哈爾濱工程大學學報，2009，8(4）:903 －908.

PANG Yong-jie，YANG Lu-chun，LI Hong-wei，et al.Approaches for predicting hydrodynamic characteristics of submarine objects，Journal of Harbin Engineering University［J］.2009，8(4）:903 －908.

［7］劉寅立，焦永芳.波浪能開發與利用研究進展［J］.中國高新技術產業，2009(2）:19－20.

LIU Ying-li，JIAO Yong-fang.Progress of research on wave energy［J］.Chinese Hi-tech Enterprises，2009(2）:19 －20.

［9］黃祥鹿，陳小紅，范菊.錨泊浮式結構波浪上運動的頻域算法［J］.上海交通大學學報，2001，45(10）:1470 －1476.

HUANG Xiang-lu，CHEN Xiao-hong，FAN Ju.Calculation of themotion of amoored floating structure by using second order frequency domain method［J］.Journal of Shanghai Jiaotong University，2001，45(10）:1470 －1476.

［10］張楠.漂浮式潮流電站葉輪與錨泊系統設計研究［D］.哈爾濱:哈爾濱工程大學，2010:38－64.

ZHANG Nan.Design and research of the turbine and mooring system of floating station formarine current energy generation［D］.Harbin:Harbin Engineering University，2010:38－64.

The effects of water depth on the hydrodynam ic coefficients of tidal current energy power p lant carrier

CHEN Si-lu1，XIE Yong-he2，XV Song-jie2
(1.School of Naval Architecture and Ocean Engineering，Jiangsu University of Science and Technology，Zhenjiang 212000，China;2.School of Naval Architecture and Ocean Engineering，Zhejiang Ocean University，Zhoushan 316000，China）

By using the 3-D linear potential theory，the radiation problem of tidal current energy power plant carrier is investigated and the hydrodynamic coefficients in various water depths are calculated numerically.It is shown that the added masses of heave，roll and pitch motions changes slowly with an increase in depth，while the added masses of surge，sway and yaw motions varies significantly.The added damping of all the motions of the carrier increased remarkably with an increase in depth，with a maximal increase slope occurs in the wave period is between 4 to 9 seconds.

tidal current energy power plant;hydrodynamic coefficient;added mass;added damping

U661.3

A

1672－7649(2014）04－0017－05

10.3404/j.issn.1672－7649.2014.04.003

2013－05－06;

2013－06－13

國家自然科學基金資助項目(51279182）

陳鷥鷺(1989－），女，碩士研究生，研究方向為海洋結構物水動力性能。