旋流式竖井泄洪洞水力要素的数值模拟与试验研究.pdfVIP

水力发电 2008年3月 为s方程常数，C =1。44，C =1。92；G为由平均流速梯度引起 2 计算结果与试验结果的比较分析 的紊动能产生项； 为紊流粘性系数，它可由紊动能k和紊 动耗散率s求出 2．1 上平段及下平段自由水面及涡室、竖井的空腔 采用水气分层两相流VOF方法的标准 双方程紊流 } (5) 数值模型能够很好地模拟旋流式竖井的上平段、涡室、竖井、 G (等+等) (6) 下平段等部位的水面线。模型试验中对上平段、涡室、竖井、 下平段．采用在间隔一定的距离布置测点．对试验数据和计 其中， 为经验常数，取C 0．09。 算值进行对比分析f见图2)。 弓【入水气分层两相流的VOF模型后，在控制体内对第q 相流体的容积分数规定为： =0表示控制体内无q相流体； ％：1表示控制体内充满q相流体；0％1表示控制体内有 部分q相流体，对所有流体的容积分数总和为 1，即 ％=1。 在水气分层两相流中， 为水的体积分数， 。为气的体积分 数，aw+a~1，P和 就是体积分数的函数，而不是一个常数。 它们可表示为 p=awpz+(1一 P (7) = k+(1—0 ) (8) 式中， 和 分别是水和气的密度； 和 分别是水和气的 分子粘性系数，通过对水的体积分数 的迭代求解。P和 都可以由式(7)、(8)求出。 1．2 计算区域、网格划分及边界条件 数学模型模拟范围从桩号0+000～0+281。5．包括旋流式 竖井泄洪洞的上平段、涡室、竖井、压坡段及下平段。对网格 的处理随体形和边界条件的复杂程度而异。一般网格越细 就越能接近客观实际情况，但占用的内存资源就越大，同 时计算所需时间就越长。所以，为了节省计算机资源，缩短 计算时间．空间网格剖分使用了非均匀网格，上平段、竖井 图2 水气两相流分布及上下游明流段水面线模拟与实测值 及下平段的网格较为稀疏而涡室内的小挑坎附近、竖井与 由图2可知：在校核流量下，设计体形仍有较大的洞顶 压坡处网格较密，网格尺度的变化范围为O．05～4．50 m。 余幅．特别是在较为敏感的涡室入口处，泄水洞的宽度变小 导致水深有所增加．但仍没有发生“水流撞顶”的现象，设计 数值模拟采用了控制体积法隐格式迭代求解，即在每个 控制体积上对基本方程进行积分得到离散方程．由于研究对 满足要求．确保了竖井旋流泄洪洞的通气量。由模拟结果和 象是带重力加速度的vOB两相流．用 PISO算法来解决速度 实测值的比较来看：从进口到涡室入口，水面线的模拟值和 场和压力场的耦合问题，效果较好。 实测值吻合良好．在涡室内部，实测水面线存在较大的波动， 整个计算区域如图1所示，水流入口采用速度入口边界 但取平均值后与模拟值符合良好 这验证了计算结果的正确 条件．在库区水位2 857．8 m时，根据模型试验得到引水道入 性 采用该方法得到的水面线分布特点与参考文献[6]相一 口断面水深及流量(O=521．o m3／s)，由此可知入口法向速度 致 下平段的水面线模拟值和实测值与上平段一致。 1O．49 m／s，出口及所有的气体边界均采用压力边界条件， 从图3可以看出：在涡室及竖井内，漩涡空腔基本位于 压力为大气压值。壁面边界条件采用无滑移边界条件，粘性