"表面に複数のディンプルが形成されたゴルフボールを打出したときの弾道をコンピューターによる演算を利用して解析する弾道シミュレーション方法であって、Aボールの打出される気流仮想空間内フィールドに、格子を生成し、B表面に複数のディンプルが形成された略球状のゴルフボールモデルを気流仮想空間内フィールドに設定し、Cその際、ゴルフボールモデルの近傍は、格子を仮想空間内フィールドの格子より細かくすると共に、ボールモデル表面から離れる方向に格子のサイズが漸次増大するように区画し、Dゴルフボールモデルに重量を設定し、Eゴルフボールモデルに初期条件初期速度、打出し角度、回転数を与え、Fゴルフボールモデルが初期条件による運動を開始し、ボールモデルの移動に追従するように、ボール近傍の格子は元に設定した格子サイズを崩さずに移動部で再度生成し、ボールが通過した後の格子は気流仮想空間内フィールドの格子に戻し、G気流仮想空間内フィールドの各々の格子区画において算出された気体流の速度、方向及び圧力を積算して、上記気流仮想空間内フィールドにおける飛行中のゴルフボールの揚力係数及び抗力係数を算出し、H飛行中のゴルフボールの高さ変化、横方向変化、速度変化、回転数変化を算出し、打出時から落下時までのゴルフボールの飛距離、左右へのばらつきを算出して、弾道を推定する際に、飛行中のゴルフボールモデル周囲のレイノルズ数Reを算出すると共に、Re=250000~170000を高速領域、Re=170000~100000を中速領域、及びRe=100000以下を低速領域とする基準に従って、ゴルフボールモデルの飛行行程を2つ又は3つの領域に分け、上記高速領域におけるゴルフボールモデル表面と接する上記格子の高さ格子高さA、中速領域におけるゴルフボールモデル表面と接する上記格子の高さ格子高さB及び低速領域におけるゴルフボールモデル表面と接する上記格子の高さ格子高さCを、格子高さA<格子高さB<格子高さCとなるように設定することを特徴とするゴルフボールの弾道シミュレーション方法。D-2上記気流仮想空間内においてゴルフボールモデルに所定速度の気流が流入する状態を設定する請求項1記載のゴルフボールの弾道シミュレーション方法。上記高速領域における上記格子高さAを0.001~0.015mmの範囲、中速領域における上記格子高さBを0.005~0.02mmの範囲、低速領域における上記格子高さCを0.01~0.025mmの範囲でそれぞれ設定する請求項1又は2記載のゴルフボールの弾道シミュレーション方法。上記気流仮想空間フィールドの下面に地面を想定した壁面を設定して、該気流仮想空間フィールドの下面を閉鎖面とする請求項1~3のいずれか1項に記載のゴルフボールの弾道シミュレーション方法。上記Eで設定する初期条件に、地面を想定した上記気流仮想空間フィールドの下面からゴルフボールモデル最下点までの距離を打撃高さとして加える請求項4記載のゴルフボールの弾道シミュレーション方法。格子形状が直交適合格子である請求項1~5のいずれか1項に記載のゴルフボールの弾道シミュレーション方法。格子形状が非構造格子である請求項1~5のいずれか1項に記載のゴルフボールの弾道シミュレーション方法。算出されるゴルフボール落下時の速度、角度、回転数から、ゴルフボールが落下した後ゴルフボールが静止するまでの距離を算出し、打出しから静止するまでの飛距離を算出する請求項1~7のいずれか1項に記載のゴルフボールの弾道シミュレーション方法。本発明は、このような問題に鑑みなされたもので、表面に複数のディンプルが形成されたゴルフボールの弾道を実物モデルによる試験評価によらずに評価して、効率的にゴルフボールを開発する方法を提供することを目的とする。本発明は、上記目的を達成するため、下記1~8のゴルフボールの弾道シミュレーション方法を提供する。1表面に複数のディンプルが形成されたゴルフボールを打出したときの弾道をコンピューターによる演算を利用して解析する弾道シミュレーション方法であって、Aボールの打出される気流仮想空間内フィールドに、格子を生成し、B表面に複数のディンプルが形成された略球状のゴルフボールモデルを気流仮想空間内フィールドに設定し、Cその際、ゴルフボールモデルの近傍は、格子を仮想空間内フィールドの格子より細かくすると共に、ボールモデル表面から離れる方向に格子のサイズが漸次増大するように区画し、Dゴルフボールモデルに重量を設定し、Eゴルフボールモデルに初期条件初期速度、打出し角度、回転数を与え、Fゴルフボールモデルが初期条件による運動を開始し、ボールモデルの移動に追従するように、ボール近傍の格子は元に設定した格子サイズを崩さずに移動部で再度生成し、ボールが通過した後の格子は気流仮想空間内フィールドの格子に戻し、G気流仮想空間内フィールドの各々の格子区画において算出された気体流の速度、方向及び圧力を積算して、上記気流仮想空間内フィールドにおける飛行中のゴルフボールの揚力係数及び抗力係数を算出し、H飛行中のゴルフボールの高さ変化、横方向変化、速度変化、回転数変化を算出し、打出時から落下時までのゴルフボールの飛距離、左右へのばらつきを算出して、弾道を推定する際に、飛行中のゴルフボールモデル周囲のレイノルズ数Reを算出すると共に、Re=250000~170000を高速領域、Re=170000~100000を中速領域、及びRe=100000以下を低速領域とする基準に従って、ゴルフボールモデルの飛行行程を2つ又は3つの領域に分け、上記高速領域におけるゴルフボールモデル表面と接する上記格子の高さ格子高さA、中速領域におけるゴルフボールモデル表面と接する上記格子の高さ格子高さB及び低速領域におけるゴルフボールモデル表面と接する上記格子の高さ格子高さCを、格子高さA<格子高さB<格子高さCとなるように設定することを特徴とするゴルフボールの弾道シミュレーション方法。2D-2上記気流仮想空間内においてゴルフボールモデルに所定速度の気流が流入する状態を設定する1記載のゴルフボールの弾道シミュレーション方法。3上記高速領域における上記格子高さAを0.001~0.015mmの範囲、中速領域における上記格子高さBを0.005~0.02mmの範囲、低速領域における上記格子高さCを0.01~0.025mmの範囲でそれぞれ設定する1又は2記載のゴルフボールの弾道シミュレーション方法。4上記気流仮想空間フィールドの下面に地面を想定した壁面を設定して、該気流仮想空間フィールドの下面を閉鎖面とする1~3のいずれかに記載のゴルフボールの弾道シミュレーション方法。5上記Eで設定する初期条件に、地面を想定した上記気流仮想空間フィールドの下面からゴルフボールモデル最下点までの距離を打撃高さとして加える4記載のゴルフボールの弾道シミュレーション方法。6格子形状が直交適合格子である1~5のいずれかに記載のゴルフボールの弾道シミュレーション方法。7格子形状が非構造格子である1~5のいずれかに記載のゴルフボールの弾道シミュレーション方法。8算出されるゴルフボール落下時の速度、角度、回転数から、ゴルフボールが落下した後ゴルフボールが静止するまでの距離を算出し、打出しから静止するまでの飛距離を算出する1~7のいずれかに記載のゴルフボールの弾道シミュレーション方法。この弾道シミュレーション方法によれば、実物モデルによる試験評価を実施することなく、表面に複数のディンプルが形成されたゴルフボールを任意の初速度、任意の回転数及び任意の角度で打出したときのゴルフボールの弾道を推定することができ、更にフォロー、アゲンスト、横風などの風の状態に応じた弾道を推定することも可能である。この場合、本発明では、上記のように、飛行中のゴルフボールモデル周囲のレイノルズ数Reを算出すると共に、Re=250000~170000を高速領域、Re=170000~100000を中速領域、及びRe=100000以下を低速領域とする基準に従って、ゴルフボールモデルの全飛行行程を2つ又は3つの領域に分け、それぞれの領域毎にゴルフボール表面と接する上記格子の高さ格子高さA~Cを設定することにより、十分な精度を担保しつつ計算時間の短縮を図るものである。例えば、ボールが高速で打出されるドライバーショットなどでは、その飛行行程を、ボールが比較的高速で飛行する前域を高速領域Re=250000~170000、中速で飛行する中盤を中速領域Re=170000~100000、比較的低速で飛行する後域を低速領域Re=100000以下の3つの領域に分け、高速領域における上記格子高さ格子高さA、中速領域における上記格子高さ格子高さB及び低速領域における上記格子高さ格子高さCを格子高さA<格子高さB<格子高さCとなるように設定することにより、計算時間の短縮を図るものである。つまり、レイノルズ数Reの大きな高速領域ではボール表面近傍を流れる境界層が薄く、レイノルズ数Reが小さくなる低速領域に移行するに従ってこの境界層は厚くなる。一方、本発明は上記のように、ボールモデル周囲に生成させた格子区画毎に気体流の速度、方向及び圧力を算出するものであるが、本発明者は、上記境界層が厚くなる上記レイノルズ数Reの小さな領域では、ゴルフボールモデル周囲の格子高さを高くして格子区画の大きさを比較的大きくすることにより格子区画の数を少なくしても実用上問題のない精度を担保し得ることを見出し、更に検討を重ね、レイノルズ数Reによって全飛行行程を上記高速領域、中速領域、低速領域の3つの領域に分けてRe値が小さい領域ほどゴルフボールモデル表面と接する格子の高さを大きく設定して、精度を損なわない範囲で格子区画の数を少なくすることにより、良好な精度を維持しつつ計算時間を効果的に短縮し得ることを見出し、本発明を完成したものである。なお、飛行行程初期のRe値が170000に満たないアイアンショットなどを想定する場合には、飛行行程は上記基準に従って中速領域Re=170000~100000と低速領域Re=100000以下の2つの領域に分けられ、格子高さB<格子高さCとなるようにそれぞれの領域における格子高さが設定される。このように、本発明によれば、実物モデルによる試験評価を実施することなく、ゴルフボールを任意の初速度、任意の回転数及び任意の角度で打出したときのゴルフボールの弾道を推定することができる。しかも、その推定シミュレーションを行う際に生成させる格子区画の数を調整して効率よく計算を行なうことにより全計算時間を短縮することができ、その結果、ゴルフボールの表面形状ディンプルの形状、配置、大きさなどの評価の時間が短縮され、また評価の正確性、客観性も十分に担保されることから、ゴルフボールの種類やグレードに合わせてより高性能の製品を効率よく開発することができる。以下、図面を参照して、本発明につき更に詳しく説明する。本発明の弾道シミュレーション方法は、表面に複数のディンプルが形成されたゴルフボールを打出したときの弾道をコンピューターによる演算を利用して解析する方法である。本発明の弾道シミュレーション方法では、まず、コンピュータによりAボールの打ち出される気流仮想空間フィールドを設定し、図1Aに示したように、その気流仮想空間内フィールドに格子を生成する。なお、特に制限されるものではないが、この気流仮想空間内フィールドに形成される格子の形状としては、直交適合格子や非構造格子とすることができる。この場合、この気流仮想空間フィールドは、ゴルフボールが打出されてから飛行し、少なくとも着地するまでの領域全体を表すものであり、後述するゴルフボールモデルはこの気流仮想空間内フィールドを移動する。また、この気流仮想空間フィールドの下面は、開放された面として設定してもよいが、地面を想定した壁面を設定して閉鎖面とすることもできる。これにより、気流仮想空間フィールドの下面から空気が出て行く要素が計算から除かれ、より精度の高いシミュレーションを行うことができる。更に、実際にゴルフをプレーする場合には、ボールをティーに載せて打撃することも行われるが、この地面を想定した壁面を設けることにより、初期条件に打撃高さとしてこの壁面からボール最下点までの距離を設定することができ、より実際の使用時に近い条件でより精度の高いシミュレーションを行うことができるものである。次に、図1に示されているように、Bこの気流仮想空間内フィールドに、表面に複数のディンプルが形成された略球状のゴルフボールモデルを設定する。この場合、このゴルフボールモデルは、3DCADなどにより作成することができ、例えば図2に示したゴルフボールモデルを例示することができる。更に、このゴルフボールモデルを気流仮想空間内フィールドに設定する際、Cゴルフボールモデルの近傍は、格子を仮想空間内フィールドの格子より細かくすると共に、ボールモデル表面から離れる方向に格子のサイズが漸次増大するように区画する。このようにゴルフボールモデル近傍の格子を区画することにより、余分な計算を省くことができ、より効率的に計算を行うことができる。具体的には、まず、図2Aに示したように、ゴルフボールモデルの表面を例えば一辺を約0.002mmに区画して、三角形、四角形等の多辺形、又は略三角形、略四角形等の略多辺形の面区画を多数設定し、この各々の面区画を一面とするゴルフボールモデル表面に隣接した格子区画を設定する。ゴルフボールモデル表面に隣接する格子区画は、略四角柱形状等の略多角柱形状や略多角錐形状に設定される。そして、図1Bに示したように、このゴルフボールモデル表面に隣接する格子区画からゴルフボールから離間する方向に格子区画の体積が漸次増大するように、気流仮想空間フィールドの残部を格子状に区画して、気流仮想領域全域を格子区画で区画する。このゴルフボールモデル表面に隣接する格子区画以外の気流仮想空間フィールドの残部に形成される格子区画の形状としては、六面体、三角柱状五面体、四角錐状五面体、三角錐状四面体などの多面体が挙げられ、これらを適宜組み合わせて設定することができる。ゴルフボール周囲の気体流は、ゴルフボールに近いものの方がより大きく影響するため、図1Bに示されるように、格子区画をゴルフボールモデル近傍は細かく、気体流が与える影響が小さいゴルフボールモデルから遠方は粗く設定するのは上述の通りであるが、この場合ゴルフボールモデルのボール表面から離間する方向への格子区画の体積の増大は、連続的なものであっても、段階的なものであってもよい。ここで、本発明では、飛行中のゴルフホールモデル周囲のレイノルズ数Reを算出して、そのRe値に応じて飛行行程を3つの領域に分割し、それぞれの領域毎にゴルフボールモデル表面に接する格子の高さを設定する。この飛行行程の領域区分及び格子高さの設定については、後で詳述する。次に、上記ゴルフボールモデルにD重量を設定し、E初期条件初期速度、打出し角度、回転数を与え、上記気流仮想空間フィールドでゴルフボールモデルを運動させる。ここで、上記のように、気流仮想空間フィールドの下面に地面を想定した壁面を設定して、該気流仮想空間フィールドの下面を閉鎖面とした場合には、この下面からゴルフボール最下点までの距離を、打撃高さとして初期条件の1つに加えることもできる。これにより、ボールをティーに載せて打撃した場合を想定したシミュレーションをより高精度に行うことができる。このようにして、F上記Dの重量を有する上記ゴルフボールモデルが、上記気流仮想空間内フィールドで上記Eの初期条件による運動を開始し、ボールモデルの移動に追従するように、ボール近傍の格子は元に設定した格子サイズを崩さずに移動部で再度生成し、ボールが通過した後の格子は気流仮想空間内フィールドの格子に戻す操作を連続的に繰り返す。そして、ゴルフボールモデルが所定回転数で回転しながら、気流仮想空間内フィールドを飛行する際、この回転しているゴルフボールモデルに気流が当たることにより生じるゴルフボールモデル表面に作用する力によって発生する気体流の気流仮想空間内フィールドの運動要素の解析を格子区画毎に行う。ゴルフボールモデルが飛行することにより生じる気流がゴルフボールモデルに当たったときに生じる運動要素は、三次元空間座標系の各軸方向の気体流の速度、気体流の方向及びボールモデル表面に対する気体流の圧力であり、計算に当たって用いる基礎方程式、即ち、以下に示される質量保存則に相当する連続の式1~3及び物体の運動保存則に相当するナビエストークの方程式4~6に数値を代入して算出することができる。上記シミュレーションでは気流仮想空間の格子区画毎に空気の流れを演算により解析することができる。この演算に上記式1~6を用い、気流仮想空間を格子区画に区切ったことに対応して上記式1~6を離散化して演算を行うことができる。シミュレーションの方法は、シミュレーションの条件等を考慮して有限差分法、有限体積法、境界要素法、有限要素法等を適宜選択して行うことができる。ここで、場合によっては、D-2上記気流仮想空間内においてゴルフボールモデルに所定速度の気流が流入する状態を設定し、風の状態に応じたシミュレーションを行うこともできる。即ち、気流仮想空間で横から、進行方向、又は逆方向の気流を作成し、フォローの風のとき、アゲンストの風のとき、横風のときの状態を考慮してシミュレーションすることもできる。具体的には、例えばアゲンストの風を考慮する場合は、気流仮想空間にゴルフボールモデル前方から所定速度の気流が流入する状態を設定し、気流仮想空間全体の各格子に上記気流が設定された状態で、ゴルフボールモデルを運動させる。アゲンストの風のとき以外の風の状態を考慮する場合は、上記流入する気流の方向を変える事でシミュレーションを行うことができる。次に、G気流仮想空間内フィールドの各々の格子区画において算出された気体流の速度、方向及び圧力を積算して、上記気流仮想空間内フィールドにおける飛行中のゴルフボールの揚力係数及び抗力係数を算出する。この場合、気体流を散乱モデルとして乱流度を加味した速度値を代入して揚力係数CL及び抗力係数CDを算出することもできる。そして、上記揚力係数と抗力係数、上記Dで設定したボール重量、及びEの初期条件初期速度、打出し角度、回転数、打撃高さからボールモデルを運動させて、H飛行中のゴルフボールの高さ変化、横方向変化、速度変化、回転数変化を算出し、打出時から落下時までのゴルフボールの飛距離、左右へのばらつきを算出して、弾道を推定する。クラブにより打撃されて飛行するゴルフボールは、図4に示したように、重力Mg、空気による抵抗抗力D、更にボールがスピンを有するためにマグナス効果による揚力Lを受けることが知らされている。なお、図中、Vは飛行方向、cはボール中心を示し、この場合、ゴルフボールbは矢印R方向に回転している。この場合、ゴルフボールに働く力は下記弾道方程式7で表される。F=FLFDMg7F:ゴルフボールに働く力FL:揚力FD:抗力Mg:重力そして、以下の運動方程式により、ゴルフボールの運動を演算することができ、Fcd=0.5×CD×ρ×A×V28Fcl=0.5×CL×ρ×A×V29FcdFclmg=m×dV/dt10m:ボール重量g:重力加速度t:時間CL:揚力係数CD:抗力係数ρ:空気密度A:ゴルフボールの断面積V:ゴルフボールの対空気速度これから、微少時間飛行した後のゴルフボールの速度が算出され、更に、その結果を放物運動に関するニュートンNewtonの運動方程式に適用することにより、微少時間飛行した後のボールの飛距離及び高さ変化を算出することができる。上記微分方程式の演算は、オイラー法やルンゲクッタ法を使用することができ、これにより微少時間に対する速度の算出が可能である。一方、飛行するゴルフボールの回転数回転速度は時間の経過によって次第に減衰するが、微少時間経過後の回転数は次の式によって計算することができる。ω=ω0×Exp{-SRD1SRD2×V×t×β}11β=π×ρ×r4/l12ω0:ゴルフボールの初期回転角速度ω:ゴルフボールの回転角速度V:ゴルフボールの対空気速度t:時間ρ:空気密度r:ゴルフボールの半径l:ゴルフボールの慣性モーメントSRD1,SRD2:実験により求められたディンプル固有の減衰係数この結果、微少時間飛行した後のボールの飛距離、高さ変化が算出され、これを第1ステップとし、更に、算出された回転数及び速度を上記微少時間経過後の回転数及び速度として、上記操作を繰り返すことにより、ボール打出時から落下時までの各々の微少時間経過毎のボール飛距離及び高さを順次算出して弾道を推定することができる。なお、設定する微少時間は、シミュレーションの効率と精度の点から0.001~0.5秒とすることが好ましい。ここで、本発明では、上記操作により微少時間経過毎のボール飛距離及び高さを算出する際、上記ゴルフボールモデル周囲のレイノルズ数Reに応じて区分した高速領域、中速領域及び低速領域の領域毎にゴルフボールモデル表面と接する格子の高さを設定する。まず、飛行する上記ゴルフボールモデル周囲のレイノルズ数Reを算出すると共に、そのレイノルズ数Reが250000~170000の飛行行程を高速領域、Reが170000~100000の飛行行程を中速領域、Reが100000以下の飛行行程を低速領域とする基準に従って全飛行行程を2つ又は3つの領域に分ける。つまり、上記飛行中のゴルフボールモデル周囲のレイノルズ数Reを計算して、打出し時からレイノルズ数Reが170000に低下するまでの飛行行程を高速領域とし、Re=170000からRe=100000までの飛行行程を中速領域とし、Re=100000以下に低下してから着地までの飛行行程を低速領域とする。この場合、上記高速領域と中速領域との境界となるRe=170000は、高速領域内としても中速領域内としてもよいが、時間短縮を優先する場合には中速領域とすることが好ましく、精度を優先する場合は、高速領域とすることが好ましい。同様に、中速領域と低速領域との境界となるRe=100000も、中速領域内としても低速領域内としてもよいが、時間短縮を優先する場合には低速領域とすることが好ましく、精度を優先する場合は中速領域とすることが好ましい。なお、アイアンショットなどの打出し直後のレイノルズ数Reが170000に満たない場合を想定するときには、上記の基準に従って全飛行行程は中速領域と低速領域の2つに分けられる。上記レイノルズ数の算出は、下記式13により算出される。Re=UL/v13Re:レイノルズ数U:特性速度ゴルフボールの対空気速度[ms]L:特性長さゴルフボールの直径[m]v:動粘度[m2s]そして、本発明では、上記レイノルズ数Reで区分した高速領域におけるゴルフボールモデル表面と接する格子の高さ格子高さA、中速領域におけるゴルフボールモデル表面と接する格子の高さ格子高さB及び低速領域におけるゴルフボールモデル表面と接する格子の高さ格子高さCを、格子高さA<格子高さB<格子高さCとなるように設定する。勿論、打出し直後のレイノルズ数がRe=170000に達せず、飛行行程を中速領域と低速領域の2つの領域に分ける場合にも、同様に格子高さB<格子高さCと設定する。これにより、ゴルフボールモデル周囲の格子を比較的大きくとって格子区画の数を少なくしても実用上問題のない精度を担保し得る領域ほど格子区画の数を少なくすることになり、精度を損なわない範囲で計算回数を少なくし、計算時間の短縮を図ることができるものである。ここで、上記格子高さA~Cのそれぞれの高さは、想定するゴルフボールモデルの初速度や打出し角度、初期回転数などに応じて十分な有用性を有する精度と優位な時間短縮が得られるように適宜設定すればよく、特に制限されるものではないが、上記高速領域における上記格子高さAは0.001~0.015mmの範囲、中速領域における上記格子高さBは0.005~0.02mmの範囲、低速領域における上記格子高さCは0.01~0.025mmの範囲でそれぞれ設定することが好ましい。この場合、上記格子高さ範囲の上限を超えると十分な精度が得られずに良好な評価を行うことができなくなる場合があり、一方下限を下回ると十分な時間短縮を図ることができない場合がある。本発明の弾道シミュレーションを実施する際のより具体的手順を示せば、例えば次の通りである。まず、図1Aのように、ボールが打ち出される気流仮想空間をコンピュータに設定すると共に、この気流仮想空間内フィルードに格子を生成する。この場合、必要に応じて気流仮想空間の下面に地面を想定した壁面を設定する。一方、例えば図2に示したゴルフボールモデルにつき、ボール重量、直径、慣性モーメントを設定すると共に、打出し初期条件として、初速度、打出し角度、初期回転数を設定し、更に必要に応じて打撃高さを設定する。このゴルフボールモデルを図1Aのように格子を生成させた気流仮想空間内フィールド内に設定し、かつ図1Bに示した細かい格子をゴルフボールモデル近傍に生成させると共にボールモデルから離れる方向に格子のサイズが漸次増大するように区画する。この場合、ゴルフボールモデル表面に接する格子の高さは、後述する高速領域に設定された格子高さAとされる。そして、このボールモデルを上記初期条件で気流仮想空間内フィールド内を運動させ、揚力係数CL及び抗力係数CDを計算し、次いで上記式8,9及び10を用いて所定の微少時間飛行した後のゴルフボールの速度を算出し、更にその結果から、放物運動に関するニュートンNewtonの運動方程式を用いて、上記微少時間飛行した後のボールの飛距離及び高さを算出し、一方で式11及び12を用いて所定時間t経過後の回転数を算出する。また、同時に上記13を用いてボールモデル周囲のレイノルズ数Reを算出する。次に、この算出された微少時間飛行した後のゴルフボールの速度及び回転数から、更に微少時間飛行した後のボールの飛距離、左右のばらつき、高さ、回転数を上記した方法で同様に算出し、これを繰り返す。そして、上記レイノルズ数Reが170000に低下した飛行点からボールが上記中速領域に達したと判断し、上記ゴルフボールモデル表面と接する格子の高さを中速領域に設定された所定の格子の高さ格子高さBへと高くする。更に、この格子高さBで上記計算を繰り返し、上記レイノルズ数Reが100000に低下した飛行点からボールが上記低速領域に達したと判断し、上記ゴルフボールモデル表面と接する格子の高さを低速領域に設定された所定の格子の高さ格子高さCへと更に高くし、ゴルフボールが落下する位置まで、この格子高さCを用いて同様に上記計算を繰り返し、図3に示したゴルフボールの弾道を得ることができる。なお、上述したように、初期のレイノルズ数Reが170000に満たない場合には、飛行行程を中速領域と低速領域との2領域に分けて同様に計算し、ゴルフボールの弾道を得る。また、必要に応じて、算出されるゴルフボール落下時の速度、角度、回転数から、ゴルフボールが落下した後ゴルフボールが静止するまでの距離ラン、転がり距離を算出し、打出しから静止するまでの飛距離を算出することもできる。ラン転がり距離は、下記式14により計算することができる。ラン転がり距離=A×LandXvelocity-B×LandYvelocityC×LandAngle-D×LandSpin14LandXvelocity水平落下速度成分m/sLandYvelocitySin2θ×垂直落下速度成分m/sLandAngleCosθLandSpin落下スピン量A,B,C,D:実際のゴルフボールの試験から逆算して得られる固有の係数である。これらの係数を変化させることで、地面の状態の違い例えば、芝とコンクリートもシミュレーションに反映させることができる。本発明の弾道シミュレーション方法では、風洞試験のような実物モデルによる試験評価を実施することなく、ゴルフボールを任意の初速度、任意の回転数及び任意の角度で打出したときのゴルフボールの弾道を推定することができ、更にフォロー、アゲンスト、横風などの風の状態に応じた弾道を推定することも可能である。しかも本発明によれば、レイノルズ数Reによってボールの全飛行行程を2つ又は3つの領域に分け、それぞれの領域毎にゴルフボールモデル表面に接する格子の高さを設定して格子区画の大きさ及び数を調整することにより、十分な精度を担保しつつ計算時間の短縮を図ることができるものである。以下、実施例及び比較例を示し、本発明を具体的に説明するが、本発明は下記実施例に制限されるものではない。実施例、比較例図2に示したゴルフボールモデルにつき、ボール重量45.3g、直径42.7mm、慣性モーメント7.8×106kgm2と設定すると共に、打出し初期条件として、初速度65m/s、打出し角度10.5°、初期回転数2800rpm、打撃高さ50mmを設定する。このゴルフボールモデルを図1Aのように格子を生成させた気流仮想空間内フィールド内に設定し、かつ図1Bに示した細かい格子をゴルフボールモデル近傍に生成させると共にボールモデルから離れる方向に格子のサイズが漸次増大するように区画する。そして、このボールモデルを上記初期条件で気流仮想空間内フィールド内で運動させ、揚力係数CL及び抗力係数CDを計算し、次いで上記式8,9及び10を用いて所定の微少時間0.001秒飛行した後のゴルフボールの速度を算出し、更にその結果から、放物運動に関するニュートンNewtonの運動方程式を用いて、上記微少時間0.001秒飛行した後のボールの飛距離及び高さを算出し、また一方で式11及び12を用いて上記微少時間0.001秒経過後の回転数を算出する。次に、この算出された微少時間0.001秒飛行した後のゴルフボールの速度及び回転数から、更に微少時間0.001秒飛行した後のボールの飛距離、左右のばらつき、高さ、回転数を上記した方法で同様に算出し、その後も同様にして、ゴルフボールが落下する位置まで、上記計算を繰り返し、ゴルフボールの弾道をシミュレートして、飛距離キャリーを推定した。このとき、実施例1及び実施例2のシミュレーションでは、飛行中のゴルフボール周囲のレイノルズ数Reを上記式13を用いて算出し、全飛行行程をRe=250000~170000の高速領域、Re=170000~100000の中速領域、及びRe=100000以下の低速領域の3つの領域に分け、各領域におけるゴルフボールモデル表面と接する格子の高さを表1の通りに設定して、上記シミュレーションを行った。一方、比較例にあっては領域の区分を行なわずに全飛行行程につき格子高さを0.001mmに設定してシミュレーションを行なった。推定された飛距離とシミュレートに要した時間計算時間を表1に示す。また、計算時間については、下記基準で時間短縮の効果を評価した。シミュレーションに用いたコンピュータの詳細は以下の通りである。なお、ボールモデル表面に接する格子区画は直交格子立方体とし、そのサイズは高さ表1記載×横0.5mm×厚さ0.5mmとした。計算時間についての評価基準◎150時間以内○250時間以内×250時間よりも長いなお、下記コンピュータを用いた今回のシミュレーションでは、このように評価基準を設定したが、計算時間の評価基準はコンピュータの性能により適宜変更する必要がある。コンピュータの詳細OS:Windows登録商標XPCPU:intelXeonX5482320GHzL2=6MB×2FSB=1600MHz2CPU8Coreメモリ:64GB4096MB×16DDR2800ECCREGFBDIMM一方、上記シミュレーションで設定したゴルフボールモデルと同様のゴルフボールにつき、下記の方法で打撃試験を行なって飛距離を測定し、得られた実測値213.2ヤードと上記シミュレーションによる結果とを比較して、下記基準により精度の評価を行った。結果を表1に示す。打撃試験打撃ロボットにドライバーW1[ブリヂストンスポーツ社製TourStageX-Drive701ロフト9°]を装着し、実際に打撃して飛距離を測定して実測値とした。このとき、打撃条件は、ボール初速が65m/s、打出し角度が10.5°、初期バックスピン量が2800rpmとなるように設定した。また、ボールをティーに載せて打撃を行い、地面からボール最下部までの高さを50mmとした。結果は、飛距離キャリー213.2ヤード、トータル距離キャリーラン227.5ヤードであった。精度の評価基準◎実測値の±5ヤード以内○実測値の±10ヤード以内×実測値の±10ヤードより大表1の通り、全飛行行程を所定のレイノルズ数Reにより3つの領域に分け、Re=250000~170000の高速領域におけるゴルフボールモデル表面と接する格子の高さ格子高さA、Re=170000~100000の中速領域におけるゴルフボールモデル表面と接する格子の高さ格子高さB、及びRe=100000以下の低速領域におけるゴルフボールモデル表面と接する格子の高さ格子高さCを格子高さA<格子高さB<格子高さCとなるように設定することにより、十分な精度を維持しつつ計算時間を効果的に短縮することができた。"