境界要素法と体積ペナルティ付き差分法を用いた頭部・耳介伝達関数のモデリング

耳の形が仮想音に与える影響

ヘッドホンを装着して後方や上方から音が来るように感じるとき、脳は頭部や耳の独特な形状が生み出す微細な音響手がかりを利用しています。本論文は、各聴取者を実験室で何時間も測定することなしに、コンピュータ上でそれらの手がかりを高い現実感で再現する方法を探ります。著者らは、音が頭部や外耳の周りでどのように屈曲・反射・回折するかをどれだけ忠実に模擬できるかを調べるため、2つの進んだ数値手法を比較します。

耳が三次元の音を符号化する仕組み

私たち一人ひとりには、頭部伝達関数（HRTF）と呼ばれる個別の「音響指紋」があります。音波が胴体、頭部、そして外耳の複雑なひだに当たると、特定の周波数が強調され、他の周波数が減衰します。これらの変化は方向によって異なります：前方と後方、上と下、左右といった違いです。脳はこれらのパターンを読み取って音源方向を識別します。説得力のある仮想現実や拡張現実のために、音響技術者は各聴取者に合わせ、空間上で非常に密にサンプリングされたHRTFを求めています。実際の人間やダミーヘッドの周囲で直接測定することは可能ですが、時間がかかり技術的に困難で、微小な位置ズレが結果に可聴的な影響を与えることがあります。

同じ聴取問題に対する二つの数学的視点

長時間の測定を回避するため、研究者たちは頭部と耳の詳細な3Dモデルの周りでの音の伝播をシミュレートします。本研究は二つの主要な戦略を比較します。一つは境界要素法と呼ばれる手法で、頭部や耳の表面だけを記述し、その表面が音をどのように散乱するかを解きます。もう一つは差分時間領域法（FDTD）と呼ばれるもので、頭部の周囲の空間を格子で満たし、時間方向に音波を進めていきます。体積法は柔軟性がありますが、大きな領域では計算コストが非常に高くなり得ます。著者らはこれに「体積ペナルティ」という工夫を加えます：耳の表面を格子上で段差状に近似する代わりに、薄い層を介して空気と固体の間を滑らかにブレンドすることで、反射や陰影の表現を大幅に改善します。

単純形状と実際の耳でのモデル検証

人頭全体に適用する前に、チームは制御されたテストケースで各手法を検証します。まず剛体球の周りでの音の回折をシミュレートし、これは教科書的な厳密解が知られています。両手法とも可聴域にわたってこの参照解に良く一致し、境界法と体積法の高密度格子は0.xデシベルの単位で一致します。次に単一の平面壁を研究し、反射音が正しいピークとディップを持つために遷移層で必要な格子点数を決めます。これらのテストから、格子間隔と最小モデル化可能な壁厚を結びつける簡単な規則を導き、可聴上のエラーを導入しないための基準を得ます。その規則を適用して、高解像度で3Dプリントされた耳介をシミュレートし、精密測定と比較します。十分に細かい格子を用いれば、シミュレーションによる耳応答は測定値と平均で約1デシベルほどしか異ならず、これは色付け（音色変化）が聴感上気づかれる閾値に近い値です。

孤立した耳から頭部全体へ

最後に、著者らは3Dオーディオ研究で一般的に用いられる全頭メッシュをシミュレートします。多くの水平方向から到来する音が閉塞された耳道でどのように変換されるかを計算し、ペナルティ付きの体積法と確立された境界法を比較します。耳の最も薄い部分を解像できるほど格子が細かければ、ほとんどの方向と周波数で両手法は非常に良く一致し、知覚される音色変化を予測する聴覚モデルで評価しても同様です。これに対し粗い格子では、耳の共鳴や反射に結びつく重要なピークやディップの周波数や強さがずれてしまい、幾何学的な詳細を犠牲にすると可聴的な影響が生じることが強調されます。

今後の仮想音響への含意

遠方場や大規模領域では、今日のコンピュータ上では境界ベースの手法が依然として効率的ですが、改良された体積法は重要な利点を提供します。小さな内部空洞、空間的に変化する材料特性、さらにはヘッドホンや耳の形状そのものをシミュレーションで最適化する将来の応用を自然に扱えます。本研究は、導き出した指針に従って格子間隔を選べば、体積ペナルティ付き体積シミュレーションが解析解や測定耳データとほぼ可聴閾付近まで一致させうることを示しています。実用的には、これによりすべての耳を実験室で測定することなく、高度に現実的で個人に合わせた3Dサウンドシーンを計算することがより現実的になります。

引用: Hölter, A.B., Lemke, M., Weinzierl, S. et al. Modeling head- and pinna-related transfer functions using boundary elements and finite differences with volume penalization. npj Acoust. 2, 16 (2026). https://doi.org/10.1038/s44384-026-00052-x

キーワード: 頭部伝達関数, 3Dオーディオ, 数値音響学, 外耳のシミュレーション, バーチャルリアリティの音響