连续弯道水流模拟中二次流修正效果评价-水科学进展.pdf

第２７卷 第２期 水 科 学 进 展 Ｖｏｌ．２７，Ｎｏ．２ ２０１６年３月　 　 ＡＤＶＡＮＣＥＳ ＩＮ ＷＡＴＥＲ ＳＣＩＥＮＣＥ Ｍａｒ．，２０１６ ＤＯＩ： １０１４０４２／ ｊｃｎｋｉ３２１３０９２０１６０２０１２ 连续弯道水流模拟中二次流修正效果评价 １ ２ ３ ４ 周　 刚 ， 姚仕明 ， 秦翠翠 ， 贾　 鹏 （１ 中国环境科学研究院国家环境保护河口与海岸带环境重点实验室， 北京　 １０００１２； ２ 长江科学院水利部江湖治理与防洪重点实验室，湖北 武汉　 ４３００１０； ３ 清华大学水沙科学与水利水电工程国家重点实验室， 北京　 １０００８４； ４ 环境保护部环境工程评估中心国家环境保护环境影响评价数值模拟重点实验室， 北京　 １０００１２） 摘要： 为了选择适用于复杂连续弯道的二次流修正方法， 选取了线性方法中两种典型的计算模式， 通过在正交曲 线坐标系二维浅水动量方程中增加扩散应力项开发了考虑二次流影响的平面二维水流模型。 基于不同复杂度的４个 连续弯道试验的水位和流速分布资料， 通过理论分析和数值试验对比定量评估了传统二维模型和两种二次流修正 方法在连续弯道水流模拟中的效果。 测试结果表明， Ｄｅｌｆｔ３Ｄ模型的二次流修正方法自由度较高， 适用于不同复杂 度的连续弯道水流模拟， 而Ｌｉｅｎ模型二次流修正方法适用于微弯或中弯， 应用于中弯时需要慎重选用， 不适用于 急弯连续弯道模拟。 两种方法比较， 建议连续弯道水流模拟中优选Ｄｅｌｆｔ３Ｄ模型的二次流修正方法。 关键词： 二次流； 扩散应力； 连续弯道； 二维模型； 急弯 中图分类号： ＴＶ１４３　 　 　 文献标志码： Ａ　 　 　 文章编号： １００１⁃６７９１（２０１６）０２⁃０２６６⁃１２ 天然河道多以连续弯道为主， 尽管三维模型能够模拟弯道水流的三维环流特征， 但水深平均二维模型在 长时期大尺度特别是弯道演变研究和工程实践中仍然是重要的辅助工具。 传统的平面二维模型通过增加动量 ［１］ 在水平方向的交换系数 考虑弯道二次流影响的方法并不理想， 有可能导致水沙运动计算结果出现定性错 ［２］ 误。 自从Ｆｌｏｋｓｔｒａ 得出平面二维模型需要考虑扩散应力项解决弯道水流模拟效果不佳的问题以来， 出现了 线性和非线性两大类方法。 然而现有的弯道水流修正模型均存在不同程度的假设且通过单一弯道单一涡推导 而来， 测试算例和所得结论仍以单一弯道为主， 其特点及适用范围也存在诸多争议。 因此， 如何准确、 高效 地实施二次流修正并能应用于现阶段工程实践具有重要的现实意义。 线性方法通过水深平均动量方程增加扩散应力项进行二次流修正。 主要有３种求解扩散应力项的方法： ① 直接假定水流纵向和横向的垂向流速分布公式［３⁃１０］ ， 通过积分平均流速与垂向流速分布差值计算， 纵向 ［３⁃４］ ［５⁃６］ 流速垂向多是对数分布 或是指数分布 ， 横向流速垂向多是线性分布， 不同研究者区别在于流速垂向分 布公式选择差异。 该类方法已应用于微弯、 急弯单弯道和微弯 “Ｓ” 型双弯道。 ② 求解二次流强度的输运方 程［１１⁃１３］ ， 通过二次流强度、 水深及平均流速等变量关系式计算， 该类方法已应用于 １８０°和 １９３°急弯单弯道 水流以及天然微弯河段水沙模拟。 ③ 求解流向涡量输运方程［１４⁃１７］ ， 通过横向流速不均匀分布产生的涡量、 水深及平均流速等变量关系式计算， 该类方法已应用于 “Ｓ”、 “Ｗ”、 “Ｕ” 型以及Ｋｉｎｏｓｈｉｔａ变曲率连续急弯 弯道。 ［１８］ ［１９］ 非线性方法相对较少， Ｊｉｎ和Ｓｔｅｆｆｌｅｒ 、 Ｙｅｈ和Ｋｅｎｎｅｄｙ 通过联立求解水深平均连续方程、 动量方程 以及两个动量矩方程来考虑