反调节水库非恒定流数值模拟.pdfVIP

第12卷第4期 水科学进展 Vol12．No4 2001年12月 INWA‘rER Dec．．200l ADVANCES SCU￡NCE 文章编号：1001—679l{2001)04．0491—08 反调节水库非恒定流数值模拟 魏根群，陈璧宏，宿晓辉 大连理工大学土木系，辽宁 大连116024) 摘要：采用二维水深平均方程组作为控制方程．利用边界拟台法进行坐标变换，采用有限体积 法计算反调节水库的流场。在计算中，由于已知上游泵站抽水流量及下游枢纽供水流量．且上 下游均找不到水位流量关系，上下游边界均为流速边界条件，没有压力边界条件，解不惟一， 所以采用全水库水量平衡方程作为补充方程来补充定解条件。 关键词：反调节水库；边界拟合法；控制体积法；有限分析法；非恒定流 中图分类号：Tv133．2 文献标识码：A 在我国某水利枢纽下游11．0km处，拟修建一低水头水利枢纽，两者之间形成反调节水 库。水库上游建大型泵站抽水蓄能发电，下游枢纽负担向其下游城市供水任务。针对水库上游 抽水和下游枢纽供水的情况，作了反调节水库水流运动的数值计算。在水库上游抽水流量及下 游枢纽供水流量已知的情况下，采用二维水深平均方程组作为控制方程，利用边界拟合法进行 坐标变换，用有限体积法离散控制方程，并利用全水库水量平衡方程作为补充方程以得到惟一 解。 1 控制方程组的坐标变换及边界拟合法 采用一种不同于11Iompson变换…的边界拟合坐标法【2-”，边界拟合法能将天然的不规则域 转化为规则的计算域，使边界条件得到准确满足。通过坐标变换将复杂物理域上的定解问题化 为规则计算域上的定解问题，同时对控制方程也进行变换。将图1的不规则域D变换到如图2 的规则计算域D’。 壁 ” f 围L不规州物理域 图2规则计算域 Fig．1IrregL】1arphysicalplan。 Fig．2R。gI山”cak山8血培pI蚰e 收稿日期：2000-∞一29；修订日期：2000．11一18 作者简介：魏根群(1975一)，男，河南开封人，大连理工大学土建学院在读硬士研究生，主要从事水 工结构及河流数值模拟方面研究。 492 水科学进展 第12卷 如图1所示，在规则的(}，口)平面中建立的变换方程式为 嘉+鲁=。 (1) 耄+骞=。 (1) 与其相对应的边界条件为 耳) (xy)∈n (z=如({， 髓滟；臻j羹；?， l，=，2(}，_) (}_)∈厂，2 日) (z fz=墨(￡， y)∈，14 7 寸) ({目)∈ra 髓豢；臻簇?， 【y=y4(∈， 对方程组(1)，(2)的求解，采用了有限分析法““，其离散表达式为 4 7 3