含生产堤等线状建筑物蓄洪区水流数值模拟.pdfVIP

2007年第5期 ·PEARLRIVER 人 民珠江 含生产堤等线状建筑物的 蓄洪区水流数值模拟 胡晓张 ，张小峰 ，何 用 ，陈文龙 (1．珠江水利科学研究院，广东广州 510611；2．武汉大学水资源与水电工程科学国家重点实验室，湖北 武昌 430072) 摘 要 ：在蓄洪区水流数值模拟中，对于含有生产堤和公路等线状建筑物的处理方法存有一定问题，根据水流在生 产堤 内外运动的不同特点，考虑生产隔堤等阻水建筑物的影响，在网格设计和模型方法上进行改进，建立适合模拟 蓄洪区内存在生产隔堤等线状建筑物时水流运动的数学模型。 关键词：蓄洪区；生产堤 ；线状建筑物；水流数学模型 中图分类号：TV32．3 文献标识码 ：B 文章编号：1001-9235(2007)05-0021-05 0 前言 方向运动方程：一 = (2) 蓄洪区内的地形大都比较复杂，有时蓄洪区内还会存在 生产隔堤和公路等线状建筑物 J，这类建筑物会对蓄洪区内 y方向运动方程：一 = (3) 的水流运动产生阻水作用 ，蓄洪区内的水流运动受到生产隔 堤等的影响将与不存在生产隔堤等建筑物时有很大的区别， 式中 —— 隔堤糙率； —— 隔堤上水流在 方向上的 因此在模型计算时必须考虑生产隔堤等阻水建筑物的影响。 速度；Pdy—— 隔堤上水流在 Y方 向上 的速度； 目前蓄洪区的水流运动过程主要采取笛卡尔坐标系下 — — 隔堤上水深；A—— 隔堤的平面面积；Q—— 的平面二维水流数学模型进行模拟计算 。当蓄洪区内存 进入隔堤内的流量； —— 隔堤的地面高程。 在生产隔堤等线状建筑物时，大多采用把生产隔堤等建筑物 1．1．2 流量公式 的高程平摊到网格单元上 ，然而由于生产隔堤等建筑物的横 隔堤布置在网格线上 ，如图 1所示，隔堤两侧流入和流 向尺寸相对于模型计算所用的网格尺寸非常小，而且水流在 出的流量采用堰流公式计算。 隔堤内外的运动规律也不相同，这种处理方法导致计算结果 如图2所示，h，为生产堤左侧网格内相对堤防水深，h 存在较大的误差。如果不针对这种情况进行处理，并改善平 为生产堤右侧网格内相对堤防水深，h为生产堤 内水深，h， 面二维水流数学模型，模型计算结果的精度难以得到保证。 为生产堤左侧 网格水深，h。为生产堤右侧 网格水深，z。、z 因此，首先在设计模型计算网格时应尽量使生产隔堤为所建 为生产堤左右两侧高程。UMU和 UMD为生产堤两侧进出的 计算网格的边，其次根据水流在生产隔堤内外运动的不同特 单宽流量。 点，对已有的蓄洪区平面二维水流数学模型进行研究和改 以生产堤东西走向为例，计算生产堤两侧的进出流量如下。 进 。 a)UMU计算公式 1 模型介绍 ① 当hl≥hd时：若2hl≥3hd，则UMU= lhl gl； 对于蓄洪区内生产隔堤外的区域 ，计算方法已非常成 若2hl3hd，则UMU= 2hd~／2g(hl—hd)； 熟 ，按照已有的蓄洪区二维水流数学模型进行计算即可，这 (当hlh时：若2hd≥3hl，贝0UMU=一 lhd√gd； 里不再赘述。重点需分析生产隔堤上 的水流运动，对于生产 若2hd3hl，则UMU=一 2hl~／2g(hd—h1)。 隔堤