CN115541006A - 基于平面piv的预测指定声源远场流致噪声的方法 - Google Patents

基于平面piv的预测指定声源远场流致噪声的方法 Download PDF

Info

Publication number
CN115541006A
CN115541006A CN202211507962.7A CN202211507962A CN115541006A CN 115541006 A CN115541006 A CN 115541006A CN 202211507962 A CN202211507962 A CN 202211507962A CN 115541006 A CN115541006 A CN 115541006A
Authority
CN
China
Prior art keywords
field
flow
information
piv
pressure
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.)
Granted
Application number
CN202211507962.7A
Other languages
English (en)
Other versions
CN115541006B (zh
Inventor
叶青青
陈朗生
陈奕宏
邵雪明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202211507962.7A priority Critical patent/CN115541006B/zh
Publication of CN115541006A publication Critical patent/CN115541006A/zh
Application granted granted Critical
Publication of CN115541006B publication Critical patent/CN115541006B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • 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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Computational Mathematics (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于平面PIV的预测指定声源远场流致噪声的方法，属于粒子图像测速领域。本发明由二维PIV获取的速度场信息出发,通过欧拉法构建并求解泊松方程获取流场压力；基于所求解的高响应压力信息，通过Curle声比拟方法，求得远场声压；再进行展向二维PIV测量，获得展向速度信息，经过展向相关性修正，获得远场噪声信息。本发明通过使用简单的二维PIV测量即可得到指定声源考虑三维效应的流致噪声信息。

Description

基于平面PIV的预测指定声源远场流致噪声的方法
技术领域
本发明涉及流体测量领域中的粒子图像测速（PIV, Particle ImageVelocimetry）技术，具体涉及一种基于平面PIV的预测指定声源远场流致噪声的方法。该方法基于平面时间解析PIV测量，通过流场压力场重构，声比拟预测，以及展向相关性修正等技术，对指定流场声源的考虑三维效应的流致噪声进行预测。
背景技术
流致噪声是20世纪50年代从流体力学和声学两个基础学科产生出来的交叉学科，主要研究由流体介质流动对声音产生和传播的影响。21世纪以来，民用航空，风力发电，海洋作业等行业快速发展，与此同时带来了不计其数的流致噪声效应，给人民的生产生活带来巨大困扰，流致噪声的机理研究与抑制需求与日俱增。
传统的噪声测量诸如空气中的远场麦克风与水下的水听器能直接获取流致噪声，但易引入除关注声源之外的其他流致噪声或环境噪声，甚至设备的噪声如风洞水洞的电机噪声会淹没掉所关注的流致噪声；消声风洞，消声水洞能消除额外的环境噪声，但却存在设备建设成本高昂的限制；麦克风阵列可有效定位特定声源但成本极高。并且以上阐述的噪声测量方式均无法将流场速度信息与噪声直接对应。若需研究流场速度与噪声的直接对应关系又需要同步测量速度场，这样会额外引入新的噪声源，如PIV的冷却器噪声等。
如专利申请号CN112763180A公开了一种声学风洞内高速列车受电弓模型气动噪声快速预测方法，虽然能快速预测受电弓模型气动噪声，但并未将其与流场结合，并且需要配合实验数据库得到测理论模型的各项系数。需要大量噪声实验测量的同时也无法将流致噪声与流场结合分析。
又如专利申请号CN111608927A公开了一种离心泵流场-压力场-声场同步测量的装置和方法。提到了需要用消声水洞抑制节流阀的流致噪声，此外还需要在水听器附近安装抑振器减小管路噪声的水声信号测试的影响。这些用于消除额外环境噪声的改进仪器构建的使用无疑会增大实验的成本与难度。
PIV测量技术具有非侵入，全局测量的优点。时间解析PIV测量则在PIV测量技术上拓展了时间维度，能捕捉到流场的高频速度信息，可从中解析出流场的高频动力学特性。通过时间解析PIV所测得的流场速度信息直接预测流致噪声，不仅能避免额外声源与环境噪声的引入，而且能将流场速度信息与流致噪声耦合，为流致噪声的抑制技术提供理论指导。
发明内容
有鉴于此，本发明旨在提供一种基于平面时间解析PIV的预测指定声源远场流致噪声的方法。该方法通过压力场重构，声比拟噪声预测，展向相关性修正等一系列技术，得到考虑三维效应的远场流致噪声。将远场流致噪声与流场耦合，可直观地从流场角度分析流致噪声的产生机理，为后续流致噪声抑制技术提供指导。并在通过平面时间解析PIV测量获得流场高频速度信息的同时，附带给出流场的高频压力信息，有利于流场的动力学分析。
为实现上述目的，本发明采用以下技术方案：
本发明提供了一种基于平面PIV的预测指定声源远场流致噪声的方法，其包括如下步骤：
1）通过平面时间解析PIV测得流场中的高频速度信息；
2）基于流场高频速度信息,通过欧拉法求得流场压力梯度，构建并求解泊松方程，重构流场高频压力信息；
3）通过Curle声比拟方法求得远场声压，获得远场流致噪声信息；
4）通过展向平面时间解析PIV测量，获得流场流向速度的展向分布信息；
5）对流向速度进行展向相关性分析，基于流向速度展向相关性修正流致噪声，得到考虑三维效应的远场流致噪声信息。
作为本发明的优选方案，所述的步骤1）具体为：利用同步器将高频双腔激光器和CMOS高速相机同步以采集和储存流场粒子图像，采集到的粒子图像通过PIV分析处理得到流场的高频速度信息。
作为本发明的优选方案，所述的步骤2）具体为：通过欧拉法求得流场压力梯度
Figure 299868DEST_PATH_IMAGE001
，对压力梯度
Figure 77856DEST_PATH_IMAGE001
取散度得到拉普拉斯算子
Figure 720058DEST_PATH_IMAGE003
构建泊松方程；给定方程的边界条件， 通过迭代法求解各个时刻速度场对应的泊松方程，得到时间解析压力信息，重构流场高频 压力信息。
作为本发明的优选方案，步骤3）所述的通过Curle声比拟方法求得远场声压为：基于时间解析PIV所测得的高频速度信息与重构的高频压力信息计算Curle声比拟模型中的声源项，得到远场声压；计算各个时刻速度场对应的声压，求得时序声压信息，进一步求得流致噪声在频域下的声压级。
与现有技术相比，本发明所具有的有益效果包括：
1.本发明通过对指定声源流场进行PIV测量、并基于PIV测量速度信息进行压力场重构、声比拟预测与展向相关性修正等技术，克服了传统流致噪声测量引入其他噪声的缺陷，获得了指定流场不含环境噪声的远场噪声信息；
2.本发明通过进行流向与展向两次平面时间解析PIV、结合展向相关性修正技术，克服了传统多物理场采集需要考虑各设备干扰、成本高的缺陷，仅仅通过平面PIV测量便预测了考虑三维效应的远场流致噪声，并同时获得了流场的高频压力信息，同步获得了流场的速度、压力与噪声信息，为后续流场多物理量耦合分析提供了数据支持；
3. 本发明通过对指定声源流场进行PIV测量、压力场重构、声比拟预测与展向相关性修正等技术，直接将流场信息与远场流致噪声耦合，可为后续流致噪声抑制技术提供理论指导。
附图说明
图1为本发明实施例提供的基于平面时间解析PIV预测流致噪声的流程示意图。
图2为本发明实施例提供的平面PIV实验总体布置图。
图3为本发明实施例提供的PIV测试区域及指定声源示意图。
图4为本发明实施例提供的瞬时压力云图。
图5为本发明实施例提供的基于流向速度的展向相关系数分布图。
图6为本发明实施例提供的基于流向速度的展向相关长度图。
图7为本发明实施例提供的未修正的声压级，基于流向速度展向相关性修正后的声压级，水听器测量声压级对比图。
其中，2-1、高频双腔激光器，2-2、PIV光学镜组，2-3、流场测量视窗，2-4、反光镜，2-5、CMOS高速相机，2-6、数据采集分析系统，2-7、同步器。
具体实施方式
下面通过具体实施例对本发明的方法技术做进一步说明，需要指出的是，以下所述实施例旨在便于对本发明的理解，而对其不起任何限制作用：
如图1所示，本实施例提供了一种基于平面时间解析PIV预测流致噪声的方法，包括以下步骤：
(1)通过平面时间解析PIV获取流场速度信息。
根据实验待测流场工况搭建相应的实验平台，其总体布置如图2所示，包括高频双腔激光器2-1、 PIV光学镜组2-2、流场测量视窗2-3、反光镜2-4、 CMOS高速相机2-5、数据采集分析系统2-6、同步器2-7。其中高频双腔激光器2-1发射出高频激光通过PIV光学镜组2-2的光路调整后，形成片光源照亮流场测量视窗2-3；待数据采集分析系统2-6向同步器2-7发出采集指令后，同步器2-7向高频双腔激光器2-1与CMOS高速相机2-5同步发出触发信号，在高频双腔激光器2-1的激光照亮下，CMOS高速相机2-5以设定采集频率同步采集测量视窗的粒子图像；最后粒子图像通过数据采集分析系统2-6上的PIV分析软件处理，得到流场的高频速度信息。
注意到，由于流场中的固体模型（在本实施例中的固体模型为双圆柱，流场为水中的双圆柱绕流，此类固体模型也可以是翼型等）对激光有遮挡作用，背光区域无法照亮粒子，导致速度信息的缺失从而无法重构压力而影响到后续的声场预测；因此需要反光镜2-4对背光区域补光，从而通过PIV捕捉该区域的速度信息。
本实施例中研究的流场为水中的双圆柱绕流。由于双圆柱绕流主要声源为下游圆柱头部，优选地，PIV测试区域以下游圆柱为中心，视窗边缘到达上游圆柱圆心，如图3所示。
通过同步器2-7将高频双腔激光器2-1和CMOS高速相机2-5同步，以采集和储存粒子图像，采集到的粒子图像通过PIV分析软件处理得到流场的高频速度信息。
(2)基于流场高频速度信息，通过欧拉法求得压力梯度，求解泊松方程后获得流场的高频压力信息,完成压力场重构；
压力梯度由不可压N-S方程导出：
Figure 823144DEST_PATH_IMAGE004
。
式中，
Figure 841784DEST_PATH_IMAGE006
为压力梯度，
Figure 815556DEST_PATH_IMAGE007
为流体密度，
Figure 894240DEST_PATH_IMAGE008
为速度矢量，
Figure 219042DEST_PATH_IMAGE009
为流体动力粘度，
Figure 654090DEST_PATH_IMAGE010
为哈 密顿算子。
对压力梯度
Figure 262794DEST_PATH_IMAGE006
取散度得到拉普拉斯算子
Figure 997532DEST_PATH_IMAGE011
构建泊松方程，考虑到实施例流场 可看作二维不可压流场，满足无散条件，泊松方程如下式所示：
Figure 58898DEST_PATH_IMAGE012
；
式中p为流场压力，
Figure 435653DEST_PATH_IMAGE007
为流体密度，u为流向速度，v为法向速度；
给定方程的边界条件，本实施例泊松方程求解边界条件分为两种：
本实施例的第一类边界条件为：如图3所示，在平行于来流的上下测量区域边界处 湍流度小于1%，可看作无旋流动，按照下式给定压力值
Figure 774230DEST_PATH_IMAGE013
；
Figure 197646DEST_PATH_IMAGE014
其中
Figure 231461DEST_PATH_IMAGE013
为当前边界节点的压力，
Figure 130016DEST_PATH_IMAGE015
为当前边界节点的平均速度，
Figure 464045DEST_PATH_IMAGE016
为当前边界节 点的速度脉动，
Figure 789853DEST_PATH_IMAGE017
为来流压强由皮托管测出，
Figure 701177DEST_PATH_IMAGE018
为来流速度由水洞设备给出。
本实施例的第二类边界条件为：如图3所示，在垂直于来流的左右测量区域边界以 及圆柱表面给定已由欧拉法求得的压力梯度
Figure 154155DEST_PATH_IMAGE006
；
结合给出的边界条件，通过迭代法求解各个时刻速度场对应的泊松方程，得到时间解析压力信息。
这里，求解给定边界条件的泊松方程，可参照现有技术中已有方法，在此不再赘述。
流场的瞬时压力分布云图如图4所示，可以看出涡与低压区一一对应，压力场的求解结果符合理论。
(3)利用重构出的时间解析压力信息与PIV测量的时间解析速度信息通过Curle声比拟方法，计算得到远场流致噪声。
本实施例为不可压，不可穿透固壁面流场，远场声压的计算公式为：
Figure 606607DEST_PATH_IMAGE020
其中等号右边第一项和第二项分别为偶极子和四级子，在本实施例中分别代表圆柱表面压力脉动与流场中湍流结构引起的噪声。
其中，
Figure 119628DEST_PATH_IMAGE021
为时序声压信号，
Figure 377303DEST_PATH_IMAGE022
、
Figure 758606DEST_PATH_IMAGE023
表示以流场声源为坐标原点，观测者为终点位移 的分量；y为以流场声源中心为坐标原点，各个声源点的坐标；
Figure 536069DEST_PATH_IMAGE024
表示观察者与声源的距离；
Figure 469259DEST_PATH_IMAGE025
为声速；t为时间，
Figure 964962DEST_PATH_IMAGE026
为声信号的产生时间，
Figure 11940DEST_PATH_IMAGE027
为当前声信号的接收时间；
Figure 909489DEST_PATH_IMAGE028
，
Figure 888946DEST_PATH_IMAGE029
、
Figure 386792DEST_PATH_IMAGE031
分别为来流压强、密度，P、
Figure 985264DEST_PATH_IMAGE032
为当地流体的压 强，密度，
Figure 721008DEST_PATH_IMAGE033
为克罗内克符号；
Figure 12312DEST_PATH_IMAGE034
、
Figure 610171DEST_PATH_IMAGE035
为速度分量；
Figure 730443DEST_PATH_IMAGE036
为固壁面法向量；S为面积,
Figure 602584DEST_PATH_IMAGE037
为固 壁面表面积；V为体积，
Figure 314057DEST_PATH_IMAGE038
表示全流场体积;
Figure 6069DEST_PATH_IMAGE040
为lighthill张量。
计算各个时刻速度场对应的声压，求得时序声压信息。
通过下式求得流致噪声在频域下的声压级：
Figure 195611DEST_PATH_IMAGE041
；
其中SPL为声压级，
Figure 49822DEST_PATH_IMAGE042
为时序声压信号
Figure 417349DEST_PATH_IMAGE043
的自功率谱密度，
Figure 845925DEST_PATH_IMAGE045
为分析频率，
Figure 589890DEST_PATH_IMAGE046
为参考压力，本实施例流体介质为水，取
Figure 420312DEST_PATH_IMAGE047
。
（4）进行展向平面时间解析PIV测量获取流场展向区域的高频速度分布信息，展向PIV测量设备与图2中相同，仅仅是调整光路与CMOS相机，将测量视窗改为上下游圆柱之间区域的展向平面。
（5）获得流场的展向速度信息后，应选取主导声源附近位置的速度展向分布信息进行展向相关性分析，基于流向速度展向相关性修正流致噪声，得到考虑三维效应的远场流致噪声信息。
优选地，由于本实施例下游圆柱头部对噪声贡献大，且流向速度能反映流场结构对下游圆柱头部的撞击作用，取下游圆柱圆心前0.75D位置处的流向速度展向分布来计算展向相关长度。
确定位置后，按时序提取出该位置的流向速度展向速度分布时序信息，按下式计算展向相关系数：
Figure 348954DEST_PATH_IMAGE048
其中，
Figure 749979DEST_PATH_IMAGE049
为展向相关性系数，
Figure 538114DEST_PATH_IMAGE050
为展向位置，
Figure 239354DEST_PATH_IMAGE051
为参考展向位置 ，
Figure 932372DEST_PATH_IMAGE052
为相应展 向位置的流向速度的功率谱密度函数，
Figure 351852DEST_PATH_IMAGE045
为分析频率。
图5为基于流向速度的展向相关系数分布图，从图5可以看出在参考展向位置附近与212.5hz,展向相关系数较大。
可选地，设定相应展向相关系数阈值
Figure 421308DEST_PATH_IMAGE049
=0.3，从参考展向位置
Figure 101688DEST_PATH_IMAGE053
向两侧取到
Figure 981919DEST_PATH_IMAGE049
= 0.3的展向位置
Figure 875314DEST_PATH_IMAGE054
、
Figure 499193DEST_PATH_IMAGE055
，展向相关长度
Figure 893134DEST_PATH_IMAGE056
为两者间距：
Figure 209846DEST_PATH_IMAGE058
。
本实施例的展向相关长度如图6所示，可以看出在212.5hz处，展向相关长度有个极大值，这也对应着上游圆柱涡脱落的一阶谐频。
考虑基于流向速度的展向分布进行展向相关性修正，修正后的声压级为：
Figure 462973DEST_PATH_IMAGE059
；
其中
Figure 139811DEST_PATH_IMAGE060
为修正后的声压级，L为模型展向长度。
在本实施例中，所选取观测者的位置离声源的距离为0.15m，在该处布置TC4013-1水听器，采样率为10kHz，采样时长60s，所采集的噪声用于与PIV预测噪声对比。
将基于PIV预测的噪声声压级、基于PIV预测并进行展向相关性修正的噪声声压级、水听器采集的噪声声压级和水听器采集的背景噪声声压级进行比较，其结果如图7所示。
由图中结果可知，在100-300hz的频段内，水听器采集的双圆柱的特征频率及其一阶谐频特征明显，与PIV预测噪声基本一致。但是在300hz以上的频段，水听器采集噪声与PIV预测噪声误差较大，这是由于本实施例使用的水洞并非消音水洞，采集信号在该频段信噪比低。也说明了传统噪声测量容易引入其他噪声的缺点。
可以从图7看出，无论是水听器采集的噪声信号，还是基于PIV预测的噪声信号，双圆柱绕流噪声都存在一个112.5Hz的单音噪声；该单音噪声频率与上游圆柱脱落的涡撞击下游圆柱的频率对应，说明双圆柱绕流噪声主要由上游圆柱脱落的涡撞击下游圆柱引起。所获得的流场信息与声场信息之间的对应关系，为后续降噪措施提供理论指导。
由此，本实施例基于简单的平面PIV测量，完成了考虑流场三维效应的特定声源——双圆柱绕流的流致噪声预测，将双圆柱流场信息与流致噪声信息耦合。同时，也附带求解出了流场的时间解析压力信息，能进一步地对该流动进行动力学分析。
以上所述，仅为本发明较佳实施例而已，并非用于限定本发明的保护范围，发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种等同变换，这些等同变换均属于本发明的保护范围。

Claims (9)

1.一种基于平面PIV的预测指定声源远场流致噪声的方法，其特征在于，所述方法包括如下步骤：
1）通过平面时间解析PIV测得流场中的高频速度信息；
2）基于流场高频速度信息,通过欧拉法求得流场压力梯度，构建并求解泊松方程，重构流场高频压力信息；
3）通过Curle声比拟方法求得远场声压，获得远场流致噪声信息；
4）通过展向平面时间解析PIV测量，获得流场流向速度的展向分布信息；
5）对流向速度进行展向相关性分析，基于流向速度展向相关性修正流致噪声，得到考虑三维效应的远场流致噪声信息。
2.根据权利要求1所述的基于平面PIV的预测指定声源远场流致噪声的方法，其特征在于，所述的步骤1）具体为：利用同步器将高频双腔激光器和CMOS高速相机同步以采集和储存流场粒子图像，采集到的粒子图像通过PIV分析处理得到流场的高频速度信息。
3.根据权利要求1所述的基于平面PIV的预测指定声源远场流致噪声的方法，其特征在于，所述的步骤2）具体为：
通过欧拉法求得流场压力梯度
Figure 656538DEST_PATH_IMAGE002
，对压力梯度
Figure 534364DEST_PATH_IMAGE002
取散度得到拉普拉斯算子
Figure 150153DEST_PATH_IMAGE004
构建 泊松方程；给定方程的边界条件，通过迭代法求解各个时刻速度场对应的泊松方程，得到时 间解析压力信息，重构流场高频压力信息。
4.根据权利要求1所述的基于平面PIV的预测指定声源远场流致噪声的方法，其特征在于，步骤3）所述的通过Curle声比拟方法求得远场声压为：基于时间解析PIV所测得的高频速度信息与重构的高频压力信息计算Curle声比拟模型中的声源项，得到远场声压；计算各个时刻速度场对应的声压，求得时序声压信息，进一步求得流致噪声在频域下的声压级。
5.根据权利要求4所述的基于平面PIV的预测指定声源远场流致噪声的方法，其特征在于，所述的远场声压的计算公式为：
Figure 630682DEST_PATH_IMAGE006
其中,等号右边第一项和第二项分别为偶极子和四级子，
Figure 143091DEST_PATH_IMAGE007
为时序声压信号，
Figure 801605DEST_PATH_IMAGE008
、
Figure 888379DEST_PATH_IMAGE010
表 示以流场声源为坐标原点，观测者为终点位移的分量；y为以流场声源中心为坐标原点，各 个声源点的坐标；
Figure 657752DEST_PATH_IMAGE011
表示观察者与声源的距离；
Figure 552895DEST_PATH_IMAGE012
为声速；t为时间，
Figure 365999DEST_PATH_IMAGE013
为声信 号的产生时间，
Figure 690801DEST_PATH_IMAGE014
为当前声信号的接收时间；
Figure 527781DEST_PATH_IMAGE015
，
Figure 621639DEST_PATH_IMAGE016
、
Figure 746590DEST_PATH_IMAGE018
分别为来流压强、密度，P、
Figure 807956DEST_PATH_IMAGE019
为当地流体的压强，密度，
Figure 653552DEST_PATH_IMAGE020
为克罗内克符号；
Figure 116763DEST_PATH_IMAGE021
、
Figure 22402DEST_PATH_IMAGE023
为速 度分量；
Figure 573994DEST_PATH_IMAGE024
为固壁面法向量；S为面积,
Figure 613494DEST_PATH_IMAGE025
为固壁面表面积；V为体积，
Figure 947524DEST_PATH_IMAGE026
表示全流场体 积;
Figure 273332DEST_PATH_IMAGE027
为lighthill张量。
6.根据权利要求4所述的基于平面PIV的预测指定声源远场流致噪声的方法，其特征在于，所述的求得流致噪声在频域下的声压级，具体为：通过下式求得流致噪声在频域下的声压级：
Figure 794443DEST_PATH_IMAGE028
；
其中SPL为声压级，
Figure 762268DEST_PATH_IMAGE029
为时序声压信号
Figure 685224DEST_PATH_IMAGE030
的自功率谱密度，
Figure 325809DEST_PATH_IMAGE032
为分析频率，
Figure 583484DEST_PATH_IMAGE033
为 参考压力。
7.根据权利要求1所述的基于平面PIV的预测指定声源远场流致噪声的方法，其特征在于，所述步骤4)具体为：基于流向PIV测试实验平台，调整光学镜组与相机，将PIV测试区域改为主要声源附近的展向区域，对该区域进行时间解析PIV测量获取流向速度的高频展向分布信息。
8.根据权利要求1所述的基于平面PIV的预测指定声源远场流致噪声的方法，其特征在于，所述的步骤5）具体为：
5.1）选取主要声源临近的流向位置；
5.2）按时序提取该位置所有展向位置速度信息，得到该流向位置的速度展向分布时序信息；
5.3）计算该流向位置基于速度的展向相关性系数；
5.4）设定相应展向相关性系数阈值
Figure 840153DEST_PATH_IMAGE034
，从参考展向位置
Figure 866883DEST_PATH_IMAGE035
向两侧取到
Figure 550805DEST_PATH_IMAGE036
的展向位置
Figure 295776DEST_PATH_IMAGE037
、
Figure 218121DEST_PATH_IMAGE038
，得到展向相干长度：
Figure 381249DEST_PATH_IMAGE039
；
5.5）将直接由Curle声比拟得到的声压级减去
Figure 485340DEST_PATH_IMAGE040
，其中L为模型展长，完成展向相 关性修正。
9.根据权利要求8所述的基于平面PIV的预测指定声源远场流致噪声的方法，其特征在于，所述的步骤5.3）中，按照下式计算该流向位置基于速度的展向相关性系数，
Figure 468339DEST_PATH_IMAGE041
其中，
Figure 50499DEST_PATH_IMAGE042
为展向相关性系数，
Figure 192768DEST_PATH_IMAGE044
为展向位置，
Figure 952913DEST_PATH_IMAGE035
为参考展向位置 ，
Figure 663687DEST_PATH_IMAGE045
为相应展向位置 的流向速度的功率谱密度函数，
Figure 65850DEST_PATH_IMAGE046
为分析频率。
CN202211507962.7A 2022-11-29 2022-11-29 基于平面piv的预测指定声源远场流致噪声的方法 Active CN115541006B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211507962.7A CN115541006B (zh) 2022-11-29 2022-11-29 基于平面piv的预测指定声源远场流致噪声的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211507962.7A CN115541006B (zh) 2022-11-29 2022-11-29 基于平面piv的预测指定声源远场流致噪声的方法

Publications (2)

Publication Number Publication Date
CN115541006A true CN115541006A (zh) 2022-12-30
CN115541006B CN115541006B (zh) 2023-04-07

Family

ID=84722141

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211507962.7A Active CN115541006B (zh) 2022-11-29 2022-11-29 基于平面piv的预测指定声源远场流致噪声的方法

Country Status (1)

Country Link
CN (1) CN115541006B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116070550A (zh) * 2023-03-07 2023-05-05 浙江大学 一种改进的基于时间解析piv的重构流场压力场方法
CN116777772A (zh) * 2023-06-07 2023-09-19 浙江大学 一种自监督的修复前视声呐图像卷帘门效应的方法
CN117556720A (zh) * 2023-10-19 2024-02-13 上海慕灿信息科技有限公司 一种基于Lighthill声类比的气动声学频域求解方法

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05282282A (ja) * 1992-04-01 1993-10-29 Fujita Corp 環境騒音の汎用予測システム
JP2004156968A (ja) * 2002-11-05 2004-06-03 Railway Technical Res Inst 空力音源探査システム及び空力音源探査方法
JP2008064692A (ja) * 2006-09-08 2008-03-21 Mitsuba Corp 音源探査方法および装置
US20080103739A1 (en) * 2006-10-20 2008-05-01 Kawasaki Jukogyo Kabushiki Kaisha Jet noise modeling method, jet noise analyzing method, and aircraft designing method
CN104239625A (zh) * 2014-09-05 2014-12-24 中国矿业大学 一种基于修正流体运动方程线性迭代的稳态求解方法
JP2015210185A (ja) * 2014-04-25 2015-11-24 国立大学法人 新潟大学 空力騒音の音源特定装置及び空力騒音の音源特定方法
US20170219998A1 (en) * 2016-02-01 2017-08-03 Regents Of The University Of Minnesota System and method for digital inline holography
US20200151730A1 (en) * 2020-01-14 2020-05-14 Joseph Carlo Pastrana Mathematical Constant Pi Dynamic-Hybrid CVV Authentication Method for Credit Cards
CN111608927A (zh) * 2020-05-19 2020-09-01 中国船舶科学研究中心 一种离心泵流场-压力场-声场同步测量的装置和方法
CN112446179A (zh) * 2020-12-10 2021-03-05 华中科技大学 一种基于变分光流模型的流体速度测量方法
CN112763180A (zh) * 2021-04-08 2021-05-07 中国空气动力研究与发展中心低速空气动力研究所 声学风洞内高速列车受电弓模型气动噪声快速预测方法
CN113378437A (zh) * 2021-06-16 2021-09-10 山东大学 一种流固耦合作用下轴流风机噪声仿真预测方法及装置
CN114216666A (zh) * 2021-12-14 2022-03-22 中国运载火箭技术研究院 一种仿生潜航器摆动推进性能一体化测试平台

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05282282A (ja) * 1992-04-01 1993-10-29 Fujita Corp 環境騒音の汎用予測システム
JP2004156968A (ja) * 2002-11-05 2004-06-03 Railway Technical Res Inst 空力音源探査システム及び空力音源探査方法
JP2008064692A (ja) * 2006-09-08 2008-03-21 Mitsuba Corp 音源探査方法および装置
US20080103739A1 (en) * 2006-10-20 2008-05-01 Kawasaki Jukogyo Kabushiki Kaisha Jet noise modeling method, jet noise analyzing method, and aircraft designing method
JP2015210185A (ja) * 2014-04-25 2015-11-24 国立大学法人 新潟大学 空力騒音の音源特定装置及び空力騒音の音源特定方法
CN104239625A (zh) * 2014-09-05 2014-12-24 中国矿业大学 一种基于修正流体运动方程线性迭代的稳态求解方法
US20170219998A1 (en) * 2016-02-01 2017-08-03 Regents Of The University Of Minnesota System and method for digital inline holography
US20200151730A1 (en) * 2020-01-14 2020-05-14 Joseph Carlo Pastrana Mathematical Constant Pi Dynamic-Hybrid CVV Authentication Method for Credit Cards
CN111608927A (zh) * 2020-05-19 2020-09-01 中国船舶科学研究中心 一种离心泵流场-压力场-声场同步测量的装置和方法
CN112446179A (zh) * 2020-12-10 2021-03-05 华中科技大学 一种基于变分光流模型的流体速度测量方法
CN112763180A (zh) * 2021-04-08 2021-05-07 中国空气动力研究与发展中心低速空气动力研究所 声学风洞内高速列车受电弓模型气动噪声快速预测方法
CN113378437A (zh) * 2021-06-16 2021-09-10 山东大学 一种流固耦合作用下轴流风机噪声仿真预测方法及装置
CN114216666A (zh) * 2021-12-14 2022-03-22 中国运载火箭技术研究院 一种仿生潜航器摆动推进性能一体化测试平台

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ALASTAIR H. MOORE: "《 Direction of Arrival Estimation in the Spherical Harmonic Domain Using Subspace Pseudointensity Vectors》" *
LIDE FANG: "《Dual-modality UDV-PIV system for measurement of solid-liquid flow in sewage facilities》" *
PAUL ROGER LEINAN: "《Comparison of ultrasound vector flow imaging and CFD simulations with PIV measurements of flow in a left ventricular outflow trackt phantom - Implications for clinical use and in silico studies 》" *
李晓辉: "《层析粒子图像测速技术研究进展》" *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116070550A (zh) * 2023-03-07 2023-05-05 浙江大学 一种改进的基于时间解析piv的重构流场压力场方法
CN116777772A (zh) * 2023-06-07 2023-09-19 浙江大学 一种自监督的修复前视声呐图像卷帘门效应的方法
CN117556720A (zh) * 2023-10-19 2024-02-13 上海慕灿信息科技有限公司 一种基于Lighthill声类比的气动声学频域求解方法

Also Published As

Publication number Publication date
CN115541006B (zh) 2023-04-07

Similar Documents

Publication Publication Date Title
CN115541006B (zh) 基于平面piv的预测指定声源远场流致噪声的方法
Ghaemi et al. PIV-based pressure fluctuations in the turbulent boundary layer
Van Oudheusden PIV-based pressure measurement
Oguma et al. Measurement of sound source distribution around a circular cylinder in a uniform flow by combined particle image velocimetry and microphone technique
Henning et al. Simultaneous particle image velocimetry and microphone array measurements on a rod-airfoil configuration
Lodermeyer et al. Aeroacoustic analysis of the human phonation process based on a hybrid acoustic PIV approach
CN103698554A (zh) 一种流场实时精确测量系统及方法
Takaishi et al. Method of evaluating dipole sound source in a finite computational domain
Gui et al. Biases of PIV measurement of turbulent flow and the masked correlation-based interrogation algorithm
CN107843406B (zh) 基于脉动压力相关函数的空腔模态波运动特征确定方法
CN116577037B (zh) 一种基于非均匀频率语谱图的风管泄漏信号检测方法
Breakey et al. Aeroacoustic source analysis using time-resolved PIV in a free jet
Wu et al. Direct numerical simulation of the self-noise radiated by the installed controlled-diffusion airfoil at transitional Reynolds number
Wang et al. Identification of engine faults based on acoustic emission signals using a 1DCNN-ViT ensemble model
Shi et al. A modified cross-correlation algorithm for PIV image processing of particle-fluid two-phase flow
Roozeboom et al. Unsteady PSP measurements on a rectangular cube
Sanders et al. Slat noise measurements in open-jet, hard-wall and hybrid wind tunnel test sections
JP2004156968A (ja) 空力音源探査システム及び空力音源探査方法
Chu et al. High-resolution localization of rotating acoustic sources: An experimental investigation and axial fan application
Hileman et al. Exploring noise sources using simultaneous acoustic measurements and real-time flow visualizations in jets
JP3970263B2 (ja) 空力音源探査システム及び空力音源探査方法
Zhang et al. Aeroacoustic analysis of an airfoil with Gurney flap based on time-resolved particle image velocimetry measurements
Kouchi et al. Extracting dominant turbulent structures in supersonic flow using two-dimensional Fourier transform
Zanon et al. Experimental study on the aeroacoustic performance of an axial fan for HVAC
Ragni et al. Jet noise predictions by time marching of single-snapshot tomographic PIV fields

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