設定
統計
0.00
最大振幅
0.00
総エネルギー
2D波動方程式を有限差分法でリアルタイム計算。クリックで波源を生成し、干渉・回折・反射パターンを観察できます。
このシミュレーションの根幹をなすのは、2次元の波動方程式です。これは、時間と空間における媒質の変位（ここでは水面の高さ）の変化を記述します。
$$ \frac{\partial^2 u}{\partial t^2}= c^2 \left( \frac{\partial^2 u}{\partial x^2}+ \frac{\partial^2 u}{\partial y^2}\right) $$ここで、$u(x,y,t)$ は位置$(x,y)$、時刻$t$における水面の変位、$c$は媒質中の波の伝播速度（波速）です。この偏微分方程式をコンピュータで解くために、有限差分法（FTCS法）が用いられています。
数値計算を安定して行うためには、時間刻み$\Delta t$と空間格子間隔$\Delta x, \Delta y$の間にCFL条件と呼ばれる制約が必要です。このシミュレーターでは以下の条件が満たされています。
$$ c \frac{\Delta t}{\Delta x}< \frac{1}{\sqrt{2}} $$この条件を破ると計算が発散し、波が無限に大きくなる非物理的な結果になってしまいます。また、減衰係数（FAQより）を各ステップで乗算することで、現実の水の抵抗に相当する減衰を表現しています。
音響・騒音解析：コンサートホールや劇場の設計では、音波の反射、干渉、回折をシミュレーションし、音響特性を最適化します。自動車や航空機の室内騒音の予測と低減にも同様の技術が用いられます。
構造物の振動解析：地震時の建物の揺れや、風による橋梁の振動は、波動方程式の応用としてモデル化できます。固有振動数や共振現象を事前に把握し、耐震・耐風設計に役立てます。
電磁波解析：アンテナの設計や電子機器内での電磁波干渉（EMI）の予測に利用されます。波動方程式は光（電磁波）の伝播も記述するため、光学デバイス設計の基礎技術でもあります。
海洋・土木工学：津波の伝播シミュレーション、港湾の波浪予測、防波堤による波の反射・散乱の評価などに応用され、防災や構造物設計に貢献しています。
このシミュレーターを使い始めるときに、特にCAE初心者が勘違いしやすいポイントがいくつかあります。まず一つ目は、「波の速さ(c)を現実の水の波そっくりに設定すれば、より正確なシミュレーションになる」という思い込みです。確かに波速は物理的なパラメータですが、数値計算の安定性を決めるCFL条件（$c \Delta t / \Delta x < 1$）と強く結びついています。例えば、波速cを大きくしすぎると、この条件を満たすために時間刻み$\Delta t$を極端に小さくする必要が出て、計算が重くなったり、不安定になったりします。実務では、現象の本質を捉えつつ計算コストを抑えるバランスが重要です。
二つ目は、減衰係数の設定についてです。このツールでは「水の抵抗」を模倣していますが、実はこの減衰の与え方にはモデル化の選択肢があります。例えば、振動が伝わる構造物の解析では、減衰を速度に比例させる「粘性減衰」や、変位に比例させる「履歴減衰」など、現象に応じてモデルを使い分けます。シミュレーターの減衰を強くしすぎると、干渉縞がはっきり見えなくなってしまうので、観察目的に合わせた調整がコツです。
最後に、「境界条件は単に壁があるだけ」という理解です。このツールでは「固定端」（壁で跳ね返る）を採用していますが、実務では「自由端」（壁が振動する）や「吸収端」（波を境界で消す）など多様です。例えば、防音室の設計では、音をできるだけ吸収する「吸収端」の条件が壁に設定されます。シミュレーション結果を解釈する際には、どんな境界条件を仮定したのかを常に意識することが、実務での落とし穴を防ぐ第一歩です。
この2D波動方程式のシミュレーションは、そのままの形、または発展形で、様々な先端工学分野の根幹を支えています。一つは非破壊検査の分野です。例えば、航空機の複合材や橋梁の鉄骨に超音波を送信し、内部のき裂からの反射波を計測・解析することで、欠陥の位置や大きさを特定します。このシミュレーターで壁からの反射を観察するのと物理的な原理は同じで、計測データの解釈にはシミュレーションによる「予習」が欠かせません。
もう一つは医用画像処理、特に超音波診断装置です。プローブから体内に送り込んだ超音波パルスが、臓器の境界で反射したり、組織中を伝わったりする様子は、まさに波動方程式で記述されます。異なる組織間の音響インピーダンス（音の伝わりにくさ）の違いを境界条件として設定することで、より現実に近いシミュレーションが可能です。これにより、新しい画像再構成アルゴリズムの開発や、診断技術の向上に貢献しています。
さらに地盤工学や資源探査にも応用されています。人工地震波を地中に送り、地層の境界で反射して戻ってくる波（地震波探査）を分析することで、地下構造を明らかにしたり、石油・天然ガスの貯留層を探したりします。この場合、媒質は均一ではなく、地層ごとに波速cが異なるため、シミュレーションはより複雑になりますが、基礎となる物理は本ツールで体験している波動の伝播と反射そのものです。
このツールに慣れてきたら、次は「なぜこの計算で波が進むのか」という数値解法の核心に踏み込んでみましょう。まず理解すべきは有限差分法の考え方です。波動方程式の二階微分 $\frac{\partial^2 u}{\partial t^2}$ や $\frac{\partial^2 u}{\partial x^2}$ を、隣接する点の値の差（差分）でどう近似するかを学びます。例えば、時間に関する二階微分は、$$ \frac{\partial^2 u}{\partial t^2} \approx \frac{u^{n+1}_{i,j} - 2u^{n}_{i,j} + u^{n-1}_{i,j}}{(\Delta t)^2} $$ と離散化されます。この式を波動方程式に代入し、未来の値 $u^{n+1}$ について解くことで、次の時間ステップの波の形が計算できる仕組みです。
次のステップとしては、3次元への拡張を考えてみることをお勧めします。現実の音波や電磁波は3次元空間を伝わります。その支配方程式は、$$ \frac{\partial^2 u}{\partial t^2}= c^2 \left( \frac{\partial^2 u}{\partial x^2}+ \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2}\right) $$ となります。計算量は格段に増えますが、本質的な解法は2Dと同じです。これが、自動車の室内騒音（3D空間内の音圧分布）のシミュレーションなどに直接繋がります。
最終的には、異なる数値解法への発展を視野に入れましょう。有限差分法は直感的ですが、複雑な形状（例えば、曲がった防音壁や不規則な地形）を扱うには不向きです。そこで、有限要素法（FEM）や境界要素法（BEM）といった手法を学ぶことになります。これらの手法は、計算領域を小さな要素（メッシュ）に分割して近似する点では共通ですが、波動問題に対してより柔軟で強力なアプローチを提供します。このツールで波動現象の物理的なイメージをしっかり掴んだことが、それらの高度な数値手法を理解するための最良の土台となるでしょう。