1.6 モード関数の内積について(TE-TM)

ストークスの定理を用いた計算

 TEモードとTMモードの場合, \begin{eqnarray} I_{[m](n)} &=& \iint_S \VEC{e}_{[m]} \cdot \VEC{e}_{(n)} \ dS \nonumber \\ &=& \iint_S \big( \VEC{a}_z \times \nabla_t \Psi_{[m]} \big) \cdot \big( -\nabla_t \Psi_{(n)} \big) dS \nonumber \\ &=& \iint_S \big( \nabla_t \Psi_{(n)} \times \nabla_t \Psi_{[m]} \big) \cdot \VEC{a}_z dS \end{eqnarray} ここで, \begin{eqnarray} \nabla_t \times \big( \Psi_{(n)} \nabla_t \Psi_{[m]} \big) &=& \nabla_t \Psi_{(n)} \times \nabla_t \Psi_{[m]} + \Psi_{(n)} \nabla _t \times \nabla_t \Psi_{[m]} \nonumber \\ &=& \nabla_t \Psi_{(n)} \times \nabla_t \Psi_{[m]} \label{eq:te-tm-i1} \end{eqnarray} これを面積分して, \begin{eqnarray} &&\iint_S \Big\{ \nabla \times \big( \Psi_{(n)} \nabla_t \Psi_{[m]} \big) \Big\} \cdot \VEC{a}_z dS \nonumber \\ &=& \iint_S \Big\{ \Big( \nabla _t + \frac{\partial}{\partial z} \VEC{a}_z \Big) \times \big( \Psi_{(n)} \nabla_t \Psi_{[m]} \big) \Big\} \cdot \VEC{a}_z dS \nonumber \\ &=& \iint_S \Big\{ \nabla_t \times \big( \Psi_{(n)} \nabla_t \Psi_{[m]} \big) \Big\} \cdot \VEC{a}_z dS \label{eq:te-tm-i2} \end{eqnarray} ストークスの定理 \begin{gather} \iint_S \Big\{ \nabla \times \big( \Psi_{(n)} \nabla_t \Psi_{[m]} \big) \Big\} \cdot \VEC{a}_z dS = \oint_C \Psi_{(n)} \nabla_t \Psi_{[m]} \cdot d\VECi{\sigma} \end{gather} より, \begin{gather} \iint_S \Big\{ \nabla_t \times \big( \Psi_{(n)} \nabla_t \Psi_{[m]} \big) \Big\} \cdot \VEC{a}_z dS = \oint_C \Psi_{(n)} \frac{\partial \Psi_{[m]}}{\partial \sigma} d\sigma \label{eq:te-tm-i3} \end{gather} ただし,面$S$は $\VEC{a}_z$ が法線方向となる平面(導波管断面), $d\VECi{\sigma}$ は周回積分路のベクトル線要素, $+\sigma$ 方向は $\VEC{a}_z$ に対して右ねじの方向である.よって,これらの結果より次式が得られる \begin{gather} I_{[m](n)} = \iint_S \big( \nabla_t \Psi_{(n)} \times \nabla_t \Psi_{[m]} \big) \cdot \VEC{a}_z dS = \oint_C \Psi_{(n)} \frac{\partial \Psi_{[m]}}{\partial \sigma} d\sigma \label{eq:Imn_one} \end{gather} 式\eqref{eq:te-tm-i1}$\sim$式\eqref{eq:te-tm-i3}において, $\Psi_{[m]}$と$\Psi_{(n)}$を交換して同様に求めると, \begin{align} &\iint_S \Big\{ \nabla_t \times \big( \Psi_{[m]} \nabla_t \Psi_{(n)} \big) \Big\} \cdot \VEC{a}_z dS = \oint_C \Psi_{[m]} \frac{\partial \Psi_{(n)}}{\partial \sigma} d\sigma \nonumber \\ &\iint_S \big( \nabla_t \Psi_{[m]} \times \nabla_t \Psi_{(n)} \big) \cdot \VEC{a}_z dS = \oint_C \Psi_{[m]} \frac{\partial \Psi_{(n)}}{\partial \sigma} d\sigma \end{align} よって, \begin{gather} -I_{[m](n)} = \oint_C \Psi_{[m]} \frac{\partial \Psi_{(n)}}{\partial \sigma} d\sigma \end{gather} また, \begin{eqnarray} \iint_S \VEC{h}_{[m]} \cdot \VEC{h}_{(n)} \ dS &=& \iint_S \big( -\nabla_t \Psi_{[m]} \big) \cdot \big( -\VEC{a}_z \times \nabla_t \Psi_{(n)} \big) \cdot dS \nonumber \\ &=& \iint_S \big( \nabla_t \Psi_{[m]} \times \nabla_t \Psi_{(n)} \big) \cdot \VEC{a}_z dS \nonumber \\ &=& -I_{[m](n)} \end{eqnarray} まとめると, \begin{eqnarray} I_{[m](n)} &=& \iint_S \VEC{e}_{[m]} \cdot \VEC{e}_{(n)} \ dS = -\iint_S \VEC{h}_{[m]} \cdot \VEC{h}_{(n)} \ dS \nonumber \\ &=& \oint_C \frac{\partial \Psi_{[m]}}{\partial \sigma} \Psi_{(n)} d\sigma = -\oint_C \Psi_{[m]} \frac{\partial \Psi_{(n)}}{\partial \sigma} d\sigma \nonumber \\ &=& -\iint_S \VEC{h}_{[m]} \cdot \VEC{h}_{(n)} \ dS = I_{(n)[m]} \end{eqnarray} モードの直交性は$I_{mn} =0$のときで,境界条件によって上式はゼロとなる.

ガウスの発散定理を用いた別の導出

 ベクトル公式 $\nabla \cdot (w \VEC{A}) = w \nabla \cdot \VEC{A} + \VEC{A} \cdot \nabla w$ を変形して, \begin{gather} \VEC{A} \cdot \nabla w = \nabla \cdot (w \VEC{A}) - w \nabla \cdot \VEC{A} \end{gather} これより, $w \equiv \Psi_{(n)}$, $\VEC{A} \equiv \VEC{a}_z \times \nabla_t \Psi_{[m]}$ とおくと, \begin{eqnarray} &&\big( \VEC{a}_z \times \nabla_t \Psi_{[m]} \big) \cdot \nabla \Psi_{(n)} \nonumber \\ &=& \nabla \cdot \big( \Psi_{(n)} \VEC{a}_z \times \nabla_t \Psi_{[m]} \big) - \Psi_{(n)} \nabla \cdot \big( \VEC{a}_z \times \nabla_t \Psi_{[m]} \big) \end{eqnarray} 上式右辺の第2項について,ベクトル公式 \begin{gather} \nabla \cdot \big( \VEC{B} \times \VEC{C} \big) = \VEC{C} \cdot \big( \nabla \times \VEC{B} \big) - \VEC{B} \cdot \big( \nabla \times \VEC{C} \big) \end{gather} より, $\VEC{B} \equiv \VEC{a}_z$, $\VEC{C} \equiv \nabla_t \Psi_{[m]}$ とおくと次式が得られる(勾配の回転はゼロ). \begin{gather} \nabla \cdot \big( \VEC{a}_z \times \nabla_t \Psi_{[m]} \big) = \big( \nabla_t \Psi_{[m]} \big) \cdot \big( \nabla \times \VEC{a}_z \big) - \VEC{a}_z \cdot \big( \nabla \times \nabla_t \Psi_{[m]} \big) = 0 \end{gather} よって, \begin{gather} \big( \VEC{a}_z \times \nabla_t \Psi_{[m]} \big) \cdot \nabla_t \Psi_{(n)} = \nabla_t \cdot \big( \Psi_{(n)} \VEC{a}_z \times \nabla_t \Psi_{[m]} \big) \end{gather} 2次元の$\nabla _t$に関するガウスの発散定理より, \begin{eqnarray} \iint _s \big( \VEC{a}_z \times \nabla_t \Psi_{[m]} \big) \cdot \nabla_t \Psi_{(n)} dS &=& \iint _s \nabla_t \cdot \big( \Psi_{(n)} \VEC{a}_z \times \nabla_t \Psi_{[m]} \big) dS \nonumber \\ &=& \oint _C \big( \Psi_{(n)} \VEC{a}_z \times \nabla_t \Psi_{[m]} \big) \cdot \VEC{n} d\sigma \nonumber \\ &=& \oint _C \Psi_{(n)} \big( \VEC{n} \times \VEC{a}_z \big) \cdot \nabla_t \Psi_{[m]} d\sigma \nonumber \\ &=& \oint _C \Psi_{(n)} (-\VEC{a}_\sigma) \cdot \nabla_t \Psi_{[m]} d\sigma \nonumber \\ &=& -\oint _C \Psi_{(n)} \frac{\partial \Psi_{[m]}}{\partial \sigma} d\sigma \end{eqnarray} ただし,$\VEC{n}$ は面$S$上(導波管断面)における周回積分路の外向き法線単位ベクトルを示し, 周回積分路に沿う方向の単位ベクトルを $\VEC{a}_\sigma \equiv \VEC{a}_z \times \VEC{n}$ とおいている.したがって, \begin{gather} I_{[m](n)} = -\iint _s \big( \VEC{a}_z \times \nabla_t \Psi_{[m]} \big) \cdot \nabla_t \Psi_{(n)} dS = \oint_C \Psi_{(n)} \frac{\partial \Psi_{[m]}}{\partial \sigma} d\sigma \end{gather}