Youngs-VOF方法模拟溃坝水流演进.pdfVIP

第 2 8卷 第 5期 水 电 能 源 科 学 V ol. 2 8 N o . 5 2 0 1 0 年 5 月 W at er Resour ces an d P ow er Ma y. 2 0 1 0 文章编号: 1000 7709( 2 010) 05 005 1 03 oungs VOF 方 法 模 拟 溃 坝 水 流 演 进 方 杰1 周 杰1 张 健2 ( 1. 华东勘测设计研究院, 浙江 杭州 3 10014; 2. 河海大学 水利水 电学院, 江 苏 南京 2 1009 8) 摘要 : 针对追踪溃坝洪水 自由面变化一直为计算流体力学 的难点问题, 采用 oun gs V OF 方法追踪流体 流场 自由面、交错网格离散求解区域、人工压缩方法求解粘性不可压缩流体控制方程, 并通过数值计算模拟 了溃坝 洪水瞬变现象。结果表 明, oung s V OF 方法追踪运动 交界面 与 MA RT IN 的试 验结果 吻合, 且略优 于 H ir t VOF 法, 为溃坝洪水的模拟提供 了新途径 。 关键词: 运动交界面; oun g s V OF 方法; 溃坝洪水; 数值模拟; 粘性不可压缩流体控制方程 中图分类号: T V122+ . 4 ; O 353 文献标志码: A 当坝体或 围堰突然决 口倒塌时, 水流由上游 急速流下, 这一非恒定流为带有运动 自由面 的瞬 变流 。自由面的追踪一直为流体力学的难点与重 点, 通常采用刚盖假定及函数拟合获得 的数值结 果显然与实际出入较大, 运动交界面 的数值模拟 效果对整个流场 的数值计算结果影响也很大 。目 图 1 交错网格单元变量 前, 追踪 自由面 的方 法 很 多, 如 标 记 点格子 法 Fig . 1 Ele me nt v ar ia b le of s t agge re d grid ( MA C) 、相域法 、零等势面方法( LEVEL SET ) 、 n+ 1 n+ 1 P i+ 1, j - P i, j t g x - - F UX - F U + VISX ( 1) 流体体 积 输送 法 ( VOF ) 等, 其 中 oun gs V OF xi+ 1/ 2 [ 1] n+ 1 n 法 原理简单、计算量小而得到了广泛应用。鉴 Vi, j + 1/ 2 = Vi, j+ 1/ 2 + 此, 本文采用 oung s V OF 方法追踪溃坝洪水 自 n+ 1 n+ 1