1次元、2次元、3次元における波動方程式

直線、平面、空間に生じた一つの波の山。それが叩かれたり弾かれたりして広がる様子を観察しましょう。波が通り過ぎたずっと後に、距離rの観測点が何を記録するのかを読み解いてみてください。

インタラクティブシミュレーションを読み込んでいます...

レッスン

理論 — 1次元、2次元、3次元における波動方程式

波動方程式 utt = c²∇²u は、各点におけるuの2階時間微分(uを変位とすれば加速度)が、ラプラシアンのc²倍になることを示しています。ラプラシアンは、その点のすぐ周囲におけるuの平均値が、点自身のuの値からどれだけ上にあるかに比例します。1次元において、すべての解は速度cで反対方向に進む二つの波形になります(ダランベール、1747年)。3次元の場合、ある領域に限定された波が遠くの点に到達するのは、その点を中心とする半径ctの球が領域を通過している間だけであり、その後その点は静止状態に戻ります。しかし2次元ではそうはなりません。残差としての裾が残り、叩かれた後の減衰は1/tに比例するのみです。

各記号の意味

u
波の値:弦や膜の変位、あるいは空間における過剰空気圧などの量
c
波の伝播速度
a
初期の波の山の半径。このページ上のすべての長さはaを単位とし、時間はa/cを単位とする
g
初期の波の山の形。中心からの距離s(単位はa)において(1 − s²)²となり、s ≥ 1でゼロ
前提
端や壁のない直線、平面、空間。そのため、観測点に反射して戻ってくるものはありません。また、変位が小さいことを仮定しているため、方程式は線形になります。つまり、2倍の強さで叩かれた波の山は、ページ上のすべての数値を2倍にします。

叩かれた弦は元に戻らない 🖖

5aで観測(叩く)を押し、青緑色の弦に注目してください。衝撃により左右に斜面が広がり、その後方の弦は0.5333 v₀a/cだけ持ち上がった平坦な状態になります。観測点の変位は波が通過する間に上昇し、ct = 20aでも80aでもその高さに留まります。5aで観測(弾く)に切り替えると、持ち上げられて放たれた同じ波の山は、半分の高さを持つ二つの波形に分かれて遠ざかります。観測点は0.5000まで上昇し、正確に0へと戻ります。ダランベールは1747年に、弦のあらゆる運動は反対方向に進む二つの波形であると示しました。初期の形fと初速度gを用いて表すと、u = ½[f(x − ct) + f(x + ct)] + (1/2c)∫x−ctx+ctg となります。形fは移動する二つの半分とともに去っていきます。速度gは点からct以内の弦の範囲で積分され、その範囲が広がって初期の領域全体を囲い込むと、積分値は全衝撃である½∫g = 8/15になります。現実の弦には固定端があり、そこでの反射が弦を引き戻します。しかし、この弦には端がないため、何かが戻ってくることは決してありません。

膜では、衝撃の到達は決して終わらない 🖖

5aで観測(叩く)を押し、二つの曲面を見てみましょう。それぞれが外側に向かって輪を運んでおり、どちらの輪もct = 6aまでに観測点を通過します。波の後端がすでに通り過ぎた点線の円の内側では、空間は平坦です。それ以降、空間の観測点の値は完全に0になり、表の最後の列はどの桁もゼロになります。一方、膜はその円の内側でも持ち上がったままです。ct = 20aで膜の観測点は0.00861を示し、これは最大値0.0867の9.93%です。5aで観測、80aまで(叩く)を押すと、値は0.00209になります。4倍の時間を待つことで、値は約4分の1に減少しています。この裾は1/tに比例して減衰し、決してゼロには到達しません。二つの解の公式がその理由を説明してくれます。空間におけるある点の解は、その点を中心とする半径ctの球面上での初期状態の平均から作られます。球が波の山を通り過ぎてしまえば、球面上には何も残りません。膜における解は、半径ctの円の内側にある円盤全体にわたる積分です。そのため、波の山は永遠にその内側に留まり続け、1/√(c²t² − ρ²) の重みがかけられます。アダマールは1923年に出版されたイェール大学での講義録の中で、この残された部分を「残差積分」と呼び、球面波とは異なり円筒波は拡散すると記しています。

空間では、放たれた波の山は山と谷になって到達する 🖖

5aで観測(弾く)を押し、青い波形を見てみましょう。初期の波の山にはゼロを下回る部分は全くありませんでしたが、空間の観測点では0.02862までの上昇に続いて同じ深さの下降が記録されます。その面積は0.01667と−0.01667であり、完全に相殺されます。r > aの場合、解は u = (r − ct) g(ct − r)/(2r) となります。これは初期の形に、波の山の中心が到達したちょうどその瞬間にゼロを横切る直線を掛けたものです。音圧もこの方程式に従うため、静止状態で放たれた圧縮空気の塊は、圧縮とその後に続く等しい希薄化として観測点に届き、最終的に押し出す力は残りません。弦は符号を変えることなく、これと同じ波を伝えます。膜は山を記録したあと、ゼロを下回る裾を描き、1/t²に比例して減衰します。これが、アダマールが2次元の残差について導き出した符号の性質です。

全プロセスの詳細解説

  1. 平坦な世界では、拍手の音はどれくらい残り続けるのか? 4 ステップ

    叩かれた膜において、5a離れた観測点は最大値0.0867 v₀a/cを記録します。裾がその1%に減衰するまで、どれくらいの時間がかかるでしょうか? 5aで観測(叩く)を押してください。ページには最大値と、裾の法則の横にct = 20aでの値が表示されます。

    この問題をツールで開く

    1. 波が通過したずっと後、膜の値は点状の衝撃に対する裾の法則に従います。ct = 20aにおいて、二つの値はすでに3桁目まで一致しています。

    2. 最大値の1パーセントは

    3. これと裾の法則を等しいと置いて、ct/aについて解きます。

    4. 空間では、同じ観測点がct = 6aから静寂になります。直径約10 cmの波、すなわちa = 0.05 mの場合、c = 343 m/sの空気中で拍手をすると、時間の単位a/cは

    解答

    約ct = 192aにおいて。これは空間が静寂になる6aの32倍です。 10 cm幅の拍手の場合、空間での1ミリ秒未満に対して28ミリ秒もかかることになります。平坦な世界では、あらゆる音が後に続く音へとその裾を引きずることになります。裾は1/tで減衰するため、しきい値を2%にすれば待ち時間はおよそ半分になります。

学習の道すじ

熱伝導方程式、波動方程式、ラプラス方程式

参考文献 (2)