明渠岸边横向取水口三维数值计算.pdfVIP

水 利 学 报 2004 年2 月 SHUILI XUEBAO 第2 期 文章编号：0559-9350(2004)02-0119-07 明渠岸边横向取水口的三维数值计算 曹继文，陈惠泉，贺益英 (中国水利水电科学研究院 水力学研究所，北京 100044) 摘要：本文是在明渠岸边横向取水口的试验研究基础上的配套数值模拟。采用雷诺平均的N—S 方程，以各向异性的 雷诺应力模型(RSM)封闭方程组，用SIMPLEC 算法对压力—速度场求解。计算结果与试验实测结果吻合较好，反映了 此项模拟的可信性，可籍以扩大、分析一些运动规律。数值计算结果启示了横向取水口水流的复杂三维流动特性， 据之对取水口附近水流各区的划分、水流的分离和反向流动、螺旋流以及床面剪应力的分布等复杂水力和泥沙输移 物理现象有了进一步认识。 关键词：明渠；横向取水口；数值模拟 中图分类号：TV671 文献标识码：A 明渠横向取水是一种常见的岸边取水方式，在火 /核电厂的冷却水取水工程中得到广泛的应用，进行 过大量的试验研究。通过对横向取水口的试验研究[1～ 3] ，得到横向取水口三维流动特性的示意性描述如图 1 所示。当上游来流接近取水口时，由于横向取水的 分流作用，在取水口门的上游产生横向流动，近岸水 体流向分槽，与连续流向下游的水流间形成分流曲 面。进入分槽的水流受剪应力和离心力的不平衡作 用，产生顺时针的二次环流，类似于弯道水流出现的 环流，形成了进水口水域复杂的三维流动。口门前下 游近底层则出现滞流区，有明显的向底流动，在滞流 区下游产生逆时针的二次环流。上述流动特点对泥沙 输移的潜在效应是来流底沙随着底部水体分流进入 图1 横向取水口水流三维流动的特征 分槽，一定条件下，产生淤积；口门的下游近岸侧则 出现冲刷。这种口门水域的冲淤现象，Barkdoll[4]在1997 年试验室试验中有详细的观测。 对此直角横向取水的水力现象，已有不少系统试验成果，在数值模拟上，近年也取得一些进展；Chen [5] [6] [7] 和Lian 、Shettar 和Murthy 用的都是二维紊流模型。Issa and Oliveira 首次用三维紊流模拟T 型分 叉结构，用壁函数的标准k～ε模型封闭三维雷诺平均的Navier Stokes 方程。Neary 和Sotiropoulos[8] [9] 进行了三维层流的数值研究；Neary 和 Sotiropoulos 采用 k～ω紊流模型对横向取水口的内部水流进行 数值模拟，计算中考虑了高低雷诺数对壁函数的影响。 前人的数值模拟基本上是基于二维或三维的各向同性简化模型，难以正确反映进水口水域复杂的三维 流动，尤其缺乏实测资料的详细验证。本项工作便是针对这些问题，作进一步的深入。首先是模拟上采用 收稿日期：2003-09-08 作者简介：曹继文(1969-)，男，辽宁沈阳人，博士生，高级工程师，主要从事水工水力学方面的研究。 119 水 利 学 报 2004 年2 月 SHUILI XUEBAO 第2 期 了更符合实际的各向异性的雷诺应力精细模型；计算结果利用文献 ［１０］提供的大量实测三维流速 资料进行检验。在获得良好验证的基础上，扩大计算参数范围，探求并认识一些试验中不易取得的横向取 水运动规律，为继后的取水防沙工程措施提供科学依据。本文可以认为是文献［10］的配套和继续。 1 模型的建立 1.1 数学模型 模化后的雷诺应力模型(RSM)的基本控制