溢流坝反弧段的动量积分方程组及其数值解.pdf

溢流坝反弧段的动量积分 方程组及其数值解 牛争鸣 李建中 (陕西机械学院) 提要 本文根据水流连续，生『原理与动量原理 ，推导 出考虑 T水深的变化 ，边幂层 的发 展 ，重力及离 力作 ．qq的溢流坝反弧段连续性方程和动量积分方程 ，并用四阶龙格库 塔法对其求解 经验证，边界层厚度与，深 的数值解与实测值基 本一致 ．说 明本文 推导 的积分方 程纽可供 实际应用 主题 词 溢流坝 边界层 动量方程 数值计算 溢流坝反弧段边界层的数值计算是近 几 年才 开始的，已有的方法[1]、[2]、[3]、[4) 都是求饵动量积分方程，只是各方法所考虑的因素不同，方程的形式不同而 已。文 献 [1] 推导了溢流坝反弧段的边界层动量积分方程，考虑了重力殛离心力的影响。值该方程在求解 时，假定流速指数 为常数，水深 沿程不变，切应力 o与陡坡上相同，边界层特忭厚度 的 定 义 与平板边界层相同。有关文献的实测资料表 明，这些假定多与实际情况不符。文 献 [2 利用赫德 (Head)的传输方程计算了溢流反弧段边界层的发展，但计算结果在反弧 中点 后 与实测值差异很大，变化趋势相反。文献 [3]利 用墨菲(Manphy)的曲面边界层动量积分方 程计算的结果也是 如此。文献 韶]、 Ca]都是根据气体在曲面上流动的方程而计算的，这 说 明以往适用于气流的曲面边界层的动量积分方程对有 自由表面，重力与离心办作．用明显， 特别是存在着双向压力梯度的溢流坝反弧段是不适应的。文献 [4]利埚赫兹 (Hinze)的极 坐标形式的雷诺方程，推导了曲面边界层动量积舟方程，但计算时同样取水探q-流球指数为 常数，所暑j压强分布的经验公式也有一定的局 性。由上可见，以上所述的动量积分方程都 未考虑水深、流速指数的变化，计算结果 与实际也不完全相符，陶此溢流坝反弧段的边界层 计算有必要进一步 的研究 1溢流坝反弧段动量积分方程组 以往 的边界层动量方程都是仅取边界层以内的水涤作为控制体推导的，在溢流反弧段， 水深与边界层厚度为同一数量级，边界层 的发艟必然引起水深的变化 ，而水深 的变化 又必然 作为边界层的外部条件影响边界层的发展。因此我们取整个愿皿段水涤内的一个微元段作为 l9 ． 控制体推导连续性方程与砖量积分方程 。溢流坝反弧 段的水流与几何参赞如图1所示 取正交曲线坐标系， 为法向坐标， 为切向坐标，沿流向为 正， 上切点 为坐标原点。 为反弧 半径，O为反弧段圆心角，a为 微元段所在断面与铅垂线的夹角，．正方向为逆时针方 向。 1．1 连续性方程 图l镒流坝反弧段水流与几何参数 设 表示反弧段沿水深的流速分布： ， o≤ ≤d “={ (1) u d《 ≤ 文献 [5] 、 _-)63， [7：指出，溢流反弧段的流速分区性质是很 显的。边界层 内柏 流 速 指数律 ≈=u ( ／d) ／’ (2) 有较大的适应性。边界层以外势流区驹流速分布： U，=U。，八l—y／R)≈u。，(1-F ／R) (3) 是近似成立的。式中U。为边界层外缘的势流流速， 为硅面势流流谜。 根据连续原理，沿任一过水断面对流速积分即为单 宽流量，且沿程不变； 』：a=』：·+u，a= ’ 将 (2) ， (3)式代入积分后： 一i (1一 一手)+h(I+一务)=口