基于OpenFOAM的数值波浪水槽研究.pdfVIP

第 54卷 第 3期 复旦 学 报 (自然科学版) Vo1．54No．3 2Ql午6月 Journa】ofFudanUniversity(NaturalScience) Jun．2015 文章编号：0427—7104(2015)03—0373—06 基于OpenFOAM 的数值波浪水槽研究 刘秀丽 ，段梦兰 ，高 攀 ，钟 超 (1．复旦大学 力学与工程科学系，上海 200433；2．北京高泰深海技术有限公司，北京 100011) 摘 要：使用设置造波边界法和质量源造波法建立了数值波浪水槽，采用阻尼消波法进行消波．设置造波边界 法根据行波的解析解或数值解来给定边界上的流速和波高值．质量源造波法是在水槽 中特定位置，将质量源项 加入到质量守恒方程中．使用该数值水槽模拟了线性波浪和二阶Stokes波，并与理论解进行了对比，结果显示 两种造波方法都能有效模拟行波．设置造波边界法必须指定速度和波面边界，不适用于波流相互作用等边界未 知的情况，而源造波法则能适用．但是由于源造波法需要在两端设定造波区域，计算区域更大，其造波效率低于 设置造波边界法． 关键词：数值造波；数值消波；波浪水槽；流体体积函数法 (VOF) 中图分类号：035；U661．1 文献标志码：A 波浪是最常见的海洋环境载荷．在港 口海岸和海洋工程中，通常会通过实尺试验、模型试验和数值计 算等方法来评估波浪诱导运动和载荷．但是，实尺试验成本很高且难以控制，另外，等待极端海况出现也 是不现实的．在工程应用中，通常采用模型试验．模型试验的缺点在于，当黏性力起主导作用时，从模型试 验的结果到实尺结果的转换有困难．同时，几何尺度和试验设备也会限制试验的可能性．随着高速计算机 和大容量存储技术的发展，数值计算在计算波浪诱导运动和载荷方面的应用越来越多．其中，波浪仿真是 计算波浪诱导运动和载荷的关键[1]． 在波浪仿真中，数值造波和消波是技术关键．模拟波浪的数值模型可以分为两类 ：基于势流理论和边 界元方法的数值模型和基于流体控制方程的数值模型．前一类方法是最早出现的一种数值水槽模拟方 法，Tanizawa和Kim等对相关理论和数值技术进行了较为全面的综述_l2]．这类方法基于流体无黏无旋 的假设，一般采用边界元法离散拉普拉斯方程和非线性 自由面边界条件，用混合欧拉拉格朗 日法来更新 自由面的运动．此类数值水池计算代价较小，但无法很好地模拟那些黏性效应不可忽略的物理现象，如波 浪远程传播和涡的演变等．后一类方法是基于Navier-Stokes方程或雷诺平均方程建立起来的黏性数值 水槽，它采用标记网格法(Marker—and—Cell，MAC)或流体体积函数法 (VolumeofFluid，VOF)方法捕捉 自由面 ]．王永学 (1994)在 国内最先开展了数值波浪水槽 的研究，采用基于N-S方程的NASA—VOF方法 建立了二维波浪水槽，造波边界处造波板的运动除了产生行进波外同时产生一个抵消反射波的局部波动 以消除二次反射的影响，出口边界应用 Sommerfeld型敷设条件防止透射波在边界的反射[5]．邹志利等 (1996)[、刘海青等 (1998)[]继续对采用 Sommerfeld辐射条件消波进行 了研究，发现单独采用 Sommerfeld辐射条件的消波效果并不理想．此后，多数研究均采用添加人工阻尼来进行消波．周勤俊等 (2005)提出了动量源造波法，该方法仅在动量方程中添加源项．该文章推导出了适用于RANS(Reynolds— AveragedNavier—Stokes)方程和VOF方法的动量源项，利用Fluent软件建立了数值波浪水槽 ，模拟了海 塘越浪过程 ]．李胜忠(2006)基于Fluent建立了数值波浪水槽，对比了推板式造波法和质量源造波法，发 现源造波法计算效率更高且更稳定_lg]．近年来随着OpenF()AM开源程序的出现和发展，国内越来越多的 人开始使用它进行计算流体力学 (ComputationalFluidDynamics，CFD)计算．查晶晶 (2011)基于 OpenFOAM，使用推板和摇板造波法建立了二维和三维波浪水池_l1．张博杰等(2012)基于OpenFOAM， 收稿 日期 ：2014-12—16 作者简介：刘秀丽(198