非正交曲线坐标下二维水流计算的SIMPLEC.pdfVIP

水 　　利 　　学 　　报 　2003 年 2 月 SHUILI 　　XUEBAO 第 2 期 文章编号 : (2003) 非正交曲线坐标下二维水流计算的 SIMPLEC 算法 1 1 2 3 吴修广 , 沈永明 , 郑永红 , 王平义 ( 1 大连理工大学 海岸和近海工程国家重点实验室　辽宁 大连　116023 ; 2 中国科学院广州能源研究所, 广东 广州　510070 ; 3 重庆交通学院河海工程系　重庆　400074) 摘要 : 本文采用Laplace 方程坐标变换方法生成正交曲线网格 , 并对浅水流动的控制方程进行坐标变换 , 方程离 散时采用 B 型交错网格。利用“水位扫描法”结合壁面函数法来处理移动边界 , 用 SIMPLEC 算法解非正交曲线 ε 坐标下的 k - 双方程紊流模型 , 修正了由网格的非正交性引起的误差。通过对美国 Colorado 洲 Fall River 的资料 进行流场验证 , 计算结果与实测资料基本符合 , 显示了本模型在不规则水域计算中的实用价值。 ε 关键词 : 坐标变换 ; k - 紊流模型 ; 水位扫描法 ; 壁面函数 ; SIMPLEC 算法 中图分类号: TV1314 文献标识码 : A 随着经济发展和社会进步 , 水利工程建设的步伐也在进一步加快 , 其中港航建设、大坝建设中的 泥沙问题以及近来倍受世人关注的水污染问题已经成为制约水利发展的瓶颈问题 , 弄清河流、湖泊、 海洋中水动力因素 , 是解决以上问题的重要基础。近年来 , 数学模型已逐步取代物理模型实验成为研 究水流的重要手段 , 而浅水流动模型是处理大区域流场的一种非常有效的模型。它属于非线性方程 组 , 在目前只能用数值方法求解 , 因此 , 有必要研究一种简单、高效的方法来求解浅水流动问题。自 Patankar 和 Spalding[1 ] 发展了 SIMPLE 算法以来 , 该方法被广泛应用于不可压缩流场的数值模拟 , 而且 [2 ] [3 ] [4 ] 该方法还得到了进一步的发展 , 主要有 SIMPLER 算法 、SIMPLEC 算法 、SIMPLEX 算法 和 SIM [5 ] PLET 算法 等。这些模型均成功地应用于速度 —压力耦合的流场计算 , 深度平均的浅水流动模型是 在静压假定下导出的 , 一般流体模型中的速度 —压力耦合也就转换成浅水流动模型中的速度 —水深 耦合[6 ] 。 天然河流、海湾的边界曲折、地形复杂 , 采用坐标变换是解决问题的途径之一。目前多数 N - S 方程的坐标变换中 , 流程全部采用逆变分量 , 这样就增加了方程的复杂程度。于是忽略掉方程中的非 正交项 , 利用正交变换下的方程进行数值求解[7 ,8 ] 。对于具有复杂边界的海湾及弯曲的河流, 坐标变 换中很难保证每个点都正交