浅海围栏牧场设施工程数值模拟分析.docxVIP

浅海围栏牧场设施工程数值模拟分析 邹海生 朱烨 洪扬 李超 摘要：为了研究浅海养殖围栏设施的可靠性、安全性，以我国海洋工程的桩基设计及渔用网片水动力性能研究为基础，应用现代仿真技术与计算机数值模拟技术，对围栏设施阻力进行模拟计算与分析。研究表明，对单一网片采用数值模拟的方法进行受力分析是可行的，其误差范围为±3 N，采用数值模拟方法对围栏设施的水动力性能进行研究，得出随着桩柱数的增加、流速的加快，围栏设施的总阻力也逐渐增大，同时在一定区间内其总阻力的变化小，适当增加桩柱数或选择合适桩柱数，对围栏设施有很好的保护作用。 关键词：养殖围栏设施;桩柱;网片;数值模拟 中图分类号： S953.2文献标志码： A 文章编号：1002-1302（2020）08-0238-04 收稿日期：20XX-03-22 基金项目：中国水产科学研究院基本科研业务费专项（编号：20XXHY-ZD0603）。 作者简介：邹海生（1981—），男，安徽合肥人，工程师，主要从事海洋渔业装备设施研究。Tel：（021E-mail：zouhaisheng@fmiri.ac.cn。 淺海牧场式围栏设施与生态养殖符合渔业转方式、调结构的发展方向，是海水养殖结构调整产业升级的重要新兴技术，有广阔的产业应用前景。围栏牧场养殖设施是渔业装备与工程集成应用、相辅相成的重要体现[1-5]。目前该模式发展迅速，但存在规模较小、成本较高、养殖对象单一、技术不够成熟等问题;同时在围栏构建中普遍缺乏关于风、浪、流的基础理论作依据，且国内相关研究较少，在实际使用中存在着结构安全性的问题，可造成养殖生产的巨大经济损失等问题。 赵云鹏等对不同流速和不同冲角组合情况下的单片网衣周围的流场进行了数值模拟[6];王磊等介绍了桩柱式铜合金围栏养殖发展现状及围栏结构强度，并提出几点建议[7];陈昌平等研究了水流作用下金属网衣的水阻力，结果发现，网线直径加大会增加水的阻力[8]。现有研究大部分都是利用数值方法对网衣受力进行3个方向上的分析，而关于桩柱和网衣的模拟研究鲜见报道。为研究立柱的抗风浪性能、围栏网绳的受力情况，本研究先采用Fluent软件对单一网片水动力性能进行数值模拟[9-10]，并修正计算过程和有关参数，验证数值模拟的可行性，再对不同流速下不同桩柱围栏设施受力情况以及在同一流速下，不同桩柱围栏设施的受力情况进行计算机模拟，探讨不同桩柱对围栏设施受力的影响，旨在为围栏设施的总体设计提供可靠的技术支持。 1流体力学控制方程 海浪运动为不可压缩性流动，因此本研究中流体运动的控制方程采用不可压缩黏性流体Navier-Stokes方程组[11-12]： ρt+xi（ρui）=0 t（ρui）+xj（ρuiuj）=-pxi+ρgi+xjμuixj+Si。（1） 式中：ρ为流体密度;t为时间;xi，xj（i=1，2，3;j=1，2，3）为3个方向上的空间坐标;ui、uj为3个方向的速度分量;p为压强;gi为3个方向上的加速度;μ为动力黏性系数;Si为源项。 RNGk-ε模型具有较高的精确度和可信度，已被广泛应用到各种湍流模拟中。本研究采用 RNGk-ε 模型进行模拟，对于不可压缩波浪，其输运方程为 kt+xi（kui）=xjakveffkxj+Gk-ε εt+xi（εui）=xjaεveffεxj+G1εεkGk-G2εε2k。（2） 式中：veff为有效黏度;Gk为由于平均速度梯度引起的湍动能k的产生项;ε为湍流耗散率;常数ak=aε=1.39，G1ε=1.42，G2ε=1.68。 2单一网片水动力性能数值模拟及试验验证 2.1网片三维模型 利用计算流体力学软件Fluent对单一网片进行水动力性能数值模拟[13-14]，网片大小为480 mm×480 mm，网绳直径为4 mm，间距为40 mm，实际网绳为上下编织形状，但该形状在进行数值计算网格划分时比较繁琐，因此将网绳简化为相互穿透齐平形状，该简化操作对于数值计算结果的影响基本可以忽略，简化后的网片三维模型见图1。 2.2网格划分 计算域选取网前1倍的网片长度，网后4倍的网片长度，如图2所示。为了得到更加精确的模拟效果以及节省计算时间，在网格划分时，将网绳附近的网格进行加密，网格由网绳处向外逐渐稀疏，以保证网格有良好的过渡性，这样既能保证计算的准确性，又能减少计算时间，整体域网格及网绳附近网格见图3。 由于要计算得到网绳对水流的阻力，而在网绳壁面附近流场变化比较剧烈，因此在划分网格时网绳周围要进行边界层网格划分，这样能