質問: 現在FDTDを用いた光導波路解析アプリケーションを作っているのですが、 吸収境界条件に用いたPMLがうまくいきません。 PMLの、解析領域と接している内側の四角から発散してしまいます。論文や参考書 をトレースして、なんども見直したりしたのですがどうもうまくいきません。 なお、PMLはYeeアルゴリズムを用いています。 どなたかご教授願います。
回答: 困りましたねえ。ご質問の内容だけではどこに問題があるか判然としません。私の経験をもとに思いつくままに原因と思われることを書いてみますので、ご参考になさってください。 1.バグがある インデクスの範囲指定に間違いがあって保持データを 破壊しているとか、PML領域内のσを与える際の距離の とり方がおかしいとか。これはもう、asamakenさんの 努力に委ねるしかありません。 2.PMLの理解、離散化に間違いがある これについては、何度も見直したということですから 大丈夫ですよね。 3.PML領域と通常領域の界面の処理に問題がある 界面における電磁界の値をPML処理と通常処理で交互 に上書きしていませんか? ちなみに、私はPML(Yeeアルゴリズムです)、MURともに用いたことがありますが、ソルバはきちんと動作しています。うーん、これでは回答とは言い難いですね。申し訳ない。(^.^)