AlxisDataBlog-数学・物理学


MathJax = {
  chtml: {
    matchFontHeight: false
  },
  tex: {
    inlineMath: [['$', '$']]
  }
};






電荷は質量ある物質です。したがって突然生み出されることもないですし消滅することもありません。任意の連結領域 $\Omega$ 内部の総電荷量の時間変化は $\Omega$ への電流の出入りとバランスをとっています。以前計算した熱エネルギーとフラックスの関係と同じですね。

\begin{align}
\partial_t\rho=-\mathrm{div}j
\end{align}

これを電荷保存則と呼びます。これを用いて、定常な電流に限定していた前回の結論を非定常な場合にも使えるように拡張してみましょう。

\begin{align}
\mathrm{rot}B=\mu_0j
\end{align}

これの発散を計算すると $\mathrm{div}\,\mathrm{rot}=0$ に注意して $\mu_0\mathrm{div}j=-\mu_0\partial_t\rho=0$ が得られます。この式が非定常時に不具合を起こすので、あらかじめ補正項を入れた

\begin{align}
\mathrm{rot}B=\mu_0j+C
\end{align}

から $C$ を求めていくことにしましょう。先程不具合を起こした式は、補正が入ることにより

\begin{align}
-\mu_0\partial_t\rho+\mathrm{div}C=0
\end{align}

となります。ここでクーロンの法則から得られた $\varepsilon_0\mathrm{div}E=\rho$ を時間微分して $\rho$ を消去すると

\begin{align}
-\varepsilon_0\mu_0\partial_t\mathrm{div}E+\mathrm{div}C=0
\end{align}

すなわち $C-\varepsilon_0\mu_0\partial_tE$ は発散のカーネルでなければなりません。発散のカーネルには回転分の不定性がありますが、実験的に $C=\varepsilon_0\mu_0\partial_tE$ となり

\begin{align}
\mathrm{rot}B=\mu_0j+\varepsilon_0\mu_0\partial_tE
\end{align}

これはアンペール マクスウェルの法則と呼ばれています。今回はここまで。次回はマクスウェル方程式の話をします。(written by J.Nishigata)


Copy Right 2020-現在 ALXIS DATA, Inc. All Right Reserved.

電荷保存則

AlxisDataBlog-数学・物理学


MathJax = {
  chtml: {
    matchFontHeight: false
  },
  tex: {
    inlineMath: [['$', '$']]
  }
};






前回、波動方程式を導出しました。

\begin{align}
\Delta f=\dfrac{1}{c^2}\partial_t^2f
\end{align}

この方程式の解には例えば $f=f_0(\langle k,x\rangle-\omega t),\omega/\|k\|=c$ 等があります。これは平面波と呼ばれる解ですが、等位相面を見てみると

\begin{align}
\langle k,x\rangle-\omega t=\mathrm{const.}\Longrightarrow\langle k/\|k\|,\partial_t x\rangle=c
\end{align}

すなわち $k$ 方向に速さ $c$ で伝搬する様子が理解できます。波動方程式はこのような伝搬する波のイメージをもちます。では早速波動方程式を解いていきましょう。方針は以前扱った熱伝導方程式と同様です。位置と時間の分離解 $f=u(x)T(t)$ を代入してみましょう。

\begin{align}
&\Delta u=\dfrac{a}{c^2}u,\\
&\partial_t^2T=aT
\end{align}

時間の方程式に着目すると、$a < 0$ の場合は単振動、$a=0$ の場合は $T=1,t$、$a > 0$ の場合は $T=e^{\pm\sqrt{a}t}$ という特殊解が得られます。この中で $T=t,e^{\sqrt{a}t}$ は時間発散してしまうため、物理学的には実現しません。


$a$ で分類される解の中で今回取り扱うものは $a < 0$ の場合です。$\omega=\sqrt{-a}$ とおきましょう。


このときの波動方程式の特殊解は位置を固定して様子を観察すると、一定の振幅で単振動をしているだけで安定しているように見えます。この特殊解を定常波、そのときの位置の方程式を定常状態の波動方程式と呼びます。


定常状態の波動方程式はヘルムホルツ方程式です。以前解いたことを覚えていますか。球面ベッセル関数を $j_k$、球面ノイマン関数を $n_k$、球面調和関数を $Y_{kl}$ とおくと

\begin{align}
u=\sum\{c_{kl}j_k(\omega \|x\|/c)+d_{kl}n_k(\omega \|x\|/c)\}Y_{kl}(x/\|x\|)
\end{align}

が一般解でしたね。$k$ は整数なので $j_k$ と $j_{-k-1}$ は同一の球面ベッセル方程式の独立する 2 解になっています。これを考慮して、改めて

\begin{align}
u=\sum\{c_{kl}j_k(\omega \|x\|/c)+d_{kl}j_{-k-1}(\omega \|x\|/c)\}Y_{kl}(x/\|x\|)
\end{align}

と書いておいた方が分り易いでしょう。今回はここまで。次回は定常状態の波動方程式の最も基本的である特殊解「球面波」の話をします。(written by J.Nishigata)


Copy Right 2020-現在 ALXIS DATA, Inc. All Right Reserved.

マクスウェル方程式と波動方程式