不同弯度水渠二维水流运动数值模拟研究.pdfVIP

第14卷第2期 宁夏工程技术 V01．14No．2 2015年6月 Technology Jun．2015 NingxiaEngineering 文章编号：1671-7244 2015 02—0132—06 不同弯度水渠二维水流运动数值模拟研究 杨程 北方民族大学土木工程学院，宁夏银川750021 摘要：为了研究水流对渠弯的影响，建立了水流运动的二维紊流数学模型，采用非结构网格有限体 积法，通过实验室连续弯道的算例，验证了数学模型的正确性，并对4种不同弯度水渠的水流运动进 行了数值模拟．结果表明，模拟值符合弯道水流运动的基本规律；在相同的水力参数下，随着渠弯角 度的增加，水流对渠弯左岸的侵蚀明显降低；验证了在渠弯设计中对渠弯角度的规定是合理的． 关键词：水渠；数值模拟；非结构网格；有限体积法 中图分类号：Tvl3l；0242．1 文献标志码：A 弯道式水渠在引水工程中应用很广，随着数值 1．2模型的离散 模拟的广泛应用，对弯道式水渠的研究也越来越 采用格心格式的非结构网格有限体积法，对式 多．尽管水渠的设计标准已非常完善，然而很多设 1 在控制单元上求体积分并利用高斯定理，则控制 计要求由经验而来，本文以水渠渠弯的角度需大于 方程可化为 90。为例，采用数值模拟手段对该设计要求进行分 + J A J J A 析和研究．二维紊流数学模型广泛应用于河流、海 』。华dA+』^鲫ndA 』。Fgradqb·ndAA优 岸等水流运动研究，国内外成熟的模型数量较多． j．s中dA， 2 由于水流运动呈现出明显的三维特性，因此许多学 式中：u为网格单元的流速矢量，几为网格界面上的 者致力于模拟三维水流运动的研究中，然而三维水 法向量． 流的复杂性，使得三维数学模型仍存在一些技术问 对式 2 进行离散，瞬态项采用一阶格式离散， 题【-1．本文采用平面二维紊流模型，利用非结构网格 源项采用线性化处理，具体过程见文献[3]． 有限体积法，首先通过1个实验室连续弯道实验， 进口边界：进口给定流量．其中k和s按照下 验证本文建立的数学模型是准确的，然后对4种弯 式计算【41 道式水渠的水流运动进行数值模拟研究． 孚； 1数学模型 而 丁3 IUi． 2，s 等等． 1．1数学模型 出口边界：出口给定水位．两岸边界：按照无滑 沿水深平均的平面二维RNGJj 一8水流数学模移处理． 型12]包括：连续性方程、动量方程、k一方程和￡一方 程，其通用方程为 2验证数学模型 优瓦,9 出 罟 +晰 通过1个实验室算例验证本文所建数学模型 hOP 十÷ 脚 +-茜- h训 -- x hF．Ox的正确性． 1 2．1模拟区域及网格划分 导 ^L罢 +Js中 z， ， ， o， o， 连续弯道示意图及弯道处3个断面的位置见 式中：通用变量少，扩散系数R和源项．s出，，， 的含 图1．弯道模型参数及实测结果引自文献[1]和文献 义及各系数的计算公式见文献【2】． L，s，糙率n 0．012，坡降五 [6]：进口流量为45．3 收稿日期：2013—12一17 作者简介：杨程 1983一 ，男，讲师，主要从事计算水力学研究， 电子信箱 yangchenglll@163．corn． 万方数据 第2期 杨程等：不同弯度水渠二维水流运动数值模拟研究 133 0．1％，出口水深h 9．4cm 流的数值模拟研究中． 图1连续弯遭不意图 采用Delaunayt5]--角化法对模拟区域进行网格 剖分 图2 ，共2031个节点，3771个单元． 2．2模拟结果分析 流场图见图3．以图l中设置的3个断面为例， 将各断面的断面垂向平均流速和水深分布的模拟 值与实测值进行对比 图4 ．约定：水流流动方向的 左侧为左岸，横轴表示距左岸距离 cm ，纵轴表示 垂向平均流速 m／s 和水深 cm ． 由图4可知：①实测值与模拟值较为一致；② 在弯道处从左岸到右岸，垂向平均流速逐渐降低， 水深则逐渐增大；③实验结果表明，本文建立的模 型、非结构网格有限体积法可用于弯道平面二维水 x／cm 图3连续弯道流场图及弯道处流场放大图 断面 伽 伽 枷 Eo≥ 伽 距左岸距离Icm 栅 断面2 瑚 瑚 姗 ；。L———一 羹7l 猢 ； _ 距左岸距离／cm 断面3 距芹卅爿，离／¨” 趴‘：岸距离／¨” ￡i：：“”为实删Ill值．“一”为艇拟值． 图4模拟值与实测值对比图 450 5 o 550 600 3渠弯水流运动数值模拟研究 x／ -nl 图2连续弯道网格剖分图及弯道处放大图 在水渠设计中，一般要求渠弯的角度大于900， 万方数据 134 宁夏工程技术 第14卷 本文设定渠弯半径为定值 r 4m ，针对4种水渠，即 900时，弯道最大流速逐渐向弯顶断面转移，且弯道