专利名称：一种突缩突扩式消能工流态的判别方法
技术领域：
本发明涉及水利水电工程泄洪洞中一种突缩突扩式消能工流态的判别方法，主要 用于判别这种消能工在不同结构体型和不同来流条件情况下，消能工内呈现的是洞塞流流 态，还是孔板流流态。
背景技术：
在高坝工程中，为了消杀库区下泄水流的巨大能量，实现库区和下游河道上下游 水流的衔接，泄水建筑物需要采取有效的消能工。传统的消能方式有底流消能，挑流消能 和面流(或戽流)消能等，这些消能方式在水利水电工程的建设和运行中发挥了巨大的作 用。随着工程建设规模的发展，同时考虑到传统形式的消能工在工程运用中的特点和存在 的局限性，如底流消能在下游需要有较长的护坦，土石方开挖和混凝土方量较大，工程造价 较高；挑流消能常常会产生强烈的雾化，影响工程安全和下游环境；面流(或戽流)消能在 下游较长距离范围内引起河道水流的不平稳，影响电站运行和下游航运。近年来，由临时 施工导流洞改建成永久泄洪洞作为一种有效的技术方案为多个工程采用，施工导流洞的利 用，大大提高了工程的经济效益，更便于合理地进行枢纽布置。在泄洪洞内采用突缩突扩式 消能工，结构简单、施工简便、消能率高，有较好的工程应用前景。目前，一般根据突缩突扩式消能工内流态的不同，将其分为洞塞式消能工和孔板 式消能工，当突缩突扩式消能工厚度较薄时，水流进入消能工突然收缩，而后直接进入下游 水体，形成这种流态的突缩突扩式消能工通常称为孔板消能工(见图1)。当突缩突扩式消 能工厚度较厚时，水流进入消能工，先在孔内收缩，然后附着孔的边壁流动，再进入下游水 体，形成这种流态的突缩突扩式消能工称为洞塞消能工。(见图2)。就突缩突扩式消能工而 言，有两个重要的特征参数，其一是孔径比0，定义为消能工孔的直径d与泄洪洞直径D的 比值，即，0 = d/D ；其二是无量纲厚度a，定义为消能工厚度T与泄洪洞直径D的比值，即， a =T/D(见下图3)。发明人发现，对于某种厚度的突缩突扩式消能工，随着工作条件的变 化，如运行的工作水头，或流速，或雷诺数(Reynolds number, Re)，当这些参数较小时，呈现 洞塞流流态；而当这些参数较大时可能出现孔板流流态，因此，对某种厚度的突缩突扩式消 能工，无论是称其为孔板消能工还是洞塞消能工，都不尽合理，因为，随着工作条件的变化， 流态将发生变化，因此，发明人认为，称这种消能工在某种工作条件下是孔板流流态或是洞 塞流流态，是更合理的。在这样的情况下，如何判别突缩突扩式消能工在某种工作条件下， 属于何种流态，即，洞塞流流态或是孔板流流态，是重要的，目前，缺乏行之有效的突缩突扩 式消能工流态的判别方法。图4和图5分别是在3 = 0. 5，a = 0. 6和3 = 0. 7，a = 0. 5的两种体型结构 情况下，当来流的雷诺数Re变化时，从洞塞流流态到孔板流流态的变化过程，其中，图4a和 图4b、图5a和图5b是洞塞流流态，图4d和图4e、图5d和图5e是孔板流流态，而图4c和图 5c则是两种流态过渡的一个临界状态。由此可见，单一从突缩突扩式消能工结构体型来区 分是孔板消能工，还是洞塞消能工是不确切的，需要综合考虑消能工结构体型和工作条件，即来流的水力参数。
发明内容
本发明在综合考虑突缩突扩式消能工的体型结构和来流条件基础上，提供一种对 这类消能工流态判别的方法。本发明采用如下技术方案—种突缩突扩式消能工流态的判别方法步骤1获取消能工孔的直径d、泄洪洞直径D以及消能工厚度T，并计算孔径比3 及无量纲厚度a 孔径比3为消能工孔的直径d与泄洪洞直径D的比值，即3 =d/D；无量纲厚度a为消能工厚度T与泄洪洞直径D的比值，即a = T/D ；同时，获取水力参数雷诺数Re，S卩Re = vD/(ii/p)，v为来流的速度，y为水的 动力粘度，P为水的密度。步骤2计算临界厚度a。a c = (3. 33 3 4_8. 01 3 3+7. 02 3 2_2. 68 3+0. 39) X l(T6Re_46. 42 3 4+l 10. 52 3 3-97. 20 3 2+37. 413—4. 94步骤3将临界厚度a。与无量纲厚度a进行比较当a > a。时，为洞塞流流态； 当a = a。时为临界流流态；当a < a。时为孔板流流态。本发明的优点是提出的临界厚度a。，概念清楚。将影响流态的体型结构参数和来流水力参数，与 临界厚度直接建立函数关系系式，流态判别方法思路清晰，简单易行，准确可靠。
图1是孔板消能工示意图。图2是洞塞消能工示意图。图3是消能工参数示意图。图4是体型结构参数为3 = 0.5, a = 0. 6时，消能工流态演变状况，其中，4a和 4b分别是雷诺数Re为0. 1 X 106和17. 53X 106时，消能工呈洞塞流流态，4c是雷诺数Re 为26. 29X 106时，消能工呈洞塞流流态到孔板流态演变过程，4d和4e分别是雷诺数Re为 35. 06X 106和52. 59X 106时，消能工呈孔板流流态。图5是体型结构参数为日=0.7，a = 0. 5时，消能工流态演变状况，其中，5a 和5b分别是雷诺数Re为0. 1 X 106和4. 38X 106时，消能工呈洞塞流流态，5c是雷诺数Re 为16. 65 X 106时，消能工呈洞塞流流态到孔板流态演变过程，5d和5e分别是雷诺数Re为 58. 72X 106和70. 11 X 106时，消能工呈孔板流流态。图6是本发明的方法流程图。图7是流态判别曲线图。
具体实施例方式一种突缩突扩式消能工流态的判别方法步骤1获取消能工孔的直径d、泄洪洞直径D以及消能工厚度T，并计算孔径比3 及无量纲厚度a 即，0 = d/D, a = T/D，同时，获取水力参数雷诺数Re，Re = vD/(u/ P)，v为来流的速度，P为水的动力粘度，P为水的密度。步骤2计算临界厚度a。，即a c = (3. 33 3 4_8. 01 3 3+7. 02 3 2_2. 68 3+0. 39) X l(T6Re-46. 42 3 4+l 10. 52 3 3
-97. 20 3 2+37. 413-4. 94步骤3将临界厚度a。与无量纲厚度a进行比较当a > a。时，为洞塞流流态； 当a = a。时为临界流流态；当a < a。时为孔板流流态。本发明对突缩突扩式消能工，提出了临界厚径a。概念和计算方法，通过消能工厚 度a与临界厚度的比较，判别突缩突扩式消能工运行中的流态，是洞塞流流态还是孔板流 流态。临界厚径a。定义为突缩突扩式消能工在工作条件下，水流流动发生从洞塞流流 态过渡到孔板流流态，或从孔板流流态过渡到洞塞流流态的一种临界流动状态时，消能工 所需要的厚度，即，ac = Tc/D0图7是不同孔径比0的突缩突扩式消能工结构，雷诺数Re与临界厚度a。的关 系曲线图。由图可以看出在各个不同的孔径比之下，临界厚度与雷诺数基本成线性关系。 对于给定的消能工厚度a和消能工的孔径比3，在确定其运行的雷诺数Re，查得若此消能 工厚度a在相应的曲线之上，则其流态为洞塞流流态，反之，在曲线之下，流态为孔板流流 态。由此，可以方便地判别突缩突扩式消能工的工作流态。发明人将这种曲线关系拟合成 一个临界厚度a。关系表达式a c = (3. 33 3 4_8. 01 3 3+7. 02 3 2_2. 68 3+0. 39) X l(T6Re-46. 42 3 4+l 10. 52 3 3-97. 20 ^2+37. 41^-4. 94当a > a。时，为洞塞流流态；当a = a。时为临界流流态；当a < a。时为孔 板流流态。
权利要求
一种突缩突扩式消能工流态的判别方法，其特征在于，步骤1获取消能工孔的直径d、泄洪洞直径D以及消能工厚度T，并计算孔径比β及无量纲厚度α孔径比β为消能工孔的直径d与泄洪洞直径D的比值，即β＝d/D；无量纲厚度α为消能工厚度T与泄洪洞直径D的比值，即α＝T/D；同时，获取水力参数雷诺数Re，即Re＝vD/(μ/ρ)，v为来流的速度，μ为水的动力粘度，ρ为水的密度。步骤2计算临界厚度αcαc＝(3.33β4-8.01β3+7.02β2-2.68β +0.39)×10-6Re-46.42β4+110.52β3 -97.20β2+37.41β-4.94步骤3将临界厚度αc与无量纲厚度α进行比较当α＞αc时，为洞塞流流态；当α＝αc时为临界流流态；当α＜αc时为孔板流流态。
全文摘要
一种突缩突扩式消能工流态的判别方法，获取消能工孔的直径d、泄洪洞直径D以及消能工厚度T，并计算孔径比β及无量纲厚度α孔径比β为消能工孔的直径d与泄洪洞直径D的比值，即β＝d/D；无量纲厚度α为消能工厚度T与泄洪洞直径D的比值，即α＝T/S；获取水力参数雷诺数Re，即Re＝vD/(μ/ρ)，v为来流的速度，μ为水的动力粘度，ρ为水的密度。计算临界厚度αcαc＝(3.33β4-8.01β3+7.02β2-2.68β +0.39)×10-6Re-46.42β4+110.52β3-97.20β2+37.41β-4.94。将临界厚度αc与无量纲厚度α进行比较当α＞αc时，为洞塞流流态；当α＝αc时为临界流流态；当α＜αc时为孔板流流态。
文档编号E02B8/06GK101858069SQ20101020781
公开日2010年10月13日 申请日期2010年6月23日 优先权日2010年6月23日
发明者吴建华, 周琦, 王静清, 艾万政, 苗沐露, 许唯临, 邱海娟 申请人:河海大学;四川大学