以前ここで、FDTD法における入射について質問(No.67248)したのですが、そのときソフトソースの励振方法として、電界の入射を
Ez(x,y)=Ez(x,y)+A*sin(ωt)-A*sin(ω(t-Δt))
と行えばよいというご指摘をTCMさんから得ました。
しかし、実際に上式のような入射方法をとったところ、電界の振幅が非常に小さく、伝播する振幅がAになりません。
また、電磁界の更新アルゴリズムも
電界の更新
電界の励振
PMLによる電界の更新
磁界の更新
PMLによる磁界の更新
のようにTCMさんと同じようにとったのですが、うまくいきません。
なにか解決方法があればご教授願います。お願いします。
A 回答 (2件)
- 最新から表示
- 回答順に表示
No.2
- 回答日時:
物理的には正しい計算になっているでしょうか?
以前の質問では単純に減衰しているとお書きになっていますが、
減衰には吸収による減衰と導波路からの放射による減衰があると思います。
・(吸収媒質で無い場合)エネルギーは全体で保たれているでしょうか(境界から出ていく量も含めて)?
・バンドギャップは十分開いているでしょうか?
・これは私も興味あるのですがPBGの導波路って
カットオフとか持っていないのでしょうか
(通常の誘電体光導波路では対称な場合はカットオフないのですが、
同じに考えていいのでしょうか?)...と書いていたら、
ご自分でも疑問に思ってらっしゃるようですね。
誘電体多層膜でサンドイッチした構造のときのフィールド条件を
考えればいいと思うのですがその場合、
減衰しながらPBGの周期性に共鳴する
フィールド(エバネッセント的フィールド)が
存在するかどうかということですよね。
線型な系ではそういう共鳴条件ってないんじゃないでしょうか?
~つまり、導波路としてリジッドな境界条件になっているのでは?
(フォトニックバンドギャップとしてではなく
単なる有効屈折率としてクラッド部分の屈折率が低いという理由で
閉じ込めが成立している場合は
共鳴しないけど減衰する場合があって導波すると思いますが
(ランダム壁の導波路になっていずれは減衰するのでしょうけど)。)
物理的に考えれば、エネルギーをつぎ込めばどんな波形であれ(式が正しければ)
必ずエネルギーは保存されるはずです(そして、伝搬モードと結合する成分があれば少しは伝搬していくはずです)。
それが計算の中でどっかの項に吸収されているか(複素数を使っている場合は特に注意)、
反射されて領域に入っていかないか(カットオフが掛かっている場合)、
あるいは計算式が全くおかしいか
だと思います。全く答えに近づいていないような気もしますがアドバイスということで。
No.1
- 回答日時:
入射がどうもうまくいかないようですね。
ここはひとつ、じっくりいきましょう。まず、振動電界のみを与えるということは、電流源による励振であると考えられます。よって参考文献の式(1.65)によって電界を更新すればよいのですが、私の経験では
Ez(x,y)=Ez(x,y)+A*sin(ωt)-A*sin(ω(t-Δt))
という形にしないと電磁界パターンがおかしなものになりました。これは、式(1.65)の形では直流成分が重畳した正弦波状の電界による励振になるためであると考えています。
ともあれ、この式の形のために電界の振幅が極端に小さくなるということはありませんので、原因は他のところにあるものと思われます。
ところで、前回の質問の内容から考えて、参考文献10ページの2次元TM_FDTD法をプログラム化し、x軸方向に電磁波を伝播させているものと思われますが、その他の境界条件や媒質の配置を教えてください。
また、媒質を空気のみ(真空でいいです)として、電磁界分布が妥当な値(Hx=0となるはず)になるかどうかも確かめてください。
参考文献:宇野亨,「FDTD法による電磁界およびアンテナ解析」,コロナ社,1998.
この回答への補足
TCMさん、ご回答ありがとうございます。何度も同じような質問をして自分の勉強不足を恥じるばかりです。
TCMさんのおっしゃるとおり、宇野さんの本の2次元TM_FDTD法を元にプログラムしています。吸収境界は16ページの図1.11にあるように、配列の一番端からPMLの計算に使用しており、入射はその境界から少なくとも10セル以上は離して入射させています。実際の入射界はハードソースで、
E[50,j]=E0*f(j)*sin(ωt) E0=1.0
としていました。f(j)は導波路やフォトニック結晶導波路のコア幅にあわせた導波分布で、最大の振幅は1になっています。これで入射させると、振幅1の波が導波モードで通してあるので、減衰することなく導波路を伝播していきます。これを
E[50,j]=E[50,j]+E0*f(j)*{sin(ωt)-sin(ω(t-Δt))} E0=1.0
とすると、x座標セル50の入射位置の電界が、電界-入射-電界吸収-磁界-磁界吸収と一回の計算の後にその値を見ると0.1くらいになり、非常に小さな値で伝播していきます。入射のすぐ後に電界の値を見ると、当然ながら入射位置では振幅1の正弦波になっています。
媒質の配置は、フォトニック結晶の場合、3.4くらいの高い屈折率の媒質を真空中に適当な間隔で配置し(形は円です)、伝播するところを一筋だけ媒質を抜いています。つまり、光が伝播するところは真空中になっています。
Hx=0となるのはどこの値ですか?一応、真空中で入射電界のすぐとなりのHxを見たところ、0ではありませんでした。どういうときにHxが0となるのですか?
まだ説明不足のところがあると思いますが、ご回答お願いします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 電子レンジが物を加熱する仕組みについて 吸収できるエネルギーに限界値はあるのか知りたい。 表題の通り 4 2022/05/06 20:10
- 物理学 磁束密度Bの一様な磁界の中に、Bと垂直に速度vで電子を入射したときの電子の軌道を求めよという問題の解 2 2022/12/05 20:35
- その他(自然科学) 電磁波モデルで疑問 4 2023/07/09 14:44
- その他(暮らし・生活・行事) 緊急!電磁波を受けて耳鳴りで悩んでますが 最適な商品を紹介してもらえませんか? 私はどこからか電磁波 1 2022/05/06 21:34
- 物理学 示すように,真空中の直交座標系を考える。y平面に平行な つ領域Iと領域Iがあり,軸上の領域Iと領域I 1 2023/06/25 14:46
- 物理学 電磁波と磁界、電界の優位性について 6 2023/02/19 05:20
- その他(自然科学) 電磁波の周波数と熱について教えて下さい。 電磁波の波長とエネルギーについて、雑学として興味があります 6 2022/04/18 20:00
- 物理学 この電気磁気学の問題が解けません。自分が解くと電界を求める時に、電位を求めてから電界を求めるのと、電 3 2022/05/26 12:50
- 大学・短大 電気電子の問題です。この問題が分からないので教えて頂きたいです 幅2mm、長さ5mm、厚さ0.2mm 1 2023/06/21 22:24
- 物理学 ごめんなさい 半径 αの円形コイルの面が地球磁界の方向に平行になるように垂直に立っている。 このコイ 1 2022/07/31 21:33
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
回折と屈折の違いは何か?
-
TM偏光とTE偏光
-
『水晶玉』を覗くと反対側の景...
-
XRDの2θ/θ法について教えてくだ...
-
屈折率の小さい方から大きい方...
-
振動数は何故変化しないの?
-
薄膜の各n,kから透過率を計...
-
光の屈折によって,水を入れた...
-
反射率+透過率 が 1 になら...
-
XRD
-
安全で安価な屈折率の高い液体...
-
プリズムについてです。 『図4...
-
波の屈折と光の屈折の違いって...
-
水の中に沈んでいるコインはな...
-
波の反射角
-
散乱と屈折の違いについて
-
水の反射率(?)を求めること...
-
大学の物理実験という授業でレ...
-
プリズムによるスペクトル
-
入射角と反射角(屈折角)
おすすめ情報