二维耦合模型在溃堤计算中的应用.docxVIP

二维耦合模型在溃堤计算中的应用 运用数值模拟法是对水库洪水风险分析的有效工具。另一方面，它可以直接为水库异常面积的防洪计划提供必要的洪水信息。另一方面，我们可以模拟人工大坝，在水库后蓄水区的洪水发展，对下游河流的防洪影响，制定系统的洪水规划方案。在长河段的洪水模拟中,为提高计算效率,一般采用一维数学模型;对于溃堤洪水的演进,多采用二维数学模型,以更好地反映物理量在平面范围的分布。将一维和平面二维模型耦合来模拟河道洪水与溃堤洪水,可以充分利用两者的优势,同时也避免了整个计算区域均采用二维模型而带来的计算量过大问题。 近年来,一维、二维耦合模型有了广泛的应用。Dhondia等、Dushmanta等分别采用荷兰代尔夫特水力学研究所开发的Sobek-Rural软件包和有限差分法耦合模型对河流漫顶洪水进行了模拟;Liang等采用动态链接技术进行一维、二维模型的耦合;张大伟等采用基于Preissmann格式的一维模型和非结构网格Roe格式的二维有限体积模型进行溃堤洪水耦合计算,同时对松花江哈尔滨段松北分洪区利用进行了模拟;槐文信等通过曲线坐标系下的有限分析法求解二维水动力方程,利用一、二维耦合模型模拟了渭河下游河道及洪泛区洪水演进;李云等将采用特征线法求解的一维模型与基于OS-FEM格式二维模型进行耦合,完成了淮河临淮岗段的洪水模拟。本文将基于黎曼近似解Godunov格式的一维、二维数学模型进行溃堤洪水的耦合计算,模型在模拟洪峰流量大,有急流、缓流及临界流等多种流态的复杂计算区域时较传统方法具有优势,张大伟等采用重叠计算区域方法建立了基于Godunov格式的一维、二维溃坝水流耦合数学模型,但是将一维河道模型和二维溃堤模型相结合的耦合溃堤模型还少见报道。模型通过堰流公式计算溃口处的流量过程,同时一维、二维模型均采用有限体积法进行求解,一维模型采用HLL格式进行通量计算,具有更好的水流条件适应性,二维模型采用Roe格式可以很好地模拟水面梯度大、水面存在间断等复杂水流流态。采用经典算例对模型进行了验证,验证后的模型对1998年松花江干流的胖头泡溃堤洪水进行了模拟,计算结果可以为此区域的防洪预警提供依据。 1 李曼模型的三维数学模型是基于模型 1.1 维控制策略 一维与沿水深平均的平面二维浅水运动控制方程的统一形式可以表示为 式中D为系数矩阵;▽·F表示变量F对空间变量的偏导数和;S为源项;Ω为控制体范围;t为时间变量;Ts为持续时间。 (1)对于一维情形 式中B为水面宽度,m;Z为河道水位,m;Q为河道断面流量,m3/s;A为过水断面面积,m2;g为重力加速度(取g=9.81 m/s2);J为沿程阻力损失,表达式为J=(n2Q|Q|)/(A2R4/3),R为水力半径,m,n为Manning糙率系数。此形式的一维控制方程在工程上有较多应用,本文考虑了旁侧入流ql,m2/s,将其作为源项,入流为正,出流为负;vx为入流沿水流方向的速度分量。 (2)对于二维情形 1.2 计算界面通量的确定 采用有限体积法,对方程(1)在控制体Ω上进行积分,其中变量定义在控制体中心,通过Gauss定理可以将控制体内的积分转化为边界积分。 (1)一维情形 式中Fi1、Fi2分别为控制体i两侧界面的数值通量;Δxi为控制体i的长度;珔S为源项在单元内的积分值。 采用HLL格式的近似Riemann解计算界面通量,此计算格式结构简易,同时有较好的系统稳定性。源项中采用水面坡度代表压力项的作用,也避免了对底坡项的复杂处理,有利于复杂地形条件下计算的稳定性。 (2)二维情形 式中m为单元边的个数;Δsi为计算单元的面积;Δlij为单元各边长度;Fij为通过第j条边的数值通量;珔S为源项在单元内的积分值。 针对分洪区内洪水流速变化剧烈以及露滩等特点,采用Roe格式的近似Riemann解计算界面通量,此方法具有良好的激波捕捉能力。底坡源项的数值离散采用特征分解方法以实现对复杂地形条件的适应性。 2 低通道二维计算单元 一维、二维数学模型通过溃口位置所在交界面的水力连接条件来实现模型的耦合。对于规模较大的堤防,溃堤洪水的出流状态与宽顶堰流非常接近,溃口流量可采用宽顶堰流公式计算(表1)。 在溃口位置二维计算单元通过多个网格点与一维计算单元进行连接,连接位置采用堰流公式计算的流量Q与口门位置的水位Z2对一维、二维模型的交换数据进行计算和衔接。在具体计算中,一维模型把旁侧入流ql作为源项(式(1)),所提供的流量值Q作为二维模型的计算边界条件,同时二维模型返回网格平均的水位Z2给一维模型来进行下一时段的计算,见图1。 由于一维、二维模型均采用了有限体积法进行计算,变量都定义在网格的中心(一维网格点即断面,见图1),因此两者连接过程中水力要素相互对应,分洪口流量可以实现一维模型与二维模型在