"冷媒を圧縮する圧縮機と、圧縮された冷媒を熱源媒体によって凝縮させる凝縮器と、凝縮された冷媒を蒸発させると共に該冷媒と熱媒とを熱交換する蒸発器と、を備えた熱源機の熱媒の流量を推定する熱媒流量推定装置であって、前記圧縮機の吸込風量を反映した第1パラメータと、前記圧縮機のヘッドを反映した第2パラメータとで表示されたマップ上に、旋回失速となる旋回失速線、及び前記圧縮機が吸い込む冷媒の音速を示す複数の機械マッハ数線が示された空力特性マップを記憶した記憶手段と、前記第2パラメータを算出し、前記空力特性マップから該第2パラメータに応じた前記第1パラメータを算出する第1パラメータ算出手段と、前記第1パラメータ算出手段によって算出された前記第1パラメータに応じた前記圧縮機の吸込風量に基づいて、前記蒸発器において冷媒と熱媒との間で交換される熱量を算出し、該熱量に基づいて熱媒の流量を算出する熱媒流量算出手段と、を備えた熱媒流量推定装置。熱媒流量算出手段は、前記第1パラメータ算出手段によって算出された前記第1パラメータに基づいた前記圧縮機の吸込風量、及び前記圧縮機に吸い込まれる冷媒の密度から、前記蒸発器を流れる冷媒の流量を算出し、該算出した冷媒の流量、及び前記蒸発器の入口側のエンタルピと出口側のエンタルピの差から、前記蒸発器において冷媒と熱媒との間で交換される熱量を算出し、該算出した熱量、及び熱媒が前記蒸発器へ流入出する温度の差に基づいて、熱媒の流量を算出する請求項1記載の熱媒流量推定装置。前記圧縮機は、回転数の制御が可能とされており、前記記憶手段は、前記圧縮機の回転数に応じて異なる複数の前記空力特性マップを記憶し、第1パラメータ算出手段は、前記圧縮機の回転数に対応した前記空力特性マップから前記第2パラメータに応じた前記第1パラメータを算出する請求項1又は請求項2記載の熱媒流量推定装置。前記圧縮機は、冷媒入口に冷媒流量を調節するベーンが備えられており、前記記憶手段は、前記ベーンの開度に応じて異なる複数の前記空力特性マップを記憶し、第1パラメータ算出手段は、前記ベーンの開度に対応した前記空力特性マップから前記第2パラメータに応じた前記第1パラメータを算出する請求項1又は請求項2記載の熱媒流量推定装置。前記凝縮器と前記蒸発器の間には、前記凝縮器内にある冷媒を前記蒸発器へと流すためのバイパス配管が設けられると共に、該バイパス配管を流れる冷媒の流量を調整するための弁が設けられ、前記記憶手段は、前記弁の開度に応じて異なる複数の前記空力特性マップを記憶し、第1パラメータ算出手段は、前記弁の開度に対応した前記空力特性マップから前記第2パラメータに応じた前記第1パラメータを算出する請求項1又は請求項2記載の熱媒流量推定装置。冷媒を圧縮する圧縮機と、圧縮された冷媒を熱源媒体によって凝縮させる凝縮器と、凝縮された冷媒を蒸発させると共に該冷媒と熱媒とを熱交換する蒸発器と、請求項1から請求項5の何れか1項に記載の熱媒流量推定装置と、を備えた熱源機。冷媒を圧縮する圧縮機と、圧縮された冷媒を熱源媒体によって凝縮させる凝縮器と、凝縮された冷媒を蒸発させると共に該冷媒と熱媒とを熱交換する蒸発器と、を備えた熱源機の熱媒の流量を推定する熱媒流量推定方法であって、前記圧縮機の吸込風量を反映した第1パラメータと、前記圧縮機のヘッドを反映した第2パラメータとで表示されたマップ上に、旋回失速となる旋回失速線、及び前記圧縮機が吸い込む冷媒の音速を示す複数の機械マッハ数線が示された空力特性マップが記憶手段に予め記憶されており、前記第2パラメータを算出することによって、前記空力特性マップから該第2パラメータに応じた前記第1パラメータを算出する第1工程と、前記第1工程によって算出された前記第1パラメータに応じた前記圧縮機の吸込風量に基づいて、前記蒸発器において冷媒と熱媒との間で交換される熱量を算出し、該熱量に基づいて熱媒の流量を算出する第2工程と、を含む熱媒流量推定方法。しかしながら、特許文献1に記載の技術では、冷却水流量を算出するために、冷水流量を計測する流量計を用いている。また、特許文献2に記載の技術では、個々の空気調和機の冷温水の流量を計測するために、全体の冷温水の流量を計測する流量センサー及び複数の差圧センサーを用いている。上記のように、特許文献1,2に記載の技術では、所定の流体の流量を算出するために、他の流体の流量を計測する流量計を用いたり、他の流体の差圧を計測する差圧計を用いたりしているため、低コストで流体の流量を把握することができない。本発明は、このような事情に鑑みてなされたものであって、流量計を用いることなく、熱媒の流量を算出することができる熱媒流量推定装置、熱源機、及び熱媒流量推定方法を提供することを目的とする。上記課題を解決するために、本発明の熱媒流量推定装置、熱源機、及び熱媒流量推定方法は以下の手段を採用する。すなわち、本発明に係る熱媒流量推定装置は、冷媒を圧縮する圧縮機と、圧縮された冷媒を熱源媒体によって凝縮させる凝縮器と、凝縮された冷媒を蒸発させると共に該冷媒と熱媒とを熱交換する蒸発器と、を備えた熱源機の熱媒の流量を推定する熱媒流量推定装置であって、前記圧縮機の吸込風量を反映した第1パラメータと、前記圧縮機のヘッドを反映した第2パラメータとで表示されたマップ上に、旋回失速となる旋回失速線、及び前記圧縮機が吸い込む冷媒の音速を示す複数の機械マッハ数線が示された空力特性マップを記憶した記憶手段と、前記第2パラメータを算出し、前記空力特性マップから該第2パラメータに応じた前記第1パラメータを算出する第1パラメータ算出手段と、前記第1パラメータ算出手段によって算出された前記第1パラメータに応じた前記圧縮機の吸込風量に基づいて、前記蒸発器において冷媒と熱媒との間で交換される熱量を算出し、該熱量に基づいて熱媒の流量を算出する熱媒流量算出手段と、を備える。本発明によれば、熱媒流量推定装置は、冷媒を圧縮する圧縮機と、圧縮された冷媒を熱源媒体によって凝縮させる凝縮器と、を備えた熱源機の熱媒の流量を推定する装置である。熱媒流量推定装置が備える記憶手段には、圧縮機の吸込風量を反映した第1パラメータと、圧縮機のヘッドを反映した第2パラメータとで表示されたマップ上に、旋回失速となる旋回失速線、及び前記圧縮機が吸い込む冷媒の音速を示す複数の機械マッハ数線が示された空力特性マップが記憶されている。なお、空力特性マップは、予め圧縮機の運転試験を綿密に行うことによって作成されるものである。第2パラメータ及び機械マッハ数は、圧縮機の運転状態に対応した値であり、第1パラメータは、第2パラメータ及び機械マッハ数で特定できるため、第2パラメータ及び機械マッハ数圧縮機が吸い込む冷媒の音速を算出することで、第1パラメータ、すなわち圧縮機の吸込風量の算出が可能である。なお、第2パラメータ及び冷媒の音速は、蒸発器内の圧力や凝縮器内の圧力から算出できる。そこで、第1パラメータ算出手段によって、まず第2パラメータが算出され、空力特性マップから該第2パラメータに応じた第1パラメータが算出される。そして、熱媒流量算出手段によって、第1パラメータ算出手段で算出された第1パラメータに応じた圧縮機の吸込風量に基づいて、蒸発器において冷媒と熱媒との間で交換される熱量が算出され、該熱量に基づいて熱媒の流量が算出される。すなわち、熱媒流量算出手段によって、蒸発器における冷媒と熱媒との熱バランスにより熱媒の流量が算出される。このように、本発明は、空力特性マップに基づいて算出した圧縮機の吸込風量を用いて、蒸発器で交換される熱量を算出し、該熱量から熱媒の流量を算出するので、流量計を用いることなく、熱媒の流量を算出することができる。また、本発明の熱媒流量推定装置は、熱媒流量算出手段が、前記第1パラメータ算出手段によって算出された前記第1パラメータに基づいた前記圧縮機の吸込風量、及び前記圧縮機に吸い込まれる冷媒の密度から、前記蒸発器を流れる冷媒の流量を算出し、該算出した冷媒の流量、及び前記蒸発器の入口側のエンタルピと出口側のエンタルピの差から、前記蒸発器において冷媒と熱媒との間で交換される熱量を算出し、該算出した熱量、及び熱媒が前記蒸発器へ流入出する温度の差に基づいて、熱媒の流量を算出してもよい。本発明によれば、冷媒や熱媒の圧力や温度を計測する計測器による計測結果等を用いて、容易に熱媒の流量を算出できる。また、本発明の熱媒流量推定装置は、前記圧縮機には、回転数の制御が可能とされており、前記記憶手段が、前記圧縮機の回転数に応じて異なる複数の前記空力特性マップを記憶し、第1パラメータ算出手段が、前記圧縮機の回転数に対応した前記空力特性マップから前記第2パラメータに応じた前記第1パラメータを算出してもよい。本発明によれば、圧縮機の回転数に対応した空力特性マップから第2パラメータに応じた第1パラメータを算出するので、熱媒の流量をより精度高く算出できる。また、本発明の熱媒流量推定装置は、前記圧縮機には、冷媒入口に冷媒流量を調節するベーンが備えられており、前記記憶手段が、前記ベーンの開度に応じて異なる複数の前記空力特性マップを記憶し、第1パラメータ算出手段が、前記ベーンの開度に対応した前記空力特性マップから前記第2パラメータに応じた前記第1パラメータを算出してもよい。本発明によれば、圧縮機の冷媒入口に設けられたベーンの開度に対応した空力特性マップから第2パラメータに応じた第1パラメータを算出するので、熱媒の流量をより精度高く算出できる。また、本発明の熱媒流量推定装置は、前記凝縮器と前記蒸発器の間には、前記凝縮器内にある冷媒を前記蒸発器へと流すためのバイパス配管が設けられると共に、該バイパス配管を流れる冷媒の流量を調整するための弁が設けられ、前記記憶手段が、前記弁の開度に応じて異なる複数の前記空力特性マップを記憶し、第1パラメータ算出手段が、前記弁の開度に対応した前記空力特性マップから前記第2パラメータに応じた前記第1パラメータを算出してもよい。本発明によれば、凝縮器と蒸発器とをバイパスするバイパス配管に設けられた弁の開度に対応した空力特性マップから第2パラメータに応じた第1パラメータを算出するので、熱媒の流量をより精度高く算出できる。また、本発明に係る熱源機は、冷媒を圧縮する圧縮機と、圧縮された冷媒を熱源媒体によって凝縮させる凝縮器と、凝縮された冷媒を蒸発させると共に該冷媒と熱媒とを熱交換する蒸発器と、請求項1から請求項5の何れか1項に記載の熱媒流量推定装置と、を備える。また、本発明に係る熱媒流量推定方法は、冷媒を圧縮する圧縮機と、圧縮された冷媒を熱源媒体によって凝縮させる凝縮器と、凝縮された冷媒を蒸発させると共に該冷媒と熱媒とを熱交換する蒸発器と、を備えた熱源機の熱媒の流量を推定する熱媒流量推定方法であって、前記圧縮機の吸込風量を反映した第1パラメータと、前記圧縮機のヘッドを反映した第2パラメータとで表示されたマップ上に、旋回失速となる旋回失速線、及び前記圧縮機が吸い込む冷媒の音速を示す複数の機械マッハ数線が示された空力特性マップが記憶手段に予め記憶されており、前記第2パラメータを算出することによって、前記空力特性マップから該第2パラメータに応じた前記第1パラメータを算出する第1工程と、前記第1工程によって算出された前記第1パラメータに応じた前記圧縮機の吸込風量に基づいて、前記蒸発器において冷媒と熱媒との間で交換される熱量を算出し、該熱量に基づいて熱媒の流量を算出する第2工程と、を含む。本発明によれば、流量計を用いることなく、熱媒の流量を算出することができる、という優れた効果を有する。以下に、本発明に係る熱媒流量推定装置、熱源機、及び熱媒流量推定方法の一実施形態について、図面を参照して説明する。第1実施形態以下、本発明の第1実施形態について説明する。図1に、本第1実施形態に係る熱源機の一例であるターボ冷凍機10の構成を示す。ターボ冷凍機10は、冷媒を圧縮する圧縮機12と、圧縮機12によって圧縮された高温高圧のガス冷媒を熱源媒体冷却水によって凝縮する凝縮器14と、凝縮器14にて凝縮された液相の冷媒液冷媒に対して過冷却を与えるサブクーラ16と、サブクーラ16からの液冷媒を膨張させる高圧膨張弁18と、高圧膨張弁18に接続されるとともに圧縮機12の中間段及び低圧膨張弁20に接続される中間冷却器22と、低圧膨張弁20によって膨張させられた液冷媒を蒸発させると共に該冷媒と熱媒冷水とを熱交換する蒸発器24を備えている。圧縮機12は、遠心式の2段圧縮機であり、電源11からの入力周波数を変更するインバータ13により回転数制御された電動モータ28によって駆動されている。圧縮機12の冷媒吸入口には、吸入冷媒流量を制御する入口ベーンIGV32が設けられており、圧縮機12の容量制御が可能となっている。また、圧縮機12は、吸い込む冷媒の温度以下、圧縮機吸込温度Tsという。を計測するための吸込温度センサー17、吸い込む冷媒の圧力以下、圧縮機吸込圧力Psという。を計測するための吸込圧力センサー19が設けられている。吸込温度センサー17及び吸込圧力センサー19の出力は、制御装置30に入力される。サブクーラ16は、凝縮器14の冷媒流れ下流側に、凝縮された冷媒に対して過冷却を与えるように設けられている。凝縮器14及びサブクーラ16には、これらを冷却するための冷却伝熱管34が挿通されている。冷却伝熱管34の冷却水の出口温水出口には、温水出口温度センサー54が設けられている。温水出口温度センサー54の出力は、制御装置30に入力される。また、熱交換器である蒸発器24には、蒸発器24内の圧力である蒸発器圧力Peを計測する圧力センサー60が設けられている。この圧力センサー60の出力は、制御装置30に入力される。蒸発器24において吸熱されることによって定格温度例えば7°Cの冷媒が得られる。蒸発器24には、外部負荷へ供給される冷水を冷却するための冷水伝熱管36が挿通されている。なお、蒸発器24よりも上流側の冷水伝熱管36には、蒸発器24内へ流入する冷水の入口温度Toを計測する冷水入口温度センサー64が設けられている。蒸発器24よりも下流側の冷水出口ノズルには、蒸発器24から流出した冷水の出口温度Tiを計測する冷水出口温度センサー62が設けられている。冷水入口温度センサー64及び冷水出口温度センサー62の出力は、制御装置30に入力される。また、凝縮器14の気相部と蒸発器24の気相部との間には、ホットガスバイパス以下、HGBP]という。配管38が設けられている。そして、HGBP配管38には、HGBP配管38内を流れる冷媒の流量を制御するためのHGBP弁40が設けられている。HGBP弁40によってHGBP流量を調整することにより、入口ベーン32では制御が十分でない非常に小さな領域の容量制御が可能となっている。制御装置30は、ターボ冷凍機10の全体の制御を司り、回転数制御部30a、冷水流量推算部30b、及び膨張弁開度制御部30cを備える。回転数制御部30aは、ターボ冷凍機10各部における状態量圧力、温度等に基づいて、電動モータ28の指示回転数に応じた指示周波数をインバータ13に出力する。冷水流量推算部30bは、冷水流量を算出し、該算出結果を膨張弁開度制御部30cへ出力する。膨張弁開度制御部30cは、ターボ冷凍機10各部の状態量圧力、温度等及び冷水流量推算部30bから入力された冷水流量に基づいて、膨張弁開度指令値を生成し、該膨張弁開度指令値を高圧膨張弁18及び低圧膨張弁20へ送信することによって、高圧膨張弁18及び低圧膨張弁20の開度を制御する。さらに、制御装置30は、入口ベーン32の開度、HGBP弁40の開度等、ターボ冷凍機10の制御に必要な各種機器も制御する。ここで、ターボ冷凍機10の冷凍能力Qは、蒸発器24を流れる冷水の入口温度To及び出口温度Ti並びに冷水流量Gwに基づいて得られる。具体的には、下式1のように、冷水の出入口温度差Ti-Toに冷水流量Gwkgs及び冷水比熱cpkJを乗じて、冷凍能力Qが得られる。この冷凍能力Qと、圧縮機12の出入口における冷媒ガスのエンタルピ差Δhとに基づいて、次式2により、蒸発器24を流れる冷媒の流量である蒸発器冷媒流量Geが得られる。ここで、kは定数である。そして、蒸発器冷媒流量Geと、飽和ガスの比体積VTem3kgと、圧縮機12の羽根車の外径Dmと、蒸発器圧力Peから算出される飽和温度Teにおける吸込冷媒音速aTemsとに基づいて、次式3により、流量変数θが得られる。この流量変数は、圧縮機12の吸込風量を反映した無次元数である。このように、流量変数θは、冷凍能力Qおよび蒸発器圧力Peから得られるものである。圧力変数Ωは、圧縮機12のヘッドを反映した無次元数であり、凝縮器圧力Pc、蒸発器圧力Pe及び蒸発器圧力Peから算出される飽和温度Teから得られる冷媒ガスのエンタルピ差ΔhTeと、蒸発器24の蒸発器圧力Peから算出される飽和温度Teにおける吸込冷媒音速aTeとから、次式4により得られる。このように、圧力変数Ωは、凝縮器圧力Pc、蒸発器圧力Peから得られるものであり、羽根車の周速と無関係に得られる。以上の流量変数θ及び圧力変数Ωによって、圧縮機12の現在の運転状態が推定される。制御装置30が有する記憶部36は、圧縮機12の空力特性マップ42を備えている。この空力特性マップ42は、予め圧縮機12の運転試験を綿密に行うことで作成されるものであり、圧力変数Ωに対する流量変数θのマップ上に、圧縮機12が旋回失速を起こす旋回失速線Lが示されたものである。例えば、図2に示すような空力特性マップ42が得られる。この空力特性マップ42において、旋回失速線Lよりも下側の領域は、旋回失速やサージングを起こさない安定領域Sとされ、旋回失速線Lよりも上側の領域は、旋回失速やサージングを起こす不安定領域NSとされる。本実施形態において、この空力特性マップ42は、入口ベーン32の開度が最大開度である100のもの最大開度時マップである。また、空力特性マップ42には、機械マッハ数圧縮機12が吸い込む冷媒の音速である吸込冷媒音速を示す複数の機械マッハ数線Mが示されている。各機械マッハ数線は同一の値の機械マッハ数を示し、上方に行くに従い機械マッハ数は大きくなる。そして、流量変数θは、圧力変数Ω及び機械マッハ数で特定されるため、圧力変数Ω及び機械マッハ数を算出することで、流量変数θ、すなわち3式を変形することで圧縮機12の吸込風量の算出が可能である。ここで、本第1実施形態に係るターボ冷凍機10は、流量を計測する流量センサーは高価である、部品点数を少なくする等の理由により、冷水や冷却水の流量を計測する流量センサーを備えていない。しかしながら、冷凍機を設計値で運転させるためには、冷水の流量の管理が必要である。そこで、本第1実施形態に係るターボ冷凍機10は、圧力変数Ωを算出し、空力特性マップから圧力変数Ωに応じた流量変数θを算出し、算出した流量変数θに応じた圧縮機12の吸込風量に基づいて、蒸発器24において冷媒と冷水との間で交換される熱量を算出し、該熱量に基づいて冷水の流量を算出する冷水流量推算処理を行う。すなわち、冷水流量推算処理は、圧縮機12の運転状態に対応した流量変数θを算出し、流量変数θから算出される圧縮機12の吸込み風量に基づいた熱量を用いて、蒸発器24における冷媒と冷水との熱バランスにより冷水の流量を算出する。図3は、冷水流量推算処理を行う場合に、制御装置30が備える冷水流量推算部30bによって実行される冷水流量推算プログラムの処理の流れを示すフローチャートであり、冷水流量推算プログラムは冷水流量推算部30bが備える記憶部の所定領域に予め記憶されている。なお、本プログラムは、例えば、所定時間間隔毎に実行される。まず、ステップ100では、吸込冷媒音速aTe、圧力変数Ω、及び吸込冷媒密度ρを算出する。吸込冷媒音速aTeは、上述したように蒸発器圧力Peから算出される飽和温度Teに基づいて算出され、圧力変数Ωは、4式に基づいて算出される。吸込冷媒密度ρは、圧縮機12に設けられた吸込温度センサー17で計測された圧縮機吸込温度Tsと吸込圧力センサー19で計測された圧縮機吸込圧力Psとから算出される。次のステップ102では、算出した圧力変数Ω及び吸込冷媒音速aTeに対応した流量変数θを空力特性マップ42から算出する。すなわち、ステップ100,ステップ102によって、圧縮機12の運転状態に対応した流量変数θが算出される。次のステップ104では、蒸発器冷媒流量Geを下記5式によって算出する。ここで、Qsは、圧縮機12の吸込風量m3sである。そして、吸込風量Qsは、ステップ102で算出した流量変数θを用いて下記6式から算出される。なお、下記6式は、吸込風量Qsを算出するために3式を変形した式であり、吸込冷媒音速aTeはステップ100で算出され、圧縮機12の羽根車の外径Dは圧縮機12の設計値から求められる。次のステップ106では、蒸発器24の入口側のエンタルピhei及び蒸発器24の出口側のエンタルピheoを算出する。次のステップ108では、蒸発器24において冷水と冷媒との間で交換される熱量である蒸発器交換熱量QekWを下記7式から算出する。次のステップ110では、冷水流量Gwを算出し、本プログラムを終了する。このように、ステップ104からステップ110によって、蒸発器24における冷媒と冷水との熱バランスにより冷水流量が算出される。そして、冷水流量推算部30bは、算出した冷水流量Gwを膨張弁開度制御部30cへ出力し、膨張弁開度制御部30cは、ターボ冷凍機10各部の状態量圧力、温度等及び冷水流量推算部30bから入力された冷水流量に基づいて、膨張弁開度指令値を生成する。以上説明したように、本第1実施形態に係る制御装置30は、圧縮機12の吸込風量を反映した流量変数θと、圧縮機12のヘッドを反映した圧力変数Ωとで表示されたマップ上に、旋回失速となる旋回失速線、及び圧縮機12が吸い込む冷媒の音速を示す複数の機械マッハ数線が示された空力特性マップ42を記憶した記憶部36を備え、冷水流量推算部30bによって、圧力変数Ωを算出し、空力特性マップ42から圧力変数Ωに応じた流量変数θを算出し、算出した流量変数θに応じた圧縮機12の吸込風量に基づいて、蒸発器24において冷媒と冷水との間で交換される熱量を算出し、該熱量に基づいて冷水の流量を算出する。従って、本第1実施形態に係る制御装置30は、流量計を用いることなく、冷水の流量を算出することができる。また、冷水流量推算部30bは、算出した流量変数θに基づいた圧縮機12の吸込風量、及び圧縮機12に吸い込まれる冷媒の密度から、蒸発器24を流れる冷媒の流量を算出し、算出した冷媒の流量、及び蒸発器24の入口側のエンタルピと出口側のエンタルピの差から、蒸発器24において冷媒と冷水との間で交換される熱量を算出し、算出した熱量、及び冷水が蒸発器24へ流入出する温度の差に基づいて、冷水の流量を算出する。従って、本第1実施形態に係る制御装置30は、冷媒や冷水の圧力や温度を計測する計測器による計測結果等を用いて、容易に冷水の流量を算出できる。第2実施形態以下、本発明の第2実施形態について説明する。本第2実施形態に係るターボ冷凍機10の構成は、図1に示す第1実施形態に係るターボ冷凍機10の構成と同様であるので説明を省略する。ただし、本第2実施形態に係る記憶部36は、インバータ13から電動モータ28に送られる指示周波数を制御することによって、圧縮機12の回転数が制御可能とされていることから、圧縮機12の回転数に応じて異なる複数の空力特性マップ42を記憶している。なお、本第2実施形態に係る空力特性マップ42は、圧縮機12の回転数が高いほど、同じ圧力変数に対する流量変数が大きくなるように表わされている。そして、本第2実施形態では、冷水流量推算プログラムのステップ102において、圧縮機12の回転数指示周波数に対応した空力特性マップ42を記憶部36から選択し、選択した空力特性マップ42から圧力変数Ωに応じた流量変数θを算出する。以上のように、本第2実施形態に係る制御装置30は、圧縮機12の回転数に対応した空力特性マップ42から圧力変数Ωに応じた流量変数θを算出するので、冷水の流量をより精度高く算出できる。第3実施形態以下、本発明の第3実施形態について説明する。本第3実施形態に係るターボ冷凍機10の構成は、図1に示す第1実施形態に係るターボ冷凍機10の構成と同様であるので説明を省略する。ただし、本第3実施形態に係る記憶部36は、ターボ冷凍機10が入口ベーン32を備えていることから、入口ベーン32の開度に応じて異なる複数の空力特性マップ42を記憶している。なお、本第3実施形態に係る空力特性マップ42は、入口ベーン32の開度が大きいほど、同じ圧力変数に対する流量変数が大きくなるように表わされている。そして、本第3実施形態では、冷水流量推算プログラムのステップ102において、入口ベーン32の開度に対応した空力特性マップ42を記憶部36から選択し、選択した空力特性マップ42から圧力変数Ωに応じた流量変数θを算出する。以上のように、本第3実施形態に係る制御装置30は、入口ベーン32の開度に対応した空力特性マップ42から圧力変数Ωに応じた流量変数θを算出するので、冷水の流量をより精度高く算出できる。第4実施形態以下、本発明の第4実施形態について説明する。本第4実施形態に係るターボ冷凍機10の構成は、図1に示す第1実施形態に係るターボ冷凍機10の構成と同様であるので説明を省略する。ただし、本第4実施形態に係る記憶部36は、ターボ冷凍機10がHGBP配管38と共にHGBP弁40を備えていることから、HGBP弁40の開度に応じて異なる複数の空力特性マップ42を記憶している。なお、本第4実施形態に係る空力特性マップ42は、HGBP弁40の開度が大きいほど、同じ圧力変数に対する流量変数が大きくなるように表わされている。そして、本第4実施形態では、冷水流量推算プログラムのステップ102において、HGBP弁40の開度に対応した空力特性マップ42を記憶部36から選択し、選択した空力特性マップ42から圧力変数Ωに応じた流量変数θを算出する。以上のように、本第4実施形態に係る制御装置30は、HGBP弁40の開度に対応した空力特性マップ42から圧力変数Ωに応じた流量変数θを算出するので、冷水の流量をより精度高く算出できる。以上、本発明を、上記各実施形態を用いて説明したが、本発明の技術的範囲は上記実施形態に記載の範囲には限定されない。発明の要旨を逸脱しない範囲で上記各実施形態に多様な変更または改良を加えることができ、該変更または改良を加えた形態も本発明の技術的範囲に含まれる。また、上記各実施形態では、凝縮器14に挿通される冷却伝熱管34内を流れる熱源媒体を冷却水とする形態について説明したが、本発明は、これに限定されるものではなく、熱源媒体を気体外気とし、凝縮器を空気熱交換器とする形態としてもよい。また、上記各実施形態では、冷凍運転を行うターボ冷凍機10に本発明を適用した場合について説明したが、これに限らず、本発明をヒートポンプ運転も可能なヒートポンプ式ターボ冷凍機に適用してもよい。また、上記各実施形態では、ターボ冷凍機10を、遠心圧縮機を用いる形態について説明したが、本発明はこれに限定されるものではなく、他の圧縮形式であっても適用することができ、例えばスクリュー圧縮機を用いたスクリューヒートポンプであってもよい。また、上記各実施形態で説明した冷水流量推算プログラムの処理の流れも一例であり、本発明の主旨を逸脱しない範囲内において不要なステップを削除したり、新たなステップを追加したり、処理順序を入れ替えたりしてもよい。"