波浪水槽シミュレーター

このシミュレーションの根幹をなすのは2D波動方程式で、波の変位の時間変化と空間的な曲率を結びつけます。

$$ \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\) は媒質中の波の伝播速度（定数）です。左辺は波の加速度、右辺は x 方向と y 方向の曲率の和を表しています。

シミュレーターはこの方程式を陽的有限差分時間領域法（FDTD法）で数値的に解きます。安定性のため、Courant–Friedrichs–Lewy（CFL）条件を満たす必要があります。

$$ \frac{c \Delta t}{\Delta x}< \frac{1}{\sqrt{2}}$$

ここで \( \Delta t \) は時間刻み、\( \Delta x \) は格子間隔です。この条件は数値計算が物理的な波の伝播より先走らないように保証し、シミュレーションの発散を防ぎます。係数 \(1/\sqrt{2}\) は2Dシミュレーション特有の値です。

音響室の設計：コンサートホールやオーディトリアムにおいて、音波がどのように反射・干渉するかをモデル化するため、エンジニアは同様の2D波動シミュレーションを使用します。さまざまな壁の形状や素材を仮想的にテストすることで、建設前に最適な音響特性を設計できます。

超音波画像診断：医療用超音波装置は波の反射と干渉の物理に依存しています。高周波音波が異なる組織を伝播する様子をシミュレーションすることで、画像分解能の向上や新しい診断技術の開発に役立ちます。

アンテナ・RF設計：アンテナアレイから放射される電波の伝播も同じ波の原理に従います。CAEシミュレーションは、携帯電話基地局や衛星通信向けに信号強度と指向性を最大化するアンテナ配置設計に不可欠です。

地震解析：地球物理学者は、地震波が地球の各層をどのように伝わるかを理解するために波動伝播モデルを使用します。これは地震リスク評価や石油・ガス探査にとって極めて重要な解析です。

このシミュレーターを使い始めたとき、特にCAE初学者が誤解しがちなポイントがいくつかあります。まず一つ目は 「波速 (c) を現実の水の波に近づけて設定すれば、より正確なシミュレーションになる」という思い込みです。確かに波速は物理的なパラメータですが、数値計算の安定性を決めるCFL条件 ($c \Delta t / \Delta x \lt 1$) と強く結びついています。例えば波速 c を大きくしすぎると、この条件を満たすために時間刻み $\Delta t$ を極端に小さくする必要が出て、計算が重くなったり不安定になったりします。実務では、現象の本質を捉えつつ計算コストを抑えるバランスが重要です。

二つ目は 減衰係数の設定についてです。このツールでは「水の抵抗」を模倣していますが、実際には減衰のモデル化には選択肢があります。例えば振動構造物の解析では、現象に応じて速度に比例する「粘性減衰」や変位に比例する「履歴減衰」などのモデルが使い分けられます。シミュレーターの減衰を強くしすぎると干渉縞が見えにくくなるので、観察目的に合わせて調整するのがコツです。

最後に、「境界条件は単に壁があるだけ」という理解です。このツールでは「固定端」（壁で波が反射）を採用していますが、実際には「自由端」（壁が振動）や「吸収端」（境界で波が消える）など、多様な条件が存在します。例えば無響室の設計では、音をできるだけ吸収するように壁に「吸収端」条件が設定されます。シミュレーション結果を解釈する際は、どんな境界条件を仮定したのかを常に意識することが、実務での落とし穴を防ぐ第一歩です。