静叶倾斜对双级轴流风机性能和噪声影响的数值研究_李优.pdf

第 卷第 期 年 月热能动力工程 ，收稿日期：；修订日期：作者简介：李 优（），男，中国大唐集团科学技术研究总院有限公司华北电力试验研究院助理工程师文章编号：（）静叶倾斜对双级轴流风机性能和噪声影响的数值研究李 优，童家麟，夏尊宇，叶学民（中国大唐集团科学技术研究总院有限公司华北电力试验研究院，北京；国网浙江省电力有限公司电力科学研究院，浙江 杭州；内蒙古大唐国际托克托发电有限责任公司，内蒙古 呼和浩特；华北电力大学 动力工程系，河北 保定）摘 要：以某 机组配套的两级轴流风机为研究对象，采用 软件对其进行三维数值计算，研究了静叶倾斜对风机性能和噪声的影响。结果表明：当正倾斜角度 时，在设计工况风机效率提高 、全压降低.；静叶正倾斜角度大于 或反倾斜时，风机效率和全压均降低。静叶正倾斜使两级静叶区的总声压级均有所降低，正倾斜 时在设计工况点第 级静叶区总声压级降低 ，且倾斜角度越大降噪效果越显著；静叶反倾斜后，第 级静叶噪声基本不变，第 级静叶区总声压级增大 。综合分析静叶倾斜对风机性能和噪声的影响，静叶正倾斜 为最优方案。关 键 词：静叶倾斜；双级轴流风机；气动性能；噪声特性；数值模拟中图分类号：文献标识码：引用本文格式李 优，童家麟，夏尊宇，等 静叶倾斜对双级轴流风机性能和噪声影响的数值研究 热能动力工程，（）：，（）：，（，：；，：；，：；，：）：，；，第 期李 优，等：静叶倾斜对双级轴流风机性能和噪声影响的数值研究 ，：，引 言风机作为流体输运的重要动力设备，在电力、农业、冶 金 和 航 空 航 天 等 领 域 广 泛 应 用。目 前，以上的火电机组已成为主力机组。轴流风机具有调节效率高、对风道系统适应性好等特点，从而广泛应用于大型火电机组中。在轴流机械中，静叶结构形式对其内部流动特征有较大影响，进而影响其气动性能及噪声特性。根据旋转机械单音噪声的产生机理可知，单音噪声的频率是叶片通过频率的倍数，其噪声强度与动静间干涉作用密切相关。程颢颐等人采用非定常雷诺时均模型和管道声类比理论模拟了倾斜 后掠静叶对风扇单音噪声的影响，研究表明：正倾斜静叶可降低噪声的声功率级，但其气动性能会降低；当倾斜角度增大时，降噪效果显著；正倾斜 时，各谐波阶次降噪量均超过 。等人利用航空发动机的低速风扇探究了静叶倾斜和掠型的降噪效果，结果显示，静叶倾斜与后掠组合的方案可有效降低风扇的单音噪声，且降噪效果优于增加动叶与静叶间距离。和 采用 提出的模型对风扇静子叶片的倾斜和掠型进行了研究，研究表明，转子尾迹与静叶间的关系是影响降低单音噪声的主要因素，且叶片数越多降噪效果越明显。等人研究了静叶后掠对风扇离散噪声和宽频噪声的影响，研究表明，转子尾迹数大于径向模态个数时静叶后掠才可以降低风扇的离散噪声，且降低宽频噪声需要静叶后掠角度足够大。等人研究了风扇出口静叶的掠型和倾斜对噪声的影响，研究表明，降低动静干涉的宽频噪声最有效的方法是减少叶片数量，静叶正倾斜和后掠是降低单音噪声最有效的方法。丁骏等人 研究了亚音速扩压叶栅端弯的积迭方式，研究表明，直叶片通道采用前缘端弯技术效果较好，而正弯叶片通道采用尾缘端弯的方式较好。等人对比了压气机采用常规导叶与倾斜导叶时的性能，研究表明，正倾斜导叶的压比和效率高于常规导叶。刘海建等人探究了静叶倾角对动静干涉噪声的影响，研究表明，总声压级中单音噪声是主要因素，静叶倾斜可降低干涉噪声，且静叶后倾（与转子转动方向相同）效果优于前倾。综上所述，压气机中的静叶多为翼型叶片，采用静叶倾斜和弯曲可降低气动噪声。对于双级轴流风机，其两级静叶结构明显不同。第 级静叶多为长短复合叶片，其内部流动复杂，第 级为等长度圆弧板叶片。因此，探究两级静叶同时倾斜对风机性能及噪声特性的影响，对于该类风机的设计优化具有重要意义。数值计算方法 几何模型轴流风机三维几何模型如图 所示。计算域包括集流器区、两级动叶区、两级静叶区和扩压器区。动叶轮直径为 ，轮毂直径 ，叶顶间隙 。两级动叶均采用截面翼型结构，每级动叶片数为。第 级静叶为长短复合圆弧板，第 级静叶为长度相同的圆弧板，每级静叶叶片数为。风机叶轮逆时针旋转（从电机侧看），设计工况点参数：体积流量 ，全压 。图 轴流风机几何模型 数值计算模型采用有限体积法进行数值计算，通过 热能动力工程 年 进行网格划分，其中集流器和扩压器区生成六面体网格来加速收敛，在动叶区和静叶区生成贴合性较好的四面体网格，并采用尺寸函数对叶片前缘、尾缘和叶顶处进行局部加密。为保证壁面边界层区域计算的准确性，需第 层网格节点位于边界层的对数率区，即壁面的无量纲距离 在 之间，经计算，第 层网格高度为 ，网格划分结果如图 所示。图 网格划分 采用 软件求解三维定常不可压缩雷诺时均方程和连续性方程，选择 湍流模型对方程组进行封闭。壁面为无滑移边界，忽略壁面粗糙度，近壁处采用标准壁面函数法求解，压力 速度耦合选择 算法。入口边界采用速度入口，出口为自由出流。采用多重参考系模型解决动静叶片之间的干涉，动静叶片交界面通过非连续 实现数据传递。当进出口流量偏差小于，、和各方向速度最大残差低于，且进、出口截面总压均不随迭代步数改变时，计算达到收敛。噪声模拟采用 混合算法，以稳态计算结果作为初场，选择 亚网格尺度模型的大涡模拟进行非定常计算，动静交界面和动叶区采用滑移网格计算，压力速度耦合采用 算法。通过设定监测点获取压力脉动时域特征，经快速傅里叶变换（）将压力时域特征转化为声压频谱特征。为排除初场对计算结果的影响，选择第 个旋转周期的计算数据进行分析。数值计算时，时间步长设为叶轮旋转 的时间 ，为确保风机叶轮旋转 周，计算步数设为 步。总声压级 如式（）所示：（）式中：功率谱密度，；参考声压，为。模型验证为验证网格尺寸对计算结果的影响，分别采用 万、万、万、万和 万网格数对风机进行模拟计算，结果如表 所示。在设计工况下当网格数达到 万，随网格数量增加，所得风机全压、效率基本不变。综合考虑计算精度、计算时长对模拟结果的影响，选择网格数量为 万的模型进行计算。表 网格无关性验证 网格数 万效率偏差 全压偏差 由表 可知，在设计流量下，模拟值与样本值的全压偏差和效率偏差分别为和。为保证数值计算结果的准确性，在 的流量范围内将模拟结果与样本值对比，如图 所示。由图 可知，数值模拟结果与风机样本性能曲线的趋势一致，全压和效率的偏差均小于。因此，所建模型和计算方法能准确反映风机的运行特性。图 模拟结果与样本值对比 静叶倾斜方案倾角 定义为静叶尾缘与轮毂法向的夹角，原结构中静叶倾角 ，如图 所示。为研究风机性能对静叶倾角的敏感性，对 ，第 期李 优，等：静叶倾斜对双级轴流风机性能和噪声影响的数值研究时的风机性能进行数值计算。其中，叶片倾斜方向与动叶旋转方向相同时定义为正倾斜（），相反则为反倾斜（）。图 倾斜静叶示意图 结果分析 性能参数图 为 静 叶 倾 斜 前、后 风 机 的 性 能 曲 线。图（）显示，静叶倾斜后风机全压降低，且倾斜角度越大降低越突出。在设计工况点，时全压下降约 ，时全压下降约 。图（）表明，倾角相同时正倾斜风机的效率高于反倾 斜。在 静 叶 正 倾 斜 的 情 形 下，时 在 范围内风机效率高于原风机，在 和 范围内效率则均略低于原风机；当 时效率显著降低。在设计工况点，时风机的效率最高，相比原风机提高了 ；在静叶反倾斜情形下，风机效率均低于原风机。静叶效率可反映出静叶区流体损失比例。图（）表明，第 级静叶效率随流量增大而降低。静叶正倾斜后，第 级静叶效率在小流量侧有小幅度提升，且静叶倾角较大时效率提升幅度也增大。在设计工况点，时风机静叶效率最高，相比原结构提高。，时，风机静叶效率提升 。在大流量侧倾斜静叶的效率低于原风机，且 时静叶效率大幅度降低，在 时降低 。静叶反倾斜时，在全流量范围内，第 级静叶效率均降低，且静叶效率降低值随倾斜角度增大而越加明显。在设计工况点，静叶效率降低，在 时降低 。图 静叶倾斜对风机性能的影响 图 不同静叶倾角下的静叶效率 热能动力工程 年图（）显示，第 级静叶效率随流量增大先增大后减小，在 时静叶效率最大。叶片倾斜后，静叶效率变化较大。这是因为静叶倾斜后，进入第 级的流场变得更加复杂，且第 级静叶为长度相同的圆弧板叶片，对第 级静叶区的流动状态影响较大。在设计工况点，仅 时静叶效率大于原结构，相比原结构提升 。倾角较小时，静叶效率变化较小，且均在 时静叶效率达到最大值。倾角达到 时，正倾斜角和反倾斜角的静叶效率均大幅度降低，且反倾斜静叶效率降低幅度更大，在设计工况点降低。表 为静叶倾斜前、后风机各级参数变化。由表 可知，静叶反倾斜后两级叶轮的全压和效率均降低，且倾斜角度越大，性能降低越明显，第 级叶轮性能参数降低幅度要大于第 级。这是因为第 级静叶倾斜后，对第 级叶轮出口处的流动状态产生了较大影响。静叶正倾斜后，当 时，第 级叶轮的性能有所提升，但第 级的性能下降。由此可知，静叶小角度正倾斜时性能降低是由第 级叶轮性能降低造成的；当 后，两级叶轮的全压和效率均降低，这是因为当叶片倾斜角度过大时增大了局部阻力损失，降低了静叶区的流通能力。表 不同静叶倾角下的级性能参数变化 方案第 级第 级相对全压 相对效率 相对全压 相对效率 原结构 熵产率分布通过熵产率分布可较直观地看出风机产生损失的部位，图 为静叶倾斜前后第 级静叶 叶高处的熵产率分布。图中显示，熵产率较高区域主要分布在静叶前缘吸力面及静叶的尾缘，这是因为来自叶轮的流体进入静叶区后，撞击到叶片前缘，其在吸力面处由于附面层分离引起较强的速度脉动，从而产生了较强的湍流耗散。由于该风机的静叶结构为等厚圆弧板，在叶片的尾缘会有流体脱离形成较多脱落涡，从而增加了湍流耗散。通过对比可知，静叶倾斜后静叶附近的熵产率增大。尤其 -时，在叶片吸力面形成大范围熵产率较高区域，且叶片的前缘和尾缘熵产率值较大。这主要是因为，在该区域内有较大范围低速区，存在较多的低能流体，从而产生较高的流动损失。图 第 级静叶 叶高处的熵产率分布 第 期李 优，等：静叶倾斜对双级轴流风机性能和噪声影响的数值研究 湍动能湍动能表示单位质量流体因湍流脉动而具有的能量，其值可表征湍流脉动的强度，同时也可以从侧面反映气动噪声的分布特征。图 为静叶倾斜前后 叶高处的湍动能分布。由图 可知，原风机第 级静叶通道中湍动能分布相对均匀，湍动能主要沿流动方向分布，在叶片的前缘吸力面处和叶片尾缘有小范围湍动能较高区域。静叶正倾斜时，湍动能较大区域分布在静叶两侧及叶片尾缘，且湍动能最大值位于叶片吸力面靠近前缘处，在流道的中间部分湍动能值较小。静叶反倾斜时，湍动能较大区域分布在静叶吸力面附近，且湍动能最大值约是原结构的 倍，而湍动能较小的流体主要分布在叶片压力面附近。这主要是因为静叶反倾斜后，叶片吸力面与机壳的夹角减小，由于气体粘性作用，使吸力面叶顶处聚集较多低能气流，从而引起了较强的湍流脉动。图 第 级静叶 叶高处的湍动能分布 噪声特性为分析静叶倾斜对风机噪声产生的影响，将压力脉动监测点设在静叶区出口通道中间位置 叶高处，且随静叶倾斜角度的变化，监测点位置也发生变化，如图 所示。图 监测点位置 时域特征图 给出了设计流量下，监测点处压力脉动的时域特征。由图可知，流道内的压力呈周期性波动，且在 个旋转周期内，压力波动的峰值有 个，与叶轮动叶个数相同。这是因为，风机内气体所获能量来自旋转叶片做功，当叶片周期性通过监测点时，使气体静压产生周期性波动。静叶倾斜角度不同时，压力的峰值所对应的时间不同。这是因为静叶倾斜后，静叶区流道的中心位置发生变化，而叶轮初始位置不变，其经过监测点所用时间不同，因此监测点处压力峰值所对应时间略有差异。图（）显示，原结构风机压力脉动的幅值约热能动力工程 年为 ，静叶正倾斜后，压力脉动的幅值明显减小，且倾斜角度越大压力波动的幅值越小。当 时，压力波动的幅值约为 。从图（）可知，静叶反倾斜后，监测点处压力波动的幅值变化较小。相比于原结构，-时压力波动的幅值略有减小。另外，静叶正