直立式防波堤迎浪面波压力计算研究.pdfVIP

第 46 卷第12 期 中国海洋大学学报 46 (12): 126~ 132 2016 年 12 月 PERIODICAL OF OCEAN UNIVERSITY OF CHINA Dec. • 2016 直立式防波堤迎浪面波压力计算研究来 于龙基1 ，董 胜IH ，段成林1 ，张华昌2 (1.中国海洋大学工程学院，山东青岛 266100; 2. 海军工程设计研究院，山东青岛 266100) 摘 要: 以 Navier-Stokes 方程为控制方程，使用 VOF方法追踪自由液面，由连续性方程和动量方程推导出源函数。基 于非线性被浪理论和 UDF对 FLUENT进行二次开发，在数值水槽的各个功能区将不同的源项添加到动量方程中，实现动 量源造波和消波的功能。此法可以有效地消除波浪在水槽出流边界的反射及披浪遇到结构物后在入射边界形成的二次反 射，从而得到稳定的波动场。本文将这种方法应用于直墙防披堤上波浪压力的数值模拟，分析了造披和1肖波的有效性，对 直立式防波堤波浪力进行了数值研究。为了得到更精确的结果，采用结构化网格和非结构化网格相结合，在直墙防波堤附 近使用三角形网格，在其他区域使用四边形结构化网格，为了较好地捕捉自由液面，在液面上下一个波高的范围内对网格 进行加密。将计算结果与二阶stokes 被理论解进行对比，该方法计算结果与理论解一致，表明该方法能够产生历时较长的 稳定的非线性规则波浪，能够很好地吸收破浪遇到建筑物产生的反射披，可以有效消除建筑物反射波的影响，为实际工程 方案的设计提供参考，有较高的实用价值。通过计算得到极端高水位，设计高水位，设计低水位时候的各个测点的最大正 向波压力和最大负向波压力，并将计算结果与物理模型试验数据和经验公式结果进行对比，结果吻合较好。并且由于该数 值方法可以考虑防波堤迎浪面不同倾角的工况，其适应范围比经验公式范围更广，因此可以为海岸结构设计提供依据。 关键词: N-S 方程; VOF方法;动量源;直墙防波堤;波压力 中固法分类号: P753 文献标志码: A 文章编号: 1672-5174(2016)12-126-07 001: 10. 16441/j. cnki. hdxb. 引用格式: 于龙基，董胜，段成林，等.直立式防披堤迎浪面被压力计算研究[1].中国海洋大学学报(自然科学版) , 2016 , 46 (12): 126-132. YU Long-Ji , DONG Sheng , DUAN Cheng-Lin , et al. Calculation of wave pressure on upright section of vertical breakwater [1]. Periodical of Ocean University of China , 2016 , 46 (1 2): 126-132. 消波功能和低通滤波功能。周勤俊等[7J 利用FLUENT 近年来，有限元数值分析软件的出现，增加了海岸 工程数值水槽计算的可行性和有效性，成为目前研究 商业软件，将人射波场作为人工的分布源项加人到动 波浪与海上结构物相互作用的重要手段。而数值波浪 量方程中，提出了一种适用于 VOF 方法的造波和消波 方法。林鹏智等[8J 通过基于 VOF 方法追踪自由面的 水槽研究的一个关键问题是要有效地消除波浪遇到建 筑物产生的反射波。有学者进行了深入探索。例如， 二维数值波浪水槽对直立防浪堤进行数值模拟研究。 Brorsen 和 Larsen[lJ 提出了适合于边界积分方程的非 董志和詹杰民[9J 基于不可压缩流体方程并采用 VOF 线性波浪的源项造波方法。 Arai 等[2