基于小波变换的核主泵惰转过程叶轮压力脉动及泵内部流动规律.pdf

书书书 年 月 第 卷 第 期 刘安林基于小波变换的核主泵惰转过程叶轮压力脉动及泵内部流动规律刘安林，叶道星，罗逸民，陈俊霖，余波，赖喜德，吴锦琛，翟凤丽（西华大学能源与动力工程学院，四川 成都 ；西华大学流体及动力机械教育部重点实验室，四川 成都 ；西华大学流体机械及工程四川省重点实验室，四川 成都 ）收稿日期：；修回日期：；网络出版时间：网络出版地址：基金项目：四川省科技厅科技计划项目（）；四川省“双一流”学科建设项目（）；德阳市开放式校市合作技术研发资助项目（）；西华大学研究生创新基金资助项目（，）第一作者简介：刘安林（），男，重庆永川人，硕士研究生（），主要从事流体机械内部流动机理研究通信作者简介：叶道星（），男，四川成都人，副教授，博士（），主要从事流体机械内部流动机理及设计方法研究摘要：为研究惰转过程核主泵内部的流动规律，应用数值模拟的方法对核主泵内部流场进行计算，得到不同惰转流量时间（）下叶轮流道内压力脉动的时频特性 对叶轮压力脉动进行小波变换分析，并采用第二代涡识别方法 准则对不同工况下各过流部件进行涡识别 研究结果表明：在惰转过程中，惰转初期压力脉动系数最大变化发生在导叶中部，随后最大变化率发生在叶轮与导叶的动静交界面，最后压力脉动变化趋于稳定；随着叶轮转速和流量逐渐减小，导叶做功区域向出口附近转移；隔舌部位由于结构特殊所受冲击较大，附近有回流产生，并有向流道中部扩散的迹象；压力脉动趋于稳定值，主要频率向低频转移，各部件流道内压力脉动规律性逐渐被破坏，其中导叶流道规律性变化趋势较为平缓；流道内的部分强涡区域呈现分离脱落迹象，并且数量逐渐减少；压水室内的涡核主要集中于左侧，与叶轮旋转方向有关 研究可为核主泵在惰转特性优化设计提供一定理论支撑关键词：核主泵；惰转过程；压力脉动；小波变换；准则；数值模拟中图分类号：；文献标志码：文章编号：（）：刘安林，叶道星，罗逸民，等 基于小波变换的核主泵惰转过程叶轮压力脉动及泵内部流动规律 排灌机械工程学报，（）：，（），（）：，（），（，；，；，）：，第 期刘安林，等基于小波变换的核主泵惰转过程叶轮压力脉动及泵内部流动规律 （），：；核主泵作为核电站一回路主要承压边界和唯一的高速旋转设备，被誉为堆芯中高温高压放射性介质环境的“心脏”，因此核主泵的稳定运行对堆芯安全性具有至关重要的作用 由于核主泵压水室等非常规设计，流动不稳定可导致严重振动，此外，反向流动、叶片后缘涡脱落等非定常流动也会引起强烈的脉动 ，其中核主泵在惰转下的压力脉动情况更为复杂，对于核主泵在惰转工况下的运行影响较大 因此，研究核主泵在惰转下的压力脉动特性，对核电站的安全稳定运行具有重要工程价值近年来国内外对核主泵压力脉动进行了大量研究 赖喜德等 、等 分别对核主泵在额定工况下的压力脉动进行了较为全面的研究，结果表明，在额定工况下，核主泵高幅值压力脉动主要出现于叶轮吸力面、导叶喉部和环形压出室左侧区域朱荣生等 、等 研究了小流量工况下的核主泵压力脉动特性，结果表明，核主泵在小流量工况下的压力脉动明显大于额定工况，相比于叶轮其他区域，叶轮叶片背面进口处压力脉动更加复杂 倪丹等 、王鹏等 研究了核主泵在大流量工况下的压力脉动特性，结果表明，在大流量工况下，叶轮叶片出口的尾迹射流导致的压力脉动增强，环形压出室右侧压力脉动非稳态性较环形压出室左侧更强断电事故工况下核主泵的惰转过渡属于瞬态流动过程，其内部流动的瞬态模拟分析相比稳态过程数值模拟分析更为紊乱，使得压力脉动特性变得复杂 等 研究了不同惰转条件下核主泵内部压力波动响应规律，然而，对于惰转条件下的压力脉动，研究还有待深入在传统的压力脉动研究过程中，一般使用傅里叶变换进行分析 但核主泵惰转过程中的压力脉动信号为非平稳信号，运用传统的傅里叶变换进行分析具有较大的局限性 小波分析能更好地对压力脉动时频特性的细节进行展示，特别是对于短时间内发生的非平稳压力脉动的分析具有较好的时频分辨率 等 对水泵水轮机从抽水模式快速过渡到发电模式期间的压力脉动进行了小波分析 等 对广义 小波进行了研究，结果表明，广义 小波本质上包含了所有其他常用的分析小波，为小波在应用中的适用性提供了框架为了研究惰转工况下核主泵内部的流动规律，文中应用数值模拟方法对核主泵内部流场进行计算，得到不同惰转流量时间下叶轮流道中压力脉动的时频特性，并对叶轮压力脉动进行小波变换分析，采用第二代涡识别方法 准则对不同工况下各过流部件进行涡识别，以期为实际核主泵安全稳定运行提供技术支撑 计算模型及数值计算方法 计算模型根据某改进压水堆核主泵，采用 软件进行等比例三维建模 为获得稳定的计算结果，对核主泵进出口延长，延长长度大于 倍进出口直径 考虑核 排灌机械工程学报第 卷主泵运行环境苛刻，设计其环形压水室为类球形 图为核主泵计算流体域三维模型，包括进口段、叶轮、导叶、环形压水室和出口段图 核主泵三维模型 核主泵主要设计参数分别为流量 ，扬程 ，转速 ，运行温度 ，出口压力 ，叶轮叶片数 ，导叶叶片数 全部参数可参考文献 惰转工况下数值计算方法 网格划分及其无关性验证采用 对整个计算流体域进行六面体结构化网格划分 对流体域近壁面、叶片进出口和环形压出室隔舌区域进行网格加密处理 网格密度对计算结果影响较大，需进行无关性验证 经过数值计算与试验数据对比，在网格数为 万时，网格的增多对扬程和效率的影响变小 综合考虑计算效率和计算机配置，最终确定计算网格数为 万 计算边界条件及惰转模型设置应用 商用软件平台进行数值计算 总计算时长为 ，核主泵在 时从正常运行工况进入惰转工况 利用文献 中的边界条件和数学模型进行计算设置 在文献 中，半流量（即流量减少到额定流量的 ）时间点为 时，综合考虑选取时间在 ，时的 个工况点（即流量分别为减小到 ，）对核主泵惰转工况下的压力脉动进行分析 为方便描述，引入惰转流量时间 的概念，即惰转过程中百分比流量减小到某一值所处的时间，如 为核主泵的惰转流量减小到设计流量的 时所处的时间即 ，则 个工况点的惰转流量时间分别为 ，和 采用 湍流模型 设置变时间步长进行计算，每旋转 为一个时间步长，时间步长随转速减小而增大 利用各工况的初始流量和初始转速作为初始条件进行定常计算，之后采用定常计算结果作为初始条件进行瞬态非定常计算，共计算 个周期，取最后 个周期数据进行分析，共 个时间步长 监测点设置为了监测不同工况下各流道的压力脉动情况，沿冷却剂流动方向，在叶轮流道设置监测点 ，其中监测点 位于叶轮进口，监测点 位于叶轮叶片流道进口，监测点 位于叶轮叶片流道中部，监测点 位于叶轮出口，如图 所示图 监测点设置 数值计算方法验证为验证数值计算方法的可靠性，在定常状态下对核主泵进行外特性试验 图 给出了不同流量工况下核主泵水力性能的数值计算结果和试验结果的比较，可以看出：计算结果与试验结果的泵水力性能曲线具有相同的变化规律；最大扬程误差发生在 处，为 ，最大效率误差发生在 处，为 ，均在工程允许误差内，这表明文中所采用的数值计算方法是可靠的图 外特性数值计算结果和试验结果对比 数据处理为更好地分析压力脉动情况，引入量纲一的压力脉动系数 ，即 ，（）式中：为监测静压与平均压力的差值；为叶轮出口圆周速度，为叶轮出口直径 第 期刘安林，等基于小波变换的核主泵惰转过程叶轮压力脉动及泵内部流动规律将后处理得到的瞬时压力信号转换为压力脉动系数，然后通过一维连续小波变换对转换后的压力脉动信号进行处理，得到压力脉动的时频特性 一维连续小波的变换原理 为（，）（）（；，），（）其中小波基函数 通过标准小波基函数平移伸缩得到，其表达式为（；，）槡()，（）式中：为伸缩因子；为平移因子文中选用 小波，小波波次为 小波表达式为（）槡 （）（）准则将第二矩阵不变量 大于 的区域定义为涡，同时要求涡区域的压强小于周围流体区域的压强，其表达式为（），（）式中：为矩阵的 范数；与 分别为速度梯度张量 的对称部分和反对称部分，表征流场中的变形与旋转由 准则的定义可以看出：旋涡结构中不仅存在反对称张量 ，而且反对称张量 能够克服对称张量 的抵消作用；当 时，流场中流体的旋转对于变形而言起主导作用由于 个工况下时间周期不同，因此，对时间周期进行量纲一化处理，即 ，其中 为叶轮旋转一周的时间步数，为数据所处时间步位置 计算结果分析 叶轮流道压力脉动分析图 为在惰转过程中叶轮流道监测点 ，各时间点的压力脉动时域特性 由图 可以看出：在惰转初期 时，监测点 和 的压力脉动系数变化规律相似，且值相近；监测点 的压力脉动系数波动较叶轮流道其他监测点更为剧烈；随着惰转过程的进行，叶轮流道内的压力脉动系数逐渐减小；监测点 处各阶段压力系数平均变化分别为 ，；监测点 处各阶段压力系数平均变化分别为 ，；监测点 处各阶段压力系数平均变化分别为 ，；监测点 处各阶 段 压 力 系 数 平 均 变 化 分 别 为 ，图 叶轮流道监测点工况点的压力脉动时域特性 图 为在惰转过程中叶轮流道监测点 ，各时间点的压力脉动小波谱，图中 排灌机械工程学报第 卷为小波强度，用以反映叶轮不同位置随惰转过程的内部能量变化图 叶轮流道监测点的压力脉动连续小波响应 由图 可以看出：由于叶轮转速随惰转过程的发展而减小，其叶轮内的通过频率 会随转速减小而减小，个工况点的叶轮初始通频分别为 ，；相同位置在不同惰转流量时间工况下，高频高小波强度区域占比逐渐减小，低频区域小波强度响应占比逐渐增大；在 工况下，监测点 ，的压力脉动强度明显强于监测点 的压力脉动强度，压力脉动系数沿叶轮流道先增大后减小；在监测点 的 频段和监测点 的 频段出现了一高一低 个信号（见图 和 箭头所指），此时区域处于 附近；叶轮叶片流道 ，在 （大于 ）的高频段响应有 个较为明显的小波波峰和波谷；在 工况下，压力脉动在监测点 的高频小波响应频段和 的高频小波响应频段近似，主要分布频段为 ，即 ，明显小于监测点 和 ；监测点 的 频段有大于 的小波强度响应（见图 箭头所指），且该频段的小波响应随流道逐渐减弱；在 工况下，叶轮进口监测点 和 的小波强度响应相较于叶轮出口将测定 和 更加的明显，其中监测点 高频脉动的小波响应最强，越接近叶轮出口越弱；叶轮流道的低频段压力脉动小波强度响应沿流道呈现先减小后增大再减小的变化趋势，此时的低小波强度响应主要频段为 ，该频段明显小 第 期刘安林，等基于小波变换的核主泵惰转过程叶轮压力脉动及泵内部流动规律于 工况时的初始叶轮通频，处于 倍转频和 倍转频之间；在 工况下，叶轮内压力脉动强度基本趋于稳定值，但在叶轮进口监测点 的高频响应区域明显大于其他监测点，此时强小波强度响应频段为 ，即大于 区域 核主泵内流场分析图 为在惰转过程中叶轮与导叶中心截面 工况下的静压和速度云图，可以看出：随着惰转发生，不同工况的压力与速度分布变化情况一致，且变化梯度一致；随着流量减小，叶轮出口速度逐渐减小，并且在导叶内进一步减小；由于在惰转过程中叶轮转速会逐渐下降，叶轮流道内部静压也呈下降趋势；随着叶轮到导叶流道的逐渐扩张，流道内的速度减小不断转化为压力，使得静压逐渐增大，在叶轮出口边与导叶进口处，对应的导叶流道扩压能力会后滞；随着惰转进行，导叶进口速度减小，流道的扩压能力逐渐减弱，扩压能力减小，导叶流道内的做功区域逐渐向出口移动，这可能是由于随着流量较小，该区域的回流减弱，而导叶做功区域转移到了出口附近图 叶轮与导叶中心截面静压与速度云图 图 为在惰转过程中环形压水室中心截面 工况下的静压和速度云图图 环形压水室截面静压和出口速度云图 由图 可以看出：环形压水室的舌部区域流动十分复杂，在 工况，泵内流量还较大，叶轮转速开始降低，舌部有冲击现象导致此处存在静高压区；工况，流量减少 以上，环形压水室中的静压分布主要集中向壁面与出口扩散，因为转速与流量的迅速下降，造成内部流动紊乱；在 工况，流量减少 以上，压水舌部的高压流体有向局部低压区流动的趋势，使得压水室舌部的冲击有 排灌机械工程学报第 卷所减少；压水室隔舌部位伴随着明显的回流现象，随着惰转进行，回流现象越发严重，扩散面积增大，有往流道中心扩散的迹象，且旋涡的个数有增多趋势，内部存在较大流动损失，流量逐渐减小 核主泵涡结构分析利用 值的大小揭示涡核的分布及强度，图 为核主泵在惰转过程中叶轮、导叶、环形压水室中间截面的 值分布图 叶轮、导叶及环形压水室中间截面 值分布 ，由图 可以看出：一个周期内不同时刻叶轮内的涡核分布逐渐较少且较为紊乱；在 工况，叶轮流道内部有较多强涡区域，主要分布在叶片吸力面中部以及叶片尾缘部分，在吸力面强涡分布主要呈现出长条形涡结构；在 工况，强涡区域面积明显减小；在 工况，叶轮流道内强涡区域基本消失，值在流道内分布均匀且趋于 ，叶轮吸力面的长条形涡结构明显，并且此长条涡结构一直存在且较稳定，由于吸力面与压力面压力相差较大，使得流体向吸力面运动由图 可以看出：值范围与叶轮一致，各工况强涡区域分布紊乱；在 工况，强涡主要分布在流道中部且涡核靠近吸力面测呈长条状，出口区域存在部分中等涡；在 工况，强涡区域面积减小，中等涡面积逐渐增大，由于叶轮转速与流量减小，流道内部分强涡脱落，并有向导叶出口移动的趋势；在 工况，流道内强涡基本消失；随着惰转进行，不同工况的导叶强涡面积与数量不断减小，流道内值趋于 由于叶轮转速与流量快速下降使得叶轮和导叶干涉剧烈，导致导叶内涡强度有所增大由图 可以看出：当 值为 ，时能够清晰观察到部分涡核，压水室的 值范围相对叶轮及导叶的 值范围大幅度增大，这是由于随惰转过程进行，叶轮速度与流量迅速减小，而经过叶轮做工，叶轮内的液流以较高速度进入导叶流道，造成了较大的冲击损失，对液流有一定的减速作用，使得液流流出时速度相对较低，所以压水室涡核 值范围增大，这也说明了叶轮与导叶间的动静干涉作用影响较大，而环形压水室与静止导叶相连接，受动静干涉作用影响较小；环形压出室出口涡 第 期刘安林，等基于小波变换的核主泵惰转过程叶轮压力脉动及泵内部流动规律核随叶轮旋转逐渐分离，且出口处有部分涡核，可能与环形压水室结构有关；叶轮旋转方向为逆时针，导致压水室左侧的涡核分布比右侧多，且强涡主要集中在左侧 结论应用数值计算的方法，基于小波变换研究了核主泵惰转下的压力脉动特性，得到结论如下：）在惰转 到 过程中，核主泵叶轮流道最大压力脉动系数变化出现在叶轮中部 到 过程中，叶轮流道最大压力脉动系数变化出现在动静交界区域，其变化率大于 ；在 到 过程中，压力脉动系数变化趋于平稳，高频小波区域占比逐渐减小）惰转过程中，随着叶轮转速降低和内部流量减小，叶轮做功区域向出口附近转移 由于叶轮与导叶之间的动静干涉作用，液流在流入压水室时速度相对减小，因此压水室内冲击较小 隔舌部位由于结构特殊所受冲击较大，附近有回流产生，并有向流道中部扩散的迹象）随着惰转进行，流道内的部分强涡区数量逐渐减小，且出现涡核分离脱落现象，由于叶轮与导叶之间的动静干涉导致导叶流道内的 值分布不均匀，随着叶轮旋转涡核向导叶出口流动，压水室内的涡核主要集中于左侧参考文献（）张恒，刘雷，刘立军 非均匀入流对 核主泵内流及性能的影响研究 西安交通大学学报，（）：，（）：（）雷明凯，陈语曦，朱宝，等 核主泵压力脉动及其诱导振动研究进展 中国核电，（）：，（）：（）赖喜德，叶道星，陈小明，等 类球形压水室的混流式核主泵压力脉动特性 动力工程学报，（）：，（）：（），：朱荣生，龙云，付强，等 核主泵小流量工况压力脉动特性 振动与冲击，（）：，（）：（），：倪丹，杨敏官，高波，等 混流式核主泵内流动结构与压力脉动特性关联分析 工程热物理学报，（）：，（）：（）王鹏，袁寿其，王秀礼，等 大流量工况下核主泵内部不稳定特性分析 振动与冲击，（）：，（）：，（），：，：，（），（）：（下转第 页）第 期王维军，等微型高速燃油泵非定常流动及振动特性 ，（）：（）赵伟国，王桂鹏，赵国寿，等 低比转速离心泵非定常空化流动特性 兰州理工大学学报，（）：，（）：（），（）：邱昊，李志鹏，肖响 基于偏置分流叶片低比转速离心泵压力脉动特征研究 中国农村水利水电，（）：，（）：（），：王维军，李泰龙 基于分流叶片技术的离心式航空燃油泵优化设计 武汉大学学报（工学版），（）：，（）：（）唐庆，单金光，周龙，等 航空燃油齿轮泵流动机理的数值模拟研究 机电工程，（）：，（）：（）（责任编辑陈建华）檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨檨（上接第 页），（）：叶道星，赖喜德，文海罡，等 核主泵断电惰转停机过渡过程性能特性研究 动力工程学报，（）：，（）：（），：赵凡超，戴石良，房华伟，等 基于改进集成多隐层小波极限学习神经网络的滚动轴承故障识别研究 机电工程，（）：，（）：（）（责任编辑陈建华）