立式曲面網板水動力及周圍流場數值仿真研究

徐琦　于文明　鄭建麗　梁晶晶

摘? 要：該文研究了立式曲面網板在不同沖角下的水動力性能變化情況，并且對網板周圍流場可視化及表面壓力分布可視化進行了探究。結果表明，當網板雷諾數在5.2×104～9.5×104區間時，升力系數、阻力系數和力矩系數隨雷諾數增大基本保持不變，而升力系數和力矩系數隨沖角增大呈先增大后減小趨勢;阻力系數隨沖角增大一直增大。升力系數的試驗值和模擬值隨沖角增大均呈先增大后減小趨勢，阻力系數的試驗值和模擬值隨沖角增大呈一直增大趨勢，升阻比的試驗值和模擬值隨沖角均呈先增大后減小趨勢。CFD結果顯示，網板表面壓力分布隨沖角增大逐漸增大，而網板壓力中心隨沖角增大基本保持不變。

關鍵詞：網板;水動力性能;模型試驗;數值模擬;流場可視化;壓力分布

中圖分類號：S971.4? ? ? ?文獻標識碼：A 文章編號：2096-4706（2020）09-0064-07

Numerical Simulation Research of Hydrodynamic Force and Surrounding

Flow Field of Vertical Curved Reticular Lamina

XU Qi1，YU Wenming2，ZHENG Jianli3，LIANG Jingjing1

（1.Zhejiang Ocean University，Zhoushan? 316022，China;2.South China Sea Fisheries Research Institute，Chinese Academy of Fishery Sciences，Guangzhou? 510300，China;3.Fishery Machinery and Instrument Research Institute，Chinese Academy of Fishery Sciences，Shanghai? 200092，China）

Abstract：This paper studies the changes in hydrodynamic performance of the vertical curved reticular lamina at different angles of attack，discusses the differences in the results of the different studies，and visualizes the flow field and surface pressure distribution around the reticular lamina. The results show that when the Reynolds number of the stencil is in the range of 5.2×104 to 9.5×104，the lift coefficient，drag coefficient and torque coefficient remain basically the same as the Reynolds number increases，while the lift coefficient and torque coefficient are first trend decreases after the increase;the drag coefficient increases with the angle of attack. The experimental value and the simulated value of the lift coefficient both increase first and then decrease as the angle of attack increases. The experimental and simulated values of drag coefficient have been increasing with increasing angle of attack. The experimental and simulated values of lift-to-drag ratio increase first and then decrease with the angle of attack. The CFD results show that the pressure distribution on the surface of the stencil gradually increases with increasing angle of attack，while the pressure center of the stencil basically remains unchanged as the angle of attack increases.

Keywords：reticular lamina;hydrodynamic performance;model test;numerical simulation;flow field visualization;pressure distribution

0? 引? 言

網板是單船拖網作業的重要屬具之一，主要作用為實現網口擴張、增大掃海面積，其水動力性能的優劣直接影響漁獲性能和經濟效益[1]。因此，開展網板水動力性能的研究尤為重要。筆者對網板在不同沖角下水動力性能的變化進行了研究，探討不同研究結果的差異，并實現網板周圍流場及表面壓力分布可視化。

在設計網板時一般須滿足幾項條件：一是在一定拖速內保證網口的擴張基本不變;二是要求網板具有較高擴張力和較低阻力;三是要求網板穩定性好;四是要求網板使用簡單，維修方便[2]。Park等[3]利用水槽試驗分析了不同彎曲度的立式曲面網板水動力性能，得到彎曲度為15%時升力系數最大，表層作業時其最適沖角范圍大于底層作業。現代拖網網板已有很多改進，包括采用新的結構形式和新型材料改進網板水動力性能，而日本漁業研究者對立式曲面網板進行了改進，采用不同的曲面形狀，使得較原型網板曲度減小，適宜沖角范圍增大，阻力降低，廣泛應用于大、中型單船拖網漁業中[4]。目前為止，在我國遠洋漁業中立式曲面網板仍有較高的使用率。因此，開展對立式曲面網板水動力性能的研究，對優化和開發新型網板具有十分重要的意義。

目前水槽模型試驗已成為研究網板水動力性能的標準方法，其具有試驗過程便于觀察、試驗條件可控的優點，研究者可以將實物網板按照尺度比例制作模型進行性能測試，而不需要制作出實物網板進行海上實測，這樣節省了大量的人力和財力。模型試驗雖是測量網板水動力性能標準的方法，但在試驗過程中消耗費用較多，并且需要精密的儀器[5]，而隨著計算機技術的發展，數值模擬技術被廣泛使用，計算流體動力學逐漸被漁業研究者應用于網板水動力性能研究中，其中FLUENT分析方法具有先進的求解技術，可以快速準確地求解出CFD結果。Xu等[6]對比了V形網板的水槽試驗和數值模擬結果，發現兩種研究結果吻合較好，最優型的網板結構為展弦比0.49，板面折角17°。Leifsson等[7]在使用CFD數值模擬技術研究立式曲面網板水動力性能時發現傳統的求解模型計算需要大量迭代次數，耗時太多，因此對模型算法進行了優化。利用模型試驗和數值模擬研究網板水動力性能已成為目前漁業學者的研究思路，通過對比兩種研究方法的結果即可互相驗證，也可彌補相互的缺點。

本文利用水槽模型試驗和CFD數值模擬方法研究立式曲面網板在不同沖角下水動力性能的變化，并對比兩種研究方法的結果差異性，提出科學的建議。同時，利用CFD數值模擬技術實現網板周圍流場可視化，直觀地分析網板水動力性能的變化，顯示出網板在水中所受到的壓力分布。

1? 材料與方法

1.1? 網板水動力特性計算公式

網板水動力特性的參數主要包括雷諾數Re、升力系數Cl、阻力系數Cd、垂直系數Cz、滾轉力矩系數Cmx、縱向力矩系數Cmz、升阻比K、翼展到壓力中心距離的壓力中心系數Cpb、翼弦到壓力中心距離的壓力中心字數Cpc，計算公式如下所示：

式中：c為網板翼弦（m），ρ為流體密度（kg/m3），b為網板翼展（m），V為流速（m/s），μ為流體動力粘性系數（kg/m/s），Fl為網板擴張力（N），Fd為網板阻力（N），Fz為網板垂直擴張力（N），Mz為縱向力矩（N·m），Mx為滾轉力矩（N·m），S為網板面積（m2），db為網板頂端到壓力中心距離（m），dc為網板翼端到壓力中心距離（m），α為網板沖角（°），β為網板傾角（°）。

水槽模型試驗中，當流速增大至某一特定值后，升力系數或阻力系數基本保持不變即進入自動模型區。本次試驗所使用升力系數、阻力系數和垂直系數均為進入自動模型區后測量值的平均值。

1.2? 水槽試驗條件設置

立式曲面網板模型如圖1所示，比例尺為1：20，展弦比為1.65，網板外側圓弧相對彎度38%，內側圓弧相對彎度23%。

水槽模型試驗在東海水產研究所循環水槽完成，試驗水槽尺寸長1.8 m，寬0.5 m，水深0.5 m，最大流速為2.5 m/s。三分力傳感器為日本共和電業制造的LSM-B-500NSA1-P型傳感器，量程500 N，精度0.1%，記錄的測量值為水流穩定后20 s內的平均值。流速計為螺旋式，由日本計測研株式會社生產。網板模型位于水槽中部，通過連接桿與三分力傳感器連接，如圖2所示。通過調節三分力傳感器下端轉軸可以改變網板模型的沖角。水槽模型試驗中流速設置為50 cm/s、60 cm/s、70 cm/s、80 cm/s、90 cm/s，網板沖角為0°～50°，其中在0°～15°和30°～50°范圍內時以5°為間隔，15°～30°之間設置為15°、18°、20°、22°、25°、27°和30°。水槽試驗中水溫為20 ℃，水密度為998.2 kg/m3，粘性密度為1.011×10-3 kg/（m*s）。

1.3? 數值模擬條件設置

數值模擬通過ANSYS 15.0軟件中的FLUENT完成[8]。數值模擬計算中使用有限體積法（Finite Volume Method，FVM）對Navier-Stokes方程進行離散。有限體積的求解方法采用SIMPLE算法，是一種壓力修正的算法，通過在交錯網格的基礎上計算壓力場，而后求解動量方程。FLUENT求解中離散格式設置為壓力二階格式、動量方程為二階迎風格式。湍流模型選擇可實現k-ε模型[9]（Realizable k-ε），其具有既滿足對雷諾應力的約束條件，又可以在雷諾應力上保持與真實情況湍流一致的優點。邊界層設置中使用增強壁面函數，加強了邊界層的網格密度，提高模型的計算精度。模擬計算區域如圖3所示，模擬計算區域與實際水槽規格一致，尺寸為1 800.00 mm*500.00 mm*500.00 mm。

網板沖角22°時數值模擬計算中局部區域和網板周圍網格劃分如圖4所示，網格劃分全部為四面體網格，最大網格尺寸為50.00 mm，最小網格尺寸為0.61 mm。網板邊界層周圍第一層網格厚度為5.1×10-2 mm（y+≈1），通過從網板表面向外膨脹15層，以1.2為增長率增強邊界層網格密度。

模擬計算區域總網格數為1.9×106，節點數為6.4×105。數值模擬計算區域邊界條件設置如圖5所示，沿X軸方向的水流速度為0.5 m/s，湍流強度為5%，湍流粘性率為10，速度入口為0.5 m/s，壓力出口為0 Pa，所有網板外壁面和所有水槽內壁面均設置為固定無滑移條件。

2? 結果與分析

2.1? 雷諾數與網板水動力系數的關系

水槽模型試驗中將測量值代入網板水動力特性計算公式（1）—（6）得到網板升力系數、阻力系數、力矩系數與雷諾數和沖角的相互關系如圖6所示。由圖可知，不同沖角條件下升力系數Cl隨雷諾數增大基本保持不變，在一定雷諾數條件下，升力系數Cl隨沖角增大呈先增大后減小趨勢，在沖角25°時達到最大。阻力系數Cd隨沖角增大呈一直增大趨勢，在沖角不變條件下隨雷諾數增大基本保持不變。在沖角0°～25°之間，力矩系數Cmx隨沖角增大逐漸增大，在25°～50°之間，隨沖角增大逐漸減小。力矩系數Cmz隨沖角增大呈先增大后減小的趨勢，在沖角30°時達到最大。升力系數、阻力系數和力矩系數在一定沖角下隨雷諾數增大基本保持不變，因此可認為當雷諾數為5.2×104～9.5×104時，在設置的流速和沖角條件下進入自動模型區。

2.2? 模型試驗和數值模擬結果對比分析

立式曲面網板的升力系數、阻力系數和升阻比與沖角相互關系的試驗值和模擬值對比結果如圖7所示。

由圖可知，試驗值和模擬值的升力系數Cl均隨沖角增大呈先增大后減小趨勢，均在臨界沖角25°時達到最大值，分別為1.45和1.59，當沖角大于臨界沖角后升力系數值急劇下降，兩種方法研究結果在沖角小于25°時差異較大，在沖角大于25°后差異較小，升力系數的模擬值始終大于試驗值。阻力系數Cd的試驗值和模擬值均隨沖角增大逐漸增大，在沖角18°～22°之間兩種方法的結果偏差較大，在沖角小于15°時吻合較好，試驗值在沖角大于27°后有明顯增大趨勢，阻力系數的模擬值也始終大于試驗值。升阻比K的試驗值和模擬值隨沖角均呈先增大后減小的趨勢，其值均在沖角15°時達到最大，分別為3.25和2.93，其中試驗值在沖角25°時又明顯增大，可能原因為試驗中出現的操作誤差。升阻比的試驗值和模擬值在沖角大于27°后吻合較好，在沖角10°～25°之間偏差較大，適宜的工作沖角范圍均為10°～22°之間，其中試驗值始終大于模擬值?？傊⑹角婢W板升力系數、阻力系數和升阻比的試驗值和模擬值偏差較小，證明了數值模擬的準確性。

2.3? 網板周圍流場可視化及表面壓力分布

立式曲面網板在不同沖角下網板中心平面周圍流場分布如圖8所示。圖中水流方向為從左至右，初始流速為0.5 m/s。由圖8可知，網板外側箭頭指向部分流速均大于0.6 m/s，內側圓形標示部分流速均小于0.5 m/s，由于外側流速增大，壓力減小，內側流速減小，壓力增大，導致在網板內外兩側形成壓力差，從而產生網板升力，圖中箭頭指向部分隨沖角增大呈先增大后減小趨勢，在沖角25°時達到最大，這與模型試驗結果相符合。網板后部矩形標示部分流速為負值，表示在該處形成渦旋，產生渦旋的主要原因為水流經過網板背部時在網板表面形成一層很薄的邊界層，隨著流速的增大水流開始逐漸離開網板表面，造成邊界層分離情況，由于邊界層分離時粘性力的存在，分離后會在后部區域產生渦旋，形成壓差阻力，隨著沖角增大該區域面積增大，表明阻力值逐漸增大，與模型試驗結果吻合。隨著沖角增大網板左側流速衰減區逐漸減小，主要原因是隨著沖角增大，水流與網板接觸面積逐漸增大，降低了流速減弱的程度。

不同沖角下網板翼端平面周圍流場分布如圖9所示。由圖9可知，網板背部產生的邊界層分離點隨著沖角增大逐漸向網板前端移動。網板翼端尾部矩形標示部分流速為負值，是水流與網板表面分離時形成的翼端渦，其隨著沖角增大逐漸增大。網板前端流速衰減區隨著沖角增大逐漸增大，后端流速衰減區隨著沖角增大明顯變寬。

在不同沖角條件下網板表面壓力分布如圖10所示，隨著沖角增大網板表面壓力分布逐漸增大，主要原因為沖角增大后網板與水流的接觸面積增大，導致網板表面壓力增大區域也逐漸增大，而網板壓力中心的壓力隨沖角增大基本保持不變。

3? 討論

3.1? 網板周圍流場可視化

本文利用水槽模型試驗和數值模擬研究了立式曲面網板水動力性能的變化，并對比分析了兩種方法得到的研究結果，最后利用數值模擬實現了網板周圍流場可視化，為今后網板的性能的研究提供了科學依據。關于網板水動力性能的研究國內外已屢見不鮮，而利用可視化手段研究網板周圍流態分布還較少，目前主要實現網板周圍流場可視化的方法主要有三種：線條法[10]、氣泡法[11]和數值模擬，其中線條法是指將輕柔線條一端固定在網板表面，根據線條另一端在水流中所呈現出的姿態來判斷網板周圍的流態;氣泡法指將水槽中的水通電，電解產生氫氣泡，氣泡隨著水流經過網板后流態發生改變進而判斷流態的分布情況。但線條法和氣泡法只能對網板周圍流場進行定性的研究，不能精確地顯示網板周圍流態的速度場。數值模擬技術具有高效、低成本等優點，可將網板周圍流場精確地描述出來，因此利用數值模擬技術研究網板周圍流態分布與網板水動力性能的關系，即可以掌握網板周圍流態分布變化的原因，也可為后續網板的優化與設計提供直觀的科學依據。

3.2? 模型試驗和數值模擬結果比較

網板的升力系數和阻力系數數值模擬結果均較水槽模型試驗結果高，其主要可能由兩個方面的原因導致：一是數值模擬在構建網板模型時未考慮摩擦力，使得數值模擬測試值偏高;二是因為在水槽模型試驗中水流的不穩定性造成網板的震蕩，使得測量水動力性能的儀器產生誤差。網板模型試驗時尺度比的大小會直接影響試驗結果的準確性，因此在今后網板模型試驗中應選擇合適的模型尺度比，盡量降低模型的尺度效應對試驗結果的影響。

3.3? 網板開發研究趨勢

網板在拖曳過程中多數情況下會存在傾斜的狀態，網板在內傾狀態下有利于增加穩定性并能夠順利越過障礙物，但不能內傾過大，會破壞平衡狀態，造成起浮。同時網板正常工作狀態也存在后傾，這樣安裝有利于越過障礙物。網板曳行時也會不同程度地受到外力的干擾，也可能會導致網板出現傾斜的情況，而不同的傾斜程度對水動力性能有一定影響，因此今后可開展不同傾斜程度對網板水動力性能的影響研究。影響網板水動力性能的因素除網板自身結構參數外，還有作業方式和作業海區底質類型等因素，在實際生產作業中，海底對網板的水動力性能有較大影響。本文利用模型試驗和數值模擬研究立式曲面網板水動力性能的過程中，由于盡量避免無關因素影響網板性能測試數據的準確性，忽略了海底底質對網板水動力性能的影響，而且模型試驗中模擬海底底質也是非常困難的，但隨著計算機技術的發展，在數值模擬中可利用離散單元法建立模擬海底，通過設置海底底質黏土的材料和粘度指數來實現不同海底底質的類型，并通過連接力模型將海底與網格化的網板整合進行求解計算[12]。因此，今后研究網板水動力性能中可引入分析不同海底底質對網板水動力性能的影響。

4? 結? 論

立式曲面網板的水動力性能參數在進入自動模型區即雷諾數為5.2×104～9.5×104區間內時變化趨勢為：升力系數的試驗值和模擬值隨沖角增大均呈先增大后減小趨勢。阻力系數的試驗值和模擬值隨沖角增大呈一直增大趨勢。升阻比的試驗值和模擬值隨沖角均呈先增大后減小趨勢。

由立式曲面網板周圍流場分布可知，網板升力隨沖角增大呈先增大后減小趨勢，在沖角25°時最大。網板壓差阻力隨沖角增大呈一直增大趨勢。網板背部產生的邊界層分離點隨沖角增大逐漸向網板前端移動。網板翼端尾部形成明顯的翼端渦，其隨沖角增大逐漸增大。網板前端流速衰減區隨沖角增大逐漸增大。網板表面壓力分布隨沖角增大逐漸增大，而網板壓力中心隨沖角增大基本保持不變。

參考文獻：

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

[2] 孫滿昌.海洋漁業技術學 [M].北京：中國農業出版社，2005.

[3] PARK C D，MATUDA K，HU F，et al. Hydrodynamic Characteristics of Cambered Plates in Free Stream and near the Bottom [J].Nippon SuisanGakkaishi，1993，59（4）：627-632.

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

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

[6] XU Q C，HUANG L Y，ZHAO F F，et al. Study on the hydrodynamic characteristics of the rectangular V-type otter board using computational fluid dynamics [J].Fisheries Science，2017，83（2）：181-190.

[7] LEIFSSON L，HERMANNSSON E，KOZIEL S. Optimal shape design of multi-element trawl-doors using local surrogate models [J].Journal of computational science，2015，10：55-62.

[8] 凌桂龍.ANSYS Workbench 15.0從入門到精通 [M].北京：清華大學出版社，2014.

[9] SHIH T H，LIOU W W，SHABBIR A，et al. A new k -? eddy viscosity model for high reynolds number turbulent flows [J].Computers Fluids，1995，24（3）：227-238.

[10] PARK C D，MATUDA K，TOKAI T. Surface Flow Visualization of Flat Plates by Tuft Method [J].Nippon Suisan Gakkaishi，1994，60（2）：193-199.

[11] PARK C D，MATUDA K，TOKAI T. Flow Visualization around Cambered Plates Using Hydrogen Bubbles [J].Nippon Suisan Gakkaishi，1994，60（4）：485-491.

[12] JUN C W，SOHN J H，PARK S H. Effect analysis of design parameter of trawl door by considering interaction with seabed [J].International Journal of Precision Engineering & Manufacturing，2016，17（11）：1515-1522.

作者簡介：徐琦（1988—），男，漢族，吉林榆樹人，漁業一級船長，中級職稱，工學學士，研究方向：漁具漁法。