闸坝水力学特性的三维数值模拟.docVIP

闸坝水力学特性的三维数值模拟 　　【摘 要】闸坝水力学特性的三维数值模拟采用标准k-ε双方程紊流模型及基于水气两相流的VOF方法，运用FLOW3D软件对闸坝进行数值模拟，得出不同水位下的闸坝泄流能力，以及不同部位的闸室及消力池水面线、流速、压力等重要水力学要素。表明数值计算一定程度上可以达到模型试验的效果。 　　【关键词】水气两相流 闸坝 数值模拟 泄流能力 水力学要素 　　1 前言 　　某闸坝共设2孔泄洪闸，1孔泄洪冲沙闸，泄洪闸及泄洪冲沙闸过流面净空尺寸均为12.0m×21.0m（净宽×净高），泄洪冲沙闸后布置深4.5m，长36m消力池。厂房段布置于枢纽中部靠左，为河床式厂房。闸坝坝顶高程494.0m，泄洪闸段闸室高程470.0m。因库区有公路桥及铁路桥，闸坝泄流能力至关重要，而坝址上下游地形较为复杂，对泄流能力有一定程度的影响，规范公式无法精确算出。鉴于模型试验周期长，成本较大，故运用大型流体计算软件FLOW3D对其进行数值模拟计算，以求得闸坝泄流能力及水面线、流速、压力等重要水力学要素。 　　2 数值模型的建立 　　2.1 基本方程 　　本计算采用k-ε紊流数学模型，引入适用于分层两相流的VOF方法求解自由水面，Hirt和Nichols[1]提出的VOF法是目前处理带自由表面分层流问题的较理想方法。鉴于闸坝水流条件并不复杂，采用标准k- e紊流模型计算，其连续方程、动量方程和k，ε方程分别为： 　　连续方程： 　　（1） 　　动量方程： 　　（2） 　　k方程： 　　（3） 　　方程： 　　（4） 　　式（2）中， 和 分别为体积分数平均的密度和分子粘性系数。P为修正压力； 为紊流粘性系数，它可由紊动能k和紊动耗散率 求出： 　　（5） 　　式（5）中， 为经验常数，取 =0.09。 和 分别是 和 的紊流普朗特数， =1.0， =1.3。 和 为 方程常数， =1.44， =1.92。G为由平均流速梯度引起的紊动能产生项，它可以由下式定义： 　　（6） 　　引入VOF模型的 紊流模型方程（1）、（2）、（3）、（4）与单相流的 模型形式完全相同，只是密度 和 的具体表达式不同，它们由体积分数加权平均值给出，即密度 和 是体积分数的函数，而不是一个常数，它们可以由下式表示出来： 　　（7） 　　（8） 　　式（7）、（8）中， 为水的体积分数， 和 分别为水和气的密度， 和 分别为水和气的分子粘性系数。通过对水的体积分数 的迭代求解， 和 的值都可以由（7）、（8）求出。 　　2.2 计算条件 　　2.2.1 闸坝模型 　　取坝轴线上游100m至坝下游160m的范围进行模拟计算，闸坝枢纽采用大型三维软件INVENTOR建模，生成STL文件后，导入Flow3D；河道模型通过AutoCAD生成的坐标点导入Flow3D。之后二者通过坐标转换对接，组合成整体模型，进行水力学的数值仿真计算。模型如图1～图2所示。 　　图1 闸坝三维模型图（下游看） 图2 闸坝三维模型图（上游看） 　　2.2.2 计算区域 　　实体地形及闸坝在FLOW 3D软件中按障碍物考虑。坐标轴规定如下： 　　X轴――顺水流向，上游指向下游为正；Y轴――垂直水流向，右岸指向左岸为正；Z轴――竖直向，竖直向上为正，符合右手螺旋定则。 　　坐标原点为坝横0+000.00与坝轴线的交点处。FLOW 3D中网格的剖分范围即计算区域。计算区域包括固体部分和流体部分。固体部分的范围：上游地形段、上游铺盖段、闸室段、下游铺盖及消力池段、下游地形段。流体部分的范围：上游地形至下游地形的稳定水体。网格剖分情况见图3和表1。因体型较为规则，故建立一个模块分析。 　　图3 闸坝及附近地形网格剖分 　　表1 计算区域剖分情况表 　　区块 网格剖分范围 网格数 网格总数 　　X向（m） Y向（m） Z向（m） X向（m） Y向（m） Z向（m） 2116400 　　block1 -100～160 -58～90 460～505 260 148 45 　　2.2.3 边界条件 　　边界条件的设定见表2。 　　表2 边界条件参数表 　　区块 X Min边界 X Max边界 Y Min边界 Y Max边界 Z Min边界 Z Max边界 　　block1 Specified pressure Specified pressure Symmetry Symmetry Symmetry Symmetry 　　边界条件说明：Specified pressure为压力边界，以符号P表示。block1： X Min边界条件为水流入口，F fraction设为1， X Max边界条件为水流出口，F fractio