方腔回流区水流运动特性三维数值分析.pdf

第 14 卷 第 1 期 水 科 学 进 展 Vol14 ,No1 2003 年 1 月 　　 ADVANCES IN WATER SCIENCE 　　J an . ,2003 　 方腔回流区水流运动特性三维数值分析 张修忠 , 王光谦 (清华大学教育部水沙科学重点实验室 , 北京 100084) 摘要 : 建立了一个三维水动力学模型 , 对方腔回流区水流运动特性进行了模拟分析 。数学模型以三维浅水环流方程 σ 为基础 , 在垂向空间引入 坐标变换 , 采用高分辨率半隐有限元离散格式 。计算结果与解析解及水槽试验数据吻 合良好 。应用该模型对方腔回流区的水流运动特性进行了计算分析 , 进一步揭示了方腔回流运动的非定常非对称 不封闭特性 、流动结构表现为竖轴环流与立面环流相叠加 、流速沿垂线分布相对均匀等流动规律 。 关 　键 　词 : 方腔回流区; 环流 ; 有限元 ; 三维数值分析 ( ) 中图分类号 : TV 1332 　　　文献标识码 : A 　　　文章编号 : 2003 0100108 自然界的水流运动通常由主流和次生流构成 , 当次生流基本上脱离主流形成相对封闭的环流时 , 便产生了 回流 。如内河港池 、船闸引航道 、盲肠河段等口门附近普遍存在着回流运动 。回流流动形态上的封闭性是相对 的 , 回流区的水流和泥沙仍然与主流区的水流和泥沙进行着质量 、动量和能量的交换 。天然河道中 , 从主流进 入方腔回流区的泥沙通常在此大量淤积[ 1] , 给水利及航运等工程的运行带来很不利的影响 。事实上 , 方腔回流 区泥沙的大量淤积 , 正是由于回流运动的特殊结构所引起的。因此 , 全面地研究回流区流动特性不仅具有理论 意义 , 而且对于三峡工程泥沙问题的研究具有实际意义 。 对于回流区水沙运动的研究主要从丁坝回流研究开始 , 由于回流区水流结构的复杂性 , 以及受量测技术的 限制 , 20 世纪 90 年代前 , 主要研究丁坝引起的回流宏观尺度和局部冲刷 , 以便工程上选择合理的丁坝群间 距 。近年来 , 关于丁坝附近水流结构的精细量测仍然罕见 , 大多数侧重于研究丁坝头的局部冲刷深度及冲刷过 程 。 文献[2 ]对盲肠河段回流区水流运动结构进行了较为系统的研究 , 定性地论述了盲肠回流的形成原因和运 动机理 。文献[3 ]对方腔回流区的水流运动开展了水槽试验研究 , 并进行了平面二维数值模拟 ,侧重于回流区水 流的非恒定过程 。文献[4 ]采用平面二维水沙模型对方腔回流区的水沙特性进行了计算分析 。当水深相对较小 , 上 、下层水体掺混相对较弱时 , 平面二维计算是可行的。但是回流区水流运动 , 不仅存在平面上的环流运动 ,还 存在着上 、下层水体之间的强烈掺混 , 表现出典型的三维特性 。因此 , 对方腔回流区开展三维计算分析是十分 必要的 , 也只有三维分析才能够比较全面地揭示回流的空间结构和时间演化 。 1 　数学模型 11 　控制方程 在静压近似及 Boussinesq 假定下 , 描述不可压流体运动的守恒方程为