JP2015179407A - 電界計算装置、電界計算方法及びプログラム - Google Patents

Abstract

【課題】計算負荷を増大させずに、荷電粒子に作用する電界が精度よく計算される電界計算装置、電界計算方法及びプログラムを提供する。
【解決手段】解析対象領域を複数の要素に分割し、電界計算に必要な各要素の物理量と境界条件とを設定し、設定された物理量と境界条件とに基づいて各要素に含まれる荷電粒子に作用する第１電界を計算し、荷電粒子間の静電相互作用力に基づいて各荷電粒子に作用する第２電界を計算し、解析対象領域に存在する荷電粒子毎に第１電界と第２電界を合算して各荷電粒子に作用する電界を求める。
【選択図】図３

Description

本発明は、電界計算装置、電界計算方法及びプログラムに関する。

特許文献１には、電子写真方式を用いた画像形成装置において顕色剤に影響する電気的状態を計算する電気的状態計算装置であって、前記計算における条件の情報である条件情報を取得する条件情報取得部と、前記取得した条件情報に基づき、前記顕色剤が含まれる第１の誘電層と前記第１の誘電層に接する第２の誘電層と前記第１の誘電層とは反対側において前記第２の誘電層に接する金属体を含む系において前記顕色剤の電荷に応じて前記金属体内に複数誘起される鏡像電荷に基づいて前記電気的状態を計算する電気的状態計算部とを有し、前記電気的状態計算部は、前記顕色剤の電荷に応じて前記第２の誘電層と前記金属体との境界を鏡面として前記金属体内に誘起される鏡像電荷である第１の金属体鏡像電荷と、前記第１の鏡像電荷が前記第１の誘電層と前記第２の誘電層との境界において生ずる影響に応じて前記第２の誘電層と前記金属体との境界を鏡面として前記金属体内に誘起される鏡像電荷である第２の金属体鏡像電荷とに基づいて前記電気的状態を計算することを特徴とする、電気的状態計算装置が開示されている。

特許文献２には、粒子挙動解析の対象の領域である解析領域内の各粒子の挙動をそれぞれ計算する粒子挙動計算部と、前記粒子挙動計算部による各粒子の挙動の計算結果に基づいて、粒子の移動量が予め定められている閾値を超える粒子を判定する移動粒子判定部と、を備え、前記粒子挙動計算部は、計算段階ごとに、前記移動粒子判定部が前記閾値を超えると判定した粒子については精度の良い第１の挙動解析処理手法を適用して粒子の挙動をそれぞれ計算し、前記移動粒子判定部が前記閾値を超えると判定しなかった粒子については前記第１の挙動解析処理手法よりも精度が劣りかつ処理時間の短い第２の挙動解析処理手法を適用して粒子の挙動をそれぞれ計算する粒子挙動解析装置が開示されている。

特開２００９−２５１４６７号公報 特開２０１０−０７９４９３号公報

本発明の目的は、計算負荷を増大させずに、荷電粒子に作用する電界が精度よく計算される電界計算装置、電界計算方法及びプログラムを提供することにある。

請求項１の発明は、解析対象領域を複数の要素に分割する分割手段と、電界計算に必要な各要素の物理量と境界条件とを設定する設定手段と、設定された物理量と境界条件とに基づいて各要素に含まれる荷電粒子に作用する第１電界を計算する第１計算手段と、荷電粒子間の静電相互作用力に基づいて各荷電粒子に作用する第２電界を計算する第２計算手段と、前記解析対象領域に存在する荷電粒子毎に前記第１電界と前記第２電界を合算して各荷電粒子に作用する電界を求める合算手段と、を有する電界計算装置である。

請求項２の発明は、前記物理量が、要素内の荷電粒子の電荷量の総和である、請求項１に記載の電界計算装置。

請求項３の発明は、前記第２計算手段が、要素内に含まれる荷電粒子間の静電相互作用力に基づいて第２電界を計算する、請求項１又は請求項２に記載の電界計算装置。

請求項４の発明は、前記境界条件により電位固定境界に隣接する領域の要素サイズを、他の領域の要素サイズより小さくする、請求項１から請求項３までのいずれか１項に記載の電界計算装置。

請求項５の発明は、前記第１計算手段が、要素内に含まれる荷電粒子間の静電相互作用力に基づく電界が除かれるように、荷電粒子に作用する第１電界を計算する、請求項１から請求項４までのいずれか１項に記載の電界計算装置。

請求項６の発明は、前記第１計算手段が、前記複数の要素を順次注目要素とし、注目要素の電荷量を０として各要素に含まれる荷電粒子に作用する第３電界を繰り返し計算し、計算された第３電界を第１電界とする、請求項５に記載の電界計算装置。

請求項７の発明は、荷電粒子を２つ以上含む注目要素の電荷量を０とする、請求項６に記載の電界計算装置。

請求項８の発明は、前記分割手段は、各要素の電荷量が解析対象領域の総電荷量に応じて予め設定された数値範囲内となるように、解析対象領域を複数の要素に分割する、請求項１から請求項７までのいずれか１項に記載の電界計算装置。

請求項９の発明は、前記設定手段は、要素毎に初期電界を設定し、計算された第１電界と初期電界との差分に基づいて前記計算された第１電界を補正する補正手段を更に備えた、請求項１から請求項７までのいずれか１項に記載の電界計算装置。

請求項１０の発明は、解析対象領域を複数の要素に分割し、電界計算に必要な各要素の物理量と境界条件とを設定し、設定された物理量と境界条件とに基づいて各要素に含まれる荷電粒子に作用する第１電界を計算し、荷電粒子間の静電相互作用力に基づいて各荷電粒子に作用する第２電界を計算し、前記解析対象領域に存在する荷電粒子毎に前記第１電界と前記第２電界を合算して各荷電粒子に作用する電界を求める、電界計算方法。

請求項１１の発明は、コンピュータを、解析対象領域を複数の要素に分割する分割手段と、設定された電界計算に必要な各要素の物理量と境界条件とに基づいて各要素に含まれる荷電粒子に作用する第１電界を計算する第１計算手段と、荷電粒子間の静電相互作用力に基づいて各荷電粒子に作用する第２電界を計算する第２計算手段と、前記解析対象領域に存在する荷電粒子毎に前記第１電界と前記第２電界を合算して各荷電粒子に作用する電界を求める合算手段と、して機能させるプログラム。

請求項１、請求項１０、請求項１１に係る発明によれば、計算負荷を増大させずに、荷電粒子に作用する電界が精度よく計算される。

請求項２に係る発明によれば、各要素の電荷量が設定される。

請求項３に係る発明によれば、要素内の荷電粒子間の静電相互作用力が考慮される。

請求項４に係る発明によれば、本構成を備えない場合に比べて、電位固定境界に隣接する領域の電界が精度よく設定される。

請求項５、請求項６に係る発明によれば、要素内に含まれる荷電粒子が電界に与える影響が重複して考慮されなくなる。

請求項７に係る発明によれば、本構成を備えない場合に比べて、計算負荷が低減される。

請求項８に係る発明によれば、解析対象領域の総電荷量に応じて要素サイズが決められる。

請求項９に係る発明によれば、計算負荷を増加させずに、要素内に含まれる荷電粒子が電界に与える影響が重複して考慮されなくなる。

電界計算装置のハードウェア構成の一例を示すブロック図である。 電界計算部の機能構成の一例を示すブロック図である。 「電界計算処理」の手順の一例を示すフローチャートである。 解析対象領域が領域分割された一例を示す模式図である。 領域分割された解析対象領域に荷電粒子を配置した例を示す模式図である。 要素毎に設定された物理量（電荷量）と境界条件の一例を示す模式図である。 （Ａ）は有限要素法により計算された各要素の電界の一例を示す模式図である。（Ｂ）は要素内の荷電粒子間の静電相互作用の一例を示す図である。 物理量（電荷量）の設定方法の他の一例を示す模式図である。 （Ａ）〜（Ｃ）は電界計算値の補正方法の一例を示す模式図である。

＜電界計算装置＞
まず、電界計算装置について説明する。
図１は電界計算装置のハードウェア構成の一例を示すブロック図である。図１に示すように、電界計算装置１０は、電界計算部２０、情報入力部２２及び情報出力部２４を備えている。電界計算部２０は、コンピュータ等の計算装置として構成されている。電界計算部２０は、情報入力部２２から入力された情報に基づいて電界計算を行い、計算結果を情報出力部２４に出力する。

電界計算部２０は、ＣＰＵ（中央処理装置; Central Processing Unit）３０、各種プログラムを記憶したＲＯＭ（Read Only Memory）３２、プログラムの実行時にワークエリアとして使用されるＲＡＭ（Random Access Memory）３４、各種情報を記憶する不揮発性メモリ３６及び入出力インターフェース（Ｉ／Ｏ）３８を備えている。ＣＰＵ３０、ＲＯＭ３２、ＲＡＭ３４、不揮発性メモリ３６及びＩ／Ｏ３８の各々は、バス４０を介して互いに接続されている。情報入力部２２及び情報出力部２４の各々は、電界計算部２０のＩ／Ｏ３８に接続されている。

図２は電界計算部２０の機能構成の一例を示すブロック図である。図２に示すように、電界計算部２０は、情報受付部４２、マッピング部４４、第１電界計算部４６、第２電界計算部４８及び合算部５０を備えている。

情報受付部４２は、情報入力部２２から入力された情報を受け付ける。情報受付部４２は、粒子情報、境界情報及び領域分割情報を受け付ける。マッピング部４４は、情報受付部４２で受け付けた情報に基づいて、解析対象領域を複数の要素に分割し、分割された解析対象領域に解析対象の荷電粒子を配置して、複数の要素の各々について各要素の物理量を設定する。

ここで、粒子情報とは、解析対象となる荷電粒子の情報である。荷電粒子の情報には、荷電粒子を識別するための粒子ＩＤ、位置情報（座標）、及び物理量（半径、質量、電荷量、誘電率）が含まれる。境界情報とは、有限要素法で境界条件を設定するための情報である。境界条件は、電極や解析対象領域と外部領域との境界（電位固定境界）における電位などである。領域分割情報とは、解析対象領域を複数の要素に分割するための情報である。領域分割情報には、分割個数、各要素のサイズが含まれる。

第１電界計算部４６は、マッピング部４４で各要素に設定された物理量と境界条件とに基づいて、有限要素法により各要素に含まれる荷電粒子に作用する第１電界を計算する。第２電界計算部４８は、荷電粒子間の静電相互作用力に基づいて各荷電粒子に作用する第２電界を計算する。合算部５０は、第１電界計算部４６で算出された第１電界と、第２電界計算部４８で算出された第２電界とを、荷電粒子毎に合算する。そして、合算して得られた各荷電粒子に作用する電界を、計算結果として情報出力部２４に出力する。

＜電界計算処理＞
次に、電界計算部２０で実行される「電界計算処理」について更に詳しく説明する。
図３は「電界計算処理」の手順の一例を示すフローチャートである。「電界計算処理」のプログラムは、ＲＯＭ３２等の記憶装置に記憶されている。ＣＰＵ３０は、「電界計算処理」のプログラムを読み出し、ＲＡＭ３４をワークエリアとしてプログラムを実行する。

まず、ステップ１００で、粒子情報、境界情報、及び領域分割情報を含む各種情報を取得する。次に、ステップ１０２で、領域分割情報に基づいて、解析対象領域を複数の要素（部分領域）に分割する（領域分割）。図４は解析対象領域が領域分割された一例を示す模式図である。以下、説明を簡単にするために二次元場で説明する。図４に示すように、解析対象領域６０は、第１電極６２と第２電極６４との間に配置された領域である。第１電極６２と第２電極６４との間に電圧が印加されることで、解析対象領域６０には電界Ｅが発生する。

分割個数、各要素のサイズ等の分割条件は、領域分割情報として与えられる。本実施の形態では、解析対象領域は、各要素が複数の荷電粒子を含むように領域分割される。各要素の形状は、三角形、四角形等の多角形とされる。複数の要素の大きさは、解析対象領域内で異なっていてもよい。例えば、電位固定境界に隣接する領域の要素サイズを、他の領域の要素径サイズより小さくしてもよい。また、各要素の電荷量が、設定された数値範囲内となるように領域分割してもよい。電荷量の数値範囲（即ち、上限値と下限値）は、例えば、解析対象領域の総電荷量に応じて設定される。

図示した例では、解析対象領域６０は、要素１から要素９までの９個の要素に等分割されている。ここでは、各要素の形状は正方形であり、正方形の一片の長さで要素サイズが規定される。空間電荷が存在しない場合の電界は、ラプラスの式で与えられる。各要素の電位、各要素に作用する電界は、有限要素法、差分法等の領域分割法により計算される。

例えば、有限要素法では、要素の電位を座標の関数で近似し、ラプラスの式がポテンシャル・エネルギー最小の原理と等価であるという前提の下で、解析対象領域全体のポテンシャル・エネルギーが最小となるように、各要素の電位、各要素に作用する電界を計算する。実際には、各要素の頂点（節点）の電位を座標の関数で近似し、節点の数と同数の連立一次方程式を境界条件の下で解くことで、各節点の電位、電界が求められる。

次に、ステップ１０４で、粒子情報に基づいて、分割された解析対象領域に荷電粒子を配置する。図５は領域分割された解析対象領域に荷電粒子を配置した例を示す模式図である。図５に示すように、解析対象領域６０には、粒子情報に含まれる位置座標に応じて、複数の荷電粒子が配置される。各荷電粒子には、粒子情報に含まれる電荷量が付与される。図示した例では、解析対象領域６０には、荷電粒子Ａから荷電粒子Ｇまでの７個の荷電粒子が配置され、各荷電粒子には電荷量Ｑ_ａから電荷量Ｑ_ｇが付与される。

次に、ステップ１０６で、複数の要素の各々について各要素内の総電荷量を計算し、要素内の総電荷量を「要素の電荷量」として設定する。ここで、要素内の総電荷量とは、要素内の荷電粒子の電荷量の総和である。図６は要素毎に設定された物理量（電荷量）と境界条件の一例を示す模式図である。図６に示すように、複数の要素の各々について「電荷量」が設定される。図示した例では、要素１から要素９までの各要素に、電荷量Ｑ_１から電荷量Ｑ_９が設定される。

例えば、要素１には電荷量Ｑ_ａの荷電粒子Ａだけが含まれるので、要素１の電荷量Ｑ_１は電荷量Ｑ_ａである。また、要素５には電荷量Ｑ_ｃの荷電粒子Ｃと電荷量Ｑ_ｄの荷電粒子Ｄとが含まれるので、要素５の電荷量Ｑ_５は電荷量（Ｑ_ｃ＋Ｑ_ｄ）である。また、要素９は荷電粒子を含まないので、要素９の電荷量Ｑ_９は０である。

また、図示した例では、境界条件として、第１電極６２上の節点６６の電位φがＶ_１に設定され、第２電極６４上の節点６６の電位φがＶ_２に設定され、解析対象領域６０と外部領域との境界上の節点６６の電位φが０に設定されている。

次に、ステップ１０８で、各要素に設定された電荷量と境界条件とに基づいて、有限要素法により各要素に含まれる荷電粒子に作用する第１電界を計算する。空間電荷が存在する場合の電界は、ポアソンの式で与えられる。各要素の電荷量や誘電率が与えられると、各要素の電位、各要素に作用する電界は、空間電荷が存在しない場合と同様に、有限要素法により計算される。本実施の形態では、誘電率は一定と仮定して「電荷量」のみ設定しているが、各要素の物理量として「誘電率」を設定してもよい。

図７（Ａ）は有限要素法により計算された各要素の電界の一例を示す模式図である。図６に示すように、要素１から要素９に電荷量Ｑ_１から電荷量Ｑ_９が付与され、境界条件が設定されると、図７（Ａ）に示すように、要素１から要素９までの各要素について、電界Ｅ_１から電界Ｅ_９が求められる。例えば、要素５に含まれる荷電粒子Ｃと荷電粒子Ｄの各々には、電界Ｅ_５が作用する。

荷電粒子の間には、静電相互作用力としてクーロン力が発生する。同じ符号の電荷を有する荷電粒子間には反発力が発生し、異なる符号の電荷を有する荷電粒子間には引力が発生する。ここで、要素５を注目要素６８とすると、注目要素６８に含まれる荷電粒子と、注目要素以外の要素に含まれる荷電粒子との間の静電相互作用力は、有限要素法で考慮されている。一方、注目要素６８に含まれる荷電粒子間の静電相互作用力は、有限要素法で考慮されていない。即ち、有限要素法だけでは、各荷電粒子に作用する電界は正確には求められない。

次に、ステップ１１０で、要素内の荷電粒子の静電相互作用力から、要素内の各荷電粒子に作用する第２電界を計算する。図７（Ｂ）は要素内の荷電粒子間の静電相互作用の一例を示す図である。図７（Ｂ）に示すように、例えば、注目要素６８に含まれる荷電粒子Ｃを注目粒子とすると、荷電粒子Ｃと荷電粒子Ｄとの間の静電相互作用力（クーロン力）から、注目粒子である荷電粒子Ｃに作用する第２電界が計算される。

次に、ステップ１１２で、ステップ１０８で算出された第１電界と、ステップ１１０で算出された第２電界とを合算して、ルーチンを終了する。これにより、注目要素に含まれる荷電粒子と注目要素以外の要素に含まれる荷電粒子との間の静電相互作用力に加え、注目要素に含まれる荷電粒子間の静電相互作用力が考慮されて、各荷電粒子に作用する電界が精度よく求められる。例えば、荷電粒子Ｃに作用する第１電界を「Ｅ_５」とし、荷電粒子Ｃに作用する第２電界を「α」とすると、荷電粒子Ｃには合計「Ｅ_５＋α」の電界が作用することになる。

＜物理量の設定＞
次に、物理量の設定方法の変形例について説明する。
上記の実施の形態では、注目要素にも要素内の総電荷量が「要素の電荷量」として設定される例について説明したが、この設定で有限要素法により第１電界を計算すると、注目要素内に含まれる荷電粒子が電界に与える影響が重複考慮されてしまう。注目要素内に含まれる荷電粒子間の静電相互作用力に基づく電界を除くために、注目要素の電荷量を「０」と設定としてもよい。

図８は物理量（電荷量）の設定方法の他の一例を示す模式図である。図８に示すように、複数の要素の各々について「要素の電荷量」を設定する際に、注目要素６８である要素５の電荷量Ｑ_５を「０」に設定する。

この場合、複数の要素を順次注目要素とし、注目要素の電荷量を「０」として有限要素法により各要素に含まれる荷電粒子に作用する第３電界を繰り返し計算し、計算された複数の第３電界の代表値を第１電界の値とする。

例えば、図６に示す例では、要素１の電荷量Ｑ_１を「０」に設定して各要素の第３電界（１）を計算し、要素２の電荷量Ｑ_２を「０」に設定して各要素の第３電界（２）を計算し、・・・、要素９の電荷量Ｑ_９を「０」に設定して各要素の第３電界（９）を計算する。そして、要素１から要素９までの各々について、計算された第３電界を各要素の各々の第１電界とする。

なお、注目要素が複数の荷電粒子を含む要素である場合にだけ、注目要素の電荷量を「０」と設定としてもよい。注目要素が荷電粒子を含まない場合は、要素内の総電荷量は元々「０」である。注目要素に含まれる荷電粒子が１つの場合は、荷電粒子間の静電相互作用力を考慮する必要がない。図６に示す例では、要素５と要素７の電荷量を「０」と設定して第３電界を２回計算し、計算された２種類の第３電界の平均値を第１電界の値とすればよい。電荷量を「０」に設定する注目要素数が減少し、計算負荷が軽減する。

＜電界計算値の補正＞
次に、電界計算値の補正方法について説明する。
上記の実施の形態では、注目要素内に含まれる荷電粒子間の静電相互作用力に基づく電界を除く方法として、注目要素の電荷量を「０」に設定する例について説明したが、有限要素法による電界計算後に注目要素の電界計算値を補正してもよい。

図９（Ａ）〜（Ｃ）は電界計算値の補正方法の一例を示す模式図である。図９（Ａ）に示すように、複数の荷電粒子を含む要素については、要素に作用する第１電界（初期値）を（Ｅ_ｘ、Ｅ_ｙ）と設定しておく。次に、図９（Ｂ）に示すように、要素内の総電荷量を「要素の電荷量」として設定して、有限要素法により第１電界（計算値）を計算する。このときの電界増加量を（２δＥ_ｘ、２δＥ_ｙ）とすると、第１電界（計算値）は（Ｅ_ｘ＋２δＥ_ｘ、Ｅ_ｙ＋２δＥ_ｙ）と表される。要素内の電荷量が「０」であれば、要素内で増加した電荷量（差分）の半分を第１電界（初期値）から減ずる補正を行えばよい。従って、図９（Ｃ）に示すように、第１電界（補正値）は（Ｅ_ｘ−δＥ_ｘ、Ｅ_ｙ−δＥ_ｙ）となる。

例えば、図６に示す例では、要素５と要素７の電荷量を「０」と設定して第３電界を２回計算する必要がある。これに対し、電界計算後に注目要素の電界計算値を補正する場合は、要素５と要素７に作用する第１電界（初期値）を各々（Ｅ_ｘ、Ｅ_ｙ）と設定して、第１電界の計算を行えば済むので、計算負荷が低減される。

なお、上記実施の形態で説明した電界計算装置、電界計算方法及びプログラムの構成は一例であり、本発明の主旨を逸脱しない範囲内においてその構成を変更してもよいことは言うまでもない。例えば、例示したフローチャートにおいて、一部のステップを省略してもよく、他のステップを追加してもよい。また、必要に応じて各ステップの順序を入れ替えてもよい。

１０ 電界計算装置
２０ 電界計算部
２２ 情報入力部
２４ 情報出力部
３０ ＣＰＵ
３２ ＲＯＭ
３４ ＲＡＭ
３６ 不揮発性メモリ
４０ バス
４２ 情報受付部
４４ マッピング部
４６ 第１電界計算部
４８ 第２電界計算部
５０ 合算部
６０ 解析対象領域
６２ 第１電極
６４ 第２電極
６６ 節点
６８ 注目要素
Ａ〜Ｇ 荷電粒子

Claims (11)

1. 解析対象領域を複数の要素に分割する分割手段と、
   電界計算に必要な各要素の物理量と境界条件とを設定する設定手段と、
   設定された物理量と境界条件とに基づいて各要素に含まれる荷電粒子に作用する第１電界を計算する第１計算手段と、
   荷電粒子間の静電相互作用力に基づいて各荷電粒子に作用する第２電界を計算する第２計算手段と、
   前記解析対象領域に存在する荷電粒子毎に前記第１電界と前記第２電界を合算して各荷電粒子に作用する電界を求める合算手段と、
   を有する電界計算装置。
2. 前記物理量が、要素内の荷電粒子の電荷量の総和である、請求項１に記載の電界計算装置。
3. 前記第２計算手段が、要素内に含まれる荷電粒子間の静電相互作用力に基づいて第２電界を計算する、請求項１又は請求項２に記載の電界計算装置。
4. 前記境界条件により電位固定境界に隣接する領域の要素サイズを、他の領域の要素サイズより小さくする、請求項１から請求項３までのいずれか１項に記載の電界計算装置。
5. 前記第１計算手段が、要素内に含まれる荷電粒子間の静電相互作用力に基づく電界が除かれるように、荷電粒子に作用する第１電界を計算する、請求項１から請求項４までのいずれか１項に記載の電界計算装置。
6. 前記第１計算手段が、前記複数の要素を順次注目要素とし、注目要素の電荷量を０として各要素に含まれる荷電粒子に作用する第３電界を繰り返し計算し、計算された第３電界を第１電界とする、請求項５に記載の電界計算装置。
7. 荷電粒子を２つ以上含む注目要素の電荷量を０とする、請求項６に記載の電界計算装置。
8. 前記分割手段は、各要素の電荷量が解析対象領域の総電荷量に応じて予め設定された数値範囲内となるように、解析対象領域を複数の要素に分割する、請求項１から請求項７までのいずれか１項に記載の電界計算装置。
9. 前記設定手段は、要素毎に初期電界を設定し、
   計算された第１電界と初期電界との差分に基づいて前記計算された第１電界を補正する補正手段を更に備えた、請求項１から請求項７までのいずれか１項に記載の電界計算装置。
10. 解析対象領域を複数の要素に分割し、
    電界計算に必要な各要素の物理量と境界条件とを設定し、
    設定された物理量と境界条件とに基づいて各要素に含まれる荷電粒子に作用する第１電界を計算し、
    荷電粒子間の静電相互作用力に基づいて各荷電粒子に作用する第２電界を計算し、
    前記解析対象領域に存在する荷電粒子毎に前記第１電界と前記第２電界を合算して各荷電粒子に作用する電界を求める、
    電界計算方法。
11. コンピュータを、
    解析対象領域を複数の要素に分割する分割手段と、
    設定された電界計算に必要な各要素の物理量と境界条件とに基づいて各要素に含まれる荷電粒子に作用する第１電界を計算する第１計算手段と、
    荷電粒子間の静電相互作用力に基づいて各荷電粒子に作用する第２電界を計算する第２計算手段と、
    前記解析対象領域に存在する荷電粒子毎に前記第１電界と前記第２電界を合算して各荷電粒子に作用する電界を求める合算手段と、
    して機能させるプログラム。

Images