孤立波在平坡海岸上爬高的数值计算.pdfVIP

第 16卷 第 6期 中 国 水 运 VoI．16 No．6 2016生 6月 China Water Transport June 2016 孤立波在平坡海岸上爬高的数值计算 张甲波 ，宫立新 ，李朋飞 。 (1．河北省地矿局秦皇岛矿产水文工程地质大队，河北 秦皇岛 066000； 2．河北省海洋地质资源调查中心，河北 秦皇岛 066000；3．中水东北勘测设计研究有限责任公司，吉林 长春 130021) 摘 要：本文建立了具有间断捕捉能力的Boussinesq波浪传播数学模型，对孤立波在平坡海岸上的爬高进行了数 值计算。对一组弱非线性和弱色散性Boussinesq水波方程采用均匀网格进行空间离散，采用MUSTA格式和高精 度状态插值方法计算控制体界面数值通量，其余项采用高精度差分公式，采用具有TVD性质的三阶龙格一库塔方法 进行时间积分。在对模型进行验证后，利用数值试验对孤立波在具有不同坡度海岸上的爬高进行了数值模拟和相关 分析。 关键词：Boussinesq波浪模型；孤立波；爬高 中图分类号：TV139．2 文献标识码：A 文章编号：1006—7973(2016)06—0097—04 一 、 引言 开数值计算研究，重点考察海岸坡度和海床摩擦对爬高的影 孤立波的传播特性和海啸波类似，其被广泛用于研究海啸 响。 波的传播和爬高过程…_。【1。众所周知，海啸波传至近岸区域， 二、数学模型 受海水变浅的影响，波高不断增大，加上增水的影响，会在海 1．Boussinesq水波方程 岸地区发生爬高，上岸波浪对海岸地区造成严重的破坏l4】。 以波面升高 和任意水深处速度 u 表达 的高阶 数值计算是研究海啸波爬高的重要手段，数值模拟是研究 Boussinesq水波方程为加【l 海啸波传播、爬高和越浪过程的重要手段。现有的海啸数值模 77，+[(+ )】 拟可分为两类：第一类是基于雷诺平均的纳维一 斯托克斯方 +[( 一 (z+ =。 程 (RANS)的数值模型，第二类非线性浅水方程模型，主要 包括全非线性浅水方程数学模型 (NLSW)和Boussinesq类 +甜“+g仉+z[ +( )]+ =0 (2) 模型。RANS详细地描述流场，且几乎适用于任何复杂形态的 ， 海底地形，但计算代价昂贵 (通常需要花费小时级别的机时模 式中，h为静水水深， 为 自由波面，d=h 77为当地水 拟秒级的实际流场)，目前很难用于工程实践。相对PANS模 深；g为重力加速度，z=一0．531h，下标t表示变量对时间的 型，NLSW模型相对简单，仅计算深度平均后的相关量，计算 偏导数，下标X和t分别表示变量对空间和时间的偏导数。t 效率更高，且已经发展出许多高精度数值格式，被诸多学者采 为依照曼宁系数的海底切应力，定义如下： 用作为控制方程2『。．5】。但致命的是，这种模型缺乏色散性， — ． (3) 仅能应用于模拟短距离海啸波传播，不能考虑如长距离传播、 绕射等短波类效应，严重限制了其应用范围。而Boussinesq 类模型同时具有色散性和非线性，能较准确的考虑海啸波的长