1. 引言
单船拖网是渔业捕捞中主要作业方式之一,在我国发展历史悠久,也是目前我国远洋渔业中主要的作业方式之一 [1] [2] [3]。网板是单船拖网作业系统中的重要构件,它具有实现网具的水平扩张、兼具网具沉降器等作用,其阻力占整个拖网系统阻力的13%~27% [4],甚至达到30% [5]。网板类型的选择、其水动力性能的优劣以及对海洋环境的适应性,是影响拖网的渔获效率和经济效益的重要因素 [6] [7] [8],因此开展网板水动力性能的研究对渔业节能减排、渔业绿色发展具有重要意义。一些发达渔业国家,如日本、丹麦、挪威等国家都对网板的水动力性能进行了深入研究,相继开发出一系列高性能网板 [9]。20世纪30年代,德国学者设计出一种大展弦比曲面网板(Suberkrub网板),具有升力系数大、效率较高等优点;20世纪50年代,前苏联研制出一种椭圆形平面开缝网板,该网板对渔场的适应性较强,适宜在高低不平的海底拖曳;20世纪70年代,法国研制出一种椭圆形曲面开缝网板,兼具扩张效率高和能越过劣质海底等特点;20世纪80年代,日本研发出一种立式V型曲面网板,它吸取了V型网板稳定性好和大展弦比曲面网板扩张力大等优点,并加以改进;20世纪90年代,丹麦基于V型网板和曲面网板成功开发了矩形V型曲面网板,有效网板的效率;20世纪末,瑞士吸取流体力学的最新研究成果研制成一种高性能喷射型网板,是当今世界上最先进的底拖网板之一;此外还有复翼型网板等 [10] [11] [12]。由于我国对网板的研究起步较晚,因此国产网板作业性能与一些先进的渔业国家还存在不小差距,目前仍以借鉴国外先进的高性能网板为主 [9] [10]。
目前,网板水动力性能的研究方法主要有理论分析、模型实验(风洞试验和水槽试验)、数值模拟与海上实测等几种方式。Kawakami等建立网板受力分析理论模型并通过模型试验进行验证 [13];Matuda等通过水槽模型试验获得传统V型网板最大升力系数为1.27 [14];王磊等通过风洞模型试验研究发现叶板尺度不仅影响网板最大升力系数,同时影响网板作业稳定性 [15];Takahashi等通过水槽模型试验验证数值模拟在网板水动力学研究的可行性 [16],Xu等利用数值模拟对V型网板展弦比、后退角等参数进行优化 [17];Sala等通过海上试验对比研究新型Clarck-Y网板与传统V型网板水动力性能,并通过海底反作用力分析网板对海底的影响 [18]。
本文研究的圆形曲面网板借鉴上述一些高效网板的优点为我国远洋渔业大型拖网船设计的网板,采用了水槽模型试验和数值模拟方法研究了该网板的水动力学性能,获得其最大升力系数、最大升阻比、临界冲角等水动力参数,并对两种方法的试验结果进行对比分析。此外,根据数值模拟得出的网板周围流场的分布,分析网板周围流场与网板水动力性能之间的关系。研究结果可为圆形曲面双导流板网板提供水动力参数,也可为网板水动力学研究方法的选择提供参考。
2. 材料与方法
2.1. 模型网板
网板原型为全钢结构,作业时拖速为3~4节,其主要结构参数为,翼展l = 0.55 m,翼弦c = 0.55 m,网板面积S = 0.24 m2,展弦比AR = l2/s = 1.25。
图1所示为模型网板结构图,D1和D2为两块导流板,P为主面板,导流板之间及导流板与主面板之间的间距为10 cm。根据弗洛德相似准则 [1] [19] [20] 设计制作网板模型,其关系式如下所示:
(1)
其中,下标“1”代表原型网板,下标“2”代表模型网板,v代表流速,s代表尺度比,本文中尺度比为s = 3.0。模型网板为全钢结构,钢板厚度为5 mm,其具体结构参数如表1所示。
Figure 1. Structure of the model otter board (The coordinate frame on it was used to calculate the center-of-pressure of the otter board)
图1. 模型网板结构图(图中坐标系用于计算网板压力中心系数)
Table 1. Detailed parameters of the model otter board
表1. 模型网板参数表
2.2. 试验方法
2.2.1. 水槽试验
模型网板的水槽试验在中国海洋大学增殖养殖工程水动力实验室循环水槽进行,水槽试验段尺寸为4 m (长) × 1.2 m (宽) × 1 m (高)。水槽配有六分力仪传感器(量程0~50 kg,精度±0.3%),以及小威龙Vectrino点式流速仪(量程0~2 m/s,精度±0.5%),分别测量网板升阻力及力矩以及流速等数据。试验过程中流速设置为0.2~0.6 m/s (间隔0.05 m/s),即雷诺数为Re = 1.09~3.28 × 105,网板冲角(AOA)为0˚~70˚,间隔5˚,分别采集各个工况下力及力矩的数据,数据采集频率为100 Hz,采集时间为90 s。
2.2.2. 数值模拟
数值模拟利用ANSYS 15.0软件进行,对水槽模型试验进行数值模拟,模拟试验条件及模拟工况与模型试验相一致。将水槽的试验段尺寸作为数值模拟的计算域 [17] [21],即4.0 m ×
1.2 m
× 1.0 m为水槽模拟计算域。湍流模型选取k-ε EARSM湍流模型,并利用可伸缩壁面函数(Scalablewallfunction)进行近壁区处理 [22] [23],数值模拟的具体参数设置如表2所示。
Table 2. Parameter settings of the numerical simulation
表2. 数值模拟参数设置
数值计算过程中,网格划分为四面体非结构化网格,对网板周围网格进行局部加密,并利用偏斜度(Skewness)检查体网格质量,其最大值小于0.83时认为网格质量较好,可以用于数值计算。对于水槽数值模拟(图2),网格数约为2.66 × 106,计算节点数约为4.81 × 105,y+值约为11.06~80.25,湍流强度为3.27%~3.75%,湍流粘度为3.54~9.25 × 10-4m2/s。
Figure 2. Computational grids of the domain during flume tank simulation
图2. 水槽试验网格划分图
2.3. 数据处理
根据所测得网板升、阻力及力矩数据,分别求得网板升力系数、阻力系数、升阻比以及力矩系数,计算公式如下:
(2)
(3)
(4)
(5)
(6)
其中,FL、FD分别为网板升力和阻力;CL、CD分别为网板升、阻力系数;ρ为流体密度;S为网板面积;Mx、Mz分别为网板沿x轴和z轴的力矩(见图1);K为网板升阻比;Cmx、Cmz为力矩系数。根据所求力矩系数及网板升阻力系数值,可计算出网板的压力中心系数,计算公式如下:
(7)
(8)
其中,Cpl、Cpc分别为网板在翼展(图1,z轴)和翼弦(图1,x轴)方向的压力中心系数;dl、dc分别为压力中心点距网板上缘和前缘的距离。
3. 结果
3.1. 网板升力系数
网板升力系数(CL)随冲角变化的关系曲线如图3所示,升力系数随冲角的增加先增大后减小,模拟结果与试验结果比较接近,且与试验结果具有相似的变化趋势。水槽模型试验结果显示,该网板最大升力系数为CL-MAX = 2.421,临界冲角为50˚;水槽数值模拟CL-MAX为2.309,临界冲角为50˚~55˚,对最大升力系数的模拟误差为4.63%。
3.2. 网板阻力系数
网板阻力系数(CD)随冲角的变化如图3所示,模拟结果与试验结果比较接近,且具有相近的变化趋势。网板阻力系数随冲角的增加逐渐增大,当冲角超过临界冲角时,水槽试验阻力系数增加趋势减缓。在0˚~70˚冲角范围内,水槽模型试验测得网板阻力系数为0.112~2.301,水槽数值模拟阻力系数为0.096~1.976。
3.3. 网板升阻比
升阻比(K)是判断网板水动力效率的重要参数,其随冲角的变化关系如图4所示。网板升阻比随冲角的增加先增大后减小,且模拟值与实测值比较接近。水槽试验测得网板最大升阻比KMAX为3.719,此时冲角为AOA = 25˚,水槽模拟KMAX为3.942,冲角为AOA = 20˚,对最大升力系数的模拟误差为6.00%。
Figure 3. The result of flume tank model experiment and numerical simulation
图3. 水槽模型试验及其数值模拟结果
Figure 4. Lift-to-drag ratio of the otter board in relation to the angle of attack
图4. 网板升阻比随冲角的变化
3.4. 网板压力中心系数
图5所示为网板压力中心系数(翼展方向压力中心系数Cpl和翼弦方向压力中心系数Cpc)随冲角的变化关系曲线,由于模型试验(水槽)时仅测得z向力矩Mz,因此试验值仅测得网板翼弦方向压力中心系数Cpc的值。水槽试验显示,模拟Cpc值与实测值比较接近,且随冲角的增加先快速减小后趋于稳定,而Cpl值则稳定在0.5附近,即随冲角的增加,网板压力中心点先往网板前缘移动,冲角超过25˚之后,压力中心点趋于稳定。
3.5. 网板周围流场
基于上述研究结果,数值模拟可以较好的反应网板水动力学性能,因此本文选取网板中心横截面分析网板周围流场状况,结果如图6所示。图中网板冲角分别设置为45˚、50˚、55˚和60˚,即该网板临界冲角附近。水槽试验流线图显示当冲角小于55˚时,网板迎流面、背流面流线均比较平缓,没有流分离现象,同时在该冲角范围内网板背流面流速逐渐增大,即网板迎流面、背流面压差逐渐增大,因此网板升、阻力系数均随冲角的增加而增大;冲角达到60˚时,网板前缘出现流分离现象并且在网板背流面形成涡流,此时网板背流面流速开始减小,网板迎流面、背流面压差逐减小,导致网板压差所产生的升、阻力均降低。
Figure 5. Center-of-pressure coefficients of the otter board in relation to the angle of attack
图5. 网板压力中心系数随冲角的变化
Figure 6. Streamlines from a horizontal plane cutting through the center section of the otter board
图6. 网板平面中心横切面流线图
4. 讨论与结论
4.1. 网板水动力性能
本文利用水槽试验及数值模拟方法研究圆形曲面网板水动力性能,结果显示,该网板最大升力系数CL-MAX高于2.0,最大升阻比KMAX超过3.5,与传统矩形网板 [24]、V型网板 [1]、V型曲面网板 [25] 相比,该网板具有良好的水动力性能。同时,该网板临界冲角(失速)为50˚,明显高于上述三种类型网板,因此,本文认为该网板曲面和导流板结构可有效减缓流分离的产生。
4.2. 数值模拟验证
本文数值模拟计算过程中,数值模型边界条件及计算域大小均与模型试验保持一致 [17] [21],结果显示,数值模拟结果均与试验结果近似,模拟误差低于10%,数值模拟可以较好的预测网板升、阻力系数及升阻比等参数随冲角的变化趋势,且可以有效预测网板临界冲角、压力中心系数等水动力学参数,同时数值模拟可以较好的反应网板周围流场流态,因此本文认为数值模拟是网板水动力学研究的一种有效的方法。由于网板是在水中进行作业,在网板水动力研究过程中应当尽量保证流体介质与实际情况一致,因此本文建议在条件允许的情况下,利用水槽试验和数值模拟方法相结合的方式研究网板水动力学性能。
基金项目
国家重点研发计划(2020YFD0901205);农业农村部远洋与极地渔业创新重点实验室开放课题(2019HY-XKQ03)。
NOTES
*通讯作者。