摘要
离心泵是一种通用水力机械，其内部流动情况一直是泵设计人员十分关注的问题，因为泵内流动的优劣直接影响泵的性能。离心泵叶轮的内部流动是很复杂的三维紊流流动，同时由于受旋转和叶片表面曲率的影响还拌有脱流、回流及二次流的现象，是流体工程中较难的试验研究和数值计算问题之一。
早在20世纪50年代，一些专家学者就开始尝试使用数值计算方法来预测叶轮中的流动情况。但具有完备形态的内流数值模拟，一般认为始于吴仲华教授的S1、S2两类相对流面理论后，叶轮机械内流数值模拟才得到了迅速发展。迄今为止，国内许多学者已经对离心泵叶轮内部三维紊流数值模拟进行了研究，也取得了一些成果，但并不理想，也缺少实验验证。近年来得益于计算机技术的高速发展，计算流体力学CFD(Computational Fluid Dynamics)发展很快，许多商用CFD软件应用非常广泛，在离心泵内部流场数值模拟上的应用也日见增多。通过CFD方法对离心泵叶轮内部流动进行数值模拟，了解液流在叶轮过流部件中的速度和压力分布，从而进行叶片的选型、设计和性能预测已成为现代泵技术的重要方法之一。
本文基于N—S方程和标准的k—ε双方程紊流模型，采用商用CFX-TASCflow计算软件对一离心泵叶轮模型进行了叶轮内部三维紊流数值模拟，得到了三个典型工况下、从轮毂到前盖板每一个断面上的速度和压力分布。通过与该叶轮内流场的激光粒子图象速度场仪——PIV(Particle Image Velocitymeter)实验结果比较，初步分析了离心泵叶轮的汽蚀和能量特性，为进一步完善泵设计理论、提高泵设计水平提供有益的参考。
本文得到国家自然科学基金项目“水泵转轮内流场激光测量研究”(No．59949010)的资助。
Centrifugal pump is a commonly hydraulic machine, whose inner flow has direct influence on pump's performance and is a considered problem by the pump designers. with the phenomena of dropping flow, backflow and secondary flow effected by circulation and surface curvature, the centrifugal pump's inner flow is very complicated 3d flow . Therefore, the centrifugal pump inner flow is one of the difficult problems both in testing research and numeric computation .
In the 1950s, some experts and scholars tried to use numeric computation method to predict the blade's inner flow. But it is thought that the mature pump inner flow numeric computation was started and then developed rapidly after the "two relative flow surface" (the S1 and S2) theory presented by professor Wuzhonghua. From then on, many domestic scholars have made much research on the numeric simulation of the centrifugal pump inner 3D turbulent, and they have made some progress, which were yet short of enough validation through experiment. In recent years, with the rapid development of the computer technology, CFD (computational fluids dynamics) develops rapidly, and many commercial computation softwares are applied widely. Some softwares have been used to compute the centrifugal pump inflow. Through the numeric simulation of the centrifugal pump inner flow by CFD, the distribution of the velocity and pressure in the pump are conducted to impeller design and hydraulic performance prediction, which has become one of the important modern pump design method.
In this paper, based on the Navier-Stokes equation and the k- turbulent model, the numerical simulation of the 3D turbulent flow in a centrifugal impeller, of which specific speed is 93. The distribution of the velocity and pressure in the blade-to-blade passage are presented at the design and off-design operating conditions. The simulation results were verifiable with the results measured by PIV. This study can be very helpful to improve and optimize the impeller design.
引文
[1] 窦国仁编．紊流力学．上、下册．北京：人民教育出版社，1982．
[2] 沈阳水泵研究所、中国农业机械化科学研究所主编．叶片泵设计手册．北京：机械工业出版社．1983，7．
[3] [美]S．V．帕坦卡著．传热与流体流动的数值计算．北京：科学出版社．1984，9．
[4] 金忠青，N-S方程的数值解和紊流模型[M]．南京：河海大学出版社，1989．
[5] 陈玉璞等编著．流体动力学．南京：河海大学出版社，1990，6．
[6] 冯汉民主编．水泵学．江苏农学院．水利电力出版社，1991，5．
[7] 赫尔姆特 舒尔茨著．[联邦德国]．泵原理、计算与结构．机械工业出版社．1991．12
[8] 是勋刚主编．湍流．天津：天津大学出版社．1994年3月．
[9] 周雪漪编．计算流体力学．北京：清华大学出版社，1995．
[10] 陈汉平主编．计算流体力学．北京：水利水电出版社．1995，11．
[11] 倪浩清，沈永明编著．工程湍流流动、传热及传质的数值模拟．北京：中国水利水电出版社．1996年7月．
[12] 日本农业土木事业协会[日]．泵站工程技术设计手册．北京：中国农业出版社．1998，12．
[13] 刘超主编．水泵及水泵站．北京：科学技术文献出版社．2003年3月．
[14] 陈池，袁寿其，金树德．离心泵叶轮内流计算方法综述．流体机械．1999．27卷第2期：35～39．
[15] 赵斌娟，王泽．离心泵叶轮内流数值模拟的现状和展望．农机化研究．2002年8月．第3期：49～52．
[16] 阎超．流体机械内部流动数值计算方法的新进展．流体机械．1994(8)：53～8．
[17] 吴玉林，曹树良等．离心泵叶轮中紊流的计算．水泵技术．1997．(1)：9～15．
[18] 李文广．离心泵叶轮内部三维紊流数值汁算．水泵技术．1997，5：20～29．
[19] 唐辉，何枫．离心泵内流场的数值模拟．水泵技术．2002年03期：3～8．
[20] 倪浩清．工程湍流模式理论综述及展望．力学进展．1996年5月25日．第26卷第2期：145～165
[21] 熊鳌魁．湍流模式理论综述．武汉理工大学学报．2001年2月．第25卷第4期：451～455．
[22] 杨华，刘超，汤方平，谷传纲，采用PⅣ研究离心泵转轮内部瞬态流场，水动力学研究与进展，2002年10月，A辑第17卷第5期：547～552．
[23] 董壮．三维水流数值模拟研究进展．水利水运工程学报．2002年9月第3期：66～73．
[24] 李文广．离心泵叶轮内部流动数值解的评价．水利学报．1999，12．第12期．
[25] 陈池，袁寿其，金树德．低比转速离心泵研究现状与展望．流体机械．第27卷第7期：29～34．
[26] 刘立军，徐忠．离心叶轮内部三维湍流流场的数值分析．工程热物理学报，1996，17(3)：296～300．
[27] 李巍，王国强，刘立军．离心叶轮内三维紊流流动数值模拟．上海交通大学学报．2000，1．第34卷第1期：143～147．
[28] 李海峰，吴玉林，赵志妹．利用三维紊流数值模拟进行离心叶轮设计比较．流体机械．2001年第29卷第9期：18～21．
[29] 李文广．采用贴体坐标计算离心泵叶轮内部流动．水泵技术．1998年6月：14～18．
[30] 齐学义，徐洁，杨从新．低比转速离心叶轮内部流动的数值计算．甘肃工业大学学报．2001，6．第27卷第2期．
[31] 刘顺隆，郑洪涛，孙海鸥，张天野．正交贴体坐标网格的生成．哈尔滨工程大学学报．1997年6月．第18卷第3期：42～46．
[32] 陈池，袁寿其，郑铭等．离心泵叶轮三维贴体网格生成．农业机械学报．2000年7月．第31卷第4期：50～53．
[33] 刘正先，曹淑珍，谷传纲，王玉明等．离心叶轮内三维湍流流场的数值计算与实验比较．流体机械．2000．第28卷第4期：9～12．
[34] 李海锋，穆忠波，吴玉林．半开式离心泵叶轮内三维紊流的数值模拟．水泵技术．2001，1：9～15．
[35] 陈池，袁寿其，郑铭．低比速离心泵叶轮内三维不可压湍流场计算．动力工程．2001，8．第21卷第4期：1346～1357．
[36] 李春，吕世和，陈康民．旋转三维紊流流场的数值分析．上海理工大学学报报．Vol．23
No．4 2001．第23卷第4期：319～321．
[37] 张扬军，陈乃祥，曹智鹏．流体机械粘性流场的数值模拟方法．清华大学学报(自然科学版)．1997．第37卷第3期：95～97．
[38] 张楚华，刘正先，谷传纲，苗永淼．旋转离心叶轮内湍流的非结构化网格数值计算．航空学报．1998(9)：556～559．
[39] 吴玉林，曹树良，戴江，孙自详，王琳．离心泵叶轮中紊流的计算．水泵技术．1997．1：9～16．
[40] 钟绍湘．离心泵叶轮分流面的数值方法．大连大学学报．2000年8月．第21卷第4期：15～20．
[41] 张士杰，袁新，叶大均．全三维粘性流场计算的快速收敛方法．清华大学学报(自然科学版)．2000．第40卷第10期：96～99．
[42] 曹立新．离心泵叶轮全三维势流的间接边界元法研究．农业机械学报．1999年7月第30卷第4期：23～25．
[43] 蔡兆麟，罗晟．用有限体积法计算叶轮机械内三维粘性流动．华中理工大学报．2000，9．第28卷第9期：72～75．
[44] 陈德江，王尚锦．应用有限元法计算离心式叶轮内部流场．应用力学学报．1999，3．第16卷第1期：27～32．
[45] 王国玉，左志钢，曹树良，赵令家等．水泵水轮机内部三维紊流流场计算与性能分析．水利水电技术．2001．第32卷第5期：28～30．
[46] 彭国义，曹树良，林汝长．水力机械转轮正问题数值研究与进展．甘肃大学学报．1997(3)：30～32．
[47] 徐宇，唐学林，吴玉林．水泵水轮机转轮内水泵工况紊流分析．水力发电学报．2000．第3期：75～83．
[48] 宋斐，刘东明，孙平．Tascflow软件包技术特点分析与应用研究．东北电力技术．2000年第4期：1～5．
[49] 杨建明，吴玉林，曹树良．流体机械中高雷诺数流动的大涡模拟．工程热物理学报．1998，3．第19卷第2期：184～188．
[50] 姚征，陈康民．CFD通用软件综述．上海理工大学学报．2002年第24卷第2期：137～144．
[51] 周晓泉，曹叔尤，瞿伦富，吴玉林，曹树良．旋转机械紊流数值模拟中的周期边界条件处理．水力发电学报．2002年第2期．总第77期：78～84．
[52] 曹树良，杨辅政，吴玉林．用代数应力紊流模型预估水轮机转轮内部三维流场．清华大学学报(自然科学版)．1998年第38卷．第4期：113～116．
[53] 吴玉林，曹树良，任燕雨，王琳．水轮机转轮内部的三维固液两相紊流计算．工程热物理学报．1997，9．第18卷第5期：580～583．
[54] 杨军虎，赵万勇，李金平．极低比转速叶轮内流体的流动分析和叶轮的设计．流体机械，2001．第29卷第7期：20～22．
[55] 王福军，刘善琨，骆大章，陈军．水泵叶轮内部的流动计算与图象模拟．农业机械学报．1997，3．第28卷第1期：58～62．
[56] 魏兵海，谢军龙，吴克启．叶轮机械设计与影响性能及内流的若干重要问题．工程热物理学报．2001，3．第22卷第2期：185～187．
[57] 杨策，蒋滋康，索沂生，老大中．用粘性体力方法计算轴流叶轮机械内部流场．应用力学学报．2000，12．第17卷第4期：40～44．
[58] 陈次昌，宋文武，杨昌明，朱兵清．离心泵三维设计的研究．农业机械学报．2002年5月．第33卷第3期34～36．
[59] 袁寿其，陈池，郑铭等．无过载离心泵叶轮内三维不可压湍流场计算．机械工程学报．2000，5．第36卷第5期：31～34．
[60] 周先进，高忠信，张世雄．水轮机转轮的三维粘性流的数值计算及其结果分析．水利水电技术．1999．第11期：18～22．
[61] 高忠信，周先进，张世雄，唐澎．水轮机转轮内部三维粘性流动计算与性能预测．水利学报．2001，7．第7期：30～35．
[62] 张武高，陈晓玲，曹广军．影响离心蜗壳内速度场分布的流量因素．流体机械．1999(12)：14～17．
[63] 蔡兆麟，罗晟．流体机械叶轮内二次流及尾迹发展．华中科技大学学报．2001年6月．第29卷增刊Ⅰ：8～10．
[64] 赖焕新，康顺，吴克启．离心式叶轮内部湍流及出口“射流—尾迹”结构的成因分析．推进技术．2000年2月第21卷第1期：45～48．
[65] 曹树良，吴玉林，杨辅政．混流式水轮机转轮内部三维紊流的数值分析．水力发电学报．1997年第4期：52～60．
[66] 齐学义，杨建明，吴玉林．混流式转轮三维湍流流动计算方法．水利学报．2000，1．第1划：67～71．
[67] 李连超，崔学明，常近时．湍流理论的发展及其在水轮机中的应用．水利水电家户技术．第33卷，2002年第10期：19～21．
[68] Jin Yan and D. G. Smith. CFD Simulation of 3-Dinentional fow in Turbomachinery Applications, FLUENT Technical Notes, Presented at Turbomachinery Flow Prediction Ⅷ. ERCOFTAC workshop, March 2000.
[69] Wu Y L, Dai J. Turbulence Flow Simulation Through Centrifugal pump impeller at Design and off-design Conditions.The 2nd Inter. Conf. Pumps Fans, Beijing, 1995: 155～167.
[70] Zhangliang, Liu shuhong, Wu yulin. THREE-DIMENSIONAL TURBULENT FLOW SIMULATION THROUFH A FRANCIS TURBINE, THE FOURTH INTERNATIONAL CONTERENCE ON PUMPS AND FANS. Beijing, China. 2002, 10: 530～534.
[71] Johnson M W, Moore J. The influence of flow rate on the wake in a centrifugal impeller[J]. J of Eng for Power, 1983, 105(1): 33～39.
[72] Launder B E, Spalding D B. The numerical computation of turbulent flows. Computer Methods in Applied Mechanics and Engineering, 1974, 3: 269～289.
[73] Wu, Y. L., Z. X. Sun, R. Oba, T. Ikohagi, Study on Flows Through Centrifugal Impeller by Turbulent Simulation, ASME FED-Vol. 222, Fluid Machinery, (1995), 63～68.
[74] Goede, E., Sebestyen, A. and Schachenmann, A. 1992, Navier-Stokes Flow Analysis for a pump impeller. Proc. of the 16th IAAR Symp: 513～523.
[75] Wiesner F J. A Review of Slip for Centrifugal Impellers. ASME J of Engineering for Power, 1967, 89(3), 558～572.