Corresponding author: LIU Chun-bao, E-mail: liuanbc@126.com
随着液力传动技术和控制技术的发展,以导叶可调式液力变矩器作为调节单元组成的液力控制系统在风力发电等领域应用日益广泛,通过调节导轮叶片的角度可以实现对整个传动系统的控制[1, 2].
在液力变矩器早期研究阶段,主要基于一维束流理论对液力变矩器的性能进行研究,通过简化流动模型结合一些经验或半经验计算式来预测液力变矩器的特性.随着液力传动技术的发展,液力变矩器的性能研究已经发展到了三维流动理论阶段[3],同时较多的学者采用CFD的方法,对液力变矩器内流场进行研究[4, 5, 6].哈尔滨工业大学毕强采用CFX软件对导叶可调LB46型液力变矩器的单流道模型进行了内流场数值计算[7];重庆大学谢清乐采用Fluent软件对LB46型液力变矩器的全流道模型进行数值模拟,并与试验测得的原始特性进行比较[8].已有文献只是对导叶可调式液力变矩器的内部流场进行了CFD数值模拟,缺乏内流场试验验证,无法证明其计算结果的准确性.目前对导叶可调式液力变矩器多开度、多工况下的内流场进行分析的文献较少.
本文利用CFD数值模拟方法,结合试验验证,对导叶可调式液力变矩器多开度、多工况下的内部流场进行分析.选择标准k-ε湍流模型,对循环圆直径为320 mm的导叶可调式液力变矩器进行CFD数值模拟,揭示不同开度、不同工况下,其内部流动状况,分析可调导轮的内部流动规律,结合PTV试验验证数值仿真结果的准确性,能够为导叶可调式液力变矩器结构参数优化和性能提高提供依据.
1 导叶可调式液力变矩器工作原理可调式液力变矩器是一种人为或自动调节叶片进出口安装角(叶片转动)而改变其外特性的一类液力变矩器.常用的调节方式有两种,一种是改变泵轮叶片的角度,另一种是改变导轮叶片的角度.导叶可调式液力变矩器是采用调节导轮叶片的方式,改变整个导叶可调式液力变矩器内部流场分布,最终改变其外特性.图 1a为导叶可调式液力变矩器结构简图,其循环圆直径为320 mm,与一般的向心涡轮式液力变矩器不同,该导叶可调式液力变矩器由离心式泵轮、离心式涡轮、轴流式固定导轮以及向心式可调导轮构成.向心式可调导轮的叶片能够通过内调节式调节机构,实现可调导叶的自由转动,如图 1b所示.
根据所设计的导叶可调式液力变矩器结构参数进行三维建模,抽取三维流道模型,如图 2a所示.非结构网格采用随机的数据结构,有利于网格自适应,容易控制网格的大小和节点的密度,生成高质量的网格.因此,本文采用非结构化网格方法划分网格,如图 2b所示.泵轮、涡轮、固定导轮和可调导轮的网格数量分别为390 228,159 019,619 152,1 907 579.
普通三元件的液力变矩器,其内部流场只随着工况的变化而变化.导叶可调式液力变矩器,其内部流场除了随着工况的变化外,还随可调导叶的角度变化而变化.目前,对导叶角度变化过程中的液力变矩器内流场进行动态数值模拟还难以实现,只能选定某一固定的导叶角度,然后对其各个工况下的内部流场进行数值模拟.设x为可调导轮的开度,θ为叶片的开启角,θmax为最大开启角,定义开度x =θ/θmax,且θmax=40°.本文选取两个开度,即θ=20°,x=0.5,与θ=40°,x=1,如图 3所示.
运用CFD数值模拟方法对导叶可调式液力变矩器进行内部流场数值模拟时,由于算法的选择对模拟的效率以及其正确性有很大的影响,因此需要选择合适的算法.雷诺平均模型中的标准k-ε模型,在工程上应用广泛,能够用于模拟导叶可调式液力变矩器三维内部流场[9, 10].
本文选用SIMPLE算法,湍流模型选择标准k-ε模型.采用滑移网格方法在导叶可调式液力变矩器的各个叶轮进、出口面上都设置为interface,各自分界面也随时间相应转动,能量通过分界面实时传递,从而实现叶轮间流动参数的实时耦合,其他边界都为wall.计算时各叶轮转速按照实际的运动情况定义.采用非稳态求解器收敛计算,时间步长为5×10-2 s,收敛精度为10-4.为便于CFD数值模拟结果与PTV试验结果对比,CFD数值模拟时液体介质为水.定义流体介质水的属性:密度ρ=998.2 kg/m3,黏度μ=0.001 003 Pa·s,参考压力为0.1 MPa.
3 内部流场分析根据CFD计算结果,重点对可调导轮在制动、设计、空载三种典型工况下的流场进行分析;在同一工况下,对比分析了开度变化对可调导轮内部流场的影响.
3.1 制动工况图 4为制动工况下,CFD数值模拟所得的不同开度可调导轮内部速度矢量图.制动工况i=0时,涡轮静止不动,液流经过固定导轮后冲击可调导叶的工作面,与叶片进口角方向形成一定的正冲角,液流从入口至出口液流方向变化程度较小.开度x=0.5时,流道入口处流动状况比较紊乱,有漩涡现象产生,而流道出口处流动分布均匀.由于液流进入叶片时冲击叶片产生冲击损失,在叶片进口处形成一个阻滞区,并且延伸到流道内部,阻滞区的高压作用使得液体在叶片工作面入口处形成漩涡现象.
随着开度的增加(开度x=1),可调导轮叶片的进出口角均增大,进口液流方向与导轮叶片骨线的夹角减小;由于开度增大使过流面积增大,液流出口速度降低,出口液流角随着开度的增加而增大.由于开度增大,流道内的液流速度降低,流动情况得到改善,漩涡区域明显减小.
3.2 设计工况图 5为设计工况下,CFD数值模拟所得的不同开度可调导轮内部速度矢量图.设计工况i=0.5时,可调导轮内部流动状态比较理想,流道内的流场分布比较均匀,流动方向基本上是从进口到出口,整个流道内没有出现漩涡等复杂流动现象.开度x=0.5时,液体冲击叶片的工作面,流速较低,压力较高,而非工作面流速较高,压力较低,速度变化大.
随着开度的增加(开度x=1),可调导轮叶片的进、出口角均增大,可调导轮进口速度变化较小;出口液流角随着叶片出口角的增大而增大,出口液流速度明显降低.
3.3 空载工况图 6为空载工况下,CFD数值模拟所得的不同开度可调导轮内部速度矢量图.空载工况i=1时,叶片非工作面出现了明显的漩涡现象,叶片工作面的流动则相对均匀,叶片入口处液流的冲角较大,液流从进口至出口的液流方向变化程度最大.开度x=0.5时,液流以较大的流速经过叶片入口处的圆头时,在叶片非工作面产生分离,使得叶片非工作面入口处的压力较小,附近液流相对压力较高从而导致漩涡的产生.在叶片中段以后,由于叶栅的作用逐渐加强,漩涡规模受到抑制不能扩展,从叶片中段到出口流道内液流趋向平缓.
随着开度的增加(x=1),可调导轮叶片的进、出口角均增大,出口液流角随着开度的增加而增大;液流进口速度变化较小,出口速度明显降低.随着开度的增加,叶片非工作表面的漩涡区域明显增大.
根据上述三种典型工况下的流场分析结果,在同一工况下,随着开度的增加,液流通过固定导轮后,进入可调导轮的速度基本不变,出口速度降低,出口速度的方向随着叶片角度的变化而变化.同一工况下,不同开度可调导轮内部流场结构类似.
4 仿真与试验对比分析 4.1 PTV试验测试粒子跟踪测速(particle tracking velocimetry,PTV)作为一种先进的非接触测量手段,已经被广泛应用于复杂元件的内部流场试验研究[11, 12, 13].试验测试不仅能了解液力变矩器流场的基本流动过程,还能为数值计算提供边界条件,并对数值计算结果的可靠性、准确性进行验证[14, 15].
对所设计的透明导叶可调式液力变矩器进行了PTV内流场试验.根据图 1所示导叶可调式液力变矩器结构图,按1∶1比例制作了相应的实物模型,如图 7所示.其循环圆直径为320 mm,壳体及叶片材料均为有机玻璃,其余为不锈钢1Cr13.导叶可调式液力变矩器运转时,液流在泵轮、涡轮、固定导轮及可调导轮间循环流动.具体测试方法、过程和结果参见文献[16].
试验测量系统由机械部分、光学部分、图像采集与显示部分组成.机械部分由交流电机、盘式制动器、齿轮泵、转速测量仪等组成.光学部分由300 mW氩离子激光器、光学组件、调速电机和脉冲盘等组成,用于产生宽为1~2 mm的连续可调激光脉冲片光源.通过调节调速电机控制脉冲盘能够调控示踪粒子在流场中的曝光时间.选用铝粉作为示踪粒子,粒径约为0.04 mm.图像采集与显示部分包括高清摄像机、图像采集卡和与之配套的采集软件.
4.2 CFD仿真结果与PTV试验对比分析由于同一工况下,不同开度可调导轮内部流动规律类似,在进行CFD仿真与PTV试验对比分析时,只选取开度x=0.5为例进行对比分析.
4.2.1 制动工况图 8为制动工况下开度x=0.5时,CFD数值模拟与PTV试验所得的可调导轮内部速度矢量图对比.比较发现其内部流动规律基本相符:入口液流冲击可调导叶的工作面,与叶片进口角方向形成一定的正冲角,流道入口处流动状况比较紊乱,而流道出口处流动分布均匀.在CFD仿真速度矢量图A处,液体流动较为紊乱,形成漩涡,对应的试验速度矢量图B处也有明显的漩涡现象.
图 9为设计工况下开度x=0.5时,CFD数值模拟与PTV试验所得的可调导轮内部速度矢量图对比.比较发现其内部流动规律基本相符:入口液流冲击可调导叶的工作面,沿着叶片弯曲方向,由叶片入口流向叶片出口,流道内速度分布比较均匀,液流的冲击较小,没有出现漩涡等复杂流动现象.在CFD数值模拟速度矢量图C处,流动分布均匀,与对应的PTV试验速度矢量图位置D处流动现象一致.
图 10为空载工况下开度x=0.5时,CFD数值模拟与PTV试验所得的可调导轮内部速度矢量图对比.比较发现其内部流动规律基本相符:入口液流冲角较大,液流以较大的速度流过叶片入口处时,在CFD仿真速度矢量图E处产生明显的漩涡流动现象.随着液流由入口流向出口,叶栅的作用逐渐加强,叶片工作面的流动则相对均匀.在对应的PTV试验速度矢量图F处同样可以看到叶片非工作面入口附近的漩涡现象.
对比结果表明,采用雷诺平均模型中标准k-ε模型的CFD方法可以准确预测导叶可调式液力变矩器内部流场分布情况.通过观察内部流场的复杂流动情况,诸如紊流、漩涡等,分析得到了流场能量损失原因.可以通过改变其结构参数、减小内部流场损失来提高液力变矩器性能.由于液力变矩器内特性决定了其外特性,CFD数值模拟方法能够为优化导叶可调式液力变矩器结构参数,提高其外特性提供指导.
5 结论1) 通过对CFD计算结果分析可知:制动工况下,随着开度增大,漩涡流动区域变小;设计工况下,开度变化对流动状态影响较小;空载工况下,随着开度增大,漩涡流动区域明显增大.
2) 将PTV试验与CFD数值模拟对比分析发现,在制动工况和空载工况下,CFD数值模拟速度矢量图和PTV试验结果矢量图均能观察到明显的漩涡现象,而设计工况流动状态比较稳定.两种方法得到的可调导轮内部流动变化趋势基本相同,故CFD数值模拟方法可以有效预测导叶可调式液力变矩器的内部流动规律,为其结构设计优化提供理论依据.
[1] | Vries E O.Drive and innovation:the DeWind D8.2 with Voith WinDrive[J]. Renewable Energy World, 2007, 10(2):55.(1) |
[2] | Muller H, Poller M, Basteck A, et al.Grid compatibility of variable speed wind turbines with directly coupled synchronous generator and hydro-dynamically controlled gearbox[C]//The Sixth International Workshop on Large-Scale Integration of Wind Power and Transmission Networks for Offshore Wind Farms.Delft:Brussels, 2006:307-315.(1) |
[3] | 马文星,何延东,刘春宝.液力传动研究现状分析与展望[J]. 农业机械学报, 2008, 39(7):51-55. (Ma Wen-xing, He Yan-dong, Liu Chun-bao.Situation and prospects of research on hydrodynamic transmission[J]. Transactions of the Chinese Society for Agricultural Machinery, 2008, 39(7):51-55.)(1) |
[4] | Kim B S, Ha S B, Lim W S, et al.Performance estimation model of a torque converter part I:correlation between the internal flow field and energy loss coefficient[J]. International Journal of Automotive Technology, 2008, 9(2):141-148.(1) |
[5] | Jung J H, Kang S, Hur N.A numerical study of a torque converter with various methods for the accuracy improvement of performance prediction[J]. Progress in Computational Fluid Dynamics, an International Journal, 2011, 11(3):261-268.(1) |
[6] | Wozniak M, Pawelski Z, Dela F P, et al.The mutual influence and blade-row interaction between pump and turbine in a hydrodynamic torque converter[J]. Archiwum Motoryzacji, 2012(1):83-93.(1) |
[7] | 毕强.LB46型可调式离心涡轮变矩器内流场数值计算及分析[D]. 哈尔滨:哈尔滨工业大学, 2011. (Bi Qiang.Numerical analysis on LB46 adjustable guide vane centrifugal turbine hydrodynamic torque converter[D]. Harbin:Harbin Institute of Technology, 2011.)(1) |
[8] | 谢清乐,李奇敏.导叶可调式液力变矩器的内流场数值分析[J]. 计算机仿真, 2013, 30(8):135-138. (Xie Qing-le, Li Qi-min.Numerical analysis on internal flow field of adjustable guide turbine torque converter[J]. Computer Simulation, 2013, 30(8):135-138.)(1) |
[9] | 刘春宝,马文星,朱喜林.液力变矩器三维瞬态流场计算[J]. 机械工程学报, 2010, 46(14):161-166. (Liu Chun-bao, Ma Wen-xing, Zhu Xi-lin.3D transient calculation of internal flow field for hydrodynamic torque converter[J]. Journal of Mechanical Engineering, 2010, 46(14):161-166.)(1) |
[10] | Fuente P D L, Stoff H, Volgmann W, et al.Numerical analysis into the effects of the unsteady flow in an automotive hydrodynamic torque converter[C]// Proceedings of the World Congress on Engineering.London:Newswood Limited, 2011:2405-2410.(1) |
[11] | Mikheev A V, Zubtsov V M.Enhanced particle-tracking velocimetry (EPTV) with a combined two-component pair-matching algorithm[J]. Measurement Science and Technology, 2008, 19(8):817-822.(1) |
[12] | Stanislas M, Okamoto K, Kahler C J, et al.Main results of the second international PTV challenge[J]. Experiments in Fluids, 2005, 39(2):170-191.(1) |
[13] | Kreizer M, Ratner D, Liberzon A.Real-time image processing for particle tracking velocimetry[J]. Experiments in Fluids, 2010, 48(1):105-110.(1) |
[14] | Adrian R J.Particle image velocimetry:a practical guide[M]. Berlin:Springer-Verlag, 2007:90-99.(1) |
[15] | Cavazzini G.The particle image velocimetry-characteristics, limits and possible applications[EB/OL]. (2012-05-23)[2016-03-15]. http://www.intechopen.com/books/the-particle-image-velocimetry-characteristics-limits-and-possible-applications.(1) |
[16] | Ma W X, Liu H, Liu C B, et al.Measurement and analysis of the internal flow field in the hydrodynamic torque converter with adjustable guide vanes[J]. Advances in Mechanical Engineering, 2015, 7(6):1-13.(1) |