|
|
| 1行目: |
1行目: |
| − | '''電磁場解析'''（でんじばかいせき、{{Lang-en-short|electromagnetic field analysis}}）とは、[[マクスウェルの方程式]]を解くことにより、対象物と[[電磁場]]の[[相互作用]]を[[解析]]することである。過去には、マクスウェルの方程式から導出される[[偏微分方程式]]を解析的に（手計算にて）計算することを指していたが、現在は専ら[[コンピュータ]]によって計算することを指し、[[積分方程式]]を解く解法もある。
| + | {{テンプレート:20180815sk}} __NOINDEX__ |
| − | | |
| − | [[工学]]分野では、'''電磁界解析'''という。電磁場解析には、[[静電場]]（静電界）解析、[[静磁場]]（静磁界）解析、[[電磁誘導]]解析、[[電磁波]]解析等が含まれる。このうち、電磁波解析は[[高周波回路]]や[[無線通信]]用回路、[[アンテナ]]や[[レーダー]]等の[[設計]]・解析、[[電磁両立性|電磁環境適合性]] ({{Lang|en|EMC}}) [[回折格子]]などに使用される。また、比較的低周波（数十Hz - 数百Hz）の[[磁界]]解析は、[[モーター]]などの回転器や{{仮リンク|リニアアクチュエータ|en|Linear actuator}}の設計などに用いられる。
| |
| − | | |
| − | == 計算電磁気学 (CEM) ==
| |
| − | コンピュータを用いて電磁場の解析を行う分野を '''計算電磁気学'''（{{Lang-en-short|computational electromagnetics}}; {{Lang|en|CEM}}) と呼ぶ<ref>電気学会編『計算電磁気学』, 2003, ISBN 978-4563067236</ref>。
| |
| − | | |
| − | [[散乱]]・[[電磁放射|放射]]・[[伝送線路|導波路]]などのようないくつかの現実の電磁界問題は、使用される物質の不規則な形状が複雑に関連しており、分析的に手計算することが不可能である。そのため、媒体と境界状態の様々なマクスウェルの方程式の補助方程式は、コンピュータにより計算される。計算電磁気学は、アンテナ・レーダー・[[衛星]]などの通信システム、{{仮リンク|ナノフォトニック|en|Nanophotonics}}デバイス、高速シリコンエレクトロニクス、[[医療画像処理]]、[[携帯電話]]のアンテナなど様々なアプリケーションの設計とモデリングにおける重要な分野である。
| |
| − | | |
| − | 計算電磁気学 (CEM) は、通常、問題の対象物について解く場合、対象物の領域の周りにある[[電界]]と[[磁界]]をコンピュータで計算することによって計算し問題を解く（言い換えると、計算電磁気学 (CEM) で任意形状の[[アンテナ]]構造について解析するには、[[アンテナ]]の[[電磁波]]放射パターンについて計算する）。
| |
| − | | |
| − | [[電磁波]]の電力進行方向（[[ポインティング・ベクトル]]）、[[導波路]]の正常モード、媒体による電磁波の[[分散 (光学)|分散]]、そして、[[散乱]]などは、電磁界解析において重要な項目であり、これらは全て[[電界]]と[[磁界]]の関係から[[コンピュータ]]を用いて計算することができる。計算電磁気学のモデルは、[[対称性]]を仮定したり、現実の構造を、筒状や球状あるいはその他の規則正しい幾何学的な形状に簡素化する。計算電磁気学のモデルは、徹底的に対称性を利用することによって、解析対象のシステムの次元を3次元から、2次元、更には 1次元にまで減少させることがある。計算電磁気学は、前述の幾つかの重要事項に挙げたような様々な問題を定式化できる。計算電磁気学の[[固有値]]定式化を用いることによって、特定の構造について、定常状態の正常モードを計算することができる。
| |
| − | | |
| − | [[過渡応答]]やインパルスによる電磁界への影響<!--"impulse field effects"-->は、[[時間領域]]の電磁界解析手法である[[FDTD法]]によって、より正確にモデル化される。曲面を含む物体は、[[有限要素法]] (FEM) や、<!--完全境界近似法（PBA法）といった-->非直交のグリッドを使用することによって、より正確に扱うことができる。{{仮リンク|ビーム伝播法|en|Beam propagation method}} (BPM)<ref>[http://pic.ee.sophia.ac.jp/html-e/bpm.html ビーム伝播法]</ref> は、導波路での電力の流れを解析する。以上のように、いくつかの手法によりある媒体の[[電磁界]]や[[電力]]の関係をモデル化できる場合でも、解析を行う応用分野ごとに計算電磁気学で使われるモデルが決まってくる。
| |
| − | | |
| − | == 解析方法 ==
| |
| − | コンピュータを用いた電磁界解析では、一般的に（直交あるいは非直交の）グリッドと呼ばれる空間の離散化（分割）し、媒体をモデル化した後、各グリッドごとにマクスウェルの方程式を解く。当然、計算に使われる空間の離散化は、コンピュータのメモリを消費するため、グリッドの数が多くなるほど、方程式を解く時間が長くなる。
| |
| − | | |
| − | {{Main2|コンピュータの仕様|並列コンピューティング}}
| |
| − | 大規模な電磁界解析では、コンピュータで使用するメモリ量や、CPU時間などにより、計算の制限を設定する。2007年現在、計算電磁気学では、[[スーパーコンピュータ]]や、高性能[[コンピュータ・クラスター|クラスタ]]、[[ベクトルプロセッサ]]、そして、[[並列コンピューティング|並列コンピュータ]]上で電磁界解析を実行することを問題としている。
| |
| − | | |
| − | 電磁界解析で用いられる定式化の典型的なものを、次に挙げる。
| |
| − | * 全時間領域の各瞬間でマクスウェルの方程式を解く。
| |
| − | * [[有限要素法]] (FEM) でモデル化されているときは、基本式の係数をまとめた逆行列を解く。
| |
| − | * {{仮リンク|転送行列法|en|transfer matrix method}}を用いる場合は、行列の積を解く。
| |
| − | * [[境界要素法|モーメント法]] (MoM)を使う場合は、積分方程式を解く。
| |
| − | * {{仮リンク|分割ステップ法|en|split-step method}}や{{仮リンク|ビーム伝播法|en|Beam propagation method}} (BPM) により計算する場合は、[[高速フーリエ変換|FFT]]と[[高速フーリエ変換|逆FFT]] (IFFT) を解く。
| |
| − | | |
| − | == 解析方法の選択 ==
| |
| − | 電磁界解析を行わなければならなくなったときに、どの手法を選ぶかは十分に検討する必要がある。積分方程式による解析を行うのか、それとも、微分方程式による解析を行うのか。そして、どのような時にどのような理由で高周波回路の近似が使われているのか理解しておく必要がある。
| |
| − | | |
| − | === 双曲型偏微分方程式のマクスウェルの方程式 ===
| |
| − | [[双曲幾何学|双曲系]]の[[偏微分方程式]]としてマクスウェルの方程式を定式化できる。これは[[双曲幾何学|双曲型]][[偏微分方程式]]を[[数値解析]]するための強力な数学理論を提供する。
| |
| − | | |
| − | 今、[[電磁波]]はx-y平面を伝播し（[[ポインティング・ベクトル|電力進行方向]]がx-y平面）、[[磁界]]がz軸方向にあると仮定する。この場合、[[電界]]はx-y平面にある。この電磁波は、'''TE波''' (Transverse Electric wave) と呼ばれる。
| |
| − | | |
| − | 2次元平面で、かつ[[誘電分極]]がないという条件が与えられている場合、マクスウェルの方程式は以下のように定式化できる。
| |
| − | {{Indent|<math>\frac{\partial}{\partial t}\bar{u}
| |
| − | + A\frac{\partial}{\partial x}\bar{u}
| |
| − | + B\frac{\partial}{\partial y}\bar{u}
| |
| − | + C\bar{u} = g</math>}}
| |
| − | ここで、u、A、BおよびCは次のように定義される：
| |
| − | {{Indent|
| |
| − | <math>\bar{u}=\left(\begin{matrix} E_x \\ E_y \\ H_z \end{matrix}\right)</math><br />
| |
| − | <math>A=\left(\begin{matrix} 0 & 0 & 0
| |
| − | \\ 0 & 0 & \frac{1}{\epsilon}
| |
| − | \\ 0 & \frac{1}{\mu} & 0 \end{matrix}\right)</math><br />
| |
| − | <math>B=\left(\begin{matrix} 0 & 0 & \frac{-1}{\epsilon}
| |
| − | \\ 0 & 0 & 0
| |
| − | \\ \frac{-1}{\mu} & 0 & 0 \end{matrix}\right)</math><br />
| |
| − | <math>C=\left(\begin{matrix} \frac{\sigma}{\epsilon} & 0 & 0
| |
| − | \\ 0 & \frac{\sigma}{\epsilon} & 0
| |
| − | \\ 0 & 0 & 0 \end{matrix}\right)</math>
| |
| − | }}
| |
| − | | |
| − | ==積分方程式による解析==
| |
| − | ===離散双極子近似 (DDA)===
| |
| − | '''{{仮リンク|離散双極子近似|en|discrete dipole approximation}}'''(discrete dipole approximation; '''DDA''')は、任意形状の対象物による[[電磁界]]の[[散乱]]と[[吸光|吸収]]を柔軟に計算できる手法である。定式化は積分形の[[マクスウェル方程式]]に基づいている。離散双極子近似(DDA)は、[[分極]]可能な箇所の有限な配列によって連続な対象物を近似する。[[分極]]された箇所は、各部分の[[電場]]に対応する[[双極子モーメント]]を持つ。もちろん、[[双極子]]は[[電場]]を介してお互いに作用するため、離散双極子近似(DDA)はしばしば結合双極子近似とも呼ばれる。結果として、線形的な連立方程式は、一般的に[[共役勾配法]]による[[反復法_(数値計算)|反復法]]を用いて解く。離散化行列には対称性があるため(積分形の[[マクスウェル方程式]]は[[畳み込み]]の形となるため)、[[共役勾配法]]による[[反復法_(数値計算)|反復]]計算を行っている間、[[行列]]と[[ベクトル]]を掛けるのに[[高速フーリエ変換]]を使用することができる。
| |
| − | | |
| − | ===モーメント法(MoM)・境界要素法(BEM)===
| |
| − | [[モーメント法|'''モーメント法'''(method of moments; '''MoM''')]]すなわち[[境界要素法|'''境界要素法'''(boundary element method; '''BEM''')]]は、[[積分方程式]] (「境界積分形」とも呼ばれる) として定式化された線形[[偏微分方程式]]を解く[[数値解析]]の手法である。[[境界要素法|境界要素法(BEM)]]は、[[流体力学]]、[[音響学]]、[[電磁気学]]、[[破壊力学]]、および[[塑性|可塑性解析]]など、工学と科学の多くの分野で適用できる。
| |
| − | | |
| − | 1980年代以降、[[境界要素法|境界要素法(BEM)]]はますます一般的に用いられるようになっている。[[境界要素法|境界要素法(BEM)]]は、[[偏微分方程式]]で定義された空間の全ての値を計算するというよりは、むしろ境界値だけについて計算するため、体積に対する表面積の比が小さい場合は、メモリの量やCPU時間といったコンピュータ・リソースの問題においてかなり効率的である。概念的に、[[境界要素法|境界要素法(BEM)]]は、モデル化された表面上で「メッシュ」を構成し問題を解く。しかしながら、多くの問題において、[[境界要素法|境界要素法(BEM)]]は体積を離散化する手法 ([[有限要素法]]、[[差分法|有限差分法]]、[[有限体積法]]より、かなり非効率的となる。[[境界要素法]]による定式化は、通常、全てにデータの入った行列<!--fully populated matrics-->を生じさせる。これは、データ量と計算時間が、解析対象物の大きさの二乗に応じて増加する傾向があることを意味する。
| |
| − | [[境界要素法|境界要素法(BEM)]]とは対照的に、[[有限要素法|有限要素]]行列は、(各要素は局所的に接続されるだけのため) 行列に必要なデータ量は 通常 解析対象物の大きさに比例して大きくなる。これらの (二乗で増えるデータ量の) 問題を改善するために、圧縮のテクニック(例えば、{{仮リンク|多重極展開|en|Multipole expansion}}やアダプティブ・クロス近似 (Adaptive Cross Approximation; ACA)<ref>[http://www.kuee.kyoto-u.ac.jp/~shamada/2004abumon_ohp01.pdf 近似品質管理付きACAによる正則化最小二乗解の高速計算法 濱田昌司(京大)]</ref> / {{仮リンク|階層行列|en|Hierarchical matrix}}(hierarchical matrices)<ref>[http://www.springer.com/math/cse/book/978-3-540-77146-3 "Hierarchical Matrices", Bebendorf, Springer], ISBN 978-3-540-77146-3</ref>)を使用できる。ただし、圧縮を行う場合は、圧縮追加の複雑さによる労力と、解析する問題の性格や幾何学的形状に大きく影響される成功率を十分に検討する必要がある。
| |
| − | | |
| − | [[境界要素法|境界要素法(BEM)]]は[[グリーン関数]]によって計算できる問題に適用可能である。通常、[[グリーン関数]]で解く問題は、[[線型性|線型]]で[[:en:Homogeneity (physics)|均質]]である媒体の領域について扱う。これは境界要素を有効に適用できる問題の範囲と一般性に関して大幅な制限を与える。[[境界要素法|境界要素法(BEM)]]の定式化で、非線形の問題も扱うこともできる。この場合、一般的に体積の積分方程式を導くことになるが、問題を解こうとする前に体積を離散化する必要があり、[[境界要素法|境界要素法(BEM)]]で最もよく引用される（均一物質を境界だけで解くという）利点が無くなることになる。
| |
| − | | |
| − | ====ソフトウェア====
| |
| − | [[境界要素法|境界要素法(BEM)]]を使ったソフトウェアの例：
| |
| − | * [http://eesof.tm.agilent.com/products/momentum_main.html Agilent Technologies社], Momentum
| |
| − | * [http://www.vectorfields.com/ Vector Fields社] Concerto
| |
| − | * [http://www.cst.com/Content/Products/MWS/IntegralEquationSolver.aspx CST社], CTS MICROWAVE STUDIO, [[積分方程式]]解析ソフト
| |
| − | * [http://www.si-list.net/swindex2.html Numerical Electromagnetic Code (NEC)]
| |
| − | * [http://www.sonnetsoftware.co.jp/ Sonnet社], [[モーメント法]]解析ソフト
| |
| − | * {{仮リンク|FEKO|en|FEKO}}, ハイブリッド[[モーメント法]]解析ソフト
| |
| − | * [http://www.sigrity.com/ Sigrity], SpeedXP
| |
| − | * [http://www.emcos.com/html/emc_studio.html EMC Studio], 全波[[電磁界]]解析ソフト, TriD, EMCoS
| |
| − | * [http://www.simion.com/cpo/ Scientific Instrument Service社], CPO (Charged Particle Optics)
| |
| − | * [http://www.ids-spa.it/page.php?f=92&id_div=7 IDS S.p.A., ADF]
| |
| − | * [http://www.efieldsolutions.com/ Efield] ESI Group
| |
| − | | |
| − | ===高速多重極法 (FMM)===
| |
| − | '''{{仮リンク|高速多重極法|en|Fast Multipole Method}}'''(fast multipole method; '''FMM''')は、[[モーメント法|モーメント法(MoM)]]([[境界要素法]])や[[エバルトの方法]]<ref>[http://staff.aist.go.jp/y-komeiji/peach/free/methods3.pdf 力とポテンシャルの計算, 古明地勇人 (産技研)]</ref>のような手法の代わりに適用されるコンピュータ上での[[電磁界]]解析手法である。高速多重極法(FMM)は、正確なシミュレーション手法であり、[[モーメント法|モーメント法(MoM)]]より計算上効率的である。メモリとプロセッサ実行時間ともに、[[モーメント法|モーメント法(MoM)]]より必要量が大幅に減る。高速多重極法(FMM)は、{{仮リンク|多重極展開|en|multipole expansion}}の手法に基づいており、{{仮リンク|Leslie Greengard|label=Greengard|en|Leslie Greengard}}と{{仮リンク|Vladimir Rokhlin|label=Rokhlin|en|Vladimir_Rokhlin_(American_scientist)}}によって最初に導入された<ref>[http://cacs.usc.edu/education/cs653/Greengard-FMM-JCP87.pdf "A Fast Algorithm for Particle Simulations", L.Greengard and V.Rokhlin, Journal of Computational Physics 73, 315-348 (1987)]</ref><!--日本語訳時にこの文献を参考とした。原文は別の文献を参考にしている可能性がある。--Kstar, 2008.10.19. -->。高速多重極法(FMM)は、[[モーメント法|モーメント法(MoM)]]の高速化に利用できる。
| |
| − | | |
| − | ==== リンク ====
| |
| − | * [http://math.nyu.edu/faculty/greengar/ Leslie Greengard] - [[ニューヨーク大学]]{{仮リンク|Leslie Greengard|en|Leslie Greengard}}のページ。幾つかの[[オープンソース]]な高速多重極法(FMM)のコードにリンク等。
| |
| − | * [http://www.umiacs.umd.edu/~ramani/ Ramani Duraiswami] - [[メリーランド大学]]Ramani Duraiswamiのページ。[http://www.umiacs.umd.edu/~ramani/cmsc878R/fmmdemo2.jar JAVA Animation of the Fast Multipole Method]、FMMの講義、コマーシャルソフトウェア等。
| |
| − | * [https://sites.google.com/site/rioyokota/home Rio Yokota] - {{仮リンク|アブデュラ王立工科大学|en|King Abdullah University of Science and Technology|label=KAUST}}横田理央のページ。ExaFMM, gemsFMM, PetFMM他。
| |
| − | | |
| − | ====ソフトウェア====
| |
| − | {{仮リンク|高速多重極法|en|Fast Multipole Method}}(FMM)とマルチレベルFMM(MLFMM)を使用するソフトウェアの例：
| |
| − | * [http://sourceforge.net/projects/puma-em/ Puma-EM] [[モーメント法|モーメント法(MoM)]]/マルチレベル[[:en:Fast Multipole Method|FMM]]による高機能で並列処理の可能なオープンソースの[[電磁界]]解析ソフト
| |
| − | * [http://www.harperlangston.com/kifmm3d/documentation/index.html KIFMM3d]
| |
| − | * [http://www.fastfieldsolvers.com FastFieldSolvers]
| |
| − | * [https://bitbucket.org/exafmm/exafmm/overview ExaFMM]
| |
| − | * [https://code.google.com/p/gemsfmm/ GemsFMM]
| |
| − | * [https://bitbucket.org/petfmm/petfmm-dev/ PetFMM]
| |
| − | | |
| − | ===漸化的T-マトリックス法 (RTMA)===
| |
| − | <!--日本語訳時参考文献：<ref>[http://ci.nii.ac.jp/naid/110003303342/ 「誘電体スラブの中に埋め込まれた柱状物体の周期的配列による電磁波の散乱」 Kushta Taras and 安元 清俊, 電子情報通信学会技術研究報告 OPE Vol.99 No.544 pp. 85-90 (2000)]</ref>-->
| |
| − | | |
| − | ===部分要素等価回路 (PEEC)===
| |
| − | '''{{仮リンク|部分要素等価回路|en|partial element equivalent circuitn}}'''(partial element equivalent circuit; '''PEEC''')は、[[電磁界]]と[[電子回路|回路]]を結合した解析に適した、三次元の全波<!--full-wave-->モデル化手法である。[[モーメント法|モーメント法(MoM)]]と異なって、部分要素等価回路(PEEC)は、[[直流|DC]]から最高[[周波数]]まで、分割<!--meshing-->の設定している全ての[[スペクトル]]で解析が可能な手法である。部分要素等価回路(PEEC)法では、[[積分方程式]]は[[キルヒホッフの法則 (電気回路)|キルヒホッフの電圧則(KVL)]]とみなして、部分要素等価回路(PEEC)の基本的なセルに適用される。そして、このことは、三次元の幾何学形状において完全な回路動作を解決する方法となる。等価回路による定式化は、[[SPICE (ソフトウェア)|SPICE]]タイプの回路素子を容易に追加することを可能にする。さらに、この方法によるモデルと解析方法は、時間領域と周波数領域の両方に適用することができる。部分要素等価回路(PEEC)モデルから生じる回路方程式は、修正{{仮リンク|ループ解析|en|Mesh analysis}}法(MLA; modified loop analysis)<ref>[http://www.ysaitoh.k.hosei.ac.jp/labor/ysaitoh/papers/Thesis/MasterThesis_Anabuki_2005.pdf 「電子機器内部の電流分布推定に関する考察」 p.12-13, 穴吹幸彦(法政大), 2005]</ref>や{{仮リンク|修正節点解析法|en|Modified nodal analysi}}(MNA; modified nodal analysis)<ref>[http://www.realize-se.co.jp/book_detail.php/RLB100142.html「回路シミュレーション技術とMOSFETモデリング」 サイペック], 2003, ISBN 4-89808-040-5</ref>による定式化を使用することで容易に構成できる。部分要素等価回路(PEEC)は、[[直流|DC]]解析を提供すること以外に、この(MLAやMNAといった)種類の問題について、様々な種類の回路素子を適切な型の行列で直接的に含めることができるため、[[モーメント法|モーメント法(MoM)]]よりもいくつもの他の利点を持っている。最近では、非直交の幾何学形状を扱うことができるように部分要素等価回路(PEEC)法が拡張されている。このモデル拡大(従来の[[直交]]座標系での定式化と一致している)は、より一般的な四角形(2次元)や6面体(3次元)の要素に加えた幾何学の[[マンハッタン距離|マンハッタン表現]]を含む。これにより、未知の数を最小限に保つことができ、その結果、非直行な幾何学形状においてもコンピュータの時間を短縮することができる<ref>[http://www.sm.luth.se/~jekman/PEEC/Program/ Partial Element Equivalent Circuit (PEEC) homepage]</ref>。
| |
| − | | |
| − | ===適応積分法(AIM)===
| |
| − | <!--原文：Adaptive Integral Method (AIM)-->
| |
| − | <!-- <ref>[http://www.egr.msu.edu/~bshanker/Fall2003%20Classes/EE%20802/aim_bleszynski.pdf "AIM: Adaptive integral method for solving large-scale electromagnetic scattering and radiation problems", E Bleszynski, M Bleszynski and T Jaroszewicz, Radio Science Vol.31 No.5 p.1225-1251, Sep.-Oct. 1996]</ref> -->
| |
| − | | |
| − | ==微分方程式による解析==
| |
| − | ===時間領域差分法 (FDTD)===
| |
| − | '''[[時間領域差分法]]'''（Finite-difference time-domain method; '''[[FDTD法]]'''）はポピュラーな計算電気力学モデリング手法である。[[FDTD法]]は分かり易くて、ソフトウェアで実装するのが簡単になるように考慮されている。
| |
| − | | |
| − | [[FDTD法]]が[[時間領域]]法であるため、一度シミュレーションを実行することにより、広い範囲の[[周波数領域]]を解くことができる。[[FDTD法]]は、一般的なグリッドを使った[[微分方程式|微分型]]の時間領域[[数値解析]]モデル化手法に属する。（[[偏微分方程式|偏微分型]]の）[[マクスウェルの方程式]]は、[[差分法|中央差分方程式]]<!--central-difference equations-->に変換され、離散化された後、ソフトウェアに実装される。方程式は{{仮リンク|リープフロッグ法|en|leapfrog integration}}で解かれる： まず、設定により与えられた時間領域の最初の時刻について、[[電場]]が解かれる。そして、時間領域の次の時刻で、[[磁場]]が解かれる。その後は、この過程が何度も繰り返される。
| |
| − | | |
| − | 基本的な[[FDTD法|FDTD]]アルゴリズムは、1966年に [http://www.ieeeaps.org/ IEEE Transactions on Antennas and Propagation] で Kane Yee が発表した論文<ref>[http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=1138693&isnumber=25435 "Numerical Solution of Initial Boundary Value Problems Involving Maxwell's Equations in Isotropic Media"], Kane S. Yee, IEEE Transactions on Antennas and Propagation Vol.AP-14,No.3 MAY,1966</ref>に遡る。「Finite-difference time-domain」という名称とその略語「FDTD」は、[http://www.emcs.org/ IEEE Transactions on Electromagnetic Compatibility] で1980年に {{仮リンク|Allen Taflove|en|Allen Taflove}}によって発表されたの論文<ref>[http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=4091372&isnumber=4091362 "Application of the Finite-Difference Time-Domain Method to Sinusoidal Steady-State Electromagnetic-Penetration Problems"], [[:en:Allen Taflove|Allen Taflove]], IEEE Transactions on Electromagnetic Compatibility, vol.EMC-22, No.3, Aug. 1980</ref>で初めて使われた。1990年頃以降、[[FDTD法]]は、物質の形状と[[電磁波]]との相互作用を扱う[[科学]]や[[工学]]の様々な問題をモデル化するための主要な方法となった。現在の[[FDTD法|FDTD]]によるモデル化の応用は、[[直流|DC]]近辺（全{{仮リンク|地球電離層導波路|en|Earth–ionosphere waveguide}}を含む超低周波の[[地球物理学]]）から[[高周波]]（レーダーシグネチャ<ref>[http://www.science.aster.ersdac.or.jp/jp/glossary/jp/si/signature.html 用語集-シグネチャ], ASTER SCIENCE PROJECT</ref>技術、[[アンテナ]]、[[無線通信]]デバイス、デジタルの相互接続<!--digital interconnects-->、生物医学画像処理<!--biomedical imaging/treatment-->）、そして[[可視光]]（[[フォトニック結晶]]、ナノ[[プラズモン|プラズモニクス]]、[[ソリトン]]、{{仮リンク|バイオフォトニクス|en|biophotonics}}）に至る範囲に及ぶ。
| |
| − | | |
| − | 現在では、FDTD法を単色波に最適化したNS-FDTD法（非標準FDTD法）も存在し、これは元のFDTD法よりも飛躍的に高い精度で計算できる。[http://www.cavelab.cs.tsukuba.ac.jp/nsfdtd/jp/index.html]
| |
| − | | |
| − | ===多重解像度時間領域法 (MRTD)===
| |
| − | '''多重解像度時間領域法'''（multiresolution time-domain; '''MRTD'''）は、[[ウェーブレット]]解析を基本として、[[FDTD法]]を適応型に変換した手法。
| |
| − | | |
| − | ===有限要素法 (FEM)===
| |
| − | '''[[有限要素法]]'''（finite element method; '''FEM'''）では、まず適切な[[偏微分方程式|偏微分方程式(PDE)]]および[[境界条件]]から境界値問題を設定し、重み付き残差法により弱形式（積分汎関数）を導出する。次に、解析領域に有限個の計算点を定義し、その点を用いて解析領域を適切に分割する（線分、三角形、四面体などの形状関数を使う）。そして、ガラーキン法などで弱形式を離散化し、線形連立方程式を構築する。そして、最後にその連立方程式を解くことで、計算点における近似解を求める。
| |
| − | | |
| − | 有限要素法は、定常状態の問題を解くことを基本とする。離散化において、スカラー場を計算する場合はスカラー形状関数を、ベクトル場を計算する場合はベクトル形状関数を用いる。しかし、ベクトル場を計算する方が煩雑になるため、スカラー場の計算で代用できる問題では、スカラー場を計算する方がよい。
| |
| − | | |
| − | 偏微分方程式の数値計算は、解析に使う近似式を作り出すことが肝要である。有限要素法は要素分割の自由度が高いため、差分法（FDTD法）などと比べて精度が高くなる。電磁界解析を行う領域が複雑な領域であるか、要求される精度が全領域で変化する場合、有限要素法は偏微分方程式を解くための良い選択肢となる。
| |
| − | | |
| − | 静電場などを計算する場合、支配方程式にはポアソン方程式が用いられ、これを定式化すると最終的には連立1次方程式ができる。また、電磁波などを計算する場合、支配方程式にはヘルムホルツ方程式が用いられ、これを定式化すると最終的には1次の固有方程式ができる。
| |
| − | | |
| − | ===有限積分法 (FIM)===
| |
| − | '''有限積分法'''(finite integration technique; '''FIT''')
| |
| − | | |
| − | ===擬似スペクトル時間領域法 (PSTD) ===
| |
| − | '''擬似スペクトル時間領域法'''（pseudospectral time domain; '''PSTD'''）<!--This classをPSTDとして訳したが、適切な訳か不明。-Kstar, 2008.11.3. -->は、一定時間内に[[マクスウェルの方程式]]の計算を完了させる<!--marching-in-time=「一定時間内に計算を完了させる」と訳したが、適切な訳か不明。-Kstar, 2008.11.3.-->ための方法として、離散[[フーリエ変換]]か離散チェビシェフ変換<ref>[http://ci.nii.ac.jp/naid/110003246316/ 「離散チェビシェフ変換を用いた画像圧縮」熱田 幸一(東海大) 他]</ref>のどちらかを使用し、2-Dまたは3-Dの単位格子に配置される電場と磁場ベクトル成分の空間的な伝播<!--spatial derivatives-->を計算する。
| |
| − | PSTD法は、[[FDTD法]]と同等程度の無視できるほど少ない位相・速度・異方性の[[数値解析]]誤差しか発生しない。したがって、PSTD法は非常に大きなサイズの問題をモデル化することができる。
| |
| − | | |
| − | [[マクスウェルの方程式]]を解くためのPSTD法に関する最近の総合的な概要は、A. Taflove, S. C. Hagness 編「Computational Electrodynamics: The Finite-Difference Time-Domain Method」<ref>「Computational Electrodynamics: The Finite-Difference Time-Domain Method」 A. Taflove, S. C. Hagness 編 Artech House, 2005, ISBN 978-1580538329</ref>内の Q. Liu と G. Zhao 著による 17章「Advances in PSTD Techniques」を参照のこと。
| |
| − | | |
| − | ===擬似スペクトル空間領域法 (PSSD) ===
| |
| − | '''擬似スペクトル空間領域法'''（Pseudo Spectral Spatial Domain; '''PSSD'''）は、[[電磁波]]の伝播する空間的な方向を選択することによって、[[マクスウェルの方程式]]を解く。したがって[[電場]]と[[磁場]]は、時間の関数として扱われ、さらに(可能であれば)特定の空間的な横軸の次元の関数としても扱われる。一時的な伝播が[[高速フーリエ変換]]を用いて周波数領域で計算されるため、この手法は擬似スペクトルである。[[電場]]と[[磁場]]が時間の関数として扱われるため、この手法により、伝播媒体の任意の分散を最小限の労力ですばやく正確にモデル化できる。J.C.A. Tyrrell他, J.Mod.Opt. 52, 973 (2005) <ref>[http://www.informaworld.com/smpp/content~content=a713736916~db=all~order=page "Pseudospectral spatial-domain: a new method for nonlinear pulse propagation in the few-cycle regime with arbitrary dispersion"], J.C.A. Tyrrell他, Journal of Modern Optics, Volume 52, Issue 7 May 2005 , pages 973 - 986</ref>を参照。
| |
| − | | |
| − | ===伝送路行列 (TLM)===
| |
| − | （[[SPICE (ソフトウェア)|SPICE]]やHSPICEなどのような<!--ala="in the manner of"=〜のような-->）[[回路シミュレータ]]が直接解析できる[[集中定数]]素子の集合や、素子を用いた固有の回路網、あるいは、[[Sパラメータ]]行列を用いるなど、いくつかの手法で'''{{仮リンク|伝送路行列|en|Transmission_line_matrix_method}}'''（Transmission Line Matrix; '''TLM'''）を定式化できる。多くの電磁界解析ソフトウェアが[[FDTD]]で解析可能であるという傾向があるが、伝送路行列（TLM）は能力において[[FDTD]]と同じく非常に柔軟性のある解析手法である。
| |
| − | | |
| − | ==そのほかの方法==
| |
| − | | |
| − | ===固有モード展開 (EME)===
| |
| − | '''[[固有モード展開]]'''(EigenMode Expansion; '''EME''')
| |
| − | | |
| − | ===物理光学(PO)===
| |
| − | '''[[物理光学]]'''（physical optics; '''PO'''）は、[[光学]]・[[電気工学]]および[[応用物理学]]では、{{仮リンク|高周波近似|en|high frequency approximation}}（短[[波長]][[近似]]）という名称で一般的に使用されている。[[物理光学]]（PO）は、[[波動]]効果を無視する[[幾何光学]]と、厳密な[[理論]]を扱う全波の[[電磁気学]]の、中間的な手法である。「物理（physical）」という言葉は、[[幾何光学]]（geometric optics）や[[:en:Optics#Classical optics|光線光学]]（ray optics）より物理的であることを意味し、正確な物理の理論であるということを意味するわけではない。
| |
| − | | |
| − | 伝播した電磁場や拡散した電磁場を計算するために、表面の電磁場と、その表面を通過する電磁場を[[積分法]]を用いて見積もり、光線光学を使用して[[近似]]を行っている。[[物理光学]]（PO）での近似は、問題の詳細が[[摂動]]として扱われるという点で、[[ボルン近似]]に類似している。
| |
| − | <!--「ボルン近似」＝「ボルン-オッペンハイマー近似」? by KStar 2008.11.16-->
| |
| − | <!--
| |
| − | ===回折幾何理論（GTD）===
| |
| − | 回折幾何理論（Geometric theory of diffraction; GTD）
| |
| − | -->
| |
| − | <!--
| |
| − | ===回折物理理論（PTD）===
| |
| − | 回折物理理論（Physical theory of diffraction; PTD）
| |
| − | -->
| |
| − | ===回折均一理論（UTD）===
| |
| − | '''{{仮リンク|回折均一理論|en|uniform theory of diffraction}}'''（Uniform theory of diffraction; '''UTD'''）は、電気的に小さい不連続や、同一点で一次元以上の不連続について、[[電磁気学]]の[[散乱]]問題を解くための高周波解析手法である。
| |
| − | | |
| − | 回折均一理論は、{{仮リンク|準光学|en|quasi optical}}として{{仮リンク|近接|en|near field}}電磁場を近似し、それぞれの物体と光源の回折関係を示す回折係数を解決するために、[[光線回折]]を使用する。これらの係数は、回折点から遠方での各方向について電磁界の強さと[[位相]]を計算するために使用される。そして、これらの電磁界は、完全な解を得るために入射電磁界と反射電磁界に加えられる。
| |
| − | | |
| − | ==検証==
| |
| − | 検証は、多くの電磁界シミュレーションソフトウェアの大きな問題の1つである。ユーザは、現在開発されているシミュレーションによって、シミュレーションできる有効な範囲を理解して習得しなければならない。シミュレーションが正しいか間違っているかということでなく、「結果は現実からどれくらい外れているか?」を確認する必要がある。
| |
| − | | |
| − | この質問に答えるために、以下の3ステップで検証する：
| |
| − | * '''シミュレーション結果と分析公式での比較'''
| |
| − | * '''解析手法間での相互比較'''
| |
| − | * '''測定とのシミュレーション結果の比較'''
| |
| − | | |
| − | ==電磁場解析ソフトウェア==
| |
| − | 電磁界シミュレータ ([[w:en:electromagnetic field solver|electromagnetic field solver]]) とも言う。解析可能な構造の次元によって、2.5次元と3次元に分けられる。
| |
| − | | |
| − | 2.5次元のシミュレータでは、例えば多層配線基板のように、複数の層で構成され各層に2次元の図形が存在する構造を解析することができる。しかし、同一の層には異なる電気物性をもつ構造を配置することはできない（異なる2つの層を接続するスルーホールを除く）。
| |
| − | | |
| − | 以下に主な電磁場解析ソフトウェアを挙げる。
| |
| − | ===2.5次元===
| |
| − | * ADS Momentum ([http://www.home.agilent.com/agilent/product.jspx?ckey=1385440&id=1385440&nid=-34333.804583.00&pageMode=OV&lc=eng&cc=US Agilent]社)
| |
| − | <!--* em ([http://www.sonnetsoftware.co.jp/ Sonnet]社) ... [[モーメント法]] ([[MOM法]]) による解析。-->
| |
| − | * EMSight (AWR社)
| |
| − | * Planar EM (Ansoft社) ... [[モーメント法]] ([[MOM法]]) による解析。
| |
| − | * [[SONNET]]（[http://www.sonnetsoftware.co.jp/ ソネット技研]）（[http://www.sonnetsoftware.com/ Sonnet Software Inc.]） ... [[モーメント法]] ([[MOM法]]) による解析。
| |
| − | | |
| − | ===3次元===
| |
| − | * [[EMS]] (EMWorks社) ... [[有限要素法]] ([[FEM法]]) による解析。[[SolidWorks]]完全アドイン。
| |
| − | * [[EMPro]] ([http://www.home.agilent.com/agilent/product.jspx?nid=-34278.0.00&cc=JP&lc=jpn Agilent]社) ... [[有限要素法]] ([[FEM法]]) および [[時間領域差分法]] ([[FDTD法]]) による解析。
| |
| − | * [[HFSS]] (Ansoft社) ... [[有限要素法]] ([[FEM法]]) による解析。
| |
| − | * [[MW Studio]] (CST社) ... [[有限積分法]] ([[FIT法]]) および [[完全境界近似法]] ([[PBA法]]) による解析。
| |
| − | * [[Microstrips]] (CST社 開発元Flomerics社の電磁波部門は2008年CST社に吸収された) ... [[伝送経路法]] ([[TLM法]]) による解析。
| |
| − | * [[ANSYS]] (ANSYS社)
| |
| − | * [[MagNet6]] (Infolytica Corporation社)
| |
| − | * [[JMAG]] ([http://www.jmag-international.com/jp/ JSOL社])
| |
| − | * [[FDTD Solutions]] ([https://www.lumerical.com/jp Lumerical Solutions]社)
| |
| − | * [[FIDERITY]] (Zeland Software社)
| |
| − | * [[IE3D]] (Zeland Software社)
| |
| − | * [[KeyFDTD]] (科学技術研究所)
| |
| − | * [[XFDTD]] (Remcom社)
| |
| − | * [[Femtet]] ([http://www.muratasoftware.com ムラタソフトウェア]社)...[[有限要素法]] ([[FEM法]]) による解析。
| |
| − | * [[EMSolution]]（サイエンスソリューションズ社）
| |
| − | * [[PHOTO-Series]]（[http://www.photon-cae.co.jp/ フォトン]）
| |
| − | * [[Poynting]]（[http://jp.fujitsu.com/solutions/hpc/app/poynting/ 富士通]) ... [[FDTD法]]による解析。
| |
| − | * [[WIPL-D]]（[http://www.wipl-d.jp/ WIPL-D(Japan), Inc.]）（[http://www.wipl-d.com/ WIPL-D d.o.o]） ... [[有限要素法]] ([[FEM法]]) による解析。
| |
| − | * [[Qm]] ([http://www.slock.co.jp シフトロック]) ... 磁気モーメント法による解析。
| |
| − | * [[CEM One]] ([https://www.esi-group.com/jp/software-solutions/virtual-environment/electromagnetics/cem-one-solution 日本イーエスアイ])
| |
| − | * [[ELF-Series]] ([http://www.elf.co.jp エルフ])
| |
| − | * [[PetaMagnetic]] ([http://petamagnetic.appspot.com ペタマグネティックソフトウェア]) ... FEM静電容量行列、FEMガルバニック腐食電流解析
| |
| − | | |
| − | == 脚注 ==
| |
| − | {{脚注ヘルプ}}
| |
| − | {{Reflist}}
| |
| − | <!-- == 参考文献 == {{Cite book}}、{{Cite journal}} -->
| |
| − | | |
| − | == 関連項目 ==
| |
| − | {{Commonscat|Computational electromagnetics}}
| |
| − | * [[数値解析]]
| |
| − | * [[CAE]]
| |
| − | * [[境界要素法]]
| |
| − | * [[有限要素法]]
| |
| − | * [[FDTD法]]（有限差分時間領域法）
| |
| − | * [[電子回路シミュレーション]]
| |
| − | * [[:en:Discrete dipole approximation codes]]
| |
| − | * [[:en:Codes for electromagnetic scattering by cylinders]]
| |
| − | * [[:en:Codes for electromagnetic scattering by spheres]]
| |
| − | | |
| − | == 外部リンク == <!-- {{Cite web}} -->
| |
| − | * [http://www.jsst.jp/j/ 日本シミュレーション学会]
| |
| − | * [http://www.jsap.or.jp/ （社）応用物理学会] [http://annex.jsap.or.jp/OSJ/ 日本光学会]
| |
| − | * [http://www.iee.or.jp/ （社）電気学会]
| |
| − | * [http://www.ieice.org/jpn/index.html （社）電子情報通信学会]
| |
| − | | |
| − | {{Physics-stub}}
| |
| − | {{Tech-stub}}
| |
| − | | |
| − | {{DEFAULTSORT:てんしはかいせき}}
| |
| − | [[category:計算物理学]]
| |
| − | [[category:物理光学]]
| |
| − | [[category:光学]]
| |
| − | [[category:電磁気学]]
| |
| − | [[category:電気回路]]
| |
| − | [[category:電気理論]]
| |