"粒子の運動の計算における任意の二つの粒子が衝突してから離れるまでの運動を計算する方法であって、時刻tpにおいて衝突している二粒子の接触面法線方向以下、法線方向という。に働く力を、弾性係数がkであり、時刻tpにおける粒子のめり込み量δpに比例する弾性反発力-k×δpと、粘性係数がηの粘性抵抗力の和で表し、一方の粒子1の、時刻tpより刻み時間Δtだけ進んだ時刻tp1における法線方向速度V1p1の予測値E1p1及び他方の粒子2の時刻tpにおける法線方向速度V2pとの差に比例する力η×V2p-E1p1を、粒子1に働く法線方向粘性力とし、さらに、他方の粒子2に働く法線方向粘性力を、同様に時刻tp1における粒子2の法線方向速度の予測値E2p1と時刻tpにおける法線方向速度V1pとの差に比例する力η×V1p-E2p1で表すことを特徴とする、二粒子が衝突してから離れるまでの法線方向の運動の計算方法。ここで、時刻指標pは整数であり、時間がΔt進行するとpの値が1増加する。また、右上添え字pは、時刻指標pにおける当該物理指標を意味する。粒子1の質量をm1及び、粒子2の質量をm2として、時刻tp1における、粒子1の予測法線方向速度E1p1及び粒子2の予測法線方向速度E2p1をそれぞれ及びで表すことを特徴とする、請求項1に記載の二粒子が衝突してから離れるまでの法線方向の運動の計算方法。時刻tpにおいて、粒子1及び粒子2が衝突から受ける法線方向の力F1p及びF2pを及びで表すことを特徴とする、請求項1又は2に記載の二粒子が衝突してから離れるまでの法線方向の運動の計算方法。多粒子系粒子の挙動を解析するための粒子挙動解析方法であって、個々の粒子の運動を運動方程式を用いて離散要素法で計算し、任意の2個の粒子が衝突してから離れるまでの運動の計算について請求項1乃至3のいずれかに記載の二粒子が衝突してから離れるまでの法線方向の運動の計算方法を用いることを特徴とする粒子挙動解析方法。シャトルコンベヤからサージホッパーに投入され,ドラムフィーダから切り出され,フルイの上を通って落下した焼結原料の粒子配置分布を求めるにあたり,請求項4に記載の粒子挙動解析方法を用いることを特徴とする、焼結原料の粒子配置分布解析方法。本発明は、大規模な設備における多数の粒子の挙動を計算するに際して、精度を落とすことなく、計算を高速化する計算方法及び、そのような計算方法を用いることで短時間での粒子の挙動の解析が可能となる粒子挙動解析方法を提供することを目的とするものである。本発明者らは、本発明が対象とする多粒子系の粒子挙動解析の場合、粒子間の接触あるいは衝突が非常に多く、このような衝突過程に関する計算を省略できれば計算負荷を低下できると考えた。精度を保ったまま、時間刻みを大きくして衝突過程を省略できる計算方法について調査研究を行った。その結果、衝突した二粒子間の接触面法線方向の、時刻tpにおける粘性力を与える、二粒子間の相対速度において、一方の粒子の速度を時刻tp1の予測速度に置き換えることにより、目的が達成されることが分かった。本発明の要旨は、以下の通りである。1粒子の運動の計算における任意の二つの粒子が衝突してから離れるまでの運動を計算する方法であって、時刻tpにおいて衝突している二粒子の接触面法線方向以下、法線方向という。に働く力を、弾性係数がkであり、時刻tpにおける粒子のめり込み量δpに比例する弾性反発力-k×δpと、粘性係数がηの粘性抵抗力の和で表し、一方の粒子1の、時刻tpより刻み時間Δtだけ進んだ時刻tp1における法線方向速度V1p1の予測値E1p1及び他方の粒子2の時刻tpにおける法線方向速度V2pとの差に比例する力η×V2p-E1p1を、粒子1に働く法線方向粘性力とし、さらに、他方の粒子2に働く法線方向粘性力を、同様に時刻tp1における粒子2の法線方向速度の予測値E2p1と時刻tpにおける法線方向速度V1pとの差に比例する力η×V1p-E2p1で表すことを特徴とする、二粒子が衝突してから離れるまでの法線方向の運動の計算方法。ここで、時刻指標pは整数であり、時間がΔt進行するとpの値が1増加する。また、右上添え字pは、時刻指標pにおける当該物理指標を意味する。2粒子1の質量をm1及び、粒子2の質量をm2として、時刻tp1における、粒子1の予測法線方向速度E1p1及び粒子2の予測法線方向速度E2p1をそれぞれ及びで表すことを特徴とする、前記1記載の二粒子が衝突してから離れるまでの法線方向の運動の計算方法。3時刻tpにおいて、粒子1及び粒子2が衝突から受ける法線方向の力F1p及びF2pを及びで表すことを特徴とする、前記1又は2に記載の二粒子が衝突してから離れるまでの法線方向の運動の計算方法。4多粒子系粒子の挙動を解析するための粒子挙動解析方法であって、個々の粒子の運動を運動方程式を用いて離散要素法で計算し、任意の2個の粒子が衝突してから離れるまでの運動の計算について上記1乃至3のいずれかに記載の二粒子が衝突してから離れるまでの法線方向の運動の計算方法を用いることを特徴とする粒子挙動解析方法。5シャトルコンベヤからサージホッパーに投入され,ドラムフィーダから切り出され,フルイの上を通って落下した焼結原料の粒子配置分布を求めるにあたり,上記4に記載の粒子挙動解析方法を用いることを特徴とする、焼結原料の粒子配置分布解析方法。[第一の実施形態]粒子の衝突の判定とめりこみ量対象とする粒子の形状は特に限定しないが、ここでは球である場合を例にとって説明する。半径r1の球である粒子1と、半径r2の粒子2の二つの粒子が接触しているかどうかの判定は、種々の方法があるが、例えば粒子1の中心と粒子2の中心の距離がr1r2以下であれば接触していると判定してもよい。粒子が球でなければ、粒子の重心をその中心とし、平均半径をその半径としてもよい。めり込み量は、前記r1r2と接触している二つの粒子の中心間距離の差によって求められる。計算に用いる場合は、そのようにして得られた差を更に、前記r1r2で除した値を用いる場合もある。粒子1と粒子2相互間のめり込み量をδとし、二つの粒子の中心の、その中心を結ぶ線上で計った位置をx1、x2とすれば、二つの粒子の中心距離は、|x1-x2|と表される。これらは、時間的に変化する。各々の量の時刻tpにおける値であることを右肩添え字pで示せば、時刻tpにおけるめり込み量δpは次のように表すことができる。めり込みの例は図1に示されている。δp={r1r2-|x1p-x2p|}/r1r2[数式1]ある時点で二つの粒子の接触判定を行った場合、接触と判定されたものについて、二つの粒子の中心距離は各々の半径の和よりもかなり小さく、即ち大きいめりこみ量である場合が多い。このようなめり込みが大きい場合、反発力も大きく、時間刻みを細かくとらないと、大きな反発力によって粒子が急速に反発してしまい、計算が不安定になってしまう場合がある。接触面法線方向に働く力接触した二つの粒子が離れるまでの衝突過程は、その接触面に対する法線方向に働く力以下、接触法線力という。が大きな影響を与える。即ち、該接触面から該法線に沿って大きく後退すれば衝突過程は早く終了する。従って、該法線方向に働く力を調整することで衝突過程を省略できるのである。また、ここで時間の推移は、刻み時間Δtずつ推移するものとする。即ち、時刻tpにおいて、時刻指標p=0からpΔtだけ時間が経過している。非特許文献1においては、時刻tpにおいて接触している二つの粒子に、その接触面の法線方向に働く力としては、粒子1に着目した場合、1めり込み量に比例する反発力-kδp、および、2二つの粒子の相対速度に比例する粘性力、ηV2p-V1pの和として与えられる。ここで右肩の添え字pは時刻指標を表す。δpは時刻tpにおける粒子1及び2のめり込み量である。また、V1pおよびV2pは各々、時刻tpにおける粒子1および2の前記接触面法線方向速度である。また、kおよびηは各々比例係数である。これらは図2に示すように、バネ定数kのバネと粘性係数ηのダッシュポットが並列につながったモデルで表される。即ち、非特許文献1にあるように、先行技術においては、粒子1に働く前記接触法線力F1pはF1p=-kδpηV2p-V1p[数式2]と表される。また、粒子2に働く該法線力は、[数式2]の下添え字の1と2を入れ替えればよい。時間刻みΔtを大きくすると、時刻tpで接触判定されても、粒子はすぐに離れ、衝突過程を省略し、計算が高速化できるが、上記[数式2]の反発力の項が大きくなり、後述するように粒子の速度が大きくなりすぎたりして、計算が不安定となったり、誤差が大きくなったりしてしまう。従って、先行技術においては、Δtを細かくとって衝突過程を細かい時間間隔で追跡して計算する必要があった。本発明者等は、上記2の粘性力の表式に着目した。そこで、粒子1の時刻tpにおける該法線方向速度V1pを、粒子1の時刻tp1における速度、即ちV1p1で置き換えたのである。ここで、時刻tp1とは、時刻tpから刻み時間Δtだけ進んだ時刻である。しかしながら、時刻tpの時点で、時刻tp1における正確な速度を把握することは困難なので、V1p1は時刻tp1における予測値E1p1で表す。即ち、本発明においては、粒子1に働く接触法線力F1pは、前記[数式2]のV1pをE1p1で置き換え、F1p=-kδpηV2p-E1p1[数式3]と表される。粒子2に働く接触法線力F2pは、[数式3]の下添え字の1と2を入れ替えればよい。即ちF2p=-kδpηV1p-E2p1[数式4]と表す。予測法線方向速度の算出方法はいくつかある。例えば、本発明になる算出方法によれば、時間刻みΔtを粗くしても、時刻tpにおいて衝突していても、時刻tp1においては粒子は離れている場合が多く、かつ反発力も大きすぎず、精度を保ったまま計算量が減少し、[数式2]を用いた先行技術の場合と比較して計算時間が1/4以下になることが分かった。ここで、注意すべきなのは、粒子2に働く法線方向粘性力は、通常なら粒子1に働く粘性力の反力となるために同じ大きさで反対方向の力となり、それぞれ、作用、反作用の関係にあり、その和はゼロとなることである。本発明においては、[数式4]となり、これは、[数式3]と比較して必ずしも同じ大きさとはならない。従って、接触している2粒子において、これら、法線方向の粘性力の和はゼロではなく、有限な値を持った残留力となり、作用、反作用の関係が満たされず、厳密にいえば物理的に正しい描像を与えない。しかし、そのために生じる誤差は、無視できるほど小さいことが分かった。更に、多数の粒子の衝突を考えた場合、このような残留力は、ランダムな方向と大きさとなり、結果的にその総和はさらに小さくなり、精度に対する影響はほとんどないことが分かったのである。予測法線方向速度粒子1の予測法線速度E1p1を、粒子1が衝突している粒子2のみから受ける力で運動した場合の時刻tp1における粒子1の法線方向速度と考える。従って、次の離散的な運動方程式から求められる。m1E1p1-V1p/Δt=F1p[数式5]上式のF1pに[数式3]を代入し、整理すれば、が得られる。[数式6]を[数式3]に再度代入すれば、時刻tpにおいて粒子1が粒子2より受ける法線方向力は、となる。同様に、粒子2について予測法線速度E2p1及び法線方向力F2pは以下の[数式8]、[数式9]のように表される。以上より、予測法線方向速度E1p1及びE2p1が時刻tpの量で示され、法線方向力も時刻tpにおける値で示されたため、効率良いプログラム作成が可能となり、更に、前述のように精度を保ったまま、計算時間が短縮できるのである。接触判定された粒子の接触後の衝突運動の計算の省略時刻tpで二つの粒子が接触しているとした場合のその後の衝突過程の進展の概念図を図3~図5に示す。従来技術である[数式2]で法線方向力を表した場合、時間刻みΔtを細かくとり、図3に示すように、粒子が離れるまでの衝突過程を細かく計算する必要がある。例えば、同じ法線方向力のままΔtを数倍に取った場合の衝突過程は、図4に示すように、接触後の粒子は次のステップで大きくめり込み、その次のステップでは反発力のために各々高速度で反発してしまい、計算が不安定になり、誤差が大きくなってしまう。次に、大きい時間刻みΔtのままで[数式3]の本発明を用いると、図5に示すように、接触してから次のステップである程度めり込むが、その後、穏やかに離れ、精度を保ったまま衝突過程の計算を省略できる。これは、[数式3]を用いた場合、めり込み量が大きくなり、反発力がかなり大きくなっても、逆方向に働く粘性力も大きくなり、ほどよくバランスされるからである。即ち、[数式3]を用いれば、精度を保ったまま時間刻みを大きくとって衝突過程を省略し、計算速度を短くできるのである。以上の計算ステップの大幅な省略のため、本発明になる計算方法を用いる場合、時間刻みΔtは、従来技術の場合の時間刻みの4倍~10倍とっても、計算結果の精度はほとんど変わらないことが分かった。即ち、従来技術による計算の精度と同一の精度のまま計算時間が1/4~1/10になることが分かった。接触粒子の運動多数の粒子からなる粒子群においては、一つの粒子が同時に複数の粒子と衝突している場合が多い。その場合も、各々の二粒子の衝突について、前記の法線方向力をもとめその合力が粒子に働くと考える。例えば、ある着目粒子iを固定し、他の対向粒子jについて粒子iと衝突しているかどうか判定する。衝突していると判定されたものについては、各々の接触面について法線方向力を[数式7]及び[数式9]において表される二粒子間の法線方向力を求める。これを繰り返し行い、粒子iについてすべての衝突粒子を求めたら、別の粒子に着目して同様な探索を行い、全ての粒子について衝突粒子と、各衝突における法線方向力を求め、これらの和を求める。但し、各々の衝突において接触面法線は異なる方向を向いているからベクトル的な和を求める必要がある。例えば着目粒子iにN個の対向粒子jが衝突しているとすれば、各々のiとjの衝突に対して異なる方向の法線速度を考え、前記[数式3]~[数式9]を適用すればよい。即ち、着目粒子iの、各対向粒子jとの各接触面法線方向の粒子iの時刻tp1での予測速度Ejip1は、[数式6]と同様に考えればよく、となる。ここで、δijpは時刻tpにおける粒子i,jの衝突によるめり込み量、Vijpは時刻tpにおける粒子jの、粒子i、jの接触面法線方向の速度成分である。Vjipは、同様に時刻tpにおける粒子iの前記接触面法線方向の速度成分である。従って、時刻tpにおいて粒子iと粒子jとの衝突から粒子iが受ける法線方向の力Fjipは、[数式7]と同様に考えて次のように表される。この力は、前記のように、時刻tpにおける粒子iと粒子jの接触面に対する法線方向の量であり、そのような法線方向の単位ベクトルをejipであらわせば、上記[数式11]の力の方向を示すベクトルとなる。時刻tpにおいて粒子iがN個の衝突から受ける各々の法線方向力の合力は、[数式11]で表される力を方向ejipを持つベクトルとして対向粒子jについて和を取ることで得る。更に、各々の衝突における接触面の接線方向の力、及び各々の衝突による、回転モーメントも求め、その合力を求める。ここで、各接触面の接線方向の力及び回転モーメントは、非特許文献1にあるように、先行技術を用いて算出することができる。即ち、接触面の接線方向の力については、せん断抗力を与える弾性スプリングと粘性ダッシュポットの並列配置を仮定し、接触点近傍のせん断変形が主として要素間の摩擦力によって生ずるとして計算を行う。回転モーメントについては、要素毎に回転モーメントの運動方程式を立て、変位の運動方程式と同様に差分法による数値解析で求めることができる。これらの合力をもとめたら、先行技術と同様にして離散的な運動方程式を解き、時刻tp1における速度と位置を求める。以上のステップのフローチャートを図6に示す。粒子挙動解析方法多粒子系粒子の挙動を解析するための粒子挙動解析方法において、上記本発明の二粒子が衝突してから離れるまでの法線方向の運動の計算方法を用いることができる。粒子挙動解析方法においては、個々の粒子の運動を運動方程式を用いて離散要素法で計算する。所定の空間領域内に存在するすべての粒子について、刻み時間をΔtとし、各粒子の運動方程式を差分方程式として数値計算を行う。他の粒子と衝突していない粒子については、外場がなければ等速直線運動、重力や電場等の外場があればそのような運動にさらに外場の力による運動が加わる。ここにおいて、任意の2個の粒子が衝突してから離れるまでの運動の計算について、上記本発明の二粒子が衝突してから離れるまでの法線方向の運動の計算方法を用いる。系の初期条件と境界条件を定めた上で、以上の手順で、多粒子系粒子を構成するすべての粒子について挙動解析を行い、時間の経過と共に対象とする多粒子系粒子が全体としてどのように運動するのかを明らかにすることができる。本発明になる計算方法、装置を用いて、大きさの異なる粒子によって構成される粒子群が、重力加速度のもとで種々の装置の中を通過し、最終的に運動を停止するまでを追跡し、最終的に粒子群が形成する塊における粒子の空間粒度分布をもとめることが可能である。この時、粒子群が通過する装置として、シャトルコンベヤや、ホッパー、ドラムフィーダ、フルイ等の機械要素を仮想的に配置し、粒子群を焼結原料の粒を模した粒径分布を持つように決めてやることによって実際の焼結原料装入装置をモデル化することができる。最終的な粒子の空間粒径分布としては、大きい粒子が下になり、小さい粒子が上にくるような分布が形成されることが作業効率上望ましい。しかしながら、前述のように、機械要素の空間的な位置関係、操業条件によってそのような望ましい分布から外れた分布ができてしまう場合が多い。本発明になる計算方法、解析方法を用いて、焼結原料装置をモデル化し、最終的に形成される空間粒径分布を求めることができる。さらに、機械要素の位置、操業条件を変えて計算を実施することによって理想的な空間粒径分布を求めることができる。即ち、焼結原料装入装置の設計、操業条件の決定を著しく高速で行うことができるのである。焼結原料の粒子配置分布解析方法即ち、上記本発明の粒子挙動解析方法を、焼結原料の粒子配置分布解析方法として用いることができる。図7に示すような焼結原料分級設備において、シャトルコンベヤ1からサージホッパー2に投入され,ドラムフィーダ3から切り出され,フルイ4の上を通って落下してパレット5上の堆積物6となるた焼結原料の粒子配置分布を、本発明の粒子挙動解析方法によって求める。まず、シャトルコンベヤ1、サージホッパー2、ドラムフィーダ3、フルイ4を空間に配置し、それぞれの装置に各粒子が接触した際の境界条件を定める。次いで、シャトルコンベア内に搭載された焼結原料の粒子配置分布を初期条件として与える。以上の準備を整えた上で、系内に存在するすべての粒子要素について、運動方程式を刻み時間Δtの差分方程式として数値計算を行う。その際、上記本発明の二粒子が衝突してから離れるまでの法線方向の運動の計算方法を用いる。このような計算を行った結果として、シャトルコンベヤからサージホッパーに投入され,ドラムフィーダから切り出され,フルイの上を通って落下した焼結原料の粒子配置分布を、時間経過と共に明らかにすることができる。接触面法線方向の力にかかる係数[数式3]に用いられる係数kとηはそれぞれ,弾性係数および粘性係数と考えることができるが,単純な定数とするよりも、非特許文献2及び非特許文献3等に開示されたように、次のように設定することで精度のよい計算を効率良く行うことができる。衝突している2粒子を粒子i及び粒子jとする。ここで,各変数は以下のように表わされる。1/E=1-νi2/Ei1-νj2/Ej[数式15]1/r=1/ri1/rj[数式16]1/m=1/mi1/mj[数式17]上記[数式13]~[数式18]の範囲内で,Eはヤング率、νはポアソン比,rは衝突粒子の半径,mは衝突粒子の質量、eは反発係数,πは円周率を表わす。画像処理用演算装置GPUを用いた解析粒子挙動解析の計算に用いる計算機として、画像処理を高速処理するための演算装置であるところのグラフィックプロセシングユニットGPUで計算を行うと、通常の中央演算装置CPUを用いる汎用計算機を用いた同じ計算と比較して、計算時間が1/50以下に短縮できることが分かった。即ち、本発明になる計算方法を用いて作成されたプログラムをGPUに適用すると、従来技術の計算方法を用いたプログラムをCPUに適用した場合に比較して計算時間を1/200以下に短縮することができるのである。ただし、GPUで実施する場合は、計算機が実施するプログラム形態がCPUに適用するプログラム形態とは異なるので、プログラム形態を変更しなければならない。しかしながら、計算のアルゴリズムは同一であり、本発明になる計算方法になるアルゴリズムを用いることは同じであるので、GPUを用いて計算しても、本発明になるものであることは明白である。また、本発明になる計算方法には、GPUによる計算の特徴である並列計算を阻害する要因は全くない。従ってGPUによる計算の特徴を効率良く引き出すことが可能である。[実施例1]汎用計算機を用いた解析数式13~[数式18]に示された係数を用いて、Δt=106秒の時間刻みで、5種類の大きさの球状粒子、2×104個からなる粒子群が、ホッパーからパレットに落下する計算を、汎用計算機を用いてシミュレーションした。従来技術で実施したものは、約50時間でシミュレーションが終了した。同様なシミュレーションを、時間刻みを50倍の5×105秒にして実施したところ、従来技術では、正確な解が得られなかった。しかしながら、時間刻みが5×105のままで本発明になる計算方法を用いると、約5時間でシミュレーションが終了し、従来技術で106の時間刻みで行った計算とほぼ同様な粒子分布を得た。即ち、前記予測速度に起因する誤差は、ほとんど無視できることが分かった。さらに、本発明になるプログラムにより、精度を保ったまま、計算時間が1/10に短縮されたことが分かった。[実施例2]汎用計算機を用いた解析と画像処理用演算装置GPUを用いた解析との比較2×105個の粒子からなる粒子群の挙動の計算を、汎用計算機を用いて1×106秒の時間刻みで従来の計算方法を用いてシミュレーションしたところ、計算終了まで約98日かかった。同様な計算を、時間刻みを5×105秒として、本発明の計算方法を用いて実施したところ、約8日かかった。時間刻み5×105秒で、本発明の計算方法を適用したプログラムをGPUを用いた装置で動作するように変換してGPUを持つ装置で計算した。その結果12時間で計算が終了し、従来法で1×106秒の時間刻みの場合とほぼ同じ結果を得た。従って、汎用計算機を用いた従来技術に対して精度を維持したまま、計算時間が約1/200に短縮された。[実施例3]図7に示す装置を用い、シャトルコンベヤ1からサージホッパー2に投入され,ドラムフィーダ3から切り出され,フルイ4の上を通って落下した焼結原料の分布をシミュレーションするにあたり,本発明による粒子運動の計算方法を、汎用計算機を用いてシミュレーションを行った。まず、シャトルコンベヤ、サージホッパー、ドラムフィーダ、フルイを空間に配置し、それぞれの装置に各粒子が接触した際の境界条件を定めた。次いで、シャトルコンベア内に搭載された焼結原料の粒子配置分布を初期条件として与えた。以上の準備を整えた上で、系内に存在するすべての粒子要素について、運動方程式を刻み時間Δtの差分方程式として数値計算を行った。その際、上記本発明の二粒子が衝突してから離れるまでの法線方向の運動の計算方法を用いた。図8に,フルイ部がプレートと複数の条材からなる装置における粒子挙動の計算結果を示す。落下した粒子は下部のパレット上に積もる。下部のパレットは図中右から左に所定の速度で移動しており,堆積した粒子は図中左方向に運ばれて次の工程に進む。計算では,3種類の大きさの粒子を考慮した大粒子最も明るい、白粒子径6mm,中粒子最も暗い、黒粒子径5mm,小粒子中間の明るさ、灰色粒子径3mm。各々の粒子は、同じ重量密度、3.3g/cm3とした。条材フルイの間隔は先端右下ほど広くなっており,大粒子ほど図右よりに落下する。落下した大粒子は堆積物斜面上を右方向に転がり,斜面右端から堆積物下部に巻き込まれ易くなる。一方小粒子は,フルイ部付け根あたりから落下し,下部バレットの堆積物の上面に降り積もる。これらにより,堆積物は,底から大粒子最も明るい、白,中粒子最も暗い、黒,小粒子中間の明るさ、灰色と分級されることが分かった。この結果を得るのに要した計算時間は10時間であった。さらに、同様な装置構成で同様な粒子挙動を、従来の離散要素法を用いて汎用計算機を用いて計算したところ約100時間かかった。本発明の手法により効率が向上したことが分かる。[実施例4]GPUを用いた画像処理専用計算機に、実施例3で用いた、本発明になるプログラムを移植し、実施例3と同じ要素配置での計算を行った。この時実施例3と同様な結果を得たが、計算に要した時間は2時間となり、汎用計算機を用いた実施例3に比較して1/50以下となった。更に、画像処理専用の計算機を用いたために、途中結果や最終結果のグラフィック表示も効率良く行うことができ、設計効率が著しく向上した。[実施例5]GPUを用いた画像処理専用計算機を用いて、フルイ部がプレートと複数の条材からなる装置において,プレートの長さと角度,条材の開き角度を変えて計算を行った。その結果,堆積した粒子の分布は,プレート長さが0.5~1.0m,角度が40~50度,条材の開き角度が1.5~4.0度の範囲が良いことが分かった。"