Processing math: 0%

2014年8月11日月曜日

マクスウェルの方程式から波動方程式の導出

マクスウェルの方程式から波動方程式の導出をする。 まず条件として真電荷と伝導電流ともにない空間で、\bm{E}\bm Hztのみの関数であるとする。

ガウスの法則 \rm{div} \ \bm D = 0

は本来

\frac{\partial D_x}{\partial x} + \frac{\partial D_y}{\partial y} + \frac{\partial D_z}{\partial z} = 0 

となるが、ztのみの関数という条件のもと始めたのでz微分の項以外は消える。よってD=\epsilon Eより

\epsilon \frac{\partial E_z}{\partial z} = 0 \tag{1}

となる。同様に磁界に関するガウスの法則

\rm{div} \ \bm B = 0

についても同様に

\mu \frac{\partial H_z}{\partial z} = 0 \tag{2}

次にファラデーの法則

\rm rot \ \bm E  = - \frac{\partial \bm B}{\partial t}

\op{E_x}{E_y}{E_z}{x}{y}{z} = \left( -\rd{B_x}{t}, \ -\rd{B_y}{t}, \ -\rd{B_z}{t} \right)

だが例によってx,yによる微分の項は消えるので

\left( -\rd{E_y}{z}, \ \rd{E_x}{z}, \ 0 \right) = \left( -\rd{B_x}{t}, \ -\rd{B_y}{t}, \ -\rd{B_z}{t} \right) \tag{3}

よって

\mu \rd{H_z}{t} = 0 \tag{4}

拡張されたアンペールの法則

\rm rot \ \bm H = \bm j + \rd{ \bm D}{t}

\op{H_x}{H_y}{H_z}{x}{y}{z} = \left(j_x,\  j_y,\  j_z \right) + \left( \rd{D_x}{t}, \ \rd{D_y}{t}, \ \rd{D_z}{t} \right)

だが、先程と同様にz微分以外を消し、また伝導電流\ \bm jがないという条件なので消え

 \left( -\rd{H_y}{z}, \ \rd{H_x}{z}, \ 0 \right) = \left( \rd{D_x}{t}, \ \rd{D_y}{t}, \ \rd{D_z}{t} \right)  \tag{5}

となり

\epsilon \rd{E_z}{t} = 0 \tag{6}

続いて(3)より

\begin{align} \rd{E_x}{z} &= -\mu \rd{H_y}{t} \tag{7} \\ \rd{E_y}{z} &= \mu \rd{H_x}{t} \tag{8} \end{align}

続いて(5)より

\begin{align} \rd{H_x}{z} &= \epsilon \rd{E_y}{t} \tag{9} \\ \rd{H_y}{z} &= -\epsilon \rd{E_x}{t}  \tag{10}\end{align}

まとめると(1)(2)よりE_zH_zは定数であることがわかる。そして電磁波の場合ともに0とおくことができる(らしい。よくわからない。) 

後はうまいこと式を変形させて波動方程式の形に持っていけば良い。例えば(7)を両辺zで微分して(10)も両辺tで微分する。積分順序は積分結果に関係ないので順番が違っても同じものとみなせば

\frac{\partial^2E_x}{\partial z^2} = \epsilon \mu \frac{\partial^2 E_x}{\partial t^2}

同様にうまいことやれば

\frac{\partial^2E_y}{\partial z^2} = \epsilon \mu \frac{\partial^2 E_y}{\partial t^2}

\frac{\partial^2H_x}{\partial z^2} = \epsilon \mu \frac{\partial^2 H_x}{\partial t^2}

\frac{\partial^2H_y}{\partial z^2} = \epsilon \mu \frac{\partial^2 H_y}{\partial t^2}

0 件のコメント:

コメントを投稿