天然河道淹没丁坝群水深平均平面二维数学模型的研究.pdfVIP

天然河道淹没丁坝群水深平均平面二维数学模型研究 皋国斌韩信 (南京水利科学研究院河港所j 摘要 本文根据平面二维水深平均水流运动方程组，利用控制体积法离散方程并运用sm舻IJ披程式求 解，研究丁天然i可道淹没丁坝群水流汁算平面二维数学模型，针对正交丁坝和菲正交丁坝提瞄了各自的处 理办法，同时分析了计算河道进口流速分布规律。此数学模型经水槽试验资料和长江八卦洲河段物理模型 实测试验资料验证，计算值和实测值符合怠好。 关键词淹没正交丁坝非正交丁坝二蛙数模控制体积法 l前言 丁坝是一种最常用的河道整治建筑物，其作用主要是整治桔水航道及河岸防护，因此 其坝顶高程通常按整冶水位确定，在汛期常处于淹没状态。丁坝根耀其轴线与岸线的夹角不 同，可分为上挑丁坝、正交丁坝和下挑丁坝。丁坝挑角不嗣及淹没程度不同，其对水流的影 响不同。乎面二维数学模型作为一种研究手段已应用予许多河段的水流及泥沙淤积计算，数 模应用中经常会遇到淹没丁坝群的水流计算问题，为精确求解此类问题，考虑淹没丁坝对水 流的影响是十分必要且有重要意义的。 2 平面二维水流计算时前入对淹没T坝的考虑及本文所作的工作 前人对天然河道淹没丁坝群的汁算．冒内外研究不多。李浩麟”’在用单元积分法计算河 口航道潮流中，将淹没丁坝过流视为堰流，使用堰流公式和二维水流运动方程组联台求解， 结果较满意，但其中仅有四条丁坝并且相距较远．因此使用堰流公式是成功的。塥建军口}在 三峡水库回水变动区九龙坡河段具有淹没丁坝的流场计算中，也采用了堰流公式，并与边界 对称点法相配合，进行了平面二维流场计算．取得r较好的效果，但其中仅有一条丁坝。然 而堰流公式对于淹没T坝群水流计算的适用性是需要进一步考虑的，众所用知，水力学中关 于堰流及分类有严格的定义，如根据堰上水头及堰顶宽度之比值丈小，可将堰分为薄壁堰、 实用堰及宽顶堰，堰与堰之间堰流系数差别很大，水流流经淹没丁坝时的流态与流过堰时的 流态尽管稍有相似之处，但毕竟有很大不同。退丽言之，假设丁坝过流可视为与堰流褶同， 但对于天然河道中拥有众多丁坝的丁坝群，堰流公式中流量系数的规律性就难以确定，其大 小与相对坝长、相对坝高、丁坝间距等一系列因素有关，尤其是当河道中丁坝个数众多、丁 坝之间互相影响时，各坝的堰流系数每每不同，数值的确定更为不易，而此数的正确与否直 接影响计算的好坏，因此以堰流公式反映淹没丁坝的作用将在地形变化复杂、丁坝众多且互 相影响之地区推广应用发生一定的困难。李国斌、李昌华c3I对天然河道淹没丁坝群水流计算 进行了研究和探讨，提出了淹没丁坝群水流计算平面二维流带模型，此模型基于绘制水流平 面图的基本思想，将天然河道水流划分为一定数量的流带，在每条流带内用差分法求解水流 运动方程式和连续方程式，在运动方程式中首次考虑了淹没T坝的局部水头损失，以概括 107 淹没丁坝对水流的影响，同时为进行数值计算提出了一系列处理办法，此方法经湘江下摄司 河段及长江八卦洲河段的实测模型试验资料验证，计算值和实测值符合良好，但此法无法计 算丁坝回流及考虑流线弯曲。Shjm】ayafq等认为淹没丁坝在水流中所受的作用力与流体力 学中一个机冀上所受的作用力类似，丁坝受到水流的推移力及上举力，推移力与水流方向一 致，上举力与水流方向垂直。根据牛顿第三定律，作用力与反作用力大小相等、方向相反． 因此相应地，水体受到丁坝的推移力及上举力，将推移力和上举力在运动方向上分解并添加 到水流运动方程组，求解此方程组即可求得考虑淹没丁坝作用的水流运动平面二维数值解， ：{：I●|：|^j●．i 但此方法同样存在推移力系数和升力系数的选取问题，而且这种方法无法考虑水流流过淹没 丁坝时的细部流态，例如丁坝坝顶流速不易计算等，这种方法相当于增加局部糙率，因而也 就只能宏观地研究丁坝的影响。 综上所述，目前对平面二维水流计算时如何考虑淹没丁坝其影响的种种方法尚不够完 善，有待进一步研究。有鉴于此，本文针对正交和非正交淹没丁坝提出了各自的处理办法， 分析了对河道水流计算有较大影响的进口流速分布规律，并将二维数模计算结果与实测资料 进行了比较。 3数学模型 3．1基本方程组 对于平面二维水深平均水流运动’其守恒型基本方程组如下： 丝+堂堕+塑尘：o