CN111832210B - 一种用于单级离心鼓风机气动噪声的仿真预测方法 - Google Patents
一种用于单级离心鼓风机气动噪声的仿真预测方法 Download PDFInfo
- Publication number
- CN111832210B CN111832210B CN202010724691.5A CN202010724691A CN111832210B CN 111832210 B CN111832210 B CN 111832210B CN 202010724691 A CN202010724691 A CN 202010724691A CN 111832210 B CN111832210 B CN 111832210B
- Authority
- CN
- China
- Prior art keywords
- noise
- centrifugal blower
- stage centrifugal
- calculation
- simulation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D27/00—Control, e.g. regulation, of pumps, pumping installations or pumping systems specially adapted for elastic fluids
- F04D27/001—Testing thereof; Determination or simulation of flow characteristics; Stall or surge detection, e.g. condition monitoring
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mechanical Engineering (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明属于仿真预测技术领域,具体公开了一种用于单级离心鼓风机气动噪声的仿真预测方法,包括以下步骤:S1:建立满足流场计算的模型,并进行模型的网格划分;S2:进行定常流场仿真计算,定常流场仿真计算收敛后得到有关的气动性能指标;S3:进行单级离心鼓风机非定常流场仿真计算,并输出与噪声激励有关的物理量;S4:建立气动噪声计算模型;S5:将S3中与噪声激励有关的物理量做声源转化处理;将转化的等效声源加载到S4中的噪声计算模型中进行噪声计算,得到噪声分析结果;S6:根据噪声分析结果预测噪声大小以及判断噪声品质,从而制定优化降噪方案。上述方法,能够解决现有的气动噪声仿真计算方法无法适用于高速单级离心鼓风机噪声仿真计算的问题。
Description
技术领域
本发明属于仿真预测技术领域,尤其涉及一种用于单级离心鼓风机气动噪声的仿真预测方法。
背景技术
单级离心鼓风机因其高效、节能等优点在污水处理曝气、电厂脱硫脱硝以及氧化工艺等领域得到广泛应用。单级离心鼓风机运行时,其内部高速气流与旋转件(叶轮)、静止件(进出口导叶以及蜗壳等)和进排气管等相互作用,产生具有鼓风机特有的叶片通过频率及其谐波频率的气动噪声,而且此噪音声级高,高频突出,传播距离远,污染范围大,特别是某些噪音的频率与人的内脏器官固有频率相接近,容易引起共振,使人产生头晕、恶心、心律过速、高血压等症状,不仅导致人们的工作和生活质量下降,而且容易引发安全事故和人际关系矛盾。
目前针对气动噪声的研究主要有理论预测、实验测试与数值模拟三种方法。理论预测虽然可以快速的对气动噪声进行评估,但对于鼓风机由于预测存在大量的假设,并且无法对声源给出细节刻画,无法针对真实结构给出定量的预测,从而受到应用限制;实验测量主要包括流场测量、声源定位和声场测量,虽然可以精确得到产品的噪声等级,但由于鼓风机产品结构复杂程度以及尺寸的增大使得噪声测试环境变得复杂,需要更为昂贵的资金与人力投入。数值方法可以精确刻画声源细节,但由于对网格尺度与数量要求高,使得多尺度问题计算量极大。
并且由于鼓风机机壳的复杂性,声波在不规则的蜗壳内表面会经过多次反射和散射,因此实际的鼓风机内的声场与自由空间中的声场有很大不同,现有常规的气动噪声仿真方法是难以考虑蜗壳的对声场的影响。
因此,对于高压比、高转速、复杂机壳的离心鼓风机,目前还没有形成较为完善的有效的气动噪声仿真预测方法。
发明内容
本发明的目的在于提供一种用于单级离心鼓风机气动噪声的仿真预测方法,以解决现有的气动噪声仿真计算方法无法适用于高速单级离心鼓风机噪声仿真计算的问题。
为了达到上述目的,本发明的技术方案为:一种用于单级离心鼓风机气动噪声的仿真预测方法,包括以下步骤:
S1:根据单级离心鼓风机机组的装配方式,建立满足单级离心鼓风机机组流场计算的模型,并进行模型的网格划分;
S2:进行单级离心鼓风机定常流场仿真计算,定常流场仿真计算收敛后得到有关的气动性能指标;
S3:根据S2中的气动性能指标,进行单级离心鼓风机非定常流场仿真计算,并输出与噪声激励有关的物理量;
S4:建立单级离心鼓风机的气动噪声计算模型;
S5:将S3中与噪声激励有关的物理量做声源转化处理;将转化的等效声源加载到S4中的噪声计算模型中进行噪声计算,得到噪声分析结果;
S6:根据噪声分析结果预测单级离心鼓风机噪声大小以及判断噪声品质,从而制定优化降噪方案。
进一步,在步骤S1中,单级离心鼓风机的叶轮、进出口导叶和蜗壳根据其几何特征线在旋转机械专用软件中进行结构网格划分。
进一步,对所有流体域划分边界层并进行网格加密处理。
进一步,在步骤S2中,气动性能指标包括压比、效率和功率。
进一步,在步骤S3中,与噪声激励有关的物理量包括单级离心鼓风机的主要发声部件表面的非定常压力脉动。
进一步,在步骤S4中,根据单级离心鼓风机机组的装配方式对气动噪声计算模型进行有限元网格划分。
进一步,具体步骤是在所述流场计算的模型的基础上建立单级离心鼓风机的空气声场传播模型,然后对空气声场传播模型进行声学网格填充;蜗壳、扩压器和管道设置全反射边界条件,出口管道处设置open duct AML surface结构,利用自动完美匹配层技术进行无反射全吸声条件处理。
进一步,主要发声部件包括叶轮、扩压器和蜗壳。
进一步,在步骤S5中,对单级离心鼓风机的叶轮表面的非定常压力脉动进行旋转偶极子声源转化,处理为旋转机械特有的扇声源;对单级离心鼓风机的扩压器和蜗壳表面的非定常压力脉动做静止偶极子声源转化,并将扩压器与蜗壳表面输出的非定常压力脉动积分插值映射到声学网格上。
进一步,在步骤S3中,对所述非定常压力脉动进行傅里叶变换处理。
本技术方案的有益效果在于:①对于常规的单级离心鼓风机噪声研究方法,对噪声源的发声原理、噪声特性以及噪声场分布都缺乏具体而直观的分析与研究,基于测试的整改方案导致过多的时间与资金成本投入。本技术方案的研究方法可以通过深入探究单级离心鼓风机噪声源(叶轮、扩压器、蜗壳等)的发声机理,挖掘流场细节与噪声之间的内在关系,形成对噪声机理较为明确的认识,因而可以从设计阶段就考虑降噪优化方案。②由于单级离心鼓风机机壳的复杂性,声波在不规则的蜗壳内表面会经过多次反射和散射,因此实际的单级离心鼓风机内的声场与自由空间中的声场有很大不同,现有常规的气动噪声仿真方法是难以考虑蜗壳的对声场的影响。本技术方案的研究方法通过对蜗壳、扩压器、管道等壁面设置全反射边界条件,并且在出口管道处设置open duct AML surface结构,可以考虑声波在不规则的蜗壳内表面会多次的反射和散射作用,同时利用噪声有限元计算方法中的自动完美匹配层技术,模拟出口管道的真实声辐射状态。
附图说明
图1为本发明一种用于单级离心鼓风机气动噪声的仿真预测方法的流程图。
具体实施方式
下面通过具体实施方式进一步详细说明:
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例基本如附图1所示:一种用于单级离心鼓风机气动噪声的仿真预测方法,包括以下步骤:
S1:根据单级离心鼓风机机组的装配方式,建立满足单级离心鼓风机机组流场计算的模型,并进行模型的网格划分;
S2:进行单级离心鼓风机定常流场仿真计算,定常流场仿真计算收敛后得到有关的气动性能指标,气动性能指标包括压比、效率和功率;
S3:根据S2中的气动性能指标,进行单级离心鼓风机非定常流场仿真计算,并输出与噪声激励有关的物理量,与噪声激励有关的物理量包括单级离心鼓风机的主要发声部件表面的非定常压力脉动,并对非定常压力脉动进行傅里叶变换处理,主要发声部件包括叶轮、扩压器和蜗壳;
S4:建立单级离心鼓风机的气动噪声计算模型,根据单级离心鼓风机机组的装配方式对气动噪声计算模型进行有限元网格划分;具体步骤是在所述流场计算的模型的基础上建立单级离心鼓风机的空气声场传播模型,然后对空气声场传播模型进行声学网格填充;蜗壳、扩压器和管道设置全反射边界条件,出口管道处设置open duct AML surface结构,利用自动完美匹配层技术进行无反射全吸声条件处理;
S5:将S3中与噪声激励有关的物理量做声源转化处理;将转化的等效声源加载到S4中的噪声计算模型中进行噪声计算,得到噪声分析结果;
S6:根据噪声分析结果预测单级离心鼓风机噪声大小以及判断噪声品质,从而制定优化降噪方案。
在步骤S1中,单级离心鼓风机的叶轮、进出口导叶和蜗壳根据其几何特征线在旋转机械专用软件中进行结构网格划分,对所有流体域划分边界层并进行网格加密处理。
在步骤S5中,对单级离心鼓风机的叶轮表面的非定常压力脉动进行旋转偶极子声源转化,处理为旋转机械特有的扇声源;对单级离心鼓风机的扩压器和蜗壳表面的非定常压力脉动做静止偶极子声源转化,并将扩压器与蜗壳表面输出的非定常压力脉动积分插值映射到声学网格上。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。
以上所述的仅是本发明的实施例,方案中公知的具体结构及特性等常识在此未作过多描述,所属领域普通技术人员知晓申请日或者优先权日之前发明所属技术领域所有的普通技术知识,能够获知该领域中所有的现有技术,并且具有应用该日期之前常规实验手段的能力,所属领域普通技术人员可以在本申请给出的启示下,结合自身能力完善并实施本方案,一些典型的公知结构或者公知方法不应当成为所属领域普通技术人员实施本申请的障碍。应当指出,对于本领域的技术人员来说,在不脱离本发明结构的前提下,还可以作出若干变形和改进,这些也应该视为本发明的保护范围,这些都不会影响本发明实施的效果和专利的实用性。本申请要求的保护范围应当以其权利要求的内容为准,说明书中的具体实施方式等记载可以用于解释权利要求的内容。
Claims (4)
1.一种用于单级离心鼓风机气动噪声的仿真预测方法,其特征在于:包括以下步骤:
S1:根据单级离心鼓风机机组的装配方式,建立满足单级离心鼓风机机组流场计算的模型,并进行模型的网格划分;
S2:进行单级离心鼓风机定常流场仿真计算,定常流场仿真计算收敛后得到有关的气动性能指标;气动性能指标包括压比、效率和功率;
S3:根据S2中的气动性能指标,进行单级离心鼓风机非定常流场仿真计算,并输出与噪声激励有关的物理量;与噪声激励有关的物理量包括单级离心鼓风机的主要发声部件表面的非定常压力脉动;主要发声部件包括叶轮、扩压器和蜗壳;
S4:建立单级离心鼓风机的气动噪声计算模型;根据单级离心鼓风机机组的装配方式对气动噪声计算模型进行有限元网格划分;具体步骤是在所述流场计算的模型的基础上建立单级离心鼓风机的空气声场传播模型,然后对空气声场传播模型进行声学网格填充;蜗壳、扩压器和管道设置全反射边界条件,出口管道处设置open duct AML surface结构,利用自动完美匹配层技术进行无反射全吸声条件处理;
S5:将S3中与噪声激励有关的物理量做声源转化处理;将转化的等效声源加载到S4中的噪声计算模型中进行噪声计算,得到噪声分析结果;具体步骤是对单级离心鼓风机的叶轮表面的非定常压力脉动进行旋转偶极子声源转化,处理为旋转机械特有的扇声源;对单级离心鼓风机的扩压器和蜗壳表面的非定常压力脉动做静止偶极子声源转化,并将扩压器与蜗壳表面输出的非定常压力脉动积分插值映射到声学网格上;
S6:根据噪声分析结果预测单级离心鼓风机噪声大小以及判断噪声品质,从而制定优化降噪方案。
2.根据权利要求1所述的一种用于单级离心鼓风机气动噪声的仿真预测方法,其特征在于:在步骤S1中,单级离心鼓风机的叶轮、进出口导叶和蜗壳根据其几何特征线在旋转机械专用软件中进行结构网格划分。
3.根据权利要求2所述的一种用于单级离心鼓风机气动噪声的仿真预测方法,其特征在于:对所有流体域划分边界层并进行网格加密处理。
4.根据权利要求1所述的一种用于单级离心鼓风机气动噪声的仿真预测方法,其特征在于:在步骤S3中,对所述非定常压力脉动进行傅里叶变换处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010724691.5A CN111832210B (zh) | 2020-07-24 | 2020-07-24 | 一种用于单级离心鼓风机气动噪声的仿真预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010724691.5A CN111832210B (zh) | 2020-07-24 | 2020-07-24 | 一种用于单级离心鼓风机气动噪声的仿真预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111832210A CN111832210A (zh) | 2020-10-27 |
CN111832210B true CN111832210B (zh) | 2022-04-08 |
Family
ID=72925919
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010724691.5A Active CN111832210B (zh) | 2020-07-24 | 2020-07-24 | 一种用于单级离心鼓风机气动噪声的仿真预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111832210B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113378437B (zh) * | 2021-06-16 | 2022-11-04 | 山东大学 | 一种流固耦合作用下轴流风机噪声仿真预测方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077292A (zh) * | 2013-01-16 | 2013-05-01 | 江苏大学 | 一种离心泵水力噪声预测方法 |
CN106503323A (zh) * | 2016-10-17 | 2017-03-15 | 江苏大学 | 一种多级离心泵流动诱导结构辐射噪声数值预测方法 |
CN109977533A (zh) * | 2019-03-22 | 2019-07-05 | 中车永济电机有限公司 | 牵引电机风扇噪声的仿真预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013134742A (ja) * | 2011-12-27 | 2013-07-08 | Daihatsu Motor Co Ltd | 空力騒音レベルの予測方法 |
-
2020
- 2020-07-24 CN CN202010724691.5A patent/CN111832210B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103077292A (zh) * | 2013-01-16 | 2013-05-01 | 江苏大学 | 一种离心泵水力噪声预测方法 |
CN106503323A (zh) * | 2016-10-17 | 2017-03-15 | 江苏大学 | 一种多级离心泵流动诱导结构辐射噪声数值预测方法 |
CN109977533A (zh) * | 2019-03-22 | 2019-07-05 | 中车永济电机有限公司 | 牵引电机风扇噪声的仿真预测方法 |
Non-Patent Citations (1)
Title |
---|
一种新型离心风机的气动噪声数值仿真;郑欣等;《自动化与仪器仪表》;20140930(第9期);正文第1-4节 * |
Also Published As
Publication number | Publication date |
---|---|
CN111832210A (zh) | 2020-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gabard et al. | Theoretical model for sound radiation from annular jet pipes: far-and near-field solutions | |
US8177496B2 (en) | Tone noise reduction in turbomachines | |
Luo et al. | Tip leakage flow and aeroacoustics analysis of a low-speed axial fan | |
Liu et al. | Numerical and experimental investigations of centrifugal compressor BPF noise | |
CN103631989A (zh) | 一种离心泵流动诱导噪声数值预测方法 | |
CN103077292A (zh) | 一种离心泵水力噪声预测方法 | |
CN111832210B (zh) | 一种用于单级离心鼓风机气动噪声的仿真预测方法 | |
Yang et al. | Numerical and experimental study on flow-induced noise at blade-passing frequency in centrifugal pumps | |
Aksenov et al. | Development of the virtual experimental bench on the basis of modernized research centrifugal compressor stage test unit with the 3D impeller. | |
CN113378437B (zh) | 一种流固耦合作用下轴流风机噪声仿真预测方法及装置 | |
Guan et al. | Experimental mode decomposition investigation on 3-stage axial flow compressor stall phenomena using aeroacoustics measurements | |
Dong et al. | Theoretical characterization and modal directivity investigation of the interaction noise for a small contra-rotating fan | |
Pietroniro et al. | Effects of boundary layer and local volumetric cells refinements on compressor direct noise computation | |
Nark et al. | Broadband liner optimization for the source diagnostic test fan | |
CN113503262B (zh) | 一种两级增压系统高低压级涡轮的仿真方法 | |
Luan et al. | Numerical study on aerodynamic noise performances of axial spacing in a contra-rotating axial fan | |
Tan et al. | Coupling bionic design and numerical simulation of the wavy leading-edge and seagull airfoil of axial flow blade for air-conditioner | |
Liu et al. | Analysis of intake silencer insertion loss in a marine diesel engine turbocharger based on computational fluid dynamics and acoustic finite element method | |
Yipeng et al. | Sensitivity analysis of impeller blade parameters to compressor performance and aerodynamic noise | |
Ketata et al. | A DFT Spectrum Acoustic Analysis for Investigating Pulse Duration Effect on Performance, Psychoacoustic Sound Level of Turbocharger Turbines Through C++ FDM Code | |
Lu et al. | Numerical optimization for radiated noises of centrifugal pumps in the near-field and far-field based on a novel MLGA-PSO algorithm | |
Zhou et al. | Multi-Objective optimization of an IGV for a large axial fan based on NSGA-II | |
Xin et al. | Numerical Investigation of the Flow Field and Aerodynamic Load on Impellers Under Radial Inlet Guide Vanes With Trailing Edge Blowing in Centrifugal Compressor | |
Zhao et al. | A RAPID PREDICTION METHOD OF BROADBAND NOISE OF SINGLE-STAGE FANS | |
Song et al. | Effect of balanced edge cutting configurations on the aeroacoustics of high-speed turbomachinery |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |