🔬 音響メタマテリアル
Arrays of subwavelength Helmholtz resonators give a material an effective negative bulk modulus near resonance, opening acoustic band gaps and enabling negative refraction of sound.
仕組み
この1D格子の各ユニットセルには、小さな局所共鳴器 — 固有(角)周波数 ω₀ = √(k/m) を持つヘルムホルツ共鳴器のような質量ばね要素 — が含まれており、本来なら一定の速度で音を伝える背景媒質と弱く結合しています。ω₀から離れた周波数では共鳴器はほとんど応答せず、格子は通常の弾性媒質のように振る舞います。しかしω₀付近では、共鳴器はエネルギーを吸収・再放射し、その際の位相遅れが材料の実効的で周波数依存の体積弾性率K_eff(ω)の符号を反転させることがあります — 光学でおなじみのローレンツ振動子分散が音響学に移植された形です。
K_eff(ω)が負になる範囲では、波動方程式に伝播解は存在しません:波数は虚数になり、入射波は格子内を振動しながら進む代わりに、距離とともに指数関数的に減衰します。この禁制周波数範囲が音響バンドギャップです。共鳴器は、それが制御する音の波長よりもはるかに小さく作ることができるため、この局所共鳴型バンドギャップは、通常であれば(通常のブラッグ散乱型フォノニック結晶のように)何倍も大きな格子が必要になる周波数で開きます — これが音響メタマテリアルの核心的な仕掛けです。
ω₀ = √(k/m) (ユニットセル共鳴)
バンドギャップ: f₀ < f < √(f₀² + fₚ²)
k(ω) = ω · √(ρ₀ / K_eff(ω)) — 複素波数、波 ∝ e^(i(kx−ωt))
よくある質問
音響メタマテリアルは通常の材料と何が違うのですか?
通常の音響材料はその特性を化学(バルク媒質の原子と結合)から得ています。音響メタマテリアルは代わりに、設計された波長以下の構造 — 共鳴するユニットセル(ヘルムホルツ共鳴器など)の周期的配列 — からその特性を得ており、共鳴付近で実効的な負の体積弾性率や負の質量密度を含む、自然の均質材料には不可能な形で音に応答します。
「負」の実効体積弾性率とは、材料が負の体積や負の剛性を持つという意味ですか?
いいえ。これは、共鳴器の固有周波数付近の狭い周波数帯域において、振動する圧力に対する材料の動的で周波数依存の応答が、駆動圧力と位相がずれていることを意味します — 波動方程式で使われる実効弾性率がその帯域で符号を反転するのです。材料自体は完全に通常のものであり、負になるのは駆動された、周波数依存の実効的な記述です。
音響メタマテリアルにおけるバンドギャップとは何ですか?
バンドギャップとは、周期構造内に伝播する波の解が存在しない周波数範囲です。バンドギャップ周波数で材料に入った音は振動波として通過できません。代わりにその振幅は距離とともに指数関数的に減衰し(エバネッセント波)、その構造はその周波数において非常に効果的な、しばしば波長以下の音響バリアとして機能します。
ヘルムホルツ共鳴器とは何で、なぜ使われるのですか?
ヘルムホルツ共鳴器は、狭い首を通じて外部とつながった空洞です。首の中の空気は質量のように、空洞内の空気はばねのように振る舞い、単一の固有共鳴周波数を持ちます。ヘルムホルツ共鳴器は、それが応答する音の波長よりもはるかに小さく作れるため、その配列を使うことで、それが制御する音の波長よりもはるかに小さいユニットセルを持つ局所共鳴型メタマテリアルを構築できます。
音の負の屈折とは何で、どのように生じるのですか?
負の屈折は、メタマテリアル内の音の実効的な位相速度が、実際にエネルギーが流れる方向(群速度)と逆を向くときに生じ、入射ビームは反対側ではなく入射ビームと同じ側の法線方向に曲がります。これは、実効体積弾性率および/または実効密度が同時に負になる共鳴付近で生じ、音に対して実効的な負の屈折率を生み出します。
音響メタマテリアルは実世界でどのような応用に使われていますか?
応用例には、質量だけでは実現できないほど低周波数をよく遮断する薄い局所共鳴型防音パネル、波長以下のイメージングのための音響レンズやスーパーレンズ、物体の周りに音波を迂回させる音響クローキングデバイス、機械類・空調ダクト・交通騒音のためのコンパクトな低周波数騒音バリアなどがあります。
音響メタマテリアルは電磁メタマテリアルとどう関係していますか?
両者は数学的な親戚です。電磁メタマテリアルはスプリットリング共鳴器やワイヤー配列を使って実効的な負の誘電率や透磁率を作り出し、ローレンツ振動子分散で記述されます。音響メタマテリアルはヘルムホルツ共鳴器や局所質量ばね共鳴器を使って実効的な負の体積弾性率や負の質量密度を作り出し、同じローレンツ振動子の数学で記述されます — 異なるのは物理的な場だけです。
局所共鳴型メタマテリアルは、質量をあまり増やさずに音を遮断できますか?
はい、これはその最も有用な特性の一つです。通常の質量則による遮音は、低周波数を遮断するために質量を追加する必要があります。局所共鳴型メタマテリアルパネルは、共鳴器の固有周波数付近で非常に大きな(または負の)実効質量密度を示すことができ、同じ物理的質量に対する質量則の予測よりもはるかに効果的に対象の低周波数帯域を遮断します — 狭い帯域における古典的な質量則の真の回避です。
実際の音響メタマテリアルには実用上どのような制限がありますか?
減衰と材料損失は、理想化された無損失の予測に比べてバンドギャップを狭く浅くします。強い効果は本質的に狭帯域であり、共鳴器の設計周波数を中心としているため、実際の装置では製造公差によってその周波数がずれます。広帯域または複数帯域の性能を実現するには、一般に異なる共鳴周波数を持つ複数の共鳴器設計を積み重ねる必要があります。
このシミュレーションについて
このシミュレーターは、周波数fの平面波でN個の局所共鳴型ユニットセルからなる1次元の鎖を駆動し、その周波数が格子を伝播するか、音響バンドギャップ内でエバネッセントに減衰するかをセルごとに計算します。実効体積弾性率K_eff(ω)は、内蔵されたヘルムホルツ様共鳴器のローレンツ振動子モデルから直接評価され、同じ式がアニメーション化された格子とその下の透過率対周波数曲線の両方を駆動します。
🔬 これが示すもの
複素波数k(f)から瞬時変位が計算されるN個のセルからなるライブ格子で、駆動された波は鎖全体を滑らかに伝播するか、数セル以内で指数関数的に減衰します — さらに、同じK_eff(ω)モデルから計算され、バンドギャップを解析的にトレースする透過率対周波数曲線も表示されます。
🎮 使い方
ユニットセル共鳴f₀、結合強度ωp(バンドギャップ幅)、減衰γ、駆動周波数f、セル数Nを設定すると、fが[f₀, √(f₀²+fₚ²)]に入るのに応じて、ステータス表示が「伝播中」と「バンドギャップ(エバネッセント)」の間で切り替わる様子を観察できます。Pで一時停止、Rでリセットします。
💡 豆知識
局所共鳴型ソニック結晶は、ユニットセル自体の数十倍も大きい波長でバンドギャップを開くことができます — 同じ周波数を遮断するために音波長の半分に匹敵するセルサイズが必要な、通常のブラッグギャップ型フォノニック結晶とは異なります。
よくある質問
表示されるバンドギャップの範囲が、減衰ではなくf₀とfₚのみに依存するのはなぜですか?
bandGapEdges()は[params.f0, Math.sqrt(f0*f0+fp*fp)]、すなわちkeffOf()の無減衰ローレンツ形式のゼロ交差を返します。減衰γはkeffOf()内部の分母の虚部にのみ入り、名目上の端そのものを動かすのではなく、端における遷移を滑らかにし浅くします。
減衰スライダーは、[f_low, f_high]の表示を変えずにどのように図を変えるのですか?
params.gammaを増やすと、keffOf(f)内部の分母の虚部が大きくなり、名目上のギャップの外側の周波数でもwaveNumberOf(f)のIm部分が大きくなります — そのため減衰を高くすると、格子と透過率曲線の「伝播」領域に実際の損失が加わりますが、統計に表示される解析的なバンドギャップの数値は固定されたままです。
ユニットセル数Nを増やすと、なぜ透過率の式自体ではなく減衰エンベロープが変わるのですか?
transmissionAt(f)はMath.exp(-waveNumberOf(f).im * params.N)を計算します。Nを増やすと、同じセルあたりの減衰率がより大きな指数で乗算されるだけなので、長い格子はバンドギャップ周波数をはるかに強く減衰させますが、keffOf()とwaveNumberOf() — セルあたりの物理 — はNによって変わりません。
なぜwaveNumberOf()ではkの虚部が非負に強制されているのですか?
waveNumberOf()はρ₀/K_eff(ω)の複素平方根を取り、ωでスケーリングする前にその虚部にMath.abs()を適用します。これは物理的に因果的な分岐を選択し、エネルギーが格子の駆動される左端から常に離れるように減衰し、誤って増大しないようにします — 受動的で損失のある媒質にとって正しい振る舞いです。
結合強度ωpを大きくすると、なぜ中心を移動させるのではなくバンドギャップを広げるのですか?
params.fpはbandGapEdges()内部でf_high = √(f0² + fp²)という平方根の中にのみ入るため、fpを大きくすると上端が外側に押し出される一方、f_lowはparams.f0に固定されたままです — これは、宿主媒質と局所共鳴器の間の結合が強くなるほど、keffOf()が負のままとなる周波数の範囲が広がるという物理的な描像とまさに一致します。
Arrays of subwavelength Helmholtz resonators give a material an effective negative bulk modulus near resonance, opening acoustic band gaps and enabling negative refraction of sound.
2D · HTML5 Canvas 2D · 60 FPS target · runs fully client-side, no install