No.1ベストアンサー
- 回答日時:
困りましたねえ。
ご質問の内容だけではどこに問題があるか判然としません。私の経験をもとに思いつくままに原因と思われることを書いてみますので、ご参考になさってください。1.バグがある
インデクスの範囲指定に間違いがあって保持データを
破壊しているとか、PML領域内のσを与える際の距離の
とり方がおかしいとか。これはもう、asamakenさんの
努力に委ねるしかありません。
2.PMLの理解、離散化に間違いがある
これについては、何度も見直したということですから
大丈夫ですよね。
3.PML領域と通常領域の界面の処理に問題がある
界面における電磁界の値をPML処理と通常処理で交互
に上書きしていませんか?
ちなみに、私はPML(Yeeアルゴリズムです)、MURともに用いたことがありますが、ソルバはきちんと動作しています。うーん、これでは回答とは言い難いですね。申し訳ない。(^.^)
この回答への補足
TCMさん、ご回答ありがとうございました。
質問の内容が漠然としていて申し訳ありません。
いろいろ試してみたところ、垂直に入射する波は吸収できている
のですが、斜めからの入射で発散しているようです。
僕が作っているのは2次元のTM波(Ez、Hx、Hy成分のみ)で、
PML層ではEzx、Ezyに分け、それぞれを計算した後に
Ez=Ezx+Ezyとしています。僕がなにか勘違いをしているかも
しれませんので、もう一度最初から見直してみるつもりです。
プログラムや式などの細かなものは、宇野亨さんの
「FDTD法による電磁界およびアンテナ解析」をかなり参考にしています。
TCMさんの回答の3番についてもう一度見直したところ、Hの計算の際に
通常処理とPML処理とでダブって計算していました。その点を修正してみる
と、見事に吸収することに成功しました。
TCMさん、どうもありがとうございました。またFDTDについての質問
をするかもしれないので、その時は良きアドバイスをお願いします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
原点修正の方法
-
savitzky-golayのフィルタについて
-
時間を100進法であらわしたい。
-
10,000百万円っていくらですか?
-
1アンペアは何ミリアンペアで...
-
合成関数の微分を使う時と、使...
-
「強度」は高い?強い?
-
穴が開く? 空く? 明く?
-
yの二乗をXで微分したら2y・y' ...
-
積分定数Cとは一体なんですか?
-
縞鋼板の曲げ応力度・たわみに...
-
組立式の家具のネジが最後まで...
-
常用対数の求め方 log10の2は約...
-
(4)でなぜcosxsin^4xが偶関数な...
-
角度の計算でdegとradがあると...
-
穴あけのあけってどんな漢字で...
-
100平方キロメートルって東京ド...
-
3割の計算
-
微分可能ならば連続ですが、 不...
-
sin角ABH=AB分のAHって公式で...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報