基于NASA—VOF模型的二维溃坝水流冲击丁坝流场模拟.pdfVIP

第 27卷第 4期 安 徽 工 程 大 学 学 报 V01．27．NO．4 2012年 12月 JournalofAnhuiPolytechnicUniversity Dec．，2012 文章编号：1672—2477(2012)04—0065—05 基于 NASA—VOF模型的 二维溃坝水流冲击丁坝流场模拟 廖 斌 ，陈善群 (安徽212程大学 建筑工程学院，安徽 芜湖 241000) 摘要：建立了一种基于VOF 自由面追踪方法的NASA-V0F数值模型，并运用该模型模拟溃坝水流的演进过 程，将数值模拟结果与实验结果进行比较，结果表明NASA—VOF数值模型可 良好地实现二维溃坝水流 自由 面的追踪．之后运用 NASA—VOF模型数值模拟二维溃坝水流冲击丁坝流场 ，并与实验结果进行 比较 ，结果表 明NASA-VOF模型可较为精确地模拟二维溃坝水流冲击丁坝流场 的演化过程． 关 键 词 ：NASA-VOF；丁坝 ；溃坝水流 ；自由面 中图分类号：TV131．2 文献标识码 ：A 多年来溃坝问题一直是研究工作者和水利工程师所关注的问题．当河道中存在码头、丁坝 、桥墩等构 筑物时，溃坝水流的演进将会受到阻碍 ，并在构筑物附近形成绕射和反射流场 ，严重影响洪水的下泄；另 外 ，构筑物本身也会受到溃坝水流的强烈冲击，构筑物 的强度和稳定性将受到严 峻的考验．丁坝是现代治 河一个重要的工程措施 ，其 自身的稳定性和安全性非常重要 ，因此模拟溃坝水流冲击丁坝流场有着非常重 要 的工程意义． 随着计算机及数值模拟技术的迅速发展，数值模拟已成为研究丁坝问题一个重要的手段．目前针对丁 坝问题的研究主要集中在非 自由面条件下的丁坝绕流 ，国内外许多研究者对此进行了研究．李国斌_1]采用 紊流三维 是．￡模型闭合雷诺方程组 ，利用刚盖假定计算非淹没丁坝的流场 ；周宜林 应用大涡模拟方法研 究了丁坝附近的水流特性 ；彭静L3等将非线性紊流模型应用于丁坝的绕流场流动模拟中，并和线性模型进 行了比较．本文拟建立一种基于 VOF 自由面追踪方法 的NASA-VOF数值模型，并运用此模型数值模拟 溃坝水流冲击丁坝过程，在含有 自由面的条件下研究丁坝绕流流场的演进过程 ，为溃坝水流对丁坝冲击的 灾害提供较为准确 的预报． 1 NASA—VOF模型 1．1 控制方程 NASA—VOF模型的控制方程为不可压缩流体 的N-S方程和连续性方程 ： 瓦+十 +十 一一一一丢-a32x+gg+十、aa3Zu2。aa3Zvvu2~，．，’ L(1) +十““+十 一一一一丢 +十g+十c嘉+十 ，，， Lc2z -o- ~(Ou)+ ()一0， (3) 式 中：t为时间；“，为流体速度分量 ；P为流体压强；i0为流体密度 ；u为流体运动粘滞系数 ； ，g 为体积加 速度分量；为部分单元体参数，其值在 0～ 1之间． NASA—VOF模型引入 VOF函数 F(x，Y，f)来表示 自由面，它的思想是追踪 自由面的流体体积分数 ， 代表网格 中某类流体 占整个网格体积的百分 比．若 F一1，则该网格单元充满流体 ；若 F=0，则该网格单 元不含流体 ；若 0 F 1，则这一网格单元必定含有 自由面．所以通过求解 F值 ，即可确定 自由面在哪些 收稿 日期 ：2012-05-03 基金项 目：安徽高校省级科研基金资助项 目(KJ2012B016、KJ2012B020)；安徽工程科技学院引进人才科研启动基金资助项 目(2OO9YQQ009)；安徽工程大学青年科研基金资助项 目(2o12YQ27) 作者简介：廖 斌(1985一)，男，江西抚州人，实验师，硕士． 安 徽 工 程 大 学 学 报