JP2011022686A - 室内温度分布モデル作成装置及び室内温度分布モデル作成方法 - Google Patents

室内温度分布モデル作成装置及び室内温度分布モデル作成方法 Download PDF

Info

Publication number
JP2011022686A
JP2011022686A JP2009165430A JP2009165430A JP2011022686A JP 2011022686 A JP2011022686 A JP 2011022686A JP 2009165430 A JP2009165430 A JP 2009165430A JP 2009165430 A JP2009165430 A JP 2009165430A JP 2011022686 A JP2011022686 A JP 2011022686A
Authority
JP
Japan
Prior art keywords
temperature
creating
zone
temperature distribution
thermal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2009165430A
Other languages
English (en)
Inventor
Hiroyasu Okuyama
博康 奥山
Yoshiaki Onishi
由哲 大西
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shimizu Construction Co Ltd
Shimizu Corp
Original Assignee
Shimizu Construction Co Ltd
Shimizu Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shimizu Construction Co Ltd, Shimizu Corp filed Critical Shimizu Construction Co Ltd
Priority to JP2009165430A priority Critical patent/JP2011022686A/ja
Publication of JP2011022686A publication Critical patent/JP2011022686A/ja
Pending legal-status Critical Current

Links

Images

Abstract

【課題】 室内の精密なモデルを作成する室内温度分布モデル作成装置及び室内温度分布モデル作成方法を提供する。
【解決手段】 モデル化する対象のデータを入力するステップと、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンＡを作成するステップと、室内外を画成する外皮が外気に接することで外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンＢを作成するステップと、発熱体ごとに発熱体ゾーンＣを設定し、発熱体から発生する熱による上昇流からなる熱プルームを発熱体ゾーンＣごとに作成するステップと、ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網２を作成するステップと、ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンス３３，３４を求め、熱回路網３を作成するステップと、モデル化した室内の温度分布を出力するステップと、を備えたことを特徴とする。
【選択図】図１

Description

本発明は、室内の温度分布をモデル化し、熱負荷等を検討できる室内温度分布モデル作成装置及び室内温度分布モデル作成方法に関する。
従来、室内の温度分布を予測する手段として、数値流体解析が用いられていた。しかしながら、数値流体解析は、計算負荷が大きく、数値的な安定性の問題等から長時間の現象を予測することが難しかった。
一方、本発明者は、建築物の熱回路網及び換気回路網モデルを考案していた。これにより実用的に熱負荷計算や平均的室温の変化を計算することができた。しかし今までは、室内の温度は一様混合の１ゾーンモデルであった。
そこで、本発明者は室内の上下温度分布を表現できるモデルとして、室空間を上下に複数に分割した温度成層化ゾーン、発熱体から上昇流が発生する部分を上下に複数に分割した熱プルームゾーン、及び外皮の内表面に接して冷却される下降流が発生する部分を上下に複数に分割した温度境界層ゾーンの複数のゾーンを設け、ゾーン間の空気の移動量や熱の伝わりやすさ等を求め、実測の温度を再現する室内の精密なモデルを作成した（特許文献１参照）。
特開２００９−５３７４９号公報
（変更？） 宇田川,石田,石野等,「熱負荷算法小委員会報告書-大空間の熱負荷計算法-」,1993年3月,空衛学会,空気調和設備委員会,熱負荷算法小委員会
ところで、上記特許文献１に記載されたものは、熱プルームを１本に想定していたので、これを加熱する複数の日射吸熱の表面が１本の熱プルームを介して平均化され、一様に近い温度になってしまうことがあった。
本発明は上記課題を解決し、室内の実用的なモデルを作成する室内温度分布モデル作成装置及び室内温度分布モデル作成方法を提供することを目的とする。
本発明は上記課題を解決するものであって、室内の温度分布をモデル化する室内温度分布モデル作成装置において、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンを作成する温度成層化領域作成手段と、室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンを作成する温度境界層領域作成手段と、発熱体ごとに発熱体ゾーンを設定し、前記発熱体から発生する熱による上昇流からなる熱プルームを発熱体ゾーンごとに作成する熱プルーム領域作成手段と、前記ゾーン底面の静圧と前記ゾーン間の空気移動量を求める換気回路網を作成する換気回路網作成手段と、前記ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンスを求め、熱回路網を作成する熱回路網作成手段と、を備えたことを特徴とする。
また、前記熱プルーム領域作成手段は、前記熱プルームの上部を束ねるヘッダーゾーンと、前記熱プルームの下部を束ねるフッターゾーンと、を作成することを特徴とする。
また、前記温度成層化ゾーン間の混合流に関する一般化熱コンダクタンスと、前記換気回路網の流路の圧力損失係数とを、予測値と実測値との誤差の評価関数を最小にするように求めるパラメータ演算手段を有することを特徴とする。
また、前記評価関数は、前記温度成層化ゾーンの各空気温度と前記発熱体の各表面温度に関して、予測値から実測値を差し引いた差の評価期間積分値を、これらの空気と表面温度それぞれで平均化し、さらに空気温度と表面温度の平均温度をとることを特徴とする。
さらに、室内の温度分布をモデル化する室内温度分布モデル作成方法において、モデル化する対象のデータを入力するステップと、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンを作成するステップと、室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンを作成するステップと、発熱体ごとに発熱体ゾーンを設定し、前記発熱体から発生する熱による上昇流からなる熱プルームを発熱体ゾーンごとに作成するステップと、前記ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網を作成するステップと、前記ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンスを求め、熱回路網を作成するステップと、モデル化した室内の温度分布を出力するステップと、を備えたことを特徴とする。
また、前記熱プルームの上部をヘッダー領域で束ねるステップと、前記熱プルームの下部をフッター領域で束ねるステップと、を有することを特徴とする。
また、前記温度成層化ゾーン間の混合流に関する一般化熱コンダクタンスと、前記換気回路網の流路の圧力損失係数とを、予測値と実測値との誤差の評価関数を最小にするように求めるステップを有することを特徴とする。
また、前記評価関数は、前記温度成層化ゾーンの各空気温度と前記発熱体の各表面温度に関して、予測値から実測値を差し引いた差の評価期間積分値を、これらの空気と表面温度それぞれで平均化し、さらに空気温度と表面温度の平均温度をとるステップを有することを特徴とする。
本発明によれば、室内の温度分布をモデル化する室内温度分布モデル作成装置において、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンを作成する温度成層化領域作成手段と、室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンを作成する温度境界層領域作成手段と、発熱体ごとに発熱体ゾーンを設定し、前記発熱体から発生する熱による上昇流からなる熱プルームを発熱体ゾーンごとに作成する熱プルーム領域作成手段と、前記ゾーン底面の静圧と前記ゾーン間の空気移動量を求める換気回路網を作成する換気回路網作成手段と、前記ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンスを求め、熱回路網を作成する熱回路網作成手段と、を備えたので、実用的なモデルを作成することができる。そして、実用的なモデルにより、冷暖房の定常・非定常の熱負荷の予測ができ、空調機器設定の目安が正確に計算できる。また、部材毎の伝熱量が概算できるので、熱負荷低減アイテムの効果が概算できる。さらに、個々の発熱体による温度変化を、対応する熱プルームが表現できる様になり、より実用的で正確なモデルを作成することができる。
また、前記熱プルーム領域作成手段は、前記熱プルームの上部を束ねるヘッダーゾーンと、前記熱プルームの下部を束ねるフッターゾーンと、を作成するので、プルーム間での循環流を抑制し、より実用的なモデルを作成することができる。
また、前記温度成層化ゾーン間の混合流に関する一般化熱コンダクタンスと、前記換気回路網の流路の圧力損失係数とを、予測値と実測値との誤差の評価関数を最小にするように求めるパラメータ演算手段を有するので、最適なパラメータを求めることができ、より正確なモデルができる。
また、前記評価関数は、前記温度成層化ゾーンの各空気温度と前記発熱体の各表面温度に関して、予測値から実測値を差し引いた差の評価期間積分値を、これらの空気と表面温度それぞれで平均化し、さらに空気温度と表面温度の平均温度をとるので、より正確な評価をすることが可能となる。
さらに、本発明は上記課題を解決する方法であって、室内の温度分布をモデル化する室内温度分布モデル作成方法において、モデル化する対象のデータを入力するステップと、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンを作成するステップと、室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンを作成するステップと、発熱体ごとに発熱体ゾーンを設定し、前記発熱体から発生する熱による上昇流からなる熱プルームを発熱体ゾーンごとに作成するステップと、前記ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網を作成するステップと、前記ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンスを求め、熱回路網を作成するステップと、モデル化した室内の温度分布を出力するステップと、を備えたので、実用的にモデルを作成することができる。そして、精密なモデルにより、冷暖房の定常・非定常の熱負荷の予測ができ、空調機器設定の目安が正確に計算できる。また、部材毎の伝熱量が概算できるので、熱負荷低減アイテムの効果が概算できる。さらに、個々の発熱体による温度変化を、対応する熱プルームが表現できる様になり、より正確なモデルを実用的に作成することができる。
また、前記熱プルームの上部をヘッダー領域で束ねるステップと、前記熱プルームの下部をフッター領域で束ねるステップと、を有するので、プルーム間での循環流を抑制し、より適切なモデルを作成することができる。
また、前記温度成層化ゾーン間の混合流に関する一般化熱コンダクタンスと、前記換気回路網の流路の圧力損失係数とを、予測値と実測値との誤差の評価関数を最小にするように求めるステップを有するので、最適なパラメータを求めることができ、より正確なモデルができる。
また、前記評価関数は、前記温度成層化ゾーンの各空気温度と前記発熱体の各表面温度に関して、予測値から実測値を差し引いた差の評価期間積分値を、これらの空気と表面温度それぞれで平均化し、さらに空気温度と表面温度の平均温度をとるステップを有するので、より正確な評価をすることが可能となる。
本実施形態の室内温度分布作成装置が作成するモデルの概要を示す図である。 本実施形態の装置構成を示す図である。 車両を図１に示す領域に分割する例を示す図である。 図３に示す車両に換気回路網を構築しモデル化した図である。 熱プルーム領域の拡大図である。 図３に示す車両に熱回路網を構築しモデル化した図である。 実測時の外気温度と相対湿度を示すグラフである。 法線面直達日射量と水平面天空拡散日射量を示すグラフである。 頭部のモデルによる計算値と実測値のグラフである。 胸部のモデルによる計算値と実測値のグラフである。 膝先のモデルによる計算値と実測値のグラフである。 足元のモデルによる計算値と実測値のグラフである。 ルーフ中央部のモデルによる計算値と実測値のグラフである。 ドア後左部のモデルによる計算値と実測値のグラフである。 右シート背面のモデルによる計算値と実測値のグラフである。 右シート座面のモデルによる計算値と実測値のグラフである。
以下、図面を参照して本発明にかかる実施形態の室内温度分布作成装置１について説明する。
図１は、本実施形態の室内温度分布作成装置１が作成するモデルの概要を示す図である。図１において、Ａは温度成層化領域、Ｂは温度境界層領域、Ｃは熱プルーム領域を示す。また、換気回路網２では、矢印２１は空気の移動を表す換気回路網の流路を示す。さらに、熱回路網３では、節点３１は例えば外気温度等の既知の温度を表す与条件温度節点、節点３２は温度が未知な未知数温度節点を示し、第１熱コンダクタンス３３は伝導・伝達・貫流・放射等の一般化熱コンダクタンス、第２熱コンダクタンス３４は移流・対流による一般化熱コンダクタンスを示し、太矢印３５は節点への発熱量を示す。なお、一般化コンダクタンスとは、伝導、伝達、放射、貫流、移流等すべての伝熱形態の熱の伝わりやすさを表す。
本実施形態では、温度成層化領域Ａ、温度境界層領域Ｂ及び熱プルーム領域Ｃに区分する。また、温度成層化領域Ａ及び温度境界層領域Ｂをそれぞれ上下に層状のゾーンＡ１，Ａ２・・・Ｂ３，Ｂ４に分割して、換気回路網２と熱回路網３を構築する。熱プルーム領域Ｃは、発熱体ごとに熱プルームＣ１，Ｃ２，Ｃ３を構築し、複数の熱プルームＣ１，Ｃ２，Ｃ３の上下端を束ね、上端にヘッダーゾーンＣｈ、下端にフッターゾーンＣｆを設けた。
温度成層化領域Ａは、室内の上下温度分布を複数の層に分割し、離散的に表現したものであり、各層内での温度は一様である。温度成層化領域Ａ間の平均的な流れは上向きまたは下向きの押し出し流れであり、こうしたピストン流れとすることが温度成層化を模擬する要件となる。また、自然対流（熱対流のみ）と強制対流（ファン作動時など）の両状態は、温度成層化領域Ａ内の各ゾーン間の混合流の一般化熱コンダクタンス３３，３４の数値を，小さくまたは大きくすることでモデル化する。
温度境界層領域Ｂは、室内外を画成する外皮が外気に接することで外皮の内表面側に下降流や上昇流を生じるゾーンである。例えば、車両の外気温度が外皮を通じて車内を冷却すると、温度境界層領域Ｂを冷却し、温度境界層領域Ｂ内で下降流が生じ，この下降流が温度成層化領域Ａに入り込み、循環流を引き起こす。
熱プルーム領域Ｃは、発熱体から発生する熱による上昇流を生じるゾーンである。例えば、日射を受けた車室内の物体から発生する上昇流を表現するものであり、熱プルーム領域Ｃは、発熱体ごとに熱プルームＣ１，Ｃ２，Ｃ３を構築し、複数の熱プルームＣ１，Ｃ２，Ｃ３の上下端を束ね、上端にヘッダーゾーンＣｈ、下端にフッターゾーンＣｆを設けた。また、上下方向を位置を設定することで，発熱体との接触位置が高い場合、低い場合の熱上昇力の違いを考慮できる様にした。
なお、ファンなどで強制対流加熱した場合には，室内は空気的に混合状態になることから，温度成層化領域Ａの全てのゾーンに発熱量は均等に分散して与えられる。これでも上下温度分布が生成されるのは押し出し上昇流が生じるからである。また、この強制対流加熱状態では熱プルーム領域Ｃ内には加熱源がなく、また熱プルーム流路の断面積が小さいことから熱プルームの流の影響はなくなる。こうして強制対流と自然対流の両状態は同じモデル構造で表現できる。
ここで、換気回路網２と熱回路網３について説明する（詳細については、非特許文献１を参照）。
換気回路網２は、各ゾーン底面の静圧と複数ゾーン間の流量を示すものである。
前提として、各ゾーンでは風量の収支が成り立つとする。また、各ゾーン底面に静圧ｐ(Pa)を持ち、ゾーンの空気密度をρ(kg/m3)とすると、この底面から上方へｈ(m)上がった点での静圧はｐ−ρ・ｈ・ｇ(Pa)とする。空気密度はそのゾーンの空気温度の関数である。流路は、圧力損失係数、指数及び流路面積、さらに基準面からの高さの属性を持ち、流路の通過風量は，前後の静圧差と圧力損失係数、指数および流路面積の関数である。そして、ゾーン数分の静圧に関する非線形連立方程式を修正ニュートンラプソン法で解くことで風量が求められる。
熱回路網３は、熱の移動を示すものである。前提として、各節点では、流入熱流から流出熱流を引いたものが、熱容量に温度変化を乗じたものに等しくなるという熱流収支が成り立つものとする。これを完全連結システムの節点方程式とする。また、伝導、表面伝達、長波長放射、移流など全ての伝熱形態の熱の伝わりやすさを一般化熱コンダクタンス３３，３４という一種類の係数で表現する。従って空気流動による一般化熱コンダクタンス３３，３４も数式記号上は同じ係数となる。
そして、全節点の温度に関する連立常微分方程式が前述の節点方程式から構成され、連立常微分方程式を後退差分などの時間積分法によって解き、毎時間ステップの温度の解を算出する。
次に、換気回路網２と熱回路網３の連成について説明する。換気回路網２の駆動条件となる各ゾーンの密度は熱回路網３の解の温度によって計算される。熱回路網３の係数となる風量は換気回路網２の解から得られる。非定常の時間ステップ毎に、一方の回路網への条件はもう一方の回路網の前時間ステップでの結果をもとにつくる。こうして両モデルの連成がとられる。
図２は本実施形態の室内温度分布作成装置１を用いたシステム構成図である。
入力手段１０１は、モデル化する対象のデータを入力するものである。記憶手段１０２は、演算式等の不変のデータを記憶するものである。入力手段１０１から入力されたデータ及び記憶手段１０２に記憶されたデータは、制御手段１１０に入力される。
室内温度分布作成装置１の制御手段１１０は、温度成層化領域作成手段１１１、温度境界層領域作成手段１１２、熱プルーム領域作成手段１１３、換気回路網作成手段１１４及び熱回路網作成手段１１５、パラメータ演算手段１１６を有する。
温度成層化領域作成手段１１１は、室内を複数のゾーンに分割し、上下温度分布を離散的に表現する温度成層化ゾーンを作成する。温度境界層領域作成手段１１２は、室内外を画成する外皮が外気に接することで外皮の内表面側に下降流や上昇流を生じるゾーンを作成する。なお、温度境界層ゾーンは、複数のゾーンに分割し作成されてもよい。熱プルーム領域作成手段１１３は、発熱体から発生する熱による上昇流を生じるゾーンを作成する。なお、熱プルームゾーンは、複数のゾーンに分割し作成されてもよい。換気回路網作成手段１１４は、ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網を作成する。熱回路網作成手段１１５は、ゾーン間の熱の移動を求める熱回路網を作成する。パラメータ演算手段１１６については後述する。
制御手段１１０により作成されたモデルは、表示装置等の出力手段１２１に出力される。
図３は、車両５０を図１に示す領域に分割する例を示す図である。車両５０の車室５１内を上下４段の温度成層化領域Ａに区分し、温度境界層領域Ｂとして、屋根温度境界層Ｂ１、窓温度境界層Ｂ２、ドア温度境界層Ｂ３及び床温度境界層Ｂ４を設定したものである。
図４は車両５０に換気回路網２を構築しモデル化した図、図５は熱プルーム領域Ｃの換気回路網２の拡大図である。
図５に示すように、本実施例では、熱プルームを合計１４本有し、各熱プルームは各日射発熱体に対応している。Ｃ１ａは左後部ピラー、Ｃ１ｂは右後部ピラー、Ｃ２ａはインスツルメントパネル、Ｃ２ｂは左ドア内側、Ｃ２ｃは右ドア内側、Ｃ２ｄはリアトレイ、Ｃ２ｅは後席左側背面、Ｃ２ｆは後席右側背面、Ｃ２ｇは左前席背面、Ｃ２ｈは右前席背面、Ｃ３ａは後席左側座面、Ｃ３ｂは後席右側座面、Ｃ３ｃは左前席座面、Ｃ３ｄは右前席座面となっている。
各発熱体に接する空気の発熱ゾーンの高さは発熱体の高さに対応し、熱プルームの上昇力は発熱ゾーンの高さに関係する。なお、発熱ゾーンの高さとはゾーンの底面高さである。これらには発熱体の実際の高さを設定している。また、流路の高さとは，その流路の出口側のゾーンに流れ込む位置での高さを意味する。流路内の空気温度は風上側のゾーンの温度に近いと考えるからである。
次に、作成するモデルのパラメータについて説明する。作成するモデルは、工学モデルなので、実験的に定めなければならない特性値がある。これらは、換気回路網２の圧力損失係数ζ及び熱回路網３の層間混合流一般化熱コンダクタンスｃijからなる二種のパラメータである。
これらパラメータをパラメータ演算手段１１６により最適化する方法について説明する。まず、パラメータを最適化するために評価関数を設定する。式（１）は、上下の平均温度は、室空気が持つ熱エネルギーに比例するので、室内の総エネルギーが一致すべきという観点から導かれる評価関数である。
Figure 2011022686
ただし、pθjは頭部又は足元の予測計算温度、
mθjは実測温度、
Ｔは検討期間
である。
ここにｊは成層化ゾーンの上下方向の位置番号を示し，一般に２以上あるものとする．
次に、頭部と足元それぞれで温度が近ければ良いという各層での温度が一致するかの評価を式（２）で行う。
Figure 2011022686
これら式（１）と式（２）の両方を平均し、総合的に考慮する式（３）の誤差評価関数Ｊを定める。
Figure 2011022686
次に、日射発熱体の表面温度が一致することの評価関数は式（４）で定める。日射発熱面は１４面存在するので、これらの表面温度は受熱体表面積で重み付け平均した面積加重平均として評価関数とするのが妥当である。
Figure 2011022686
 ただし、Ｓｊは受熱体の表面積である。
以上の空気温度と表面温度の誤差評価関数の平均の温度である式（５）を最終的な評価関数とする。
Figure 2011022686
次に、未定パラメータである換気回路網の圧力損失係数ζと温度成層化ゾーン間の混合流熱コンダクタンスｃijを求める。まず、空気温度についての評価関数である式（３）を調べる。
圧力損失係数ζは１，１２，１００の３通りについて、層間混合流一般化熱コンダクタンスｃijは０．０３５，０．３５，３．５の３通りについて計算した。
Figure 2011022686
次に、同様なパラメータの組み合わせにより、発熱体表面温度の評価関数式(４)を計算する。
Figure 2011022686
そして、空気温度と表面温度の平均値を表す式（５）について計算する。
Figure 2011022686
以上のような検討により、最適なパラメータを表（４）の様に定める。ここに、対流熱伝達率は予測表面温度が実測値に近付く様に調整したものである。
Figure 2011022686
図６は、車両５０に熱回路網３を構築しモデル化した図である。本実施例のモデルについては、室内表面の放射伝達と表面対流伝達を分離したモデルなので、総合伝達率のモデルと比較して、室内部位間の長波長放射伝達の様子を描いた分だけ複雑になっている。
次に、このような室内温度分布作成装置１により作成したモデルによる計算値と実測値との比較結果について説明する。図７は実測時の外気温度と相対湿度を示すグラフ、図８は法線面直達日射量と水平面天空拡散日射量を示すグラフである。
図７において、薄線は外気温度、濃線は相対湿度である。１日目及び２日目は、日中に、高温且つ低湿度であり、３日目及び４日目は、日中に、気温が上がらず、高湿度であった。
また、図８において、濃線は法線面直達日射量、薄線は水平面天空拡散日射量を示す。１日目及び２日目は、日中に、法線面直達日射量が多く、３日目及び４日目は、日中に、法線面直達日射量が少ない。
図７及び図８から、１日目及び２日目は、日中に晴天であり、３日目及び４日目は、日中に曇天又は雨天であったことがわかる。
図９から図１２は、４日分の室内の空気についてモデルによる計算値と実測値とで比較したものである。
図９は頭部のモデルによる計算値と実測値のグラフ、図１０は胸部のモデルによる計算値と実測値のグラフ、図１１は膝先のモデルによる計算値と実測値のグラフ、図１２は足元のモデルによる計算値と実測値のグラフである。
図９から図１２のグラフをみると、４日間の全期間について、頭部、胸部、膝先、足元の４層での空気温度に対して、モデルによる計算値は、実測値に近い結果が得られている。
図１３から図１６は、１日分の表面温度についてモデルによる計算値と実測値とで比較したものである。
図１３はルーフ中央部のモデルによる計算値と実測値のグラフ、図１４はドア後左部のモデルによる計算値と実測値のグラフ、図１５は右シート背面のモデルによる計算値と実測値のグラフ、図１６は右シート座面のモデルによる計算値と実測値のグラフである。
図１３から図１６のグラフをみると、１日間の全期間について、ルーフ、ドア、シート背面、シート座面の表面温度に対して、モデルによる計算値は、実測値に近い結果が得られている。これは、熱プルームを発熱体ごとに複数設けたことにより、個々の発熱体による温度変化を、対応する熱プルームが表現できる様になったからであると考えられる。
このように、本実施形態によれば、室内の温度分布をモデル化する室内温度分布モデル作成装置１において、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンＡ１，Ａ２，Ａ３，Ａ４を作成する温度成層化領域作成手段１１１と、室内外を画成する外皮が外気に接することで外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンＢ１，Ｂ２，Ｂ３，Ｂ４を作成する温度境界層領域作成手段１１２と、発熱体ごとに発熱体ゾーンＣ１，Ｃ２，Ｃ３を設定し、発熱体から発生する熱による上昇流からなる熱プルームを発熱体ゾーンＣ１，Ｃ２，Ｃ３ごとに作成する熱プルーム領域作成手段１１３と、ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網２を作成する換気回路網作成手段１１４と、ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンス３３，３４を求め、熱回路網３を作成する熱回路網作成手段１１５と、を備えたので、実用的なモデルを作成することができる。そして、実用的なモデルにより、冷暖房の定常・非定常の熱負荷の予測ができ、空調機器設定の目安が正確に計算できる。また、部材毎の伝熱量が概算できるので、熱負荷低減アイテムの効果が概算できる。さらに、個々の発熱体による温度変化を、対応する熱プルームが表現できる様になり、実用的で正確なモデルを作成することができる。
また、熱プルーム領域作成手段１１３は、熱プルームの上部を束ねるヘッダーゾーンＣｈと、熱プルームの下部を束ねるフッターゾーンＣｆと、を作成するので、プルーム間での循環流を抑制し、より実用的なモデルを作成することができる。
また、温度成層化ゾーンＡ１，Ａ２，Ａ３，Ａ４間の一般化熱コンダクタンス３３，３４と、換気回路網２の流路の圧力損失係数ζとを、予測値と実測値との誤差の評価関数を最小にするように求めるパラメータ演算手段１１６を有するので、最適なパラメータを求めることができ、より正確なモデルができる。
さらに、本実施形態は、室内の温度分布をモデル化する室内温度分布モデル作成方法において、モデル化する対象のデータを入力するステップと、室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンＡ１，Ａ２，Ａ３，Ａ４を作成するステップと、室内外を画成する外皮が外気に接することで外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンＢ１，Ｂ２，Ｂ３，Ｂ４を作成するステップと、発熱体ごとに発熱体ゾーンＣ１，Ｃ２，Ｃ３を設定し、発熱体から発生する熱による上昇流からなる熱プルームを発熱体ゾーンＣ１，Ｃ２，Ｃ３ごとに作成するステップと、ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網２を作成するステップと、ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンス３３，３４を求め、熱回路網３を作成するステップと、モデル化した室内の温度分布を出力するステップと、を備えたので、実用的にモデルを作成することができる。そして、精密なモデルにより、冷暖房の定常・非定常の熱負荷の予測ができ、空調機器設定の目安が正確に計算できる。また、部材毎の伝熱量が概算できるので、熱負荷低減アイテムの効果が概算できる。さらに、個々の発熱体による温度変化を、対応する熱プルームが表現できる様になり、より正確なモデルを実用的に作成することができる。
また、熱プルームの上部をヘッダー領域Ｃｈで束ねるステップと、熱プルームの下部をフッター領域Ｃｆで束ねるステップと、を有するので、プルーム間での循環流を抑制し、より適切なモデルを作成することができる。
また、温度成層化ゾーンＡ１，Ａ２，Ａ３，Ａ４間の一般化熱コンダクタンス３３，３４と、換気回路網２の流路の圧力損失係数ζとを、予測値と実測値との誤差の評価関数を最小にするように求めるステップを有するので、最適なパラメータを求めることができ、より正確なモデルができる。
なお、本実施形態では、車両５０のモデル作成について説明したが、これに限らず、建築物や構造物等の室内に適用してもよい。
１…室内温度分布作成装置、２…換気回路網、２１…流路、３…熱回路網、３１…与条件温度節点、３２…未知数温度節点、３３…第１熱コンダクタンス、３４…第２熱コンダクタンス、３５…節点への発熱量、３６…熱橋熱コンダクタンス、５０…車両、５１…車室、５２…熱橋、１０１…入力手段、１１０…制御手段、１１１…温度成層化領域作成手段、１１２…温度境界層領域作成手段、１１３…熱プルーム領域作成手段、１１４…換気回路網作成手段、１１５…熱回路網作成手段、１１６…パラメータ演算手段、１２１出力手段、Ａ…温度成層化領域、Ｂ…温度境界層領域、Ｃ…熱プルーム領域、Ｃｈ…ヘッダーゾーン、Ｃｆ…フッターゾーン、Ｃ１，Ｃ２，Ｃ３…発熱体ゾーン

Claims (8)

  1. 室内の温度分布をモデル化する室内温度分布モデル作成装置において、
    室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンを作成する温度成層化領域作成手段と、
    室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンを作成する温度境界層領域作成手段と、
    発熱体ごとに発熱体ゾーンを設定し、前記発熱体から発生する熱による上昇流からなる熱プルームを発熱体ゾーンごとに作成する熱プルーム領域作成手段と、
    前記ゾーン底面の静圧と前記ゾーン間の空気移動量を求める換気回路網を作成する換気回路網作成手段と、
    前記ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンスを求め、熱回路網を作成する熱回路網作成手段と、
    を備えたことを特徴とする室内温度分布モデル作成装置。
  2. 前記熱プルーム領域作成手段は、
    前記熱プルームの上部を束ねるヘッダーゾーンと、
    前記熱プルームの下部を束ねるフッターゾーンと、
    を作成することを特徴とする請求項１に記載の室内温度分布モデル作成装置。
  3. 前記温度成層化ゾーン間の混合流に関する一般化熱コンダクタンスと、前記換気回路網の流路の圧力損失係数とを、予測値と実測値との誤差の評価関数を最小にするように求めるパラメータ演算手段を有することを特徴とする請求項１又は請求項２に記載の室内温度分布モデル作成装置。
  4. 前記評価関数は、前記温度成層化ゾーンの各空気温度と前記発熱体の各表面温度に関して、予測値から実測値を差し引いた差の評価期間積分値を、これらの空気と表面温度それぞれで平均化し、さらに空気温度と表面温度の平均温度をとることを特徴とする請求項３に記載の室内温度分布モデル作成装置。
  5. 室内の温度分布をモデル化する室内温度分布モデル作成方法において、
    モデル化する対象のデータを入力するステップと、
    室内を複数のゾーンに分割し、上下温度分布を表現する温度成層化ゾーンを作成するステップと、
    室内外を画成する外皮が外気に接することで前記外皮の内表面側に下降流や上昇流を生じる温度境界層ゾーンを作成するステップと、
    発熱体ごとに発熱体ゾーンを設定し、前記発熱体から発生する熱による上昇流からなる熱プルームを発熱体ゾーンごとに作成するステップと、
    前記ゾーン底面の静圧とゾーン間の空気移動量を求める換気回路網を作成するステップと、
    前記ゾーン間の熱の伝わりやすさを表す一般化熱コンダクタンスを求め、熱回路網を作成するステップと、
    モデル化した室内の温度分布を出力するステップと、
    を備えたことを特徴とする室内温度分布モデル作成方法。
  6. 前記熱プルームの上部をヘッダー領域で束ねるステップと、
    前記熱プルームの下部をフッター領域で束ねるステップと、
    を有することを特徴とする請求項５に記載の室内温度分布モデル作成方法。
  7. 前記温度成層化ゾーン間の混合流に関する一般化熱コンダクタンスと、前記換気回路網の流路の圧力損失係数とを、予測値と実測値との誤差の評価関数を最小にするように求めるステップを有することを特徴とする請求項５又は請求項６に記載の室内温度分布モデル作成方法。
  8. 前記評価関数は、前記温度成層化ゾーンの各空気温度と前記発熱体の各表面温度に関して、予測値から実測値を差し引いた差の評価期間積分値を、これらの空気と表面温度それぞれで平均化し、さらに空気温度と表面温度の平均温度をとるステップを有することを特徴とする請求項７に記載の室内温度分布モデル作成方法。
JP2009165430A 2009-07-14 2009-07-14 室内温度分布モデル作成装置及び室内温度分布モデル作成方法 Pending JP2011022686A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009165430A JP2011022686A (ja) 2009-07-14 2009-07-14 室内温度分布モデル作成装置及び室内温度分布モデル作成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009165430A JP2011022686A (ja) 2009-07-14 2009-07-14 室内温度分布モデル作成装置及び室内温度分布モデル作成方法

Publications (1)

Publication Number Publication Date
JP2011022686A true JP2011022686A (ja) 2011-02-03

Family

ID=43632725

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009165430A Pending JP2011022686A (ja) 2009-07-14 2009-07-14 室内温度分布モデル作成装置及び室内温度分布モデル作成方法

Country Status (1)

Country Link
JP (1) JP2011022686A (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102012201538A1 (de) 2011-02-04 2012-08-09 Fuji Jukogyo K.K. Antriebsvorrichtung für ein Fahrzeug

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102012201538A1 (de) 2011-02-04 2012-08-09 Fuji Jukogyo K.K. Antriebsvorrichtung für ein Fahrzeug

Similar Documents

Publication Publication Date Title
Wang et al. Advances in building simulation and computational techniques: A review between 1987 and 2014
Hajdukiewicz et al. Formal calibration methodology for CFD models of naturally ventilated indoor environments
Cholewa et al. On the heat transfer coefficients between heated/cooled radiant floor and room
Khaoua et al. Analysis of greenhouse ventilation efficiency based on computational fluid dynamics
Molina-Aiz et al. Comparison of finite element and finite volume methods for simulation of natural ventilation in greenhouses
Villar et al. Numerical 3-D heat flux simulations on flat plate solar collectors
Fuliotto et al. Experimental and numerical analysis of heat transfer and airflow on an interactive building facade
Teodosiu et al. CFD modeling of buoyancy driven cavities with internal heat source—Application to heated rooms
Guo et al. Simulation and measurement of air temperatures and mean radiant temperatures in a radiantly heated indoor space
Qi et al. The effects of non-uniform temperature distribution on neutral plane level in non-adiabatic high-rise shafts during fires
Gan General expressions for the calculation of air flow and heat transfer rates in tall ventilation cavities
Taieb et al. Natural convection in an asymmetrically heated vertical channel with an adiabatic auxiliary plate
Mahyuddin et al. Computational fluid dynamics modelling of the air movement in an environmental test chamber with a respiring manikin
Larsen et al. Calculation methods for single-sided natural ventilation: Now and ahead
Singh et al. A literature review of building energy simulation and computational fluid dynamics co-simulation strategies and its implications on the accuracy of energy predictions
Wurtz et al. An equation-based simulation environment to investigate fast building simulation
Mullen et al. Automatic extraction of reduced-order models from CFD simulations for building energy modelling
Saraz et al. A CFD approach to assess the effects of different opening combinations in poultry houses
JP2011022686A (ja) 室内温度分布モデル作成装置及び室内温度分布モデル作成方法
Sun et al. Investigation of the thermal characteristics of condensers in nuclear power plant by simulation with zoning model
Stavridou et al. Unsteady CFD simulation in a naturally ventilated room with a localized heat source
Heidarinejad et al. Urban scale modeling of campus building using Virtual PULSE
Desta et al. Combining CFD and data-based mechanistic (DBM) modelling approaches
Boulet et al. Sequential inverse method implemented into CFD software for the estimation of a radiation boundary
CN115344919A (zh) 建筑内部环境分析方法、装置、设备及介质