基于螺旋桨尾流压力脉动计算的空化噪声特征估计方法技术

本发明专利技术公开了基于螺旋桨尾流压力脉动计算的空化噪声特征估计方法，属于水声目标特征提取领域。基于螺旋桨尾流压力脉动计算的空化噪声特估计取方法，包括以下步骤：(1)备用网格生成并导入计算程序后生成算例文件，(2)空化模型和湍流模型设定，(3)数值计算参数设定，(4)数值计算，(5)数值方法可靠性验证及网格确定，(6)空化尾流压力脉动非定常数值计算，(7)压力脉动信号功率谱变换及低频线谱幅值提取，(8)线谱特征估计及分析。本发明专利技术将现代流体力学、空泡动力学和信号处理领域中相关研究成果引入水下目标的噪声特征分析，体现多学科和多领域的交叉性。

【技术实现步骤摘要】
基于螺旋桨尾流压力脉动计算的空化噪声特征估计方法
本专利技术涉及水声目标特征提取领域，具体地说，涉及基于螺旋桨尾流压力脉动计算的空化噪声特征估计方法。
技术介绍
螺旋桨噪声是船舶三大噪声源之一，包含了目标推进器种类信息和结构特征，这些特征宽容性强，具有较好的可分性，是识别水下目标的主要特征和重要依据。而空化一旦出现，空化噪声就成为螺旋桨主要噪声。这些目标源噪声由于被海洋环境噪声干扰和在复杂的水声信道传播中而产生畸变，使得被动声纳所接收到的噪声信号特征不明显，信噪比降低。因此传统的以信号处理方法提取噪声特征，进行水下目标识别越来越困难。进一步挖掘螺旋桨噪声本质特征是水下目标识别急待解决的问题。对于采用信号处理技术对实测的螺旋桨噪声进行特征提取方面的研究，国外的学者很早就已经开始了。早在1971年Whalen就已经提出了最大似然调制接收机。随着这一技术的发展，高阶谱、AR谱、双重谱和小波分析等时频处理方法，以及分形、混沌、极限环和模态分解等非线性处理方法，都在螺旋桨噪声特征提取中广为尝试。近年来，李启虎等学者采用理论分析和数值仿真研究了强干扰背景噪声下单频信号分量检测方法和检测系统性能。南京大学鲍菲等将经验模型分解法(empiricalmodedecomposition)和奇异值分解法(singularvaluedecomposition)相结合，从强干扰背景噪声中提取螺旋桨的空化噪声调制成分。现代信号处理方法对背景噪声下的实测噪声信号特征进行提取，取得了不错的效果。但是对强干扰背景噪声，由于实测信号中缺乏机理特征，这一方法适应能力较不高。由此，一些学者开展了基于模型的噪声特征分析及方法研究。陶笃纯将噪声调制包络作为有相同形状、相等重复周期、随机幅度，具有成组结构的脉冲性随机过程处理。并从舰船辐射噪声调制包络的功率谱密度和自相关函数中提取与舰船各种物理属性有关的丰富的节奏信息。蒋国健和林建恒等人利用指数衰减形随机脉冲序列的理论模型来分析舰船螺旋桨空泡噪声，得到螺旋桨空化噪声谱。近年来，史广智等学者针对螺旋桨叶片数识别问题，建立空化噪声信号模型。并对双螺旋桨舰船噪声包络建模，研究双桨目标调制谱谐波族特征的结构问题，进一步采用模型特征提取技术，研究基于模型匹配的噪声特征精细分析方法。除了叶频特征，这些模型没有考虑螺旋桨几何形状和工况等参数，很难体现空化噪声的机理特征。螺旋桨空化是空化噪声的直接声源，并且螺旋桨空化尾流是空化噪声重要的传播途径。由于尾流受到螺旋桨周期性转动节拍的作用，具有周期性脉动特征。这些特征反映了螺旋桨工况和几何形状等特征信息。同时，螺旋桨旋转节拍对其辐射的空化噪声有明显的振幅调制作用，其功率谱的线谱特征也反映包括螺旋桨工况和几何形状等特征信息在内的螺旋桨节奏信息。因此，由于同样受到螺旋桨桨叶的节拍作用，螺旋桨空化尾流与空化噪声具有特征相关性，其特征都反映了螺旋桨工况参数和几何形状参数。由于目前螺旋桨空化噪声的声学机理研究还很不完善，因此本专利技术从空化噪声的本源即空化尾流入手来阐述其噪声特征的一种预报方法。对于螺旋桨空化尾流国内外有不少学者进行了研究。意大利船模水池实验室的Francesc等利用RANS、LES和BEM方法分别对空化和非空化条件下E779A螺旋桨尾流场进行数值模拟。瑞典Rickard和Goran基于混合两相流模型，利用隐式LES方法和Kunz空化模型模拟了E779A螺旋桨在非均匀流场中空化的动态行为，对中小尺度的流场结构和螺旋桨梢涡空化的模拟较为成功。清华大学季斌等学者利用Rayleigh–Plessete方程和k-ωShearStressTransport(SST)湍流模型对高侧斜螺旋桨均匀和非均匀入流的空化尾流进行了数值模拟。片空化和梢涡空化被较好地预报，同时空化诱导的尾流场压力脉动特征与螺旋桨轴频叶频特征一致。海军工程大学杨琼方对空化模型和湍流模型在螺旋桨空化模拟进行评估分析，选择改进Sauer空化模型和修正SSTk-ω湍流模型，较准确地预报螺旋桨空化斗图谱。对与七叶大侧斜桨的非均匀进流，分析了其空化引起的推力和力矩崩溃性能以及对叶背梢涡空化初生的影响，描述了空化推力和力矩的脉动特征、桨叶空化面积和空化形态随周向位置的变化，并给出了伴流中螺旋桨是否出现叶面片空化的区间划分。目前，国内外对空化尾流的研究主要侧重于某种桨模片空化的数值预报，对于空化与螺旋桨工况与几何形状之间的特征关系方面的研究较少，而用空化尾流来研究空化噪声特征则更少。另外，中国专利申请号ZL201310538724.7，文件也公开了一种基于非均匀入流中螺旋桨空化噪声数值预报的特征提取方法，步骤包括：首先，对螺旋桨计算域进行网格划分，检查网格质量并定义边界条件；接下来，在CFD软件中，设置计算模型，进行稳态迭代计算淌水性能参数和入流口速度验证模型准确性；然后，在CFD软件中，将稳态计算作为非稳态计算的初始值进行非稳态迭代计算，并通过后处理显示螺旋桨片空化周期形态及记录片空化面积变化；最后，根据单空泡辐射噪声理论由螺旋桨片空化面积计算螺旋桨空化辐射噪声，进行特征提取。该申请文件中所用方法将空化区域折算成球形体积，并得到球形体积半径，再将半径变化带入球形单空泡噪声辐射模型中，来预报空化噪声及其特征。由于螺旋桨空化与球形单空泡有很大不同，这种折算方法的准确性有待于进一步检验。本专利技术的方法则是利用螺旋桨空化尾流压力脉动与空化噪声之间的特征相关性来估计噪声特征。具体来说，空化尾流压力脉动信息含有螺旋桨工况和几何形状参数特征，而空化噪声也具有这一属性。因此，它们具有相同的本源关系，即是螺旋桨的旋转导致空化尾流并产生噪声，同时空化噪声还受到旋转桨叶的调制作用。空化尾流和噪声产生的根本原因是螺旋桨在流体中的转动。
技术实现思路
螺旋桨空化噪声主要特征有：1.螺旋桨空化是空化噪声的来源，螺旋桨空化噪声总是伴随着螺旋桨空化的出现而出现；2.螺旋桨空化尾流不仅是空化噪声的声源，还是空化噪声传播的重要载体；3.空化噪声的声强与空泡体积变化密切相关，特别是其在溃灭瞬间的体积变化最大，其辐射噪声也最强；4.螺旋桨空化体积变化和位置分布随着桨叶的旋转而具有周期性特征，使得空化噪声也具有周期性特征，并会反映到其噪声频谱分布上；5.空化尾流和空化噪声同时会受到螺旋桨桨叶转动节拍的调制作用；6.上述特征使得螺旋桨空化与其噪声的特征具有严密的本质相关性；7.片空化噪声一般分布在低频段，其频谱呈现线谱特征，而梢涡空化发出的噪声一般分布在中高频段，其频谱呈现连续特征。本专利技术的原理就是依据上述螺旋桨空化的主要特征，基于黏性多相流理论，利用现代计算流体力学方法对水下螺旋桨尾流场构建N-S方程，并结合湍流模型和空化模型对方程组进行数值求解，从而得到水下螺旋桨叶面周围汽相体积分数和尾流场中压力脉动等相关信息；再利用功率谱等信号处理方法对数值计算的流场信息数据的低频特征进行提取和分析；最后利用流场压力脉动与噪声之间的特征相关性对水下目标螺旋桨噪声特征进行估计和判断。虽然压力脉动是尾流场中力学参数而噪声声压是声学参数，它们的物理概念不同，但它们在某些方面的特征，如低频线谱幅值分布特征，又有一些共同点。这些共同点实质是反映了螺旋桨几何和工况参数特征，这里把它本文档来自技高网...

【技术保护点】
基于螺旋桨尾流压力脉动计算的空化噪声特征估计取方法，其特征在于：包括以下步骤：(1)备用网格生成并导入计算程序后生成算例文件：利用专业建模软件制作螺旋桨三维几何模型后导入网格生成软件，在网格划分软件中建立三种备选网格，这三种备选网格的计算域相同，速度入流边界距离螺旋桨中心为1D，D为螺旋桨直径，下游压力出口边界距离为5D，螺旋桨中心至侧面外围距离为2.5D，这三个网格的网格单元数量按照倍数逐渐增加，对网格中相邻边界的网格单元尺寸在边界点合理过渡，使得网格中所有体网格单元的skew都限定在0.9以内；(2)空化模型和湍流模型设定：采用全空化模型和重整化群湍流模型，并对其重要参数进行修正；(3)数值计算参数设定：对工况条件、边界条件和数值算法的相关参数进行设定；(4)数值计算：采用逐级分步骤的计算过程，在螺旋桨工况参数中，环境压力和入流速度可以直接设定到工况值，而螺旋桨转速采用分级增加，直到增加到预定工况值；先计算无空化模型流场分布，等到计算稳定后再打开空化模型；先对压力、密度、动量和汽相分数等参数进行一阶精度离散格式计算，计算稳定后，再将离散精度提高到二阶或QUCIK等，并采用并行计算技术来进行计算；(5)数值方法可靠性验证及网格确定：将典型工况下对螺旋桨桨的水动力参数和空化的数值计算结果与相关实验结果进行比较，以验证网格无关性和所采用数值方法的可靠性；并对数值计算结果中水动力参数和空化进行比较，当结果随着网格数量的增加而趋于稳定并与实验结果一致时，则选定满足条件中网格单元数量最少的网格作为下面数值计算的选定网格；否则适当增加网格数量，重复步骤1重新开始；(6)空化尾流压力脉动非定常数值计算：采用步骤5中的选定网格，对螺旋桨的尾流场在所需工况条件下进行非定常数值计算，在计算程序中对尾流场中某一特定位置(A点)压力脉动检测并保存其检测数据，同时对数据进行无量纲化；(7)压力脉动信号功率谱变换及低频线谱幅值提取：采用信号处理中快速傅立叶变换方法对流场中压力脉动等物理量和噪声信号数据进行功率谱变换，并对低频线谱幅值进行提取；再利用尾流场压力脉动特征与空化噪声被桨叶调制特征的相似性，建立从压力脉动的低频线谱幅值到噪声的低频线谱幅值的特征对应关系；(8)线谱特征估计及分析：将步骤7中低频线谱幅值一一对应到噪声信号功率谱的低频线谱幅值，作为对空化噪声信号低频线谱幅值分布特征的估计。...

【技术特征摘要】
1.基于螺旋桨尾流压力脉动计算的空化噪声特征估计方法，其特征在于：包括以下步骤：(1)备用网格生成并导入计算程序后生成算例文件：利用专业建模软件制作螺旋桨三维几何模型后导入网格生成软件，在网格划分软件中建立三种备选网格，这三种备选网格的计算域相同，速度入流边界距离螺旋桨中心为1D，D为螺旋桨直径，下游压力出口边界距离为5D，螺旋桨中心至侧面外围距离为2.5D，这三个网格的网格单元数量按照倍数逐渐增加，对网格中相邻边界的网格单元尺寸在边界点合理过渡，使得网格中所有体网格单元的skew都限定在0.9以内；(2)空化模型和湍流模型设定：采用全空化模型和重整化群湍流模型，并对其重要参数进行修正；(3)数值计算参数设定：对工况条件、边界条件和数值算法的相关参数进行设定；(4)数值计算：采用逐级分步骤的计算过程，在螺旋桨工况参数中，环境压力和入流速度可以直接设定到工况值，而螺旋桨转速采用分级增加，直到增加到预定工况值；先计算无空化模型流场分布，等到计算稳定后再打开空化模型；先对压力、密度、动量和汽相分数参数进行一阶精度离散格式计算，计算稳定后，再将离散精度提高到二阶或QUCIK，并采用并行计算技术来进行计算；(5)数值方法可靠性验证及网格确定：将典型工况下对螺旋桨桨的水动力参数和空化的数值计算结果与相关实验结果进行比较，以验证网格无关性和所采用数值方法的可靠性；并对数值计算结果中水动力参数和空化进行比较，当结果随着网格数量的增加而趋于稳定并与实验结果一致时，则选定满足条件中网格单元数量最少的网格作为下面数值计算的选定网格；否则适当增加网格数量，重复步骤1重新开始；(6)空化尾流压力脉动非定常数值计算：采用步骤5中的选定网格，对螺旋桨的尾流场在所需工况条件下进行非定常数值计算，在计算程序中对尾流场中某一特定位置A点压力脉动检测并保存其检测数据，同时对数据进行无量纲化；(7)压力脉动信号功率谱变换及低频线谱幅值提取：采用信号处理中快速傅立叶变换方法对流场中压力脉动物理量和噪声信号数据进行功率谱变换，并对低频线谱幅值进行提取；再利用尾流场压力脉动特征与空化噪声被桨叶调制特征的相似性，建立从压力脉动的低频线谱幅值到噪声的低频线谱幅值的特征对应关系；(8)线谱特征估计及分析：将步骤7中低频线谱幅值一一对应到噪声信号功率谱的低频线谱幅值，作为对空化噪声信号低频线谱幅值分布特征的估计。2.如权利要求1所述的基于螺旋桨尾流压力脉动计算的空化噪声特征估计方法，其特征在于：所述的步骤6中的空化尾流压力脉动非定常计算包括以下步骤：(6-1)导入步骤5中的选定网格生成算例文件；(6-2)空化模型和湍流模型设定；(6-3)数值计算参数设定；(6-4)数值计算；(6-5)压力脉动信号提取：对尾流场中某一特定位置A点压力脉动检测并保存其检测数据。3.如权利要求1所述的基于螺旋桨尾流压力脉动计算的空化噪声特征估计方法，其特征在于：所述的步骤1中的备选网格采用分区域混合网格划分方法：螺旋桨周围流场区域采用非结构网格方法划分，网格由桨毂到叶梢逐渐减小，叶梢处面网格为三角形，网格单元边长大小为0.001D，桨榖处单元为0.02D；在桨叶表面建立边界层网格；采用结构网格划分螺旋桨外围规则形状的计算域；对梢涡区域网格进行加密，同时桨叶表面采用边界层网格以提高对梢涡空化的预报精度，梢涡区域网格单元尺寸为0.001D，边界层网格共有4层，其相邻两层高度比为1.1，第一层网格单元高度为0.001D，使得无量纲参数20<y+<300。4.如权利要求1所述的基于螺旋桨尾流压力脉动计算的空化噪声特征估计方法，其特征在于：所述的步骤2中的全空化模型设定及其参数修正为：当p<pv时，蒸汽产生率为：

【专利技术属性】
技术研发人员：朱志峰，
申请(专利权)人：安徽工业大学，
类型：发明
国别省市：安徽;34