长江弯曲型分汊河道分流区水流泥沙运动特性.pdf

三、泥沙运动基本理论及泥沙输移 ‘325’ 长江弯曲型分汊河道分流区水流泥沙运动特I生 姚仕明12余文畴2 (1．清华大学水利系 北京 100084；2．长江科学院 武汉430010) 摘要：长江中下游分汉河道分流区水流泥沙运动较为复杂。本文运用数值模拟的手段计算分析了弯曲 型分汊河道分流区的水位、流速及挟沙因子的变化特点。结果表明，在不同流量情况下，分汉河道分流 区主汉侧的水位低于支汉侧，存在横比降，主汉侧流速及挟沙因子大于支汉侧；分流区横向水位差值随 流量增加而增大，流速与挟沙因子也随流量的增加而增大；横向水位沿程变化差值增大，流速与挟沙因 子沿程变化则由大流量情况下的减小到小流量情况下的增大。 关键词：分汊河道 分流区 数值模拟水沙运动 l前言 分汊河道分流区的水沙运动较为复杂。正确认识其水沙运动规律具有重要的意义。 目前，国内外对分汊河道分流区水沙运动特性进行系统研究的尚不多，研究方法和手段 基本为概化模型试验与原型资料分析。由于概化模型与天然分汉河道的实际地形存在差 异，且分汊河道原型观测资料也较有限，给深人认识天然河流分汊河道的水沙运动特性 带来困难。因此，结合天然分汊河道的地形特点，运用数值模拟的手段来研究分汊河道 分流区的水沙运动规律是非常必要的。本文运用二维紊流数学模型，对弯曲型分汊河道 分流区水流运动及挟沙因子分布特性进行了研究，为进一步认识分汊河道特性与有效治 理分汊河道奠定了基础。 2数学模型的建立 2．1 正交曲线网格的生成 天然河道往往具有复杂的几何边界，一般难以直接用直角坐标系下的水流运动方程 去求解。通常采用坐标变换的方法将天然河道较为复杂的边界变换为贴体正交曲线网格， 便于适应复杂天然河道边界条件的水流计算，本文就是采用贴体正交曲线网格，其中网 格横向节点数为50，纵向为150，如图1(a)所示。 2．2水流控制方程的离散 在正交曲线坐标系下，水流控制的通用方程可变为： 堕％笋+堕铲=善(LH苦簧]+南[LH旦C．塑Oft)1+ ㈣a善 a玎 a善L‘ G a毒／a，7 L。 、’ 作者简介：姚仕明，男，1974年出生，博士研究生。 ．326． 第六届全国泥沙基本理论研究学术讨论会论文集 式中：“、v分别为正交曲线坐标系下毒、，7方向的速度分量；H为水深；几为扩散系数； 毋为待求变量；为广义源项，由三部分组成，包括方程标准化后对流项剩余项、扩散 项剩余项与笛卡尔坐标下的源项，通过方程负坡线性化后可表示为=Sc+鄙咖。不同变 量的求解差别主要体现在上。 0．395E+06 芒 k 530000 535 000 540000 545 000 550000 x(m) (a) (b) 图1 计算河段网格布置与三维地形示意图 本文计算水流所采用的紊流模型为标准的拈占模型，采用控制体积法对式(1)进行离 散，不同方程的求解变量按交错网格模式布置，即在交错网格模式中，压力、日与k、占布 置在主节点上，而U、v分别位于善、枋向的界面节点上。在计算域(善、啪中，令△乒 A rl=l，因此一般性控制方程的离散形式可表示为 aP西P-E口6西H6+ap‘aI)p‘+6 (2) 式中：郇=∑％+西一品△弘叩；口；=旦旦兰At塑；6=△丛叩；4柏为对应节点位置的值。 3水位一流速的耦合求解与定解条件 3．1 水位一流速的耦合求解 求解动量方程的难点在于未知的压力(水位)场，而压力场又是连续方程间接给定的， 也就是对水位一流速进行耦合求解。本文采用压力校正法，即先假设一预测水位值矿， 由动量方程求解出与之相应的而不一定满足连续性方程的速度U+、v+。为了能同时满足 水流连续方程与动量方程，需要通过连续方程修正压力。由连续方程在△豁玎控制体上离 散并整理得： aeZe=芝：tin6乙6+b ap=口E+aw+aN+口c 其中，f71￡= ．g．．．(．．C。．．．c．．H．．．．A．．．．善．．．)．—2—． ％ ， g(CcHA善)： as 2———————-。 ％ 三、泥沙运动基本理论及泥沙输移 ·327· 6_[(q吼+W--(C，Hu+)。]△刁+[(Q胁’3m(Q协+)。]△孝+竺墅掣 以上数学模型控制方程及定解条件的数值离散和求解基于有限体积法和SIMPLEC 算法。在求解离散方程时采用Patanker和Spalding给出的欠松弛技术、块校正技术及逐 行迭代法。 3．2定解条件 在计算区域进口断面处给定计算变量(“，v，k，占)的分布；出口断面则给定计算变量 U，v，k，占)的沿程梯度善=o，并给定控制水位Z=ZOUT；边界为无滑动边界：“=v=o， 娑：害：0。在七一占模型中，固壁条件要特殊处理，本文给定_Ok：0，占在临近岸的节 ／1j／qLj／Z点上给定：Cp=兰兰，修正岸边壁紊动黏性系数”。 t,：y