パラメータ設定
波の種類
媒質1（入射側）
媒質2（透過側）
⚠ 全反射発生！ θᵢ ≥ θcr
理論メモ
垂直入射：
$R = \dfrac{Z_2 - Z_1}{Z_2 + Z_1}$，$T = \dfrac{2Z_2}{Z_2 + Z_1}$
スネルの法則：
$\dfrac{\sin\theta_t}{c_2}= \dfrac{\sin\theta_i}{c_1}$
臨界角（$c_2 > c_1$ のとき）：
$\theta_{cr}= \arcsin\!\left(\dfrac{c_1}{c_2}\right)$
波の反射・透過シミュレーターとは
🧑🎓
このシミュレーターで「インピーダンス」って何ですか？鋼と水とかで値が全然違うみたいですが。
🎓
ざっくり言うと、波が通りにくさを表す「媒質のキャラ」みたいなものだね。音響インピーダンス $Z$ は密度 $\rho$ と波の速さ $c$ を掛けた $Z = \rho c$ で計算するよ。例えば、上のスライダーで「媒質1」を鋼（Z₁≈45 MRayl）、「媒質2」を水（Z₂≈1.5 MRayl）に設定してみて。インピーダンスの差が大きいと、反射係数がほぼ1になって、ほとんど全てのエネルギーが跳ね返されちゃうんだ。
🧑🎓
え、そうなんですか！反射係数が1ってことは、波が全部反射するんですね。でも、入射角を変えると結果も変わりますか？
🎓
もちろん変わるよ！垂直（入射角0度）じゃなくて斜めから波が当たると、反射の仕方が複雑になる。シミュレーターで「入射角」のスライダーをゆっくり大きくしてみて。反射係数と透過係数のグラフが滑らかに変化するのがわかるかな。実務では、非破壊検査で欠陥の角度を探る時とか、この斜め入射の計算がすごく重要になるんだ。
🧑🎓
「臨界角」って表示が出てきました！これが「全反射」が起こる角度なんですか？どういう時に起こるんですか？
🎓
その通り！波が遅い媒質から速い媒質へ進む時（c₁ < c₂）に起こる現象だよ。例えば、水中（c₁）から鋼中（c₂）に超音波を送る場合だね。入射角を臨界角以上にすると、透過波は界面に沿って進むようになり、全てのエネルギーが反射される「全反射」が起こる。シミュレーターで波速c₂をc₁より大きく設定して、入射角を臨界角の表示値まで引き上げてみよう。透過係数が一気にゼロに落ちるのが確認できるはずだよ。
物理モデルと主要な数式
垂直入射（入射角0度）の時の反射・透過係数は、両媒質のインピーダンスのみで決まります。エネルギーの反射率・透過率は係数の2乗で表されます。
$$R = \frac{Z_2 - Z_1}{Z_2 + Z_1}, \quad T = \frac{2Z_2}{Z_2 + Z_1}$$
$R$: 反射係数（振幅比）, $T$: 透過係数（振幅比）, $Z_1, Z_2$: 媒質1, 2の音響インピーダンス ($Z = \rho c$)
斜め入射の場合、波の進行方向はスネルの法則で決まり、反射・透過係数は入射角と波の種類（P波、SH波）に依存するより複雑な式で計算されます。媒質2の波速が媒質1より大きい場合、全反射が起こる臨界角が存在します。
$$\frac{\sin\theta_t}{c_2}= \frac{\sin\theta_i}{c_1}, \quad \theta_{cr}= \arcsin\!\left(\frac{c_1}{c_2}\right) \quad (c_2 > c_1)$$
$\theta_i$: 入射角, $\theta_t$: 透過角, $\theta_{cr}$: 臨界角, $c_1, c_2$: 媒質1, 2中の波の速度
実世界での応用
超音波非破壊検査：金属部材内部のき裂や空洞を検出するために広く用いられます。検査対象（鋼）とカプラント（水やゲル）のインピーダンス差を考慮し、効率的に内部に超音波を送り込むための入射角度や周波数を設計します。
地震波解析：地下の地層構造を探る際、地層境界で反射・透過する地震波を観測します。異なる地層の密度と弾性波速度（インピーダンス）の違いから反射波の強さが決まり、石油探査や断層調査に活用されています。
音響デバイス設計：超音波センサーや医用超音波プローブでは、圧電素子から被検体（人体組織）へ効率的に音響エネルギーを伝達するため、中間層（整合層）のインピーダンスを最適化します。反射を最小化する設計が重要です。
建築・防音設計：音が壁や窓を透過する際の透過損失を評価します。コンクリートと空気のようにインピーダンスが大きく異なる材料の組み合わせでは、音の大部分が反射され、遮音性能が高まります。
よくある誤解と注意点
このシミュレーターを使い始めるときに、特に気をつけてほしいポイントがいくつかあるよ。まず「インピーダンスが同じなら反射はゼロ」というのは垂直入射の時の話だ。斜め入射では、たとえインピーダンスが同じ媒質でも、波の速度が違えばスネルの法則で屈折が起こるし、P波とSH波で結果が分かれる。例えば、密度と速度の組み合わせが違う二つの媒質で、たまたまインピーダンス値が一致することはあり得る。そんな時に「反射しないはず」と思い込んで斜め入射の結果を無視すると、実務で痛い目を見る。
次に、シミュレーター上の「反射係数」は振幅の比であって、エネルギー（強度）の比ではない点だ。実際に反射で失われるエネルギーの割合を知りたい時は、反射係数を2乗する必要がある。例えば、振幅の反射係数が0.5（50%）なら、エネルギー反射率は0.25（25%）になる。残りの75%のエネルギーは透過と他のモード変換に分配されるんだ。グラフを見て「半分くらい反射してるな」と直感的に考えがちだけど、エネルギー的にはもっと少ない場合があるから注意してね。
最後に、臨界角は「波が遅い媒質から速い媒質へ進む時」だけ出現するという根本条件を忘れないで。鋼（P波速度約5900 m/s）から水（約1500 m/s）へ波が進む場合は、そもそも臨界角は存在しない。逆に、水から鋼への入射では、透過角が90度を超える（全反射が起こる）臨界角が計算される。この条件を間違えると、非破壊検査で最適な探触子角度を決められなくなっちゃうから、まずはシミュレーターでc₁とc₂の大小関係を確認するクセをつけよう。
関連する工学分野
この反射・透過の理論は、CAEの世界を超えて、様々な「波」を扱う工学分野の根底にある共通言語だ。まず挙げるのは「電子工学・電磁波工学」。光ファイバーやアンテナ設計では、電磁波が異なる誘電率を持つ媒体の境界で反射・透過する。このときのインピーダンスは「特性インピーダンス」と呼ばれ、式の形は音響インピーダンスとそっくり（$Z=\sqrt{\mu/\epsilon}$）で、同じく整合が重要になるんだ。
もう一つは「機械工学における振動・騒音解析」。自動車のボディや建築構造物の振動は、弾性波として部材間を伝わる。溶接部やボルト結合部のような「界面」で振動エネルギーがどう反射・透過するかを理解することで、振動の伝わり方を制御し、騒音を低減する設計が可能になる。例えば、防振ゴムは意図的にインピーダンスを変化させて、振動の透過を抑えているんだ。
さらに「地球物理学探査」も深く関連している。地中を伝わる地震波の反射を利用した反射法地震探査では、地層ごとの密度と弾性波速度（つまりインピーダンス）の違いから得られる反射波を解析し、地下構造を3次元で可視化する。石油・天然ガスの貯留層探査や、活断層の調査に不可欠な技術で、ここで学ぶ反射係数の計算がその基礎データを作り出していると言えるよ。
発展的な学習のために
このシミュレーターに慣れたら、次は「モード変換」の概念を学ぶのがおすすめだ。実は、斜め入射で界面にP波（縦波）が当たると、反射波や透過波としてSH波（横波）が発生することがある。これは入射波のエネルギーが別の種類の波に変換される現象で、非破壊検査では欠陥の形状を詳細に判別するのに利用される。シミュレーターでP波を選んだ時の計算式の背後には、このモード変換を記述するより複雑な連立方程式（ゾエプリッツの方程式）が隠れているんだ。
数学的な背景を深めたいなら、波動方程式とその境界条件を勉強しよう。反射・透過係数の式はすべて、「界面で変位と応力が連続する」という物理的な境界条件を波動方程式の解に適用することで導出される。例えば、垂直入射の簡単なケースから自分で導出してみると、式の意味が手に取るようにわかるようになる。数式的には、微分方程式と連立一次方程式の組み合わせだね。
最後に、実務に近い学習ステップとしては、シミュレーターのパラメータを実際の材料データベースに基づいて設定してみること。例えば、アルミニウム（Z≈17 MRayl）とエポキシ樹脂（Z≈3 MRayl）の界面、あるいは人体組織（肝臓 Z≈1.7 MRayl）と骨（Z≈6-8 MRayl）の界面など、具体的な組み合わせで反射率を計算してみよう。そうすることで、教科書的な知識が、実際の設計や解析の現場でどう使われるのか、その感覚を養うことができるはずだ。