缓坡方程的有限元解法及应用 - 水科学进展.pdf

第 20 卷 第 2 期 水 科 学 进 展 Vol20 ,No2 2009 年 3 月 　　 ADVANCES IN WATER SCIENCE 　 Mar . , 2009 　 缓坡方程的有限元解法及应用 1 1 ,2 2 1 郑 　俊 , 李瑞杰 , 江森汇 , 罗 　锋 ( 1 海岸灾害及防护教育部重点实验室 , 江苏 南京 　210098 ; 2 河海大学物理海洋研究所 , 江苏 南京 　210098) 摘要 : 利用有限元方法离散椭圆型缓坡方程 , 能适用于复杂区域 , 并很好地拟合不规则边界 ; 采用改进共轭梯度法 求解离散方程组 , 可以大大降低计算内存要求 , 提高计算效率 。利用结合上述两种方法的模式对规划的日照港区水 域进行了波浪数值计算 , 并将计算结果与物理模型试验值进行比较 , 结果表明: 该模式能适用于较大区域的波浪场 计算 , 并可以得到较好的计算结果 。 关 　键 　词 : 波浪 ; 数值计算 ; 缓坡方程 ; 有限元法 ( ) 中图分类号 : TV1392 　　　文献标识码 : A 　　　文章编号 : 2009 Berkhoff [ 1] 提出的缓坡方程考虑了波浪的折射 、绕射 、反射等联合效应 , 对复杂地形下的波浪计算具有较 满意的精度 , 已被广泛运用于近岸波浪传播计算 。该方程为椭圆型 , 可采用有限差分方法和有限元方法对其进 行离散求解 。李瑞杰等[2 ,3 ] 利用有限差分方法对椭圆型缓坡方程进行求解 , 得到了较好的结果 。有限差分方法 具有形式直观 、理论成熟 、精度可选 、易于编程 、易于并行等优点 , 但对不规则区域的处理较为繁琐 , 且对复 杂不规则边界拟合较差 。有限元方法相对有限差分方法的优越性主要体现在 : 能适应不规则区域 , 拟合复杂不 [4 ] [5 ] [6 ] [7 ] 规则边界等方面 。赵明等 、柳淑学等 、Bellotti 等 、Khellaf 等 利用有限元方法离散椭圆型缓坡方程对试 验地形或理想地形上的波浪场进行过求解 , 都得到了较理想的结果 。 在波浪计算中 , 一般要求每个波长内要求至少有 8 个节点 , 在大区域上对椭圆型缓坡方程进行离散将得到 庞大的稀疏系数矩阵 。由于受计算机内存和速度的限制 , 大型稀疏系数矩阵方程的求解比较困难 。经验表明 , ( ) [8 ] ( 直接法 如 Gauss 消去法 只能适用 10 个波长以内的计算域 。迭代方法 如J acobi 、GaussSeidel 和超松弛迭代 )