"表面に少なくとも1つの凹部又は凸部が成形されたゴルフボールが回転しながら飛行する際の該ゴルフボール周囲の気流を解析してゴルフボールの性能を評価するシミュレーション方法であって、Aコンピューターにおいて設定される仮想空間内に、ゴルフボールモデルと、該ゴルフボールモデルの周囲を取り囲む気流仮想空間を設定し、B該気流仮想空間に格子を生成すると共に、その際にゴルフボールモデルの表面近傍の格子が細かく、表面から離れる方向に格子のサイズが漸次増大するように設定し、C上記気流仮想空間にゴルフボールモデルの前方から所定速度の気流が流入する状態を設定し、D上記気流仮想空間に流入した気流が、初期気流の影響が無くなり安定するまで計算を進行させる気流安定時間を設定し、E該気流安定時間を過ぎた時点からゴルフボールモデルを所定回転数で回転させ、F上記気流安定時間経過後から上記ゴルフボールモデルが1/2回転以上する間、上記気流仮想空間内で上記ゴルフボールモデルの前方から該ボールモデルの周囲を通過し後方へと流れる上記気流の速度、方向及び圧力を上記格子による区画毎に算出し、G各々の格子区画毎に算出された気流の速度、方向及び圧力を積算して、上記気流仮想空間内における上記ゴルフボールモデルの揚力係数及び抗力係数を算出することを特徴とするゴルフボール評価シミュレーション方法。上記気流安定時間を、1≦V×T/X≦20[V:気流速度m/s、T:気流安定時間s、X:気流仮想空間の気流の流れ方向の長さm]を満足するように設定する請求項1記載のゴルフボール評価シミュレーション方法。上記気流仮想空間内に生成させる格子の形状がテトラメッシュ形状である請求項1又は2記載のゴルフボール評価シミュレーション方法。上記気流仮想空間内に生成させる格子の形状がポリゴンメッシュ形状である請求項1又は2記載のゴルフボール評価シミュレーション方法。請求項1~4のいずれか1項に記載のシミュレーション結果から、上記ゴルフボールモデル周囲の気体の流れ方向及び速度を、ベクトルの方向及び長さにより可視化することを特徴とするゴルフボール評価シミュレーション方法。請求項1~4のいずれか1項に記載のシミュレーション結果から、上記ゴルフボールモデル周囲の気体の圧力分布を、等圧線又は等値面により可視化することを特徴とするゴルフボール評価シミュレーション方法。請求項1~4のいずれか1項に記載のシミュレーション結果から、上記ゴルフボールモデル周囲の気体の渦度分布を、渦度の等値線又は等値面により可視化することを特徴とするゴルフボール評価シミュレーション方法。請求項1~4のいずれか1項に記載のシミュレーション結果から、上記ゴルフボールモデル周囲の気体を、気体の流線、流跡線、パーティクルトレース、又はボリュームレンダリングにより可視化することを特徴とするゴルフボール評価シミュレーション方法。本発明は、このような問題に鑑みなされたもので、コンピューターにゴルフボールモデルと、該ゴルフボールモデルの周囲を取り囲む気流仮想空間を設定し、該気流仮想空間に気流が流入する状態を設定してゴルフボールの揚力係数及び抗力係数を算出する際に、正確性を低下させることなく簡易性を向上させ、その算出時間を効果的に短縮することができるゴルフボール評価シミュレーション方法を提供することを目的とする。本発明は、上記目的を達成するため、下記1~8のゴルフボール評価シミュレーション方法を提供する。1表面に少なくとも1つの凹部又は凸部が成形されたゴルフボールが回転しながら飛行する際の該ゴルフボール周囲の気流を解析してゴルフボールの性能を評価シミュレーション方法であって、Aコンピューターにおいて設定される仮想空間内に、ゴルフボールモデルと、該ゴルフボールモデルの周囲を取り囲む気流仮想空間を設定し、B該気流仮想空間に格子を生成すると共に、その際にゴルフボールモデルの表面近傍の格子が細かく、表面から離れる方向に格子のサイズが漸次増大するように設定し、C上記気流仮想空間にゴルフボールモデルの前方から所定速度の気流が流入する状態を設定し、D上記気流仮想空間に流入した気流が、初期気流の影響が無くなり安定するまで計算を進行させる気流安定時間を設定し、E該気流安定時間を過ぎた時点からゴルフボールモデルを所定回転数で回転させ、F上記気流安定時間経過後から上記ゴルフボールモデルが1/2回転以上する間、上記気流仮想空間内で上記ゴルフボールモデルの前方から該ボールモデルの周囲を通過し後方へと流れる上記気流の速度、方向及び圧力を上記格子による区画毎に算出し、G各々の格子区画毎に算出された気流の速度、方向及び圧力を積算して、上記気流仮想空間内における上記ゴルフボールモデルの揚力係数及び抗力係数を算出することを特徴とするゴルフボール評価シミュレーション方法。2上記気流安定時間を、1≦V×T/X≦20[V:気流速度m/s、T:気流安定時間s、X:気流仮想空間の気流の流れ方向の長さm]を満足するように設定する上記1のゴルフボール評価シミュレーション方法。3上記気流仮想空間内に生成させる格子の形状がテトラメッシュ形状である上記1又は2のゴルフボール評価シミュレーション方法。4上記気流仮想空間内に生成させる格子の形状がポリゴンメッシュ形状である上記1又は2のゴルフボール評価シミュレーション方法。5上記1~4のいずれかのシミュレーション結果から、上記ゴルフボールモデル周囲の気体の流れ方向及び速度を、ベクトルの方向及び長さにより可視化することを特徴とするゴルフボール評価シミュレーション方法。6上記1~4のシミュレーション結果から、上記ゴルフボールモデル周囲の気体の圧力分布を、等圧線又は等値面により可視化することを特徴とするゴルフボール評価シミュレーション方法。7上記1~4のシミュレーション結果から、上記ゴルフボールモデル周囲の気体の渦度分布を、渦度の等値線又は等値面により可視化することを特徴とするゴルフボール評価シミュレーション方法。8上記1~4のシミュレーション結果から、上記ゴルフボールモデル周囲の気体を、気体の流線、流跡線、パーティクルトレース、又はボリュームレンダリングにより可視化することを特徴とするゴルフボール評価シミュレーション方法。この評価シミュレーション方法は、上記仮想気流空間内に気流を流入させて回転するゴルフボールモデルの周囲を流れる気流の速度、方向及び圧力を格子区画毎に算出する際、上記D気流仮想空間に流入した気流が、初期気流の影響が無くなり安定するまで計算を進行させる気流安定時間を設定し、次いで上記E該気流安定時間を過ぎた時点からゴルフボールモデルを所定回転数で回転させて、上記気流の速度、方向及び圧力を格子区画毎に算出するようにしたことにより、正確性を低下させることなく算出時間を大幅に短縮することができ、これら格子区画毎の気流の速度、方向及び圧力から容易に揚力係数と抗力係数を算出することができるものである。即ち、図1Aのフローチャートに示したように、従来方法では、1仮想空間内にゴルフボールモデルと気流仮想空間を設定し、2気流仮想空間内に格子を生成させると共に、3ボールモデル表面近傍の格子を細かく、表面から離間するに従って格子を大きくし、この状態で4気流仮想空間にゴルフボール前方から所定速度の気流が流入する状態を設定すると共に、ゴルフボールモデルを所定回転数で回転させて解析をスタートさせ、5ゴルフボールモデルが所定回転する間の解析を進めて気流速度、方向、圧力を求め、これらから揚力係数と抗力係数を算出するのに対し、本発明では、図1Bに示したように、1~3までは従来と同様の手順に従って設定を行った後、4-1気流仮想空間にゴルフボール前方から所定速度の気流が流入する状態を設定してボールを回転させずに解析をスタートさせ、4-2初期気流の影響が無くなり気流が安定するまでの間気流安定時間解析を進め、4-3気流安定時間が経過して気流が安定した時点からゴルフボールモデルを所定回転数で回転させ、5ゴルフボールモデルが1/2回転以上回転する間解析を進めて気流速度、方向、圧力を求め、これらから揚力係数と抗力係数を算出するものである。つまり、従来は初期気流が安定しないままボールを回転して計算するのに対し、本発明は初期気流が安定してからボールを回転させて計算するものである。そしてこのように、本発明のシミュレーション方法では、設定した気流仮想空間に気流を流入させて、該気流仮想空間内のゴルフボールモデル周囲の気流の状態を解析する際、気流仮想空間に流入した気流がボールに当たって方向、圧力、速度などが不規則に変化しながらボールの後方へと流れていく不安定な初期気流状態の間は、ボールを回転させずに解析を進め、この初期気流の影響が無くなり気流が安定した後気流安定時間経過後、ゴルフボールモデルを所定回転数で回転させて1/2回転以上の所定数回転する間解析を進めて気流の速度、方向、圧力を算出し、更にこれらから揚力係数及び抗力係数を算出することにより、初期気流の影響による不安定な気流のデータの混入を排除して安定な気流によるデータのみから揚力係数及び抗力係数の算出が行われ、正確な解析結果を短時間で得ることができるものである。本発明のシミュレーション方法では、風洞試験のような実物モデルによる試験評価を実施することなく、回転しながら飛行するゴルフボール周囲の気流を解析してゴルフボールの揚力係数及び抗力係数を算出することができ、これによりゴルフボールの性能評価を行うことができる。しかも本発明によれば、上記ゴルフボール周囲の気流の解析時間を効果的に短縮して正確な揚力係数及び抗力係数を短時間で算出することができ、コンピューターを用いたシミュレーションによるゴルフボールの性能評価の実施を現実的なものとすることができるものである。以下、図面を参照して、本発明につき更に詳しく説明する。本発明のシミュレーション方法は、表面に少なくとも1つの凹部又は凸部が成形されたゴルフボールが回転しながら飛行する際の該ゴルフボール周囲の気流を解析し、ボールの揚力係数及び抗力係数をコンピューターによる演算を利用して算出する方法である。本発明のシミュレーション方法では、まず、Aコンピューターにおいて設定される仮想空間内に、ゴルフボールモデルと、該ゴルフボールモデルの周囲を取り囲む気流仮想空間を設定する図1Bのフローチャートの1。このコンピューターによるゴルフボールモデルと気流仮想空間の設定に関し、図2A、B及びCにその一例を示す。図2Aは仮想空間内にゴルフボールモデルと気流仮想領域を設定した状態の一例を示す概念斜視図、図2Bは図2Aのゴルフボールモデル1の中心を通る横断面図であり、後述する格子区画の一例を示す概念図、図2Cは図2Bのゴルフボールモデル1の表面部及び表面近傍の格子区画の一例を示す拡大概念図である。本発明では、上記の通り、例えば図2Aに示されるように、仮想空間内にゴルフボールモデル1と該ボールモデル1の周囲を取り囲む気流仮想空間2を設定する。このゴルフボールモデル1は、3DCADなどにより作成することができる。また、気流仮想空間2は、このゴルフボールモデル1を中心にした所定の大きさの直方体形状とすることができ、特に立方体形状とすることが好ましく、より具体的には、例えば3m×3m×3m~15m×15m×15mの立方体とすることができる。更に、この気流仮想空間2の気流の流れ方向の長さXは、特に制限されるものではないが、3~15mとすることが好ましい。この気流仮想空間2は、ゴルフボールの運動に影響を与えるゴルフボール周囲の気体流をほぼ網羅する範囲とすることが必要であるが、ゴルフボールから遠方の気体流はゴルフボールの運動への影響が小さく、また、気流仮想空間の大きさが小さすぎるとシミュレーション精度が低下するため、気流仮想空間2の大きさはシミュレーションの効率又は精度の点から上記範囲が好適である。次に、B該気流仮想空間に格子を生成すると共に、その際にゴルフボールモデルの表面近傍の格子が細かく、表面から離れる方向に格子のサイズが漸次増大するように設定する図1Bのフローチャートの2及び3。具体的には、まず、図3Aに示されるように、ゴルフボールモデル1の表面を例えば一辺を約0.002mmに区画して、三角形、四角形等の多辺形、又は略三角形、略四角形等の略多辺形の面区画11を多数設定し、図2Cに示されるように、この各々の面区画を一面とするゴルフボールモデル表面10に隣接した格子区画21を設定する。ゴルフボールモデル1表面に隣接する格子区画21は、略四角柱形状等の略多角柱形状や略多角錐形状に設定される。そして、このゴルフボールモデル表面10に隣接する格子区画からゴルフボールモデル1から離間する方向に格子区画21の体積が漸次増大するように、気流仮想空間2の残部を格子状に区画して、気流仮想領域2全域を格子区画21で区画する。ここで、上記気流仮想空間2に形成される格子区画の形状は、ポリゴンメッシュ多面体、テトラメッシュ四面体、プリズムメッシュ三角柱、ヘキサメッシュ六面体、又はこれらが混在した形状など、適宜な形状とすることができる。なお、これらの中でも特にポリゴンメッシュ形状又はテトラメッシュ形状が好ましく採用される。ゴルフボール周囲の気体流は、ゴルフボールに近いものの方がより大きく影響するため、図2Bに示されるように、格子区画をゴルフボールモデル1近傍は細かく、気体流が与える影響が小さいゴルフボールモデル1から遠方は粗く設定するのは上述のとおりであるが、この場合ゴルフボールモデル1のボール表面から離間する方向への格子区画の体積の増大は、連続的なものであっても、段階的なものであってもよい。次に、C上記気流仮想空間2にゴルフボールモデル1の前方から所定速度の気流が流入する状態を設定し図1Bのフローチャートの4-1、D上記気流仮想空間2に流入した気流が、初期気流の影響が無くなり安定するまで計算を進行させる気流安定時間を設定し図1Bのフローチャートの4-2、E該気流安定時間を過ぎた時点からゴルフボールモデル1を所定回転数で回転させる図1Bのフローチャートの4-3。上記気流の速度は、特に制限されるものではなく、想定されるゴルフボールの飛行速度等に応じて適宜設定され、通常は5~90m/sの範囲の任意の速度とすることができる。上記気流安定時間は、初期気流の不安定な状態が解消し気流仮想空間2内が安定した気流に満たされるまでの時間とすることが好ましく、気流の速度Vや気流仮想空間2の気流の流れ方向の長さX図2A参照に応じて適宜設定される。具体的には、気流速度Vm/s、気流安定時間Ts、気流仮想空間2の気流の流れ方向の長さXmが、1≦V×T/X≦20の関係を満足するように設定すればよい。即ち、気流の流入開始から気流が気流仮想空間2の下流面に到達するまでの時間乃至その20倍の時間の範囲で適宜設定すればよく、より好ましくは1≦V×T/X≦8の範囲に設定すればよい。上記気流安定時間T経過後、ゴルフボールモデル1を回転させるが、その回転数も使用クラブの番手やヘッドスピードなどから想定される飛行中の回転数に応じて適宜設定され、特に制限されるものではないが、通常は15000rpm以下の任意の回転数とすることが好ましい。次に、F上記気流安定時間経過後から上記ゴルフボールモデルが1/2回転以上する間、上記仮想空間内で上記ゴルフボールモデルの前方から該ボールモデルの周囲を通過し後方へと流れる上記気流の速度、方向及び圧力を上記格子による区画毎に算出する図1Bのフローチャートの5の前段。即ち、気流が気流仮想空間2に流入してゴルフボールモデル1に当たったときに生じる運動要素は、三次元空間座標系の各軸方向の気流の速度、気流の方向及びボールモデル1表面に対する気流の圧力であり、計算に当たって用いる基礎方程式、つまり、以下に示される質量保存則に相当する連続の式1~3及び物体の運動保存則に相当するナビエストークの方程式4~6に数値を代入して算出することができる。上記ゴルフボールモデル1の周囲に、図2Aに示すように、気体を矢印方向に流すシミュレーションでは、気流仮想空間2の格子区画毎に空気の流れを演算により解析することができる。この演算に上記式1~6を用い、気流仮想空間2を格子区画に区切ったことに対応して上記式1~6を離散化して演算を行うことができる。シミュレーションの方法は、シミュレーションの条件等を考慮して有限差分法、有限体積法、境界要素法、有限要素法等を適宜選択して行うことができる。ここで、本発明では、上記のように、格子による区画毎の気流の速度、方向及び圧力をゴルフボールモデルが1/2回転以上する間計算する必要がある。解析時間内のゴルフボールモデル1の回転数が1/2回転未満であると正確な気流の速度、方向及び圧力が得られない場合があり、正確性が低下して本発明の目的を達成することが出来ない。一方、解析時間内の回転数が多ければ、それだけ正確性を向上させることができるが計算時間が長くなる不利があり、コンピューターの性能によっては十分な時間短縮を達成することができない場合がある。このため、解析時間内のゴルフボールモデル1の回転数は1/2以上であればよく特に制限されるものではないが、好ましくは1/2回転~10回転、より好ましくは1/2回転~7回転である。次に、G気流仮想空間内2の各々の格子区画において算出された気流の速度、方向及び圧力を積算して、上記気流仮想空間内におけるゴルフボールモデルの揚力係数及び抗力係数を算出する図1Bのフローチャートの5の後段。更に本発明では、上記シミュレーション方法により揚力係数及び抗力係数を算出し、その結果から、下記1~4の通り、ゴルフボールモデル1の周囲の気体の流れ状態、圧力、渦度などを可視化することもできる。1図4のとおり、上記ゴルフボールモデル1周囲の気体の流れ方向及び速度を、ベクトルの方向及び長さにより可視化すること。2図5のとおり、上記ゴルフボールモデル1周囲の気体の圧力分布を、等圧線又は等値面により可視化すること。3図6のとおり、上記ゴルフボールモデル1周囲の気体の渦度分布を、渦度の等値線又は等値面により可視化すること。4図7のとおり、上記ゴルフボールモデル1周囲の気体を、気体の流線、流跡線、パーティクルトレース、又はボリュームレンダリングにより可視化すること。本発明のシミュレーション方法では、風洞試験のような実物モデルによる試験評価を実施することなく、回転しながら飛行するゴルフボール周囲の気流を解析してゴルフボールの揚力係数及び抗力係数を算出することができ、これによりゴルフボールの性能評価を行うことができる。しかも本発明によれば、上記ゴルフボール周囲の気流の解析時間を効果的に短縮して正確な揚力係数及び抗力係数を短時間で算出することができ、コンピューターを用いたシミュレーションによるゴルフボールの性能評価の実施を現実的なものとすることができるものである。以下、実施例及び比較例を示し、本発明を具体的に説明するが、本発明は下記実施例に制限されるものではない。実施例1~5及び比較例1,2ディンプルの設計が異なる2種類のゴルフボールモデルモデル1、モデル2につき、図1に示された解析手法A又はBに従って、下記表1に示した条件で揚力係数CL及び抗力係数CDを算出し、両ゴルフボールモデルの性能を比較評価した。各シミュレーションによる揚力係数CL評価/抗力係数CD評価の精度、及び計算時間を下記基準で評価した。なお、シミュレーションに用いたコンピューターの詳細は下記の通りである。コンピューターの詳細OS:Windows登録商標XPCPU:intelXeonX5482320GHzL2=6MB×2FSB=1600MHz2CPU8Coreメモリ:64GB4096MB×16DDR2800ECCREGFBDIMMCD評価/CL評価の精度判定の基準◎:安定した値で比較評価ができる。○:十分に比較評価ができる。△:比較評価はできる。×:比較評価ができない。計算時間◎:30時間以内○:60時間以内△:90時間以内×:90時間超表1のとおり、図1のフローチャートBに示された本発明にかかるシミュレーション方法に従って行われた実施例1~5の性能評価によれば、ゴルフボールの正確な性能評価を短時間で行うことができる。それに対し、図1のフローチャートAに示された従来のシミュレーション方法に従って行われた比較例1の性能評価では、240時間もの長時間を要した。この比較例1と解析手法以外は同一の条件で行われた実施例4が1/4の60時間で同等の精度で評価が行われており、本発明による時間短縮の効果は絶大である。また、ゴルフボールモデルが1/4回転の間でCD、CLを算出した比較例2では、ボールの回転が足りないために安定した正確なCD値、CL値を得ることができず、例えばCL値は、他の例とは逆にモデル1よりもモデル2が大きな値となってしまった。このため、比較例2では、両ボールモデルの性能を比較評価することができなかった。"