時間領域FEMで完全整合層(PML:Perfectly Matched Layer)を用いる場合の定式化にはいくつかあるようです。 いままで用いてきたのは、Jiaoらによる定式化です。彼らの定式化による結果は、時間領域FEMで解く場合もっとも良好な結果が得られています。 今回、下…

Givoli,Neta,Patlashenkoが提出した高次吸収境界条件は、導波管のスカラ波の吸収には効果がありました。 今回は境界上の媒質が変化する場合について試してみました。文献では媒質は一様として導出されているためそのまま使用できません。先ずは媒質が変化す…

前回Givoli-Neta-Patlashenkoの高次ABCを時間領域有限要素法(TD-FEM/ FETD)に用いてH面導波管直角コーナーベンドの散乱係数を計算しました。 前回記事： 時間領域FEMによるH面導波管直角コーナーベンドの散乱係数の計算(Givoli-Neta-Pattashenko's High orde…

前回の記事で試したGivoli-Neta-Patlashenkoの高次吸収境界条件(High order Absorbing Boundary Conditions)を用いた時間領域FEMをDelFEM4Netで実装してみました。 今回は、前回の記事で用いた反復計算は行わないで、電界Ez、及び補助変数Φ1～ΦJ-1に関する方…

時間領域FEMに適用可能な高次吸収境界条件(Higer oder absorbing boundary condition)の構成方法が下記の文献に記されていたので試してみました。 Dan Givoli , Beny Neta & Igor Patlashenko "Finite element analysis of time-dependent semi-infinite wav…

時間領域FEMで完全整合層(PML)を用いる方法は計算精度が良いのですが、畳み込み積分を評価する関係で時間ステップ毎にFEM行列の組み立てを行う必要があり時間がかなりかかりました。 それに比べると、吸収境界条件(ABC)を用いる方法はFEM行列の組み立てはあ…