冲刷漏斗纵剖面形态数值模拟技术研究.pdfVIP

冲刷漏斗纵剖面形态数值模拟技术研究1 1 2 2 1 冯小香 ，张小峰 ，韦直林 ，李建兵 1. 交通部天津水运工程科学研究所，天津 (300456) 2.武汉大学水资源与水电工程科学国家重点实验室，湖北武汉 (430072) E-mail：kathyfxx@163.com 摘 要：提出σ 坐标系垂向分层的优化方法，较好地解决了传统σ 坐标变换用于底孔泄流 数值模拟时存在网格边界变化的问题；采用一维水流模型为立面二维模型提供自由水面边 界，即“动刚盖法”确定自由表面；导出σ 坐标系下床面泥沙边界条件和底孔泄流初期河床滑 塌的判别式。在此基础上建立底孔泄流条件下立面二维水沙数学模型，并将之应用于有刚性 护坦条件下冲刷漏斗形成的数值模拟。模拟结果表明，漏斗冲刷平段长度、漏斗纵坡等数值 试验结果与水槽试验基本一致，验证了模型的合理性和可靠性。 关键词：冲刷漏斗，底孔泄流，河床滑塌，立面二维 中图分类号：TV8 1 引言 在多沙河流上修建水库或引水工程，经常会设置底孔以排水排沙。当底孔闸门开启后， 近孔水流向孔口汇集，孔口周围的淤积物甚易发生冲刷，特别是深水底孔前，形成局部冲刷 坑，形似漏斗，称之为冲刷漏斗。冲刷漏斗的形态和变化过程同底孔前水沙运动特性密切相 关，其形态与规模在一定程度上反映了底孔排沙的效果，同时又反过来影响坝前水沙分布， 进而影响枢纽的正常运行。从这一意义上讲，漏斗形态直接关系到泄水建筑物的正常运行乃 至整个枢纽工程的寿命。 目前对冲刷漏斗形态研究的技术手段主要有类比分析，物理模型和数学模型。由于原型 中漏斗的成因十分复杂，受模型试验技术无法严格相似的限制，物理模型试验成果与实际情 况有较大差异。尽管如此，由于原型观测的漏斗规模较模型规模大，原型漏斗的纵坡和横坡 较模型平缓，模型试验成果偏于安全。对于重要的水库建设，用模型试验研究冲刷漏斗形态， [1] [2] 为工程设计及调度方式提供依据和验证，是十分必要的。三峡水库 、小浪底水库 等大型 水库建设工程均采用物理模型试验研究了冲刷漏斗形态，对模型试验成果进行整理和分析， 极大地丰富了冲刷漏斗形态的研究成果。 相比上述资料分析和模型试验的研究成果，采用数值模拟对冲刷漏斗形态进行研究工作 相对薄弱。曹志先等[3]根据提出的水沙两相流立面二维数学模型，研究了水库坝前冲刷漏斗 平衡形态及相应水、沙流动特性。周新民[4]基于两相流理论和随机理论，提出挟沙水流两相 流湍流模式，对排沙底孔前冲刷漏斗的形态进行了数值模拟。 鉴于数值模拟在多方案比选试验中具有独特的优势，且其试验结果不受模型比尺及观测 设备的影响。本研究将在已建立面二维悬浮物分布数学模型[5] 的基础上，针对底孔前水沙运 动规律的特殊性，对模型功能进行针对性的优化和完善，建立能够反映复杂边界形状、准确 模拟孔口泄流条件下局部水沙运动特性及河床形态变化过程的立面二维数学模型，其研究成 果将具有重要的理论研究意义和工程应用价值。 1本课题得到教育部高校博士点基金项目(水电站厂前水沙运动三维数学模型研究，20030486030) 的资助。 - 1 - 2 σ坐标系下立面二维水流泥沙数学模型的建立 2.1 基本方程 σ 坐标系下立面二维水沙运动方程的通用形式为 ∂ ∂ ∂ ∂ ⎛ ∂Φ ∂⎞ ⎛ΓΦ ∂Φ⎞ ρH Φ + ρHuΦ + ρHωΦ H =Γ + +S (1) ( ) ( ) ( ) ⎜ Φ ⎟ ⎜