
有限差分時間領域法( FDTD ) またはYee 法( 1934 年生まれの中国系アメリカ人応用数学者Kane S. Yeeにちなんで命名) は、数値解析手法であり、計算電気力学のモデリング (関連する微分方程式系の近似解の探索)に使用されます。時間領域法であるため、FDTD ソリューションは1 回のシミュレーション実行で広い周波数範囲をカバーでき、非線形材料特性を自然な方法で処理できます。
FDTD 法は、グリッドベースの微分数値モデリング法 (有限差分法)の一般的なクラスに属します。時間依存のマクスウェル方程式(偏微分形式 )は、空間および時間の偏微分に対する中心差分近似を使用して離散化されます。結果として得られる有限差分方程式は、ソフトウェアまたはハードウェアのいずれかでリープフロッグ方式で解かれます。つまり、空間ボリューム内の電界ベクトル成分が特定の瞬間に解かれ、次に同じ空間ボリューム内の磁界ベクトル成分が次の瞬間に解かれます。このプロセスは、目的の過渡的または定常状態の電磁場動作が完全に展開されるまで何度も繰り返されます。
歴史
時間依存偏微分方程式(PDE) の差分スキームは、計算流体力学の問題で長年使用されてきました。 [1]これには、空間と時間の交互に配置されたグリッド上で中心化された差分演算子を使用して 2 次精度を達成するというアイデアが含まれます。[1] 1966 年の独創的な論文[2]で発表された Kane Yee の FDTD スキームの斬新さは、 Maxwell の回転方程式の各電気ベクトル場と磁気ベクトル場の成分に対して、空間と時間の交互に配置されたグリッド上で中心化された差分演算子を適用することでした。「有限差分時間領域」という記述子と、それに対応する「FDTD」という頭字語は、1980 年にAllen Tafloveによって考案されました。 [3] 1990 年頃から、FDTD 技術は、材料構造と電磁波の相互作用 を扱う多くの科学的および工学的問題を計算的にモデル化するための主要な手段として登場しました。現在のFDTDモデリングの用途は、近直流(地球電離層導波管全体を含む超低周波地球物理学)からマイクロ波(レーダーシグネチャ技術、アンテナ、無線通信デバイス、デジタル相互接続、生物医学的画像/治療)、可視光(フォトニック結晶、ナノプラズモニクス、ソリトン、バイオフォトニクス)まで多岐にわたります。[4] 2006年には、科学および工学の文献に推定2,000件のFDTD関連の出版物が掲載されました(人気を参照)。2013年の時点で、少なくとも25の商用/独自のFDTDソフトウェアベンダー、13のフリーソフトウェア/オープンソースソフトウェアのFDTDプロジェクト、および2つのフリーウェア/クローズドソースのFDTDプロジェクトがあり、一部は商用利用不可です(外部リンクを参照)。
FDTD法とマクスウェル方程式の開発
マクスウェル方程式の FDTD 数値手法の基礎、技術開発、および将来性を理解するには、まずその歴史を検討する必要があります。次に、この分野の主要な出版物の一部を示します。
FDTDモデルと手法
マクスウェルの微分方程式を調べると、時間における電界の変化(時間微分)は空間を横切る磁界の変化(回転)に依存していることがわかります。これにより、空間内の任意の点において、時間における電界の更新値は、電界の保存値と空間における磁界の局所分布の数値回転に依存するという、基本的なFDTD時間ステップ関係が生まれます。[2]
H フィールドも同様の方法で時間ステップ化されます。空間のどの点でも、時間における H フィールドの更新値は、H フィールドの保存値と空間における E フィールドのローカル分布の数値回転に依存します。E フィールドと H フィールドの更新を繰り返すと、時間進行プロセスが発生し、対象とする連続電磁波のサンプル データ アナログが、コンピュータ メモリに保存された数値グリッド内を伝播します。

この説明は、1次元、2次元、3次元のFDTD法に当てはまります。複数の次元を考慮すると、数値回転の計算は複雑になる可能性があります。Kane Yeeの1966年の独創的な論文では、直交座標の計算グリッドの長方形の単位セルの周囲に電界と磁界のベクトル成分を空間的にずらして配置し、各電界ベクトル成分が磁界ベクトル成分のペアの中間に位置するようにし、その逆も行うことを提案しました。[2]現在Yee格子 として知られるこの方式は非常に堅牢であることが証明されており、現在の多くのFDTDソフトウェア構造の中核となっています。
さらに、Yee は、E フィールドと H フィールドの更新をずらして、E フィールドの更新が各タイムステップの途中で H フィールドの更新の間に実行され、その逆も行われるようにする、時間を進めるためのリープフロッグ方式を提案しました。[2]プラス面としては、この明示的なタイムステップ方式により、同時方程式を解く必要がなくなり、さらに散逸のない数値波動伝播が得られます。マイナス面としては、この方式では数値安定性を確保するためにタイムステップの上限が必須です。[9]その結果、特定の種類のシミュレーションでは、完了までに何千ものタイムステップが必要になる場合があります。
FDTD法の使用
マクスウェル方程式の FDTD ソリューションを実装するには、まず計算領域を確立する必要があります。計算領域とは、シミュレーションが実行される物理的な領域です。計算領域内の空間のすべてのポイントで E フィールドと H フィールドが決定されます。計算領域内の各セルの材質を指定する必要があります。通常、材質は自由空間 (空気)、金属、または誘電体のいずれかです。透磁率、誘電率、導電率が指定されている限り、任意の材質を使用できます。
表形式の分散性材料の誘電率は、FDTD法に直接代入することはできません。代わりに、複数のデバイ項、ドルーデ項、ローレンツ項、または臨界点項を使用して近似することができます。この近似は、オープンフィッティングプログラム[70]を使用して取得できますが、必ずしも物理的な意味を持つわけではありません。
計算領域とグリッド材料が確立されると、ソースが指定されます。ソースは、ワイヤ上の電流、印加電界、または衝突する平面波です。最後のケースでは、FDTDを使用して、任意の形状のオブジェクトからの光散乱、さまざまな入射角での平面周期構造、[71] [72]、および無限周期構造のフォトニックバンド構造をシミュレートできます。[73] [74]
E フィールドと H フィールドは直接決定されるため、シミュレーションの出力は通常、計算領域内の 1 点または一連の点における E フィールドまたは H フィールドです。シミュレーションは、E フィールドと H フィールドを時間的に前進させます。
シミュレーションによって返された E フィールドと H フィールドに対して処理が行われる場合があります。シミュレーションの実行中にデータ処理が行われる場合もあります。
FDTD法はコンパクトな空間領域内の電磁場を計算するが、散乱および/または放射された遠方場は近距離場から遠方場への変換によって得ることができる。[14]
FDTDモデリングの強み
すべてのモデリング手法には長所と短所があり、FDTD 法も例外ではありません。
- FDTD は、マクスウェル方程式を解くために使用される多目的モデリング手法です。直感的なので、ユーザーは使い方を簡単に理解でき、特定のモデルから何を期待できるかを知ることができます。
- FDTD は時間領域技術であり、広帯域パルス (ガウス パルスなど) をソースとして使用すると、広範囲の周波数にわたるシステムの応答を 1 回のシミュレーションで取得できます。これは、共振周波数が正確にわからないアプリケーションや、広帯域の結果が必要な場合に役立ちます。
- FDTD は、計算領域内のあらゆる場所の E フィールドと H フィールドを時間の経過に従って計算するため、モデル内の電磁場の動きをアニメーションで表示できます。このタイプの表示は、モデル内で何が起こっているかを理解し、モデルが正しく動作していることを確認するのに役立ちます。
- FDTD 技術により、ユーザーは計算領域内のすべてのポイントで材料を指定できます。さまざまな線形および非線形の誘電体および磁性材料を自然かつ簡単にモデル化できます。
- FDTD を使用すると、開口部の影響を直接的に判定できます。シールド効果を検出でき、構造の内側と外側の両方のフィールドを直接的または間接的に検出できます。
- FDTD は E フィールドと H フィールドを直接使用します。ほとんどの EMI/EMC モデリング アプリケーションは E フィールドと H フィールドを対象としているため、シミュレーションを実行してこれらの値を取得した後に、変換を行う必要がないのは便利です。
FDTDモデリングの弱点
- FDTD法では計算領域全体をグリッド化する必要があり、グリッド空間の離散化はモデル内の最小の電磁波波長と最小の幾何学的特徴の両方を解くのに十分細かくなければならないため、非常に大きな計算領域が開発され、その結果、解析時間が非常に長くなります。細長い特徴(ワイヤなど)を持つモデルは、必要となる計算領域が大きすぎるため、FDTD法でモデル化するのは困難です。固有モード展開などの方法は、Z方向に沿って細かいグリッドを必要としないため、より効率的な代替手段となります。[75]
- 材料界面における誘電率と透磁率の一意の値を決定する方法はありません。
- 空間ステップと時間ステップはCFL 条件を満たす必要があります。そうでないと、偏微分方程式を解くために使用されるリープフロッグ積分が不安定になる可能性があります。
- FDTD法は計算領域内のどこでもE/H場を直接検出します。ある距離での場の値が必要な場合、この距離によって計算領域が過度に大きくなる可能性があります。FDTD法では遠距離場の拡張が可能ですが、ある程度の後処理が必要です。[4]
- FDTD シミュレーションは計算領域内のすべてのポイントで E および H フィールドを計算するため、計算領域はコンピュータ メモリ内に保存できるように有限である必要があります。多くの場合、これはシミュレーション空間に人工的な境界を挿入することによって実現されます。このような境界によって生じるエラーを最小限に抑えるように注意する必要があります。無限の境界のない計算領域をシミュレートするために、非常に効果的な吸収境界条件 (ABC) が多数用意されています。[4]最新の FDTD 実装のほとんどは、代わりに吸収境界を実装するために完全整合層(PML) と呼ばれる特殊な吸収「材料」を使用しています。[42] [47]
- FDTD は時間領域でフィールドを前方に伝播させることによって解かれるため、媒体の電磁時間応答を明示的にモデル化する必要があります。任意の応答の場合、計算コストの高い時間畳み込みが必要になりますが、ほとんどの場合、媒体 (または分散 (光学) ) の時間応答は、再帰畳み込み (RC) 手法、補助微分方程式 (ADE) 手法、または Z 変換手法のいずれかを使用して適切かつ簡単にモデル化できます。任意の分散を簡単に処理できるマクスウェル方程式を解く別の方法は、擬似スペクトル空間領域 (PSSD)です。これは、代わりにフィールドを空間で前方に伝播します。
グリッド切り捨て技術
オープン領域FDTDモデリング問題で最も一般的に使用されるグリッド切断技術は、Mur吸収境界条件(ABC)[13] 、 Liao ABC、[16]、およびさまざまな完全整合層(PML)定式化です。[4] [43] [42] [47] MurおよびLiao技術はPMLよりも単純です。ただし、PML(技術的には境界条件自体ではなく吸収領域)は、桁違いに低い反射を提供できます。PMLの概念は、1994年にJournal of Computational Physicsで発表された独創的な論文でJ.-P. Berengerによって導入されました。[42] 1994年以来、Berengerのオリジナルの分割フィールド実装は修正され、単軸PML(UPML)、畳み込みPML(CPML)、および高次PMLに拡張されました。後者の 2 つの PML 定式化では、エバネッセント波を吸収する能力が向上しているため、原理的には、Berenger の元の定式化よりもシミュレートされた散乱構造または放射構造に近づけることができます。
PMLからの不要な数値反射を減らすために、追加の背面吸収層技術を使用することができる。[76]
人気
同じ期間の学術出版スループットの全体的な増加と、すべての計算電磁気学 (CEM) 技術に対する関心の全体的な拡大にもかかわらず、マクスウェル方程式の FDTD 計算ソリューション アプローチに対する関心が大幅に拡大した主な理由は 7 つあります。
- FDTD では行列反転は必要ありません。FDTD は完全に明示的な計算であるため、周波数領域積分方程式や有限要素電磁気学モデルのサイズを一般に 10 9未満の電磁場未知数に制限する行列反転の困難を回避できます。[4] 10 9もの未知数を持つ FDTD モデルが実行されましたが、この数に固有の上限はありません。[4]
- FDTD法は正確かつ堅牢です。FDTD法の計算における誤差の原因は十分に理解されており、非常に多様な電磁波相互作用問題に対して正確なモデルを作成できるように制限することができます。[4]
- FDTDはインパルス的な振る舞いを自然に扱います。時間領域技術であるFDTDは、電磁システムのインパルス応答を直接計算します。したがって、単一のFDTDシミュレーションで、励起スペクトル内の任意の周波数で超広帯域の時間波形または正弦波定常応答のいずれかを提供できます。[4]
- FDTDは非線形挙動を自然に扱う。時間領域技術であるFDTDは、電磁システムの非線形応答を直接計算する。これにより、古典的または半古典的な観点から非線形性を記述する補助微分方程式のセットとFDTDを自然にハイブリッド化することができる。[4]研究の最前線の一つは、FDTD古典電気力学モデルと量子電気力学から生じる現象、特にカシミール効果 などの真空変動を結合するハイブリッドアルゴリズムの開発である。[4] [77]
- FDTDは体系的なアプローチです。FDTDでは、モデル化される新しい構造を指定することは、潜在的に複雑な積分方程式の再定式化ではなく、メッシュ生成の問題に簡略化されます。たとえば、FDTDでは構造に依存するグリーン関数の計算は必要ありません。[4]
- 並列処理コンピュータアーキテクチャはスーパーコンピューティングの主流になってきました。FDTDは並列処理CPUベースのコンピュータ上で高い効率でスケーリングし、最近開発されたGPUベースのアクセラレータ技術でも非常に優れています。[4]
- コンピュータの視覚化能力は急速に向上しています。この傾向はすべての数値手法に良い影響を与えますが、特にFDTD法には有利です。FDTD法は、フィールドのダイナミクスを示すカラービデオに適したフィールド量の時間進行配列を生成します。[4]
タフロブは、これらの要因を組み合わせると、FDTDが計算電気力学の主要な手法の1つであり続けるだろうと主張している(潜在的には他のマルチフィジックス問題も同様である)。[4]
参照
参考文献
- ^ ab J. von Neumann; RD Richtmyer (1950年3月). 「流体衝撃の数値計算法」. Journal of Applied Physics . 21 (3): 232–237. Bibcode :1950JAP....21..232V. doi :10.1063/1.1699639.
- ^ abcdef Kane Yee (1966). 「等方性媒体におけるマクスウェル方程式を含む初期境界値問題の数値解」. IEEE Transactions on Antennas and Propagation . 14 (3): 302–307. Bibcode :1966ITAP...14..302Y. doi :10.1109/TAP.1966.1138693. S2CID 122712881.
- ^ ab A. Taflove (1980). 「正弦定常電磁波浸透問題への有限差分時間領域法の適用」(PDF) . IEEE Trans. Electromagn. Compat. 22 (3): 191–202. Bibcode :1980ITElC..22..191T. doi :10.1109/TEMC.1980.303879. S2CID 39236486.
- ^ abcdefghijklmno Allen TafloveおよびSusan C. Hagness (2005)。計算電気力学:有限差分時間領域法、第 3 版。Artech House Publishers。ISBN 978-1-58053-832-9。
- ^ TafloveとHagness(2005)の許可を得て改変。
- ^ リチャード・クーラント;クルト・オットー・フリードリヒス。ハンス・レヴィ (1928)。 「数学物理学の重要な部分」。Mathematische Annalen (ドイツ語)。100 (1): 32–74。ビブコード:1928MatAn.100...32C。土井:10.1007/BF01448839。JFM 54.0486.01。MR 1512478。S2CID 120760331 。
- ^ GG O'Brien、M. A Hyman、S. Kaplan (1950)。「偏微分方程式の数値解法の研究」。Journal of Mathematical Physics。29 ( 1): 223–251。doi : 10.1002 /sapm1950291223。MR 0040805 。
{{cite journal}}: CS1 maint: 複数の名前: 著者リスト (リンク) - ^ Dong-Hoa Lam (1969). 「電磁散乱問題に対する有限差分法」(PDF) .ミシシッピ州立大学、Interaction Notes . 44 .
- ^ ab A. Taflove; ME Brodwin (1975). 「時間依存マクスウェル方程式を使用した定常電磁散乱問題の数値解法」(PDF) . IEEE Transactions on Microwave Theory and Techniques . 23 (8): 623–630. Bibcode :1975ITMTT..23..623T. doi :10.1109/TMTT.1975.1128640.
- ^ A. Taflove; ME Brodwin (1975). 「マイクロ波照射を受けた人間の目のモデルにおける電磁場と誘導温度の計算」(PDF) . IEEE Transactions on Microwave Theory and Techniques . 23 (11): 888–896. Bibcode :1975ITMTT..23..888T. doi :10.1109/TMTT.1975.1128708.
- ^ R. Holland (1977). 「Threde: 自由場EMP結合および散乱コード」. IEEE Transactions on Nuclear Science . 24 (6): 2416–2421. Bibcode :1977ITNS...24.2416H. doi :10.1109/TNS.1977.4329229. S2CID 35395821.
- ^ KS Kunz; KM Lee (1978). 「複雑な過渡電磁環境に対する航空機の外部応答の 3 次元差分法」IEEE Trans. Electromagn. Compat. 20 (2): 333–341. doi :10.1109/TEMC.1978.303727. S2CID 31666283.
- ^ ab G. Mur (1981). 「時間領域電磁場方程式の有限差分近似のための吸収境界条件」. IEEE Trans. Electromagn. Compat. 23 (4): 377–382. doi :10.1109/TEMC.1981.303970. S2CID 25768246.
- ^ ab KR Umashankar; A. Taflove (1982). 「複雑な物体の電磁散乱を解析する新しい方法」(PDF) . IEEE Trans. Electromagn. Compat. 24 (4): 397–405. Bibcode :1982ITElC..24..397U. doi :10.1109/TEMC.1982.304054. S2CID 37962500.
- ^ A. Taflove; KR Umashankar (1983). 「一般的な3次元散乱体のレーダー断面積」(PDF) . IEEE Trans. Electromagn. Compat. 25 (4): 433–440. doi :10.1109/TEMC.1983.304133. S2CID 40419955.
- ^ ab ZP Liao; HL Wong; BP Yang; YF Yuan (1984). 「過渡波解析のための透過境界」Scientia Sinica、シリーズ A。27 : 1063–1076。
- ^ W. Gwarek (1985). 「任意形状の平面回路の解析 - 時間領域アプローチ」. IEEE Transactions on Microwave Theory and Techniques . 33 (10): 1067–1072. Bibcode :1985ITMTT..33.1067G. doi :10.1109/TMTT.1985.1133170.
- ^ DH Choi; WJ Hoefer (1986). 「有限差分時間領域法と固有値問題への応用」. IEEE Transactions on Microwave Theory and Techniques . 34 (12): 1464–1470. Bibcode :1986ITMTT..34.1464C. doi :10.1109/TMTT.1986.1133564.
- ^ GA Kriegsmann; A. Taflove; KR Umashankar (1987). 「表面放射境界条件アプローチを使用した電磁波散乱の新しい定式化」(PDF) . IEEE Transactions on Antennas and Propagation . 35 (2): 153–161. Bibcode :1987ITAP...35..153K. doi :10.1109/TAP.1987.1144062.
- ^ TG Moore; JG Blaschak; A. Taflove; GA Kriegsmann (1988). 「放射境界演算子の理論と応用」(PDF) . IEEE Transactions on Antennas and Propagation . 36 (12): 1797–1812. Bibcode :1988ITAP...36.1797M. doi :10.1109/8.14402.
- ^ KR Umashankar; A. Taflove; B. Beker (1987). 「任意の形状の空洞内の結合ワイヤ上の誘導電流の計算と実験的検証」(PDF) . IEEE Transactions on Antennas and Propagation . 35 (11): 1248–1257. Bibcode :1987ITAP...35.1248U. doi :10.1109/TAP.1987.1144000.
- ^ A. Taflove; KR Umashankar; B. Beker; FA Harfoush; KS Yee (1988). 「厚い導電性スクリーンの狭いスロットと重ね合わせたジョイントを貫通する電磁場の詳細な FDTD 解析」(PDF) . IEEE Transactions on Antennas and Propagation . 36 (2): 247–257. Bibcode :1988ITAP...36..247T. doi :10.1109/8.1102.
- ^ TG Jurgens; A. Taflove; KR Umashankar; TG Moore (1992). 「曲面の有限差分時間領域モデリング」(PDF) . IEEE Transactions on Antennas and Propagation . 40 (4): 357–366. Bibcode :1992ITAP...40..357J. doi :10.1109/8.138836.
- ^ DM Sullivan; OP Gandhi; A. Taflove (1988). 「人体モデルにおける電磁波吸収の計算における有限差分時間領域法の使用」(PDF) . IEEE Transactions on Biomedical Engineering . 35 (3): 179–186. doi :10.1109/10.1360. PMID 3350546. S2CID 20350396.
- ^ X. Zhang; J. Fang; KK Mei; Y. Liu (1988). 「時間領域有限差分法によるマイクロストリップの分散特性の計算」. IEEE Transactions on Microwave Theory and Techniques . 36 (2): 263–267. Bibcode :1988ITMTT..36..263Z. doi :10.1109/22.3514.
- ^ T. Kashiwa; I. Fukai (1990). 「FDTD法による電子分極に伴う分散特性の解析」.マイクロ波・光技術レター. 3 (6): 203–205. doi :10.1002/mop.4650030606.
- ^ R. Luebbers; F. Hunsberger; K. Kunz; R. Standler; M. Schneider (1990). 「分散材料の周波数依存有限差分時間領域定式化」. IEEE Trans. Electromagn. Compat. 32 (3): 222–227. doi :10.1109/15.57116.
- ^ RM Joseph; SC Hagness; A. Taflove (1991). 「吸収を伴う線形分散媒体におけるマクスウェル方程式の直接時間積分によるフェムト秒電磁パルスの散乱と伝播」(PDF)。Optics Letters。16 (18): 1412–4。Bibcode :1991OptL...16.1412J。doi :10.1364/OL.16.001412。PMID 19776986。
- ^ JG Maloney; GS Smith; WR Scott Jr. (1990). 「有限差分時間領域法を用いた単純なアンテナからの放射の正確な計算」. IEEE Transactions on Antennas and Propagation . 38 (7): 1059–1068. Bibcode :1990ITAP...38.1059M. doi :10.1109/8.55618. S2CID 31583883.
- ^ DS Katz; A. Taflove ; MJ Piket-May ; KR Umashankar (1991). 「ホーンアンテナを含むシステムからの電磁波放射の FDTD 解析」(PDF) . IEEE Transactions on Antennas and Propagation . 39 (8): 1203–1212. Bibcode :1991ITAP...39.1203K. doi :10.1109/8.97356.
- ^ PA Tirkas; CA Balanis (1991). 「ホーンアンテナによる放射の有限差分時間領域法」。アンテナおよび伝播協会シンポジウム 1991 ダイジェスト。第 3 巻。pp. 1750–1753。doi : 10.1109 /APS.1991.175196。ISBN 978-0-7803-0144-3. S2CID 122038624。
- ^ E. Sano; T. Shibata (1990). 「ピコ秒光伝導スイッチの全波解析」. IEEE Journal of Quantum Electronics . 26 (2): 372–377. Bibcode :1990IJQE...26..372S. doi :10.1109/3.44970.
- ^ SM El-Ghazaly; RP Joshi; RO Grondin (1990). 「サブピコ秒光伝導スイッチモデリングにおける電磁気と輸送の考慮」. IEEE Transactions on Microwave Theory and Techniques . 38 (5): 629–637. Bibcode :1990ITMTT..38..629E. doi :10.1109/22.54932.
- ^ PM Goorjian; A. Taflove (1992). 「フェムト秒電磁ソリトンの伝播と散乱のための非線形分散媒体におけるマクスウェル方程式の直接時間積分」(PDF)。Optics Letters。17 (3): 180–182。Bibcode :1992OptL...17..180G。doi :10.1364/OL.17.000180。PMID 19784268。
- ^ RW Ziolkowski; JB Judkins (1993). 「有限応答時間を示す非線形カー媒体における超短光パルスの自己集束の全波ベクトルマクスウェル方程式モデリング」。Journal of the Optical Society of America B 。10 ( 2 ): 186–198。Bibcode :1993JOSAB..10..186Z。doi :10.1364/JOSAB.10.000186。
- ^ RM Joseph; PM Goorjian; A. Taflove (1993). 「フェムト秒電磁ソリトンの伝搬と散乱のための 2 次元誘電体導波路におけるマクスウェル方程式の直接時間積分」(PDF) . Optics Letters . 18 (7): 491–3. Bibcode :1993OptL...18..491J. doi :10.1364/OL.18.000491. PMID 19802177.
- ^ RM Joseph; A. Taflove (1994). 「FDTD マクスウェル方程式モデリングによる空間ソリトン偏向メカニズム」(PDF) . IEEE Photonics Technology Letters . 2 (10): 1251–1254. Bibcode :1994IPTL....6.1251J. doi :10.1109/68.329654. S2CID 46710331.
- ^ W. Sui; DA Christensen; CH Durney (1992). 「2次元FDTD法のアクティブおよびパッシブ集中素子によるハイブリッド電磁システムへの拡張」IEEE Transactions on Microwave Theory and Techniques . 40 (4): 724–730. Bibcode :1992ITMTT..40..724S. doi :10.1109/22.127522.
- ^ B. Toland; B. Houshmand; T. Itoh (1993). 「FDTD法による非線形アクティブ領域のモデリング」. IEEE Microwave and Guided Wave Letters . 3 (9): 333–335. doi :10.1109/75.244870. S2CID 27549555.
- ^ Aoyagi, PH; Lee, JF; Mittra, R. (1993). 「ハイブリッド Yee アルゴリズム/スカラー波動方程式アプローチ」. IEEE Transactions on Microwave Theory and Techniques . 41 (9): 1593–1600. Bibcode :1993ITMTT..41.1593A. doi :10.1109/22.245683.
- ^ VA Thomas; ME Jones; MJ Piket-May ; A. Taflove; E. Harrigan (1994). 「FDTD 高速電子回路設計のサブグリッド モデルとしての SPICE 集中回路の使用」(PDF) . IEEE Microwave and Guided Wave Letters . 4 (5): 141–143. doi :10.1109/75.289516. S2CID 32905331.
- ^ abcd J. Berenger (1994). 「電磁波吸収のための完全整合層」(PDF) . Journal of Computational Physics . 114 (2): 185–200. Bibcode :1994JCoPh.114..185B. doi :10.1006/jcph.1994.1159.
- ^ DS Katz; ET Thiele; A. Taflove (1994). 「FDTD メッシュの Berenger PML 吸収境界条件の検証と 3 次元への拡張」(PDF) . IEEE Microwave and Guided Wave Letters . 4 (8): 268–270. doi :10.1109/75.311494. S2CID 10156811.
- ^ CE Reuter; RM Joseph; ET Thiele; DS Katz; A. Taflove (1994). 「FDTD シミュレーションにおける導波構造の終端に対する超広帯域吸収境界条件」(PDF) . IEEE Microwave and Guided Wave Letters . 4 (10): 344–346. doi :10.1109/75.324711. S2CID 24572883.
- ^ WC Chew; WH Weedon (1994). 「引き伸ばされた座標を持つ修正マクスウェル方程式による 3D 完全整合媒体」.マイクロ波および光技術レター. 7 (13): 599–604.書誌コード:1994MiOTL...7..599C. doi :10.1002/mop.4650071304.
- ^ abc SD Gedney (1996). 「FDTD 格子の切り捨てのための異方性完全整合層吸収媒体」. IEEE Transactions on Antennas and Propagation . 44 (12): 1630–1639. Bibcode :1996ITAP...44.1630G. doi :10.1109/8.546249.
- ^ ZS Sacks、DM Kingsland、R. Lee、JF Lee ( 1995)。「吸収境界条件として使用する完全に整合した異方性吸収体」。IEEE Transactions on Antennas and Propagation。43 ( 12): 1460–1463。Bibcode :1995ITAP...43.1460S。doi : 10.1109/8.477075。
- ^ QH Liu (1997). 「擬似スペクトル時間領域 (PSTD) 法: マクスウェル方程式の解の新しいアルゴリズム」IEEE アンテナおよび伝搬学会国際シンポジウム 1997. ダイジェスト. 第 1 巻. pp. 122–125. doi :10.1109/APS.1997.630102. ISBN 978-0-7803-4178-4. S2CID 21345353。
- ^ OM Ramahi (1997). 「FDTDシミュレーションにおける相補演算子法」. IEEE Antennas and Propagation Magazine . 39 (6): 33–45. Bibcode :1997IAPM...39...33R. doi :10.1109/74.646801.
- ^ JG Maloney、MP Kesler (1998)。「周期構造の解析」。『計算電気力学の進歩: 有限差分時間領域法』第 6 章、A. Taflove 編、Artech House、出版社。
- ^ AS Nagra; RA York (1998). 「非線形吸収および利得媒体における波動伝搬の FDTD 解析」. IEEE Transactions on Antennas and Propagation . 46 (3): 334–340. Bibcode :1998ITAP...46..334N. doi :10.1109/8.662652.
- ^ SC Hagness; A. Taflove; JE Bridges (1998). 「乳がん検出用パルスマイクロ波共焦点システムの2次元FDTD解析:固定焦点およびアンテナアレイセンサー」(PDF) . IEEE Transactions on Biomedical Engineering . 45 (12): 1470–1479. doi :10.1109/10.730440. PMID 9835195. S2CID 6169784.
- ^ JB Schneider; CL Wagner (1999). 「FDTD 分散の再考: 光より速い伝搬」. IEEE Microwave and Guided Wave Letters . 9 (2): 54–56. CiteSeerX 10.1.1.77.9132 . doi :10.1109/75.755044.
- ^ F. Zhen; Z. Chen; J. Zhang (2000). 「3次元無条件安定有限差分時間領域法の開発に向けて」IEEE Transactions on Microwave Theory and Techniques . 48 (9): 1550–1558. Bibcode :2000ITMTT..48.1550Z. doi :10.1109/22.869007.
- ^ F. Zheng; Z. Chen (2001). 「無条件に安定した 3-D ADI-FDTD 法の数値分散解析」. IEEE Transactions on Microwave Theory and Techniques . 49 (5): 1006–1009. Bibcode :2001ITMTT..49.1006Z. doi :10.1109/22.920165.
- ^ JA Roden; SD Gedney (2000). 「畳み込み PML (CPML): 任意の媒体に対する CFS-PML の効率的な FDTD 実装」. Microwave and Optical Technology Letters . 27 (5): 334–339. doi :10.1002/1098-2760(20001205)27:5<334::AID-MOP14>3.0.CO;2-A. 2013-01-05 にオリジナルからアーカイブされました。
- ^ T. Rylander; A. Bondeson (2000). 「マクスウェル方程式のための安定したFDTD-FEMハイブリッド法」. Computer Physics Communications . 125 (1–3): 75–82. doi :10.1016/S0010-4655(99)00463-4.
- ^ M. Hayakawa、T. Otsuyama (2002)。「不均質な電離層下導波管モデルにおけるELF波伝搬のFDTD解析」。ACESジャーナル。17 :239-244。2012年5月27日時点のオリジナルよりアーカイブ。
- ^ JJ Simpson; A. Taflove (2002). 「地球の反対側の ELF 伝搬とシューマン共鳴の 2 次元 FDTD モデル」(PDF) . IEEE Antennas and Wireless Propagation Letters . 1 (2): 53–56. Bibcode :2002IAWPL...1...53S. CiteSeerX 10.1.1.694.4837 . doi :10.1109/LAWP.2002.805123. S2CID 368964. 2010-06-17 に オリジナル(PDF)からアーカイブ。
- ^ H. De Raedt、K. Michielsen 、 JS Kole、MT Figge (2003)。「チェビシェフ法によるマクスウェル方程式の解法:ワンステップ有限差分時間領域アルゴリズム」IEEE Transactions on Antennas and Propagation。51 ( 11 ) : 3155–3160。arXiv : physics/0208060。Bibcode : 2003ITAP...51.3155D。doi : 10.1109/TAP.2003.818809。S2CID 119095479 。
- ^ I. Ahmed; EK Chua; EP Li; Z. Chen (2008). 「3次元無条件安定LOD-FDTD法の開発」. IEEE Transactions on Antennas and Propagation . 56 (11): 3596–3600. Bibcode :2008ITAP...56.3596A. doi :10.1109/TAP.2008.2005544. S2CID 31351974.
- ^ 谷口 勇; 馬場 勇; 長岡 直樹; 雨谷 明 (2008). 「FDTD 計算のための改良された細線表現」. IEEE Transactions on Antennas and Propagation . 56 (10): 3248–3252. Bibcode :2008ITAP...56.3248T. doi :10.1109/TAP.2008.929447. S2CID 29617214.
- ^ RMS de Oliveira; CLSS Sobrinho (2009). 「有限差分時間領域法による変電所の落雷シミュレーションのための計算環境」. IEEE Transactions on Electromagnetic Compatibility . 51 (4): 995–1000. doi :10.1109/TEMC.2009.2028879.
- ^ FI Moxley III; T. Byrnes; F. Fujiwara; W. Dai (2012). 「N 体相互作用ハミルトニアンに対する一般化有限差分時間領域量子法」. Computer Physics Communications . 183 (11): 2434–2440. Bibcode :2012CoPhC.183.2434M. doi :10.1016/j.cpc.2012.06.012.
- ^ FI Moxley III; DT Chuss; W. Dai (2013). 「非線形シュレディンガー方程式を解くための一般化された有限差分時間領域スキーム」. Computer Physics Communications . 184 (8): 1834–1841. Bibcode :2013CoPhC.184.1834M. doi :10.1016/j.cpc.2013.03.006.
- ^ Frederick Moxley 他 (2014)。現代数学: 連続および離散動的システムの数学。アメリカ数学会。ISBN 978-0-8218-9862-8。
- ^ RMS de Oliveira; RR Paiva (2021). 「最小二乗有限差分時間領域」。IEEE Transactions on Antennas and Propagation . 69 (9): 6111–6115. Bibcode :2021ITAP...69.6111D. doi :10.1109/TAP.2021.3069576. S2CID 234307029.
- ^ 「誘電関数のフィッティング」。
- ^ I. Valuev; A. Deinega; S. Belousov (2008). 「有限差分時間領域法による斜入射周期構造解析のための反復手法」. Opt. Lett . 33 (13): 1491–3. Bibcode :2008OptL...33.1491V. doi :10.1364/ol.33.001491. PMID 18594675.
- ^ A. Aminian、Y. Rahmat-Samii (2006)。「スペクトルFDTD:周期的構造上の斜め入射平面波の解析のための新しい手法」IEEE Transactions on Antennas and Propagation。54 ( 6): 1818–1825。Bibcode : 2006ITAP ...54.1818A。doi :10.1109/tap.2006.875484。S2CID 25120679。
- ^ A. Deinega; S. Belousov; I. Valuev (2009). 「層状周期構造のためのハイブリッド転送マトリックスFDTD法」. Opt. Lett . 34 (6): 860–2. Bibcode :2009OptL...34..860D. doi :10.1364/ol.34.000860. PMID 19282957. S2CID 27742034.
- ^ Y. Hao; R. Mittra (2009). メタマテリアルのFDTDモデリング:理論と応用。Artech House Publishers。
- ^ D. Gallagher (2008). 「フォトニクス CAD の成熟」(PDF) . LEOS ニュースレター.
- ^ A. Deinega; I. Valuev (2011). 「層状周期構造のPML吸収境界の長時間挙動」Comput. Phys. Commun . 182 (1): 149–151. Bibcode :2011CoPhC.182..149D. doi :10.1016/j.cpc.2010.06.006.
- ^ SG Johnson、「カシミール相互作用を計算するための数値的手法」、Casimir Physics (D. Dalvit、P. Milonni、D. Roberts、F. da Rosa 編)、Lecture Notes in Physicsの第 834 巻、第 6 章、pp. 175–218、ベルリン: Springer、2011 年 6 月。
さらに読む
Nature Milestones: Photonsの次の記事は、マクスウェル方程式に関連した FDTD 法の歴史的重要性を説明しています。
- David Pile (2010 年 5 月)。「マイルストーン 2 (1861) マクスウェル方程式」。Nature Milestones: Photons。doi : 10.1038 /nmat2639。2010年6 月 17日閲覧。
マクスウェル方程式の出版 150 周年を記念した 2015 年 1 月号のNature Photonics特集号に掲載された Allen Taflove のインタビュー「数値解法」。このインタビューでは、FDTD の開発がマクスウェルの電気力学理論の 150 年の歴史とどのように結びついているかについて触れています。
- ネイチャーフォトニクスのインタビュー
以下の大学レベルの教科書は、FDTD 法の一般的な入門書として最適です。
- Karl S. Kunz、Raymond J. Luebbers (1993)。電磁気学のための有限差分時間領域法。CRC Press。ISBN 978-0-8493-8657-2. 2007年12月10日時点のオリジナルよりアーカイブ。2006年8月5日閲覧。
- Allen Taflove、Susan C. Hagness (2005)。計算電気力学: 有限差分時間領域法、第 3 版。Artech House Publishers。ISBN 978-1-58053-832-9。
- Wenhua Yu 、 Raj Mittra、Tao Su、Yongjun Liu、Xiaoling Yang (2006)。並列有限差分時間領域法。Artech House Publishers。ISBN 978-1-59693-085-8。
- John B. Schneider (2010)。FDTD 法の理解。オンラインで入手可能。
外部リンク
フリーソフトウェア/オープンソースソフトウェアFDTD プロジェクト:
- FDTD++: 高度なフル機能の FDTD ソフトウェア、洗練された材料モデル、定義済みのフィット、ディスカッション/サポート フォーラム、電子メール サポート
- openEMS (C++ で記述され、 Matlab / Octaveインターフェイスを使用する、完全な 3D 直交座標および円筒形の段階的メッシュ EC-FDTD ソルバー)
- pFDTD (Se-Heon Kim が開発した 3D C++ FDTD コード)
- JFDTD (Jeffrey M. McMahon がナノフォトニクス用に開発した 2D/3D C++ FDTD コード)
- WOLFSIM 2008-07-02 にWayback Machineでアーカイブ(NCSU) (2-D)
- Meep ( MIT、2D/3D/円筒並列 FDTD)
- (ジオ)レーダーFDTD
- bigboy (メンテナンスされていないため、リリース ファイルはありません。ソースは CVS から取得する必要があります)
- C++ の並列 (MPI&OpenMP) FDTD コード (Zs. Szabó によって開発)
- Fortran 90 の FDTD コード
- 2D EM 波シミュレーション用の C FDTD コード
- Angora (Ilker R. Capoglu が管理する 3D 並列 FDTD ソフトウェア パッケージ)
- GSvit (グラフィック カード コンピューティングをサポートする 3D FDTD ソルバー、C 言語で記述、グラフィカル ユーザー インターフェイス XSvit が利用可能)
- gprMax (オープンソース (GPLv3)、GPR 用に開発された Python/Cython の 3D/2D FDTD モデリング コードですが、一般的な EM モデリングにも使用できます。)
フリーウェア/クローズドソースFDTD プロジェクト (一部は商用利用不可):
- EMTL (電磁テンプレート ライブラリ) (電磁シミュレーション用の無料の C++ ライブラリ。現在のバージョンでは主に FDTD を実装しています)。
