曲线坐标系二维带自由表面强紊动水流数值模拟 - 水科学进展.pdf

第 15 卷 第 6 期 水 科 学 进 展 Vol15 ,No6 2004 年 11 月 　　 ADVANCES IN WATER SCIENCE 　 Nov. , 2004 　 曲线坐标系二维带自由表面强紊动水流数值模拟 戴会超1 ,2 , 魏文礼3 ( 1 中国长江三峡工程开发总公司 , 湖北 宜昌　443002 ; 2 河海大学 , 江苏 南京 　210098 ; 3 西安理工大学水利水电学院 , 陕西 西安 　710048) ( ) 摘要 : 将追踪自由表面的流体体积 VOF 法应用于曲线坐标系下水流控制方程的求解中 ,计入流线弯曲对水流紊动 特性的影响 , 建立了垂向二维强紊动水流的曲率修正的紊流模型 , 并对溢流坝反弧段的紊流特性进行了数值模拟 。 κ ε γ 数值计算时 , 采用有限体积法离散水流的控制方程 ; 物理变量 , 如 : 压力 P 、紊动参量 、 、 等 , 采用交错方 t ( ) 式排列 交错网格布置 ,用 SIMPLEC 算法求解离散方程 。计算结果表明 , 得到的溢流坝反弧段的自由表面位置 、速 度场 、压力场 、剪应力分布和紊动能分布与实验结果吻合良好 。 关 　键 　词 : 曲线坐标系 ; 有限体积法 ; 自由表面 ; 强紊动水流 ; 数值模拟 ; VOF 法 ( ) 中图分类号 : TV13121 　　　文献标识码 : A 　　　文章编号 : 2004 随着工程实践的发展 , 高水头 、大流量水工建筑物 日益增多 , 与之相关的诸如空化 、空蚀等一系列高速水 力学问题 日益突出。泄水建筑物的反弧段易产生负压 , 导致水流不稳定或水流空化数过低 , 反弧段由于具有自 由表面和受显著曲率壁面的影响 , 离心力和重力都不容忽视 , 形成了十分复杂的水流现象 , 高速水流流经反弧 段往往造成严重的空蚀破坏 , 是工程界极为关注的研究课题 。传统的方法是以实验为主的 , 随着数值计算技术 的发展 , 数值模拟逐步成为强有力的研究手段 , 但是由于泄水建筑物流场复杂 , 并且具有不规则的几何边界及 自由水面 , 进行数值计算比较困难 。以势流理论为基础采用有限元法 , 文献[ 1 ,2 ]对于溢流坝面的流场进行了 数值分析 , 文献 [ 3 ] 则研究了三维过坝水流 问题 ; 随着紊流理论的发展 , 紊流模型在水力学中广泛应用 , ε 文献[4 ]采用极坐标系下 k 的模型对反弧段紊流场进行了模拟 ; 文献[ 5 ]采用紊流模型 , 并辅以VOF 法追踪 自 由表面对淹没水跃进行了