空力音響解析（FW-H方程式）
理論と物理
概要 — FW-H方程式とは
先生、FW-H方程式って何ですか？ CFDの結果から騒音を計算するって聞いたんですけど、どういう仕組みなんですか？
いい質問だね。FW-H方程式は「Ffowcs Williams-Hawkings方程式」の略で、1969年にJohn Ffowcs WilliamsとDavid Hawkingsが発表した音響アナロジーの手法だ。ざっくり言うと、Lighthillの音響アナロジーを、運動する物体表面に拡張したものなんだ。
なるほど…。具体的にはどう使うんですか？ CFDで全部の音を直接計算するんじゃないんですか？
それが問題でね。音は流体の微小な圧力変動だから、遠方まで直接CFDで解こうとすると、数値的な減衰や反射の影響で精度が出ない。しかもメッシュを観測点まで広げなきゃいけないから、計算コストが爆発する。
そこでFW-Hの出番だ。CFDでは物体近傍の流れ場だけを解いて、遠方場の音圧はFW-H積分方程式でポスト処理的に計算する。この2段階アプローチのおかげで、航空機のジェットノイズ、ファン騒音、自動車のサイドミラーの風切り音といった実務的な問題が、現実的な計算コストで解けるようになったんだ。
つまりCFDの「後処理」として使うんですね！ それなら計算領域をコンパクトにできるのか…。でも、なぜ積分で遠くの音がわかるんですか？
音の本質は波の伝播だからね。物体表面（またはその周囲の仮想面）で圧力と速度の変動がわかっていれば、線形波動方程式のグリーン関数を使って任意の遠方点の音圧を計算できる。これがまさにFW-H積分の原理だ。
Lighthillの音響アナロジー
さっき「Lighthillのアナロジーを拡張」って言いましたよね。まずLighthillのほうから教えてもらえますか？
Lighthillの音響アナロジー（1952年）は、空力音響学の出発点だ。ナビエ・ストークス方程式を変形して、次の波動方程式を導く：
ここで $\rho' = \rho - \rho_0$ は密度変動、$c_0$ は一様場の音速、$T_{ij}$ はLighthillの応力テンソルで次のように定義される：
$T_{ij}$って何を表しているんですか？ 右辺を見ると運動量とか圧力とか粘性応力とか、いろいろ入ってますけど…。
ポイントは、左辺が「静止した一様媒質中の波動方程式」の形をしていること。右辺の$T_{ij}$は等価音源として解釈できる。つまり、「実際の乱流をすべて右辺の音源項に押し込んで、左辺はあたかも静止流体中の音の伝播と見なす」——これがアナロジー（類推）の本質だ。
実際の支配的な項は $\rho u_i u_j$ で、これがジェット騒音の主因。レイノルズ応力テンソルとも呼ばれる。低マッハ数では他の項は小さくなるから、$T_{ij} \approx \rho_0 u_i u_j$ と近似できる。
なるほど、乱流を「音源」に読み替えるから「アナロジー」なんですね。でもこれだと物体表面の効果が入ってなくないですか？
鋭いね。そのとおりで、Lighthillの元の定式化は自由空間中の乱流を対象にしている。回転するファンブレードや飛行する航空機の表面効果を扱うには、境界条件を積分に組み込む必要がある。それがFW-H方程式だ。
FW-H方程式の定式化
いよいよFW-H方程式本体ですね。どんな形なんですか？
FW-H方程式は、一般化関数（Heaviside関数、Diracのデルタ関数）を使ってLighthillの方程式を運動する物体表面 $f=0$ に拡張したものだ。結果として遠方場の音圧変動 $p'(\mathbf{x},t)$ は3つの項の和として表される：
ここで $[\cdot]_{\mathrm{ret}}$ は遅延時間（retarded time）での評価を意味し、$r$ は音源点から観測点までの距離、$M_r$ は音源のマッハ数の観測方向成分、$v_n$ は面の法線速度、$\ell_i = (p\delta_{ij} - \tau_{ij})n_j + \rho u_i(u_n - v_n)$ は荷重ベクトルである。
3つの音源項の物理的意味
うわ、3つも項があるんですね…。それぞれ具体的にどんな音なんですか？
よし、1つずつ具体例で説明しよう。
1. 厚み音源（Thickness Noise）：物体が運動すると、その前方で空気が押しのけられ、後方で引き寄せられる。ヘリコプターのブレードが回転するとき、ブレードの「厚み」が空気を周期的に押す効果だ。回転数×ブレード枚数の整数倍の周波数（BPF: Blade Passing Frequency）に離散的なトーンとして現れる。
2. 荷重音源（Loading Noise）：物体表面に作用する圧力（揚力・抗力の変動）が音源になる。ファンブレードが後流と干渉して揚力が変動するとき、これが音になる。自動車のサイドミラーでは、ミラー表面の圧力変動が「ヒュー」という風切り音として車内に伝わるんだ。
3. 四重極音源（Quadrupole Noise）：物体の外側の空間にある乱流そのものが発する音。ジェットエンジンの排気噴流から出る「ゴー」という広帯域騒音が典型例。体積積分が必要で計算コストが高い。ただし低マッハ数（$M < 0.3$程度）では厚み・荷重音源に比べて小さいので、実務的には省略することが多い。
なるほど！ つまり低速の自動車やファンなら面積分の2項だけでだいたいOKで、ジェットエンジンみたいな高マッハ数の問題になると四重極も効いてくるってことですね。
そのとおり。実際にAnsys FluentやSTAR-CCM+のFW-H実装でも、デフォルトでは厚み＋荷重の面積分のみが有効で、四重極はオプション扱いになっている。
Lighthillの第8乗則と音源スケーリング
Lighthillが導いた重要な結果の一つが、低マッハ数ジェットの音響パワーが噴流速度$U$の8乗に比例するという第8乗則だ：
$$ W_a \propto \frac{\rho_0 U^8 D^2}{c_0^5} $$ここで$D$はジェット直径。これは実験と非常によく一致し、空力音響学の基礎となった。物理的には、四重極音源の放射効率がマッハ数の8乗に依存することを意味する。一方、面上の荷重音源（双極子）は$U^6$、厚み音源（単極子）は$U^4$にスケールする。
FW-H方程式の微分形（一般化関数形式）
FW-H方程式の元の微分形は、Heaviside関数$H(f)$を用いて次のように書ける：
$$ \Box^2 [p' H(f)] = \frac{\partial^2}{\partial x_i \partial x_j}[T_{ij} H(f)] - \frac{\partial}{\partial x_i}[\ell_i \delta(f)] + \frac{\partial}{\partial t}[Q_n \delta(f)] $$ここで $\Box^2 = \frac{1}{c_0^2}\frac{\partial^2}{\partial t^2} - \nabla^2$ はダランベール演算子、$Q_n = \rho_0 v_n + \rho(u_n - v_n)$ は質量流束項（厚み音源）、$\delta(f)$ はDiracのデルタ関数である。この微分方程式を自由空間のグリーン関数で解くと、先述の積分形が得られる。
各音源の指向性パターン
| 音源タイプ | 多重極次数 | 速度依存性 | 指向性 | 代表例 |
|---|---|---|---|---|
| 厚み（Thickness） | 単極子（Monopole） | $\propto U^4$ | 全方位 | ブレード回転、ピストン |
| 荷重（Loading） | 双極子（Dipole） | $\propto U^6$ | $\cos\theta$ パターン | ファン騒音、サイドミラー |
| 四重極（Quadrupole） | 四重極 | $\propto U^8$ | $\cos^2\theta$ パターン | ジェット噴流騒音 |
透過面（Permeable Surface）定式化
「透過面」ってたまに聞くんですけど、物体表面じゃなくて仮想的な面を使うんですか？
これは非常に重要な実務テクニックだ。元のFW-H方程式では積分面を物体表面（固体面）に取るけど、di Francescantonio（1997）が提案した透過面定式化では、物体を囲む仮想的な閉曲面（permeable surface）を積分面に使う。
メリットは2つある。まず、四重極の体積積分を避けられること。透過面を十分大きくとって音源領域を包含すれば、面積分だけで四重極の寄与も捕捉できる。次に、非線形の流れ効果（衝撃波を含むような高マッハ数流れ）も面上のデータとして取り込める。
じゃあ透過面を使えば万事解決？ デメリットはないんですか？
残念ながら、透過面を通過する渦が「虚偽音源（spurious noise）」を生むことがあるんだ。特にジェット流れが透過面の下流端を横切ると、面の内側と外側で計算される音がキャンセルしきれなくなって、非物理的な低周波ノイズが乗る。これをend-cap問題と呼ぶ。面の配置は職人芸的なノウハウが要るんだよ。
FW-H方程式を生んだ2つの天才
John Ffowcs Williams（ウェールズ出身、ケンブリッジ大学教授）は「航空音響学の父」と呼ばれ、コンコルドのジェットノイズ研究で名を馳せた。共同研究者のDavid Hawkingsは当時の大学院生で、1969年のPhilosophical Transactions誌に掲載された論文は、その後50年以上にわたって空力音響解析の基盤となっている。一方、理論の土台を築いたSir James Lighthillは流体力学全般の巨人で、音響アナロジーの1952年の論文は引用数10,000回を超える。面白いのは、Lighthillが1952年に提案した理論が、ジェットエンジン1基の騒音試験に数億円かかっていた時代に「計算で騒音を予測する」という革命的なアイデアだったことだ。
数値解法と実装
時間領域積分法（Farassat 1A）
FW-H方程式を実際にコンピュータで解くにはどうするんですか？ あの積分を直接数値的に…？
実務で最も使われているのがFarassat 1A定式化だ。NASAのFarassatが1980年代に導出した、FW-H方程式の時間領域における解の閉形式表現だよ。
厚み音源の寄与は：
荷重音源の寄与は：
遅延時間（retarded time）の計算が大変そうですね。各面要素からオブザーバーまでの音の到達時間を計算して…。
そこがFW-Hコードの実装でキモになるところだ。各面要素について、遅延時間$\tau$を次式から求める：
この陰的方程式をニュートン法で反復的に解く。一様流中であれば解析的に解けるが、平均流の効果がある場合は修正が必要になる。商用ソフトの多くは一様流を仮定した高速アルゴリズムを採用している。
周波数領域法
時間領域以外のアプローチもあるんですか？
周波数領域でFW-H方程式を解く方法もある。面上の圧力・速度データをFFTで周波数成分に分解し、各周波数ごとにヘルムホルツ方程式のグリーン関数で積分する。
メリットは特定周波数帯だけ計算すればいいときに効率がいいこと。デメリットは大量の非定常データを時間方向に保存する必要があること。実務では回転機械のトーン騒音（BPF成分）の解析に向いているね。
CFDとの連成フロー
CFDとFW-Hの連成って、具体的にはどういうデータのやり取りをするんですか？
典型的なフローはこうだ：
- CFD非定常計算を実行し、FW-H面上の圧力$p$、密度$\rho$、速度$\mathbf{u}$を各時間ステップで記録する
- 十分な統計量が得られたら（通常は数十〜数百回のフロー通過時間）、記録したデータをFW-Hソルバーに渡す
- FW-Hソルバーが各オブザーバー点の音圧時刻歴$p'(t)$を計算する
- 音圧時刻歴をFFTしてスペクトル（SPL vs. 周波数）やOASPL（Overall SPL）を算出する
Ansys FluentやSTAR-CCM+では、CFD計算と同時にFW-H積分を実行する「オンザフライ」モードも使える。メモリ節約になるが、後からオブザーバー位置を変更できないのが欠点だ。
乱流モデルの選択
CFD側の乱流モデルって何を使えばいいんですか？ k-εモデルとかでも大丈夫ですか？
RANS（k-ε, k-ω SST等）は空力音響解析には使えない。これは非常に重要なポイントだ。RANSは時間平均した流れ場を解くので、音源となる瞬時の圧力変動を解像できない。
空力音響に必要なのは：
- LES（Large Eddy Simulation）：大スケールの渦を直接解き、小スケールはSGSモデルで近似。最も信頼性が高いが計算コストも高い
- DES/DDES（Detached Eddy Simulation）：壁近傍はRANS、剥離領域はLESで解くハイブリッド。実務的なバランスが良い
- SAS（Scale-Adaptive Simulation）：Menter-Egorovモデル。URANSベースだがLES的な解像が可能
例えば自動車のサイドミラー風切り音であれば、DDES + FW-Hが実務的な標準になっている。ファン騒音ではLES、ジェット噴流ではLESが望ましいが計算コストの関係でDDESも使われる。
メッシュ・時間刻みの要件
音を計算するとなると、メッシュの要件も普通のCFDより厳しいんですか？
そのとおり。空力音響に必要な解像度の目安：
| パラメータ | 推奨値 | 備考 |
|---|---|---|
| 1波長あたりの格子数 | 15〜25点以上 | 対象最大周波数での音波長で判定 |
| CFL数（音響CFL） | $\leq 1.0$ | $\Delta t \leq \Delta x / c_0$ |
| サンプリングレート | $f_s > 2 f_{\max}$ | ナイキスト条件。通常は$f_s \geq 20 f_{\max}$推奨 |
| FW-H面解像度 | 10〜20点/波長 | 面の法線方向に十分な分解能 |
| データ記録時間 | 10〜50 フロー通過時間 | 低周波まで解像するには長い記録が必要 |
10kHzまでの音を解像しようとすると、波長が3.4cmだから格子幅は1.5mm程度…。かなり細かいですね。
だからこそ、FW-H法の価値がある。CFDのメッシュ解像が必要なのは音源近傍だけで、遠方場はFW-H積分が担当する。もしDirect Noise Computation（DNC、直接数値音響計算）でやろうとしたら、音源から観測点まで全域に微細メッシュが必要になる。それでは計算コストが2〜3桁以上増える。
実践ガイド
解析ワークフロー
空力音響解析を最初からやるとして、全体の流れを教えてください。
標準的なCFD + FW-Hワークフローはこうだ：
Step 1: 定常RANS計算（初期場の作成）
- k-ω SSTやSST等で定常解を求める。これが非定常計算の初期条件になる
- この段階でメッシュ品質、y+、収束性を確認
Step 2: 非定常LES/DES計算（音源データ生成）
- 定常解を初期値にLES/DDESに切り替える
- 「遷移時間」（数フロー通過時間）を捨てる。定常解の影響が消えるまで待つ
- 遷移後にFW-H面のデータ記録を開始
Step 3: FW-Hポスト処理（遠方場音圧の計算）
- 記録データからFW-H積分を実行
- オブザーバー点での音圧時刻歴を取得
Step 4: 音響後処理
- FFT → SPLスペクトル、1/3オクターブバンド、OASPL、A特性補正
- 実験データとの比較検証
FW-H面の配置戦略
透過面をどこに置くかって、結構重要なんですよね？ 何かガイドラインはありますか？
FW-H面の配置は結果の精度に直結するから、慎重にやる必要がある。主なルール：
- 音源を完全に包含する：面の内側に全ての音源領域（渦放出域、剥離域、ジェット混合域）が収まること
- 非線形効果が十分減衰する位置に置く：強い渦や衝撃波が面を横切ると虚偽音源が発生
- CFDメッシュが十分に細かい領域内に置く：メッシュ粗化領域に面があると数値散逸で音源データが劣化
- 面の法線方向に3〜5セル分のバッファを確保：スポンジ層や出口境界に近すぎない
- 複数の面でロバスト性を確認：面を少し動かして結果が大きく変わらないことを検証
面を1つ動かすだけで結果が変わるなら、怖いですね…。実務ではどうやって信頼性を担保してるんですか？
経験的なベストプラクティスとしては、少なくとも3つの異なるサイズの透過面で計算して、結果が±1〜2 dB以内に収束していることを確認する。また、固体面でのFW-H結果と透過面でのFW-H結果を比較して、面積分のみで四重極の寄与が捕捉できているかをチェックするのも定番だ。
音響後処理とスペクトル解析
FW-Hで音圧の時刻歴が得られた後、どうやってSPLに変換するんですか？
音圧時刻歴$p'(t)$からのスペクトル解析手順：
- 窓関数の適用：Hanning窓やHamming窓を適用してスペクトルリーケージを低減
- FFT：離散フーリエ変換で$P(f)$を得る
- PSD（パワースペクトル密度）の計算：$S_{pp}(f) = |P(f)|^2 / (\Delta f \cdot T)$
- SPLへの変換：
ここで$p_{\mathrm{ref}} = 20 \; \mu\text{Pa}$は空気中の基準音圧。総合音圧レベル（OASPL）は全帯域の積分で得る：
A特性補正って何ですか？ 規制値で出てきますけど…。
人間の耳の感度は周波数によって違う。低周波（100 Hz以下）と高周波（10 kHz以上）は聞こえにくく、1〜4 kHzが最も敏感だ。A特性（A-weighting）は、この人間の聴覚特性を模擬した周波数補正フィルターで、騒音規制では「dBA」が使われる。FW-Hの計算結果をdBAに変換して評価するのが一般的だよ。
よくある失敗と対策
初心者がハマりやすい落とし穴ってありますか？
| 症状 | 原因 | 対策 |
|---|---|---|
| SPLが実験より20 dB以上高い | RANS計算を使っている（人工的な変動） | LES/DDESに切り替える |
| 低周波に非物理的なピーク | 透過面のend-cap問題 | 面形状の見直し、固体面結果と比較 |
| 高周波が実験より減衰 | メッシュ解像度不足、数値散逸 | FW-H面近傍のメッシュ細分化 |
| データ記録時間が短い | 低周波成分が解像できない | 記録時間を延長（$\Delta f = 1/T$） |
| 結果がFW-H面位置に強く依存 | 面がメッシュ粗化域や渦通過域にある | 面をCFDの高解像度領域内に移動 |
| 双極子パターンが崩れている | 面上の圧力サンプリング不足 | 面の空間分解能を上げる |
ドローンの騒音問題 — FW-Hが使われる最前線
都市部でのドローン配送が現実になりつつある今、最大の壁の一つが「騒音」だ。小型マルチロータの回転翼は数千RPMで回転し、ブレード通過周波数（BPF）とその高調波が耳障りな純音成分として放射される。各国規制当局が定める85 dBA以下という目標値に対して、何枚翼にするか・ピッチ角をどう設定するか・翼形状をどう最適化するかをCFD+FW-Hで仮想評価している。NASAのARM（Advanced Air Mobility）プロジェクトでは、OpenFOAMのLES + FW-H functionObjectを使ったeVTOLの騒音予測が標準的なワークフローとして確立されつつある。
初心者が陥りやすい落とし穴：RANSで空力音響
「k-ε で非定常計算を回して、圧力変動を記録すれば騒音が計算できるのでは？」——これは空力音響解析で最も多い誤りだ。RANS（Reynolds-Averaged Navier-Stokes）は乱流の時間平均を解いているだけで、音源となる渦スケールの変動を解像していない。RANSの「非定常」計算で出てくる圧力変動は数値的なアーティファクトであり、物理的な音源ではない。たとえFW-Hに通しても出てくるSPLは全く意味がない。空力音響にはLES/DES/DDESが必須——これは鉄則だ。
ソフトウェア比較
Ansys Fluent
FluentのFW-H機能について教えてください。設定は難しいですか？
Fluentは商用CFDの中でも最も成熟したFW-H実装を持っている。主な特徴：
- Acoustics Modelメニューから「Ffowcs Williams and Hawkings」を選択するだけ
- 固体面（Solid）と透過面（Permeable）の両方に対応
- オブザーバー点を複数設定可能（円弧配置、グリッド配置）
- 計算と同時にFW-H積分を実行する「オンザフライ」モード
- 面データを書き出して後処理する「ファイルベース」モードも選択可能
設定で注意すべきは、Source Zoneの選択とExport Intervalだ。Export Intervalが粗いとナイキスト周波数が下がり、高周波が折り返してしまう。
Simcenter STAR-CCM+
STAR-CCM+はどうですか？ 回転機械が得意って聞きましたけど。
STAR-CCM+のFW-H実装は回転機械との統合が特に優れている：
- Sliding Mesh/Overset Meshとシームレスに連携し、回転するブレード面を自動的にFW-H音源面として扱える
- 「FW-H impermeable」（固体面）と「FW-H permeable」（透過面）を明示的に分けて設定
- ポスト処理でSPLスペクトル、指向性パターン、A特性補正が一画面で完結
- ファン騒音のトーン＋広帯域の分離も容易
GUIベースの設定が充実しているので、FW-H解析の学習コストは比較的低い。航空宇宙や自動車OEMでの採用実績も多いね。
OpenFOAM
OpenFOAMにもFW-Hがあるんですか？ 無料で使えるのはいいですけど…。
OpenFOAMにはnoiseModel functionObjectとしてFW-H実装がある。ESI版（.com）とFoundation版（.org）で実装が若干異なるが、基本的にはFarassat 1A定式化を採用している。
設定例（controlDictに追記）：
functions
{
noise
{
type noiseModel;
libs ("libfieldFunctionObjects.so");
noiseModel FfowcsWilliamsHawkings;
patchName "fwhSurface";
observers
{
obs1 (0 10 0); // x, y, z [m]
obs2 (10 0 0);
}
rhoRef 1.225;
cRef 340;
writeInterval 1;
}
}注意点として、OpenFOAMのFW-H実装は商用ソフトに比べて検証事例が少ないので、ベンチマーク問題（NACA0012翼型のBPM騒音、SUNYジェットノイズ等）で自分の設定の妥当性を確認してから実務に入ることを強く勧める。
Actran / Virtual.Lab
専用の音響ソフトもあるんですか？
FFT（Free Field Technologies）社のActranとSiemensのVirtual.Lab Acoustics（現Simcenter 3D Acoustics）は、空力音響の専門ツールだ。CFDソルバーとは別に、FW-Hよりも高度な手法（BEM：境界要素法、FEM-BEM連成）を提供する。
これらは特に散乱効果（音波が胴体やナセルに当たって反射・回折する効果）を高精度に扱いたい場合に使われる。例えばエンジンファンの騒音がインレットやエキゾーストダクトを通って放射されるケースでは、FW-H単体では不十分で、BEMやFEMとの連成が必要になる。
| ソフトウェア | FW-H実装 | 透過面 | 回転機械 | 散乱効果 | コスト |
|---|---|---|---|---|---|
| Ansys Fluent | Farassat 1A | 対応 | 対応 | 非対応 | 商用 |
| STAR-CCM+ | Farassat 1A | 対応 | 優秀 | 非対応 | 商用 |
| OpenFOAM | Farassat 1A | 一部対応 | 要設定 | 非対応 | 無償（GPL） |
| Actran | BEM/FEM連成 | N/A | 対応 | 対応 | 商用（高額） |
| Simcenter 3D Acoustics | BEM/FEM | N/A | 対応 | 対応 | 商用 |
「ソルバーより上流のCFDが大事」という現場の声
空力音響の分野で長年仕事をしているエンジニアに聞くと、異口同音に「FW-Hソルバーの差よりも、上流のCFD計算の品質が結果の9割を決める」と言う。同じFW-Hコードを使っても、メッシュ解像度、乱流モデルの選択、時間刻みの設定、記録時間の長さで結果が10 dB以上変わることは珍しくない。逆に言えば、CFDの設定が適切であれば、Fluent、STAR-CCM+、OpenFOAMのどれを使っても、同じ結果が出るべきなのだ。ツール選びで迷う時間があったら、メッシュ収束性の検証に時間を使おう。
先端技術
広帯域騒音モデル
LES+FW-Hは計算コストが高いですよね。もっと手軽に騒音を見積もる方法はないんですか？
RANS結果から騒音レベルを推定する広帯域騒音モデルがある。これはFW-Hとは異なるアプローチで、乱流統計量（k, ε, ω等）から半経験的な式でSPLを見積もる。
- Proudman法：等方性乱流の仮定のもと、$T_{ij}$の自己相関からSPLを推定。Fluentに実装済み
- BPMモデル（Brooks-Pope-Marcolini）：翼型の後縁騒音を境界層パラメータから予測。風車の騒音設計で標準
- Amiet理論：乱流の入射による前縁騒音。ファンのローター・ステータ干渉騒音に適用
精度はLES+FW-Hに及ばないが、設計の初期段階でのスクリーニングやパラメトリックスタディには十分使える。
機械学習との融合
最近は機械学習も使われているんですか？
空力音響でもML（機械学習）の活用が急速に進んでいる：
- サロゲートモデル：LES+FW-Hの結果をニューラルネットワークで学習し、形状パラメータからSPLを瞬時に予測。設計最適化のループに組み込める
- PINN（Physics-Informed Neural Network）：波動方程式の物理則をロス関数に組み込んだニューラルネット。データが少ない場合の汎化に有効
- 超解像（Super-Resolution）：粗いLESの結果から細かいスケールの乱流場を推定し、FW-Hに渡す。計算コストを大幅に削減
ただし現時点ではまだ研究段階のものが多く、規制認証が必要な航空分野では「LES+FW-Hの結果でMLモデルを検証する」という使い方が現実的だ。
UAM・eVTOLへの適用
空飛ぶクルマ（eVTOL）の騒音問題にもFW-Hが使われるんですか？
UAM（Urban Air Mobility）はまさにFW-H解析のキラーアプリケーションだ。eVTOLは複数のロータを住宅地の上空で運用するから、騒音問題は型式認証の最大のハードルの一つになっている。
技術的な課題：
- 複数ロータの干渉：各ロータのウェイクが隣のロータに入射してBVI（Blade-Vortex Interaction）騒音を生む。全ロータを同時にLESで解く必要があり、計算規模は数億セルに達する
- 遷移飛行：ホバーからクルーズへの遷移中は流れ場が刻々と変わるため、定常的なFW-H面設定が困難
- 地面反射：離着陸時に地面からの反射が住民の暴露騒音レベルに大きく影響
NASAは「革命的垂直揚力技術（RVLT）」プロジェクトで、LES+FW-HによるeVTOL騒音予測の大規模ベンチマークを実施し、その結果を公開している。
トラブルシューティング
虚偽音源と数値ノイズ
FW-Hの結果に非物理的なノイズが乗ることがあるって聞きました。どう対処すればいいですか？
虚偽音源（spurious noise）は透過面定式化で最も頻繁に発生する。原因と対策を整理しよう：
- 渦が面を横切る：面の内側で発生した渦が下流端を通過するとき、面の内外での音源キャンセルが不完全になる。→ 面の下流端をジェット・ウェイク領域から離すか、end-cap補正を適用
- メッシュの粗い領域に面がある：数値散逸で面上のデータが劣化し、非物理的な変動が混入。→ 面をメッシュが十分細かい領域に移動
- 時間刻みが粗い：エイリアシング（折り返し）による虚偽の高周波成分。→ サンプリング定理を満たす時間刻みに設定
実験とのSPL不一致
シミュレーション結果と実験のSPLが全然合わないんですが…。何を見直せばいいですか？
SPL不一致のデバッグフロー：
- まず上流CFDを疑う：壁面y+は適切か？ LES/DDESの解像度は十分か？ 時間平均のCp, Cd は実験と合っているか？ → CFDの空力特性が合っていないのに音だけ合うことはない
- スペクトルの形状を比較：レベルのオフセットだけなら音源強度の問題。特定周波数でピークがずれるなら流れ場の特性（剥離点位置、渦放出周波数）の問題
- 周波数依存性をチェック：低周波が合わないなら記録時間が短い。高周波が落ちるならメッシュ解像度不足
- マイク位置・方向を確認：指向性がある音源では、観測角度の1°のずれで数dBの差が出る
- 実験条件との整合：風洞のバックグラウンドノイズ、反射面の有無、マイクの周波数応答を確認
End-cap問題
End-cap問題って先ほど出てきましたけど、もう少し詳しく教えてください。
透過面FW-Hでは、面は閉じた曲面でなければならない。しかしジェットの下流方向に面を閉じると、ジェット流がその「蓋」（end-cap）を横切り、虚偽の低周波ノイズが発生する。
対処法としては：
- 面を閉じない：下流を開いたままにするが、理論的には面が閉じていないので四重極の寄与を完全には捕捉できない。低マッハ数では許容される場合が多い
- 時間平均のend-cap補正：面のend-cap部分での時間平均流束を計算して引き去る手法。FluentとSTAR-CCM+にオプションとして実装されている
- 面をジェットの外に配置：ジェット混合域を完全に包含せず、ジェット音源は固体面FW-Hの荷重項で捕捉する妥協案
実務では「複数の面位置で計算して、end-capの影響が小さいことを確認する」のが最もロバストだ。
デバッグの心構え — 空力音響は「精度の階段」
空力音響のCFD+FW-H解析で±3 dB以内の精度を出すのは、実はかなり高い目標だ。流れ場のCFDで平均量（Cd, Cp）が±5%以内に収まっていても、音圧の変動成分はさらに微細なスケールに依存するため、追加の解像度と注意が必要になる。「CFDの定常解が合ったから空力音響も合うだろう」と思い込むのは危険だ。段階的に品質を上げていくことが重要：まず空力特性 → 次に壁面圧力変動のPSD → そしてFW-HのSPL、という順番で検証しよう。
なった
詳しく
報告