三维边界元方法的波浪变形计算.pdfVIP

第四届全国海事技术研讨会文集 三维边界元方法的波浪变形计算 孙大鹏 李玉成 (大连理工大学海岸及近海工程国家重点实验室。 大连116024) 数为未知量的数值积分方程，进而对空间区域内理想流体有势运动的波浪变形问题进行数值计算．由于三维空问区域 的表面包括入射面，自由水面，底面以及岸壁等具有不同类型边界条件的边界面．因此，这种三维边界元模式涵盖了 折射、绕射及反射等诸多波浪变形因子，可适用于求解具有复杂岸壁形状。任意反射系数的非恒定水深区域波浪变形 问题． 关■髑Laplace方程．波浪变形，三绝边界元 前言 近岸区域，由于受水深地形变化的影响，加之岛屿、岸壁的屏障，使得自外海传来的波浪在这里 发生诸如折射、绕射及反射等一系列的变形，自由水面呈现为～种复杂的表面波动。若从有势运动的 角度来考虑和解析这种近岸波浪现象，无疑基于质量守恒原理建立起来的Laplace方程是其最基本的 控制方程。但由于数值求解．Laplace方程在计算机容量及计算机对方面存在着诸多限制．使得以往的 研究工作常常在引入某些特定波场条件的基础上，对三维Laplace-方程做一些适当的简化。如建立在 水深变化非剧烈区域上的缓坡方程(Berkhoff 射方程等。为克服上述方程的数值模式在评价任意空间区域波浪变形时的局限性，探讨和建立近岸区 域波浪变形计算的三维数值模式，具有着实际工程意义和广泛的应用前景． 1 基本方程和边界条件 波动水体应遵循质量守恒原理(鄙连续性方程)，而对于有势线性波动，波场中任意一点盯良只圳， 可以采用对空闻、时间可以分离变量的形式，给出三维波势函数中“月邑∥的表达式： 巾(t．y，乙f)；』圣i巾(_y，力exp(一f∞，) (1．1) Z∞ 其中：霄，为入射波波高； m为圆频率；F为重力加速度 因此，求解波势函数÷(x，Y，：)的基本方程即为如下的三维Laplace方程： 黝萨 +塑矿 塑铲 =O (1|2) ·19· 若三维空间区域Q由闭曲面r所围成，任意点肌幺以∥和点聪，Tl，‘)位于边界面r上，其间的 距离 则依Green定理，点彤伍月∥处的势函数值可由边界上点P(1，T1，‘)处的势函数值及其法线方向 导数值以下列积分方程式表示： 刮c加=小国羞(9丢警卜∽吐坩， ms， 著把近岸波浪变形区域分成等水深的外海区域和非恒定水深的计算区域，如图l所示，显然·在 等水深外海区域。波势函数应由入射波~|和岸壁影响造成的散乱波税这两部分组成，即 40u《x，y，z)；再(x，Y，z)+虬@，y，z) (1．4) 微小振幅波·在水深为^的波动区域内，(1．1)式中的三维波势函数÷(x，y，：)可分离出水深因子 鄙(1．4)式入射波孓和散乱波瓦可分剐展开成为： (1．5) 玩(z，y，力=fo(x，，)cosch。k＆b(z肋+h) ÷。(x，圳=^(x，)，．，—cosh磊k函(z+-h) (1．6) 式中t女为波数； 兀@，y)、A(x，y)为入射波和散乱波的平面波势函数；且 矗(x，Y)=-iexp{i(kcos，+如sing)} (1．7) 式中：p为入射波波向和x轴之间的夹角。 关于散乱波平面波势函数^0，力的二维Ht，lmholtz方程， 尝+磐“^：o (1．8) 国2。却2…。1” 对于散乱波．若加之以无远处的Soramcrreid辐射条件，即 缨．占(等嗡Ⅲ