現在FDTD法によるフォトニック結晶光波回路の解析を行っているのですが、入射が思うようにいきません。
僕が行っている解析はEz、Hx、Hy成分で構成されるTM波で、入射のさせかたは、解析領域におけるx=0の位置でEz成分のみについて、
Ez(x,y) = sin((y-y0)π/W) * sin(2πft)
y0:波の中心
W :波の幅
f :入射波の周波数
として、時間的にも位置的にも正弦波の連続波として入射させています。Hx,Hyについてはなにも行っていません。
この入射の計算は、FDTDの電磁界を計算する前に入れており、入射の計算を行ったあとに電界と磁界をそれぞれ計算するようにしています。
しかしこのままでは発散していまってので、FDTDの電界の計算は入射位置を除いた(入射位置がEz(1,j)とするとEz(2,j)から)位置から行っています。
一応このやり方で導波路に波は伝搬するのですが、入射初期の波頭は減衰してなくなってしまいます。しかも遠くに伝搬するにしたがって波はどんどん減衰していき、定常応答を得るまでにかなりの時間を要します。波の振る舞いについてはこれで正しいのでしょうか?
あと、時間と位置的にガウス波の孤立パルスも入射させてみたのですが、すぐに減衰して広がってしまって最初の振幅を保ったまま伝搬しません。なにかおかしい所があるのでしょうか?それともこれで正しいのでしょうか?FDTDに関する論文や本などには、入射についてあまり詳しく書いてないので確かめようがありません。どなたかアドバイスしてくれれば幸いです。
質問にまとまりがなく、わかりにくくて申し訳ありませんが、なにかちょっとしたことでもかまいませんので、ご回答をよろしくお願いします。
No.2ベストアンサー
- 回答日時:
いくつか気になるところがありますので、以下に挙げてみます。
よろしければ参考になさってください。1.励振方法に問題あり
これは、今回の不具合には関係ないかもしれませんが、励振方法に問題があります。
Ez(x,y) = sin((y-y0)π/W) * sin(2πft)
とのことですが、これでは励振面の電界に他の部分から伝播してきた電磁界の影響を盛り込むことができません。ちなみに、このような励振をhard sourceといいますが、やはり、
Ez(x,y)=Ez(x,y)+A*sin(ωt)-A*sin(ω(t-Δt))
とすべきだと思います。また、私の場合、電磁界の更新アルゴリズムは以下のようにしています。
電界の更新
電界の励振
PMLによる電界の更新
磁界の更新
PMLによる磁界の更新
なお、上式で-A*sin(ω(t-Δt))とした意味については、よーく考えてみてください。磁界の励振も行うとこの項は必要なくなるんですけどね。
2.励振面の位置がx=0である
x=0の位置はPML吸収境界面になっていると思いますが、私の場合は吸収境界面と励振面を数セル離しています。まあ、問題ないような気もしますが。
3.光導波路には遮断周波数のようなものがある?
光導波路のことはよくは知らないのですが、導波管の遮断周波数fcのように、それ以下の周波数の電磁波を伝播させないというような特性はないのでしょうか? この部分、自信なしです。
参考文献:Allen Taflove,"COMPUTATIONAL ELECTRODYNAMICS The Finite-Difference Time-Domain Method",ARTECH HOUSE,1995.
TCMさん、貴重なご助言ありがとうございます。
さっそくTCMさんのご回答をもとに、プログラムを見直してみます。
TCMさんには何度もお世話になっていますね。本当にありがとうございます。
自分でもいろいろ文献を調べて勉強しているつもりなのですが、まだまだのようです。もっと勉強します!
どうもありがとうございました。
No.1
- 回答日時:
素人がこんな事を言って申し訳ありませんが,どの様な状況での解析でしょうか。
御質問文から受ける印象では,大学か企業での研究の一環のように思えるのですが。そうだとすると,ここで御質問されるよりは,周りに居るであろう先輩なり指導者なりに相談されるのが先だと思います。
この回答への補足
ご回答ありがとうございます。
僕は大学生で、FDTDによる解析は卒業研究のテーマにしています。
rei00さんのおっしゃる通り、先生などに相談すべきなのですが、FDTDのことを知っていても使ったことがなく、詳しいところまで説明できる先生がいません。
そこで神にもすがる思いで質問した次第です。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 高校物理です。 (1)がよく分かりません。 Oから出た波が壁で反射する時、位相は変わらない。また、入 1 2023/07/28 13:57
- 都市伝説・怖い話 イルミナティ・ディープステイト・ㇵザールマフィアの回し者である関暁夫君がやりすぎ都市伝説 1 2022/12/23 22:26
- その他(自然科学) 電磁波の周波数と熱について教えて下さい。 電磁波の波長とエネルギーについて、雑学として興味があります 6 2022/04/18 20:00
- 物理学 有限の大きさの物質では、周期的境界条件を満たすように格子振動が発生する。もし、満たさない場合、物質の 1 2022/07/05 18:37
- 物理学 電子レンジが物を加熱する仕組みについて 吸収できるエネルギーに限界値はあるのか知りたい。 表題の通り 4 2022/05/06 20:10
- 工学 高周波工学について 3 2022/05/11 11:04
- その他(暮らし・生活・行事) 緊急!電磁波を受けて耳鳴りで悩んでますが 最適な商品を紹介してもらえませんか? 私はどこからか電磁波 1 2022/05/06 21:34
- その他(自然科学) この2つはどこが間違いでしょうか? 2. 植物の葉に含まれるクロロフィルは可視赤色域の波長帯を強く吸 2 2022/05/17 16:55
- 物理学 物理 7 2023/08/05 11:51
- 物理学 開管共鳴と定常波 2 2023/03/23 11:15
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報