(57)【要約】【課題】  室内の精密なモデルを作成する室内温度分布モデル作成装置及び室内温度分布モデル作成方法を提供する。【解決手段】  モデル化する対象のデータを入力するステップと、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンＡを作成するステップと、室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンＢを作成するステップと、発熱体から発生する熱による上昇流を生じる熱プルームゾーンＣを作成するステップと、ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網２を作成するステップと、ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンス３３，３４を求め、熱回路網３を作成するステップと、モデル化した室内の温度分布を出力するステップと、を備えたことを特徴とする。【選択図】図１
2009053749.tif 000002
【特許請求の範囲】
【請求項１】
室内の温度分布をモデル化する室内温度分布モデル作成装置において、
室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンを作成する温度成層化領域作成手段と、
室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンを作成する温度境界層領域作成手段と、
発熱体から発生する熱による上昇流を生じる熱プルームゾーンを作成する熱プルーム領域作成手段と、
前記ゾーン底面の静圧と前記ゾーン間の空気移動量を求める換気回路網を作成する換気回路網作成手段と、
前記ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンスを求め、熱回路網を作成する熱回路網作成手段と、
を備えたことを特徴とする室内温度分布モデル作成装置。
【請求項２】
前記温度境界層領域作成手段は複数の温度境界層ゾーンを作成し、前記熱プルーム領域作成手段は複数の熱プルームゾーンを作成することを特徴とする請求項１に記載の室内温度分布モデル作成装置。
【請求項３】
前記熱回路網作成手段は、前記一般化熱コンダクタンスの数値を小さくまたは大きくすることで、前記温度成層化領域内の自然対流と強制対流の両状態をモデル化することを特徴とする請求項１又は請求項２に記載の室内温度分布モデル作成装置。
【請求項４】
前記熱回路網作成手段は、熱橋による熱橋熱コンダクタンスを作成し、前記一般化熱コンダクタンスと並列の熱流となる経路としてモデル化することを特徴とする請求項１乃至請求項３のいずれかに記載の室内温度分布モデル作成装置。
【請求項５】
前記熱橋による熱橋熱コンダクタンスと前記一般化熱コンダクタンスとの比を示す熱橋割合と、温度成層化ゾーン間の混合流熱コンダクタンスと、換気回路網の流路の圧力損失係数とを、予測値と実測値との誤差の評価関数を最小にするように求めるパラメータ演算手段を有することを特徴とする請求項１乃至請求項４のいずれかに記載の室内温度分布モデル作成装置。
【請求項６】
室内の温度分布をモデル化する室内温度分布モデル作成方法において、
モデル化する対象のデータを入力するステップと、
室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンを作成するステップと、
室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンを作成するステップと、
発熱体から発生する熱による上昇流を生じる熱プルームゾーンを作成するステップと、
前記ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網を作成するステップと、
前記ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンスを求め、熱回路網を作成するステップと、
モデル化した室内の温度分布を出力するステップと、
を備えたことを特徴とする室内温度分布モデル作成方法。
【請求項７】
前記温度境界層ゾーン及び前記熱プルームゾーンが、複数のゾーンにより作成されることを特徴とする請求項６に記載の室内温度分布モデル作成方法。
【請求項８】
前記一般化熱コンダクタンスの数値を小さくまたは大きくすることで、前記温度成層化領域内の自然対流と強制対流の両状態をモデル化することを特徴とする請求項６又は請求項７に記載の室内温度分布モデル作成方法。
【請求項９】
熱橋による熱橋熱コンダクタンスを作成し、前記一般化熱コンダクタンスと並列の熱流となる経路としてモデル化することを特徴とする請求項６乃至請求項８のいずれかに記載の室内温度分布モデル作成方法。
【請求項１０】
前記熱橋による熱橋熱コンダクタンスと前記一般化熱コンダクタンスとの比を示す熱橋割合と、温度成層化ゾーン間の混合流熱コンダクタンスと、換気回路網の流路の圧力損失係数とを、予測値と実測値との誤差の評価関数を最小にするように求めることを特徴とする請求項６乃至請求項９のいずれかに記載の室内温度分布モデル作成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、室内の温度分布をモデル化し、熱負荷等を検討できる室内温度分布モデル作成装置及び室内温度分布モデル作成方法に関する。
【背景技術】
【０００２】
従来、室内の温度分布を予測する手段として、数値流体解析が用いられていた。しかしながら、数値流体解析は、計算負荷が大きく、数値的な安定性の問題等から長時間の現象を予測することが難しかった。
【０００３】
そこで、本発明者は、建築物の熱回路網及び換気回路網をモデル化し、熱負荷計算を実行することで、適度な室内環境を実現するものを発表した（非特許文献１参照）。
【非特許文献１】宇田川,石田,石野等,「熱負荷算法小委員会報告書-大空間の熱負荷計算法-」,1993年3月,空衛学会,空気調和設備委員会,熱負荷算法小委員会
【発明の開示】
【発明が解決しようとする課題】
【０００４】
ところで、上記非特許文献１に記載されたものは、計算するモデルを、１室に対して１節点で設定しているため、室内の温度分布に対しては、考慮されていなかった。
【０００５】
本発明は上記課題を解決し、室内の精密なモデルを作成する室内温度分布モデル作成装置及び室内温度分布モデル作成方法を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明は上記課題を解決するものであって、室内の温度分布をモデル化する室内温度分布モデル作成装置において、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンを作成する温度成層化領域作成手段と、室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンを作成する温度境界層領域作成手段と、発熱体から発生する熱による上昇流を生じる熱プルームゾーンを作成する熱プルーム領域作成手段と、前記ゾーン底面の静圧と前記ゾーン間の空気移動量を求める換気回路網を作成する換気回路網作成手段と、前記ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンスを求め、熱回路網を作成する熱回路網作成手段と、を備えたことを特徴とする。
【０００７】
また、前記温度境界層領域作成手段は複数の温度境界層ゾーンを作成し、前記熱プルーム領域作成手段は複数の熱プルームゾーンを作成することを特徴とする。
【０００８】
また、前記熱回路網作成手段は、前記一般化熱コンダクタンスの数値を小さくまたは大きくすることで、前記温度成層化領域内の自然対流と強制対流の両状態をモデル化することを特徴とする。
【０００９】
また、前記熱回路網作成手段は、熱橋による熱橋熱コンダクタンスを作成し、前記一般化熱コンダクタンスと並列の熱流となる経路としてモデル化することを特徴とする。
【００１０】
また、前記熱橋による熱橋熱コンダクタンスと前記一般化熱コンダクタンスとの比を示す熱橋割合と、温度成層化ゾーン間の混合流熱コンダクタンスと、換気回路網の流路の圧力損失係数とを、予測値と実測値との誤差の評価関数を最小にするように求めるパラメータ演算手段を有することを特徴とする。
【００１１】
さらに、本発明は上記課題を解決する方法であって、室内の温度分布をモデル化する室内温度分布モデル作成方法において、モデル化する対象のデータを入力するステップと、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンを作成するステップと、室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンを作成するステップと、発熱体から発生する熱による上昇流を生じる熱プルームゾーンを作成するステップと、前記ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網を作成するステップと、前記ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンスを求め、熱回路網を作成するステップと、モデル化した室内の温度分布を出力するステップと、を備えたことを特徴とする。
【００１２】
また、前記温度境界層領域及び前記熱プルーム領域が、複数のゾーンにより作成されることを特徴とする。
【００１３】
また、前記一般化熱コンダクタンスの数値を小さくまたは大きくすることで、前記温度成層化領域内の自然対流と強制対流の両状態をモデル化することを特徴とする。
【００１４】
また、熱橋による熱橋熱コンダクタンスを作成し、前記一般化熱コンダクタンスと並列の熱流となる経路としてモデル化することを特徴とする。
【００１５】
また、前記熱橋による熱橋熱コンダクタンスと前記一般化熱コンダクタンスとの比を示す熱橋割合と、温度成層化ゾーン間の混合流熱コンダクタンスと、換気回路網の流路の圧力損失係数とを、予測値と実測値との誤差の評価関数を最小にするように求めることを特徴とする。
【発明の効果】
【００１６】
本発明によれば、室内の温度分布をモデル化する室内温度分布モデル作成装置において、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンを作成する温度成層化領域作成手段と、室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンを作成する温度境界層領域作成手段と、発熱体から発生する熱による上昇流を生じる熱プルームゾーンを作成する熱プルーム領域作成手段と、前記ゾーン底面の静圧と前記ゾーン間の空気移動量を求める換気回路網を作成する換気回路網作成手段と、前記ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンスを求め、熱回路網を作成する熱回路網作成手段と、を備えたので、精密なモデルを作成することができる。そして、精密なモデルにより、冷暖房の定常・非定常の熱負荷の予測ができ、空調機器設定の目安が正確に計算できる。また、部材毎の伝熱量が概算できるので、熱負荷低減アイテムの効果が概算できる。
【００１７】
また、前記温度境界層領域作成手段は複数の温度境界層ゾーンを作成し、前記熱プルーム領域作成手段は複数の熱プルームゾーンを作成するので、加熱高さによる上昇力の違いが表現できて，より精密なモデルができる。
【００１８】
また、前記熱回路網作成手段は、前記一般化熱コンダクタンスの数値を小さくまたは大きくすることで、前記温度成層化領域内の自然対流と強制対流の両状態をモデル化するので、強制対流と自然対流の両状態を同じモデル構造で表現できる。
【００１９】
また、前記熱回路網作成手段は、熱橋による熱橋熱コンダクタンスを作成し、前記一般化熱コンダクタンスと並列の熱流となる経路としてモデル化するので、より精密なモデルができる。
【００２０】
また、前記熱橋による熱橋熱コンダクタンスと前記一般化熱コンダクタンスとの比を示す熱橋割合と、温度成層化ゾーン間の混合流熱コンダクタンスと、換気回路網の流路の圧力損失係数とを、予測値と実測値との誤差の評価関数を最小にするように求めるパラメータ演算手段を有するので、最適なパラメータを求めることができ、より精密なモデルができる。
【００２１】
さらに、本発明は上記課題を解決する方法であって、室内の温度分布をモデル化する室内温度分布モデル作成方法において、モデル化する対象のデータを入力するステップと、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンを作成するステップと、室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンを作成するステップと、発熱体から発生する熱による上昇流を生じる熱プルームゾーンを作成するステップと、前記ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網を作成するステップと、前記ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンスを求め、熱回路網を作成するステップと、モデル化した室内の温度分布を出力するステップと、を備えたので、精密なモデルを作成することができる。そして、精密なモデルにより、冷暖房の定常・非定常の熱負荷の予測ができ、空調機器設定の目安が正確に計算できる。また、部材毎の伝熱量が概算できるので、熱負荷低減アイテムの効果が概算できる。
【００２２】
また、前記温度境界層領域及び前記熱プルーム領域が、複数のゾーンにより作成されるので、より精密なモデルができる。
【００２３】
また、前記一般化熱コンダクタンスの数値を小さくまたは大きくすることで、前記温度成層化領域内の自然対流と強制対流の両状態をモデル化するので、強制対流と自然対流の両状態を同じモデル構造で表現できる。
【００２４】
また、熱橋による熱橋熱コンダクタンスを作成し、前記一般化熱コンダクタンスと並列の熱流となる経路としてモデル化するので、より精密なモデルができる。
【００２５】
また、前記熱橋による熱橋熱コンダクタンスと前記一般化熱コンダクタンスとの比を示す熱橋割合と、温度成層化ゾーン間の混合流熱コンダクタンスと、換気回路網の流路の圧力損失係数とを、予測値と実測値との誤差の評価関数を最小にするように求めるので、最適なパラメータを求めることができ、より精密なモデルができる。
【発明を実施するための最良の形態】
【００２６】
以下、図面を参照して本発明にかかる実施形態の室内温度分布作成装置１について説明する。
【００２７】
図１は、本実施形態の室内温度分布作成装置１が作成するモデルの概要を示す図である。図１において、Ａは温度成層化領域、Ｂは温度境界層領域、Ｃは熱プルーム領域を示す。また、換気回路網２では、矢印２１は空気の移動を表す換気回路網の流路を示す。さらに、熱回路網３では、節点３１は例えば外気温度等の既知の温度を表す与条件温度節点、節点３２は温度が未知な未知数温度節点を示し、第１熱コンダクタンス３３は伝導・伝達・貫流・放射等の一般化熱コンダクタンス、第２熱コンダクタンス３４は移流・対流による一般化熱コンダクタンスを示し、太矢印３５は節点への発熱量を示す。なお、一般化コンダクタンスとは、伝導、伝達、放射、貫流、移流等すべての伝熱形態の熱の伝わりやすさを表す。
【００２８】
本実施形態では、温度成層化領域Ａ、温度境界層領域Ｂ及び熱プルーム領域Ｃに区分すると共に、それぞれ上下に層状のゾーンＡ１，Ａ２・・・Ｃ３，Ｃ４に分割して、換気回路網２と熱回路網３を構築する。
【００２９】
温度成層化領域Ａは、室内の上下温度分布を複数の層に分割し、離散的に表現したものであり、各層内での温度は一様である。温度成層化領域Ａ間の平均的な流れは上向きまたは下向きの押し出し流れであり、こうしたピストン流れとすることが温度成層化を模擬する要件となる。また、自然対流（熱対流のみ）と強制対流（ファン作動時など）の両状態は、温度成層化領域Ａ内の各ゾーン間の混合流の一般化熱コンダクタンス３３，３４の数値を，小さくまたは大きくすることでモデル化する。
【００３０】
温度境界層領域Ｂは、室内外を画成する外皮が外気に接することで外皮の内表面側に下降流や上昇流を生じるゾーンである。例えば、車両の外気温度が外皮を通じて車内を冷却すると、温度境界層領域Ｂを冷却し、温度境界層領域Ｂ内で下降流が生じ，この下降流が温度成層化領域Ａに入り込み、循環流を引き起こす。
【００３１】
熱プルーム領域Ｃは、発熱体から発生する熱による上昇流を生じるゾーンである。例えば、日射を受けた車室内の物体から発生する上昇流を表現するものであり、熱プルーム領域Ｃ全体で煙突状のものを想定し、やはり上下方向を分割することで，発熱体との接触位置が高い場合、低い場合の熱上昇力の違いを考慮できる様にした。
【００３２】
なお、ファンなどで強制対流加熱した場合には，室内は空気的に混合状態になることから，温度成層化領域Ａの全てのゾーンに発熱量は均等に分散して与えられる。これでも上下温度分布が生成されるのは押し出し上昇流が生じるからである。また、この強制対流加熱状態では熱プルーム領域Ｃ内には加熱源がなく、また熱プルーム流路の断面積が小さいことから熱プルームの流の影響はなくなる。こうして強制対流と自然対流の両状態は同じモデル構造で表現できる。
【００３３】
ここで、換気回路網２と熱回路網３について説明する（詳細については、非特許文献１を参照）。
【００３４】
換気回路網２は、各ゾーン底面の静圧と複数ゾーン間の流量を示すものである。
【００３５】
前提として、各ゾーンでは風量の収支が成り立つとする。また、各ゾーン底面に静圧ｐ(Pa)を持ち、ゾーンの空気密度をρ(kg/m3)とすると、この底面から上方へｈ(m)上がった点での静圧はｐ−ρ・ｈ・ｇ(Pa)とする。空気密度はそのゾーンの空気温度の関数である。流路は、圧力損失係数、指数及び流路面積、さらに基準面からの高さの属性を持ち、流路の通過風量は，前後の静圧差と圧力損失係数、指数および流路面積の関数である。そして、ゾーン数分の静圧に関する非線形連立方程式を修正ニュートンラプソン法で解くことで風量が求められる。
【００３６】
熱回路網３は、熱の移動を示すものである。前提として、各節点では、流入熱流から流出熱流を引いたものが、熱容量に温度変化を乗じたものに等しくなるという熱流収支が成り立つものとする。これを完全連結システムの節点方程式とする。また、伝導、表面伝達、長波長放射、移流など全ての伝熱形態の熱の伝わりやすさを一般化熱コンダクタンス３３，３４という一種類の係数で表現する。従って空気流動による一般化熱コンダクタンス３３，３４も数式記号上は同じ係数となる。
【００３７】
そして、全節点の温度に関する連立常微分方程式が前述の節点方程式から構成され、連立常微分方程式を後退差分などの時間積分法によって解き、毎時間ステップの温度の解を算出する。
【００３８】
次に、換気回路網２と熱回路網３の連成について説明する。換気回路網２の駆動条件となる各ゾーンの密度は熱回路網３の解の温度によって計算される。熱回路網３の係数となる風量は換気回路網２の解から得られる。非定常の時間ステップ毎に、一方の回路網への条件はもう一方の回路網の前時間ステップでの結果をもとにつくる。こうして両モデルの連成がとられる．
図２は本実施形態の室内温度分布作成装置１を用いたシステム構成図である。
【００３９】
入力手段１０１は、モデル化する対象のデータを入力するものである。記憶手段１０２は、演算式等の不変のデータを記憶するものである。入力手段１０１から入力されたデータ及び記憶手段１０２に記憶されたデータは、制御手段１１０に入力される。
【００４０】
室内温度分布作成装置１の制御手段１１０は、温度成層化領域作成手段１１１、温度境界層領域作成手段１１２、熱プルーム領域作成手段１１３、換気回路網作成手段１１４及び熱回路網作成手段１１５、パラメータ演算手段１１６を有する。
【００４１】
温度成層化領域作成手段１１１は、室内を複数のゾーンに分割し、上下温度分布を離散的に表現する温度成層化ゾーンを作成する。温度境界層領域作成手段１１２は、室内外を画成する外皮が外気に接することで外皮の内表面側に下降流や上昇流を生じるゾーンを作成する。なお、温度境界層ゾーンは、複数のゾーンに分割し作成されてもよい。熱プルーム領域作成手段１１３は、発熱体から発生する熱による上昇流を生じるゾーンを作成する。なお、熱プルームゾーンは、複数のゾーンに分割し作成されてもよい。換気回路網作成手段１１４は、ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網を作成する。熱回路網作成手段１１５は、ゾーン間の熱の移動を求める熱回路網を作成する。パラメータ演算手段１１６については後述する。
【００４２】
制御手段１１０により作成されたモデルは、表示装置等の出力手段１２１に出力される。
【００４３】
図３は、車両５０を図１に示す領域に分割する例を示す図である。車両５０の車室５１内を上下４段の温度成層化領域Ａに区分し、温度境界層領域Ｂとして、屋根温度境界層Ｂ１、窓温度境界層Ｂ２、ドア温度境界層Ｂ３及び床温度境界層Ｂ４を設定したものである。
【００４４】
図４は、図３に示す車両５０に換気回路網２と熱回路網３を構築しモデル化した図である。この例では、熱プルーム領域ＣのＣ３ゾーンに発熱体としてヒーターを設置している。
【００４５】
次に、構造物内部の熱回路網３について説明する。図５は、熱橋５２を示す図、図６は車両５０の熱橋５２をモデル化した図である。熱橋５２とは、車室５１から外気への熱損失性能が実験値と合うようにするために考慮する板状部材内に存在するリブ等で熱が伝わり易い部分を示す。
【００４６】
熱橋５２による熱橋熱コンダクタンス３６は、図６の屋根部分に示すように、一般化熱コンダクタンス３３と並列の熱流となる経路として構成される。熱橋熱コンダクタンス３６をＣｂとし、一般化熱コンダクタンス３３をＣｓとすると、熱コンダクタンスの合成は、Ｃｓ＋Ｃｂとなる。
【００４７】
次に、作成するモデルのパラメータについて説明する。作成するモデルは、工学モデルなので、実験的に定めなければならない特性値がある。これらは、熱橋割合Ｃｂ／Ｃｓ、温度成層化ゾーン間の混合流熱コンダクタンスｃｉｊ及び換気回路網の圧力損失係数ζからなる三種のパラメータである。
【００４８】
これらパラメータをパラメータ演算手段１１６により最適化する方法について説明する。まず、パラメータを最適化するために評価関数を設定する。実験と予測計算で対応を検討したのは、図７に示す頭部温度θ１と足下温度θ２である。実測と予測は頭部と足元の平均温度が近いほど良いという評価を式（１）で行う。式（１）は、上下の平均温度は、室空気が持つ熱エネルギーに比例するので、一致すべきという観点から導かれる。
【００４９】
【数１】
2009053749.tif 000003
・・・（１）
ただし、pθjは頭部又は足元の予測計算温度、
mθjは実測温度、
Ｔは検討期間
である。
【００５０】
ここにjは成層化ゾーンの上下方向の位置番号を示し，一般に２以上あるものとする．
次に、頭部と足元それぞれで温度が近ければ良いという評価を式（２）で行う。
【００５１】
【数２】
2009053749.tif 000004
・・・（２）
これら式（１）と式（２）の両方を平均し、総合的に考慮する式（３）の誤差の評価関数Ｊを定める。また、誤差は発熱量Ｈの大きさに比例すると考えられるので、Ｈで除して正規化する。
【００５２】
【数３】
2009053749.tif 000005
・・・（３）
次に、式（３）を最小化する三種のパラメータを探索法により求める。まず、熱橋割合Ｃｂ／Ｃｓを０とする。次に、温度成層化ゾーン間の混合流熱コンダクタンスｃｉｊと換気回路網の圧力損失係数ζとをそれぞれ何通りか変化させ、評価関数を計算して、図８に示すような２次元的数表で最小化する組み合わせを求める。さらに、図９に示すように、熱橋割合Ｃｂ／Ｃｓを３次元方向として変化させ、評価関数を計算し、最適な三種のパラメータを求める。
【００５３】
次に、このような室内温度分布作成装置１により作成したモデルによる計算値と実測値との比較結果について説明する。図１０は実測時の外気温と日射量を示すグラフ、図１１は室温推移の比較を示すグラフである。
【００５４】
図１０において、点線は外気温、正方形を含む実線は水平面全日射量、菱形を含む点線は直達日射量、×を含む実線は拡散日射量である。実測時は晴天で、徐々に温度が上昇する高温の状態であった。
【００５５】
図１１において、○を含む実線は頭部付近のモデルによる計算値、菱形を含む実線は頭部付近の実測値、正方形を含む実線は足元付近のモデルによる計算値、×を含む実線は足元付近の実測値である。図１０に示すように、頭部と足元で約２０度の大きな温度差を有する状態において、モデルによる計算値は、実測値にかなり近い値を算出している。
【００５６】
このように、室内の温度分布をモデル化する室内温度分布モデル作成装置において、室内を複数のゾーンＡ１，Ａ２・・・に分割し、上下温度分布を表現する温度成層化領域Ａを作成する温度成層化領域作成手段１１１と、室内外を画成する外皮が外気に接することで外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンＢを作成する温度境界層領域作成手段１１２と、発熱体から発生する熱による上昇流を生じる熱プルームゾーンＣを作成する熱プルーム領域作成手段１１３とゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網２を作成する換気回路網作成手段１１４と、ゾーン間の熱の移動を表す一般化熱コンダクタンス３３，３４を求め、熱回路網３を作成する熱回路網作成手段１１５と、を有するので、精密なモデルを作成することができる。そして、精密なモデルにより、冷暖房の定常・非定常の熱負荷の予測ができ、空調機器設定の目安が正確に計算できる。また、部材毎の伝熱量が概算できるので、熱負荷低減アイテムの効果が概算できる。
【００５７】
また、温度境界層領域作成手段１１２は複数の温度境界層ゾーンＢ１，Ｂ２・・・を作成し、前記熱プルーム領域作成手段１１３は複数の熱プルームゾーンＣ１，Ｃ２・・・を作成するので、加熱高さによる上昇力の違いが表現できて，より精密なモデルができる。
【００５８】
また、熱回路網作成手段１１５は、一般化熱コンダクタンス３３，３４の数値を小さくまたは大きくすることで、温度成層化領域Ａ内の自然対流と強制対流の両状態をモデル化するので、強制対流と自然対流の両状態を同じモデル構造で表現できる。
【００５９】
また、熱回路網作成手段１１５は、熱橋５２による熱橋熱コンダクタンス３６を作成し、一般化熱コンダクタンス３３，３４と並列の熱流となる経路としてモデル化するので、より精密なモデルができる。
【００６０】
また、熱橋５２による熱橋熱コンダクタンス３６と一般化熱コンダクタンス３３，３４との比を示す熱橋割合Ｃｂ／Ｃｓと、温度成層化ゾーンＡ１，Ａ２・・・間の混合流熱コンダクタンスｃｉｊと、換気回路網の流路の圧力損失係数ζとを、予測値と実測値との誤差の評価関数Ｊを最小にするように求めるパラメータ演算手段１１６を有するので、最適なパラメータを求めることができ、より精密なモデルができる。
【００６１】
さらに、本発明は上記課題を解決する方法であって、室内の温度分布をモデル化する室内温度分布モデル作成方法において、モデル化する対象のデータを入力するステップと、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化領域Ａを作成するステップと、室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンＢを作成するステップと、発熱体から発生する熱による上昇流を生じる熱プルーム領域Ｃを作成するステップと、ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網２を作成するステップと、ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンス３３，３４を求め、熱回路網３を作成するステップと、モデル化した室内の温度分布を出力するステップと、を備えたので、精密なモデルを作成することができる。そして、精密なモデルにより、冷暖房の定常・非定常の熱負荷の予測ができ、空調機器設定の目安が正確に計算できる。また、部材毎の伝熱量が概算できるので、熱負荷低減アイテムの効果が概算できる。
【００６２】
また、温度境界層領域Ｂ及び前記熱プルーム領域Ｃが、複数のゾーンＢ１，Ｂ２，・・・Ｃ１，Ｃ２・・・により作成されるので、より精密なモデルができる。
【００６３】
また、一般化熱コンダクタンス３３，３４の数値を小さくまたは大きくすることで、温度成層化領域Ａ内の自然対流と強制対流の両状態をモデル化するので、強制対流と自然対流の両状態を同じモデル構造で表現できる。
【００６４】
また、熱橋５２による熱橋熱コンダクタンス３６を作成し、一般化熱コンダクタンス３３，３４と並列の熱流となる経路としてモデル化するので、より精密なモデルができる。
【００６５】
また、熱橋５２による熱橋熱コンダクタンス３６と一般化熱コンダクタンス３３，３４との比を示す熱橋割合Ｃｂ／Ｃｓと、温度成層化ゾーンＡ１，Ａ２・・・間の混合流熱コンダクタンスｃｉｊと、換気回路網の流路の圧力損失係数ζとを、予測値と実測値との誤差の評価関数Ｊを最小にするように求めるので、最適なパラメータを求めることができ、より精密なモデルができる。
【００６６】
なお、本実施形態では、車両５０のモデル作成について説明したが、これに限らず、建築物や構造物等の室内に適用してもよい。
【図面の簡単な説明】
【００６７】
【図１】本実施形態の室内温度分布作成装置が作成するモデルの概要を示す図である。
【図２】本実施形態の装置構成を示す図である。
【図３】車両を図１に示す領域に分割する例を示す図である。
【図４】図３に示す車両に換気回路網と熱回路網を構築しモデル化した図である。
【図５】熱橋を示す図である。
【図６】車両の熱橋をモデル化した図である。
【図７】車両内の頭部温度と足下温度を示す図である。
【図８】２次元的数表で最小化する組み合わせを示す図である。
【図９】熱橋割合を３次元方向として変化させることを示す図である。
【図１０】実測時の外気温と日射量を示すグラフである。
【図１１】室温推移の比較を示すグラフである。
【符号の説明】
【００６８】
１…室内温度分布作成装置、２…換気回路網、２１…流路、３…熱回路網、３１…与条件温度節点、３２…未知数温度節点、３３…第１熱コンダクタンス、３４…第２熱コンダクタンス、３５…節点への発熱量、３６…熱橋熱コンダクタンス、５０…車両、５１…車室、５２…熱橋、１０１…入力手段、１１０…制御手段、１１１…温度成層化領域作成手段、１１２…温度境界層領域作成手段、１１３…熱プルーム領域作成手段、１１４…換気回路網作成手段、１１５…熱回路網作成手段、１２１出力手段
【図１】
2009053749.tif 000006
【図２】
2009053749.tif 000007
【図３】
2009053749.tif 000008
【図４】
2009053749.tif 000009
【図５】
2009053749.tif 000010
【図６】
2009053749.tif 000011
【図７】
2009053749.tif 000012
【図８】
2009053749.tif 000013
【図９】
2009053749.tif 000014
【図１０】
2009053749.tif 000015
【図１１】
2009053749.tif 000016