前回、Split-TypeのPMLと磁界ベクトルを用いた時間領域FEMの計算を行いましたが、あまり良い結果ではありませんでした。 今回はUniaxial PMLを使った場合の計算結果を記します。 引用文献はこちら。 Dan Jiao, Jian-Ming Jin, Eric Michielssen, and Douglas…

TMzモードの場合の界成分分離型(Split-Type)のPMLでは、x方向PMLのみとすると、 ε dEz1 / dt + σx Ez1 = dHy / dx ε dEz2 / dt = - dHx / dy μ dHx / dt = - dEz / dy μ dHy / dt + σxμε-1 Hy = dEz / dx ここで、 Ez = Ez1 + Ez2 です。 Hx, Hyを消去して、…