基于大涡模拟贯流式水轮机尾水涡带研究.docVIP

基于大涡模拟贯流式水轮机尾水涡带研究 　　摘 要 基于大涡模拟方法并采用Smagorinsky-Lilly动态亚格子模型，分别计算了灯泡贯流式水轮机的6个不同工况的非定常流场，得到了每个工况尾水管段尾水管涡带信息。研究发现：尾水管中心区域的脉动强度大于外围。水头对尾水涡带形态影响显著，随着水头的降低，涡带形状由双螺旋变成单螺旋，直至平直状。 　　关键词 大涡模拟；贯流式水轮机；压力脉动；涡带 　　中图分类号 TK7 文献标识码 A 文章编号 1674-6708（2015）150-0147-02 　　1963年，美国气象学者 Smagorinsky[1]为研究天气预报问题首次提出了大涡模拟（Large Eddy Simulation）方法[2]。该方法认为湍流运动是由不同尺度的涡组成，通过滤波计算，将湍流的瞬时运动分解成小尺度涡和大尺度涡两部分。其中，大尺度涡拥有主导的湍流动能，通过控制方程直接求解；小尺度（亚格子尺度）涡主要起耗散作用，在高雷诺数下趋于各向同性，受边界条件影响较小，可用通用模型进行模拟。因此，亚格子模型是大涡数值模拟控制方程封闭的前提，是大涡模拟方法的核心。Smagorinsky提出了第一个亚格子模型，Lilly[3]改进并发展了Smagrinsky-Lilly动态亚格子模型，其中Smagorinsky模型系数根据解析尺度提供的信息动态计算，由当时当地的速度场动态处理，故而能动态地反映湍流信息，且与实验结果符合良好。 　　目前，LES在湍流的数值计算中获得了大量的关注，我国的研究者也取得了丰硕的成果。茅媛婷[4]将LES方法应用于混流式水轮机的研究，取得了与实验吻合较好的结果，夏晨宇[5]将大涡模拟方法应用于贯流式水轮机压力脉动的研究，然而大涡模拟方法用于贯流式水轮机尾水涡带的研究并不多见。本文以灯泡贯流式水轮机为研究对象，重点针对尾水涡带进行了研究，为大涡模拟的应用和贯流式水轮机尾水涡带的研究做了积极 　　探索。 　　1 几何模型与计算方法 　　1.1 计算区域与计算工况 　　基于圆柱坐标系建立三维叶片，进而建立灯泡贯流式水轮机全流道模型，转轮直径D1=7.4m，轮毂直径d=3.1m，叶片数Z=5，导叶数16。计算控制域即水轮机的全流道流场空间，在进行了网格无关性验证的基础上，采用非结构三维网格技术划分网格。水轮机发电机额定转速n=83.33r/min，旋转方向为向下游顺时针。本文按照水头高低，选择了运行范围内从最高到最低水头的6个不同工况（详见表1），进行了三维紊流非定常流场 　　计算。 　　1.2 计算求解方法 　　1.2.1 滤波方法 　　大涡模拟第一步就是滤波，把流动分成小尺度量和大尺度量。滤波运算通过消除湍流中小尺度脉动，滤掉比过滤网格小的漩涡，从而得到可解尺度运动的控制方程。过滤方法如下： 　　（1） 　　其中V为计算单元的体积，是计算变量，是经过过滤后的变量，G为决定过滤尺寸的函数，是实际流动区域中的坐标，x为过滤后的大尺度空间坐标。过滤过滤函数选择盒式滤波器。 　　1.2.2 计算方法 　　为了节省计算时间，首先用RANS方法定常计算出一个初始流场，再改用LES方法计算非定常的流动状态，LES计算直到湍流变得统计地稳定。 　　非定常阶段采用大涡模拟计算方法，亚格子模型采用Smagorinsky-Lilly动态模型，用以模拟小尺度涡的运动。由于水轮机的进口压力比较稳定且易于获得，因此采用压力进出口边界条件，压差为水头，考虑吸出高度。固壁采用无滑移边界，动静干涉选择滑移网格模型。考虑到转轮频率的关系，采用 0.002s 为时间步长，即每个时间步转轮转过1°。离散化方法采用有界中心差分格式，压力差值选择PRESTO！格式，压力速度耦合选择PISO算法。 　　2 尾水管涡带特性分析 　　通过设置等压面，选择合适的压强，显示尾水涡带形态。 　　（a） 工况1 （b） 工况2 　　（c） 工况3 （d） 工况4 　　（e） 工况5 （f） 工况6 　　图1 各工况下尾水管涡带形态图 　　图1（a）显示工况1尾水涡带以泄水锥为起点，由互相对称的两条螺旋线交错形成“麻花”状。该工况下转轮出流带有环量，水流旋转着进入尾水管，为涡带的形成创造了条件，但涡带处于尾水管中心区域，并没有发生偏心。此时旋转的水流在尾水管内受离心力的作用，而附着在尾水管壁面，导致流速放缓、湍流粘滞性上升，形成“贴壁效应”[6]，故而尾水管内涡带外围区域压力脉动不强，尾水管壁面也不会受到强烈的水力振动。工况2也是双螺旋交织形成的涡带，与工况1不同的是，两条螺旋的对称性不明显，大小相差较大形成主次之分。大螺旋与工况1的相似，较粗壮，随着流场向下游发展而衰减变细直至消失。小螺旋与大螺旋交织在一