水轮机尾水管涡带的模拟论文.pdfVIP

水轮机尾水管涡带的模拟 ThomasAschenbrenner (福伊特西门子水电集团，上海200240) 【摘要】本文研究部分负荷工况尾水管内部的流动。对尾水管的涡带进行预测及由此产生的压力脉动，并与 实测值进行了比较。计算结果在幅值和频率方面均与模型试验数据吻合较好。但是，在模拟中涡带被严重衰 减了，在下游的涡带强度没有完全预测出来。这主要是由所使用的紊流模型造成的。因此。后面讨论了改进 紊流模型并对效果进行了评价。特别需要指出的，本文还考虑了水电站流道系统对尾水管涡带造成的压力脉 动响应。 【关键词】部分负荷涡带数值模拟特征线法紊流模型 1引言 部分负荷工况下的运行是水力机械领域的一大难题，特别是对于混流式水轮机组在部分负荷下可能发 生严重的压力脉动，这是由尾水管内部非稳态漩涡特性(即众所周知的涡带)造成的。通常。涡带以 30％一50％倍转频旋转，形成尾水管内旋转的压力脉动。除了该压力脉动以外，另外一个同步的压力脉动 可能发生在尾水管肘管段内。这意味着在尾水管进口存在一个总的压力水平的波动。该同步部分对整个 水电站系统进行激振，尤其是当激振频率接近压力钢管的特征频率的时候，因此。需要对运行范围进行 必要的限制。 本文的目标是对涡带进行数值预测，并研究涡带对水电站整个流道的影响。所以采用CFD模拟尾水 管的内部流动。因为该模拟要求更加精确的紊流模型，而标准的扣8模型不能预测该现象，所以本文中讨 论了各种不同的紊流模型。 2数值方法 or Maih6fel周。 Stokes方程的各种紊流模型。关于算法的详细信息，读者可参考Ruprechttll 3紊流模型 对非稳态漩涡运动进行模拟比进行稳态模拟计算要求更加复杂的紊流模型。否则漩涡运动会被严 重衰减甚至所预测的流动变成了稳态的。比如采用标准的k-e模型(该模型经常被用于稳态计算)会 好一点的结果。该双尺度模型额外考虑了大涡的耗散率，在8方程中增加了额外的乘积项。改造后的8 方程为 鲁+奶酱2击[(v+vt，I酱卜-。詈Pk-c：睾+竺(譬)·Pk， 采用该模型可以获得比采用扣￡模型较好的结果。图1中显示了直扩散段的流动模拟结果。在进12给 定了有旋的速度边界条件，该速度分布是通过部分负荷下转轮内的流动模拟得到的。从图1中可以看出． 标准的如￡模型获得的仅仅是稳态的流动，然而Kim和Chen提出的扩展的扣￡模型可以预测非稳态的旋 转涡带。 i86 第一届水力发电技术国际会议论文集(第一卷) 图1采用标准扣8模型及KimChen提出的k-8扩展模型预测直扩散段内漩涡流动 4弯肘型尾水管内涡带 弯肘型是水力机械中常用的尾水管型式。本文所研究的尾水管有两个支墩。计算模型的进口转轮泄 水锥处，从转轮泄水锥到尾水管进El之间简化为一直管段。计算网格和进口边界条件见图2，网格有约 250000节点。从三维稳态计算得到的部分负荷下转轮出口的速度的圆周方向平均值作为计算的进口边界 条件。在尾水管出口假定了一个常压。计算采用了Kim和Chen提出的扩展的扣8模型。 图2计算网格及弯肘型尾水管的进I=I边界条件 计算结果显示出尾水管内存在很强的非稳态流动特性。图3显示了两个不同时刻以等压面表示的涡带 形状，可以看出涡带具有旋转的螺旋形的形状。涡带尺寸的改变表示存在很强的压力变化，也就意味着出 现同步的、不旋转的压力脉动。该同步的分项往往造成水电站里面出现的某些问题。 图3两不同时刻尾水管内涡带 为了对非稳态流动模拟结果与试验结果进行对比，所取每一 时间步长处压力值的位置与模型试验中压力脉动测量位置相同， 如图4所示。作为一个例子图5中显示了点l、2处压力的实测值 和计算值沿着时间的变化对比情况。 可以明显地看出。在点1位置模拟得到的压力值与实测值符 合较好。预测的压力脉动幅值比实测值稍小，但是两者的频率非 图4测