复式断面明渠二次流的数值模拟.pdfVIP

第16卷2期 应用基础与工程科学学报 Vo1．16，No．2 2008年4月 JOURNAL OF BASIC SCIENCE AND ENGINEERING April 2008 文章编号：1005-0930(2008)02-0296-08 中图分类号：TV131．1 文献标识码：A 复式断面明渠二次流的数值模拟 ‘ 梁爱国 ， 槐文信2 (1．广东省电力设计研究院，广东广州510600；2．武汉大学水资源与水rgn-程科学国家重点实验室，湖北武汉 430072) 摘要：采用Launder—Ying雷诺应力模型对复式断面明渠二次流流动进行数值模 拟，和实测资料对比，验证了数值计算的合理可靠性，并进一步分析了不同滩地 水深时复式断面明渠二次流运动的流动结构和特性．结果表明，二次流的阻滞导 致主流流动减缓；二次流主要产生在明渠固壁拐角和近岸区域；二次流对主流速 度等值线有强烈影响，使其沿二次流流动方向凸起，在滩槽结合部表现最明显； 随着滩地水深的增加，主槽涡强度相对减小，滩地涡强度相对增大，滩槽结合部 主流速度等值线在更大范围向自由水面凸起，主流高速度核心区域也逐渐由主 槽中部向滩槽结合方向移动；在滩地水深h／H=0．5左右时二次流最强． 关键词：复式断面明渠：二次流；应力模型；纵向旋涡 几乎所有的天然河道都存在浅滩，这种由主槽和浅滩构成的渠道称为复式断面明渠． 当洪水发生时，滩地上流动较慢的流体会向主槽内输运，这就产生了大量的流体横向交 换，由于速度不一致而形成的剪切层影响到湍流的流动结构并形成许多纵向旋涡．这一现 象已经在许多实验中被观测到，Tominaga和Nezu 、Naot和Nezu_2 J、Shiono和Feng_3 等 测量发现，二次流对动量输运交换和流动结构有重要影响． 近年来，许多学者对复式断面明渠二次流运动进行了数值计算和研究．工程中常用的 标准k-8模型是建立在湍流各向同性假定基础上，不能准确反映二次流运动，其主要缺陷 在于假定雷诺应力与时均速度梯度成正比，涡粘性系数由当地湍流结构特性确定，而与流 体的物理性质无关．为了克服这一难题，改善计算效果，人们采用不同的湍流模式进行了 大量数值实验．早期，Launder和Ying_4 J、Naot和Rodi_5 采用雷诺应力模式模拟了矩形断 面明渠二次流运动；在此基础上，Kawahara和Tamai。。 、Naot和Nezu 分别采用修正雷诺 应力模式模拟复式断面明渠二次流运动．Pezziga_8 J、Lin和Shiono_9 采用简化Speziale_l 0_ 非线形k-8模型成功模拟了复式断面明渠中二次流运动，并分析了二次流对浓度分布的 影响．Cokljkat和Younis【11]运用完全雷诺应力输运模式、Thomas和willi锄s_l 运用大涡 模拟分别对复式断面明渠中二次流运动进行了数值计算，研究了采用不同的数学模型模 拟二次流的准确程度，并得到了二次流动的一些特性．Kearney¨ 采用三种模式(Launder 和Ying模型、Naot和Rodi模型和Speziale非线形k-8模型)对比研究复式断面明渠中二 收稿日期：2006-08-02：修订日期：2008-01—15 基金项目：国家自然科学基金项目资助 作者简介：梁爱国(1975一)，男，博士． No．2 梁爱国等：复式断面明渠二次流的数值模拟 297 次流运动对浓度分布的影响，认为Launder和Ying模型能够更好模拟二次流运动．Shiono 和Kearney【J 也采用Launder和Ying模型对复式断面明渠中二次流运动对浓度分布的影 响进行计算，得到较好的结果．SHAO Xuejun等¨ 也采用三种不同的Reynolds应力计算 方法，包括采用Launder和Ying、Naot和Rodi的代数应力模型及运用非线性k-e模型计算 紊动黏性系数的SY模型，针对复式断面顺直明渠流动进行模拟，结果表明Launder和 Ying模型与SY模型数值模拟计算得到的二次流结构和紊动应力分布与试验数据均能达 到趋势上的符合． 本文亦采用Launde