吉牛水电站冲击式水轮机水斗沙水流动特性及泥沙磨损研究
Study on Internal Flow Characteristics and Sediment Wear of Pelton Turbine Bucket in Jiniu Hydropower Station
Corresponding author: LIU Xiaobing,liuxb@mail.xhu.edu.cn
-
摘要:
冲击式水轮机由于水头高,其高速射流引起的水轮机水斗泥沙磨损非常严重。为探究冲击式水轮机在高速水流下泥沙颗粒对水斗的磨损,基于CFD技术对吉牛水电站冲击式水轮机全流域进行液气固三相流动数值计算和水斗泥沙磨损预估,获得转轮水斗的内部流动特性和泥沙磨损情况。结果表明:水斗根部的出水边、水斗工作面靠近根部处以及切口和分水刃磨损严重,分水刃及切口的磨损大于水斗工作面的磨损,分水刃到出水边的磨损率总体上先上升后下降;水斗磨损的计算结果与实测结果较为一致,计算结果仅比实测结果偏大6%左右,说明泥沙磨损预估方法有效可行。本文的研究结果为冲击式水轮机泥沙磨损预估和电站运行维护提供了技术方法和依据。
Abstract:Due to the high head of the Pelton turbine, the high-speed jet causes very serious sediment wear to the turbine bucket. In order to explore the wear of sediment particles on the bucket of Pelton turbine under high-speed water flow, the liquid-gas-solid three-phase flow numerical calculation and sediment wear prediction of the bucket were carried out based on CFD technology for the whole flow field of Pelton turbine of Jiniu Hydropower Station, and the internal flow characteristics and sediment wear of a runner bucket were obtained. The results show that the trailing edge at the root of the bucket, the leading face of the bucket near the root as well as the notch and the splitter are severely eroded. The wear of the splitter and the notch is greater than that of the working face of the bucket. The wear rate from the splitter to the trailing edge generally increases first and then decreases. The calculated results of bucket wear are in good agreement with the measured results, and the calculated results are only about 6% larger than the measured results, which shows that the method of sediment wear prediction is effective and feasible. This study provides a technical method and basis for sediment wear prediction of Pelton turbine and operation and maintenance of power station.
-
水轮机泥沙磨损是水电机组常见的现象。泥沙磨损常导致水电站水轮机效率下降,过流部件损坏,检修频繁,难以维持正常运行,造成巨大的经济损失。冲击式水电站由于水头高,一旦水轮机处于沙水运行,高速射流对水轮机过流部件的磨损更是非常严重,因此水电工程建设非常重视冲击式水电站的选址和抗泥沙磨损设计。由于近年来水土流失严重,不可避免地造成一些冲击式水电站在浑水中运行,因此冲击式水轮机泥沙磨损研究引起了国内外学者的高度重视。
Guo等[1-2]采用Eulerian-Lagrangian方法分析了冲击式水轮机水斗内部气液固流动特性,发现涡旋脱落引起的二次流,颗粒分离增加,从而增强了脱落区域的磨损。Han等[3-6]研究发现在喷嘴的收缩段之后,颗粒更多地集中在喷嘴的内侧,并分析环管及喷嘴处的涡旋结构和二次流动,还发现颗粒干扰水的分布,从而降低了总扭矩,水力效率降低了9%左右。肖业祥等[7]提出一种新的算法预测了粒子的冲击行为,对水斗的磨损进行了预估,与实测的磨损分布较为一致。葛新峰等[8-9]发现速度的增加会加大水斗的磨损率,而喷嘴开度的变化对磨损区域和最大磨损率的影响不大,喷射机构受到的最大磨损率随泥沙浓度的增加而增大。Tarodiya等[10]发现喷嘴喷针的磨损存在不对称现象,不同开度下喷嘴的磨损分布相似。Messa等[11]研究结果表明喷针顶角的减小可能会增加磨损的风险。Thakur等[12]通过改变泥沙参数及射流速度进行多工况的冲击式水轮机磨损试验,获得了磨损率公式。Din等[13]通过测量喷针和喷嘴在磨损过程中的材料去除量,得到了喷针和喷嘴分别磨损了3.71%和5%。Padhy等[14-15]通过显微镜观察磨损实验前后的试样,分析磨损机制,研究发现泥沙尺寸是产生磨损的一个重要参数。Pachón等[16]发现经过等离子渗氮工艺或抗磨材料喷涂处理后的喷嘴耐磨性显著提高。Abgottspon等[17]用3D光学扫描仪、超声波测厚仪测量出电站冲击式水轮机磨蚀厚度。Rai等[18-19]将电站水斗按比例缩小,进行不同泥沙大小、浓度等工况的实验,开发了预测水斗磨损的磨损模型。
由于冲击式水轮机速度高、内部流动复杂,目前水轮机泥沙磨损研究成果并不多,尤其是水轮机过流部件表面泥沙磨损分布情况。本研究针对吉牛水电站冲击式水轮机建立了全流域几何物理模型,基于CFD技术进行水轮机内部沙水(含气)流动计算,并采用磨损模型预估水斗的磨损情况。
1. 吉牛水电站工程及运行概况
吉牛水电站位于四川甘孜州丹巴县革什扎河干流上,是革什扎河“一库四级”水电开发方案的最后一级,为低闸引水式电站。水库正常蓄水位2 378.00 m,总库容197.5万m3,具有日调节性能。引水线路全长22.42 km,额定引用流量60.28 m3/s。电站装机2台,总装机容量240 MW,年利用小时4929 h,多年平均发电量11.829亿kW·h。电站在系统负荷低谷期安排电站停机冲沙,冲沙历时6 h。当入库日平均流量大于190 m3/s时,电站停机避峰,全闸打开,停机冲沙。在汛期出现流量不大但含沙量达到一定量时,为了避免其对水轮机的磨损,也将停机冲沙。据电站运行数据分析,水轮机多年在暴雨期运行的平均最大含沙量约0.212 kg/m3,累计泥沙过机时间每年不足15天,其他时间过机泥沙量很小。
革什扎河流泥沙特性参数见表1。沙样矿物成分中,莫氏硬度大于等于5的矿物有石英、长石、千枚岩屑、火山岩屑、绿帘石、磷灰石等,其含量约占81.5%。
表 1 水文泥沙参数Table 1. Hydrological sediment parameters参数 数值 泥沙中值粒径/mm 0.1 泥沙密度/(kg/m3) 2650 多年暴雨期过机平均最大含沙量/(kg/m3) 0.212 2. 数学模型
2.1 多相流模型
冲击式水轮机内部沙水流动实际上是液气固三相流动。VOF模型可以通过求解单一的动量方程,并跟踪区域内每个流体的体积分数来模拟两种或两种以上的非混溶流体,这里用它来捕捉水气交界面。VOF模型要求在每个控制体中,所有相的体积分数之和为1,因此在本模拟中,要保证液气固三相的体积分数之和始终为1。这里对液气固三相流体的物理量进行线性平均:
$$ \rho=\alpha_i \rho_i $$ (1) $$ \mu=\alpha_i \mu_i $$ (2) 式中:α为体积分数;ρ为密度;μ为动力黏性系数;下标i为张量坐标(表示液气固各相)。
2.2 湍流模型
SST k-ω模型在近壁面采用k-ω模型,而在边界层外采用k-ε模型,它包含了修正的湍流黏性公式,考虑了湍流剪切应力的效应,适用于冲击式水轮机的数值模拟。SST k-ω模型的运输方程如下:
$$ \frac{\partial }{{\partial t}}\left( {\rho k} \right) + \frac{\partial }{{\partial {x_i}}}\left( {\rho k{u_i}} \right) = \frac{\partial }{{\partial {x_j}}}\left( {{{{\varGamma}} _k}\frac{{\partial k}}{{\partial {x_j}}}} \right) + {G_k} - {Y_k} + {S_k} $$ (3) $$ \frac{\partial }{{\partial t}}\left( {\rho \omega } \right) + \frac{\partial }{{\partial {x_i}}}\left( {\rho \omega {u_i}} \right) = \frac{\partial }{{\partial {x_j}}}\left( {{\varGamma _\omega }\frac{{\partial \omega }}{{\partial {x_j}}}} \right) + {G_\omega } - {Y_\omega } + {D_\omega } + {S_\omega } $$ (4) 式中:t为时间;u为流体速度;x为坐标;Gk为湍动能k产生项;Gω为ω产生项;Yk和Yω分别为k和ω的有效扩散项;Dω为正交发散项;Sω和Sk为用户自定义源项;下标i、j为张量坐标。其中:
${\varGamma _k} = \mu + \dfrac{{{\mu _t}}}{{{\sigma _k}}}$ ,${\varGamma _\omega } = \mu + \dfrac{{{\mu _t}}}{{{\sigma _\omega }}}$ ,μt为湍涡流动力黏性系数,σk和σω分别为k和ω的湍流普朗特数,${\sigma _k} = \dfrac{1}{{0.85{F_1} + (1 - {F_1})}}$ ,$ {\sigma _\omega } = \dfrac{1}{{0.5{F_1} + 0.856(1 - {F_1})}} $ ,F1为壁面函数值;$ {G_k} = {\mu _t}{S^2} $ ,$ S = \sqrt {2{S_{ij}}{S_{ij}}} $ ,${S_{ij}} = \dfrac{1}{2}\left( {\dfrac{{\partial {u_j}}}{{\partial {x_i}}} + \dfrac{{\partial {u_i}}}{{\partial x_j}}} \right)$ ;$ {G_\omega } = \dfrac{{\rho {G_k}}}{{{\mu _t}}} $ ;Yk=ρβ*kω,β*=0.09;Yω=ρβkω2,β=0.075F1+0.0828(1−F1);$ {D_\omega } = 2(1 - {F_1})\rho {\sigma _{\omega ,2}}\dfrac{1}{\omega }\dfrac{{\partial k}}{{\partial {x_j}}}\dfrac{{\partial \omega }}{{\partial {x_j}}} $ ,σω,2=1.168。2.3 颗粒轨迹模型
多相流中的离散相颗粒采用颗粒运动拉格朗日方程模型[20]来追踪颗粒轨迹。本模型适用于平均颗粒相体积分数小于0.001(含沙量2.5 kg/m3左右)的情况。离散项颗粒的控制方程如下:
$$\begin{split} \frac{{{\text{d}}{u_{pi}}}}{{{\text{d}}t}} =& \frac{1}{{{K_m} + \bar \rho }}\left[ {\frac{3}{{4{d_p}}}{C_D}\left| {\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\rightharpoonup}$}} {u} - {{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\rightharpoonup}$}} {u} }_p}} \right|\left( {{u_i} - {u_{pi}}} \right)}+ \right.\\& \frac{3}{{2{d_p}}}{K_B}\sqrt {\frac{\nu }{{\text{π }}}} \int_{ - \infty }^t {\left( {\frac{{{\text{d}}{u_i}}}{{{\text{d}}\tau }} - \frac{{{\text{d}}{u_{pi}}}}{{{\text{d}}\tau }}} \right)\frac{{{\text{d}}\tau }}{{\sqrt {t - \tau } }}}+\\& \frac{6}{{{\text{π }}{d_p}}}{K_S}{\left| {\nu \frac{{\partial {u_j}}}{{\partial {x_i}}}} \right|^{\frac{1}{2}}}\left( {{u_j} - {u_{pj}}} \right){{\rm{sgn}}} \left( {\frac{{\partial {u_j}}}{{\partial {x_i}}}} \right) +\\& \frac{3}{4}{C_M}{\Omega _i} \cdot \left( {{u_i} - {u_{pi}}} \right) + {K_m}\frac{{{\text{d}}{u_i}}}{{{\text{d}}t}} \left. { - \frac{1}{\rho }\frac{{\partial P}}{{\partial {x_i}}} + \bar \rho {g_i}} \right] \end{split} $$ (5) 式中:up为颗粒速度;Km为虚拟质量力系数;
$ \bar \rho $ 为颗粒密度ρp与流体密度ρ的比值;CD为颗粒阻力系数;dp为颗粒粒径;KB为Basset力系数;ν为流体运动黏性系数;KS为Saffman升力系数;CM为Magnus升力系数;${{{\varOmega }}_i} = {\omega _{pi}} - 0.5\nabla \cdot {u_i}$ ,ωp为颗粒自身旋转角速度;P为压力;g为重力加速度;sgn为符号函数。2.4 磨损模型
含沙水射流撞击水斗,对转轮水斗表面造成磨损,采用Generic模型预估转轮的磨损分布及磨损率变化,表达式如下:
$$ {R_e} = \sum\limits_{p = 1}^{{N_P}} {\frac{{{m_p}c\left( {{d_p}} \right)f\left( \alpha \right){w^{b\left( w \right)}}}}{{{A_f}}}} $$ (6) 式中:Re为磨损率;NP为颗粒总数;mp为颗粒质量流率;c(dp)为粒径函数;f(α)为冲角函数(α为颗粒冲击壁面的角度);b(w)为颗粒相对速度的函数(w为颗粒与壁面的相对速度);Af为壁面面积。
3. 几何模型及计算设置
3.1 过流部件三维几何模型建立
根据吉牛水电站冲击式水轮机设计等资料,建立水轮机过流部件计算几何模型。水轮机基本设计参数见表2,其中过流部件金属材质如下:喷嘴材质为00Cr13Ni5Mo,喷针材质为00Cr13Ni5Mo,转轮整体材质为X3CrNi13-4。转轮如图1所示。
表 2 水轮机基本设计参数Table 2. Design parameters of hydraulic turbine参数 数值 最大水头/m 506.5 额定水头/m 457.0 最小水头/m 456.5 额定流量/(m³/s) 30.14 额定转速/(r/min) 300 额定出力/MW 123 额定效率/% 91.79 水斗数 21 喷嘴数 6 转轮节圆直径/mm 2890 所建水轮机计算全流域三维几何模型如图2所示,主要由以下4部分组成:1)包含水斗的旋转域,2)围着旋转域外一圈的静止域,3)喷嘴出口的射流域,4)包含喷嘴喷针的射流机构。
3.2 几何模型网格划分
冲击式水轮机的结构较为复杂,对射流机构及射流域采用结构化网格,而对旋转域采用非结构化网格,并对水斗表面网格进行了加密处理,准确地捕捉水斗表面的流态,图3为网格划分示意图。由于网格对数值模拟的结果影响较大,对全流域的网格进行了无关性验证,见表3。网格数量增加,预测效率值逐渐逼近设计效率值,综合考虑计算资源及数值模拟的准确性,最终选择方案2,全流域模型网格数约为1254万个。
表 3 网格无关性验证Table 3. Grid independence verification方案 网格数/万 预测效率/% 设计效率/% 相对误差/% 1 677 88.27 91.79 3.52 2 1254 91.48 91.79 0.31 3 1843 91.67 91.79 0.12 3.3 计算设置
计算方法采用SIMPLEC算法,时间步长设置为1.11×10–4 s,每个时间步迭代20次。将水和空气视为连续相,加入离散相颗粒进行液气固三相流计算,离散相采用双向耦合,考虑离散相与连续相的相互作用,每10步连续相迭代进行一次轨迹计算。
采用速度进口,流速为7.51 m/s,出口设置为压力出口,压力为0 。旋转域与静止域采用交界面连接,固定壁面采用无滑移边界条件,用标准壁面函数来模拟近壁面区域的流动。离散相颗粒简化为球形颗粒,垂直于进口面射入,入射速度与水流速度保持一致,与壁面的接触方式均设置为反弹。根据吉牛电站提供的水文泥沙资料,泥沙参数为:泥沙中值粒径dp=0.1 mm,多年暴雨期过机平均最大含沙量0.212 kg/m3,额定工况过机平均最大泥沙质量流量6.39 kg/s,泥沙颗粒密度2 650 kg/m3。
4. 数值计算及结果分析
4.1 转轮水斗内部流动特性
水流经射流机构,从喷嘴出口射流到大气与空气结合。图4为全流域平面的水气体积分布图,红色为水相,蓝色则为气相,喷嘴内部完全充满了水,在流出喷嘴出口后,水相与射流域、静止域以及旋转域中的气相形成了一层薄的空气 – 水界面。由于各水斗的扭矩都具有周期性,现以单个水斗为研究对象,在计算过程中监测了扭矩随时间的变化情况,并提取了一时间段的扭矩 – 时间曲线,如图5所示。在该周期内选取t1(扭矩上升阶段)、t2(高扭矩阶段)、t3(扭矩下降阶段)3个时间点,对应的水相分布情况在图5中标明(图上到下分别为1#水斗、2#水斗、3#水斗),具体分析这3个时间点的流动特性。
4.1.1 水斗工作面速度分布
水斗工作面的速度及流线分布如图6所示。t1时刻,射流只冲击到分水刃头部,该处流速较大,流线方向与分水刃呈锐角,出水边处流速较小,分水刃一侧的流线与该侧靠近出水边处的流线均朝着一个方向汇集。t2时刻,部分射流较为垂直地冲击水斗基圆处,该处的流线方向与分水刃呈直角,分水刃处速度最大,沿流线向两侧出水边呈减小的趋势,水流朝出水边方向充分延展。t3时刻,仍有部分流体在水斗内,速度沿分水刃向两侧出水边减小,在水斗工作面靠近出水边处达到了最小值,单侧的流线朝着一个方向汇集。
4.1.2 水斗工作面压力分布
水斗工作面的压力分布情况如图7所示。t1时刻,射流部分被2#水斗的分水刃头部拦截,造成了局部的高压区,在分水刃两侧也有小部分射流产生的压力区,随着射流越来越多地被2#水斗拦截,2#水斗的扭矩持续上升。t2时刻,射流被1#水斗部分拦截,部分射流朝向2#水斗基圆附近冲击,此刻2#水斗扭矩达到该周期内的最大值,此时的高压区域对称分布于分水刃两侧,且高压区域在水斗底处往水斗根部偏移,压力值从水斗最底处至出水边逐渐减小。t3时刻,已无射流冲击2#水斗,扭矩迅速下降,水斗工作面的水膜延展开来,此时高压力区域往出水边移动。
4.1.3 水斗工作面泥沙浓度分布
对t1、t2以及t3时刻下的水斗工作面泥沙浓度分布进行分析,如图8所示。t1时刻,分水刃头部泥沙浓度较高。t2时刻,分水刃中部泥沙浓度较高,分水刃两侧泥沙分布范围较广,泥沙较为对称地分布于分水刃两侧。t3时刻,分水刃处无明显泥沙分布,泥沙分布范围由分水刃向出水边逐渐扩大。取3个水斗基圆截面观察泥沙浓度分布,如图9所示,流体中携带的泥沙颗粒沿分水刃向两侧出水边分布,同时2#水斗表面泥沙含量及浓度均比1#水斗及3#水斗大。此时,1#水斗虽拦截了部分射流,但射流未冲击到基圆处,1#水斗分水刃处泥沙浓度较低,泥沙主要位于分水刃附近两侧。与此同时,3#水斗中由于之前射流进水斗的流体还未完全排出,因此水斗底处及出水边仍聚集了较多泥沙颗粒,且离工作面壁越近泥沙浓度越高。
4.2 转轮水斗磨损情况
4.2.1 水斗工作面磨损率分布
水斗受到高速射流的冲击会造成较为严重的磨损,计算得到的水斗磨损率云图如图10所示。从磨损率的分布可以预测磨损较为严重的部位位于靠近水斗根部的出水边、水斗工作面靠近根部处、切口以及分水刃。
为便于分析磨损率与磨损相对位置的关系,以分水刃为中心(即横坐标的0点),−1~0代表左半边截线弦长,0~1代表右半边截线弦长,横坐标定义为该点在基圆截线(基圆截线示意图如图11所示)上的相对位置,绘制出如图12所示的关系曲线。从图12中可以观察到,分水刃处的磨损最为严重,水斗分水刃两边的磨损分布较为对称,磨损率总体上都有先上升后下降的趋势,最大值略有变化,出水边处存在一定磨损。
4.2.2 水斗工作面磨损量的预估与实测
选取分水刃、切口及水斗工作面靠近根部部位进行泥沙磨损分析。根据吉牛水电站多年水文泥沙资料以及运行情况,计算水轮机运行7年后转轮水斗的泥沙磨损量,计算结果见表4。从表4的数据可以发现分水刃及切口处的磨损大于水斗工作面的磨损,分水刃的磨损大于切口处,水斗工作面靠近根部部位的磨损也较严重。
表 4 水斗不同部位最大磨损率及磨损量Table 4. Maximum wear rate of different parts of the bucket部位 分水刃 切口 水斗工作面 计算磨损率/(mm/s) 5.20×10−7 5.08×10−7 4.19×10−7 计算磨损量/mm 4.31 4.21 3.47 实测磨损量/mm 4.04 3.94 3.28 对吉牛水电站运行7年后的实际磨损情况进行现场实测和拍照,水斗磨损量的实测结果见表4。磨损较为严重的区域(即分水刃与切口)位于图13中椭圆圈示意处,水斗工作面磨损分布的计算结果与实测结果较一致。同时,现场观察还发现有个别区域出现较大的刮痕、冲击坑,可能是电站运行时,个别较大的高硬度泥沙颗粒或石块冲击所致。另外,水斗表面还存在一些“蜂窝结构”磨蚀坑,应是磨损与空化的联合作用加剧了过流部件的磨蚀破坏。
从表4的数据可以发现,水轮机水斗磨损计算结果比实测结果偏大6%左右。我们认为这个偏差应该与电站获取的水文泥沙和运行工况资料的准确度以及计算的磨损模型精度有关,不过这个水轮机泥沙磨损预估偏差应该在工程的允许范围内,计算结果已能较准确地反映出水轮机水斗的磨损情况。
5. 结论
针对吉牛水电站冲击式水轮机建立全流域模型,采用SST k-ω湍流模型、VOF多相流模型和颗粒动力学模型进行液气固三相流动模拟,并采用Generic模型预估转轮水斗的磨损分布及磨损率变化,得到以下主要结论。
1)水轮机水斗根部的出水边易受到磨损,水斗工作面靠近根部处、切口以及分水刃磨损严重,分水刃到出水边的磨损率总体上先上升后下降。
2)计算所得水斗工作面最大磨损率为4.19×10–7 mm/s,切口最大磨损率为5.08×10–7 mm/s,分水刃最大磨损率为5.20×10–7 mm/s,分水刃及切口的磨损大于水斗工作面的磨损。
3)水轮机水斗磨损量的计算结果比实测结果偏大6%左右,计算的水轮机磨损率和磨损量分布与电站水轮机实测磨损分布较为一致,表明其泥沙磨损预估方法是可行的。尽管吉牛电站的水质一般情况下很好,仅在暴雨期存在短暂含沙量较大的情况,也非常注意避沙峰运行,但研究发现水轮机的磨损还是很厉害,因此冲击式水轮机要高度重视抗泥沙设计和避沙峰运行。本研究为冲击式水轮机泥沙磨损预估和电站运行维护提供了技术方法和依据。
-
表 1 水文泥沙参数
Table 1 Hydrological sediment parameters
参数 数值 泥沙中值粒径/mm 0.1 泥沙密度/(kg/m3) 2650 多年暴雨期过机平均最大含沙量/(kg/m3) 0.212 表 2 水轮机基本设计参数
Table 2 Design parameters of hydraulic turbine
参数 数值 最大水头/m 506.5 额定水头/m 457.0 最小水头/m 456.5 额定流量/(m³/s) 30.14 额定转速/(r/min) 300 额定出力/MW 123 额定效率/% 91.79 水斗数 21 喷嘴数 6 转轮节圆直径/mm 2890 表 3 网格无关性验证
Table 3 Grid independence verification
方案 网格数/万 预测效率/% 设计效率/% 相对误差/% 1 677 88.27 91.79 3.52 2 1254 91.48 91.79 0.31 3 1843 91.67 91.79 0.12 表 4 水斗不同部位最大磨损率及磨损量
Table 4 Maximum wear rate of different parts of the bucket
部位 分水刃 切口 水斗工作面 计算磨损率/(mm/s) 5.20×10−7 5.08×10−7 4.19×10−7 计算磨损量/mm 4.31 4.21 3.47 实测磨损量/mm 4.04 3.94 3.28 -
[1] GUO B, XIAO Y, RAI A K, et al. Analysis of the air-water-sediment flow behavior in Pelton buckets using a Eulerian-Lagrangian approach[J]. Energy, 2021, 218: 119522. doi: 10.1016/j.energy.2020.119522
[2] GUO B, XIAO Y, RAI A K, et al. Sediment-laden flow and erosion modeling in a Pelton turbine injector[J]. Renewable Energy, 2020, 162: 30 − 42. doi: 10.1016/j.renene.2020.08.032
[3] HAN L, WANG Y, ZHANG G F, et al. The particle induced energy loss mechanism of Pelton turbine[J]. Renewable Energy, 2021, 173: 237 − 248. doi: 10.1016/j.renene.2021.03.136
[4] HAN L, ZHANG G F, WANG Y, et al. Investigation of erosion influence in distribution system and nozzle structure of pelton turbine[J]. Renewable Energy, 2021, 178: 1119 − 1128. doi: 10.1016/j.renene.2021.06.056
[5] HAN L, WANG Y, ZHANG G F, et al. Investigation of the erosion characteristics in the distribution parts of pelton turbine[J]. Journal of Applied Fluid Mechanics, 1999, 14(6): 1809 − 1816.
[6] HAN L, DUAN X L, GONG R Z, et al. Physic of secondary flow phenomenon in distributor and bifurcation pipe of Pelton turbine[J]. Renewable Energy, 2019, 131: 159 − 167. doi: 10.1016/j.renene.2018.06.116
[7] 肖业祥, 郭豹, 张瑾, 等. 冲击式水轮机含沙三相流与磨损特性数值研究[J]. 水电与抽水蓄能, 2021, 7(2): 4 − 10. XIAO Y X, GUO B, ZHANG J, et al. Numerical analysis of the three-phase flow and sand abrasion of a pelton turbine[J]. Hydropower and Pumped Storage, 2021, 7(2): 4 − 10.
[8] GE X F, SUN J, ZHOU Y, et al. Experimental and numerical studies on opening and velocity influence on sediment erosion of pelton turbine buckets[J]. Renewable Energy, 2021, 173: 1040 − 1056. doi: 10.1016/j.renene.2021.04.072
[9] 葛新峰, 孙洁, 李阳, 等. 冲击式水轮机喷射机构泥沙磨损特性的数值模拟[J]. 水利学报, 2020, 51(12): 1486 − 1494. doi: 10.13243/j.cnki.slxb.20200430 GE X F, SUN J, LI Y, et al. Numerical simulation of silt erosion characteristics of an injector of Pelton turbine[J]. Journal of Hydraulic Engineering, 2020, 51(12): 1486 − 1494. doi: 10.13243/j.cnki.slxb.20200430
[10] TARODIYA R, KHULLAR S, LEVY A. Particulate flow and erosion modeling of a Pelton turbine injector using CFD-DEM simulations[J]. Powder Technology, 2022, 399: 117168. doi: 10.1016/j.powtec.2022.117168
[11] MESSA G V, MANDELLI S, MALAVASI S. Hydro- abrasive erosion in Pelton turbine injectors: A numerical study[J]. Renewable Energy, 2019, 130: 474 − 488. doi: 10.1016/j.renene.2018.06.064
[12] THAKUR R, KUMAR A, KHURANA S, et al. Correlation development for erosive erosion rate on pelton turbine buckets[J]. Int J Mech Prod Eng Res Dev, 2017, 7: 259 − 274.
[13] DIN M Z U, HARMAIN G A. Assessment of erosive wear of Pelton turbine injector: Nozzle and spear combination–A study of Chenani hydro-power plant[J]. Engineering Failure Analysis, 2020, 116: 104695. doi: 10.1016/j.engfailanal.2020.104695
[14] PADHY M K, SAINI R P. Study of silt erosion mechanism in Pelton turbine buckets[J]. Energy, 2012, 39(1): 286 − 293. doi: 10.1016/j.energy.2012.01.015
[15] PADHY M K, THATOI D N, ACHARYA A K. Effect of shape of silt particles on erosive erosion of pelton turbine bucket[C]//IEEE-International Conference on Advances in Engineering, Science and Management (ICAESM-2012). [S.l.]: IEEE, 2012: 19 − 24.
[16] PACHÓN I F, LOBOGUERRERO J, MEDINA J A. Development of a test rig to evaluate abrasive erosion on pelton turbine nozzles. A case study of Chivor hydropower[J]. Erosion, 2017, 372: 208 − 215.
[17] ABGOTTSPON A, STAUBLI T, FELIX D, et al. Hydro-abrasive erosion of Pelton buckets and suspended sediment monitoring[C]// Proceedings of the Hydro Vision International. Denver, USA: [s.n.] , 2013.
[18] RAI A K, KUMAR A, STAUBLI T. Effect of concentration and size of sediments on hydro-abrasive erosion of Pelton turbine[J]. Renewable Energy, 2020, 145: 893 − 902. doi: 10.1016/j.renene.2019.06.012
[19] RAI A K, KUMAR A, STAUBLI T. Analytical modelling and mechanism of hydro-abrasive erosion in Pelton buckets[J]. Wear, 2019, 436/437: 203003. doi: 10.1016/j.wear.2019.203003
[20] 刘小兵, 曾永忠, 华红, 等. 水轮机沙水流动及磨损[M]. 北京 : 中国水利水电出版社, 2020. LIU X B, ZENG Y Z, HUA H, et al. Sand flow and wear of turbine [M]. Beijing: China Water Power Press, 2020.
-
期刊类型引用(2)
1. 詹怡阳,汪勇杰,栗正新. 金刚石耐磨涂层在水轮机叶片上的应用研究进展. 超硬材料工程. 2024(04): 25-31 . 百度学术
2. 范文睿,孙龙刚,郭鹏程,李琛玺,孙红. 冲击式水轮机内部流动与泥沙磨损特性研究进展. 水动力学研究与进展A辑. 2023(05): 774-793 . 百度学术
其他类型引用(4)