数值模拟孤立波在平直海岸上的传播和爬坡.pdfVIP

第 14卷 第 11期 中 国 水 运 VoI．14 No．11 2014钲 11月 ChinaWaterTransport November 2O14 数值模拟孤立波在平直海岸上的传播和爬坡 金 红 ’，苟大旬 ，房克照 。 (1．大连市环境科学设计研究院，辽宁 大连 116023；2．辽宁省交通厅港航管理局，辽宁 沈阳 ll0003 3．大连理工大学 海岸和近海工程国家重点实验室，辽宁 大连 116023) 摘 要：对文献 【lJ所建立的非线性浅水方程数值模型进行了改进 ，采用水位梯度法 (SGM，SurfaceGradient Method)处理底坡源项以提高计算效率 ，采用半隐格式处理水底摩擦项以提高模型的稳定性。通过解析解对模型 进行验证后，利用改进后的模型 ，对破碎和非破碎孤立波在平直海岸上的传播、爬坡和回落过程进行了数值模拟研 究。数值计算结果同实验数据、解析解以及文献中相关数值结果进行比较和分析，结果表明模型可用于计算孤立波 的传播和爬高问题。 关键词：孤立波；爬坡；非线性浅水方程 中图分类号：P731 文献标识码：A 文章编号：1006-7973 (2014)11-0137-04 一 一 、 引言 维守恒形式的非线性浅水方程为 海啸通常由于海底地震或者海底滑坡引起，产生伊始，海 +( 1 0 (1) 啸波波高可能较小，但其在向岸传播过程中发生的浅水变形致 ～ 使波高大幅度增加，最终导致海啸上岸后对海岸地区造成严重 a(du)+( ：+ ：) 跳 +I／-U了2 (2) ： 的破坏。譬如，1998年7月发生在巴布亚新几内亚的海啸对 。 u | ‘ “ 该地区造成严重的破坏，并导致2，000余人死亡，据 目击者称 式中：d=-h+ 水深；h为静水深；釉 波面；g为重 海啸波有4-15m 高l2J。2004年印度洋海底地震诱发的海啸， 力加速度；毛为底摩擦系数。 致使 l1个国家的225，000人死亡，海啸波达 30多米高。 Fang等采用高精度 Godunov类有限体积法对上述方程 海啸上岸后造成的损失程度同海啸爬高密切相关，许多学 进行求解。其中采用Roe格式计算控制体边界的流量，通过 者致力于这方面的研究，而孤立波和基于非线性浅水方程的数 应用状态插值法 (MUSCL)达到空间的二阶精度，同时使用 学模型通常被用于模拟海啸的传播和爬坡过程。孤立波是以固 非线性限制因子使得格式无震荡。时间上则采用二阶精度的 定速度传播的水体，水质点运动方向一致，并且传播过程中形 预报矫正格式进行积分。此外 ，通过对水底坡度项进行特征 状变化不大，这同海