非结构网格上的三维浅水流动数值模型-水科学进展.PDFVIP

第 17 卷 第 5 期 水 科 学 进 展 Vol17 ,No5 2006 年 9 月 　　 ADVANCES IN WATER SCIENCE 　 Sep . , 2006 　 非结构网格上的三维浅水流动数值模型 1 2 3 4 赖锡军 , 曲卓杰 , 周 　杰 , 申满斌 ( 1 中国科学院南京地理与湖泊研究所 , 江苏 南京 　210008 ; 2 山东省南水北调工程建设管理局 , 山东 济南 　250013 ; 3 河海大学环境科学与工程学院 , 江苏 南京 　210098 ; 4 清华大学水利水电工程系 , 北京 　100084) 摘要 : 针对当前复杂环境水流模拟的需求 , 建立了新型的基于特征型高分辨率数值算法的三维非结构网格浅水动 力模型 。模型采用有限体积法离散 sigma 坐标下的三维浅水方程 , 运用 Roe 黎曼近似解评估水平界面通量 。模型网 格拟合边界能力强 , 可根据需要局部加密 ; 格式数值性能优 良, 具有守恒性 、单调迎风性 、高数值分辨率等特性 。 同时 , 应用干湿判别法处理动边界 , 以适应浅滩地形漫/ 露过程模拟的需要 。封闭水池内部风生环流 、干河床上溃 坝过程和长江口实际潮流场的模拟从不同侧面展示了模型的特点 , 结果表明它能够准确地预测水流的三维流动结 构 , 而且计算简单高效 , 具有良好的数值稳定性 。 关 　键 　词 : 三维 ; 浅水 ; 非结构网格 ; 黎曼近似解 ; sigma 坐标 ; 动边界 ; 水动力模型 ( ) 中图分类号 : TV1312 　　　文献标识码 : A 　　　文章编号 : 2006 大型湖泊 、河口、海岸等水域的水动力特征研究是工程建设 , 区域资源 、环境开发 , 环境保护与规划等的 基础 。数值模拟技术的发展 , 为研究大型水体的水动力特征提供了快速方便的手段 , 当前 , 二维 、三维水流数 学模型已经开始在这类区域广为应用 。尤其是二维模型 , 在 20 世纪 90 年代以来 , 建立在非结构网格有限体积 法基础上的新一代特征型的高分辨率数值算法得到了发展 , 由于其优 良的数值性能 , 得到了迅速的推广应 用[ 1] 。在实际水域中 , 流动三维特性明显 , 尤其是具有密度分层的水体 。但是由于精细模拟大区域非恒定流动 仍存在计算代价过高等问题 , 基于静水压强假定的三维浅水流动的数学模型在这类水体计算中被广为采用 。这 类三维模型大多是在水平方向采用基于结构网格的有限差分法和三角形网格的有限单元法[2 ] , 而垂向则采用 sigma 坐标或垂向固定分层 。由于有限单元法计算量较大 , 为了拟合复杂的自然几何边界 , 又发展了曲线拟合 坐标 , 但是曲线拟合坐标对真实环境的拟合能力仍然有限。而且 , 采用的数值格式较为古典 , 未能充分应用当 前浅水流动计算的众多高性能格式 。 为此 , 本文在非结构网格上 , 采用有限体积法离散方程 , 运用 Roe 黎曼近似解评估界面通量 , 建立了 sig ma 坐标系下的特征