[发明专利]基于KIVA和SYSNOISE耦合内燃机的燃烧噪声预测分析方法

权利要求书

1.一种基于KIVA和SYSNOISE耦合内燃机的燃烧噪声预测分析方法，其特征在于，包括：

采用KIVA计算第N时间步长当前燃烧点源以及第N时间步长至初始时间步长内声波再传播对内燃机缸内网格单元的第一压力更新量；

采用SYSNOISE计算第N时间步长所述内燃机缸内壁面反射对内燃机缸内网格单元的第二压力更新量，包括：将内燃机缸内整体的流场网格与激励力所处的壁面作局部流固固体耦合，得到流固耦合的有限元模型；采用SYSNOISE分别计算各时间步长的每一个当前激励力对内燃机缸内所有网格单元产生的压力改变量；将各时间步长的当前激励力的所述压力改变量存储至本地；根据第N时间步长当前激励力在所述内燃机缸内网格单元的声波传播距离在所述压力改变量中提取第一压力改变量；根据所述第N时间步长至初始时间步长内各时间步长对应的声波继续传播距离在所述压力改变量中提取声波继续传播的第二压力改变量；迭加所述第一压力改变量与所述第二压力改变量得到第N时间步长内燃机缸内壁面反射对内燃机缸内网格单元的第二压力更新量；

迭加所述第一压力更新量与所述第二压力更新量得到第N时间步长内燃机缸内网格单元的压力总升高量；

根据所述各时间步长的所述压力总升高量确定所述内燃机缸内的压力曲线。

2.根据权利要求1所述的方法，其特征在于，所述采用KIVA计算第N时间步长当前燃烧点源以及第N时间步长至初始时间步长内声波再传播对内燃机缸内网格单元的第一压力更新量，包括：

处理器采用KIVA分别计算各时间步长的每一个当前燃烧点源对内燃机缸内所有网格单元产生的压力改变量；

所述处理器将各时间步长的当前燃烧点源的所述压力升高量存储至本地；

所述处理器根据第N时间步长当前燃烧点源在内燃机缸内网格单元的声波首次传播的距离在所述压力改变量中提取声波首次传播的第三压力改变量；

所述处理器根据第N时间步长至初始时间步长内各时间步长对应的声波继续传播的距离在所述压力改变量中提取声波继续传播的第四压力改变量；

所述处理器迭加所述第三压力改变量与所述声波继续传播的第四压力改变量得到第N时间步长内燃机缸内的第一压力更新量；

所述处理器根据各时间步长的所述压力更新量确定内燃机缸内的压力曲线。

3.根据权利要求2所述的方法，其特征在于，所述处理器根据第N时间步长当前燃烧点源在内燃机缸内网格单元的声波首次传播的距离在所述压力改变量中提取声波首次传播的第三压力改变量，包括：

处理器根据KIVA的化学反应子程序确定第N时间步长内当前燃烧点源的总个数；

所述处理器根据所述当前燃烧点源对应的网格单元中各个节点的空间坐标计算所述当前燃烧点源的第N时间步长内的空间位置坐标；

所述处理器根据小振幅三维声波方程计算所述第N燃烧点源第N时间步长的声速，并根据所述声速计算所述第N燃烧点源在第N时间步长内在空间传播的距离；

所述处理器根据所述当前燃烧点源的空间位置坐标以及所述声波传播的距离确定所述第N时间步长内所述当前燃烧点源影响的声内单元；

所述处理器根据所述声内单元在已存储的压力改变量中提取声波首次传播的第三压力改变量。

4.根据权利要求2所述的方法，其特征在于，所述处理器根据第N时间步长至初始时间步长内各时间步长对应的声波继续传播的距离在所述压力改变量中提取声波继续传播的第四压力改变量，包括：

所述处理器根据声速确定第N时间步长至初始时间步长中各时间步长声波再传播的距离；

所述处理器根据所述各时间步长当前燃烧点源空间位置坐标、所述当前燃烧点源传播的距离以及所述声波再传播距离确定各时间步长对应的声波再传播的声内单元；

所述处理器根据所述声波再传播的声内单元在已存储的压力改变量中提取声波继续传播的第四压力改变量。

5.根据权利要求2所述的方法，其特征在于，所述处理器迭加所述第三压力改变量与所述声波继续传播的第四压力改变量得到第N时间步长内燃机缸内的第一压力更新量之后，还包括：

所述处理器将第N时间步长的所述压力更新量对应的文件存储至本地，所述文件包括：所述第N时间步长对应的所述当前燃烧点源的声传播距离、所述声波传播的声内单元总个数、所述燃烧点源的空间位置坐标、所述内燃机缸内网格单元压力升高量、所述内燃机缸内网格单元的起始号码和所述内燃机缸内网格单元的终止号码。