被冷却体と間隔を開けて配置された複数のノズルから冷却媒体を供給することにより前記被冷却体を冷却する冷却装置で、少なくとも1の前記ノズルから供給される前記冷却媒体が前記被冷却体に当たる領域を単位解析領域とし、該単位解析領域を平均化した数値解析を行って前記被冷却体の温度を予測する際に、前記冷却媒体が、前記ノズルの先端よりも前記被冷却体側、且つ、前記被冷却体の表面よりも前記ノズル側から供給されていると仮定し、前記冷却媒体の少なくとも質量および運動量をソース項として取り扱い、前記質量のソース項を質量保存式へ、前記運動量のソース項を運動量保存式へと加えて、少なくとも質量場、速度場、および、圧力場を計算する数値解析を行って、前記被冷却材の温度を予測することを特徴とする、被冷却体の温度予測方法。さらに、前記冷却媒体のエネルギーをソース項として取り扱い、前記エネルギーのソース項をエネルギー保存方程式へと加えて、少なくとも質量場、速度場、圧力場、および、温度場を計算する数値解析を行って、前記被冷却体の温度を予測することを特徴とする、請求項1に記載の被冷却体の温度予測方法。前記被冷却体および前記冷却媒体の間の熱流束と、前記被冷却体の表面温度と、前記被冷却体表面近傍の冷却媒体温度と、前記被冷却体の表面性状との関係式を予め導出し、前記単位解析領域について、前記関係式を用いて熱流束を算出し伝熱解析を行うことで、各単位解析領域の前記被冷却体の温度を予測することを特徴とする、請求項1又は2に記載の被冷却体の温度予測方法。請求項1~3のいずれか1項に記載の被冷却体の温度予測方法を用いて、冷却装置により冷却される被冷却体の温度分布を予測する温度分布予測工程と、予測された前記被冷却体の温度分布を用いて、前記冷却装置の動作を制御する冷却制御工程と、を有することを特徴とする、金属板の製造方法。特許文献1に開示されている手法では、冷却媒体がどのように流動するか、また冷却媒体が液体である場合、被冷却体周囲にどの程度溜まるかを予測する手法が不明である。冷却媒体は被冷却体から熱を受けながら、すなわち冷却媒体自身の温度を高くしながら、被冷却体周囲および冷却装置内から流出していく。そのため、冷却媒体の流動状態が分からなければ、冷却媒体の温度分布は分からず、したがって、冷却媒体の温度に影響を受ける被冷却体の冷却状況を正確に予測することはできない。また、非特許文献1に開示されている技術では、鋼板のごく局所の解析に100万点程度の解析格子を要しており、当該手法を冷却装置全体といった広範囲な領域に適用し、工学的に実用性のある解析を行うことは、解析負荷の点から非現実的である。そこで本発明は、解析負荷を抑制しつつ簡便的に被冷却体、すなわち高温の金属板の温度分布の予測することが可能な、被冷却体の温度予測方法およびこれを用いる金属板の製造方法を提供することを課題とする。数値解析は、科学分野の数学的問題を、数値の演算によって解く方法である。被冷却体の温度予測温度分布予測における数値解析では、支配方程式として少なくとも質量保存方程式および運動量保存方程式を含み、それらを離散化し、少なくとも質量場、圧力場、速度場を数値的に求める流動数値解析、またはこれらに加え、支配方程式としてエネルギー保存方程式をさらに含み、それらを離散化し、温度場を数値的に求める熱流動数値解析が行われる。流動数値解析または熱流動数値解析以下において、これらをまとめて熱流動数値解析という。を用いれば、冷却装置に複数備わるノズルから供給される、ラミナー流、噴流、スプレー流等といった冷却水冷却媒体ともいう。以下において、冷却媒体を冷却水ということがあるが、純水に限らず、工業用水であってもよく、油、その他の液体、およびそれらに不可避的に混入する異物が含まれていてもよい。の流動状態や温度分布、また被冷却体の温度分布を、その流動を解像可能な程度に詳細な解析格子を用いて詳細に知ることが可能となる。ところが、このような熱流動数値解析を用いて、例えば鉄鋼製造プロセスの冷却装置全体といった、長さ数十mにもおよぶ広範囲な領域の熱流動状態を詳細に予測することは、解析負荷が過大になるため非現実的である。金属板の製造で用いられる、複数のノズルを備えた通常の冷却装置は、個々のノズルからのラミナー流、噴流、スプレー流等よりも非常に大きい。したがって、個々のノズルから供給される冷却水は装置全体の熱流動状態に直接は影響しないと考えられ、それらの簡略化、モデル化が可能と考えられる。本発明者らは、その方法として、ノズルから供給される個々のラミナー流、噴流、スプレー流等の冷却水を、その熱流動を精緻に解析できるほどに詳細な解析格子では解像せず、一つまたは複数のラミナー流、噴流、スプレー流等を含む領域を単位解析領域とし、これを平均化して取り扱う方法について検討した。このような方法を用いることにより、解析格子数を減少させることが可能となり、現実的な解析負荷で冷却装置全体の熱流動数値解析が可能となる。その結果、これまで困難であった冷却媒体の流動状態およびその温度分布予測が可能となり、被冷却体の温度予測が可能となる。単位解析領域を平均化して温度予測するには、平均化する領域単位解析領域から供給される冷却水の流量、流速、温度が分かればよいことが一般的に知られている。鉄鋼製造プロセスにおける冷却装置のように供給する冷却水の水量密度が大きい場合、被冷却体周囲に、排水されるまでの冷却水が溜まる。このような冷却水は、滞留水と呼ばれる。冷却装置内は、気体空気や蒸気と液体冷却水とが混合する流動現象となっており、滞留水の流動状態もその相互作用で決定される。したがって、滞留水の流動状態を得るためには、少なくとも気液二相流を取り扱える数値解析法が必要となる。本発明者らは、検討の結果、数値解析において平均化したラミナー流、噴流、スプレー流等の冷却水を与えるとき、実際にノズルが設置されている場所に与えるより、被冷却体表面近傍に与えた方例えば、冷却媒体が、ノズル先端より下方かつ被冷却体の上方から供給されていると仮定して数値解析することが、より正確な予測が可能であることを知見した。その理由は、実際にノズルが設置されている場所に冷却水を与えると、平均化して与えた冷却水の動きが滞留水に阻害されるためである。個々のノズルからラミナー流、噴流、スプレー流等として供給された冷却水は、通常、滞留水を貫通しつつ被冷却体表面へと達する。ところが数値解析において、例えば滞留水が多くなる条件で解析を行った場合、平均化した冷却水を供給すると滞留水を貫通できず、滞留水にまぎれてしまう場合があった。個々のノズルからのラミナー流、噴流、スプレー流等の直下における被冷却体表面近傍では、冷却水が供給され続けることで水温が低く保たれ、冷却が進む。ところが熱流動数値解析において、平均化することで供給された冷却水が滞留水を貫通できないと、被冷却体表面近傍の冷却水温度が上昇して冷却が進まなくなり、実装置とは異なる状況に陥る。そこで数値解析において、ノズルが設置されている場所からではなく、被冷却体表面近傍の空間から平均化された冷却水を供給すること、すなわち、数値流体解析的にいえば、ソース湧き出し項により冷却水を系解析する空間に供給することで、このような問題がなく、数値解析できることがわかった。実際の冷却装置では、個々のノズルから供給される冷却水は質量、運動量、エネルギー温度を持っている。そのため、本発明による数値解析においてもこれらをソース項として、それぞれ質量保存方程式、運動量保存方程式、エネルギー保存方程式に組み込めばよい。ただし、滞留している冷却水の量や冷却水の流速を知ることができれば、経験則で温度や温度ムラを推測できる場合もある。計算の負荷を低減するという観点からは、特にエネルギー温度をソース項として考慮せず、冷却水の数値解析のみを行ってもよい。以上のような数値計算を行えば、十分に精度の高い数値解析ができ、結果として精度の高い被冷却体の温度予測、ひいては温度分布予測が可能となる。より精度の高い数値解析を行うには、被冷却体-冷却水間の運動量の相互作用を考慮してもよい。鉄鋼製造プロセスに見られる高温鋼材の冷却を取り扱う場合は、鋼材近傍の冷却水は沸騰状態にあるとして取り扱う。本発明は、以上の知見に基づいてされたものであって、本発明の趣旨は以下のとおりである。以下、本発明について説明する。なお、本発明の理解を容易にするため、添付図面の参照符号を括弧書きにて付記するが、それにより本発明が図示の形態に限定されるものではない。本発明の第1の態様は、被冷却体1と間隔を開けて配置された複数のノズルから冷却媒体3、3、...を供給することにより被冷却体を冷却する冷却装置で、少なくとも1のノズルから供給される冷却媒体が被冷却体に当たる領域を単位解析領域とし、該単位解析領域を平均化した数値解析を行って被冷却体の温度を予測する際に、冷却媒体が、ノズルの先端よりも被冷却体側、且つ、被冷却体の表面よりもノズル側から供給されていると仮定し、冷却媒体の少なくとも質量および運動量をソース項として取り扱い、質量のソース項を質量保存式へ、運動量のソース項を運動量保存式へと加えて、少なくとも質量場、速度場、および、圧力場を計算する数値解析を行って、被冷却体の温度を予測することを特徴とする、被冷却体の温度予測方法である。また、上記本発明の第1の態様において、さらに、冷却媒体3、3、...のエネルギーをソース項として取り扱い、該エネルギーのソース項をエネルギー保存方程式へと加えて、少なくとも質量場、速度場、圧力場、および、温度場を計算する数値解析を行って、被冷却体1の温度を予測してもよい。また、上記本発明の第1の態様において、被冷却体1および冷却媒体3、3、...の間の熱流束と、被冷却体の表面温度と、被冷却体表面近傍の冷却媒体温度と、被冷却体の表面性状との関係式を予め導出し、上記単位解析領域について、上記関係式を用いて熱流束を算出し伝熱解析を行うことで、各単位解析領域の被冷却体の温度を予測してもよい。本発明の第2の態様は、上記本発明の第1の態様にかかる被冷却体の温度予測方法を用いて、冷却装置により冷却される被冷却体1の温度分布を予測する温度分布予測工程と、予測された被冷却体の温度分布を用いて、冷却装置の動作を制御する冷却制御工程と、を有することを特徴とする、金属板の製造方法である。本発明によれば、多大な解析時間や費用を要さず、冷却装置内全体の冷却水流動状況を予測することが可能になるため、被冷却体1と冷却水との間の熱流束に大きな影響を及ぼす冷却水の温度分布や冷却水の流動状態を求めることが可能になる。それにより、冷却装置の設計をより高い信頼性をもって行うことが可能となる。また、稼働中の冷却装置の内部で生じている現象を推測することが可能になるので、推測結果を冷却装置の改良に活かすことが可能となる。特に、鉄鋼冷却プロセスにおいては冷却装置に鋼材が挿入される段階で、鋼材に温度分布が生じている場合がほとんどであるが、本発明を適用することで、鋼板温度分布を低減させ、冷却停止温度を均一化させる冷却装置の制御方法および改良方法を提案することや、冷却停止温度を均一化させて品質を均一化させた金属板の製造方法を提供することが可能になり、鋼材品質の均一化や熱応力等による変形を抑制することも可能になる。図1は、被冷却体表面温度、および、被冷却体と冷却媒体との間の熱流束熱輸送量を表す一般的な図である。図1に示したように、一般に沸騰現象は、冷却媒体に接する物質被冷却体の温度が上昇するにしたがって、核沸騰、遷移沸騰、膜沸騰と、沸騰の形態が変化する。ここで、核沸騰は、蒸気気泡が固体表面から立ち上る状態であり、蒸気気泡が発生する場所以外では固体および冷却媒体は接する状態となる。膜沸騰は、活発な沸騰が生じた結果固体と冷却媒体との間に蒸気の膜ができた状態であり、ほとんどの場所で固体および冷却媒体は直接には接していない状態となる。一方、遷移沸騰は核沸騰および膜沸騰が入り混じった状態とされる。数値解析を行うに当たり、沸騰形態と被冷却体-冷却媒体間の運動量の境界条件との関係を考えると、核沸騰時は、被冷却体と冷却媒体とがほとんどの部分で接しているため、熱流動数値解析で一般的に用いられる滑りなし壁条件、または乱流状態では壁関数を用いた壁面境界条件を与えればよい。膜沸騰時は、被冷却体と冷却媒体との間に摩擦が小さいと考えれば、熱流動数値解析で一般的に用いられる滑り壁境界条件、または蒸気の膜が持つ粘性によるせん断応力を壁面境界条件として与えればよい。遷移沸騰時は、核沸騰に近いか膜沸騰に近いかを判別し、これらの中間の値を与えればよい。なお、被冷却体表面の状態が非沸騰、核沸騰、遷移沸騰、膜沸騰のどの状態であるかも、被冷却体表面温度、被冷却体表面近傍の冷却水温度、滞留水の流動状態といったパラメータにより分類できる。したがって、数値解析により前記パラメータを求めれば、沸騰の状態を判別可能である。また、冷却装置で一般に用いられる冷却水の温度は、室温程度である、すなわち沸点よりも十分に低い。このような場合、高温の被冷却体表面において冷却水は沸騰するものの、生じた蒸気は滞留水により冷却され、凝縮し、水に戻る。したがって、熱流動数値解析において、このような相変化現象は解析してもよく、解析しなくてもよい。数値解析で被冷却体の温度を予測する際には、被冷却体と冷却水との間の熱流束を考慮する必要がある。被冷却体と冷却水との間の熱流束は、被冷却体の表面温度、被冷却体の表面性状、被冷却体表面近傍の冷却水温度、および冷却水の流動状態をパラメータとして数式化することができる。このため、熱流束とこれらの関係式を予め導出し、当該関係式を用いて熱流束を算出し伝熱解析を行うことで各解析格子の被冷却体の温度を決定し、被冷却体と冷却水との間の熱流束を各時刻各位置において与え、被冷却体の温度を予測する。ここで、被冷却体の表面性状は、少なくとも被冷却体の表面粗さを含み、特に鉄鋼製造プロセスに見られる高温鋼材の冷却を取り扱う場合には、被冷却体である鋼材表面に生成するスケール酸化物の表面粗さ、厚み、物性、空隙率といった、スケールの物理的性状を含むものである。被冷却体の表面性状については、予め与えたものを使用してもよく、また、特に鉄鋼製造プロセスに見られる高温鋼材の冷却を取り扱う場合、例えば冷却中にスケール厚みが変化するといったモデルを組み込み、熱流動数値解析中に変化する値として与えてもよい。また、被冷却体表面近傍の冷却水温度の被冷却体表面近傍とは、流体力学において一般的に知られているように、冷却媒体が被冷却体表面に形成する温度境界層の外縁付近かそれよりも遠方を意味する。また、冷却水の流動状態とは、例えば冷却媒体の攪拌状態が、噴流やスプレー流が被冷却体に直接衝突する領域と、噴射された冷却媒体が滞留水となって流動する領域とでは大きく異なり、このような差異を意味する。他にも滞留水の流速による差異や乱流状態も流動状態として考慮すべきである。被冷却体の表面の各単位解析領域における温度を求めれば、被冷却体全体の温度分布を予測することができる。以上の計算を時間進展させつつ繰り返すことで、冷却装置中の冷却水の流動状態や温度分布を、現実的な解析負荷範囲で精度および効率良く推定することができ、さらに被冷却体の温度ひいては温度分布を精度および効率良く推測することが可能となる。以下に、複数の柱状噴流により鋼板を冷却する装置に対して、本発明の被冷却体の温度予測方法を適用した例を示す。図2に、数値解析を行う際の被冷却体表面近傍における解析格子の例を示す。図2aは、本発明を用いずに噴流を解像できる解析格子を用いて解析を行う場合を説明する概念図であり、図2bは、本発明を用いて噴流を解像しない解析格子を用いて解析を行う場合を説明する概念図である。図2a及び図2bの各解析格子には数値計算時の解析イメージを合わせて着色を施した。図2a及び図2bの紙面下方が鋼板被冷却体であり、鋼板における着色は黒いほど低温であることを示す。また、図2a及び図2bの紙面上方は冷却水が供給される領域であり、着色が濃いほど水の濃度が高いことを示す。なお、図2bでは冷却媒体がノズル先端より下方かつ被冷却体の上方から供給されることをイメージしやすいように冷却水の湧き出しと表記し、1の解析格子を単一濃度にせず濃淡があるように表示したが、実際の計算では、冷却水の質量等は解析格子ごとに平均化して計算される。図2aは詳細な解析格子メッシュとした例であり、この例の数値解析では、被冷却体の上方に位置する冷却媒体の供給ノズルからの冷却水を考慮して数値解析がされる。この場合、柱状噴流を解析格子で解像しないと数値解析が行えず、解析格子数が膨大になるため、計算負荷が大きくなり冷却装置全体といった広範囲な領域を解析することは困難である。一方、図2bは解析格子を粗くした例であり、解析格子数を低減し、計算負荷を小さくすることで現実的な解析が可能となる。また、この例ではノズルからの柱状噴流を考慮しなくてよく、この分も計算負荷の軽減に寄与する。続いて、実際の試験結果と、図2bに示したような解析格子に本発明を適用した数値解析結果との比較について説明する。図3に、実際の試験に用いた装置の例を示す。本装置は、柱状噴流を被冷却体へ供給可能な多数のノズルを有する冷却ヘッダを具備し、該冷却ヘッダよりも幅が広く側方の一面が開いた水槽に冷却水が噴射できるようになっている。噴射された冷却水は滞留水となって水槽内に溜まり、水槽側方の開いた一面より流出する。その流出流量分布を測定した。一方、この装置を想定した解析領域について、流動数値解析を行った。図4に、図3に示した装置を用いた試験結果と、本発明を用いた解析結果との比較を示す。図4bは図4aに示す例の2倍の冷却水を用いた例である。図4a及び図4bでは、冷却ヘッダ中央からの距離を横軸とし、冷却水の質量流量を縦軸とし、縦軸は冷却水が流出する出口の平均流量により無次元化し、横軸は冷却水が流出する出口の水槽の幅により無次元化した。図4a及び図4bに示したように、異なる流量に対しても試験結果と解析結果とがよく一致している。この結果から、本発明における数値解析は滞留水の流動を妥当に表現可能であることが分かる。本発明者らは、本発明の妥当性、有用性についてさらに検証するため、同一の解析対象に対する金属板の温度予測について、本発明を用いて行った結果と、本発明を用いず、詳細な熱流動解析を用いて行った結果とを比較した。図5に、計算のための解析領域10を示す。ここでは、12本のノズルから柱状噴流3、3、...が厚鋼板1に噴射される領域を仮定した。厚鋼板1の上面には滞留水2が存在している。本発明の温度予測方法による計算では、この解析領域10を水平方向に20×20の解析格子に分割してその1を単位解析領域としたのに対し、詳細計算の温度予測方法では、この解析領域をさらに細かい解析格子に分割し、その解析格子を単位解析領域として計算を行った。図6に、本発明の温度予測方法で求めた温度分布結果を示す。図6aは、厚鋼板の上面に滞留した滞留水の流動も含めた、計算結果の全体を示す図であり、図6bは、厚鋼板の表面温度分布を示す図である。また、図7に、詳細に計算した温度予測方法で求めた温度分布結果を示す。図7aは、厚鋼板の上面に滞留した滞留水の流動も含めた、計算結果の全体を示す図であり、図7bは、厚鋼板の表面温度分布を示す図である。図7aでは、滞留水と噴射される柱状噴流が滞留水を貫通している様子も示されている。なお、本発明の温度予測方法で温度分布を求める際には、冷却水が、ノズルの先端よりも下方、且つ、厚鋼板の表面よりも上方から供給されていると仮定し、冷却水の少なくとも質量および運動量をソース項として取り扱い、質量のソース項を質量保存式へ、運動量のソース項を運動量保存式へと加えて、質量場、速度場、および、圧力場を計算する数値解析を行って、厚鋼板の温度を予測した。なお、本発明の温度予測方法で温度分布を求める際には、厚鋼板および冷却水間の熱流束と、厚鋼板の表面温度と、厚鋼板表面近傍の冷却水温度と、厚鋼板の表面性状との関係式を予め導出し、この関係式を用いて熱流束を算出し伝熱解析を行うことで、温度を導出した。図8に、表面粗さなし、鋼材表面近傍の冷却媒体の平均温度が37°Cおよび56°Cの場合における、鋼材表面温度と鋼材-冷却水間の熱流束との関係を示す。図8aは鋼板表面近傍の平均水温が37°Cの場合の関係を示す図であり、図8bは鋼板表面近傍の平均水温が56°Cの場合の関係を示す図である。これらは伝熱学において沸騰曲線として一般的に知られる線図であり、被冷却体の温度を予測する上で必要な線図である。このような線図から導き出される多数の関係式を多数用意することで、温度の導出が可能になる。図6bおよび図7bを見る限り、両者はほぼ同じであり、柱状噴流を詳細に捉えなくても、本発明による温度予測方法で求めた表面温度分布は、詳細計算による温度予測方法で求めた表面温度分布と差異はない。図9に、図6および図7の数値解析における厚鋼板の平均温度の経時変化を示す。図9の縦軸は厚鋼板の平均温度[°C]であり、同横軸は冷却開始からの経過時間[s]である。図9に示したように、両者の結果はよく一致した。すなわち、本発明の温度予測方法でも詳細に計算した温度予測方法で求めた温度分布結果と同様の結果が得られ、本発明による温度予測方法の妥当性が示されたことがわかる。なお、本発明では、冷却装置に供給される個々のノズルからのラミナー流、噴流、スプレー流等の冷却水の質量、運動量、エネルギーをソース項として解析に導入し、単位解析領域を平均化して温度予測する。これは、以下の理由による。まず熱流動数値解析において、冷却水の質量、運動量、エネルギーを系内に導入するには、境界条件として導入するか、ソース項として導入する必要がある。また本発明で扱うような金属板の製造プロセス例えば、鉄鋼製造プロセスで用いられる冷却装置は、長さが数十mにもおよぶものが多く、ノズルを一本ずつ解像できるような細かい解析格子を用いたのでは、冷却装置全体の流動状態や温度分布を、現実的な計算負荷で解析することは困難である。そこで個々のノズルからの流動を平均化する必要がある。ここで、ノズルから供給されるラミナー流、噴流、スプレー流を平均化して与える場合、実際の冷却装置におけるノズル先端位置から境界条件として導入すると、その挙動を捉えることができず、妥当な結果が得られない。したがって、本発明では、これらをソース項として導入する必要がある。また、本発明に関する上記説明では、厚鋼板および冷却水間の熱流束と、厚鋼板の表面温度と、厚鋼板表面近傍の冷却水温度と、厚鋼板の表面性状との関係式を予め導出し、この関係式図8に示したような沸騰曲線から導出される関係式を用いて熱流束を算出し伝熱解析を行うことで温度を予測する形態を例示したが、本発明は当該形態に限定されない。ただし、被冷却体から冷却媒体への熱流束は、被冷却体の表面温度、被冷却体の表面性状、被冷却体表面近傍の冷却媒体温度、および冷却媒体の流動状態と大きな相関がある。そのため、被冷却体から冷却媒体への熱流束と、被冷却体の表面温度、被冷却体の表面性状、被冷却体表面近傍の冷却媒体温度、および冷却媒体の流動状態との関係式を用いて数値解析を行うことにより、これらの被冷却体各所における分布を常に把握し、常に適切な熱流束を与えながら冷却の予測を行うことが可能になる結果、被冷却体の温度予測の精度を高めやすくなる。したがって、本発明では、沸騰曲線から導出される関係式を用いて数値解析を行うことが好ましい。また、本発明に関する上記説明では、被冷却体の上方に配置された複数のノズルから被冷却体の上面へと冷却媒体を供給することにより、被冷却体を冷却する形態について主に説明したが、本発明の形態は、これに限定されない。本発明は、被冷却体の下方に配置された複数のノズルから被冷却体の下面へと冷却媒体を供給することにより、被冷却体を冷却する形態であっても良い。このほか、本発明は、被冷却体の上方、および、被冷却体の下方にそれぞれ配置された複数のノズルから、被冷却体の上面、および、被冷却体の下面へと冷却媒体をそれぞれ供給することにより、被冷却体を冷却する形態であっても良い。以上、現時点において実践的であり、且つ好ましいと思われる実施形態に関連して本発明を説明したが、本発明は本願明細書中に開示された実施形態に限定されるものではなく、特許請求の範囲および明細書全体から読み取れる発明の要旨或いは思想に反しない範囲で適宜変更可能であり、そのような変更を伴う被冷却体の温度予測方法および金属板の製造方法も本発明の技術的範囲に包含されるものとして理解されなければならない。