-PR-
解決済み

FDTDにおけるPML吸収境界条件について

  • すぐに回答を!
  • 質問No.48459
  • 閲覧数3200
  • ありがとう数2
  • 気になる数0
  • 回答数1
  • コメント数0

お礼率 55% (52/94)

現在FDTDを用いた光導波路解析アプリケーションを作っているのですが、
吸収境界条件に用いたPMLがうまくいきません。
PMLの、解析領域と接している内側の四角から発散してしまいます。論文や参考書
をトレースして、なんども見直したりしたのですがどうもうまくいきません。
なお、PMLはYeeアルゴリズムを用いています。
どなたかご教授願います。
通報する
  • 回答数1
  • 気になる
    質問をブックマークします。
    マイページでまとめて確認できます。

質問者が選んだベストアンサー

  • 回答No.1
レベル10

ベストアンサー率 44% (81/181)

 困りましたねえ。ご質問の内容だけではどこに問題があるか判然としません。私の経験をもとに思いつくままに原因と思われることを書いてみますので、ご参考になさってください。

1.バグがある
 インデクスの範囲指定に間違いがあって保持データを
 破壊しているとか、PML領域内のσを与える際の距離の
 とり方がおかしいとか。これはもう、asamakenさんの
 努力に委ねるしかありません。
2.PMLの理解、離散化に間違いがある
 これについては、何度も見直したということですから
 大丈夫ですよね。
3.PML領域と通常領域の界面の処理に問題がある
 界面における電磁界の値をPML処理と通常処理で交互
 に上書きしていませんか?

 ちなみに、私はPML(Yeeアルゴリズムです)、MURともに用いたことがありますが、ソルバはきちんと動作しています。うーん、これでは回答とは言い難いですね。申し訳ない。(^.^)
補足コメント
asamaken

お礼率 55% (52/94)

TCMさん、ご回答ありがとうございました。
質問の内容が漠然としていて申し訳ありません。

いろいろ試してみたところ、垂直に入射する波は吸収できている
のですが、斜めからの入射で発散しているようです。
僕が作っているのは2次元のTM波(Ez、Hx、Hy成分のみ)で、
PML層ではEzx、Ezyに分け、それぞれを計算した後に
Ez=Ezx+Ezyとしています。僕がなにか勘違いをしているかも
しれませんので、もう一度最初から見直してみるつもりです。

プログラムや式などの細かなものは、宇野亨さんの
「FDTD法による電磁界およびアンテナ解析」をかなり参考にしています。
投稿日時 - 2001-03-08 16:34:54
お礼コメント
asamaken

お礼率 55% (52/94)

 TCMさんの回答の3番についてもう一度見直したところ、Hの計算の際に
通常処理とPML処理とでダブって計算していました。その点を修正してみる
と、見事に吸収することに成功しました。
 TCMさん、どうもありがとうございました。またFDTDについての質問
をするかもしれないので、その時は良きアドバイスをお願いします。
投稿日時 - 2001-03-11 18:28:24
関連するQ&A
-PR-
-PR-
このQ&Aで解決しましたか?
関連するQ&A
-PR-
-PR-
このQ&Aにこう思った!同じようなことあった!感想や体験を書こう
このQ&Aにはまだコメントがありません。
あなたの思ったこと、知っていることをここにコメントしてみましょう。

その他の関連するQ&A、テーマをキーワードで探す

キーワードでQ&A、テーマを検索する
-PR-
-PR-
-PR-

特集


新大学生・新社会人のパソコンの悩みを解決!

いま みんなが気になるQ&A

関連するQ&A

-PR-

ピックアップ

-PR-
ページ先頭へ