Nystrom法计算水波绕射波高的数值模型.pdfVIP

NQ．1 陕西科技大学学报 Feb．2009 · 154 · J0URNALOFSHAANXIUNIVERSITYOFSCIENCE TECHN0L0GY V0L27 文章编号：1000—5811(2009)01—0154—04 Nystrom法计算水波绕射波高的数值模型 薛 晓，王连堂 (西北大学数学系，陕西 西安 710127) 摘 要：研究了一般等截面柱体的水波绕射问题中波高的计算．应用线性小振幅波理论将水波 绕射 问题转化为一个二维的Helmholtz方程 ，再采用Nystrom方法来求解二维的Helmholtz 方程．通过求得的数值解与解析解的对比，说明作者给出的数值模型计算简单且有较高的精 度． 关键词 ：波浪绕射；Helmholtz方程；Nystrom方法 中图分类号 ：TV139．2；024 文献标识码 ：A 0 引言 在海港建设中，常常建造防波堤来掩护港域免受波浪的直接侵袭，港区波浪是决定防波堤的主要条 件．波浪在传播过程中遇到建筑物和地形变化时会发生绕射与折射，对港区波浪的分布产生明显的影响， 因此研究波浪的折射和绕射的特性对海港的设计具有十分重要的意义[1≈]．本文的研究假设入射波为线性 微小振幅波，海底平缓近似为等水深，即波数为常数，波浪传播过程中只有单纯绕射，不发生折射现象．当 结构物(防波堤，码头等)为近似露出水面的等截面柱体时，数学模型可以简化为二维的Helmholtz方程． 在文献[3，4]中应用边界元方法给出了数值模型，而本文则应用Nystrom方法给出了数值模型，通过这种 数值方法求得的数值解与解析解在结构物为圆柱时的对比验证，说明该方法计算简单且有很高的精度，并 且本文给出了一般等截面柱体的水波绕射模型． 1 基本方程的建立 在流体为无粘，无旋，不可压缩条件下，可用势流理论得到域 内任一点的速度势为： (z，Y，，￡)一6(x，Y，z)e-~ 线性水波绕射问题可由下式表示： △ 一 0 一d≤ ≤ 0 地～ + 一0 2—0 Ot2十 g 一 0 一 7一 7一一去 z_一uo 一 0 ～ =-- d 一 一 0 物面上 辐射条件 无穷远 收稿 日期：2OO8—1O一19 作者简介：薛 晓(1983-)，女，陕西省渭南市人，在读硕士生，研究方向：数学物理方程反问题 第 1期 薛 晓等：Nystrom法计算水波绕射波高的数值模型 ·155 · 其中刁为自由表面在静水以上的高度，d为静水深度，，z为固壁上的单位法向矢量，方向指向流体．由 线性假定，水域内各点总波势庐(z，y，)可分解为人射波势 和由障碍物扰动而引起的向外散射的散射波 势 ，即： (z，y，)一 。(z，Y，)+ (．z，Y，) (7) hk z~ h ( 入射波势 可写为 (，，)===z() (，)，其 中z()一A—cos[ d)]A _一 g ，