动网格技术数值模拟-中国科技论文在线.pdfVIP

A辑第 16卷第 2期 水 动 力 学 研 究 与 进 展 S A，vd．16，N。．2 2oo1年 5月 JOURNAL OF HYDRODYNAM ICS J ，2001 文章鳊号 2001)02—0156—08 动网格技术数值模拟 挑流冲刷过程 杨建明， 吴建华 (南京水利科学研究院，江苏南京 210029) 捕 要： 实际挑流泄洪消能时冲坑及 堆丘的形成是一动态过程 ；已有对挑 流冲刷的数值模 拟多考虑平衡冲剧坑情况。本文采用动网格技术模拟挑流冲刷时固河床边界变化引起的流场变化 过程 ，控制方程组由非正交 同位 网格有 限体积法求解。计算结果表明，本文方法是有效 的，正确地 反映了执流冲刷过程中流场结构及速度、压力等参数的变化 ： 关 键 词： 执流冲刷 ；数值模拟 ；体积守恒定律；动网格 中围分类号： TV863 文献标识码 ：A 1 引言 挑流消能是高坝泄洪消能常用 的一种方式 。对于我国在建及计划中的大批深山峡谷高坝 工程，其泄洪水头高达 2∞～300m，此时挑流方式仍可有效消煞余能；但挑流水舌对下游河床 的巨大冲击力 ，则有可能危及大坝和岸坡安全。故而，挑流的冲刷作用一直是人们关心的重要 问题 。 。。 对于挑流冲刷 的研究 由来 已久 。陈椿庭 。提 出了关于高坝挑流消能和局部冲刷深度 的估 算方法，主要由水垫消能率的观点，根据冲坑单位体积消能率及其几何形状二者来估计冲坑深 度 。这是一种经验性很强的计算方法 ，目前在工程规划与初步设计阶段仍有应用 。柴华一应用 热膜技术对三峡溢流坝下游冲坑流场进行了量测 。近年来开始采用数值模拟方法来研究冲刷 坑 内的水流结构。文献 [3 采用局部流线迎风有限元法数值模j!：{了冲坑水流 王艳明等 运 用正交变换技术将不规则冲坑流场区域变换为矩形区域，对文献 Ez_q所做的试验工作进行了 数值模拟验证。刘沛清等 应用大涡模拟对冲坑内的水流结构进行了研究。但这些工作一般均 着眼于平衡冲刷坑 ，对于挑流冲刷过程 中，因边界运动而引起流场结构流速及压力等参数变化 的数值研究工作似乎 尚未有报道。 由于冲刷过程中河床边界的运动，计算区域是随时问而变化的，则贴体计算网格必须跟随 边界一起运动 如果坐标系统 的单位矢量保持不变 ，唯一 的变化是守恒方程 中的对流项 中会出 收稿 日期 2000—07—24 作者简介： 杨建 明(1974~)，男 ．助理工程师 ，硕 士。 转载 中国科技论文在线 杨建明等：动网格技术数值模拟挑流冲刷过程 lS7 现相对速度 。在移动边界 问题的计算 中，这个相对速度如何定义是一个关键 ，处理不好时常会 引人误差 ，产生人为质量源 ，最终使计算失败 。但若在计算过程中始终保证所谓的体积守恒定 律得到满足 的话 ，则可避免这一情况的发生 Thomas和 Lombard 。将体积守恒定律与质量、动量及能量输运方程列在一起 ，确定为在 动 网格 中数值求解流动问题时的基本控制方程组 ，并认识到必须将其与其它守恒方程 同时求 解 以避免产生人为质量源 Demirdzic和 Peric：q提 出了一种将体积守恒定律结合到有 限体积 法中的方法．井在文献 E83中计算了一个移动边界的层流流动问题，但对时间的离散采用了 一 阶精度的臆式 Euler格式 本文引入这一方法，计算挑流冲刷过程 的紊流流动问题，井对时 间和空间的离散都采用二阶精度格式。 2 基本控制方程与数值求解方法 考虑由