Skip to content
Published on

工業数学シリーズ 第21回:熱方程式、波動方程式、ラプラス方程式

Authors

工業数学シリーズ 第21回:熱方程式、波動方程式、ラプラス方程式

偏微分方程式の入門で最初に親しくなるべき3つの式があります。この3つは名前だけ覚えるのではなく、どのような現象を説明するのかを先に理解する方がはるかに重要です。

熱方程式

代表形は

ut=kuxxu_t = k u_{xx}

です。

この式は温度や濃度のようなある値が高いところから低いところへ広がり、徐々に均一になる現象をモデリングします。急激な変化が時間とともに滑らかになるのが核心です。

直観

熱い部分と冷たい部分があれば差が徐々に縮まります。熱方程式はまさにこの「平坦化」プロセスを数学的に書いたものです。

波動方程式

代表形は

utt=c2uxxu_{tt} = c^2 u_{xx}

です。

この式は振動が空間に沿って伝播する現象を扱います。弦を一度弾くと形が広がっていく絵をイメージしてください。

直観

熱方程式は鋭い変化が広がって消えますが、波動方程式は形が移動し反射され、振動を維持することができます。

ラプラス方程式

代表形は

uxx+uyy=0u_{xx}+u_{yy}=0

です。

時間が全く入っていないため定常状態の問題を表します。すでに平衡に達した電位分布や温度分布を扱うときによく登場します。

直観

ラプラス方程式は「内側の値が周囲の値と調和的に合わされた状態」を述べています。

3つの式の比較

  • 熱方程式:広がりと緩和
  • 波動方程式:振動と伝播
  • ラプラス方程式:平衡と静的分布

この違いを頭に入れておけば、PDEを初めて見るときも問題の性格をはるかに速く読み取れます。

手で見る短い例題

長さ1の棒の定常状態温度分布を見てみましょう。両端の温度がそれぞれ0と100だとします。定常状態では時間変化がないので

uxx=0u_{xx}=0

です。

したがって

u(x)=Ax+Bu(x)=Ax+B

で、境界条件

u(0)=0,u(1)=100u(0)=0, \quad u(1)=100

を入れると

B=0,A=100B=0, \quad A=100

なので

u(x)=100xu(x)=100x

です。

この例題はPDEが複雑に見えても、状況によっては非常に単純で意味のある結果を与えることを示しています。

工学応用

半導体とチップ冷却

熱方程式はパッケージングと冷却設計の基礎です。

通信と音響

波動方程式は伝送線路、音響、電磁波伝播の基本的直観を提供します。

電場と静電ポテンシャル

ラプラス方程式は電位分布とポテンシャル流れの問題に繰り返し登場します。

よくある間違い

3つの方程式を同じ感覚で見る

外見は似て見えても解の振る舞いはかなり異なります。熱はsmoothing、波動はpropagation、ラプラスはequilibriumです。

境界条件なしで式だけ見る

特にPDEは境界条件が一緒にあってこそ物理的な問題が完成します。

分離解法を早すぎる段階で機械的に適用する

入門段階ではまず物理的意味と代表形態を理解する方がより重要です。

一行まとめ

熱方程式は広がり、波動方程式は伝播、ラプラス方程式は平衡を扱います。

次回予告

次の記事では実数軸を超えて、回路と信号処理で核心的に登場する複素数と解析関数に入ります。

参考資料

  • Erwin Kreyszig, Advanced Engineering Mathematics, 10th Edition
  • Richard Haberman, Applied Partial Differential Equations
  • Lawrence C. Evans, Partial Differential Equations