- 締切済み
FDTD法でのソフトソース入射
motsuanの回答
- motsuan
- ベストアンサー率40% (54/135)
物理的には正しい計算になっているでしょうか？ 以前の質問では単純に減衰しているとお書きになっていますが、 減衰には吸収による減衰と導波路からの放射による減衰があると思います。 ・（吸収媒質で無い場合）エネルギーは全体で保たれているでしょうか（境界から出ていく量も含めて）？ ・バンドギャップは十分開いているでしょうか？ ・これは私も興味あるのですがＰＢＧの導波路って カットオフとか持っていないのでしょうか （通常の誘電体光導波路では対称な場合はカットオフないのですが、 同じに考えていいのでしょうか？）．．．と書いていたら、 ご自分でも疑問に思ってらっしゃるようですね。 誘電体多層膜でサンドイッチした構造のときのフィールド条件を 考えればいいと思うのですがその場合、 減衰しながらＰＢＧの周期性に共鳴する フィールド（エバネッセント的フィールド）が 存在するかどうかということですよね。 線型な系ではそういう共鳴条件ってないんじゃないでしょうか？ ～つまり、導波路としてリジッドな境界条件になっているのでは？ （フォトニックバンドギャップとしてではなく 単なる有効屈折率としてクラッド部分の屈折率が低いという理由で 閉じ込めが成立している場合は 共鳴しないけど減衰する場合があって導波すると思いますが （ランダム壁の導波路になっていずれは減衰するのでしょうけど）。） 物理的に考えれば、エネルギーをつぎ込めばどんな波形であれ（式が正しければ） 必ずエネルギーは保存されるはずです（そして、伝搬モードと結合する成分があれば少しは伝搬していくはずです）。 それが計算の中でどっかの項に吸収されているか（複素数を使っている場合は特に注意）、 反射されて領域に入っていかないか（カットオフが掛かっている場合）、 あるいは計算式が全くおかしいか だと思います。全く答えに近づいていないような気もしますがアドバイスということで。
関連するQ&A
- ＦＤＴＤ法における入射について
現在ＦＤＴＤ法によるフォトニック結晶光波回路の解析を行っているのですが、入射が思うようにいきません。 僕が行っている解析はＥｚ、Ｈｘ、Ｈｙ成分で構成されるＴＭ波で、入射のさせかたは、解析領域におけるx＝０の位置でＥｚ成分のみについて、 Ez(x,y) = sin((y-y0)π/W) * sin(2πft) y0:波の中心 W :波の幅 f :入射波の周波数 として、時間的にも位置的にも正弦波の連続波として入射させています。Ｈｘ，Ｈｙについてはなにも行っていません。 この入射の計算は、ＦＤＴＤの電磁界を計算する前に入れており、入射の計算を行ったあとに電界と磁界をそれぞれ計算するようにしています。 しかしこのままでは発散していまってので、ＦＤＴＤの電界の計算は入射位置を除いた（入射位置がEz(1,j)とするとEz(2,j)から)位置から行っています。 一応このやり方で導波路に波は伝搬するのですが、入射初期の波頭は減衰してなくなってしまいます。しかも遠くに伝搬するにしたがって波はどんどん減衰していき、定常応答を得るまでにかなりの時間を要します。波の振る舞いについてはこれで正しいのでしょうか？ あと、時間と位置的にガウス波の孤立パルスも入射させてみたのですが、すぐに減衰して広がってしまって最初の振幅を保ったまま伝搬しません。なにかおかしい所があるのでしょうか？それともこれで正しいのでしょうか？ＦＤＴＤに関する論文や本などには、入射についてあまり詳しく書いてないので確かめようがありません。どなたかアドバイスしてくれれば幸いです。 質問にまとまりがなく、わかりにくくて申し訳ありませんが、なにかちょっとしたことでもかまいませんので、ご回答をよろしくお願いします。
- ベストアンサー
- 物理学
- FDTDによるマイクロストリップ共振器の解析
FDTDでマイクロストリップライン共振器の解析をおこないたいのですが、うまくいきません。 モデリングで、X,Y,Z方向の媒質IDを指定し、完全導体の場合、電界は０、誘電体の場合はその誘電体の誘電率で計算をおこなっています。 励振にはガウシャンパルスを使っています。 空気と完全導体や、誘電体と完全導体の境界等で何か処理が必要なんですか？
- 締切済み
- 物理学
- マクスウェル方程式？？？
高周波磁界は高周波電界よりも９０°位相が進んでいることを証明したいのです。 条件として電界はｘ方向のみ、磁界はｙ方向のみの成分を有していて 電界 Ｅ＝Ｅ0 exp(jwt) Ｅ0は振幅 磁界 Ｈ＝Ｈ0 exp(jwt) Ｈ0は振幅 で振動しているとする。 直交座標で解きたく電磁波関係の本を見ているのですがイマイチ理解ができないので、よろしくお願いします。
- ベストアンサー
- 物理学
- 正弦波の合成について
二つの振幅と位相が異っている、同周期の振動波があり、 これを合成する方法がわかりません。 二つの振動波は、Y1=A1sinωt、Y2=A2sin（ωt＋φ）です。 導き方を詳しく教えて下さい。どうか宜しくお願いします。
- ベストアンサー
- 物理学
- FDTD法を用いた電磁界解析にOpenMPを導入
FDTD法を用いた電磁界解析にOpenMPを導入するにはどうしたらよいのでしょうか？ ちなみに３次元のプログラムを作成しており、電磁界解析のエンジン部分は下記の様になっています。 for( kt=0; kt<=it ; kt++) { /* ※※※forの括弧開始※※※ */ /*********** 1ステップ前の電磁界[jt-1]の値の定式 ***********/ for( jx=0; jx<=ix; jx++ ) for( jy=0; jy<=iy; jy++ ) for( jz=0; jz<=iz; jz++ ) { exg[jx][jy][jz] = ex[jx][jy][jz]; eyg[jx][jy][jz] = ey[jx][jy][jz]; ezg[jx][jy][jz] = ez[jx][jy][jz]; } /*********** 磁界(hx)の計算 ***********/ for( jx=0; jx<=ix; jx++ ) for( jy=0; jy<=iy-1; jy++ ) for( jz=0; jz<=iz-1; jz++ ) { double dey = (ey[jx][jy][jz+1]-ey[jx][jy][jz])/dz; double dez = (ez[jx][jy+1][jz]-ez[jx][jy][jz])/dy; hx[jx][jy][jz] = hx[jx][jy][jz]+dt/amuy[jx][jy][jz]*(dey-dez); } /*********** 磁界(hy)の計算 ***********/ for( jx=0; jx<=ix-1; jx++ ) for( jy=0; jy<=iy; jy++ ) for( jz=0; jz<=iz-1; jz++ ) { double dez = (ez[jx+1][jy][jz]-ez[jx][jy][jz])/dx; double dex = (ex[jx][jy][jz+1]-ex[jx][jy][jz])/dz; hy[jx][jy][jz] = hy[jx][jy][jz]+dt/amuz[jx][jy][jz]*(dez-dex); } /*********** 磁界(hz)の計算 ***********/ for( jx=0; jx<=ix-1; jx++ ) for( jy=0; jy<=iy-1; jy++ ) for( jz=0; jz<=iz; jz++ ) { double dex = (ex[jx][jy+1][jz]-ex[jx][jy][jz])/dy; double dey = (ey[jx+1][jy][jz]-ey[jx][jy][jz])/dx; hz[jx][jy][jz] = hz[jx][jy][jz]+dt/amuz[jx][jy][jz]*(dex-dey); } /*********** 励振 ***********/ for( jx=1; jx<=ix-1; jx++ ) for( jz=1; jz<=iz-2; jz++ ) { ez[jx][80][jz] = sin( 2.0*pi*(kt%100)/100 )*sin(pi*float(jx)/float(ix)); } /*********** 電界(Ez)の計算 ***********/ for( jx=1; jx<=ix-1; jx++ ) for( jy=1; jy<=iy-1; jy++ ) for( jz=0; jz<=iz-1; jz++ ) { double term1 = (1.0-sigmax[jx][jy][jz]*dt/(2.0*epsirx[jx][jy][jz])) /(1.0+sigmax[jx][jy][jz]*dt/(2.0*epsirx[jx][jy][jz])); double term2 = 1.0/(1.0+sigmax[jx][jy][jz]*dt/(2.0*epsirx[jx][jy][jz])); double dhy = (hy[jx][jy][jz]-hy[jx-1][jy][jz])/dx; double dhx = (hx[jx][jy][jz]-hx[jx][jy-1][jz])/dy; ez[jx][jy][jz] = term1*ez[jx][jy][jz]+dt/epsirx[jx][jy][jz]*term2*(dhy-dhx); } ・ ・ ・ この後電界の計算が同様に続きます。 どなたかこのプログラムにOpenMPを導入するにはどうすればよいか教えていただきたいです。
- 締切済み
- 物理学
- 磁界の式による電界の式の導出
電磁気学の問題を解いていたところ以下の問題が分かりませんでした。自分なりに答えを出したものの、これで正解なのかどうか分からないため、皆さんよろしくお願いします。 真空中をz方向に伝播する電磁波の磁界がH(z,t)=Asin(ωt-bz)iであったとする。ただしi,j,kをそれぞれx,y,z方向の単位ベクトルとする。また、真空中の特性インピーダンスをηとする。 ・電界E(z,t)の式を導出せよ。 ｜E｜=｜H｜η=Aη 電界Eは磁界Hに対して直交なので、 E(z,t)=Aηsin(ωt-bz)j
- 締切済み
- 地学
- 磁界の式による電界の式の導出について
電磁気学の問題を解いていたところ以下の問題が分かりませんでした。自分なりに答えを出したものの、これで正解なのかどうか分からないため、皆さんよろしくお願いします。 真空中をz方向に伝播する電磁波の磁界がH(z,t)=Asin(ωt-bz)iであったとする。ただしi,j,kをそれぞれx,y,z方向の単位ベクトルとする。また、真空中の特性インピーダンスをηとする。 ・電界E(z,t)の式を導出せよ。 ｜E｜=｜H｜η=Aη 電界Eは磁界Hに対して直交なので、 E(z,t)=Aηsin(ωt-bz)j
- ベストアンサー
- 物理学
- 波のグラフについて
固定された反射板による 音波の反射を考える 図は入射した音波の進行方向をｘ軸 ｘ方向の変位を正としてt=0 における入射波を示したもの入射波は正弦曲線であらわされ周期Ｔとする、反射板で固定端反射されるものとする。 (1) 入射波に対する反射波の波形をかけ (2) 図の状態から時間が経過して入射波と反射波の合成波の変位がどのｘについても０になる最初の時刻を求めよ。 以上が問題文に記されています 反射波の作図は下記のy-xグラフにて太線が入射波なのですが、薄い線が反射波。 あるxにおいて反射波は振幅をＡとしてy=Asin{2π(t/T-x/λ)+π/3}で 入射はのxにおける波の式はy'=Acos2π(t/T-x/λ)で 合成波Y=y+y'=・・・・・・・・・・・ 考え方が正しいか教えていただけると助かります
- 締切済み
- 物理学
- アンテナの電界強度について
今、自分は大学でアンテナについて研究し始めました。 その中で、半波長ダイポールアンテナについて質問させていただきたいと思い、投稿いたしました。 (質問) 半波ダイポールアンテナをi=I sin ωt という電流で励振させたとき、中心Oから両端へ向かって、それぞれxだけはなれた距離の点A,Bに線素ｄｘを考える。A点の電流iaとB点の電流ibは、 ia=I cos mx sin ωt ib=I cos mx sin ωt m=2π/λ となる。 半波長ダイポールから十分遠い距離の点Pにia,ibによって生じる電界強度をｄea,debとして求め、その電界強度を合成して、合成電界強度deを求めます。そのdeを[0,λ/4]の範囲で積分してタブレット全体によりP点に生じる電界強度eを求めました。その結果、 e=j(60I/r)〔｛cos(π/2 cosθ）｝/sinθ〕sin ωt となるところまでは 分ったのですが・・・ 本論はここからなんです 半波長ダイポールの放射抵抗を求めるとき、電界強度Eθを求めるって書いてあって、その教科書には Eθ＝| e | となり Eθ＝(60I/r)〔｛cos(π/2 cosθ）｝/sinθ〕 だと書いてあるのですが、イマイチ意味がわかりません。説明が下手で申し訳ないのですが、もし詳しいかたいらっしゃいましたら教えてください。自分的にはEθはsin ωtの振幅を意味してるので、そこらへんから来てるのかと思うのですが。是非ともアドバイスお願いします。
- ベストアンサー
- 物理学