完全频散性、非线性波浪传播数值模式.pdf

完全频散性、非线性波浪传播数值模式1 吴中 河海大学交通学院，南京（210098) E-mail ： zwu@ ，wuhohai@ 摘 要：本文以洪广文提出的具有完全频散性非线性波浪传播方程为基础，建立了一个包括 线性与非线性模型、能够适用于深水至浅水、包纳流场、地形与水位缓慢变化、可以满足开 敞与多连通水域波浪传播的统一模式和计算方法。论文对深水至极浅水非线性波进行了较全 面的验证，计算模拟了波浪浅化、立波、三维干涉波、潮波、孤立波以及波浪越过潜堤、单 突堤绕射、双突堤绕射、岛式堤、Berkhoff 实验和多连通水域波浪传播等问题，还对特殊地 形条件下波浪发生港湾共振、布拉格共振进行了仿真；论文研究了水流中波浪传播折射变形 问题，Arthur 流场上非线性折射-绕射、水位缓变的转向潮流中波浪传播等缓变水流中波浪 问题，各类问题的较好解决验证和表明了本文模式与方法具有适用于深浅水的可计算性与非 线性。 关键词：非线性；频散性；波浪数值模拟；波流相互作用；波浪共振 中图分类号：O353 1. 引言 [1～5] 波浪模式包括波能平衡方程和波浪传播方程两大类 ，传播方程的发展经历了四个阶 段：线性波折射模型[1～19](1940)、联合折射绕射模型[20～32] (“缓坡方程”1972)、非线性长波模 型[33～37] (Boussinesq 方程、1976)和完全频散性非线性模型[38～40](近十余年) 。线性折射模型适 用于大范围深水水域，不能描述浅水非线性波，不适应浅水复杂地形产生波向线相交情况； 缓坡方程具有完全频散、可以描述波浪折射与绕射行为，但无法描述波浪的非线性；缓坡方 程抛物近似模式忽略了反射波、且绕射角受到限制[41～53] ；缓坡方程双曲近似模式引进虚拟 、 、 中间量[7 53 54] ，近似模式在非线性、准确性上都有一定的不足；Boussinesq 非线性长波方程 只具有弱频散性，并且无法在深水计算；Boussinesq 改进型方程虽然频散性有提高，但深水 计算仍受到一定限制。目前国内外在近海、海岸工程中，大多应用缓坡方程的等价方程，抛 物方程在商业软件中得到最为广泛的应用，绕射角受到限制的问题也未得到彻底的解决。从 以往波浪传播模式看，完全频散性与非线性未能较好地统一到一个方程中。1994 年， Nadaoka[38]提出了一个无流中无能量交换项的完全频散非线性方程，同年 Isobe[39]也提出一 个无流水域中无能量交换项的完全频散非线性方程，但形式极为复杂。它们方程的性能与普 适性也都在研究中。近年来洪广文[40]提出了缓变水流、地形、水位条件、且具有能量交换 因子的完全频散非线性方程，本文依据洪广文方程，建立了能够适用于深水至极浅水不同大 小计算范围、包纳流场、地形与水位缓慢变化、绕射角无限制、时步上计算结果准确合理、 具有波浪传播过程中能量交换因子、可以满足开敞与多连通复杂地形水域时域波浪传播统一 计算的模式和计算方法。 2. 计算方法 数值计算采用时间前差与空间中心差、采用无条件稳定的Crank-Nicolson 格式。基本解 法为预测-校正-迭代解法，对强非线性问题采用本文提出的“ 同胚线性的预测-校正-迭代方 1本课题得到国家自然科学基金40676053) 的资助。 - 1 - 法”和“ 同胚线性的交替推进迭代方法” 。它们是依据滤波与通量控制思想，利用低阶方程做 自适应滤波器与限制器的求解非线性方程的两种解法，尤其适用于求解强非线的细节问题。 边界条件主要采用洪广文[55]改进的Sommerfeld 型反射与辐射边界条件。 3. 数值模型验证 针对完全频散性、非线性方程，进行计算普适性研究。从深水(h/λ=1.0)至极浅水(h/λ=0.01) 验证非线性方程与算法的深水波(h/λ≤1.0)～极浅水长波表达