二维异重流数学模型的研究.pdfVIP

中国水利学会．7．003学术年会论文集 二维异重流数学模型研究 祁伟黄永健毛蛀新 (中国水利水电科学研究皖) 摘 要 本文主要研究因含沙量的差异而形成的发生在河道或水库底部的浑水异重流。对 它的运动规律进行研究．在指导水库异重游捧沙和碱轻航道淤积等方面都具有重要的理论意 义和应用价值。本文针对异重流运动的数值模拟问题，提出了一套二维异重流上下两层水体 的双层计算模式下的数学模型基本方程。同时在对两层水体基本方程简化的基础上。运用控 制体积法对方程进行了离散．并建立了二雏异重漉教学模型。最后运用建立舯二维异重流数 学模型对黄河流壤刘家峡水库中的浑水异重流运动和异重流捧沙进行了二雏数值模拟．结果 令人较为满意。 关■词 异重流二维数学模型双层计算模式刘家峡水库 一、前言 异重流是指两种密度相差不大、可以相混的流体，因为密度的差异而发生的一种相对运动。异重流按 照它的形成方式可分为温度异重流、盐水异重流和浑水异重流等。本论文中主要研究因含沙量的差异而 形成的这种在河道或水库底部的浑水异重流。 浑水异重流是一种经常发生于水库、引航道和河口的特殊水流泥沙运动形式，容易引起泥沙淤积，造 成库容减少和航道维护困难等问题。对异重流的运动规律进行研究和数值模拟，在揖导水库异重流排沙 和减轻航道淤积中具有重要的理论意义和应用价值。 二、模型的建立 由于浑水异重流在接近稳定的状态下，含沙量沿垂线分布曲线在异重流交界面附近有明显的拐点，在 拐点以上含沙量一般较小。可以近似于清水；而 在拐点以下含沙量则较大，且一般沿垂线分布 较为均匀。因此从异重流含沙量分布的这一特 性出发，本文中以异重流的交界面作为分界线， 把水体分为上、下两层(见图1)．然后从三维的 水流基本方程出发，在考虑界面上的掺混与交 换的情况下，把上、下两层水体分别沿垂线进行 积分，得到了双层模式下的二维水流基本方程， 分别称之为上层水流基本方程和下层异重流基 本方程。得到了双层模式下的这两层基本方 图1异重流分层运动示意图 二维异重流教学模型研究 程，理论上可以把它们进行耦合求解，但在实际操作与程序实践中，这种双层耦合求解方法却存在着诸多 的困难与缺点，主要表现在：①上下两层之问相互影响，导致模型极其敏感，误差容易扩大；②模型的数值 收敛性难以控制，稳定性不好；③计算程序庞大而复杂，难以维护。正是由于双层耦合求解模型在实际中 不便于实现，所以将双层模式下的上、下两层基本方程进行了简化，即假设水体的自由水面在达到稳定后 是一个不随时间变化的恒定值，简化结果为以下两个子模型，分别称为二维水流子模型和二维异重流子模 型。 二维水流子模型方程组 鲁+警+盟as=o (1) 警+未(譬)+砉(％})+一鲁+寿Q，师一他=0。(z) (s) 警‘+圭(马乎)．+昙(譬)+曲舅；南Q，4-eJ-cf+恕=o 式中Q。Q，在这里指在整个水体(上层水体、下层水体合在一起)中沿*方向和y方向的单宽流量，h 指整个水体的水深，互为水体自由水面的水面高程，c为谢才系数，g为重力加速度，为科氏力系数。’ 二维水流子模型采用“交错方向隐格式”的数值算法(ADI法)，动边界采用“窄缝法”进行处理。它的 二维异重流子模型的方程组之中，使二维异重流子模型方程组封闭‘。 二维异重流子模型方程组 誓+豢+誓=一咋 ∽’