低比转速半螺旋吸水室双吸泵的CFD分析.docVIP

万　方数据 排灌 机械 第２６卷 扬程Ⅳ为２９．５ｍ，转速，ｌ为１为５８，叶片数彳为６． ４５０ ｒ／ｌＩｌｉｎ，比转速ｎ．３．２出口边界条件 由于出口压力和速度分布未知，出口边界采用第二类边界条件．３．３壁面条件 在叶片和轮毂等各固体壁面，均采用无滑移条件，对于近壁区的流体流动则采用标准壁面函数法．３．４压力边界条件 为了保证计算的稳定性，在计算域中，选择某一点的压力值作为参考压力，边界上压力法向导数为零‘５－６１． 所选用双吸泵及半螺旋形吸水室为对称结构，其对称面满足对称边界条件，所以计算区域选择吸水室和叶轮、泵体流道的一半‘１｜，见图１． ４计算结果及分析 通过对０．６倍、０．８倍、１．ｏ倍、１．２倍和１．４倍 图ｌ计算域与吸水室 Ｆｉｇ．１ 设计流量这５个流量点分别进行模拟分析，可以得 ｃｈ锄ｂｅｒ 到这５个流量点下的内部流场分布结果，同时可以对整体流道水力性能进行预测．４．１蜗壳内叶轮压力 图２和图３分别为设计流量点时对称面上叶轮和蜗壳流道内的静压和总压分布图． ＣａＪｃｕｌａｔｉｏｎａｒｅａ彻ｄｓｕｃｔｉｏｎ ２ 网格划分和计算模型的选取 采用Ｐｒｏ／Ｅ软件对半螺旋形吸水室、叶轮和蜗 壳部件内部流道进行三维造型，利用ＧＡＭＢＩＴ软件对其进行网格划分．由于半螺旋形吸水室内部流道的结构造型比较复杂，进水流道和叶轮部分均采用适应性较强的非结构四面体单元进行划分，而蜗壳部分则采用混合网格划分，以便得到质量较好的网格条件．经过前期计算校核，最后确定网格总数为 ３２０ ０００左右时，可得到稳定的计算结果． 半螺旋形吸水室流道弯曲变化较大，计算采用 修正的尼一占湍流模型，压力速度耦合方式选用 ＳＩＭＰＬＥ算法，迭代过程采用欠松弛因子进行［２。Ｊ． Ｆｉｇ．２ 图２对称面叶轮蜗壳流道静压分布 Ｓｔａｔｉｃｐｒｅｓｓｕｒｅｉｎ ｄｉｓｍｂｕｔｉｏｎｏｆ ｖｏｌｕｔｅ ８ｙ咖ｅｔｒｙ ｐｌａｉｎ ｉｍｐｅｌｌｅｒ肌ｄ ３边界条件 ３．１进口边界条件 进口采用第一类边界条件，在进口边界上给定初始速度．假定速度方向垂直于进口面厂舻速度： ‰：孚：墨 式中Ｑ为进口体积流量，ｒ为进口处流道半径． 进口湍动能Ｉ｜｝和耗散率占由下列经验公式确定： 图３对称面叶轮蜗壳流道总压分布 ｜｜｝观００５Ｈ：，占：守学 式中Ｚ为进口混合长度；巳取值０．０９． Ｆｉｇ．３ ＴｏｔａｌｐｒｅｓｓｕｒｅｄｉｓｔＴ曲ｕｔｉｏｎｏｆｉｎｉｍｐｅｌｌｅｒ鲫ｄｖｏｌｕｔｅ ｓｙⅡ皿ｅｔｒｙｐｌａｉｎ 由图２可见，在叶轮进口有明显的低压区，并且 万方数据　 第４期朱荣生等：低比转速半螺旋吸水室双吸泵流场的数值分析 ４１ 在叶片吸力面处存在一个最低压区，这与双吸泵最容易发生汽蚀的位置相吻合；沿着叶轮流道，随着半径的增大，静压值逐渐升高，在叶轮出口处压力达到最高．由图３可见，在叶轮流道内动压值随着叶轮半径的增大而增大，由于隔舌处过流面积较小，在隔舌后的蜗壳流道小区域内流速明显增大；在扩散管流道内，靠近叶轮一侧管壁附近的的流体动压，比另一侧要大，流速相对较高．４．２蜗壳内叶轮速度 图４为设计流量点时，在对称面蜗壳和叶轮流道内的绝对速度矢量图．由图４可见，流速从叶片进口到出口逐渐增大，叶轮流道内流速分布不均．在叶轮流道进口到中间位置，靠近叶片压力面的流体流速明显大于吸力面的；在流道中间处，压力面和吸力面附近流体流速差变小；在叶轮出口，受叶轮内部流道流速影响，压力面和吸力面附近流体流速都较高．在蜗壳流道内，流速沿着蜗壳内流道逐渐平稳降低．流体流速在扩散管内也随着管径的增大而逐渐降低．总体速度分布符合蜗壳设计（由于从蜗壳第八断面到扩散管前包括隔舌的一段使用了加密网格，所以在矢量图中，这一段速度矢量显示较密）． 图４对称面蜗壳叶轮流道速度矢量 Ｆｉｇ．４ Ｖｅｌｏｃｉｔｙｖｅｃｔｏｒ ｃｈａｎｏｆ８ｙ姗曲ｒｙ ｐｌａｉｎ ｉｎ ｉｍｐｅｌｌｅｒ加ｄ ｖｏｌｕｔｅ ４．３半螺旋形吸水室流道内的压力和速度 图５为设计流量点时，半螺旋形吸水室进水流道沿进口平面绝对速度矢量图． 图５吸水室流道绝对速度矢量 Ｆｉｇ．５ Ａｂ∞ｌｕｔｅｖｅｌｏｃｉＩｙｖｅｃｔｏｒｃＩｌ昌Ｉｎｏｆｓｕｃｔｉ∞ｃｌｌ丑【Ｉｎｌ）盯 由图５可见，在进水流道进口靠近蜗壳外壁面的分水区域，由于流道弯曲变化太大，受蜗壳外壁面 万　 方数据的阻挡，形成一个小的低速区．在流道中间区域，随着过水面逐渐变小，流速逐渐增大，符合进水流道的设计．图６为设计流量点时，半螺旋形吸水室进水流道垂直出口方向截面的绝对速度矢量图，由图可见流速从流道进口到出口的变化，以及出口处流速沿轴线截面的分布．图６中，由于受密封体挡水板和隔舌（见图１右下图所示）的影响，在设计流量点，进口流道出口挡水板前面流道内出现了不明显的漩涡，但是在整