抛物型缓坡方程的变分及数值模拟.pdfVIP

维普资讯 第 22卷 第 1期 海 洋 学 报 VoI_22．No．1 2O00—0l ACTA OCEANOLoGICA SINICA January，2000 抛物型缓坡方程的变分及数值模拟 ／7歹／2‘z 一 |30 林 钢 ． (大连理工大学海岸与近海工程 国家重点实验室，大连 116023) (华南理工大学交通学院，广州． 510641) 邱 大 洪 (大连理工大学海岸与近海工程国家重点实验室，大连 116023) 摘 要 对线性水波的折射 一绕射 问题应用变分原理，对非等深、具有缓玻和不连续 的底坡导出了一种修改的抛物型缓坡方程近似模型，可预测三维地形上波浪的折 射 一绕射．同抛物型缓坡方程的线性方程进行 了对 比．通过数值模拟方法进行数值 求解，表 明本方法可用于地形条件下的波浪折射 一绕射 问题． 关键词 銮坌堡堡．盔业 垣盟 绕射 修改地塑型堡垫左整 中圉法分类号 ：P731．2 1 引言 BerklIoff1【】使用近似方程结合经典Galerkin方法，将垂 向积分引入研究线性 水波折射 一 绕射的近似模型．许多学者相继提出了各种模型来模拟不同的地形引起的波浪折射 一绕射． 由于地形变化的捕获问题对线性水波可 由拉普拉斯方程结合边界条件和辐射条件作为控制方 程 ．对边值 问题的近似 已经提出了大量的方法 ． Kirby[·]对速度势使用WKB方法导出了抛物型缓坡方程，提出了 由缓慢变化 (缓坡 )地 形组成的附加迅速变化的小振幅分量处理的扩展缓坡方程 ． Mi1es[]表 明了对非线性 自由表面流动如何存在变分原理，并 由Luke_5提 出的变分原理 出发将其修改应用到线性 自由表面水渡问题，导出了缓坡方程 Chamberlain等[6]应用变分原理得到了修改的椭圆型缓坡方程和扩展椭圆型缓坡方程，指 出了在底坡不连续情况下，为使质量流量连续，必须满足跳跃条件 ． 本文应用变分原理对线性 自由表面水渡问题推导 出一种修改的抛物型缓坡方程，通过对 具体三维地形上的波浪折射一绕射问题的研究同抛物型缓坡方程比较，得到修改的抛物型缓 本文于 1998—05—24收到．修改稿于 1998419—24收到 第一作者茼舟：林 钢 ．男。43岁．高级工程师，博士，从事港 lit程和海岸及近梅工程研究 维普资讯 海洋学报 22卷 坡方程具有形式简单、计算简便、适应地形广等特点． 2 修改的抛物型缓坡方程 假设流体不可压缩、无旋流动，取正交笛卡儿坐标系，坐标系建在未扰 自由面 =O上，z， Y分别为水平方 向，=轴正向向上 ． 我们选择拉格朗 日密度函数 为 ￡=『[++号(V k 式中，中为速度势；f为 自由表面位移． 使用哈密顿原理，假定在时间间隔r内对函数，的小的变分是稳定的。则 最，=ff8『￡ d o， (2) 由此可得到描述水波的控制方程．将式(1)代入式 (2)。得 s，=一!盯f[号cV ++]f：s+[号cV中++=1：一a一 [中·f一 +1：s中一[中·一1一s中一 V ]8中d 曲d 0． (3) 我们假定 中 ： (z，Y，=)o0s ，