格子ボルツマン法を基とした流体音の直接数値計算手法の確立

基于格子玻尔兹曼法的流体声直接数值计算方法的建立

• 批准号: 17760143
• 负责人: 片岡 武
• 金额: $ 2.18万
• 依托单位: Kobe University
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for Young Scientists (B)
• 财政年份: 2005
• 资助国家: 日本
• 起止时间: 2005 至 2006
• 项目状态: 已结题
• 来源: https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-17760143/
• 关键词: 流体音; 直接計算; 数値計算法

昨年度は、差分格子ボルツフン法による空力音の直接計算の計算格子として新たに有限体積法を導入することで精度向上を測った。つまり、流体の圧縮性を考慮した格子ボルツマン法である熱流体モデルに対して、有限体積法を適用し、エオルス音の発生メカニズムについて検証するとともに、層流流れの基本的な音場を直接計算した。本年度は、流れが乱流である場合をとり扱った。実際、空力騒音の抑制をはじめとする工学的問題において乱流場の音の現象を捉えることは非常に重要になってきている。ただ最初から直接計算によって音波の乱流流れを捉えることは難しいため、まずは乱流モデルを導入することにした。スマゴリンスキーモデルを差分格子ボルツマン法に組み込んだモデルを提案し、乱流および空力音の数値計算を行った。計算対象としては、角柱周りの乱流計算および空力音の計算、自動車のフロントピラー部をモデル化したウェッジボックス周りの流れの数値計算をおこなった。その結果、前者に対しては平均流速や変動速度などの乱流統計量に関しては実験結果と若干異なる結果となったが、角柱から発生する音波の周波数特性については、変動のピークにおけるストローハル数や広帯域における変動において実験値との定性的な一致が確認できた。また後者のウェッジボックス周りの乱流計算に対しては、流れの再現性について乱れの非等方性を有する旋回流の再現ができるなど、良好な結果が得られた。しかし、ウェッジボックス側面部の圧力変動の周波数特性は特に高周波数俗において再現できなかった。この原因としては格子解像度の不足が考えられ、今後格子数を増やして計算を行うことで改善できると考えられる。以上の結果より、乱流場における音の直接数値計算コードを完成させる土台ができたものと考えている。

In the past, the difference lattice method was introduced into the calculation of spatial force and sound. The lattice method is suitable for thermal fluids, and the finite volume method is suitable for laminar flow. This year, the flow is turbulent. The problem of noise suppression in real time and space is very important. The initial calculation of the acoustic turbulence was carried out. The difference between the number of points and the number of points Calculation of the number of objects, the turbulence of the prism cycle, the calculation of the air force sound, the calculation of the number of objects in the automatic vehicle cycle The results of the former are related to the average velocity, dynamic velocity and turbulence statistics. The results of the former are related to the number of cycles of the acoustic wave generated by the prism. The results of the former are related to the number of cycles of the acoustic wave generated by the prism. The results of the former are related to the number of cycles of the acoustic wave generated by the prism. The results of the former are related to the turbulence statistics of the acoustic wave generated by the prism. The results of the calculation of turbulence in the latter part of the equation are satisfactory. The frequency characteristics of the pressure variation at the bottom of the tank are different from those at the high frequency. The reason for this is that the grid resolution is insufficient, and the number of grids in the future is increased. The above results show that the direct numerical value calculation of chaotic flow is completed.