連続体力学編 第9章

多次元波動方程式の解

この章では、2次元および3次元の波動方程式について、初期値問題の解を扱う。1次元の波動方程式については、第7章で扱った。線型方程式なので、デルタ関数による畳み込みを考える。

9.13次元波動方程式

3次元波動方程式は である。ただし、3次元のラプラス演算子 である。

この方程式()を、初期値 のもとで解きたい。この節では、この解を与える。また、ソース項 がある場合 についての解も求める。

初期値問題の解:式()

1次元波動方程式()の場合、初期値問題の解は、常微分方程式の解に微分演算子を形式的に代入して求めることができた(7.2節)。3次元波動方程式の場合も同様に、形式的な解は と書ける。これは実際に計算することができ、以下の【9.1-注1】の式()のようになる。

【9.1-注1】3次元波動方程式の初期値問題の解:キルヒホッフの解の公式

3次元の波動方程式 の初期値問題の解は、初期値 のもとで、以下のようになる:(キルヒホッフの解の公式)

証明

(未)

波動が伝わる速度は

式()を見ると、波動が伝わる速度が であることが分かる。これは1次元の波動方程式の場合と同じである。特に、 は、 を中心とする半径 の球上の初期値だけで決まっている。これは、3次元の場合の特徴である。

体積力がある場合の初期値問題の解:式()

ソース項 がある場合の波動方程式 の初期値問題の解は、デュアメルの原理により となる。

シミュレーション

3次元波動方程式の数値計算を行うと、右図のようになる。同図左は初期変位 のみがある場合、同右は初期速度 のみがある場合である。 の値は、球の体積で表されている。が正、が負の値である。

周期的なソース項がある場合、右図のようになる。

9.22次元波動方程式

2次元波動方程式は、式()において としたものである。

初期値問題の解:式()

2次元波動方程式の初期値問題の解は、以下の【9.2-注1】で与えられる。この場合もやはり、波動が伝わる速度は である。

【9.2-注1】2次元波動方程式の初期値問題の解:ポアソンの解の公式

2次元の波動方程式 の初期値問題の解は以下のようになる:(ポアソンの解の公式)

証明

(未)

体積力がある場合の初期値問題の解:式()

ソース項 がある場合の波動方程式 の初期値問題の解は、デュアメルの原理により