三种湍流模型在阶梯溢流坝数值模拟中的比较.pdfVIP

第32卷第 7期 人 民 黄 河 VoI．32．No．7 2010年 7月 YELL0W RIVER Ju1．，2010 【水利水 电工程 】 三种湍流模型在阶梯溢流坝数值模拟中的比较 补舒棋 ，唐新军 (新疆农业大学 水利与土木工程学院，新疆 乌鲁木齐 830052) 摘 要：采用 standardk—s、RNGk— 、Realizablek—s3种不同的湍流模型并结合VOF法，对阶梯溢流坝气液两相流动 进行了数值模拟，将不同湍流模型所模拟的水面线、台阶速度场分布及台阶壁面压力分布与物理模型试验实测结果进行 了对比分析。结果表明：RNGk—s湍流模型能够较好地模拟阶梯溢流坝的水流特性。 关 键 词 ：阶梯溢流坝水流 ；VOF方法 ；standardk一8模型；RNGk一￡模型；Realizablek一￡模型 中图分类号 ：TV135．2 文献标识码：A doi：10．3969／j．issn．1000—1379．2010．07．063 随着计算流体力学的发展 ，国内许多学者开展了对阶梯溢 人 口边界条件采用速度入 口边界条件，入 口边界处压力均为大 流坝的数值模拟研究 ，并取得了较好的研究成果，如陈群、陈永 气压。由于出口水深未知，水气边界未知，因此出口边界采用 明 -3]等利用 standard 一s湍流模型结合VOF方法对阶梯水 压力出口边界，引水槽底部、WES堰面及阶梯面为固壁边界。 流进行了数值模拟。目前对阶梯水流进行数值模拟，较多的是 给定无滑移边界，近壁区的黏性底层用壁函数处理。 选择 standardk一 、RNGk—s或Realizablek— 模型中的一 阶梯溢流坝网格划分采用平铺非结构化网格的方式，网格 种结合VOF方法对阶梯水流进行数值模拟 ，但对这3种湍流模 采用四边形网格，在每一级台阶面上适 当加密网格 ，从而更好 型的模拟效果缺乏综合比较分析。笔者结合 VOF法，通过实例 地模拟台阶间的旋涡。压力、速度的耦合采用压力隐式算子分 来分析和探讨上述 3种湍流模型用于阶梯溢流坝数值模拟的 割法 (PISO算法)进行计算，PISO算法在计算过程中进行了相 适用性和计算效率，以期为今后开展阶梯水流数值模拟研究提 邻校正和偏斜校正 ，减少了迭代次数，可提高计算效率。笔者 供参考。 在模拟计算中，对于不同的湍流模型均使用相同的网格，以避 免网格对模拟结果的影响。 1 数学模型及计算方法 2 计算结果与对比 1．1 V0F法 VOF法是Hirt和Nichols于 1981年提出的处理复杂 自由表 通过VOF法和3种湍流模型的结合，模拟了实例模型水头 面的有效方法 。阶梯溢流坝面上水流是～种复杂的带有 自由 为9．7em、流量为0．02m ／s的情况，分别得到了水面线 、阶梯 表面的气液两相流动，因此可采用VOF法通过求解体积分数的 速度分布、速度矢量分布图及阶梯压力分布图。 连续方程来求解其 自由表面。F：g(x，，。，t)表示计算区域内 2．1 水面线对比 流体体积 占计算区域的相对比例 ，对于某个计算单元 ，F=1表 流量为 0．02m ／s时，3种不同湍流模型得到的模拟水面 示该单元被流体充满 ，F=0表示它是一个空单元，0 ’1表 线和实测水面