某水库溢洪道明渠扩散段三角翼体型参数设计优化

黄智峰，鄢慧欢

(江西华钦建设工程有限公司，江西 南昌 330038)

水利水电工程中，扩散过渡段的位置和形态会直接影响到工程的安全运行，必须认真设计和施工[1]。当前，在水利水电工程设计中，为了满足国家和行业标准中关于扩散段扩散角的要求，一般会采取增加过渡段长度的方式。但是，在部分水利水电工程建设中，受坝址区地质以及空间限制，不得不使用较短的扩散段长度，而选择大扩散角过渡[2]。但是，相关学者的理论研究和工程实践也证明，如果扩散段的扩散角过大，边界对水流的影响作用就会显著增强，进而导致水流断面的能量分布不均，导致显著的侧壁回流、边界脱流以及主流偏流等问题，从而影响水工建筑物的正常使用，情况严重的还会危及水工建筑物的安全[3]。由此可见，对于没有足够空间布置小扩散角过渡段的水利水电工程，如何对大扩散角条件下的水流进行控导和干预是亟待解决的工程问题，而这也正是当前水利工程建设领域的薄弱之处，该领域的研究成果也必将具有巨大的理论和实践应用价值。相关工程经验表明，在扩散段设置潜没式三角翼能够对该部位的水体流动起到有效的控制作用，属于比较有效的控导手段，但是如果三角翼设置不合理，仍会产生比较明显的回流区[4]。因此，有必要就三角翼对扩散段的水流控导效果进行深入研究。

1 计算模型的构建

1.1 几何模型的选择

本次研究的原型选择某水库工程溢洪道明渠扩散段，该水库位于江西省，下游距离县城约35km。该水库的主要任务是城市供水，同时兼具防洪和下游农业灌溉条件的改善。水库的设计库容为1.6亿m3，工程等别为Ⅱ等。水库的永久性建筑物主要包括大坝、副坝、溢洪道和输水设施。考虑本文研究的重点与计算的复杂程度，几何模型选择的是溢洪道扩散段以及上游的直线窄段50m和下游直线宽段120m。

1.2 模型的构建

模型创建是流体数值模拟研究的基础和关键环节。结合本次研究对象的实际特点以及相关研究经验，采用ANSYS三维建模软件进行研究对象的几何建模，利用FLUENT软件进行数值计算模型的构建[5]。一般情况下，模型的网格划分越密集，计算精度越高，但是计算量也会急剧增大，因此，本文对三角翼部分采用混合网格划分[6-8]，最终获得294110个网格单元，244268个计算节点，过渡段混合网格示意图如图1所示。

图1 过渡段混合网格示意图

1.3 边界条件

结合本次研究的实际情况以及平面大扩散角衔接过渡段流动数值模拟研究的相关成果，设定如下边界条件：对于模型的入口选取液相速度入口边界条件与气相压力入口边界条件[9-10]。具体而言，根据模型入口部位的设计流量，计算入口边界均匀来流条件。对于模型的出口，设置压力出口边界条件，根据某水库下游的实测水位设置水流自由面的高程，本次研究的模拟工况全部为平面大扩散段两相流问题，模型下游出口处的水流可以自由发展，没有回流现象出现；模型设置为无滑移的固壁边界条件。

2 三角翼数值模拟优化研究方案

2.1 优化研究的关键参数

扩散段三角翼的几何体型特征主要受如下参数控制：①分流度(α)，一般情况下该角应设定为单侧扩散角的2倍，同时要保证与来流相适应；②三角翼的潜没度(h)，增大三角翼的高度有利于提高水流控制能力，但同时会增大水力损失，从而在三角翼的后部产生回流区；③三角翼的长度，也就是控导度(η)，该参数主要影响三角翼的控导范围，可以对扩散段的水流流态以及水流动量的调整强度产生显著影响；④三角翼的攻角，也就是攻流度(θ)，也就是三角翼的翼弦与扩散段来流之间的夹角，其大小可以对匹流和控导效果产生直接影响。此外，三角翼的丰满度(b)、回流度(φ)以及位置度(λ)也是影响三角翼几何体型的重要因素。

2.2 模拟研究方案设计

目前，水利工程研究中三角翼的体型参数优化数值模拟方案并不明确，因此，本次研究基于上节提出的主要控制参数，探索三角翼体型优化设计的数值模拟方案，为研究的顺利进行打下良好基础。在工程研究中，控制变量法是针对多因素问题的有效研究方法，理论上，我们可以通过上述方法寻求每个三角翼控制参数的最优值，并最终确定出三角翼体型设计的最佳方案。在具体的研究方案设计中，结合相关演技成果与方便后期结果分析的考虑，仅对控导度、攻流度2个无量纲控制变量进行优化研究，研究方案见表1。

表1 计算方案设计

3 模拟计算结果分析

3.1 攻流度模拟计算结果

利用上节构建的模型，对不同攻流度条件下的溢洪道扩散段流场进行模拟计算，结果显示：利用a型三角翼体型设计时，水流的流态和整体分布比无控导条件下表现更为稳定，回流区长度明显减小，强度也有所降低。但是在边壁的部分区域产生一定的边界层分离现象，两侧均出现了范围较大的不对称型漩涡，主流位置存在一定程度的偏移，流速的不均匀性仍旧较大；在利用b型三角翼体型设计时，水流的流态和整体分布与a型基本相同，但是两侧漩涡的范围区强度有所增加且中心有向下推移的趋势，下游主流受两侧剪切挤压，摆动不稳定现象比较明显；在利用b型三角翼体型设计时，两侧边壁部位的动量明显增加，边界层分离不复存在，下游主流偏移现象消失，但是在翼后方存在范围不大的回流区，使单宽动量不均匀度增大，影响范围也显著增加，呈现出矫枉过正的情况。

为了对扩散段三角翼控导效果进行定量分析，根据模型计算结果绘制出不同计算工况下的各个计算断面的水流单宽动量不均匀度变化曲线，结果如图2所示。由图2可知，b型三角翼的控导效果最优，因此，就本次研究的某水库溢洪道扩散段而言，攻流度的优化取值为0.75，也就是取67.5°的攻流角。

图2 不同攻流度断面单宽动量不均匀度

3.2 控导度模拟计算结果分析

利用上节构建的模型，对不同控导度条件下的溢洪道扩散段流场进行模拟计算，结果显示：利用d型三角翼体型设计时，扩散段边壁位置不存在明显的边界分离，主流也不存在偏移现象。但是，在三角翼的后方存在一定范围的汇流区，水流流速呈现出两侧大中间小的特征，由于两侧水流速度较大，造成左右激荡作用明显；利用e型三角翼体型设计时，在溢洪道扩散段的左侧存在呈条带状分布的回流区，并且随着三角翼长度的增加，在三角翼的后部和右侧出现了一定范围的弱流区，在一定程度上降低了扩散段的过流能力。

为了对扩散段三角翼控导效果进行定量分析，根据模型计算结果绘制出b、d、e型工况下各个计算断面的水流单宽动量不均匀度变化曲线，结果如图3所示。由图3可知，d型三角翼的控导效果最优，因此，就本次研究的某水库溢洪道扩散段而言，控导度的最优取值为0.4，也就是三角翼的长度为4.0m。

图3 不同控导度断面单宽动量不均匀度

4 结论

本次研究以某水库溢洪道扩散段为工程背景，利用CFD和FLUENT软件对大扩散角条件下的三角翼优化设计进行数值模拟研究。结论对推动大扩散角过渡段水流调整以及边界层控制流域的理论完善，促进复杂边界条件下的输水工程设计具有重要理论意义和实践价值。但是，由于三角翼轮廓和造型的复杂性，使数值模型研究中的网格剖分仍存在一些问题。例如混合网格剖分虽然可以较好体现计算域的完整性，但是会造成交界面网格质量的降低。今后研究中需要探寻更好的网格剖分方式，避免网格剖分质量带来较大的计算误差。