海岸碎波拍的计算15t∫5x-力学学报.PDF

第 1 期 董国海等 : 海岸碎波拍的计算 海岸碎波拍的计算1) 董国海 　　　　　邹志利 (大连理工大学海岸和近海工程国家重点实验室 , 大连 116023) 　　摘要 　研究了一维碎波拍的计算方法. 波浪控制方程由欧拉方程在短波周期上平均而得到. 为了考 虑低频波浪在海岸上的破碎 , 采用了 WAF ( Weighted Average Flux) 方法求解长波方程. 对由辐射 应力和海底摩擦力引起的方程非齐次项采用了时间算子分裂法进行了处理 , 并对由此产生的误差进 行了修正. 计算了单色波和双色波垂直入射海岸的情况 , 并与实验结果进行了对比. 关键词 　碎波拍 , 破碎波浪 , 低频波浪 引 　　言 碎波拍 ( surf beat) 是海岸处伴随入射波群而产生的低频波浪 (周期为 30 秒以上达几分 ) [1 ] 钟 , 大量研究表明其对海岸泥沙输运 , 港口波浪和系泊船运动都有重要作用 . 但关于其产 生机理仍不完全清楚. 60 年代 Longuet - Higgins 和 Stewart 提出的“辐射应力”的概念[2 ] , 可 用来解释波群中产生的非线性长波 Set - down. Set - down 在海岸处反射是碎波拍产生的一个 重要原因. Symonds 等[3 ]提出了碎波拍产生的另一机制 , 即波浪破碎点在空间和时间上的周期 性变化 , 其作用类似于造波板造波 , 从而产生低频波浪. Kostense[4 ] 的模型实验定性地证实了 以上现象的存在. Roelvink[1 ]求解以流量表达的一维低频波浪方程 , 研究了碎波拍对海岸变形 的作用. 通过坐标变换将计算域变换为固定域 , 处理了海岸处动边界. 本文将用于求解短波方 程的 WAF 方法[5 ]用来求解一维长波方程. 对方程中与辐射应力和底摩擦力有关的非齐次项 , 采用了时间算子分裂法进行处理 , 并对由此产生的误差进行了修正. 该方法将破碎波处理为 “激波”, 采用了 TVD 差分格式 , 对自由表面的陡然变化 (波浪破碎) 有很强的分辨能力 , 因 而可以考虑低频波浪在海岸处的破碎现象. 计算对单色波和双色波进行 , 并与实验结果进行了 对比. 1 　波浪方程 考虑波群垂直传向海岸的情况 , 如图 1 所示. 假定沿海岸方向地形变化是均匀的 , 这样方 程沿水深平均后所研究的问题是一维的. 对欧拉方程沿水深和短波周期平均 , 得 ζ ζ 2 d ρ ( ρ ) τ ( ) u d z + u + p d z - p b h + b = 0 1 t ∫ x ∫ d x - h - h ρ ζ τ 式中 , , h 和 u 分别为流体密度、静水水深和流体速度 ,