基于VOF 法的风暴潮洪水演进三维数值模拟 - 中国工程科学.PDFVIP

基于 VOF 法的风暴潮洪水演进三维数值模拟 １ １ ２ ２ １ 王晓玲 ，孙宜超 ，陈华鸿 ， 杨丽美 ， 孙蕊蕊 （１．天津大学水利工程仿真与安全国家重点实验室，天津 ３０００７２； ２．天津市水利勘测设计院，天津 ３００２０４） ［摘要］　为了提高风暴潮灾害应急处置能力，确保防潮安全，文章建立了耦合流体体积函数（ＶＯＦ）法的三维 非稳态水气两相流k －ε模型，采用等效糙率的方法处理城市密集建筑群，既考虑了其阻水作用，又考虑了其 蓄水作用。 针对天津市滨海新区海河与永定新河之间区域的风暴潮洪水演进数值模拟与分析，对１００年一遇 风暴潮洪水淹没情况进行分析，并对不同频率的风暴潮洪水的严重性进行了比较。 结果表明，随着风暴潮发 生频率的增加，风暴潮淹没范围逐渐减小，水深随着频率的增加是逐渐减少的，该研究为海堤安全管理、风暴 潮灾害的快速科学评估提供了理论依据和技术支持。 ［关键词］　风暴潮洪水演进；流体体积函数（ＶＯＦ）法；三维k －ε紊流数学模型；水深分布 ［中图分类号］　ＴＶ１２２＋．１　［文献标识码］　Ａ　［文章编号］　１００９－１７４２（２０１１）１２－００５６－０７ 产生的影响［６］ ；姚志坚等提出在溃坝洪水演进计算 1　前言 中用“等效糙率”模拟建筑群的方法及“等效糙率” 风暴潮是指由于大风以及气压急剧变化等因素 取值的水槽试验手段，并成功应用在某城市水库溃 造成的沿海或河口水位的异常升降现象。 风暴潮往 坝洪水演进研究中［７］ ；张大伟通过考虑社区和楼房 往导致沿海地区水位暴涨，造成严重的淹没损失。 内部的容水性，引入侵入水量的概念并给出其合理 ［１］ 风暴潮数值模拟的研究开始于２０ 世纪５０ 年代 ， 的估算方法，建立了能够模拟建筑物密集城区溃堤 国外一些学者采用数值模拟的手段研究了密集城区 水流运动的二维浅水方程数学模型，对哈尔滨市可 的溃堤洪水问题，Ｈａｉｄｅｒ 等和 Ｉｎｏｕｅ 等采用二维浅 能发生的溃堤洪水进行了计算研究［８］ 。 水方程分别模拟了法国Ｎｉｍｅｓ 城和日本Ｏｓａｋａ 城的 以上研究多是针对水库溃坝后洪水的数值模 溃堤大洪水，较好地模拟了阻水区域和蓄水区域，均 ［９］ 拟，而且多是一、二维 ，国内外基于风暴潮洪水演 取得了理想的计算结果［２，３］ ；Ｍｉｇｎｏｔ 等和 Ｓｏａｒｅｓ 等 进三维数值模拟的相关研究成果较少，缺乏对含有 通过城区溃坝物理模型试验对二维浅水方程数学模 复杂地貌的大范围洪水演进的数值模拟。 文章建立 型的精确性进行了验证，通过对试验数据和模拟数 了耦合流体体积函数（ｖｏｌｕｍｅ ｏｆ ｆｌｕｉｄ，ＶＯＦ）法的三 据的比较，发现洪水淹没水深、流速和演进时间均较 维非稳态水气两相流k －ε模型，采用等效糙率的方 吻合，然而由于地形条件的不确定性和数学模型模 法处理城市密集建筑群，重点对天津市滨海新区海 拟的局限性，仍然存在一些偏差［４，５］ 。 对于密集城 河与永定新河之间区域１００年一遇的风暴潮洪水淹 区内溃堤洪水计算方面，国内涉及较少。 张大伟等 没情况进行了数值模拟与分析，并对不同频率的风 在堤坝溃决水流数学模型中，采用真实地形法对建 暴潮洪水的严重性做了比较分析。 筑物进行处理，较真实地模拟了建筑物对溃决水流 ［收稿日期］　２０１１