二维流场模拟矩阵标识法.pdfVIP

第 30 卷第 3 期 河 海 大 学 学 报 Vol . 30 No . 3 2002 年 5 月 JOURNAL OF HOHAI UNIVERSITY May 2002 二维流场模拟的矩阵标识法 1 1 2 李光炽 ,王船海 ,周晶晏 ( 1. 河海大学水文水资源及环境学院 ,江苏 南京 　210098 ; 2 . 河海大学交通与海洋工程学院 ,江苏 南京 　210098) 摘要 :采用二分步法 ,从积分型方程出发 ,在有限控制体上建立守恒型差分格式 , 对二维浅水波方 程进行求解. 求解线性代数方程组采用矩阵标识法 ,提高了计算效率 ,特别适合于流域大尺度水流 的模拟. 计算网格采用无序编码方法 ,灵活方便. 算例表明 ,本方法能够模拟宽广湖泊一类大尺度 的二维水流. 关键词 :二维水流 ;数值模型 ;矩阵标识法 ( ) 中图分类号 :TV131. 3 　　　文献标识码 :A 　　　文章编号 :1000 1980 2002 030080 05 二维流场的数值模拟研究由来已久 ,取得了十分丰富的成果. 数值模拟方法主要有三大类 :有限差分法 、 有限单元法和有限分析法. 其中有限差分法应用最为普遍. 常用的有限差分法有 ADI 法和剖开算子法 ,这些 方法又可分为显式和隐式两种格式. 显式方法编程简单 ,但稳定性受到 Courant 数的限制 ,时步长较小 , 当地 形变化剧烈时 ,还会出现“轴化”现象 ,使计算结果违背物理规律. 采用隐式方法求解 ,将形成巨大的高阶方程 [ 1] 组 ,对于河道型的计算区域 , 由于有一个方向的尺度相对较小 ,可采用矩阵追赶法 求解. 但对于大范围的计 算区域 ,当两个方向的尺度都较大时 ,这种方法也难以实施. 特别是对于宽广的计算区域 ,通常采用显式的方 法. 剖开算子法求解使得方法简化 ,易于实现. 为保证计算的稳定性 ,可采用隐式的方法. 隐式方法首要问题 是采用高效的存储和简捷的求解. 本文试图探讨这种方法. 1 　基 本 方 程 描述二维水流运动的基本方程是浅水波方程 : Z U V + + = 0 t x y U U U Z u2 + v2 2 ν + u + v + g h + g U - c W W - f V = U ( ) t x y x c2 h w x 1 V V V Z u2 + v2 2