二阶入射边界条件在完全非线性数值波浪水槽中的应用.doc

?完全非线性深水波的数值模拟 周斌珍1, 宁德志1 2, 滕斌1, 宋伟华1 (1. 大连理工大学 海岸和近海工程国家重点实验室，辽宁 大连 116024；2. 上海交通大学 海洋工程国家重点实验室，上海， 200030) 摘 要：基于势流理论，并结合深水波质点运动从水面向下呈e指数衰减的特性，建立了完全非线性数值变深水槽模型，通过实时模拟活塞式造波机运动来产生波浪。采用时域高阶边界元法进行模拟，利用混合欧拉-拉格朗日方法和四阶Runge-Kutta方法追踪流体瞬时水面，应用镜像格林函数消除了水槽两个侧面的积分，在水槽末端布置人工阻尼层消除反射波浪。利用GMRES加速算法，提高了计算效率。利用所建模型对深水波进行了模拟研究，得到了稳定的波形，在造波板运动幅值较小时，与线性解析结果吻合良好，当造波板运动幅值较大时，体现出波浪的非线性特性。同时利用所建模型研究了造波板所在位置（上部台阶）水深对生成波浪高度的影响，由此可以选择合适的造波板所在位置水深及造波板运动幅值来得到所需要的深水波浪。 关键词：造波板运动；数值变深水槽；完全非线性；高阶边界元；深水波 中图文分类号： U661　 文献标识码：A 文章编号：0253-4193（2009）00-0000-00 1 引 言 近年来，由于中等和浅水深的开发技术已经成熟，以及当今世界大部分地区浅海油气资源已日渐减少，因而深海资源开发已经成为当今世界能源开发的热点和必然趋势。许多专家和学者对深水问题开展了各种理论研究工作，Newman[1]通过水平台阶理论分析了波浪从浅水到无限水深和无限水深到浅水的传播变化；Fenton[2] 通过解析研究得到了五阶深水波浪速度势和波面的解析解；Ferront[3]通过频域和时域方法研究了深水中非线性波浪辐射问题。然而，由于理论分析适用范围的局限性，开展深水问题的数值模型试验是很重要的。国内外专家开展了许多相关研究工作，且都希望建立一个与物理波浪水槽功能相当的数值模型。譬如Lee等[4]通过截取有限长度计算域来研究深水中运动物体的荷载；Fonseca 和 Soares[5]利用满足自由水面条件的格林函数对无限水深中船舶做大振幅运动问题进行了数值模拟研究；宁德志等[6, 7]利用源造波技术对无限水深的完全非线性波浪进行了数值模拟研究。然而，大多数数值模型中都是通过在入射边界处给定特定的理论波浪特性（如速度或波高）来实现造波，尽管这样可以令入射边界条件固定而避免网格的重新划分，使问题求解相对简单，但在入射边界给定一个与实际波浪非线性相匹配的理论速度或波高通常是很难的[8]，特别是对于造波板做大振幅运动的情况。而通过模拟造波板实时运动产生波浪则可以解决入射边界不匹配问题，最终真实模拟物理波浪水槽中的自然现象。如Bai和Eatock Taylor[9-11] 采用高阶边界元方法，在数值波浪水槽中模拟造波板实时运动产生波浪，分别对完全非线性规则波、聚焦波与垂直圆柱相互作用问题及完全非线性波浪对固定和漂浮结构物作用问题进行了模拟；周斌珍等[12]和Ning等[13]采用类似的方法建立了完全非线性数值波浪水槽模型，并对规则波、不规则波进行了模拟。 以往的波浪水槽大多以平底水槽为主，在水槽一端设置造波机，通过给定造波机的运动形式产生所需要的波浪。对于深水波而言，质点运动从水面向下呈e指数衰减，因此只需要把造波机设置在水面附近一定深度，只在表面造波，由此建立台阶式变深水槽。Zhou等[14]等利用特征函数展开法对台阶式变深水槽中推板式造波机造波问题进行了解析研究。 本文利用时域高阶边界元方法建立了一个完全非线性数值变深水槽，通过实时模拟造波板在上部台阶运动来产生深水波浪，且自由水面满足完全非线性边界条件。在每一时间步内，利用四阶Runge-Kutta法和混合欧拉-拉格朗日法对瞬时水面和入射边界面进行更新；为了避免由于网格运动变形而引起数值不稳定问题，在一定时间步内对网格进行重新划分，通过采用二次形状函数计算新网格的物理量和几何量；将镜像格林函数运用到所建立的边界积分方程中，进而消除水槽两侧面的积分；通过在入射边界面和台阶立面上网格沿水深方向逐渐变粗的方法，减少了计算所需要的网格；利用GMRES加速算法，提高了计算效率。利用所建模型对深水波进行了模拟，并与线性解析结果进行了对比验证，同时利用所建模型研究了造波板所在位置（上部台阶）水深对生成波浪高度的影响。 2 控制方程及初边值条件 考虑如图1所示的台阶式变深水槽，在静水面上建立笛卡尔坐标系，z轴向上为正，并于台阶立面重合，x轴向右为正。波浪沿x轴正方向传播。在水槽的左端装有造波机，水槽右端为吸