[发明专利]基于非均匀入流中螺旋桨空化噪声数值预报的特征提取方法

摘要

本发明公开了一种基于非均匀入流中螺旋桨空化噪声数值预报的特征提取方法，步骤包括首先，对螺旋桨计算域进行网格划分，检查网格质量并定义边界条件；接下来，在CFD软件中，设置计算模型，进行稳态迭代计算淌水性能参数和入流口速度验证模型准确性；然后，在CFD软件中，将稳态计算作为非稳态计算的初始值进行非稳态迭代计算，并通过后处理显示螺旋桨片空化周期形态及记录片空化面积变化；最后，根据单空泡辐射噪声理论由螺旋桨片空化面积计算螺旋桨空化辐射噪声，进行特征提取。

主权项

一种基于非均匀入流中螺旋桨空化噪声数值预报的特征提取方法，其特征在于，步骤包括：首先，对螺旋桨计算域进行网格划分，检查网格质量并定义边界条件；接下来，在CFD软件中，设置计算模型，进行稳态迭代计算淌水性能参数和入流口速度验证模型准确性；然后，在CFD软件中，将稳态计算作为非稳态计算的初始值进行非稳态迭代计算，并通过后处理显示螺旋桨片空化周期形态及记录片空化面积变化；最后，根据单空泡辐射噪声理论由螺旋桨片空化面积计算螺旋桨空化辐射噪声，进行特征提取；所述对螺旋桨计算域进行网格划分，检查网格质量并定义边界条件；具体步骤如下：第一步，对螺旋桨计算域进行网格划分：设定坐标原点在螺旋桨中心点处，螺旋桨的旋转轴为x轴，螺旋桨的直径为D，假定x正方向为下游，x负方向为上游；建立计算动域和计算静止域；进行网格划分，检查网格质量；定义流体计算域和边界条件；第二步，采用一阶差分插值计算非均匀入流速度面的速度值：实测非均匀入流速度面的三向速度值，测试点的位置分布在r/R＝0.2,0.4,0.6,0.8,1.0…0.2Mtest，其中R为螺旋桨半径，r为测试点的半径，Mtest为测试半径的个数；入流速度面1的面网格划分采用结构网格划分，径向的节点数为M，一般情况下，Mtest＜M，在径向方向进行一阶插值计算，径向插值的点数为Mr；第三步，导入计算网格并进行网格相关操作：将网格文件读入到CFD软件中进行网格检查，确保最小网格体积大于零，否则重新划分网格；进行计算域尺寸的调整，调整网格的尺寸比例，使最后尺寸符合实际模型的大小；将网格进行交换和光滑处理，依次按交换系数从小到大进行光滑和交换操作，直到每个交换系数的交换网格数目为零；读入非均匀入流面速度的文本文件；第四步，计算模型设置：在CFD软件中进行计算模型设置，定义求解器，按照默认的设置；选择K‑epsilon模型，在K‑epsilon模型下保存默认选项Stadard，在Near‑Wall Treatment选项下选择标准壁面函数；流体介质的选取，选取water‑liquid和water‑vapor，根据实际值设置流体介质的密度和粘性系数；设置操作环境，在Operating Pressure中根据实际值写入环境压力；将螺旋桨的旋转速度的单位选为rpm；设置边界条件，包括选择流体动域的坐标系，稳态计算中选取MRF，非稳态计算中选取Moving Mesh，定义动域旋转轴为x轴，根据实际值设定选准速度，速度方向遵循右手定则，默认流体静止域为静止坐标系，非均匀入流速度面中速度采用笛卡尔坐标系，分别选择profile读入文本文件中的vx，vy，vz，静止域外围速度采用笛卡尔坐标系，x方向的速度值根据实际的入流速度进行设定，压力出口选择静压力为零，定义螺旋桨的桨叶和桨毂壁面速度无滑移，粗糙系数为0；定义交界面；第五步，设置求解参数、初始化以及收敛条件：设置松弛因子，其中vapor，Turbulent Kinetic Energy，Specifiec Dissipation Rate，Turbulent Viscosity的松弛因子为λ，其它保持不变，定义差分方程形式，其中Pressure为Standard，其它的均为First Order Upwind形式；将入流口速度面1的速度作为初始值；设置各参数的残差。