JP2004086248A - 分散媒質の電磁界の見積もり方法 - Google Patents
分散媒質の電磁界の見積もり方法 Download PDFInfo
- Publication number
- JP2004086248A JP2004086248A JP2002242447A JP2002242447A JP2004086248A JP 2004086248 A JP2004086248 A JP 2004086248A JP 2002242447 A JP2002242447 A JP 2002242447A JP 2002242447 A JP2002242447 A JP 2002242447A JP 2004086248 A JP2004086248 A JP 2004086248A
- Authority
- JP
- Japan
- Prior art keywords
- electric field
- field vector
- vector
- storage device
- electromagnetic field
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Complex Calculations (AREA)
Abstract
【課題】ローレンツ媒体やデバイ媒体用の混成媒体についても、同じアルゴリズムで扱うことが可能であり、したがって同時の計算が可能となり、それぞれの媒体の相互作用も取りいれられるため、正確な見積りができる分散媒質の電磁界の見積もり方法を提供する。
【解決手段】対象とする領域を微小領域に分け、それぞれに微小領域について、電磁界のマックスウェル方程式に時間領域有限差分法（ＦＤＴＤ）を適用して、１）磁界ベクトルを、時間的に前の磁界ベクトルと、時間的に前の電界ベクトルとを用いて演算装置で計算し、２）電荷成分を、時間的に前の電界ベクトルと、時間的に前の電荷成分とを用いて演算装置で計算し、３）電界ベクトル、を磁界ベクトル、時間的に前の電界ベクトル、電荷成分、時間的に前の電荷成分とを用いて演算装置で計算するステップにより、逐次、時間にそって数値解法により電界を求める。
【選択図】図１
【解決手段】対象とする領域を微小領域に分け、それぞれに微小領域について、電磁界のマックスウェル方程式に時間領域有限差分法（ＦＤＴＤ）を適用して、１）磁界ベクトルを、時間的に前の磁界ベクトルと、時間的に前の電界ベクトルとを用いて演算装置で計算し、２）電荷成分を、時間的に前の電界ベクトルと、時間的に前の電荷成分とを用いて演算装置で計算し、３）電界ベクトル、を磁界ベクトル、時間的に前の電界ベクトル、電荷成分、時間的に前の電荷成分とを用いて演算装置で計算するステップにより、逐次、時間にそって数値解法により電界を求める。
【選択図】図１
Description
【０００１】
【発明の属する技術分野】
この発明は、予め与えられた境界条件から電磁界の分布あるいは電磁界の伝播を、入力装置と記憶装置と演算装置と出力装置とを備えた計算装置を用いて見積もることのできる分散媒質の電磁界の見積もり方法に関するものである。
【０００２】
【従来の技術】数値解析により電磁界を見積る方法としては、時間領域有限差分法（ＦＤＴＤ：Ｆｉｎｉｔｅ Ｄｉｆｆｅｒｅｎｃｅ Ｔｉｍｅ Ｄｏｍａｉｎ）がよく知られており、これは文献１（Ｒ．Ｊ． Ｌｕｅｂｂｅｒｓ， Ｆ． Ｈｕｎｓｂｅｒｇｅｒ， Ｋ． Ｓ． Ｋｕｎｚ， Ｒ． Ｂ． Ｓｔａｎｄｌｅｒ， ａｎｄ Ｍ． Ｓｃｈｎｅｉｄｅｒ， “Ａ ｆｒｅｑｕｅｎｃｙ−ｄｅｐｅｎｄｅｎｔ ｆｉｎｉｔｅ−ｄｉｆｆｅｒｅｎｃｅ ｔｉｍｅ−ｄｏｍａｉｎ ｆｏｒｍｕｌａｔｉｏｎ ｆｏｒ ｄｉｓｐｅｒｓｉｖｅｍａｔｅｒｉａｌｓ，” ＩＥＥＥ Ｔｒａｎｓ． Ｅｌｅｃｔｒｏｍａｇｎ． Ｃｏｍｐａｔ．，ｖｏｌ． ３２， ｐｐ． ２２２−２２７， Ａｕｇ． １９９０．）で報告されている。また、ＦＤＴＤを拡張して分散媒質を取り扱いできるような方法として、文献２（Ｍ．Ｄ． Ｂｕｉ， Ｓ． Ｓ． Ｓｔｕｃｈｌｙ， ａｎｄ Ｇ． Ｉ． Ｃｏｓｔａｃｈｅ， “Ｐｒｏｐａｇａｔｉｏｎ ｏｆ ｔｒａｎｓｉｅｎｔｓ ｉｎｄｉｓｐｅｒｓｉｖｅ ｄｉｅｌｅｃｔｒｉｃ ｍｅｄｉａ，” ＩＥＥＥ Ｔｒａｎｓ． Ｍｉｃｒｏｗａｖｅ Ｔｈｅｏｒｙ Ｔｅｃｈ．， ｖｏｌ． ３９， ｐｐ．１１６５−１１７２， Ｊｕｌｙ １９９１．）に記載された再帰的畳込（ＲＣ：Ｒｅｃｕｒｓｉｖｅ Ｃｏｎｖｏｌｕｔｉｏｎ）法や、文献３（Ｔ． Ｋａｓｈｉｗａ， Ｎ．Ｙｏｓｈｉｄａ， ａｎｄ Ｉ． Ｆｕｋａｉ， “Ａ ｔｒｅａｔｍｅｎｔ ｂｙ ｔｈｅ ｆｉｎｉｔｅ−ｄｉｆｆｅｒｅｎｃｅ ｔｉｍｅ−ｄｏｍａｉｎｍｅｔｈｏｄ ｏｆ ｔｈｅ ｄｉｓｐｅｒｓｉｖｅ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ａｓｓｏｃｉａｔｅｄ ｗｉｔｈ ｏｒｉｅｎｔａｔｉｏｎｐｏｌａｒｉｚａｔｉｏｎ，” Ｔｒａｎｓ． ＩＥＩＣＥ， ｖｏｌ． Ｅ７３， ｎｏ． ８， ｐｐ． １３２６−１３２８， １９９０．）補助微分方程式（ＡＤＥ：ＡｕｘｉｌｉａｒｙＤｉｆｆｅｒｅｎｔｉａｌ Ｅｑｕａｔｉｏｎ）法、あるいは、文献４（Ｄ． Ｍ． Ｓｕｌｌｉｖａｎ， “Ｆｒｅｑｕｅｎｃｙ−ｄｅｐｅｎｄｅｎｔ ＦＤＴＤｍｅｔｈｏｄｓ ｕｓｉｎｇ Ｚ ｔｒａｎｓｆｏｒｍｓ，” ＩＥＥＥ Ｔｒａｎｓ． Ａｎｔｅｎｎａｓ Ｐｒｏｐａｇａｔ．， ｖｏｌ． ４０， ｐｐ．１２２３−１２３０， Ｏｃｔ． １９９２．）に記載されたＺ変換（ＺＴ：Ｚ−ｔｒａｎｓｆｏｒｍ）法等が知られている。これらは、入力装置と記憶装置と演算装置と出力装置とを備えたコンピュータに適用される。
【０００３】
この中で、ＲＣ法が用いるメモリ容量がもっとも少なくてすむ方法であるが、正確さにおいては、ＡＤＥ法やＺＴ法には及ばない。しかし、ＲＣ法でも、最近になって文献５（Ｄ．Ｆ． Ｋｅｌｌｅｙ ａｎｄ Ｒ． Ｊ． Ｌｕｅｂｂｅｒｓ， “Ｐｉｅｃｅｗｉｓｅ ｌｉｎｅａｒ ｒｅｃｕｒｓｉｖｅ ｃｏｎｖｏｌｕｔｉｏｎ ｆｏｒｄｉｓｐｅｒｓｉｖｅ ｍｅｄｉａ ｕｓｉｎｇ ＦＤＴＤ，” ＩＥＥＥ Ｔｒａｎｓ． Ａｎｔｅｎｎａｓ Ｐｒｏｐａｇａｔ．， ｖｏｌ． ４４， ｐｐ．７９２−７９７， Ｊｕｎｅ １９９６．）で報告されているように、部分的に線形の再帰的畳込（ＰＬＲＣ：Ｐｉｅｃｅｗｉｓｅ−Ｌｉｎｅａｒ ＲｅｃｕｒｓｉｖｅＣｏｎｖｏｌｕｔｉｏｎ）法を適用することにより、ＡＤＥ法やＺＴ法に匹敵する精度で法取り扱うことができるようになってきた。また、逆にＡＤＥ法やＺＴ法でもＲＣ法に匹敵する程度に必要なメモリ容量が少なくなってきている。他の特徴としては、ＡＤＥ法では、異なる分散媒質（ローレンツ媒体とデバイ媒体）を取り扱う時には、異なる取り扱いが必要であるが、ＲＣ法やＰＬＲＣ法では、同じ取り扱いで見積ることができる。また、数学的にはＡＤＥ法がもっとも簡素な取り扱いですむことが知られている。
【０００４】
ＡＤＥ法に関しては、多重極のローレンツ媒体やデバイ媒体用の３種類のアルゴリズムが文献６（Ｍ． Ｏｋｏｎｉｅｗｓｋｉ， Ｍ． Ｍｒｏｚｏｗｓｋｉ， ａｎｄＭ． Ａ． Ｓｔｕｃｈｌｙ， “Ｓｉｍｐｌｅ ｔｒｅａｔｍｅｎｔ ｏｆ ｍｕｌｔｉ−ｔｅｒｍ ｄｉｓｐｅｒｓｉｏｎ ｉｎ ＦＤＴＤ，” ＩＥＥＥＭｉｃｒｏｗａｖｅ Ｇｕｉｄｅｄ Ｗａｖｅ Ｌｅｔｔ．， ｖｏｌ． ７， ｐｐ． １２１−１２３， Ｍａｙ １９９７．）に報告されている。１）ＤＡＤＥＳ（ＤｅｂｙｅＡＤＥ Ｓｙｎｃｒｏｎｉｚｅｄ）法、２）ＬＡＤＥＳ（Ｌｏｒｅｎｔｚ ＡＤＥ Ｓｙｎｃｒｏｎｉｚｅｄ）法、３）ＬＡＤＥＰ（Ｌｏｒｅｎｔｚ ＡＤＥＰａｒｔｉａｌｌｙ Ｓｙｎｃｒｏｎｉｚｅｄ）法である。ＤＡＤＥＳ法においては、有限差分表現で半陰性の形式を用いて完全な時間同期が取られている。ＬＡＤＥＳ法においては、全時間での同期は行われているが、これには、さらに以前の電界成分を保存しておく必要がある。ＬＥＤＥＰ法においては、さらに以前の電界成分を保存しておく必要はないが、計算はＬＥＤＥＳ法より正確でない。これは、同期が部分的であるためである。
【０００５】
従来の数値解析により電磁界を見積る方法においては、電界ベクトルと磁界ベクトルと電流ベクトルとに対応した漸化式が用いられてきた。これに対し、本発明では、電界ベクトルと磁界ベクトルと電荷成分とに対応した漸化式を用いることにより、完全な同期をとりながら使用するメモリ領域を減少させることができ、また、デバイ媒体とローレンツ媒体とを、同じアルゴリズムで取り扱うことができるようになった。
【０００６】
【発明が解決しようとする課題】
従来の電磁界の計算方法では、ローレンツ媒体やデバイ媒体用の混成媒体については、同じアルゴリズムで扱うことはできなかったため、別々の計算が必要であり、正確な見積りができなかった。
【０００７】
この発明は上記に鑑み提案されたもので、ローレンツ媒体やデバイ媒体用の混成媒体についても、同じアルゴリズムで扱うことが可能であり、したがって同時の計算が可能となり、それぞれの媒体の相互作用も取りいれられるため、正確な見積りができる分散媒質の電磁界の見積もり方法を提供することを目的とする。
【０００８】
【課題を解決するための手段】
上記目的を達成するために、本発明における第１の発明は、予め与えられた境界条件から電磁界の分布あるいは電磁界の伝播を、入力装置と記憶装置と演算装置と出力装置とを備えた計算装置を用いて見積もる方法であって、見積もり対象とする領域を微小領域に分割し、その座標を上記の入力装置に入力し記憶装置に保存するステップと、各微小領域の誘電率と透磁率を設定し、それらの値を上記の入力装置に入力し記憶装置に保存するステップと、入射電磁界の波長を設定し、その値を上記の入力装置に入力し記憶装置に保存するステップと、入射電磁界を与える微小領域を設定し、それらの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存するステップと、を含む処理を行い、見積もり対象とする微小領域の、電界ベクトルと磁界ベクトルと電荷成分とを用いて、各微小領域について、電磁界のマックスウェル方程式に時間領域有限差分法（ＦＤＴＤ：ＦｉｎｉｔｅＤｉｆｆｅｒｅｎｃｅ Ｔｉｍｅ Ｄｏｍａｉｎ）を適用して、１）磁界ベクトルを、時間的に前の磁界ベクトルと、時間的に前の電界ベクトルとを用いて演算装置で計算し、２）電荷成分を、時間的に前の電界ベクトルと、時間的に前の電荷成分とを用いて演算装置で計算し、３）電界ベクトル、を磁界ベクトル、時間的に前の電界ベクトル、電荷成分、時間的に前の電荷成分とを用いて演算装置で計算するステップにより逐次、電界を時間にそって数値解法により求め、最後の結果を出力装置で出力することを特徴としている。
【０００９】
また、本発明における第２の発明は、ローレンツ媒体の電磁界を見積る方法で、後に述べる数２３、数２４および数２５の漸化式に基づいて計算することを特徴としている。
【００１０】
また、本発明における第３の発明は、デバイ媒体の電磁界を見積る方法で、後に述べる数２３、数２５および数３５の漸化式に基づいて計算することを特徴としている。
【００１１】
また、本発明における第４の発明は、デバイ媒体の電磁界を見積る方法で、後に述べる数２５、数４０および数４１の漸化式に基づいて計算することを特徴としている。
【００１２】
また、本発明における第５の発明は、ローレンツ媒体とデバイ媒体とが混在する場合の電磁界を見積る方法で、ローレンツ媒体部分については、後に述べる数２３、数２４および数２５を適用し、デバイ媒体については、数２３、数２５および数３３、あるいは、数２５、数４０および数４１の漸化式を用いて計算することを特徴としている。
【００１３】
【発明の実施の形態】
以下にこの発明の原理を説明し、次にコンピュータに適用する場合のアルゴリズムについて説明する。
【００１４】
以下の様に、電荷を内部に持たない場合を想定する。この電磁界のマックスウェルの方程式は、次の様になる。
【００１５】
【数１７】
【００１６】
【数１８】
【００１７】
［ローレンツ媒体］
ここで、ローレンツ分散特性をもったローレンツ媒体については、次の様に計算する。ωｐをｐ番目の分極対の共鳴周波数、Ａｐをその振幅、δｐを減衰率、εｓを静電気的な誘電率、ε∞を周波数無限大での誘電率、また、Δε＝εｓ−ε∞とするとき、ｐ対の分極対があるローレンツ媒質の誘電率は周波数座標で次の様になることが知られている。
【００１８】
【数１９】
【００１９】
角周波数ωでは、
【数２０】
とするとき、周波数座標でのアンペアの法則は次の様になる。
【００２０】
【数２１】
ここで、
【００２１】
【数２２】
である。
【００２２】
ここで、数２１あるいは数２２は、逆フーリエ変換をとり、時間についての差分表現にすると、次のようになる。
【００２３】
【数２３】
【００２４】
【数２４】
本発明のＦＤＴＤ法では、数１７から磁場は次の差分漸化式になる。
【００２５】
【数２５】
【００２６】
これらの数２３から数２５を用いてベクトルＥｎからベクトルＥｎ＋１を求める手順を、図１を用いて説明する。
１）ベクトルＨｎ−１／２とベクトルＥｎからベクトルＨｎ＋１／２を求める。
２）ベクトルＱｐ ｎ−１とベクトルＱｐ ｎとベクトルＥｎからベクトルＱｐ ｎ＋１を求める。
３）ベクトルＨｎ＋１／２とベクトルＥｎとベクトルＱｐ ｎ＋１とベクトルＱｐ ｎからベクトルＥｎ＋１を求める。
【００２７】
これに対し、ＬＡＤＥＳ法による従来の方法では、数２１に対応する磁場の差分漸化式を次のようにしている。
【数２６】
また、数２２に対応する電流の差分漸化式は、
【００２８】
【数２７】
である。
完全に時間同期を取ったＦＤＴＤ法では、漸化式は、次のようになる。
【００２９】
【数２８】
電流に関しては、
【００３０】
【数２９】
ここで、
【００３１】
【数３０】
【００３２】
【数３１】
である。
【００３３】
これらの数２３から数２５を用いてベクトルＥｎからベクトルＥｎ＋１を求める手順を、図６を用いて説明する。
１）ベクトルＨｎ−１／２とベクトルＥｎからベクトルＨｎ＋１／２を求める。
２）ベクトルＨｎ−１／２とベクトルＥｎとベクトルＥｎ−１とベクトルＪｐ ｎ−１とベクトルＪｐ ｎとからベクトルＥｎ＋１を求める。
３）ベクトルＥｎ＋１とベクトルＥｎ−１とベクトルＪｐ ｎ−１とベクトルＪｐ ｎとからベクトルＪｐ ｎ＋１を求める。
【００３４】
この様に従来のＦＤＴＤ法では８種の変数を用いるが、本発明のＦＤＴＤ法の漸化式では７種の変数を用いることが分かる。
【００３５】
また、数２３と数２７が同等なものであることは次の様にして理解することができる。
数２１と数２６から、次の数３２が導かれる。
【００３６】
【数３２】
また、ＦＴＤＴ法では、次の関係がある。
【００３７】
【数３３】
この数３３と数２９を数２３に代入することにより数２８が得られる。
【００３８】
［デバイ媒体］
次にデバイ分散特性を示すデバイ媒体の場合について説明する。τｐをｐ番目の分極対の緩和時間とするとき、ｐ対の分極対があるデバイ媒質の誘電率は周波数座標で次の様になることが知られている。
【００３９】
【数３４】
ここで、本発明で用いる数２４に対応する漸化式は、次の様になる。
【００４０】
【数３５】
【００４１】
この式の他に、数２３と数２５を用いて、図２に示す様に処理を進める。
１）ベクトルＨｎ−１／２とベクトルＥｎからベクトルＨｎ＋１／２を求める。
２）ベクトルＱｐ ｎ−１とベクトルＱｐ ｎとベクトルＥｎからベクトルＱｐ ｎ＋１を求める。
３）ベクトルＨｎ＋１／２とベクトルＥｎとベクトルＱｐ ｎ＋１とベクトルＱｐ ｎからベクトルＥｎ＋１を求める。
【００４２】
これは、上記のＦＴＤＴの場合と同様であるため、ローレンツ媒体とデバイ媒体との演算処理を同時に進められることがわかる。このため、これらの媒体の混成物についての演算処理を容易に行うことができることも明らかである。
【００４３】
一方文献（）に記載された従来のＤＡＤＥＳ法では、数２８と数２９とに対応する漸化式は、次の様になる。
【００４４】
【数３６】
【００４５】
【数３７】
ここで、
【００４６】
【数３８】
【００４７】
【数３９】
である。
【００４８】
これらの数３６から数３９を用いて演算処理を進めるには、図７に示す次の手順に従う。
１）ベクトルＨｎ−１／２とベクトルＥｎからベクトルＨｎ＋１／２を求める。
２）ベクトルＨｎ＋１／２とベクトルＥｎとベクトルＪｐ ｎからベクトルＥｎ＋１を求める。
３）ベクトルＥｎとベクトルＥｎ＋１とベクトルＪｐ ｎからベクトルＪｐ ｎ＋１を求める。
【００４９】
ちなみに、数２３と数３６とが同等のものであることは、数３７を数３２に代入し、この結果を数２３に代入することによって確認することができる。
【００５０】
さらに本発明の方法として、電流ベクトルＪｐ ｎ換えて電荷ベクトルＱｐ ｎを用いることにする。この場合、ＦＴＤＴ法では、数３５は数４０となる。
【００５１】
【数４０】
前記の数４０を数２３に代入して次の漸化式を得る事ができる。
【００５２】
【数４１】
ここで、
【００５３】
【数４２】
【００５４】
【数４３】
である。
【００５５】
これらの数４０から数４３を用いて演算処理を進めるには、図３に示す次の手順に従う。
１）ベクトルＨｎ−１／２とベクトルＥｎからベクトルＨｎ＋１／２を求める。
２）ベクトルＨｎ＋１／２とベクトルＥｎとベクトルＱｐ ｎからベクトルＥｎ＋１を求める。
３）ベクトルＥｎとベクトルＥｎ＋１とベクトルＱｐ ｎからベクトルＱｐ ｎ＋１を求める。この手順は、電流ベクトルＪｐ ｎを電荷ベクトルＱｐ ｎと入れ換えただけのものであるので、この場合も、ローレンツ媒体とデバイ媒体との演算処理を同時に行なうことは容易である。
【００５６】
［アルゴリズム］
ここで想定するコンピュータは、入力装置と記憶装置と演算装置と出力装置とを備えた計算装置である。このフローチャートの例を図４、図５に示し、以下に説明する。
Ｆ１）見積もり対象とする領域を微小領域に分割し、その座標を上記の入力装置に入力し記憶装置に保存する。
Ｆ２）各微小領域の誘電率と透磁率を設定し、それらの値を上記の入力装置に入力し記憶装置に保存する。
Ｆ３）入射電磁界の波長を設定し、その値を上記の入力装置に入力し記憶装置に保存する。
Ｆ４）入射電磁界を与える微小領域を設定し、それらの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存する。
Ｆ５）各微小領域の電界Ｅ、磁界Ｈ、あるいは、演算の途中で用いられる中間変数Ｑ１からＱｍまでをリセットし、前記の各微小領域の電界Ｅ、磁界Ｈの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存する。
Ｆ６）前記の電界Ｅ、磁界Ｈ、あるいは、演算の途中で用いられる中間変数Ｑ１からＱｍまでを一時的に保存する変数をそれぞれａＥ、ａＨ、ａＱ１〜ａＱｍとし、中間変数Ｑ１からＱｍまでを一時的に保存する他の変数をｂＱ１〜ｂＱｍとするとき、これらの変数をリセットする。
Ｆ７）その他の変数の初期値を設定する。
【００５７】
Ｆ８）第１の繰り返し演算の指標を当初値に設定する。
Ｆ９）繰り返し演算の第１のフラグを設定し、計算対象とする微小領域を決定する。
Ｆ１０）それぞれの微小領域について、それが入射電界を与える領域かどうかの正否を判断する。
Ｆ１４）この判断が正の場合に、見積り対象となる微小領域の磁界ベクトルＨと電界ベクトルＥに入射電磁界の値を引き渡し、後述の第２のフラグの位置から処理を行う。
Ｆ１１）この判断が否の場合には、各微小領域について、ｐ＝１、２、・・・ｍ、とするとき、以下の処理を行なう。まず、ＨをａＨあるいはａＥから計算する。ここで、ａＨあるいはａＥは、時間的に前の磁界ベクトルあるいは電界ベクトルである。より具体的には、数２５を用いて演算する。
Ｆ１２）ＱｐをａＥ、ａＱｐあるいはｂＱｐから計算する。ここで、ａＱｐやｂＱｐは、時間的に前の電荷成分である。より具体的には、ローレンツ媒体については数２４、デバイ媒体については数３５あるいは数４０を用いて演算する。
Ｆ１３）ＥをＨ、ａＥ、Ｑ１からＱｍまで、ａＱ１〜ａＱｍから計算する。より具体的には、数２４、数３５あるいは数４０を用いて演算する。ローレンツ媒体については数２３、デバイ媒体については数４１を用いて演算する。
Ｆ１５）第２のフラグを設定し、ａＨをＨに保存する。
Ｆ１６）ｂＱｐにａＱｐを保存する。
Ｆ１７）ａＱｐにＱｐを保存する。
Ｆ１８）ａＥにＥを保存する。
Ｆ１９）全微小領域についての計算が終了したかどうかの正否を判断する。否ならば第１のフラグの位置から再び計算する。
Ｆ２０）正ならば、次の処理を行なう。微小領域におけるＨ、Ｅ、Ｑ１からＱｍまでを上記の記憶装置に保存する。
Ｆ２１）すべての計算が終了したかどうかの正否を判断する。
Ｆ２２）否ならば、繰り返し演算の指標を１段階進めて第１のフラグの位置から再び計算する。
Ｆ２３）正ならば、この結果を出力装置に出力する。
Ｆ２４）計算を終了する。
【００５８】
【発明の効果】
本発明では、電界ベクトルの漸化式と磁界ベクトルの漸化式と電荷成分の漸化式とを用いることにより、完全な同期をとりながら使用するメモリ領域を減少させることができ、また、デバイ媒体とローレンツ媒体とを、同じアルゴリズムで取り扱うことができるようになり、これらの媒体の混在した媒体については、これを分割することなく取り扱うことができる。
【図面の簡単な説明】
【図１】本発明においてローレンツ媒体を取り扱う場合の漸次求める変数の流れを示すフローチャートである。
【図２】本発明においてデバイ媒体を取り扱う場合の漸次求める変数の流れを示すフローチャートである。
【図３】本発明においてデバイ媒体を取り扱う第２の場合の漸次求める変数の流れを示すフローチャートである。
【図４】本発明のフローチャートである。
【図５】本発明のフローチャートである。
【図６】従来のローレンツ媒体を取り扱う場合の漸次求める変数の流れを示すフローチャートである。
【図７】従来のデバイ媒体を取り扱う場合の漸次求める変数の流れを示すフローチャートである。
【発明の属する技術分野】
この発明は、予め与えられた境界条件から電磁界の分布あるいは電磁界の伝播を、入力装置と記憶装置と演算装置と出力装置とを備えた計算装置を用いて見積もることのできる分散媒質の電磁界の見積もり方法に関するものである。
【０００２】
【従来の技術】数値解析により電磁界を見積る方法としては、時間領域有限差分法（ＦＤＴＤ：Ｆｉｎｉｔｅ Ｄｉｆｆｅｒｅｎｃｅ Ｔｉｍｅ Ｄｏｍａｉｎ）がよく知られており、これは文献１（Ｒ．Ｊ． Ｌｕｅｂｂｅｒｓ， Ｆ． Ｈｕｎｓｂｅｒｇｅｒ， Ｋ． Ｓ． Ｋｕｎｚ， Ｒ． Ｂ． Ｓｔａｎｄｌｅｒ， ａｎｄ Ｍ． Ｓｃｈｎｅｉｄｅｒ， “Ａ ｆｒｅｑｕｅｎｃｙ−ｄｅｐｅｎｄｅｎｔ ｆｉｎｉｔｅ−ｄｉｆｆｅｒｅｎｃｅ ｔｉｍｅ−ｄｏｍａｉｎ ｆｏｒｍｕｌａｔｉｏｎ ｆｏｒ ｄｉｓｐｅｒｓｉｖｅｍａｔｅｒｉａｌｓ，” ＩＥＥＥ Ｔｒａｎｓ． Ｅｌｅｃｔｒｏｍａｇｎ． Ｃｏｍｐａｔ．，ｖｏｌ． ３２， ｐｐ． ２２２−２２７， Ａｕｇ． １９９０．）で報告されている。また、ＦＤＴＤを拡張して分散媒質を取り扱いできるような方法として、文献２（Ｍ．Ｄ． Ｂｕｉ， Ｓ． Ｓ． Ｓｔｕｃｈｌｙ， ａｎｄ Ｇ． Ｉ． Ｃｏｓｔａｃｈｅ， “Ｐｒｏｐａｇａｔｉｏｎ ｏｆ ｔｒａｎｓｉｅｎｔｓ ｉｎｄｉｓｐｅｒｓｉｖｅ ｄｉｅｌｅｃｔｒｉｃ ｍｅｄｉａ，” ＩＥＥＥ Ｔｒａｎｓ． Ｍｉｃｒｏｗａｖｅ Ｔｈｅｏｒｙ Ｔｅｃｈ．， ｖｏｌ． ３９， ｐｐ．１１６５−１１７２， Ｊｕｌｙ １９９１．）に記載された再帰的畳込（ＲＣ：Ｒｅｃｕｒｓｉｖｅ Ｃｏｎｖｏｌｕｔｉｏｎ）法や、文献３（Ｔ． Ｋａｓｈｉｗａ， Ｎ．Ｙｏｓｈｉｄａ， ａｎｄ Ｉ． Ｆｕｋａｉ， “Ａ ｔｒｅａｔｍｅｎｔ ｂｙ ｔｈｅ ｆｉｎｉｔｅ−ｄｉｆｆｅｒｅｎｃｅ ｔｉｍｅ−ｄｏｍａｉｎｍｅｔｈｏｄ ｏｆ ｔｈｅ ｄｉｓｐｅｒｓｉｖｅ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ａｓｓｏｃｉａｔｅｄ ｗｉｔｈ ｏｒｉｅｎｔａｔｉｏｎｐｏｌａｒｉｚａｔｉｏｎ，” Ｔｒａｎｓ． ＩＥＩＣＥ， ｖｏｌ． Ｅ７３， ｎｏ． ８， ｐｐ． １３２６−１３２８， １９９０．）補助微分方程式（ＡＤＥ：ＡｕｘｉｌｉａｒｙＤｉｆｆｅｒｅｎｔｉａｌ Ｅｑｕａｔｉｏｎ）法、あるいは、文献４（Ｄ． Ｍ． Ｓｕｌｌｉｖａｎ， “Ｆｒｅｑｕｅｎｃｙ−ｄｅｐｅｎｄｅｎｔ ＦＤＴＤｍｅｔｈｏｄｓ ｕｓｉｎｇ Ｚ ｔｒａｎｓｆｏｒｍｓ，” ＩＥＥＥ Ｔｒａｎｓ． Ａｎｔｅｎｎａｓ Ｐｒｏｐａｇａｔ．， ｖｏｌ． ４０， ｐｐ．１２２３−１２３０， Ｏｃｔ． １９９２．）に記載されたＺ変換（ＺＴ：Ｚ−ｔｒａｎｓｆｏｒｍ）法等が知られている。これらは、入力装置と記憶装置と演算装置と出力装置とを備えたコンピュータに適用される。
【０００３】
この中で、ＲＣ法が用いるメモリ容量がもっとも少なくてすむ方法であるが、正確さにおいては、ＡＤＥ法やＺＴ法には及ばない。しかし、ＲＣ法でも、最近になって文献５（Ｄ．Ｆ． Ｋｅｌｌｅｙ ａｎｄ Ｒ． Ｊ． Ｌｕｅｂｂｅｒｓ， “Ｐｉｅｃｅｗｉｓｅ ｌｉｎｅａｒ ｒｅｃｕｒｓｉｖｅ ｃｏｎｖｏｌｕｔｉｏｎ ｆｏｒｄｉｓｐｅｒｓｉｖｅ ｍｅｄｉａ ｕｓｉｎｇ ＦＤＴＤ，” ＩＥＥＥ Ｔｒａｎｓ． Ａｎｔｅｎｎａｓ Ｐｒｏｐａｇａｔ．， ｖｏｌ． ４４， ｐｐ．７９２−７９７， Ｊｕｎｅ １９９６．）で報告されているように、部分的に線形の再帰的畳込（ＰＬＲＣ：Ｐｉｅｃｅｗｉｓｅ−Ｌｉｎｅａｒ ＲｅｃｕｒｓｉｖｅＣｏｎｖｏｌｕｔｉｏｎ）法を適用することにより、ＡＤＥ法やＺＴ法に匹敵する精度で法取り扱うことができるようになってきた。また、逆にＡＤＥ法やＺＴ法でもＲＣ法に匹敵する程度に必要なメモリ容量が少なくなってきている。他の特徴としては、ＡＤＥ法では、異なる分散媒質（ローレンツ媒体とデバイ媒体）を取り扱う時には、異なる取り扱いが必要であるが、ＲＣ法やＰＬＲＣ法では、同じ取り扱いで見積ることができる。また、数学的にはＡＤＥ法がもっとも簡素な取り扱いですむことが知られている。
【０００４】
ＡＤＥ法に関しては、多重極のローレンツ媒体やデバイ媒体用の３種類のアルゴリズムが文献６（Ｍ． Ｏｋｏｎｉｅｗｓｋｉ， Ｍ． Ｍｒｏｚｏｗｓｋｉ， ａｎｄＭ． Ａ． Ｓｔｕｃｈｌｙ， “Ｓｉｍｐｌｅ ｔｒｅａｔｍｅｎｔ ｏｆ ｍｕｌｔｉ−ｔｅｒｍ ｄｉｓｐｅｒｓｉｏｎ ｉｎ ＦＤＴＤ，” ＩＥＥＥＭｉｃｒｏｗａｖｅ Ｇｕｉｄｅｄ Ｗａｖｅ Ｌｅｔｔ．， ｖｏｌ． ７， ｐｐ． １２１−１２３， Ｍａｙ １９９７．）に報告されている。１）ＤＡＤＥＳ（ＤｅｂｙｅＡＤＥ Ｓｙｎｃｒｏｎｉｚｅｄ）法、２）ＬＡＤＥＳ（Ｌｏｒｅｎｔｚ ＡＤＥ Ｓｙｎｃｒｏｎｉｚｅｄ）法、３）ＬＡＤＥＰ（Ｌｏｒｅｎｔｚ ＡＤＥＰａｒｔｉａｌｌｙ Ｓｙｎｃｒｏｎｉｚｅｄ）法である。ＤＡＤＥＳ法においては、有限差分表現で半陰性の形式を用いて完全な時間同期が取られている。ＬＡＤＥＳ法においては、全時間での同期は行われているが、これには、さらに以前の電界成分を保存しておく必要がある。ＬＥＤＥＰ法においては、さらに以前の電界成分を保存しておく必要はないが、計算はＬＥＤＥＳ法より正確でない。これは、同期が部分的であるためである。
【０００５】
従来の数値解析により電磁界を見積る方法においては、電界ベクトルと磁界ベクトルと電流ベクトルとに対応した漸化式が用いられてきた。これに対し、本発明では、電界ベクトルと磁界ベクトルと電荷成分とに対応した漸化式を用いることにより、完全な同期をとりながら使用するメモリ領域を減少させることができ、また、デバイ媒体とローレンツ媒体とを、同じアルゴリズムで取り扱うことができるようになった。
【０００６】
【発明が解決しようとする課題】
従来の電磁界の計算方法では、ローレンツ媒体やデバイ媒体用の混成媒体については、同じアルゴリズムで扱うことはできなかったため、別々の計算が必要であり、正確な見積りができなかった。
【０００７】
この発明は上記に鑑み提案されたもので、ローレンツ媒体やデバイ媒体用の混成媒体についても、同じアルゴリズムで扱うことが可能であり、したがって同時の計算が可能となり、それぞれの媒体の相互作用も取りいれられるため、正確な見積りができる分散媒質の電磁界の見積もり方法を提供することを目的とする。
【０００８】
【課題を解決するための手段】
上記目的を達成するために、本発明における第１の発明は、予め与えられた境界条件から電磁界の分布あるいは電磁界の伝播を、入力装置と記憶装置と演算装置と出力装置とを備えた計算装置を用いて見積もる方法であって、見積もり対象とする領域を微小領域に分割し、その座標を上記の入力装置に入力し記憶装置に保存するステップと、各微小領域の誘電率と透磁率を設定し、それらの値を上記の入力装置に入力し記憶装置に保存するステップと、入射電磁界の波長を設定し、その値を上記の入力装置に入力し記憶装置に保存するステップと、入射電磁界を与える微小領域を設定し、それらの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存するステップと、を含む処理を行い、見積もり対象とする微小領域の、電界ベクトルと磁界ベクトルと電荷成分とを用いて、各微小領域について、電磁界のマックスウェル方程式に時間領域有限差分法（ＦＤＴＤ：ＦｉｎｉｔｅＤｉｆｆｅｒｅｎｃｅ Ｔｉｍｅ Ｄｏｍａｉｎ）を適用して、１）磁界ベクトルを、時間的に前の磁界ベクトルと、時間的に前の電界ベクトルとを用いて演算装置で計算し、２）電荷成分を、時間的に前の電界ベクトルと、時間的に前の電荷成分とを用いて演算装置で計算し、３）電界ベクトル、を磁界ベクトル、時間的に前の電界ベクトル、電荷成分、時間的に前の電荷成分とを用いて演算装置で計算するステップにより逐次、電界を時間にそって数値解法により求め、最後の結果を出力装置で出力することを特徴としている。
【０００９】
また、本発明における第２の発明は、ローレンツ媒体の電磁界を見積る方法で、後に述べる数２３、数２４および数２５の漸化式に基づいて計算することを特徴としている。
【００１０】
また、本発明における第３の発明は、デバイ媒体の電磁界を見積る方法で、後に述べる数２３、数２５および数３５の漸化式に基づいて計算することを特徴としている。
【００１１】
また、本発明における第４の発明は、デバイ媒体の電磁界を見積る方法で、後に述べる数２５、数４０および数４１の漸化式に基づいて計算することを特徴としている。
【００１２】
また、本発明における第５の発明は、ローレンツ媒体とデバイ媒体とが混在する場合の電磁界を見積る方法で、ローレンツ媒体部分については、後に述べる数２３、数２４および数２５を適用し、デバイ媒体については、数２３、数２５および数３３、あるいは、数２５、数４０および数４１の漸化式を用いて計算することを特徴としている。
【００１３】
【発明の実施の形態】
以下にこの発明の原理を説明し、次にコンピュータに適用する場合のアルゴリズムについて説明する。
【００１４】
以下の様に、電荷を内部に持たない場合を想定する。この電磁界のマックスウェルの方程式は、次の様になる。
【００１５】
【数１７】
【数１８】
［ローレンツ媒体］
ここで、ローレンツ分散特性をもったローレンツ媒体については、次の様に計算する。ωｐをｐ番目の分極対の共鳴周波数、Ａｐをその振幅、δｐを減衰率、εｓを静電気的な誘電率、ε∞を周波数無限大での誘電率、また、Δε＝εｓ−ε∞とするとき、ｐ対の分極対があるローレンツ媒質の誘電率は周波数座標で次の様になることが知られている。
【００１８】
【数１９】
角周波数ωでは、
【数２０】
【００２０】
【数２１】
【００２１】
【数２２】
【００２２】
ここで、数２１あるいは数２２は、逆フーリエ変換をとり、時間についての差分表現にすると、次のようになる。
【００２３】
【数２３】
【数２４】
【００２５】
【数２５】
これらの数２３から数２５を用いてベクトルＥｎからベクトルＥｎ＋１を求める手順を、図１を用いて説明する。
１）ベクトルＨｎ−１／２とベクトルＥｎからベクトルＨｎ＋１／２を求める。
２）ベクトルＱｐ ｎ−１とベクトルＱｐ ｎとベクトルＥｎからベクトルＱｐ ｎ＋１を求める。
３）ベクトルＨｎ＋１／２とベクトルＥｎとベクトルＱｐ ｎ＋１とベクトルＱｐ ｎからベクトルＥｎ＋１を求める。
【００２７】
これに対し、ＬＡＤＥＳ法による従来の方法では、数２１に対応する磁場の差分漸化式を次のようにしている。
【数２６】
【００２８】
【数２７】
完全に時間同期を取ったＦＤＴＤ法では、漸化式は、次のようになる。
【００２９】
【数２８】
【００３０】
【数２９】
【００３１】
【数３０】
【数３１】
【００３３】
これらの数２３から数２５を用いてベクトルＥｎからベクトルＥｎ＋１を求める手順を、図６を用いて説明する。
１）ベクトルＨｎ−１／２とベクトルＥｎからベクトルＨｎ＋１／２を求める。
２）ベクトルＨｎ−１／２とベクトルＥｎとベクトルＥｎ−１とベクトルＪｐ ｎ−１とベクトルＪｐ ｎとからベクトルＥｎ＋１を求める。
３）ベクトルＥｎ＋１とベクトルＥｎ−１とベクトルＪｐ ｎ−１とベクトルＪｐ ｎとからベクトルＪｐ ｎ＋１を求める。
【００３４】
この様に従来のＦＤＴＤ法では８種の変数を用いるが、本発明のＦＤＴＤ法の漸化式では７種の変数を用いることが分かる。
【００３５】
また、数２３と数２７が同等なものであることは次の様にして理解することができる。
数２１と数２６から、次の数３２が導かれる。
【００３６】
【数３２】
【００３７】
【数３３】
【００３８】
［デバイ媒体］
次にデバイ分散特性を示すデバイ媒体の場合について説明する。τｐをｐ番目の分極対の緩和時間とするとき、ｐ対の分極対があるデバイ媒質の誘電率は周波数座標で次の様になることが知られている。
【００３９】
【数３４】
【００４０】
【数３５】
この式の他に、数２３と数２５を用いて、図２に示す様に処理を進める。
１）ベクトルＨｎ−１／２とベクトルＥｎからベクトルＨｎ＋１／２を求める。
２）ベクトルＱｐ ｎ−１とベクトルＱｐ ｎとベクトルＥｎからベクトルＱｐ ｎ＋１を求める。
３）ベクトルＨｎ＋１／２とベクトルＥｎとベクトルＱｐ ｎ＋１とベクトルＱｐ ｎからベクトルＥｎ＋１を求める。
【００４２】
これは、上記のＦＴＤＴの場合と同様であるため、ローレンツ媒体とデバイ媒体との演算処理を同時に進められることがわかる。このため、これらの媒体の混成物についての演算処理を容易に行うことができることも明らかである。
【００４３】
一方文献（）に記載された従来のＤＡＤＥＳ法では、数２８と数２９とに対応する漸化式は、次の様になる。
【００４４】
【数３６】
【数３７】
【００４６】
【数３８】
【数３９】
【００４８】
これらの数３６から数３９を用いて演算処理を進めるには、図７に示す次の手順に従う。
１）ベクトルＨｎ−１／２とベクトルＥｎからベクトルＨｎ＋１／２を求める。
２）ベクトルＨｎ＋１／２とベクトルＥｎとベクトルＪｐ ｎからベクトルＥｎ＋１を求める。
３）ベクトルＥｎとベクトルＥｎ＋１とベクトルＪｐ ｎからベクトルＪｐ ｎ＋１を求める。
【００４９】
ちなみに、数２３と数３６とが同等のものであることは、数３７を数３２に代入し、この結果を数２３に代入することによって確認することができる。
【００５０】
さらに本発明の方法として、電流ベクトルＪｐ ｎ換えて電荷ベクトルＱｐ ｎを用いることにする。この場合、ＦＴＤＴ法では、数３５は数４０となる。
【００５１】
【数４０】
【００５２】
【数４１】
【００５３】
【数４２】
【数４３】
【００５５】
これらの数４０から数４３を用いて演算処理を進めるには、図３に示す次の手順に従う。
１）ベクトルＨｎ−１／２とベクトルＥｎからベクトルＨｎ＋１／２を求める。
２）ベクトルＨｎ＋１／２とベクトルＥｎとベクトルＱｐ ｎからベクトルＥｎ＋１を求める。
３）ベクトルＥｎとベクトルＥｎ＋１とベクトルＱｐ ｎからベクトルＱｐ ｎ＋１を求める。この手順は、電流ベクトルＪｐ ｎを電荷ベクトルＱｐ ｎと入れ換えただけのものであるので、この場合も、ローレンツ媒体とデバイ媒体との演算処理を同時に行なうことは容易である。
【００５６】
［アルゴリズム］
ここで想定するコンピュータは、入力装置と記憶装置と演算装置と出力装置とを備えた計算装置である。このフローチャートの例を図４、図５に示し、以下に説明する。
Ｆ１）見積もり対象とする領域を微小領域に分割し、その座標を上記の入力装置に入力し記憶装置に保存する。
Ｆ２）各微小領域の誘電率と透磁率を設定し、それらの値を上記の入力装置に入力し記憶装置に保存する。
Ｆ３）入射電磁界の波長を設定し、その値を上記の入力装置に入力し記憶装置に保存する。
Ｆ４）入射電磁界を与える微小領域を設定し、それらの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存する。
Ｆ５）各微小領域の電界Ｅ、磁界Ｈ、あるいは、演算の途中で用いられる中間変数Ｑ１からＱｍまでをリセットし、前記の各微小領域の電界Ｅ、磁界Ｈの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存する。
Ｆ６）前記の電界Ｅ、磁界Ｈ、あるいは、演算の途中で用いられる中間変数Ｑ１からＱｍまでを一時的に保存する変数をそれぞれａＥ、ａＨ、ａＱ１〜ａＱｍとし、中間変数Ｑ１からＱｍまでを一時的に保存する他の変数をｂＱ１〜ｂＱｍとするとき、これらの変数をリセットする。
Ｆ７）その他の変数の初期値を設定する。
【００５７】
Ｆ８）第１の繰り返し演算の指標を当初値に設定する。
Ｆ９）繰り返し演算の第１のフラグを設定し、計算対象とする微小領域を決定する。
Ｆ１０）それぞれの微小領域について、それが入射電界を与える領域かどうかの正否を判断する。
Ｆ１４）この判断が正の場合に、見積り対象となる微小領域の磁界ベクトルＨと電界ベクトルＥに入射電磁界の値を引き渡し、後述の第２のフラグの位置から処理を行う。
Ｆ１１）この判断が否の場合には、各微小領域について、ｐ＝１、２、・・・ｍ、とするとき、以下の処理を行なう。まず、ＨをａＨあるいはａＥから計算する。ここで、ａＨあるいはａＥは、時間的に前の磁界ベクトルあるいは電界ベクトルである。より具体的には、数２５を用いて演算する。
Ｆ１２）ＱｐをａＥ、ａＱｐあるいはｂＱｐから計算する。ここで、ａＱｐやｂＱｐは、時間的に前の電荷成分である。より具体的には、ローレンツ媒体については数２４、デバイ媒体については数３５あるいは数４０を用いて演算する。
Ｆ１３）ＥをＨ、ａＥ、Ｑ１からＱｍまで、ａＱ１〜ａＱｍから計算する。より具体的には、数２４、数３５あるいは数４０を用いて演算する。ローレンツ媒体については数２３、デバイ媒体については数４１を用いて演算する。
Ｆ１５）第２のフラグを設定し、ａＨをＨに保存する。
Ｆ１６）ｂＱｐにａＱｐを保存する。
Ｆ１７）ａＱｐにＱｐを保存する。
Ｆ１８）ａＥにＥを保存する。
Ｆ１９）全微小領域についての計算が終了したかどうかの正否を判断する。否ならば第１のフラグの位置から再び計算する。
Ｆ２０）正ならば、次の処理を行なう。微小領域におけるＨ、Ｅ、Ｑ１からＱｍまでを上記の記憶装置に保存する。
Ｆ２１）すべての計算が終了したかどうかの正否を判断する。
Ｆ２２）否ならば、繰り返し演算の指標を１段階進めて第１のフラグの位置から再び計算する。
Ｆ２３）正ならば、この結果を出力装置に出力する。
Ｆ２４）計算を終了する。
【００５８】
【発明の効果】
本発明では、電界ベクトルの漸化式と磁界ベクトルの漸化式と電荷成分の漸化式とを用いることにより、完全な同期をとりながら使用するメモリ領域を減少させることができ、また、デバイ媒体とローレンツ媒体とを、同じアルゴリズムで取り扱うことができるようになり、これらの媒体の混在した媒体については、これを分割することなく取り扱うことができる。
【図面の簡単な説明】
【図１】本発明においてローレンツ媒体を取り扱う場合の漸次求める変数の流れを示すフローチャートである。
【図２】本発明においてデバイ媒体を取り扱う場合の漸次求める変数の流れを示すフローチャートである。
【図３】本発明においてデバイ媒体を取り扱う第２の場合の漸次求める変数の流れを示すフローチャートである。
【図４】本発明のフローチャートである。
【図５】本発明のフローチャートである。
【図６】従来のローレンツ媒体を取り扱う場合の漸次求める変数の流れを示すフローチャートである。
【図７】従来のデバイ媒体を取り扱う場合の漸次求める変数の流れを示すフローチャートである。
Claims (5)
- 予め与えられた境界条件から電磁界の分布あるいは電磁界の伝播を、入力装置と記憶装置と演算装置と出力装置とを備えた計算装置を用いて見積もる方法であって、
見積もり対象とする領域を微小領域に分割し、その座標を上記の入力装置に入力し記憶装置に保存するステップと、
各微小領域の誘電率と透磁率を設定し、それらの値を上記の入力装置に入力し記憶装置に保存するステップと、
入射電磁界の波長を設定し、その値を上記の入力装置に入力し記憶装置に保存するステップと、
入射電磁界を与える微小領域を設定し、それらの座標と入射電磁界の振幅、あるいは、強度と位相とを上記の入力装置に入力し記憶装置に保存するステップと、を含む処理を行い、
見積もり対象とする微小領域の、電界ベクトルと磁界ベクトルと電荷成分とを用いて、
各微小領域について、電磁界のマックスウェル方程式に時間領域有限差分法（ＦＤＴＤ：Ｆｉｎｉｔｅ Ｄｉｆｆｅｒｅｎｃｅ Ｔｉｍｅ Ｄｏｍａｉｎ）を適用して、
１）磁界ベクトルを、時間的に前の磁界ベクトルと、時間的に前の電界ベクトルとを用いて演算装置で計算し、
２）電荷成分を、時間的に前の電界ベクトルと、時間的に前の電荷成分とを用いて演算装置で計算し、
３）電界ベクトル、を磁界ベクトル、時間的に前の電界ベクトル、電荷成分、時間的に前の電荷成分とを用いて演算装置で計算するステップにより、
逐次、電界を時間にそって数値解法により求め、最後の結果を出力装置で出力することを特徴とする分散媒質の電磁界の見積もり方法。 - 予め与えられた境界条件から電磁界の分布あるいは電磁界の伝播を、入力装置と記憶装置と演算装置と出力装置とを備えた計算装置を用いて見積もる方法であって、
見積もり対象とする領域を微小領域に分割し、その座標を上記の入力装置に入力し記憶装置に保存するステップと、
各微小領域の誘電率と透磁率を設定し、それらの値を上記の入力装置に入力し記憶装置に保存するステップと、
入射電磁界の波長を設定し、その値を上記の入力装置に入力し記憶装置に保存するステップと、
入射電磁界を与える微小領域を設定し、それらの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存するステップと、
分極対を、ｐ＝１、２、・・・ｍ、とするとき、各微小領域の電界Ｅ、磁界Ｈ、あるいは、演算の途中で用いられる中間変数Ｑ１からＱｍまでをリセットし、前記の各微小領域の電界Ｅ、磁界Ｈの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存するステップと、を含む処理を行い、
第１の繰り返し演算の指標を当初値に設定し、繰り返し演算の第１のフラグを設定するステップと、
計算対象とする微小領域を決定するステップと、
それぞれの微小領域について、それが入射電界を与える領域かどうかの正否を判断するステップと、
この判断が正の場合に、見積り対象となる微小領域の磁界ベクトルＨと電界ベクトルＥに入射電磁界の値を引き渡し、後述の第２のフラグの位置から処理を行うステップと、
この判断が否の場合には、
各微小領域で、分極対を、ｐ＝１、２、・・・ｍ、とするとき、ωｐをｐ番目の分極対の共鳴周波数、Ａｐをその振幅、δｐを減衰率、εｓを静電気的な誘電率、ε∞を周波数無限大での誘電率、また、Δε＝εｓ−ε∞とし、
１）磁界ベクトルＨを、電界ベクトルＥと磁界ベクトルＨを用いた次の時間についての漸化式；
２）電荷成分Ｑｐを、電界ベクトルＥと電荷成分Ｑｐを用いた次の時間についての漸化式；
３）電界ベクトルＥを、磁界ベクトルＨと電界ベクトルＥと電荷成分Ｑｐを用いた次の時間についての漸化式；
第２のフラグを設定するステップと、
全微小領域についての計算が終了したかどうかの正否を判断するステップと、
否ならば第１のフラグの位置から再び計算するステップと、
正ならば、微小領域における磁界ベクトルＨと電界ベクトルＥと電荷成分Ｑｐを上記の記憶装置に保存するステップと、
すべての計算が終了したかどうかの正否を判断するステップと、
否ならば、繰り返し演算の指標を１段階進めて第１のフラグの位置から再び計算するステップと、
正ならば、計算を終了させ、この結果を出力装置に出力させるステップと、を含む処理を行うことを特徴とする分散媒質の電磁界の見積もり方法。 - 予め与えられた境界条件から電磁界の分布あるいは電磁界の伝播を、入力装置と記憶装置と演算装置と出力装置とを備えた計算装置を用いて見積もる方法であって、
見積もり対象とする領域を微小領域に分割し、その座標を上記の入力装置に入力し記憶装置に保存するステップと、
各微小領域の誘電率と透磁率を設定し、それらの値を上記の入力装置に入力し記憶装置に保存するステップと、
入射電磁界の波長を設定し、その値を上記の入力装置に入力し記憶装置に保存するステップと、
入射電磁界を与える微小領域を設定し、それらの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存するステップと、
各微小領域で、分極対を、ｐ＝１、２、・・・ｍ、とするとき、各微小領域の電界Ｅ、磁界Ｈ、あるいは、演算の途中で用いられる中間変数Ｑ１からＱｍまでをリセットし、前記の各微小領域の電界Ｅ、磁界Ｈの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存するステップと、を含む処理を行い、
第１の繰り返し演算の指標を当初値に設定し、繰り返し演算の第１のフラグを設定するステップと、
計算対象とする微小領域を決定するステップと、
それぞれの微小領域について、それが入射電界を与える領域かどうかの正否を判断するステップと、
この判断が正の場合に、見積り対象となる微小領域の磁界ベクトルＨと電界ベクトルＥに入射電磁界の値を引き渡し、後述の第２のフラグの位置から処理を行うステップと、
この判断が否の場合には、
ωｐをｐ番目の分極対の共鳴周波数、Ａｐをその振幅、δｐを減衰率、εｓを静電気的な誘電率、ε∞を周波数無限大での誘電率、また、Δε＝εｓ−ε∞とし、各微小領域について、
１）磁界ベクトルＨを、電界ベクトルＥと磁界ベクトルＨを用いた次の時間についての漸化式；
２）電荷成分Ｑｐを、電界ベクトルＥと電荷成分Ｑｐを用いた次の時間についての漸化式；
３）電界ベクトルＥを、磁界ベクトルＨと電界ベクトルＥと電荷成分Ｑｐを用いた次の時間についての漸化式；
第２のフラグを設定するステップと、
全微小領域についての計算が終了したかどうかの正否を判断するステップと、
否ならば第１のフラグの位置から再び計算するステップと、
正ならば、微小領域における磁界ベクトルＨと電界ベクトルＥと電荷成分Ｑｐを上記の記憶装置に保存するステップと、
すべての計算が終了したかどうかの正否を判断するステップと、
否ならば、繰り返し演算の指標を１段階進めて第１のフラグの位置から再び計算するステップと、
正ならば、計算を終了させ、この結果を出力装置に出力させるステップと、を含む処理を行うことを特徴とする分散媒質の電磁界の見積もり方法。 - 予め与えられた境界条件から電磁界の分布あるいは電磁界の伝播を、入力装置と記憶装置と演算装置と出力装置とを備えた計算装置を用いて見積もる方法であって、
見積もり対象とする領域を微小領域に分割し、その座標を上記の入力装置に入力し記憶装置に保存するステップと、
各微小領域の誘電率と透磁率を設定し、それらの値を上記の入力装置に入力し記憶装置に保存するステップと、
入射電磁界の波長を設定し、その値を上記の入力装置に入力し記憶装置に保存するステップと、
入射電磁界を与える微小領域を設定し、それらの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存するステップと、
各微小領域で、分極対を、ｐ＝１、２、・・・ｍ、とするとき、各微小領域の電界Ｅ、磁界Ｈ、あるいは、演算の途中で用いられる中間変数Ｑ１からＱｍまでをリセットし、前記の各微小領域の電界Ｅ、磁界Ｈの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存するステップと、を含む処理を行い、
第１の繰り返し演算の指標を当初値に設定し、繰り返し演算の第１のフラグを設定するステップと、
計算対象とする微小領域を決定するステップと、
それぞれの微小領域について、それが入射電界を与える領域かどうかの正否を判断するステップと、
この判断が正の場合に、見積り対象となる微小領域の磁界ベクトルＨと電界ベクトルＥに入射電磁界の値を引き渡し、後述の第２のフラグの位置から処理を行うステップと、
この判断が否の場合には、
ωｐをｐ番目の分極対の共鳴周波数、Ａｐをその振幅、δｐを減衰率、εｓを静電気的な誘電率、ε∞を周波数無限大での誘電率、また、Δε＝εｓ−ε∞とするとき、
１）磁界ベクトルＨを、電界ベクトルＥと磁界ベクトルＨを用いた次の時間についての漸化式；
２）電荷成分Ｑｐを、電界ベクトルＥと電荷成分Ｑｐを用いた次の時間についての漸化式；
３）電界ベクトルＥを、磁界ベクトルＨと電界ベクトルＥと電荷成分Ｑｐを用いた次の時間についての漸化式；
第２のフラグを設定するステップと、
全微小領域についての計算が終了したかどうかの正否を判断するステップと、
否ならば第１のフラグの位置から再び計算するステップと、
正ならば、微小領域における磁界ベクトルＨと電界ベクトルＥと電荷成分Ｑｐを上記の記憶装置に保存するステップと、
すべての計算が終了したかどうかの正否を判断するステップと、
否ならば、繰り返し演算の指標を１段階進めて第１のフラグの位置から再び計算するステップと、
正ならば、計算を終了させ、この結果を出力装置に出力させるステップと、を含む処理を行うことを特徴とする分散媒質の電磁界の見積もり方法。 - 予め与えられた境界条件から電磁界の分布あるいは電磁界の伝播を、入力装置と記憶装置と演算装置と出力装置とを備えた計算装置を用いて見積もる方法であって、
見積もり対象とする領域を微小領域に分割し、その座標を上記の入力装置に入力し記憶装置に保存するステップと、
各微小領域の誘電率と透磁率を設定し、それらの値を上記の入力装置に入力し記憶装置に保存するステップと、
入射電磁界の波長を設定し、その値を上記の入力装置に入力し記憶装置に保存するステップと、
入射電磁界を与える微小領域を設定し、それらの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存するステップと、
各微小領域で、分極対を、ｐ＝１、２、・・・ｍ、とするとき、各微小領域の電界Ｅ、磁界Ｈ、あるいは、演算の途中で用いられる中間変数Ｑ１からＱｍまでをリセットし、前記の各微小領域の電界Ｅ、磁界Ｈの座標と入射電磁界の強度と位相とを上記の入力装置に入力し記憶装置に保存するステップと、を含む処理を行い、
第１の繰り返し演算の指標を当初値に設定し、繰り返し演算の第１のフラグを設定するステップと、
計算対象とする微小領域を決定するステップと、
それぞれの微小領域について、それが入射電界を与える領域かどうかの正否を判断するステップと、
この判断が正の場合に、見積り対象となる微小領域の磁界ベクトルＨと電界ベクトルＥに入射電磁界の値を引き渡し、後述の第２のフラグの位置から処理を行うステップと、
この判断が否の場合には、
ωｐをｐ番目の分極対の共鳴周波数、Ａｐをその振幅、δｐを減衰率、εｓを静電気的な誘電率、ε∞を周波数無限大での誘電率、また、Δε＝εｓ−ε∞とするとき、
１）磁界ベクトルＨを、電界ベクトルＥと磁界ベクトルＨを用いた次の時間についての漸化式；
ローレンツ媒体部分については、
２−１）電荷成分Ｑｐを、電界ベクトルＥと電荷成分Ｑｐを用いた次の時間についての漸化式；
３−１）電界ベクトルＥを、磁界ベクトルＨと電界ベクトルＥと電荷成分Ｑｐを用いた次の時間についての漸化式；
デバイ媒体部分については、
２−２）電荷成分Ｑｐを、電界ベクトルＥと電荷成分Ｑｐを用いた次の時間についての漸化式；
３−２）電界ベクトルＥを、磁界ベクトルＨと電界ベクトルＥと電荷成分Ｑｐを用いた次の時間についての漸化式；
あるいは、デバイ媒体部分について、
２−３）電荷成分Ｑｐを、電界ベクトルＥと電荷成分Ｑｐを用いた次の時間についての漸化式；
３−３）電界ベクトルＥを、磁界ベクトルＨと電界ベクトルＥと電荷成分Ｑｐを用いた次の時間についての漸化式；
第２のフラグを設定するステップと、
全微小領域についての計算が終了したかどうかの正否を判断するステップと、
否ならば第１のフラグの位置から再び計算するステップと、
正ならば、微小領域における磁界ベクトルＨと電界ベクトルＥと電荷成分Ｑｐを上記の記憶装置に保存するステップと、
すべての計算が終了したかどうかの正否を判断するステップと、
否ならば、繰り返し演算の指標を１段階進めて第１のフラグの位置から再び計算するステップと、
正ならば、計算を終了させ、この結果を出力装置に出力させるステップと、を含む処理を行うことを特徴とする分散媒質の電磁界の見積もり方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002242447A JP3829191B2 (ja) | 2002-08-22 | 2002-08-22 | 分散媒質の電磁界の見積もり方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002242447A JP3829191B2 (ja) | 2002-08-22 | 2002-08-22 | 分散媒質の電磁界の見積もり方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004086248A true JP2004086248A (ja) | 2004-03-18 |
| JP3829191B2 JP3829191B2 (ja) | 2006-10-04 |
Family
ID=32051528
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002242447A Expired - Lifetime JP3829191B2 (ja) | 2002-08-22 | 2002-08-22 | 分散媒質の電磁界の見積もり方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3829191B2 (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005274233A (ja) * | 2004-03-23 | 2005-10-06 | Fujitsu Ltd | 電磁界シミュレータ、電磁界解析装置、電磁界シミュレートプログラム、および電磁界解析プログラム |
| KR20160024633A (ko) * | 2014-08-26 | 2016-03-07 | 국방과학연구소 | 플라즈마 매질에서 전자기파 해석 방법 |
| CN107944113A (zh) * | 2017-11-17 | 2018-04-20 | 西北工业大学 | 一种计算三维高速平动目标电磁散射场的方法 |
| JP2018136753A (ja) * | 2017-02-22 | 2018-08-30 | 富士通株式会社 | 磁性材料シミュレーションプログラム、磁性材料シミュレーション方法および磁性材料シミュレーション装置 |
| CN111783339A (zh) * | 2020-06-30 | 2020-10-16 | 西安理工大学 | 电磁波在随机色散介质中传播的pce-fdtd方法 |
-
2002
- 2002-08-22 JP JP2002242447A patent/JP3829191B2/ja not_active Expired - Lifetime
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005274233A (ja) * | 2004-03-23 | 2005-10-06 | Fujitsu Ltd | 電磁界シミュレータ、電磁界解析装置、電磁界シミュレートプログラム、および電磁界解析プログラム |
| KR20160024633A (ko) * | 2014-08-26 | 2016-03-07 | 국방과학연구소 | 플라즈마 매질에서 전자기파 해석 방법 |
| KR101671189B1 (ko) | 2014-08-26 | 2016-11-01 | 국방과학연구소 | 플라즈마 매질에서 전자기파 해석 방법 |
| JP2018136753A (ja) * | 2017-02-22 | 2018-08-30 | 富士通株式会社 | 磁性材料シミュレーションプログラム、磁性材料シミュレーション方法および磁性材料シミュレーション装置 |
| CN107944113A (zh) * | 2017-11-17 | 2018-04-20 | 西北工业大学 | 一种计算三维高速平动目标电磁散射场的方法 |
| CN107944113B (zh) * | 2017-11-17 | 2021-01-05 | 西北工业大学 | 一种计算三维高速平动目标电磁散射场的方法 |
| CN111783339A (zh) * | 2020-06-30 | 2020-10-16 | 西安理工大学 | 电磁波在随机色散介质中传播的pce-fdtd方法 |
| CN111783339B (zh) * | 2020-06-30 | 2024-04-16 | 西安理工大学 | 电磁波在随机色散介质中传播的pce-fdtd方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3829191B2 (ja) | 2006-10-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fu et al. | Transient eddy-current driver pickup probe response due to a conductive plate | |
| Giani et al. | Adaptive finite element methods for computing band gaps in photonic crystals | |
| Rekanos et al. | Approximation of Grünwald–Letnikov fractional derivative for FDTD modeling of Cole–Cole media | |
| Pinchuk et al. | Error estimation for automatic adaptive finite element mesh generation | |
| Yu et al. | Implementation of material interface conditions in the radial point interpolation meshless method | |
| JP2004086248A (ja) | 分散媒質の電磁界の見積もり方法 | |
| Xu et al. | Simulation of nonlinear superconducting rf losses derived from characteristic topography of etched and electropolished niobium surfaces | |
| Sboui et al. | Global modelling of microwave active circuits by an efficient iterative procedure | |
| Fu et al. | A compact higher‐order ADI‐FDTD method | |
| Williams et al. | Causality and waveguide circuit theory | |
| Razmjoo et al. | Electromagnetic time domain modeling using an improved meshless method | |
| Rogier et al. | Convergence behavior and acceleration of the Berenger and leaky modes series composing the 2-D Green's function for the microstrip substrate | |
| Sakamoto et al. | Multi-Objective Motor Design Optimization with Physics-Assisted Neural Network Model | |
| Hassan Hosseini et al. | Comparison of high frequency detailed generator models for partial discharge localization | |
| Liu et al. | Shift operator method and Runge-Kutta exponential time differencing method for plasma | |
| Wilhelm et al. | Exploiting the T (x) function in fast hysteresis models for transient circuit simulations | |
| Ma et al. | Improved finite difference time-domain method for anisotropic magnetised plasma based on shift operator | |
| Naamen et al. | Anisotropy analysis of the 3D-radial point interpolation method in lossy media | |
| Honarbakhsh et al. | Scattering by a 2D Crack: The meshfree collocation approach | |
| Shams et al. | Unconditionally stable vector-based meshless method for transient electromagnetic analysis | |
| Becerra et al. | Stable simulation of nonlinearly loaded lossy transmission lines with time marching approach | |
| Huang et al. | A Fast Time Domain Electric Field Integral Equations for Predicting Transient Responses of Dielectric Objects | |
| Kalantari et al. | Optical cloak design exploiting efficient anisotropic adjoint sensitivity analysis | |
| Zhang et al. | Recent advance in neuro-transfer function-assisted yield-driven EM optimization | |
| Harada et al. | PML absorbing boundary condition for noncubic cell time-domain method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060221 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060424 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20060606 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 3829191 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |