"解析空間を要素で分割して解析メッシュを設定する解析メッシュ設定手段と、上記解析空間の予め設定する条件下での上記各要素の流体力学に関するパラメータを反復計算する計算実行手段と、上記反復計算により得られた上記各要素の上記パラメータが所定に収束したか否か判定する収束判定手段と、上記各要素の上記パラメータが所定に収束したと判定された場合、上記反復計算の結果を出力する出力手段と、を備えた流体解析装置において、上記反復計算により得られた上記各要素の上記パラメータが予め設定しておいた閾値を超えた場合、該閾値を超えた上記パラメータを有する要素を特定して分離する要素分離手段と、上記分離した要素に流れが入り込まない予め設定する抵抗値を設定する抵抗値設定手段と、上記抵抗値設定手段で設定した抵抗値を上記パラメータに応じて可変する抵抗値可変手段と、を備えたことを特徴とする流体解析装置。上記流体力学に関するパラメータは、少なくとも圧力であることを特徴とする請求項1記載の流体解析装置。解析空間を要素で分割して解析メッシュを設定する解析メッシュ設定ステップと、上記解析空間の予め設定する条件下での上記各要素の流体力学に関するパラメータを反復計算する計算実行ステップと、上記反復計算により得られた上記各要素の上記パラメータが所定に収束したか否か判定する収束判定ステップと、上記各要素の上記パラメータが所定に収束したと判定された場合、上記反復計算の結果を出力する出力ステップと、を備えた流体解析方法において、上記反復計算により得られた上記各要素の上記パラメータが予め設定しておいた閾値を超えた場合、該閾値を超えた上記パラメータを有する要素を特定して分離する要素分離ステップと、上記分離した要素に流れが入り込まない予め設定する抵抗値を設定する抵抗値設定ステップと、上記抵抗値設定ステップで設定した抵抗値を上記パラメータに応じて可変する抵抗値可変ステップと、を備えたことを特徴とする流体解析方法。上記流体力学に関するパラメータは、少なくとも圧力であることを特徴とする請求項3記載の流体解析方法。しかしながら、上述の特許文献1に開示される方法では、メッシュモデルの生成には、予め解析結果から解析誤差を算出し、その解析誤差を評価してメッシュモデルの粗密分布制御の対象となる領域を判定し、その領域に存在する有限要素群を結合して形状モデルを創成し、その形状モデルを分割しなければならないという問題があった。また、上述の特許文献2に開示される方法では、本来、流れ場に流路等があるところを個体領域として閉塞してしまうため、計算結果が実際の現象と異なったものとなってしまう虞がある。本発明は上記事情に鑑みてなされたもので、予め解析誤差を評価してメッシュモデルの粗密分布を作成するような計算の工数も必要なく、また、実際の流れ場を忠実に再現して実現象に近い解析結果を安定して得ることができる流体解析装置、及び、流体解析方法を提供することを目的としている。本発明の流体解析装置の一態様は、解析空間を要素で分割して解析メッシュを設定する解析メッシュ設定手段と、上記解析空間の予め設定する条件下での上記各要素の流体力学に関するパラメータを反復計算する計算実行手段と、上記反復計算により得られた上記各要素の上記パラメータが所定に収束したか否か判定する収束判定手段と、上記各要素の上記パラメータが所定に収束したと判定された場合、上記反復計算の結果を出力する出力手段とを備えた流体解析装置において、上記反復計算により得られた上記各要素の上記パラメータが予め設定しておいた閾値を超えた場合、該閾値を超えた上記パラメータを有する要素を特定して分離する要素分離手段と、上記分離した要素に流れが入り込まない予め設定する抵抗値を設定する抵抗値設定手段と、上記抵抗値設定手段で設定した抵抗値を上記パラメータに応じて可変する抵抗値可変手段とを備えた。また、本発明の流体解析方法の一態様は、解析空間を要素で分割して解析メッシュを設定する解析メッシュ設定ステップと、上記解析空間の予め設定する条件下での上記各要素の流体力学に関するパラメータを反復計算する計算実行ステップと、上記反復計算により得られた上記各要素の上記パラメータが所定に収束したか否か判定する収束判定ステップと、上記各要素の上記パラメータが所定に収束したと判定された場合、上記反復計算の結果を出力する出力ステップとを備えた流体解析方法において、上記反復計算により得られた上記各要素の上記パラメータが予め設定しておいた閾値を超えた場合、該閾値を超えた上記パラメータを有する要素を特定して分離する要素分離ステップと、上記分離した要素に流れが入り込まない予め設定する抵抗値を設定する抵抗値設定ステップと、上記抵抗値設定ステップで設定した抵抗値を上記パラメータに応じて可変する抵抗値可変ステップとを備えた。本発明による流体解析装置、及び、流体解析方法によれば、予め解析誤差を評価してメッシュモデルの粗密分布を作成するような計算の工数も必要なく、また、実際の流れ場を忠実に再現して実現象に近い解析結果を安定して得ることが可能となる。以下、図面に基づいて本発明の実施の形態を説明する。図1において、符号1は、ROM,RAM等のメモリ及びCPUを有するコンピュータで構成された流体解析装置を示し、この流体解析装置1には、入力側にキーボード、マウス等を含む入力装置2が接続され、出力側にモニタ等の表示装置3が接続されている。流体解析装置1は、解析モデル設定処理部11、モデル解析処理部12、メモリ部13、及び、表示処理部14を備えて主要に構成されており、各処理部11,12,14の機能を実現するためのプログラムをメモリ部13にロードして実行することにより解析処理が実行される。解析モデル設定処理部11は、入力装置2からの入力信号により、例えば、予めメモリ部13に記憶した解析空間Aに対し、予め用意された流体解析用のメッシュ生成ソフトを用いて、所定形状例えば、各要素が所定の大きさの三角形形状で表現される三角メッシュ等の解析メッシュを設定し、解析に係る様々な条件例えば、初期状態における圧力P0、流速V0、密度ρ0等を各要素に設定するようになっている。そして、これら解析メッシュの各要素にはそれぞれを識別するIDが付される。尚、解析メッシュの形状としては、3角形形状に限るものではなく、他に、4角形形状の解析メッシュ等であっても良く、解析空間Aが3次元空間のものでは、4面体、6面体等の立体形状の解析メッシュであっても良い。このように、解析モデル設定処理部11は、解析メッシュ設定手段として設けられている。本実施の形態では、解析空間Aとして、図3に示すような、長方形形状の空間に、二つの部品固体領域C1、C2が、空間の中央下部に近接して設けられた解析モデルを例に、左端から右端に流れが生じる場合を説明する。そして、この解析空間内には、三角メッシュが、部品固体領域C1、C2を除いた全ての領域に生成される。ここで、図4に示すように、解析空間A内に設定される三角メッシュは、部品固体領域C1、C2の間では、この間が極めて狭くなっているため、他の領域とは異なった、歪な形状のメッシュが形成されることになる。モデル解析処理部12は、上述の解析空間Aに対し、予め設定する条件下での各要素の流体力学に関するパラメータ圧力P、流速V等を公知のCFDComputationalFluidDynamicsの理論に基づいて反復計算し、この反復計算により得られた各要素のパラメータが所定に収束したか否か判定し、各要素のパラメータが所定に収束したと判定された場合、反復計算の結果を表示処理部14に出力する。ここで、計算結果が収束しない場合には、反復計算により得られた各要素の圧力Pが予め設定しておいた閾値Pcを超えた場合、この圧力Pが予め設定しておいた閾値Pcを超えた要素を特定して分離しIDを特定し、この分離した要素に流れが入り込まない予め設定する抵抗値Rmaxを設定し、この設定した抵抗値を、以後繰り返される計算によって、圧力Pに応じて可変設定するようになっている。このように、モデル解析処理部12は、計算実行手段、収束判定手段、出力手段、要素分離手段、抵抗値設定手段、抵抗値可変手段として設けられている。メモリ部13は、ハードディスク装置や光磁気ディスク装置、フラッシュメモリ等の読み書き自在な不揮発性記録媒体で構成されており、解析モデル設定処理部11及びモデル解析処理部12における処理内容等を適宜記憶する。表示処理部14は、解析処理を実行している際の画像、解析結果の画像等を所定にポスト処理して表示装置3へ出力する。尚、解析モデル設定処理部11、モデル解析処理部12、表示処理部14は、コンピュータによって処理される必要はなく、それぞれ専用のハードウェアにより実現されるものであってもよい。次に、流体解析装置1で実行される流体解析プログラムを、図2に示すフローチャートで説明する。まず、ステップ以下、Sと略称101では、解析モデル設定処理部11は、入力装置2からの入力信号により、解析空間Aを読み込んで、メモリ部13に記憶する。次に、S102に進み、解析モデル設定処理部11は、S101で読み込んだ解析空間Aに対し、予め用意された流体解析用のメッシュ生成ソフトを用いて、所定形状例えば、各要素が所定の大きさの三角形形状で表現される三角メッシュ等の解析メッシュを設定する。この際、解析メッシュの各要素には、それぞれを識別するIDが付される。次いで、S103に進み、解析に係る様々な条件例えば、初期状態における圧力P0、流速V0、密度ρ0等が各要素に設定され、解析空間Aに対する流れの環境例えば、流速の方向等、収束条件例えば、所定領域における前回と今回の流速の差が、予め設定する閾値以下になった場合に収束と判定する等の設定がなされる。尚、収束条件は、前回と今回の流速の差以外にも、前回と今回の圧力の差が、予め設定する閾値以下になった場合に収束と判定する等の条件としても良い。そして、S104に進み、モデル解析処理部12は、設定された解析空間Aに対し、公知のCFDComputationalFluidDynamicsの理論に基づいて反復計算を実行し、解析メッシュを構成する全ての要素について計算が終了した後、S105に進んで、例えば、解析空間Aにおける所定領域における前回と今回との流速の差が予め設定する閾値以下であって収束条件が満足されているか否か判定する。このS105の判定の結果、収束条件が満足されているのであれば、S106に進み、表示処理部14は、ポスト処理を行って、表示装置3に結果を出力する。一方、モデル解析処理部12は、S105の判定の結果、収束条件が満足されない場合は、S107に進み、全ての解析メッシュを構成する要素で、圧力Pと予め設定しておいた閾値Pcとが比較され、圧力Pが閾値Pc以上となっている要素がある場合は、S108に進み、圧力Pが閾値Pc以上P≧Pcとなっている要素のIDを区別することで分離処理する。そして、このP≧Pcのメッシュを分離処理した後は、S110に進み、分離処理したメッシュに流れが入り込まない予め設定する抵抗値R=Rmaxを設定して、再び、S104に戻って、反復計算を実行させる。また、S107の判定の結果、全てのメッシュで、P<Pcとなった場合はS109に進み、分離されたメッシュが有るか否か判定され、分離されたメッシュがある場合は、S110に進み、設定された抵抗値Rを、例えば、図5に示すような、圧力Pの関数の値に応じて可変設定し、再び、S104に戻って、反復計算を実行させる。可変設定される抵抗値Rは、図5の関数では、P≧Pcの領域では、メッシュに流れが入り込まない予め設定する抵抗値Rmaxが設定され、P<Pcの領域では、圧力Pが低下するに従ってすなわち、計算が安定してくるに従って、低い値が設定され、最終的には、0に設定されるような特性となっている。尚、抵抗値の特性は、図5の特性に限るものではない。また、S109の判定の結果、分離されたメッシュが無い場合ば、再び、S104に戻って、反復計算を実行させる。このように本発明の実施の形態によれば、メッシュを設定した解析空間Aに対し、予め設定する条件下での各要素の流体力学に関するパラメータ圧力P、流速V等を公知のCFDComputationalFluidDynamicsの理論に基づいて反復計算し、計算結果が収束しない場合に、反復計算により得られた各要素の圧力Pが予め設定しておいた閾値Pcを超えた場合、この圧力Pが予め設定しておいた閾値Pcを超えた要素を特定して分離しIDを特定し、この分離した要素に流れが入り込まない予め設定する抵抗値Rmaxを設定し、この設定した抵抗値を、以後繰り返される計算によって、圧力Pに応じて可変設定するようになっている。このため、予め解析誤差を評価してメッシュモデルの粗密分布を作成するような計算の工数も必要なく、解析計算を実行することが可能である。また、分離メッシュに設定される抵抗値Rは、該分離メッシュの圧力Pの値に応じて、解析結果が安定するに従って、次第に0に設定されるので、実際の流れ場を忠実に再現して実現象に近い解析結果を安定して得ることが可能となる。図6aは、本実施の形態の方法により得られた解析結果所定圧力の等圧線分布を示し、図6bは、二つの部品固体領域C1、C2の間にも十分な品質のメッシュが生成されるように、解析空間全体に細かいメッシュを設定し、計算を実行して得た解析結果を示すものである。図6bによる計算結果は、メッシュの目が細かいため、図6aの解析計算よりも収束するまでに、長い時間がかかっている。本発明の実施の形態による解析方法では図6aでは、図6bと同等の計算結果を得るまでに、はるかに短い時間で収束させることが可能となっている。尚、本実施の形態では、長方形形状の空間に、二つの部品固体領域C1、C2が、空間の中央下部に近接して設けられた解析モデルを例に、説明したが、このような解析モデルは、他の形状であっても良く、また、三次元モデルであっても適用できることは云うまでもない。"