边界条件谏壁泵站的进水流道在立面方向采用了对称的对拼肘形的形式，因此，进水流道的计算只在排涝、灌溉中选一种就可以，选择计算区域时，要包括所有的进水区域。在整个计算域内，只要确定了进口边界条件、固体壁面边界条件就可以计算44中国农村水利水电得到整个流道和出口处的流场。
（ 1）进口边界条件。计算区域的进口为流道的进口（只有一侧） ，出口为泵的进口区域，其他为固体壁面。对于进口，认为来流在垂直方向上为对数分布，同时还要假定进口的湍流动能k和湍流动能耗散率< 3>。
（ 2）固壁边界。在固体壁面附近为湍流的粘性底层，粘性的作用远远超过惯性力的作用，因此采用壁面定律求解近壁处流速。壁面定律就是在靠近壁面的地方，流动主要由粘性作用控制。平均流动参数只与到壁面的距离y、流体密度、粘性和壁面切应力w有关。即：u = f （y， ， ，w）应用量纲分析得到：u + = U u = f （u y） = f （y +）（ 5）式（ 5）就称为壁面定律，其中包括2个重要的无量纲组u +和y +的定义。速度尺度u = （w / ）1/ 2，称为摩擦速度。
计算结果及分析计算对象以及评估流场的目标函数谏壁泵站效率偏低，分析其原因，除选用的2. 8CJ-70泵的运行工况点偏离其高效区以外，流道方面也具有重要的原因。在不是自排（泵运行打水）的情况下，现行双向流道在泵的进出口处形成了不规则的三通，因此损失较大。同时考虑到自排的情况，联系到单侧钟形进水流道的优点，这里尝试在进水流道处加近似钟形流道的导水锥。
进水流道是对称肘形的，在一侧流道进水时，另外一侧为固壁边界。进水流道的设计就是为了使轴流泵获得最佳的进口流场，具体地说，也就是使进水流道计算区域的出口流速分布达到理想状态，据此，引入了2个评价进水流道计算区域出口情况的目标函数：流速分布均匀度V u = < 1 - 1 u a（u ai - u a）2 m > 100%（ 6）式中： u a为计算域出口断面的平均轴向速度； u ai为计算域出口断面各单元的轴向速度； m为流场数值计算时该断面所划分的单元数量。
速度加权平均角度- = u ai < 90 - arctg（u ti u ai） > u ai（ 7）式中： u ti为计算区域出口各单元的横向速度。
可以看出，当速度分布均匀度越接近1、速度加权平均角度越接近90时，流动状态越理想。
计算结果及分析显示截面示意图。由于最优流场的目标是使计算区域的出口（泵进口）的流态达到最优，因此这里列出了接近泵进口的几个断面的流场图， 5.单个流道的设计工况点流量为27 m 3 / h.列出了优化目标的计算结果。
泵站进水流道出水断面目标函数计算结果是否设导水锥最大流速u max / （ m s - 1）最小流速u min / （ m s - 1）平均流速u/ （ m s - 1）均匀度V u / %最大角度max / 最小角度min / 平均角度- / 不设5. 76 4. 58 5. 24 94. 50 89. 98 79. 63 86. 56设5. 82 4. 61 5. 24 94. 84 89. 83 78. 29 86. 15显示截面示意图 X向断面流场图（不设导水锥）由流场图可以看到：在进水流道内，一部分水流平顺地从前方直接进入喇叭管，由于流道宽度较大且未在平面方向收 X向断面流场图（设导水锥） Z 1断面流场图（不设导水锥） Z 1断面流场图（设导水锥）缩，故而一部分水流从两侧进入喇叭管，还有少量水流绕至后张-何法则是分别给出/ P = 1， 2， 3， 4下流量系数曲线图。将张-何法中4张曲线图进行叠加表明，流量关系受/ P影响并非完全按（ / P）/ （ / P + 0. 18）变化的。因此，在实际工程使用中虽然USBR法和张-何法差别不大，但在适用范围方面以及准确程度上张-何法略显优势，更具实用价值。
综上所述，本文认为，在实际工程中，直线形迷宫堰单宫为三角形或0. 8 max梯形时的流量系数计算应推荐使用张志军-何建京法。
迷宫堰流量系数的公式拟合选用张志军-何建京法中m W（本文中的m W0）曲线关系图作为公式拟合的对象。对曲线进行均匀采样取点后，利用商业软件origin进行数据的曲线公式拟合，可以得出对于固定的/ P和H / P值所有曲线均满足：m W 0 = x（1 - e - y L/ W）式中： x和y为待定的参数。
拟合曲线与采样数据点的相关系数R > 0. 99，说明拟合的精确程度较高，是可信的。
由因次分析可知，待定参数x、y应同时与H / P和/ P有关。为简化拟合的维数，仍将/ P作为已知参数。这样，对应于一个/ P，参数x、y就只有与H / P有关。当L/ W 1时，迷宫堰将趋近于传统直线堰，流量系数m W 0就将逼近直线堰时的流量系数表达形式。这里选取！经验公式m = 0. 402 + 0. 054（H / P）作为参考直线堰流量系数公式，由此可得y与x的函数关系式：y = - ln< 1 - 1 x（0. 402+ 0. 054 H P） >再通过寻找x与对应H / P的关系并整理，可最终得到流量系数m W0关于L/ W和H / P的表达式：m W 0 = x< 1 - e L W ln（1- 0. 41 x）> x = 0. 9 + 4. 2 e - H 0. 17 P（P = 1）（ 7）m W 0 = x< 1 - e L W ln（1- 0. 42 x）> x = 0. 85+ 4 e - H 0. 29 P（P = 2）（ 8）m W0 = x< 1 - e L W ln（ 1- 0. 45 x）> x = 0. 82 + 3. 7 e - H 0. 21 P（P = 3）（ 9）m W 0 = x< 1 - e L W ln（1- 0. 45 x）> x = 0. 8 + 3. 7 e - H 0. 34 P（P = 4）（10）上述公式的拟合范围是： 2 L/ W 8， 0. 1 H/ P 0. 9.
为验证拟合公式的准确程度，对原始曲线上的采样点得出的流量系数值与拟合公式计算的计算值进行回归统计分析，结果如下：式（7）U = 6. 864 1， Q = 0. 483 1， R = 0. 966 6；式（8）U = 8. 866 2， Q= 0. 430 7， R = 0. 976 6；式（9）U = 9. 071 8， Q = 0. 542 4， R = 0. 971 4；式（10）U = 9. 434 4， Q = 0. 590 5， R = 0. 970 1.
U为回归平方和， Q为残差， R为相关系数。由此可见，拟合公式计算的流量系数m W0值与张-何法查图所得值非常接近，实际应用中可以用上述拟合公式来计算流量系数值。
|
