有压管道充水过程数值模拟.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 有压管道充水过程数值模拟 杨敏 李强 李琳 郑杨 天津大学建筑工程学院 天津 摘要 基于明渠非恒定流的运动方程和连续方程 采用 四点隐式差分格式 结合虚拟流量法 对有压管道 充水过渡过程进行了数值仿真模拟 该充水数学模型能模拟管道从无水变为无压流 再过渡到有压流的整个水流 运动过程 其中包括漫流和明满过渡流 将该模型应用于南水北调天津干线有压输水系统充水仿真计算了充水过 程中出现的水力过渡现象 最不利的情况出现在已经充满后的波动阶段 结果可为具体工程的充水操作预测水锤 规避风险和加强防护提供科学的依据 关键词 充水 数值模拟 隐式差分法 明满流 中图分类号 文献标识码 近十几年来 供水和需水矛盾已经严重地制约着经济的发展 社会的进步 为此 我国相继兴建了 一系列大规模 长距离的调水工程 已建和在建的如引滦入津 引黄入晋 南水北调等工程 有压输水方 式由于可以提高输水规模和保证水质而往往被采用 但其水力控制引起的非恒定流非常复杂 特别是在 有压管道频繁检修时进行的充水 操作中 复杂的水力过渡现象还未被人们所认知 由于管道容水体 积巨大 流量的大小与管道承载有着直接的联系 鉴于此 建立有压管道充水的数学模型显得非常必 要 充水过程 近年来修建的长距离调水工程多采用梯级 的有压输水方式 图 特点是在不输水时也能保证 管道处于有压状态 缩短输水响应时间 如昆明掌鸠河调水工程 南水北调中线天津干线 图 有压管道输水单元 充水前 整个管道处于无水状态 水从初始时刻开始流入管道 由于重力等作用逐渐往下游流动 引入一组名词 水流到达指定位置所经过的时间称为流达时间 而水流前沿 即有水与无水的交界面 称 作波前 为了叙述方便 本文把水流波前到达管内某处 使得管道部分有水而另一部分无水的状态称为 漫流 水流到达末端后水面开始上壅 随着水不断流入管道 水位持续涨高 水面接触管顶后 出现明 满过渡流 在管道为正坡的情况下 明满流分界点将随着水位的不断壅高而逐渐往进口处移动 分界点 收稿日期 基金项目 国家自然科学基金 新世纪优秀人才支持计划资助 作者简介 杨敏 男 吉林长春人 教授 主要从事水工水力学研究 到达首端时整个管道都充满水 此后全段处于有压状态 当达到正常运行所需要的水位时 进行开闸 泵启动等一系列水力控制操作 充水过程结束 充水期间 管道内水流将经过以下 种过渡流态的转换 无水状态 漫流 无压流 明满流 有压流 正常运行 数学模型 建立数学模型需要解决两个问题 漫流和明满流的模拟 漫流 当波前未到达末端 整个管道分成明渠流动的有水部分和无水部分 有水部分的明渠流动 可采用明渠控制方程求解 而水深 流量均为零的无水段无法与有水段采用统一的微分方程来描述 为 了模拟漫流阶段水流逐渐往管道末端流动的过程 对无水部分采用虚拟水深 流量法 即认为自波前 开始的无水部分有一股流量 水深均非常小的水流 其流量为入流流量的百分之一至千分之一 具体编 制程序实现时 在每一时间步长的循环迭代过程中跟踪水流波前的到达位置 并记录下来 而自波前开 始的无水段又重新赋予虚拟的水深及流量 即保持无水部分的虚拟流量 水深为恒定值 直至管道内水 流从漫流过渡到无压流 明满流 明渠非恒定流控制方程如下 运动方程 ? ? 连续方程 式中 为流量 为水深 为管道底坡 为水面宽 为过水断面面积 为湿周 为曼宁粗糙系 数 为有压流动的波速 明满流已有较多文献 作了深入研究 本文根据