《Water》：Energy Dissipation Mechanism of Draft Tube Vortex Rope Suppression by Self-Driven Jets from the Discharge Cone
编辑推荐：
本研究考察了水轮机尾水管在部分负荷工况下轴向水注入所引起的解析流结构重分布及基于速度梯度的熵产。采用稳态SST k–ω计算对6种射流半径（4~9 mm）进行筛选，随后对无射流工况与选定的7 mm工况开展非定常嵌入大涡模拟（ELES）。7 mm半径对应于锥体半径
本研究考察了水轮机尾水管在部分负荷工况下轴向水注入所引起的解析流结构重分布及基于速度梯度的熵产。采用稳态SST k–ω计算对6种射流半径（4~9 mm）进行筛选，随后对无射流工况与选定的7 mm工况开展非定常嵌入大涡模拟（ELES）。7 mm半径对应于锥体半径的0.175倍及名义射流流量份额1.6%，该候选方案系基于稳态流动拓扑筛选而非已证实的熱力学最优值。在所报道的数值设置下，选定射流消除了稳态计算中的中心回流区，并在非定常可视化中破坏了相干螺旋结构。所报道的1.6 Hz压力幅值在DT1处由249 Pa降至167 Pa，在DT2处由110 Pa降至71 Pa，降幅分别约为33%和36%。对于假设的10%旁通水头损失，射流动能通量估算为水轮机水力输入功率的1.44%；该值并非实测净效率惩罚。速度梯度分析表明解析耗散位置发生改变，但未包含独立的亚格子尺度贡献及壁面函数熵产计算。因此，结果支持对该控制机制的条件性解释，而非完整的损失平衡或已验证的自驱动系统设计。
研究背景方面，部分负荷工况下水轮机转轮出口残余旋流与尾水管中心轴向速度亏损相互作用，径向平衡要求压力沿径向向外增加（近似正比于ρuθ2/r），从而形成低压中心区。扩张尾水管中的减速可增强逆压梯度并促进中心回流或涡带破裂，低动量核心与外围旋流之间的剪切层可支撑螺旋全局不稳定性，其进动偏心核心产生周向非均匀壁面压力脉动。进动频率通常报道为转轮旋转频率fn的0.17~0.33倍。空化取决于局部绝对压力是否低于汽化压力，结构振动取决于激励与结构响应，二者均非仅由旋流决定。缓解策略分为被动几何修改与主动流动控制两类：尾水管导流栅、导流装置及锥体修改可削弱旋流或壁面载荷，但固定几何在偏离目标工况时可能损害性能；补气可提高涡核压力并降低脉动，但需可控气源且可能引入效率损失或气液两相流复杂性；轴向水注入通过增加轴向动量、削弱维持旋转涡带的环形剪切层，更直接地作用于低动量核心。自驱动轴向射流概念由Muntean、Susan-Resiga、Bosioc、T?nas?等开创，少量旁通流从蜗壳经水轮机轴路由至锥体顶点喷出，由水轮机压差驱动而无需辅助泵。前人已在旋流扩散器及水泵水轮机透平工况中证实压力脉动衰减与涡带缓解。然而三个问题尚未充分解决：注入能量相对于原始涡带损失在何处耗散；孔口半径如何平衡抑制效果与附加射流混合损失；控制后相干涡拓扑与压力谱如何变化。
针对上述问题，研究人员开展研究：采用稳态RANS计算筛选射流半径以供非定常ELES对比，考察涡形态与截面计数，并解读低频压力响应。参照Herwig与Kock的平均/脉动速度分解，利用瞬时解析速度评估平均与脉动速度梯度对熵产的贡献。该诊断用于解读核心区附近解析耗散空间分布的变化，但不确立总热力学损失，因为未评估未解析尺度耗散及独立的近壁熵产处理。研究所用主要关键技术方法包括：稳态SST k–ω湍流模型计算用于六种射流半径筛选；非定常嵌入大涡模拟（ELES）在RANS区采用SST k–ω、LES区采用Smagorinsky型模型，RANS–LES界面应用ANSYS Fluent中的Vortex Method生成解析速度脉动；基于Herwig–Kock平均/脉动速度梯度分解的熵产诊断（解析梯度LEGR）；基于Liutex准则（MLTX=20）的涡形态可视化与45个截面（P1–P45，间距5 mm，跨度220 mm）涡分量计数；以及压力脉动监测与频谱分析。研究对象为某电站模型试验关联的水轮机，计算域含14个固定导叶、28个活动导叶、15长15短叶片转轮及肘形尾水管，转轮进/出口直径分别为0.630 m与0.349 m。锥体流道半径40 mm，考虑4~9 mm六种孔口半径。采用H=11.87 m、hl/H=0.10、g=9.81 m s?2、Q=0.13962 m3 s?1得Vjet≈14.48 m s?1，六种半径对应射流流量份额约0.52%、0.81%、1.17%、1.60%、2.08%、2.64%。时间步采用转轮角度法（每步2°），fn≈5.55 Hz时Δt≈1.00×10?3 s，共计算10000步约10 s，前2160步排除于平均与谱分析，剩余名义记录约7.84 s，原始Fourier频率间隔约0.128 Hz。LES Index of Quality（LES_IQ）用于估计解析湍动能份额，LES_IQ>0.8表示解析动能超过解析加模化总量的80%。与匿名电站模型试验测量对比，转轮力矩、净水头与水力效率的最大相对偏差分别为4.02%、9.84%、7.73%。
研究结果方面，第一，稳态RANS初步筛选：稳态SST k–ω计算仅用于通过稳态轴向流型与模化耗散分布筛选六种孔口半径。无射流稳态计算中锥体下方出现Vz≤0的核心，相邻轴向剪切层与高模化LEGR区域重合，该重合标识稳态解中强剪切与混合的可能位置，但不表明熵产维持相干进动模式。Rjet由4增至6 mm逐步填充轴向速度亏损；7 mm时图示回流区消失、肘部入流更均匀，但直锥混合区仍存在强模化LEGR区；8与9 mm显示射流下游影响更大。这些观察促成选择7 mm进行非定常对比，但等值线视觉范围或峰值颜色不能确立积分熵损失或净效率收益。7 mm半径系按所述稳态流动准则选出的候选方案，而非已证实的全局或热力学最优值。第二，非定常结果与分析：（1）涡形态演化：采用Liutex准则MLTX=20可视化，无射流工况显示相干螺旋结构；在所述可视化阈值下，7 mm工况未显示连续螺旋结构，解析涡图案因注入而改变，此观察与相干核心组织改变一致，但不表明所有结构在肘部之前衰减，下游分量计数需单独且更有限的解读。截面计数在MLTX=20阈值下：无射流曲线在P5附近升至约29后总体下游下降，肘部有小残余值，扩散段趋近于零；7 mm射流曲线直锥段计数约1~3，但在肘部后段与扩散段增至约5~16。因此图示上游核心附近检测分量更少而下游检测分量更多，不支持下游减少或肘部前完全耗散。上游计数降低本身不能证明螺旋全局模式被抑制，下游计数增大本身也不能确立更强压力激励。（2）涡数量变化：计数代表阈值化截面分量，不度量三维相干涡的能量或尺寸；须与形态及压力响应分别考虑。（3）压力对比：所报道主导频率约1.6 Hz（0.29 fn），处于常报涡带范围，3.2 Hz峰与二次谐波相容但非独立验证。7 mm工况下1.6 Hz幅值DT1由249降至167 Pa、DT2由110降至71 Pa，约33%与36%，残余分量仍有限。图10中表观平均压力变化须与振荡幅值分离，局部平均静压偏移可反映运行/参考压力变化、速度重分布或物理压力恢复变化，仅凭两个测点不能识别为总水头损失；两面板显示不同时间跨度，共同基准比较前须核查出口与运行压力、静水基准、主与旁通流量、采样窗口及压力定义。显示0~20 Hz频带内fn附近未报道显著新峰，但图示不能证明更高频激励不存在；名义1 ms时间步若每步记录Nyquist限约500 Hz，该采样限不证明该频率以下数值精度；至少200 Hz的谱将涵盖名义15 fn≈83.3 Hz与30 fn≈166.5 Hz通道尺度，此类谱及其收敛性未在此评估，亦不能从低频图推断结构安全或空化抑制。（4）水力惩罚与运行范围适用性：名义旁通假设下7 mm射流Qjet/Q≈1.6%，动能通量约1.44% ρgQH；物理旁通在转轮前抽取水力输入，部分在供应路径耗散，其余再注入尾水管，净单位效率变化需总供应流量、转轮力矩及一致的总压或水头平衡，局部静压偏移与解析梯度LEGR等值线均不能提供该平衡。该概念针对部分负荷运行，此处无受控BEP或HL计算，旁通在这些工况保持开启的影响未知；全范围适用性需额外工况对比与耦合旁通特性及射流交互区充分解析，当前结果未确立可投产的自驱动供应系统或全运行范围效率收益。
讨论与结论方面，稳态六半径筛选与非定常无射流/7 mm对比表明，在所报道部分负荷工况下轴向注入可改变中心速度亏损与相干涡形态。7 mm半径为按稳态流动筛选准则选出的候选方案，名义Qjet/Q约1.6%，尚未确立为热力学最优值。所报道1.6 Hz压力幅值降幅在现有数值设置下DT1约33%、DT2约36%。将Herwig–Kock平均/脉动梯度分解应用于瞬时ELES速度为解析梯度熵产提供诊断，观测到的空间变化表明注入相关的解析耗散重分布；未包含独立SGS熵项与壁面函数熵修正，故分析未闭合总损失平衡，也未证明任一省略贡献可忽略。图9按其图例解读显示射流工况上游分量更少而下游分量更多，这阻止将计数解释为上游完全耗散的证据。名义射流动能通量1.44%以假设10%旁通水头损失为条件，并非实测净效率惩罚。模型试验对比评估整体性能，而受控非定常流动与旁通损失特性的验证需进一步证据。RANS–LES界面采用Vortex Method，射流流量由供应至锥体压差决定，这些实现选择本身不确立对局部网格分辨率、时间步或界面调整长度的独立性。后续工作需检验速度梯度统计、在寻求总预算处量化省略损失贡献、考察Liutex阈值敏感性、核查压力基准与采样程序并扩展压力谱；在确立设计优化或工程适用性之前，仍需非定常半径对比、完整效率平衡及额外运行工况。
