买专利卖专利找龙图腾,真高效! 查专利查商标用IPTOP,全免费!专利年费监控用IP管家,真方便!
申请/专利权人:中国空气动力研究与发展中心设备设计与测试技术研究所
摘要:本发明公开了一种风洞内部结构件噪声辐射强度计算方法,涉及噪声辐射强度计算领域,所述方法包括:步骤1:建立湍流运动混合模拟模型,利用湍流运动混合模拟模型进行流场仿真模拟获得风洞回路流场仿真结果;步骤2:基于风洞回路流场仿真结果,利用格林函数积分对闵可夫斯基时空中的能量与动量关系方程进行求解获得每一声源时刻流场的气动数据;步骤3:基于给定的接收点时刻t,计算获得接收点时刻t对应的声源时刻τ;步骤4:基于声源时刻τ和坐标信息,从每一声源时刻流场的气动数据中查询获得对应的噪声数据;本方法能够准确且高效的计算获得风洞内部结构件噪声辐射强度。
主权项:1.一种风洞内部结构件噪声辐射强度计算方法,其特征在于,所述方法包括:步骤1:建立湍流运动混合模拟模型,利用湍流运动混合模拟模型进行流场仿真模拟获得风洞回路流场仿真结果;步骤2:基于风洞回路流场仿真结果,利用格林函数积分对闵可夫斯基时空中的能量与动量关系方程进行求解获得每一声源时刻流场的气动数据;步骤3:基于给定的接收点时刻t,计算获得接收点时刻t对应的声源时刻τ;步骤4:基于声源时刻τ和坐标信息,从每一声源时刻流场的气动数据中查询获得对应的噪声数据;所述湍流运动混合模拟模型的建立方式为:基于RANS长度尺度和LES长度尺度混合构成混合湍流长度尺度;基于混合湍流长度尺度建立湍流运动混合模拟模型;利用格林函数积分对闵可夫斯基时空中的能量与动量关系方程进行求解,具体为: ;其中,为声压辐射强度,t为接收点时刻,为接收点在t时刻位置矢量,为时刻t和位置矢量对应的单极子声源项,为时刻t和位置矢量对应的偶极子声源项,为时刻t和位置矢量对应的四极子声源项; 的计算方式为: ;其中,;其中,f为计算频率,r为计算半径,ρ0为空气密度,Ui为质点振速,为矢量方向,Un为特定方向上的质点振速,为Un对声源时刻τ求导数,为Ui对声源时刻τ求导数,M为介质的质量矩阵,Mr为半径r球面上介质质量矩阵,为Mr对声源时刻τ求导数,为离散点i对应的计算半径的矢量,为Mi对声源时刻τ求导数,Mi为不同频率下的刚度矩阵,S为积分域,c为声速; 的计算方式为: ;其中,;其中,f为计算频率,c为声速,r为计算半径,Mr为半径r球面上介质质量矩阵,为Lr对声源时刻τ求导数,Li为不同频率的刚度矩阵,Lr为不同计算半径上的刚度矩阵,S为计算域,LM为第十一中间量,为Mr对声源时刻τ求导数,M为介质的质量矩阵,为离散点i对应的计算半径的矢量,为Li对声源时刻τ求导数,Mi为不同频率下的风度矩阵; 的计算方式为: ; ; ; ;其中,;其中,f为计算频率,c为声速,r为计算半径,K1为第十二中间量,K2为第十三中间量,K3为第十四中间量,V为速度域,为Trr对声源时刻τ求二阶导数,Mr为半径r球面上介质质量矩阵,为Mi对声源时刻τ求二阶导数,Trr为第十五中间量,为Mr对声源时刻τ求导数,为Trr对声源时刻τ求导数,为Tii对声源时刻τ求导数,Tii为在空间对角线上的总温向量,为TMr对声源时刻τ求导数,TMr为第十六中间量,为Mi对声源时刻τ求导数,Mi为不同频率下的风度矩阵,M为介质的质量矩阵,TMM为第十七中间量,i为和j为离散点的记号,Tij为总温空间分布矩阵,Mj为空间某一截面上的介质质量矩阵,为离散点j对应的计算半径的矢量,为Tij对声源时刻τ求导数,为离散点i对应的计算半径的矢量。
全文数据:
权利要求:
百度查询: 中国空气动力研究与发展中心设备设计与测试技术研究所 一种风洞内部结构件噪声辐射强度计算方法
免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。