ライトヒルの音響アナロジー

カテゴリ: 連成解析・音響連成 | 2026-10-01 改稿

理論:流れの方程式を音の方程式に書き直す

概要

🙋

ジェット機の音は、なぜあんなに大きいんですか?

🎓

高速の噴流が周りの空気と混ざるときの乱れが、音を出しているからだ。ライトヒルは1952年、流れの方程式を変形して、左辺を静止した空気の中の音の波の式、右辺を乱れた流れが作る音源の項に書き直した。この右辺の音源は四重極と呼ばれる形で、次元を調べると、出てくる音のパワーは速度の8乗に比例することがわかる。速度が2倍なら音は256倍、つまり24 dBも大きくなる。だから噴流の速度をわずかに下げるだけで、騒音は大きく下がる。現代のジェットエンジンは、大きなファンで多くの空気をゆっくり流して推力を出し、噴流の速度を下げて静かにしている。

ライトヒルの方程式

$$ \frac{\partial^2\rho'}{\partial t^2} - c_0^2\nabla^2\rho' = \frac{\partial^2 T_{ij}}{\partial x_i\,\partial x_j},\qquad T_{ij} = \rho u_i u_j + (p' - c_0^2\rho')\delta_{ij} - \tau_{ij} $$

ρ′ は密度の変動、c₀ は静止空気の音速、T_ij はライトヒル応力テンソル(流れのレイノルズ応力、エントロピーの変動、粘性応力からなる)。右辺の2階微分が四重極の音源を表す。

8乗則と音響効率

$$ W = K\,\frac{\rho_0 U^8 D^2}{c_0^5},\quad K \approx 3\times10^{-5},\qquad L_W = 10\log_{10}\frac{W}{10^{-12}\,\mathrm{W}},\qquad \eta = \frac{W}{\tfrac12\rho_0 U^3 A} \propto M^5 $$

K は実験で決まる係数、D は噴流の直径、A は噴流の断面積、M は噴流のマッハ数。亜音速の冷たい噴流に当てはまる近似。

Coffee Break よもやま話

ジェット旅客機の時代の始まりに

ライトヒルが音響アナロジーを発表した1952年は、世界初のジェット旅客機コメットが就航した年でもあった。ジェットエンジンの騒音は空港の周りで大きな問題となり、8乗則は、噴流の速度を下げることが最も効く対策であることを示した。これがのちに、多くの空気をファンで送る高バイパス比のターボファンエンジンにつながり、同じ推力でも騒音は大きく下がった。音響アナロジーは今もCFDで騒音を予測する方法の土台になっている。

計算例

例1:直径0.5 mの噴流の音響パワー

噴流速度マッハ数音響パワー音響パワーレベル噴流の運動エネルギー音響効率100 mでの音圧レベル
100 m/s0.290.02 W102.8 dB0.12 MW1.6×10⁻⁷51.8 dB
200 m/s0.584.9 W126.9 dB0.94 MW5.1×10⁻⁶75.9 dB
300 m/s0.87124 W140.9 dB3.18 MW3.9×10⁻⁵90.0 dB

100 mの音圧レベルは、すべての方向に均等に広がると仮定した値(実際の噴流音は後方斜め方向に強い)。速度を1割下げると−3.7 dB、2倍にすると+24 dB。

例2:音を計算するための格子(300 m/s)

ストローハル数周波数波長格子(1波長に10セル)
0.2120 Hz2.86 m286 mm
1.0600 Hz0.57 m57 mm
🙋

流れのエネルギーのたった0.004%しか音にならないのに、こんなに大きな音になるんですね。

🎓

噴流が運ぶエネルギーが3 MWと非常に大きいので、その0.004%でも124 Wになる。人の耳は10⁻¹² W/m²(1 W/m²の1兆分の1)という弱い音から聞き取れるので、124 Wは非常に大きな音だ。この小ささは、音の予測を難しくしてもいる。音の変動は流れの変動よりずっと小さいので、流れと音を同じ計算で直接求めると、数値の誤差に音が埋もれてしまう。そこで、乱れた流れの部分だけをLESなどで細かく計算して音源を求め、遠くへの伝わりはライトヒルの式やFW-H法などの積分で計算する分業(ハイブリッド法)がよく使われる。流れの計算でも、関心のある最も高い周波数の音源を解像できる格子と時間刻みが必要になる。

騒音予測の進め方

  1. 8乗則などの簡単な式で、騒音の大きさと速度の効果の目安をつかむ。
  2. 関心のある周波数から、格子の大きさと時間刻みを決める。
  3. 乱れの領域をLESなどの非定常計算で解き、音源を求める。
  4. FW-H法やライトヒルの積分で、遠くの観測点の音を計算する。
  5. 実測の騒音と、周波数ごと・方向ごとに比べる。
Coffee Break よもやま話

「排気の速さを1割下げたら」

ある工場の排気ダクトの出口から出る気流の騒音が、敷地の境界で基準をわずかに超えていた。消音器を追加する案もあったが、まず8乗則で見積もると、出口の速度を1割下げれば約3.7 dB下がり、基準を満たせることがわかった。出口の断面を約1割広げて速度を下げると、測定でもほぼ計算どおり騒音が下がった。大がかりな対策の前に、速度の効果を見積もることが役に立った例だ。

よくある間違い

間違いと対策

間違い起きること対策
流れと音を同じ格子で直接解く音が誤差に埋もれるハイブリッド法
RANSで音源を求める変動が得られないLESなど非定常
格子が粗い高い周波数が欠ける周波数から格子を決める
指向性を無視観測点の音を誤る方向ごとに評価
8乗則を超音速に使う過大評価適用範囲を確認
🙋

関連する内容も知りたいです。

関連シミュレーター

この分野のインタラクティブシミュレーターで理論を体感しよう

シミュレーター一覧

関連する分野

構造解析電磁界解析熱解析
この記事の評価
ご回答ありがとうございます!
参考に
なった
もっと
詳しく
誤りを
報告
参考になった
0
もっと詳しく
0
誤りを報告
0