矩形曲面網板水動力性能的數值模擬

劉圣聰+劉佳茗

摘 要：拖網網板為網衣提供水平張力，是拖網作業系統非常重要的漁具構件之一。該文利用CFD軟件FLUENT研究了展弦比對矩形曲面網板其水動力性能的影響。研究中設計了2種不同展弦比λ的網板，在流速1.54m/s，迎流沖角α = 0°～50°時，建立數值水槽進行數值模擬，得到網板的阻力系數Cd、升力系數Cl、俯仰力矩系數Cm和升阻比K，對比不同網板的水動力性能差異。結果表明，1號網板（λ = 2.0）與2號網板（λ = 1.5）的最大升力系數分別為0.85（α =15°）和0.92（α = 42°）；最大升阻比分別為9.59和8.35，俯仰力矩系數的絕對值分別為0.013和0.001。可見，展弦比越大，臨界沖角和Cl值越小，最大升阻比越高，但穩定性越差。研究結果可為拖網網板的結構優化設計提供參考。

關鍵詞：拖網網板；展弦比；水動力學性能；數值模擬

中圖分類號 S141 文獻標識碼 A 文章編號 1007-7731（2017）18-0092-04

Numerical Simulation on Hydrodynamic Performance of Rectangular Otter Boards with Different Aspect Ratios

Liu Shengcong et al.

（Dalian Tianzheng Industrial Co. Ltd，Dalian 116001，China）

Abstract：The otter board is a vital component for the efficient expansion performance of an otter trawl system. The numerical CFD code ANSYS fluent was conducted for rectangular otter board model to study the effect of aspect ratio on hydrodynamic performance. 2 different kinds of aspect ratio of otter boards were numerical simulated using the FLUENT with the flow velocity at 1.54m/s and the angle of attack measured from 0° to 50°. The cal-culated drag coefficient Cd，lift coefficient Cl，pitching moment coefficient Cm，lift to drag ratio K and the relations curve of these value and angle of attack were compared. The results showed that the maximum lift coefficient of otter boards 1 and 2 was 0.85（α =15°） and 0.92（α =42°），the maximum lift to drag ratio was 9.59 and 8.35，and the absolute value of pitching moment coefficient was 0.013 and 0.001 respectively. The result indicated that the greater the aspect ratio，the smaller the critical angle of attack and lift coefficient，the higher the maximum lift to drag ratio，but the worse of stability. The results would provide reference for the structural optimization design of trawl otter board.

Key words：Otter boards；Aspect ratio；Hydrodynamic performance；CFD analysis

拖網作業具有靈活、高效、適應性強等特點，是最常見的一種捕撈方式，也是目前我國遠洋捕撈中最主要的作業方式之一。網板是單船拖網系統中的重要組成部分，其水動力性能的優劣，關系到網具能否得到充分的擴張，從而直接影響到拖網漁具的漁獲性能和經濟效益[1-2]。因此，對網板水動力性能的研究以及對新型網板的開發，一直是國內外漁業工程研究者的關注重點。

挪威、日本、法國等漁業發達國家對網板水動力性能的改進研究較早，相繼開發出大展弦比矩形曲面網板、綜合型網板（橢圓形曲面開縫網板）、立式V型曲面網板和矩形V型曲面開縫網板等各種高性能網板[3]。A.Sala等[4]通過模型與原型試驗相對比的方法，研究了2種矩形V型曲面開縫網板的水動力性能及其對海底棲息地的破壞影響。Yuki等[5]用計算流體動力學方法研究了復翼式曲面網板的水動力特性，并通過軟件可視化將網板周圍的流場形態展現出來，與模型試驗進行對比，得到較好的吻合效果。Balash等[6]研究了2種柔性蝠翼式網板，分析了網板的膨脹和扭曲對網板水動力特性及穩定性的影響，發現膨脹可以提高升力，扭曲可以減小阻力從而提高升阻比，適當的扭曲對網板穩定性是有利的。

中國對拖網網板水動力性能的研究起步較晚，網板的結構形式較國外比較單一。中國的網板模型試驗主要在風洞設備中完成。郭根喜等[7]較系統的對網板進行了研究，包括網板的幾何結構、開縫數量、導流翼相對高度對不同類型網板水動力性能的影響。張勛等[8]、王錦浩等[9]分別研究了矩形V型曲面網板開縫口導流板曲率、展弦比、縫口位置和縫口寬度等參數對網板水動力性能的影響。王磊等[10]研究了葉板尺度比例變化對單縫立式曲面網板水動力性能的影響，發現導流板尺度過小會增加網板阻力，過長則會降低網板升阻比。劉健等[11]通過水槽模型試驗，研究了展弦比為1.0的矩形曲面網板在不同水平傾角和迎流沖角下的水動力性能。endprint

國內目前對網板水動力性能的研究手段多數是在風洞和水槽中進行模型試驗，對設備要求較高，同時需要大量的人力財力，且誤差較大。隨著計算機模擬仿真技術的發展，利用數值模擬的方法對網板水動力性能進行研究，將會更加省時省力，研究內容也會更加靈活。本文利用CFD軟件FLUENT對矩形曲面網板進行數值模擬，研究展弦比變化對網板水動力性能的影響，以期為優化拖網網板性能提供參考。

1 數值模型

1.1 控制方程 本文假設流體不可壓縮、各向同性，控制方程包括連續性方程和N-S方程：

[??V=0] （1）

[?V?t+（V??）?V=-1ρ?p+μρ?2V+F] （2）

式中：[V]為速度；[t]為時間；[p]為壓力；[F]為體積力；[ρ]為密度；[μ]為動力粘度。

1.2 參數定義 描述網板水動力特性的參數主要有雷諾數Re、阻力系數Cd、升力系數Cl、升阻比K、縱向力矩系數Cm。計算公式如下：

[Re=Vbρμ] （3）

[Cd=Fd（0.5?ρV2S）] （4）

[Cl=Fl（0.5?ρV2S）] （5）

[K=ClCd] （6）

[Cm=Mz（0.5?ρV2Sb）] （7）

式中：V為來流速度（m·s-1）；b為網板弦長（m）；ρ為流體密度（kgm-3）；μ為流體動力粘度（kgm-1s-1）；Fd為網板阻力（N）；S為網板面積（m2）；Fl為網板擴張力（N）；Mz為俯仰力矩（N·m）。

1.3 網板結構參數 網板的基本尺寸如圖1所示，L為網板翼展，b為網板翼弦，f為網板的最大厚度，展弦比λ=L/b，迎流沖角α為網板翼弦與來流方向的夾角，網板模型如圖2所示。因為矩形曲面網板為對稱結構，因此在計算時，可模擬網板的一半，將大大節省計算時間。

圖1 網板基本尺寸

圖2 網板模型

本文研究對象為2種不同展弦比的矩形曲面網板，迎流沖角α均為0°～50°范圍，網板的具體規格見表1。

表1 網板規格

[編號 翼展

（L/m） 翼弦

（b/m） 展弦比

（λ） 面積

（S/m2） 最大厚度

（f/m） 1 2.14 1.07 2.0 2.3 0.107 2 1.86 1.24 1.5 2.3 0.124 ]

1.4 模型參數設置 網板周圍的網格劃分如圖3所示，靠近網板的區域選擇六面體單元，遠場區域選擇四面體單元來進行計算，第一層網格高度為b/100。計算區域如圖4所示，水的密度ρ為998.2kg·m-3，動力粘度μ為0.001003kg·m-1·s-1，上游入口處設為速度入口邊界（velocity inlet），速度為1.54m/s，方向沿x軸正向的均勻流，湍流強度為2.6%，湍流長度尺度為0.07b。下游出口設為完全發展自由出流邊界（outflow）。計算區域的底面設為對稱邊界（symmetry），網板表面及計算域其他邊界均設為不可滑移壁面（no-slip wall）。

圖3 網板周圍網格劃分（α=20°）

圖4 計算域整體網格劃分

1.5 數值計算方法 經驗證，湍流模型選擇標準[k-ω]模型，控制方程采用有限體積法進行離散，對流項采用二階迎風格式，擴散項采用二階中心差分格式，計算步長為0.004s。

2 結果與分析

2.1 升力系數Cl和阻力系數Cd 網板的升力系數Cl和阻力系數Cd隨迎流沖角α的變化曲線如圖5所示。從圖5可以看出，網板的阻力系數Cd隨迎流沖角α的增大均呈上升趨勢。1號網板的升力系數Cl在α=0°～15°范圍內隨著α的增大而增大，在α=15°時達到最大值Clmax=0.85，此時的Cd=0.17。在α=15°～20°范圍時，Cl有所下降，在α=20°時Cl=0.73。在α=20°～35°之間，隨著α的增大，升力系數有略微增加，當α>35°后，隨著α的增大，升力系數呈下降趨勢。2號網板的Cl在α=0°～27°范圍內先增后降，在α= 15°左右達到較大值為0.80，在α=27°～42°范圍內，Cl持續上升，在α=42°時達到最大值Clmax=0.92，此時的Cd= 0.82。當α>42°后，升力系數開始下降。圖中Cl-exp和Cd-exp分別為2號網板的物理模型試驗所得到的升力系數和阻力系數曲線[7]，從圖中可以看出，數值模擬得到的Cd曲線與物模試驗相差不大，吻合較好。試驗得到的最大升力系數Clmax=1.08（α=27°），數值模擬為Clmax=0.92（α=42°），兩者在數值上相差不大，但是所對應的沖角有較大差距。產生誤差的原因有以下：（1）數值模擬沒有考慮摩擦力；（2）由于物模試驗時間較早，且在風洞中進行，試驗裝置和條件有限制，只在迎流沖角為5°～30°范圍內進行了試驗，數據有限。從阻力系數曲線的擬合效果與升力系數曲線的變化范圍和趨勢來看，數值模擬的結果是正確的。

圖5 升力系數Cl和阻力系數Cd隨沖角α的變化曲線

2.2 升阻比K和俯仰力矩系數Cm 如圖6所示，升阻比K隨著迎流沖角α的增加呈先升后降趨勢。從圖7可以看出，1號網板的升阻比K較高，最大為K=9.59（α=5°），2號網板的最大升阻比K=8.35（α=5°），兩種網板最大升阻比時對應的迎流沖角大致相同。但是，在最大升阻比對應的沖角α=5°時，網板的升力系數并不高，無法充分擴張網口，最大升力系數對應沖角附近升阻比的大小也很重要。1號網板在α=15°（最大升力系數對應沖角）時，升阻比K=5.05，2號網板在α=42°時，升阻比K=1.12。2號網板升阻比K的模擬結果與試驗相比較，兩者在α=5°時完全吻合，隨著α的增大，K均呈下降趨勢，數值模擬結果與試驗結果的變化趨勢是一致的。

圖6 升阻比K和俯仰力矩Cm隨迎流沖角α的變化曲線

俯仰力矩也稱為縱向力矩，是指作用在網板的流體動力對其重心所產生的力矩沿橫軸的分量，一般轉換為俯仰力矩系數Cm來分析比較。俯仰力矩系數有正負之分，其絕對值大小表示俯仰力矩的高低，一般Cm趨于0表示該網板的穩定性越優。在α=5°左右時，兩塊網板均達到最大升阻比K，此時1號網板Cm的絕對值為0.013，2號網板Cm的絕對值為0.001。在α=0°～20°范圍內，1號網板Cm的絕對值始終比2號網板大，意味著2號網板的穩定性較好。在α>25°后，兩塊網板Cm的絕對值均小于0.10，比較穩定。

3 結論

本文利用數值模擬的方法，研究了展弦比對矩形曲面網板水動力學性能的影響。模擬結果與物理模型試驗結果吻合較好，驗證了模擬的正確性。有以下結論：

（1）展弦比為2.0的1號網板的臨界沖角α=15°，Clmax=0.85，Cd=0.17。展弦比為1.5的2號網板的臨界沖角α=42°，Clmax=0.92，Cd=0.82。因此展弦比越大，臨界沖角越小，對應的Clmax和Cd越小。

（2）1、2號網板的最大升阻比K分別為9.59和8.35，在迎流沖角α=5°左右時達到，此時的俯仰力矩系數Cm的絕對值分別為0.013和0.001。在各自最大升力系數處，1、2號網板的升阻比K分別為5.05和1.12，俯仰力矩系數的絕對值分別為0.133和0.094。因此，展弦比越大，升阻比K越大，但穩定性越差。

參考文獻

[1]張勛，王明彥，徐寶生.拖網網板型式、結構與性能的研究與應用進展[J].中國水產科學，2004，11：107-113.

[2]徐寶生，張勛，王明彥.單船拖網網板的現狀及發展趨勢[J].福建水產，2010（1）：86-90.

[3]劉健，黃洪亮，吳越，等.2種立式曲面縫翼式網板水動力學性能的實驗研究[J].南方水產科學，2015，11（1）：68-74.

[4]Sala A，Farran JDP，Antonijuan J，et al. Performance and impact on the seabed of an existing- and an experimental-otterboard： Comparison between model testing and full-scale sea trials[J]. Fisheries Research，2009，100（2）：156-166.

[5]Takahashi Y，Fujimori Y，Hu F，et al. Design of trawl otter boards using computational fluid dynamics[J]. Fisheries Research，2015，161：400-407.

[6]Balash C，Sterling D，Lustica M，et al. Twist-and-camber effects on the performance of simple hydro-sails for efficiently spreading penaeid trawls[J]. Ocean Engineering，2015，109：161-168.

[7]郭根喜，劉同渝，黃小華.拖網網板動力學理論研究與實踐[M]. 廣州：廣東省出版集團，2008.

[8]張勛，王錦浩，王明彥，等.矩形V型曲面網板水動力性能的研究——Ⅰ導流板曲率與網板水動力性能的關系[J]. 中國水產科學，2004，11（s1）：5-8.

[9]王錦浩，張勛，王明彥，等.矩形V型曲面網板水動力性能的研究——Ⅱ展弦比、縫口位置、縫口寬度對網板水動力性能的影響[J]. 中國水產科學，2004，11（z1）：9-13.

[10]王磊，王魯民，馮春雷，等.葉板尺度比例變化對單縫曲面網板水動力性能的影響[J]. 漁業現代化，2015，42（6）：55-60.

[11]劉健，黃洪亮，陳帥，等.小展弦比立式曲面網板的水動力性能[J]. 水產學報，2013，37（11）：1742-1749.

（責編：張宏民）