动网格技术数值模拟挑流冲刷过程.pdfVIP

　 辑第 16 卷第 2 期　　　　　 水 动 力 学 研 究 与 进 展 　　　 . , . 16, . 2 A Ser A V o l N o 　　2001 年 6 月　　　　　　　JOU RNAL O F H YDROD YNAM ICS　　　　　June, 2001 文章编号: (2001) 动网格技术数值模拟 挑流冲刷过程 杨建明, 　　吴建华 (南京水利科学研究院, 江苏南京 210029) 　　摘　要: 　实际挑流泄洪消能时冲坑及堆丘的形成是一动态过程。已有对挑流冲刷的数值模 拟多考虑平衡冲刷坑情况。本文采用动网格技术模拟挑流冲刷时因河床边界变化引起的流场变化 过程, 控制方程组由非正交同位网格有限体积法求解。计算结果表明, 本文方法是有效的, 正确 地反映了挑流冲刷过程中流场结构及速度、压力等参数的变化。 关　键　词: 　挑流冲刷; 数值模拟; 体积守恒定律; 动网格 中图分类号: 　TV 863　　　文献标识码: A 1　引言 挑流消能是高坝泄洪消能常用的一种方式。对于我国在建及计划中的大批深山峡谷高坝 工程, 其泄洪水头高达 200～ 300 , 此时挑流方式仍可有效消煞余能; 但挑流水舌对下游河 m 床的巨大冲击力, 则有可能危及大坝和岸坡安全。故而, 挑流的冲刷作用一直是人们关心的 重要问题。 [ 1 ] 对于挑流冲刷的研究由来已久。陈椿庭 提出了关于高坝挑流消能和局部冲刷深度的估 算方法, 主要由水垫消能率的观点, 根据冲坑单位体积消能率及其几何形状二者来估计冲坑 [ 2 ] 深度。这是一种经验性很强的计算方法, 目前在工程规划与初步设计阶段仍有应用。柴华 应 用热膜技术对三峡溢流坝下游冲坑流场进行了量测。近年来开始采用数值模拟方法来研究冲 刷坑内的水流结构。文献 [3 ] 采用局部流线迎风有限元法数值模拟了冲坑水流。王艳明等[ 4 ] 运用正交变换技术将不规则冲坑流场区域变换为矩形区域, 对文献 [2 ] 所做的试验工作进行 [ 5 ] 了数值模拟验证。刘沛清等 应用大涡模拟对冲坑内的水流结构进行了研究。但这些工作一般 均着眼于平衡冲刷坑, 对于挑流冲刷过程中, 因边界运动而引起流场结构流速及压力等参数 变化的数值研究工作似乎尚未有报道。 由于冲刷过程中河床边界的运动, 计算区域是随时间而变化的, 则贴体计算网格必须跟 随边界一起运动。如果坐标系统的单位矢量保持不变, 唯一的变化是守恒方程中的对流项中 收稿日期: ( ) 作者简介: 　杨建明 1974～ , 男, 助理工程师, 硕士。 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 杨建明等: 动网格技术数值模拟挑流冲刷过程 157 会出现相对速度。在移