模拟孤立波在斜坡上爬高的新模型.pdfVIP

第六届夸国流体力学会议论文集 2Q01 上海 模拟孤立波在斜坡上爬高的新模型，) 刘桦 王本龙 何友声 张近‘ 上海变通太学t上海200030)‘(香避抽I市大学．香港 摘要针对斜坡上的波浪爬高计算问题，讨论并建立了适合于Bou8sinesq模型的严格的水边 界条件，进而引^水平坐标的伸缩变换．建立了数学模型．本文提出的计算模型不仅在动边 界处理上具有合理的理论基础，而且大量的数值实验表明本模型是有效的。应用本模型的数 值结果，讨论了色散对孤立波爬高的影响。 关键词孤立渡，哂而商矗面j}程，爬高 ————勺i 引言 。 为了认识和预报海啸这一海洋灾害对海岸带的影响，人们在过去十多年内对孤立波在海 滩上的爬高问题进行了较广泛的研究，涉及解析、数值模拟和物理模型实验等方面。尽管美 国和日本等国家已建立了较完善的海啸预警系统，包括编制以孤立波爬高理论计算为基础的 et 沿岸风险图，但是，有关的基础研究工作似乎还方兴未艾．如Liual(1994)o]、Kanoglu＆ 何条件)和水位梯度为零的提法，严格地来讲在理论上是不完备的。 在海岸工程界，为了 et 模型中实现，并且可推广至平面二维的爬高计算，这些模型毕竟在理论上回避了处理水边线 的难点，且得到的最大爬高值偏小。Zhang(1996)t61蛤出了水边线的几何条件、运动学条件和 动力学条件，但是仍有待进行广泛的计算验证，特别是推广至平面二维的工作有待继续深入。 此外，斜坡式防波堤仍然是一种广泛使用的海岸工程结构形式，研究浅水太波，如孤立波， 在斜坡上的爬高对海岸工程结构物的设计具有重要的实际意义。 etal et 在处理动边界的困难。例如，针对带自由表面的水动力学问题，Kim f1983)口】并口Liu a1(1992)t 8】采用基于样条插值的边界积分方法和混台的Euler-Lagranga自由面跟踪技术，对自 由表面的晃荡、飞溅以及斜坡上孤立波的变形和破碎进行了模拟。依据以下原理．位于自由 表面与固壁相交处的水质点一直沿固壁运动且占踞动边界位置，我们可以根据得到的固壁上 波的传播及其在直墙上的爬高问题。当然，若在动边界上考虑粘性效应，如引进接触线模型， et 问题会变得复杂，见Liual(1998)㈣。 本文将讨论应用经典Boussinesq方程计算斜坡上孤立波的爬高问题。首先使用坐标伸缩 得到一组新坐标系下的Boussinesq方程。进而采用有限差分方法，建立了数值模型，对陡坡 上海市科委青年科技启明星计划【98∞14。44)和上海市重点学科基金蕾叻项目 253 上的孤立波爬高进行了全过程模拟．并与现有的BIEM计算结果进行比较。 1基本方程 不妨取平底部分的水深h。为特征长度，√丙为特征流速，√ii万为特征时间。无量纲 31 的Boussinesq方程为，见Whitham(1974)11 等+兰[(^+f)“】=0 (1) ∽ ％ 加v-a“ar1，a3“ 研栅夏+云2i扩孬 (2) 其中，“为水深平均流速，f为水位，h为水深。地形如图1所示，Lo为孤立波峰到角点的 距离；h=h(x)为水深．平底部分水深为l：防波堤的斜率为1／L。 右侧自由边边界条件为 ^(z(f))+f(xO)，t)=0