(57)【要約】【課題】三次元空間における熱移動を考慮した高精度の温度予測を、顧客の目の前で計算可能な程度の計算負荷で実現した、断熱容器内の温度算出装置を提供する。【解決手段】蓄熱材２１と断熱容器２０とを含む解析対象となる空間をＸ軸、Ｙ軸及びＺ軸に沿って複数のメッシュに区画した場合に、温度算出対象メッシュＭ’を通るＸ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群それぞれに、メッシュを占める物質に応じた伝熱式を設定した伝熱式データ１４を記憶する記憶部１１と、予め設定された解析条件と伝熱式データを用い、Ｘ軸、Ｙ軸及びＺ軸上のメッシュ群に沿って温度算出対象メッシュＭ’に対して出入りする単位時間あたりの熱量総和Ｑを算出する熱移動量算出部１２と、熱量総和Ｑを用いて、温度算出対象メッシュＭ’の単位時間経過後の温度を算出する温度算出部１３と、を有する。【選択図】図１
2015187827.tif 000002
【特許請求の範囲】
【請求項１】
蓄熱材を収容した断熱容器内温度を算出する装置であって、
前記蓄熱材と前記断熱容器とを含む解析対象となる空間をＸ軸、Ｙ軸及びＺ軸に沿って複数のメッシュに区画した場合に、温度算出対象となるメッシュを通るＸ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群それぞれに、メッシュを占める物質に応じた伝熱式を設定した伝熱式データを記憶する記憶部と、
予め設定された解析条件と前記伝熱式データを用い、前記Ｘ軸上のメッシュ群、前記Ｙ軸上のメッシュ群及び前記Ｚ軸上のメッシュ群に沿って前記温度算出対象となるメッシュに対して出入りする単位時間あたりの熱量総和を算出する熱移動量算出部と、
前記熱量総和を用いて、前記温度算出対象となるメッシュの単位時間経過後の温度を算出する温度算出部と、
を備える断熱容器内の温度算出装置。
【請求項２】
前記伝熱式には、熱伝導及び境膜熱伝達に関する式を含む請求項１に記載の装置。
【請求項３】
前記伝熱式には、放射伝熱に関する式を更に含む請求項２に記載の装置。
【請求項４】
前記蓄熱材の比熱を温度毎に関連付けた温度依存比熱データを記憶する比熱記憶部を備え、
前記蓄熱材に対応するメッシュでの熱量計算には、前記比熱記憶部に記憶されている前記温度依存比熱データのうち現時点の蓄熱材の温度に対応する比熱を読みだして用いる、請求項１〜３のいずれかに記載の装置。
【請求項５】
前記断熱容器内において前記蓄熱材を除く収容用の内部空間が前記温度算出対象としての１つのメッシュに対応するように、前記解析対象となる空間が複数のメッシュに区画されている請求項１〜４のいずれかに記載の装置。
【請求項６】
前記温度算出対象となるメッシュを通るＸ軸、Ｙ軸及びＺ軸について、軸に沿った物質配置パターンは、各軸において単一のパターンが設定されている請求項５に記載の装置。
【請求項７】
前記温度算出対象となるメッシュを通るＸ軸、Ｙ軸及びＺ軸のうち、少なくともいずれかの軸に沿った物質配列パターンが複数存在する場合に、各々の物質配列パターンに対応する伝熱式データが前記記憶部に記憶されており、
前記熱移動量算出部は、前記各々の物質配列パターンに対応する伝熱式データと、前記各々の物質配列パターンが前記温度算出対象となるメッシュに占める面積とに基づき、各々の物質配列パターン毎に前記温度対象となるメッシュに出入りする単位時間あたりの熱量を算出し、各計算結果を合計することで前記温度対象となるメッシュに出入りする単位時間あたりの熱量総和を算出する請求項５に記載の装置。
【請求項８】
前記温度算出対象となるメッシュを高さ方向に積層する複数の小領域に分割した場合に、前記小領域間の単位時間あたりの対流熱量に関する対流熱式データが前記記憶部に記憶されており、
前記熱移動量算出部により算出される温度算出対象となるメッシュに出入りする単位時間あたりの熱量総和と前記対流熱式データとに基づき、前記各々の小領域に出入りする単位時間あたりの熱量総和を算出する小領域熱量算出部と、
前記小領域熱量算出部により算出された前記熱量総和に基づき単位時間経過後の小領域の温度を算出する小領域温度算出部と、を備える請求項５〜７のいずれかに記載の装置。
【請求項９】
蓄熱材を収容した断熱容器内温度をコンピュータが算出する方法であって、
前記蓄熱材と前記断熱容器とを含む解析対象となる空間をＸ軸、Ｙ軸及びＺ軸に沿って複数のメッシュに区画した場合に、温度算出対象となるメッシュを通るＸ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群それぞれに、メッシュを占める物質に応じた伝熱式を設定した伝熱式データを記憶部から読み出すステップと、
予め設定された解析条件と前記伝熱式データを用い、前記Ｘ軸上のメッシュ群、前記Ｙ軸上のメッシュ群及び前記Ｚ軸上のメッシュ群に沿って前記温度算出対象となるメッシュに対して出入りする単位時間あたりの熱量総和を算出するステップと、
前記熱量総和を用いて、前記温度算出対象となるメッシュの単位時間経過後の温度を算出するステップと、
を含む、断熱容器内の温度算出方法。
【請求項１０】
前記伝熱式には、熱伝導及び境膜熱伝達に関する式を含む請求項９に記載の方法。
【請求項１１】
前記伝熱式には、放射伝熱に関する式を更に含む請求項１０に記載の方法。
【請求項１２】
前記蓄熱材に対応するメッシュの熱量を計算する際に、前記蓄熱材の比熱を温度毎に関連付けた温度依存比熱データのうち、現時点の蓄熱材の温度に対応する比熱を読みだして熱量計算に用いる、請求項９〜１１のいずれかに記載の方法。
【請求項１３】
前記断熱容器内において前記蓄熱材を除く収容用の内部空間が前記温度算出対象としての１つのメッシュに対応するように、前記解析対象となる空間が複数のメッシュに区画されている請求項９〜１２のいずれかに記載の方法。
【請求項１４】
前記温度算出対象となるメッシュを通るＸ軸、Ｙ軸及びＺ軸について、軸に沿った物質配置パターンは、各軸において単一のパターンが設定されている請求項１３に記載の方法。
【請求項１５】
前記温度算出対象となるメッシュを通るＸ軸、Ｙ軸及びＺ軸のうち、少なくともいずれかの軸に沿った物質配列パターンが複数存在する場合に、各々の物質配列パターンに対応する伝熱式データを記憶部から読み出すステップと、
前記各々の物質配列パターンに対応する伝熱式データと、前記各々の物質配列パターンが前記温度算出対象となるメッシュに占める面積とに基づき、各々の物質配列パターン毎に前記温度対象となるメッシュに出入りする単位時間あたりの熱量を算出するステップと、
各計算結果を合計することで前記温度対象となるメッシュに出入りする単位時間あたりの熱量総和を算出するステップと、を含む請求項１３に記載の方法。
【請求項１６】
前記温度算出対象となるメッシュを高さ方向に積層する複数の小領域に分割した場合に、前記小領域間の単位時間あたりの対流熱量に関する対流熱式データを記憶部から読み出すステップと、
前記温度算出対象となるメッシュに出入りする単位時間あたりの熱量総和と前記対流熱式データとに基づき、前記各々の小領域に出入りする単位時間あたりの熱量総和を算出するステップと、
前記各々の小領域に出入りする単位時間あたりの熱量総和に基づき単位時間経過後の小領域の温度を算出するステップと、を含む請求項１３〜１６のいずれかに記載の方法。
【請求項１７】
請求項９〜１６のいずれかに記載の方法をコンピュータに実行させるコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、断熱容器に収容された保冷剤等の蓄熱材によって断熱容器内の温度がどの程度の時間、所望温度範囲内に保持されるかを予測するために利用可能な、断熱容器内の温度算出装置、方法及びコンピュータプログラムに関する。
【背景技術】
【０００２】
食品、医薬品等を梱包して保管・輸送する際に、内容物を一定温度に保持することが求められる。輸送容器として用いられる断熱容器の設計は、内容物の保持温度や時間によって、容器自体の断熱性能、使用する蓄熱材の種類、個数、配置を変更する必要がある。顧客のニーズに対応するためには、断熱容器及び蓄熱材の試作品を実際に作成し、実際に温度を計測することで、容器性能を評価しており、時間とコストがかかる。このため、断熱容器は、短時間で製品化するのが非常に困難な製品である。
【０００３】
近年では、ＣＡＥ（Computer Aided Engineering）システムの一種として、コンピュータ上で断熱容器内温度の経時変化をシミュレーションするための設計システム及びコンピュータプログラムが開発されている。例えば、特許文献１には、断熱容器内部の二次元平面について温度予測を行い、容器内部の食品の鮮度評価を行う設計システムが記載されている。しかしながら、このシステムは、二次元平面での温度予測しかできないので、計算負荷が少ないものの、断熱容器内に配置される蓄熱材の種類、個数、配置を考慮した正確な温度予測を行うことはできない。
【０００４】
正確な温度予測を行うためには、三次元空間における熱移動量を算出する必要がある。そのためには、有限要素法や有限体積法といった解析手法を用い、断熱容器及び蓄熱材を含む空間を三次元の微小メッシュに区画し、メッシュ毎に熱量計算を行い、熱量に基づき温度変化を算出することが一つの有効な手段として考えられる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００２−１０８９５１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、上記有限要素法や有限体積法といった解析手法では、計算負荷が莫大となり、結果が得られるまで数時間〜数日必要になってしまう。これでは、営業ツールとして顧客の目の前で計算結果を出して検証することはできない。
【０００７】
本発明は、このような課題に着目してなされたものであって、その目的は、三次元空間における熱移動を考慮した高精度の温度予測を、顧客の目の前で計算可能な程度の計算負荷で実現した、断熱容器内の温度算出装置、方法及びコンピュータプログラムを提供することである。
【課題を解決するための手段】
【０００８】
本発明は、上記目的を達成するために、次のような手段を講じている。
【０００９】
すなわち、本発明に係る断熱容器内の温度算出装置は、蓄熱材を収容した断熱容器内温度を算出する装置であって、前記蓄熱材と前記断熱容器とを含む解析対象となる空間をＸ軸、Ｙ軸及びＺ軸に沿って複数のメッシュに区画した場合に、温度算出対象となるメッシュを通るＸ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群それぞれに、メッシュを占める物質に応じた伝熱式を設定した伝熱式データを記憶する記憶部と、予め設定された解析条件と前記伝熱式データを用い、前記Ｘ軸上のメッシュ群、前記Ｙ軸上のメッシュ群及び前記Ｚ軸上のメッシュ群に沿って前記温度算出対象となるメッシュに対して出入りする単位時間あたりの熱量総和を算出する熱移動量算出部と、前記熱量総和を用いて、前記温度算出対象となるメッシュの単位時間経過後の温度を算出する温度算出部と、を備える。
【００１０】
本発明に係る断熱容器内の温度算出方法は、蓄熱材を収容した断熱容器内温度をコンピュータが算出する方法であって、前記蓄熱材と前記断熱容器とを含む解析対象となる空間をＸ軸、Ｙ軸及びＺ軸に沿って複数のメッシュに区画した場合に、温度算出対象となるメッシュを通るＸ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群それぞれに、メッシュを占める物質に応じた伝熱式を設定した伝熱式データを記憶部から読み出すステップと、予め設定された解析条件と前記伝熱式データを用い、前記Ｘ軸上のメッシュ群、前記Ｙ軸上のメッシュ群及び前記Ｚ軸上のメッシュ群に沿って前記温度算出対象となるメッシュに対して出入りする単位時間あたりの熱量総和を算出するステップと、前記熱量総和を用いて、前記温度算出対象となるメッシュの単位時間経過後の温度を算出するステップと、を含む。
【００１１】
このような装置構成又は方法によれば、熱量計算が、温度算出対象となるメッシュを通るＸ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群に限られるので、解析空間にある全メッシュについて計算する場合に比べて、著しく計算コストを低減することが可能となる。それでいて、温度算出対象となるメッシュに影響を与える三次元方向の熱量を計算しているので、二次元平面だけの算出結果に比して高精度な結果を得ることが可能となる。
【００１２】
本発明の具体例としては、前記伝熱式には、熱伝導及び境膜熱伝達に関する式を含むことが好ましい。
【００１３】
精度を向上させるためには、前記伝熱式には、放射伝熱に関する式を更に含むことが好ましい。
【００１４】
本発明の装置において、蓄熱材の比熱の温度依存性を考慮するためには、前記蓄熱材の比熱を温度毎に関連付けた温度依存比熱データを記憶する比熱記憶部を備え、前記蓄熱材に対応するメッシュでの熱量計算には、前記比熱記憶部に記憶されている前記温度依存比熱データのうち現時点の蓄熱材の温度に対応する比熱を読みだして用いるように構成することが好ましい。
【００１５】
本発明の方法において、前記蓄熱材に対応するメッシュの熱量を計算する際に、前記蓄熱材の比熱を温度毎に関連付けた温度依存比熱データのうち、現時点の蓄熱材の温度に対応する比熱を読みだして熱量計算に用いることが好ましい。
【図面の簡単な説明】
【００１６】
【図１】本発明の第１実施形態における断熱容器内の温度算出装置を模式的に示すブロック図。
【図２】解析対象となる断熱容器及び蓄熱材の配置の一例を示す斜視図及び断面図。
【図３Ａ】従来解析方法において計算対象となるメッシュ及び熱移動方向を示す斜視図。
【図３Ｂ】本発明において計算対象となるメッシュ及び熱移動方向を示す斜視図。
【図４】解析空間における各物質とメッシュの対応関係、並びに温度分布の一例を示す図。
【図５】Ｘ軸、Ｙ軸及びＺ軸に沿って温度算出対象メッシュに対して出入りする熱量に関する説明図。
【図６Ａ】外気と断熱材間の境界部分の熱計算式に関する説明図。
【図６Ｂ】断熱材内部又は蓄熱材内部の熱計算式に関する説明図。
【図６Ｃ】断熱材と蓄熱材の間、又は、蓄熱材同士の間の境界部分の熱計算式に関する説明図。
【図７】比熱記憶部に記憶される温度依存比熱データに関する説明図。
【図８】本発明の装置が実行する温度算出処理ルーチンを示すフローチャート。
【図９】第１実施形態において、断熱容器、蓄熱材及び内部空間とこれらに対応するメッシュを示す平面図、側面図。
【図１０Ａ】第２実施形態における断熱容器、蓄熱材及び内部空間を示す平面図、右側面図。
【図１０Ｂ】第２実施形態において、各々の物質配列パターンがメッシュに占める面積を示す説明図。
【図１１】第３実施形態における断熱容器内の温度算出装置を模式的に示すブロック図。
【図１２Ａ】温度算出対象となるメッシュを複数の小領域に分割した様子を示す説明図。
【図１２Ｂ】図１２Ａに示す小領域における周辺壁面の定義に関する説明図。
【図１３】第３実施形態の装置が実行する内部空間温度算出処理ルーチンを示すフローチャート。
【発明を実施するための形態】
【００１７】
以下、本発明の一実施形態について、図面を参照して説明する。
【００１８】
［第１実施形態］
［温度算出装置］
本実施形態の温度検出装置は、断熱容器内部温度をどれぐらいの間、所望温度範囲内に保持できるかをシミュレーションするために、保冷剤などの蓄熱材を収容した断熱容器の内部温度を算出する装置である。具体的には、断熱容器の大きさ、断熱容器の内外温度、蓄熱材の数とその配置、各物質の物性値、外気温の経時変化などの解析に必要な解析パラメータを入力値として予め設定しておき、断熱容器の外気温度及び蓄熱材の保温能力を考慮して単位時間毎の断熱容器内部の温度の計算を繰り返す。
【００１９】
図１に示すように、装置１は、設定部１０と、記憶部１１と、熱移動量算出部１２と、温度算出部１３と、を有する。これら各部１０〜１３は、ＣＰＵ、メモリ、各種インターフェイスを備えたパソコン等の情報処理装置においてＣＰＵが予め記憶されている図８に示す温度算出処理ルーチンを実行することによりソフトウェア及びハードウェアが協働して実現される。
【００２０】
設定部１０は、キーボードやマウス、タッチパネル等の既知の操作部１５を介してユーザからの操作を受け付け、指定されたファイルから解析パラメータ及び伝熱式データ１４を記憶部１１に記憶させる。解析パラメータには、断熱容器に関する値（形状、物性、配置、初期温度）、蓄熱材に関する値（形状、物性、配置、初期温度）、断熱容器の内部の初期温度などが挙げられる。物性には、熱伝達率ｋ[W/(m・℃)]、密度ρ[kg/m3]、比熱Ｃｐ[J/(Kg・℃)]、境膜熱伝達係数ｈ[W/(m・℃)]、放射率εなどが挙げられる。
【００２１】
詳細を説明する前に、本発明の温度算出法の概要について説明する。
【００２２】
図２は、解析対象となる蓄熱材２１と断熱容器２０の一例を示す斜視図及び断面図である。本実施形態では、断面Ａ１と断面Ａ２と断面Ａ３の交差点となる内部空間の温度を算出対象とする。
【００２３】
従来のような有限要素法や有限体積法を用いた温度解析では、図３Ａに示すように、断熱容器及びその周囲の外気を含む解析対象となる空間をＸ軸、Ｙ軸及びＺ軸に沿って複数のメッシュに区画し、一つのメッシュ毎に３方向の熱の出入り（計６つの式）を解く。図３Ａに示すように、一辺のメッシュの数が例えば１０個である場合には、１０×１０×１０＝１０００個のメッシュについて計算が必要となる。一つのメッシュあたりの式が６つあるので、計算量は、１０００×６式＝６０００となる。
【００２４】
これに対し、本発明では、図３Ｂに示すように、図２に示す断面Ａ１，Ａ２，Ａ３の交差点となる内部空間の温度を算出対象とした場合に、この温度算出対象となるメッシュＭ’を通るＸ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群を用いる。具体的には、熱移動量算出部１２が、Ｘ軸上のメッシュ群に沿って温度算出対象メッシュＭ’に対して出入りする熱量総和、Ｙ軸上のメッシュ群に沿って温度算出対象メッシュＭ’に対して出入りする熱量総和、及びＺ軸上のメッシュ群に沿って温度算出対象メッシュＭ’に対して出入りする熱量総和を算出する。
【００２５】
温度算出部１３（図１参照）が、温度算出対象メッシュＭ’に対する熱量総和を用いて、当該メッシュＭ’の単位時間経過後の温度を算出する。具体的には、図５に示すように、各方向からメッシュＭ’に対する熱量総和Ｑ[W]は、次の式（１）で表現できる。
Ｑ＝Ａｘ１（ｑｈ＋ｑｒ）ｘ１＋Ａｘ２（ｑｈ＋ｑｒ）ｘ２＋Ａｙ１（ｑｈ＋ｑｒ）ｙ１＋Ａｙ２（ｑｈ＋ｑｒ）ｙ２＋Ａｚ１（ｑｈ＋ｑｒ）ｚ１＋Ａｚ２（ｑｈ＋ｑｒ）ｚ２    …（１）
ただし、総熱量はＱ[W]、表面積はＡ[ｍ２]、熱流束はｑ[Ｗ・ｍ−２]、境膜熱伝達による熱量はｑｈ、放射伝熱による熱量はｑｒである。
【００２６】
メッシュＭ’の単位時間経過後の温度Ｔｎ＋１は、次の式（２）で算出できる。
Ｔｎ＋１＝Ｔｎ＋Δｔ・（Ｑ／ρ・Ｃｐ・Ｖ）    …（２）
ただし、空間温度がＴ[K]、現時点の空間温度がＴｎ、空間密度がρ[kg・ｍ−３]、空間比熱がＣｐ[J・Kg−１・Ｋ−１]、空間体積がＶ[ｍ３]、単位時間がｔ[s]である。
【００２７】
上記熱量総和を算出するために必要となる伝熱式データ１４が、図１に示す記憶部１１に記憶されている。
【００２８】
熱量計算に用いる伝熱式データ１４は、３つのデータ（１４ｂ、１４ｃ、１４ｄ）を有する。温度差により同一物質中を移動する熱を表現する熱伝導式のデータ１４ｂ、温度差により異なる物質の境界を超えて移動する熱を表現する境膜熱伝達式のデータ１４ｃ、及び、温度差により物質表面間を移動する熱を表現する放射伝熱式のデータ１４ｄである。これらのデータ（１４ｂ、１４ｃ、１４ｄ）は、伝熱式データ１４として記憶部１１に記憶されている。これらのどの式（１４ｂ、１４ｃ、１４ｄ）を用いるかは、メッシュを占める物質に応じて定まる。そのため、図１に示す伝熱式データ１４は、図４に示すように各メッシュに、当該メッシュを占める物質を識別するためのメッシュに関するデータ１４ａ（図１参照）を有する。図４の上段部は、断面の概略図であり、図４の中段部は、断面とメッシュを対応付けた図であり、図４の下段部は、各メッシュの温度分布を示すイメージ図である。一つのメッシュは、一種類の物質に対応している。また、同図には、物質に応じて利用する伝熱式を図示している。熱伝導はｑｃであり、境膜熱伝達はｑｈであり、放射伝熱はｑｒを示す。
【００２９】
各軸上のメッシュ群に沿って温度算出対象メッシュＭ’に対して出入りする熱量総和を計算するために、各軸上のメッシュ群を構成するメッシュ一つ毎に、単位時間あたりの熱量の出入りを計算する。
【００３０】
外気と断熱材間の境界部分の熱計算式について図６Ａを用いて説明する。同図に示すように、メッシュをｉ０〜ｉｎで表現している。同図に示すように、外気メッシュ、境膜メッシュ、断熱材メッシュの温度をＴｉ−１、Ｔｉ、Ｔｉ＋１と表記して説明する。外気から境膜への入力は、境膜熱伝達ｑｈと放射伝熱ｑｒであり、境膜から断熱材への出力は、熱伝導ｑｃとなる。熱移動量は等しいので、ｑｃ＝ｑｈ＋ｑｒとなる。
熱伝導による熱移動量：
ｑｃ＝ｋ（∂Ｔ／∂ｘ）      …（３）
境膜熱伝達による熱移動量：
ｑｈ＝ｈ（Ｔｉ＋１−Ｔｉ）      …（４）
放射伝熱による熱移動量：
ｑｒ＝εｉ＋１σ（Ｔ４ｉ−１−Ｔ４ｉ＋１）    …（５）
ただし、メッシュ長さがｘ[m]、ステファンボルツマン定数がσ[W/m2K4]である。
すなわち、温度Ｔｉは、次の式により算出できる。
【数１】
2015187827.tif 000003
ただし、ｒ＝ｋ／ｈ・Δｘである。
【００３１】
断熱材内部又は蓄熱材内部の熱計算式について図６Ｂを用いて説明する。同図に示すように、３つの断熱材の時間ｍ時点の温度をＴｉ−１，ｍ、Ｔｉ，ｍ、Ｔｉ＋１，ｍと表記して説明する。所定の単位時間経過後の温度は、Ｔｉ−１，ｍ＋１、Ｔｉ，ｍ＋１、Ｔｉ＋１，ｍ＋１となる。勿論、蓄熱材と境膜間（境膜、蓄熱材及び境膜）についても同様である。いずれも熱伝導ｑｃとなる。温度Ｔｉ，ｍ＋１は、次の式により算出できる。
【数２】
2015187827.tif 000004
ただし、α＝ｋ／ρ・Ｃｐである。
【００３２】
断熱材と蓄熱材の間、又は、蓄熱材同士の間の境界部分の熱計算式について図６Ｃを用いて説明する。同図に示すように、断熱材、境膜、空間、境膜及び蓄熱材の時間ｍ時点の温度をＴｉ−１，ｍ、Ｔｉ，ｍ、Ｔｉ＋１，ｍ、Ｔｉ＋２，ｍ、Ｔｉ＋３，ｍと表記して説明する。勿論、蓄熱材同士の間（蓄熱材、境膜、空間、境膜及び蓄熱材）についても同様である。断熱材からの境膜への入力は熱伝導ｑｃとなり、境膜から空間への出力は、境膜熱伝達ｑｈと放射伝熱ｑｒである。
放射伝熱による熱移動量：
【数３】
2015187827.tif 000005
ただし、放射伝熱面間の形態係数はＦ、ζ＝１−εである。
すなわち、温度Ｔｉは、次の式により算出できる。
【数４】
2015187827.tif 000006
ただし、ｒ＝ｋ／（ｈ・Δｘ）である。
【００３３】
すなわち、上記式（３）〜（９）を用いることにより、熱移動量算出部１２は、温度算出対象メッシュＭ’を通るＸ軸上のメッシュ群それぞれについて、温度差に起因する熱移動量及び熱移動による温度変化量を計算することによって、温度算出対象メッシュＭ’に出入りする熱量総和を計算する。温度算出対象メッシュＭ’を通るＹ軸上のメッシュ群、温度算出対象メッシュＭ’を通るＺ軸上のメッシュ群についても同様である。
【００３４】
本実施形態では、蓄熱材の能力を適切に再現するために、図７及び図１に示すように、蓄熱材の比熱を温度毎に関連づけた温度依存比熱データ１７ａを記憶する比熱記憶部１７を有している（図１参照）。熱移動量算出部１２は、式（７）を用いた蓄熱材に対応するメッシュでの熱量計算には、比熱記憶部１７に記憶されている温度依存比熱データ１７ａのうち、現時点の蓄熱材２１の温度に対応する比熱Ｃｐを読み出して熱移動計算に用いている。勿論、断熱材の比熱は常に一定である。
【００３５】
同様に、解析パラメータとして、経過時間と外気温を経過時間毎に関連づけたデータを記憶しており、現時点に対応する外気温を読み出して、熱移動計算に用いている。勿論、外気温を変化させずに一定に設定してもよい。
【００３６】
熱移動量算出部１２及び温度算出部１３による計算は、所定の終了条件が成立するまで、単位時間毎に繰り返し実行され、各メッシュの温度及び温度算出対象メッシュＭ’の経時温度が算出される。
【００３７】
［温度算出方法］
上記装置１を用いた断熱容器内温度を算出する方法について、図８を参照しつつ説明する。
【００３８】
まず、操作部１５に対して設定操作がなされると、ステップＳＴ１において、設定部１０が、指定されたファイルを読み込み、解析パラメータ及び伝熱式データ１４を記憶部１１に設定する。ここでは、初期条件（外気温度、初期温度等）、放射伝熱設定（放射率）などの解析に必要なパラメータが設定される。
【００３９】
次のステップＳＴ２において、熱移動量算出部１２は、解析に必要なデータ（伝熱式データ１４及び解析パラメータ）を記憶部１１から読み出す。次のステップＳＴ３において、熱移動量算出部１２は、予め設定された解析条件と伝熱式データを用い、Ｘ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群に沿って温度算出対象となるメッシュＭ’に対して出入りする単位時間あたりの熱量総和を算出する。具体的には、図５に例示する外気メッシュからメッシュＭ’までの間にある各メッシュについて、熱移動量計算及び熱移動による単位時間経過後の新たな温度を算出する。
【００４０】
次のステップＳＴ４において、温度算出部１３は、熱移動量算出部１２が算出した熱量総和（図５参照）に基づき、温度算出対象メッシュＭ’の単位時間経過後の温度を算出する。
【００４１】
次のステップＳＴ５において、装置１は、所定の終了条件を満たすか否かを判定する。判定基準として、例えば断熱容器内部空間（温度算出対象メッシュＭ’）の温度が所定の範囲内を超えた場合、又は、計算ステップ数（解析時間でもよい）が所定値を超えた場合に、終了条件を満たすとすることが挙げられる。ステップＳＴ５において、所定の終了条件を満たしていないと判定された場合（ＳＴ５：ＮＯ）には、次のステップＳＴ６において、解析時点を所定の単位時間進め、ステップＳＴ３の実行に移行する。
【００４２】
一方、ステップＳＴ５において、所定の終了条件を満たしたと判定した場合（ＳＴ５：ＹＥＳ）には、次のステップＳＴ７において、演算結果を出力部１６（ディスプレイやファイル）を介して外部に出力し、処理を終了する。
【００４３】
以上のように、本実施形態における断熱容器内の温度算出装置は、蓄熱材２１を収容した断熱容器２０内温度を算出する装置１であって、
蓄熱材２１と断熱容器２０とを含む解析対象となる空間をＸ軸、Ｙ軸及びＺ軸に沿って複数のメッシュに区画した場合に、温度算出対象となるメッシュＭ’を通るＸ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群それぞれに、メッシュを占める物質に応じた伝熱式を設定した伝熱式データ１４を記憶する記憶部１１と、
予め設定された解析条件と伝熱式データ１４を用い、Ｘ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群に沿って温度算出対象となるメッシュＭ’に対して出入りする単位時間あたりの熱量総和Ｑを算出する熱移動量算出部１２と、
熱量総和Ｑを用いて、温度算出対象となるメッシュＭ’の単位時間経過後の温度を算出する温度算出部１３と、
を備える。
【００４４】
また、本実施形態における断熱容器内の温度算出方法は、蓄熱材２１を収容した断熱容器２０内温度をコンピュータが算出する方法であって、
蓄熱材２１と断熱容器２０とを含む解析対象となる空間をＸ軸、Ｙ軸及びＺ軸に沿って複数のメッシュに区画した場合に、温度算出対象となるメッシュＭ’を通るＸ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群それぞれに、メッシュを占める物質に応じた伝熱式を設定した伝熱式データ１４を記憶部１１から読み出すステップ（ＳＴ２）と、
予め設定された解析条件と伝熱式データ１４を用い、Ｘ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群に沿って温度算出対象となるメッシュＭ’に対して出入りする単位時間あたりの熱量総和Ｑを算出するステップ（ＳＴ３）と、
熱量総和Ｑを用いて、温度算出対象となるメッシュＭ’の単位時間経過後の温度を算出するステップ（ＳＴ４）と、
を含む。
【００４５】
このような装置構成又は方法によれば、熱量計算が、温度算出対象となるメッシュＭ’を通るＸ軸上のメッシュ群、Ｙ軸上のメッシュ群及びＺ軸上のメッシュ群に限られるので、解析空間にある全メッシュについて計算する場合に比べて、著しく計算コストを低減することが可能となる。それでいて、温度算出対象となるメッシュＭ’に影響を与える三次元方向の熱量を計算しているので、二次元平面だけの算出結果に比して高精度な結果を得ることが可能となる。
【００４６】
特に、本実施形態では、伝熱式には、熱伝導ｑｃ及び境膜熱伝達ｑｈに関する式を含んでいるので、適切に熱移動計算を実施できる。
【００４７】
さらに、本実施形態では、伝熱式には、放射伝熱ｑｒに関する式を更に含んでいるので、放射伝熱ｑｒを考慮しない計算方法に比べて精度を向上させることが可能となる。それでいて、放射伝熱ｑｒの計算式を導入しても、上記のように計算対象となるメッシュ数が適切に抑制されているので、計算コストが乗数オーダーで増大することもない。
【００４８】
さらに、本実施形態の装置では、蓄熱材２１の比熱を温度毎に関連付けた温度依存比熱データ１７ａを記憶する比熱記憶部１７を備え、蓄熱材２１に対応するメッシュでの熱量計算には、比熱記憶部１７に記憶されている温度依存比熱データ１７ａのうち現時点の蓄熱材の温度に対応する比熱を読みだして用いるように構成されている。
同様に、本実施形態の方法では、蓄熱材２１に対応するメッシュの熱量を計算する際に、蓄熱材２１の比熱を温度毎に関連付けた温度依存比熱データ１７ａのうち、現時点の蓄熱材２１の温度に対応する比熱を読みだして熱量計算に用いる。
【００４９】
このような装置構成又は方法によれば、蓄熱材の比熱の温度依存性を考慮した計算が可能になるので、従来ように外気温及び熱流束が一定で蓄熱材の比熱が一定としていたシミュレーションではできなかった、条件変動を伴う解析を精度良く行うことが可能となる。
【００５０】
本実施形態に係るコンピュータプログラムは、上記断熱容器内の温度算出方法を構成する各ステップをコンピュータに演算させるプログラムである。
これらプログラムを実行することによっても、上記方法の奏する作用効果を得ることが可能となる。言い換えると、装置１は、上記方法を使用しているとも言える。
【００５１】
以上、本発明の実施形態について図面に基づいて説明したが、具体的な構成は、これらの実施形態に限定されるものでないと考えられるべきである。本発明の範囲は、上記した実施形態の説明だけではなく特許請求の範囲によって示され、さらに特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれる。
【００５２】
例えば、本実施形態では、熱伝導、境膜伝達及び放射伝熱の３つの伝熱式を用いているが、放射伝熱式を省略することも可能である。これにより、精度が多少落ちるものの、計算コストを有効に低減することが可能となる。
【００５３】
上記は装置として説明しているが、方法及びコンピュータプログラムとしても特定可能である。
【００５４】
第１実施形態では、図２及び図９に例示するように、蓄熱材２１を除く収容用の内部空間Ａｒ１が温度算出対象としての１つのメッシュＭ’に対応するように、解析対象となる空間が複数のメッシュに区画されている。
【００５５】
このようにすれば、収容用の内部空間Ａｒ１に注目したメッシュを構成でき、収容用の内部空間Ａｒ１を複数のメッシュに区画する場合に比べて計算コストを低減させることが可能となる。勿論、計算コストの増大を或る程度許容できるのであれば、収容用の内部空間Ａｒ１に、温度算出対象となるメッシュを複数設定してもよい。
【００５６】
また、第１実施形態では、温度算出対象となるメッシュＭ’を通るＸ軸、Ｙ軸及びＺ軸について、軸に沿った物質配置パターンは、各軸において単一のパターンが設定されている。すなわち、メッシュＭ’を通るＸ軸上の物質配置パターンは、「断熱材、空間、蓄熱材、空間、蓄熱材、内部空間、蓄熱材、空間、蓄熱材、空間、断熱材」という１つのパターンのみが存在する。メッシュＭ’を通るＹ軸上の物質配置パターンは、「断熱材、内部空間、断熱材」という１つのパターンのみが存在する。メッシュＭ’のいずれの部位を通るＺ軸上の物質配置パターンは、「断熱材、内部空間、断熱材」という１つのパターンのみが存在する。
【００５７】
この構成によれば、収容用の内部空間Ａｒ１に対応する１つのメッシュＭ’に対し、Ｘ軸２つ、Ｙ軸２つ、Ｚ軸２つの計６軸に沿った伝熱式を用いればよいので、有限要素法又は有限体積法を用いた全メッシュ計算に比べて、計算コストを必要最小限に抑制することが可能となる。
【００５８】
［第２実施形態］
第１実施形態では、比較的小さな断熱容器で、収容用の内部空間の全面を一方（例えばＸ軸）から１つの蓄熱材で覆う例である。しかし、図１０Ａに示す第２実施形態のように、断熱容器が大きくなり、内部空間の全面を一方（例えばＸ軸）から１つの蓄熱材で覆うことができず、蓄熱材をＸ軸に沿って内部空間に投影した場合に蓄熱材が配置されていない領域が生じる場合や使用する蓄熱材の種類が異なる場合がある。例えば、図１０Ａの例では、ｘ１１軸には、蓄熱材２１が配置されているが、ｘ１２軸には蓄熱材が配置されていないため、ｘ１１軸を通ってメッシュＭ’に至る熱量と、ｘ１２軸を通ってメッシュＭ’に至る熱量が異なる。また、ｘ１１軸に配置される蓄熱材２１と、ｘ１３軸に配置される蓄熱材２１’は種類が異なり、比熱などの各種パラメータが異なるので、ｘ１１軸を通ってメッシュＭ’に至る熱量と、ｘ１３軸を通ってメッシュＭ’に至る熱量が異なる。
【００５９】
このように、温度算出対象となるメッシュＭ’を通るＸ軸、Ｙ軸及びＺ軸のうち、少なくともいずれかの軸に沿った物質配列パターンが複数存在する場合に、各々の物質配列パターンに対応する伝熱式データが記憶部１１に記憶されている。図１０Ａの例では、ｘ１１軸、ｘ１２軸、ｘ１３軸に沿った各物質配列パターンに対応する伝熱式データが記憶部１１に記憶されている。また、記憶部１１には、図１０Ｂに示すように、各々の物質配列パターンがメッシュＭ’に占める面積に関する面積データを記憶する。例えば、図１０Ｂにて蓄熱材２１が占める領域は右上がりで間隔の広い斜線で示す領域であり、その面積はＡｘ１１である。同図にて蓄熱材２１’が占める領域は左上がり斜線であり、その面積はＡｘ１３である。同図にて蓄熱材が配置されていない領域は十字状に配置された右上がりで間隔の狭い斜線であり、その面積がＡｘ１２である。この場合、メッシュＭ’の表面積Ａｘ１＝Ａｘ１１＋Ａｘ１２＋Ａｘ１３となる。
【００６０】
熱移動量算出部１２は、各々の物質配列パターン（ｘ１１軸、ｘ１２軸、ｘ１３軸）に対応する伝熱式データと、各々の物質配列パターンが温度算出対象となるメッシュＭ’に占める面積（Ａｘ１１、Ａｘ１２、Ａｘ１３）とに基づき、各々の物質配列パターン（ｘ１１軸、ｘ１２軸、ｘ１３軸）毎に温度対象となるメッシュＭ’に出入りする単位時間あたりの熱量「Ａｘ１１（ｑｈ＋ｑｒ）ｘ１１、Ａｘ１２（ｑｈ＋ｑｒ）ｘ１２、Ａｘ１３（ｑｈ＋ｑｒ）ｘ１３」を算出する。熱移動量算出部１２は、物質配列パターン毎に算出した算出結果を合計することで、温度算出対象となるメッシュＭ’に出入りする単位時間あたりの熱量総和を算出する。図１０の例では、次の式（１０）を用いる。
Ａｘ１（ｑｈ＋ｑｒ）ｘ１＝Ａｘ１１（ｑｈ＋ｑｒ）ｘ１１＋Ａｘ１２（ｑｈ＋ｑｒ）ｘ１２＋Ａｘ１３（ｑｈ＋ｑｒ）ｘ１３    …（１０）
【００６１】
上記構成の装置１の動作は次の通りである。図８に示すように、ステップＳＴ２において、熱移動量算出部１２（図１参照）は、温度算出対象となるメッシュＭ’を通るＸ軸、Ｙ軸及びＺ軸のうち、少なくともいずれかの軸に沿った物質配列パターンが複数存在する場合に、各々の物質配列パターンに対応する伝熱式データ１４を記憶部１１から読み出す。
【００６２】
次のステップＳＴ３において、熱移動量算出部１２は、各々の物質配列パターン（ｘ１１軸、ｘ１２軸、ｘ１３軸）に対応する伝熱式データ１４と、各々の物質配列パターンが温度算出対象となるメッシュＭ’に占める面積（Ａｘ１１、Ａｘ１２、Ａｘ１３）とに基づき、各々の物質配列パターン（ｘ１１軸、ｘ１２軸、ｘ１３軸）毎に温度対象となるメッシュＭ’に出入りする単位時間あたりの熱量「Ａｘ１１（ｑｈ＋ｑｒ）ｘ１１、Ａｘ１２（ｑｈ＋ｑｒ）ｘ１２、Ａｘ１３（ｑｈ＋ｑｒ）ｘ１３」を算出する。次に、熱移動量算出部１２は、各計算結果を合計することで温度対象となるメッシュに出入りする単位時間あたりの熱量総和を算出する。
【００６３】
この構成又は方法によれば、種々の蓄熱材の配置や複数種類の蓄熱材を用いた場合でも、これら物質の配列パターンを考慮した熱量計算が可能となるので、温度を精度良く予測することが可能となる。なお、図１０Ａ及び図１０Ｂに示す例は一例であり、これに限定されるものではない。
【００６４】
［第３実施形態］
上記第１及び第２実施形態では、収容用の内部空間に対応する温度算出対象となるメッシュＭ’の温度が一様であるとしているが、内部空間又は内部空間に収容された収容物が自然対流を許容する流動体（液体や気体）である場合には、温度が一様に分布しない。すると、算出結果と実際の温度とに誤差が生じてしまう。例えば、蓄熱材を保冷剤として場合には、内部空間の上方の温度が高く、下方の温度が低くなることが考えられる。
【００６５】
そこで、内部空間又は内部空間に収容された収容物の自然対流を考慮するために、図２及び図９に示すように、断熱容器２０内において蓄熱材２１を除く収容用の内部空間Ａｒ１が温度算出対象としての１つのメッシュＭ’に対応するように、解析対象となる空間が複数のメッシュに区画されていることを前提とする。そして、図１２Ａに示すように、温度算出対象となるメッシュＭ’を、高さ方向に積層する複数（本実施形態では３つ）の小領域（Ｆ１、Ｆ２、Ｆ３）に分割する。
【００６６】
各々の小領域（Ｆ１、Ｆ２、Ｆ３）に出入りする単位時間あたりの熱量（Ｑ１、Ｑ２、Ｑ３）は、メッシュＭ’に対してＸ軸、Ｙ軸及びＺ軸に沿って出入りする単位時間あたりの熱量と、小領域間の単位時間あたりの対流熱量（Ｑｎ１、Ｑｎ２）とに基づき算出可能である。鉛直方向となるＺ軸に沿ってメッシュＭ’に出入りする熱量は、鉛直方向の上方及び下方にある小領域にそれぞれ出入りする。例えば、図１０に示すように、ｚ１軸であれば、小領域Ｆ１に出入りし、ｚ２軸であれば、小領域Ｆ３に出入りする。水平方向となるＸ軸及びＹ軸に沿ってメッシュＭ’に出入りする熱量は、各々の小領域（Ｆ１、Ｆ２、Ｆ３）に均等分割されて入力する。例えば、図１０に示すように、ｘ１軸及びｘ２軸に沿った熱量は、小領域（Ｆ１、Ｆ２、Ｆ３）に三等分される。
【００６７】
図１２Ａの例では、各々の小領域（Ｆ１、Ｆ２、Ｆ３）に対する単位時間あたりの熱量総和（Ｑ１、Ｑ２、Ｑ３）[W]は、次の式（１１）〜（１３）で表現できる。
Ｑ１＝｛１／３｝｛Ａｘ１（ｑｈ＋ｑｒ）ｘ１＋Ａｘ２（ｑｈ＋ｑｒ）ｘ２＋Ａｙ１（ｑｈ＋ｑｒ）ｙ１＋Ａｙ２（ｑｈ＋ｑｒ）ｙ２｝＋Ａｚ１（ｑｈ＋ｑｒ）ｚ１−Ｑｎ１    …（１１）
Ｑ２＝｛１／３｝｛Ａｘ１（ｑｈ＋ｑｒ）ｘ１＋Ａｘ２（ｑｈ＋ｑｒ）ｘ２＋Ａｙ１（ｑｈ＋ｑｒ）ｙ１＋Ａｙ２（ｑｈ＋ｑｒ）ｙ２｝−Ｑｎ２＋Ｑｎ１    …（１２）
Ｑ３＝｛１／３｝｛Ａｘ１（ｑｈ＋ｑｒ）ｘ１＋Ａｘ２（ｑｈ＋ｑｒ）ｘ２＋Ａｙ１（ｑｈ＋ｑｒ）ｙ１＋Ａｙ２（ｑｈ＋ｑｒ）ｙ２｝＋Ａｚ２（ｑｈ＋ｑｒ）ｚ２＋Ｑｎ２    …（１３）
【００６８】
各々の小領域（Ｆ１、Ｆ２、Ｆ３）の単位時間経過後の温度は、上述した式（２）を用いて算出できる。
【００６９】
小領域間の単位時間あたりの対流熱量（Ｑｎ１、Ｑｎ２）[W]は、次の式（１４）〜（１８）で表現できる。
【数５】
2015187827.tif 000007
ただし、熱量はＱ[W]、熱伝達係数はｈ[Ｗ・ｍ−２・Ｋ−１]、表面積はＡ[ｍ２]、温度はＴ[K]、小領域の境界の底面からの距離はＤ[ｍ]、内部空間又は収納物の動粘度はｖ[ｍ２・ｓ−１]、内部空間又は収納物の体膨張係数はαｖ[Ｋ−１]、重力加速度はｇ[ｍ・ｓ−２]、内部空間又は収納物の比熱はＣｐ[Ｊ・Ｋｇ−１・Ｋ−１]、内部空間又は収納物の粘度はη[Ｐａ・ｓ−１]、内部空間又は収納物の熱伝導度はλ[Ｗ・ｍ−１・Ｋ−１]、内部空間又は収納物の温度はＴＦ[K]、周辺壁面の平均温度はＴＷ[K]である。ここで、周辺壁面とは、内部空間Ａｒ１に面した側面を意味し、周辺壁面の平均温度は、図１２Ｂ及び式（１９）で表される。
【数６】
2015187827.tif 000008
図１２Ｂに示す側面ｘｗ１の温度はＴｘｗ１であり、側面ｘｗ２の温度はＴｘｗ２であり、側面ｙｗ１の温度はＴｙｗ１であり、側面ｙｗ２の温度はＴｙｗ２である。
【００７０】
プラントル数Ｐｒは、式（１７）と予め設定されたパラメータにより算出可能である。グラスホフ数Ｇｒは、式（１８）と予め設定されたパラメータにより算出可能である。上記式（１４）〜（１８）で表現される、小領域間の単位時間あたりの対流熱量に関する対流熱式データ１４ｅは、図１１に示すように、伝熱式データ１４の一部として記憶部１１に記憶されている。
【００７１】
図１１に示すように、装置１は、小領域熱量算出部１８と小領域温度算出部１９とを有する。
【００７２】
小領域熱量算出部１８は、熱移動量算出部１２により算出される温度算出対象となるメッシュＭ’に出入りする単位時間あたりの熱量総和と対流熱式データ１４ｅとに基づき、式（１１）〜（１８）を用い、各々の小領域（Ｆ１、Ｆ２、Ｆ３）に出入りする単位時間あたりの熱量総和（Ｑ１、Ｑ２、Ｑ３）を算出する。
【００７３】
小領域温度算出部１９は、小領域熱量算出部１８により算出された熱量総和（Ｑ１、Ｑ２、Ｑ３）に基づき、式（２）を用いて単位時間経過後の小領域の温度を算出する。
【００７４】
上記構成の装置１の動作について説明する。第１実施形態の図８におけるステップＳＴ４の処理が、図１３に示す内部空間温度算出処理に置き換わる。すなわち、小領域熱量算出部１８は、ステップＳＴ１０１において、温度算出対象となるメッシュＭ’を高さ方向に積層する複数の小領域（Ｆ１、Ｆ２、Ｆ３）に分割した場合に、小領域間の単位時間あたりの対流熱量に関する対流熱式データを含む前記伝熱式データを記憶部から読み出す。
【００７５】
次のステップＳＴ１０２において、小領域熱量算出部１８は、温度算出対象となるメッシュに出入りする単位時間あたりの熱量総和（熱移動量算出部１２が算出）と対流熱式データ１４ｅとに基づき、各々の小領域（Ｆ１、Ｆ２、Ｆ３）に出入りする単位時間あたりの熱量総和（Ｑ１、Ｑ２、Ｑ３）を算出する。
【００７６】
次のステップＳＴ１０３において、小領域温度算出部１９は、小領域熱量算出部１８が算出した結果に基づき単位時間経過後の小領域の温度を算出する。
【００７７】
上記構成又は方法によれば、小領域（Ｆ１、Ｆ２、Ｆ３）の間の対流を考慮して小領域（Ｆ１、Ｆ２、Ｆ３）毎に温度が算出可能となり、実際の温度分布に精度を近づけることが可能となる。
【００７８】
第３実施形態では、温度算出対象となるメッシュＭ’を高さ方向に積層する３つの小領域に分割しているが、少なくとも２つ以上に分割すれば、その数は適宜変更可能である。
【００７９】
上記の各実施形態で採用している構造を他の任意の実施形態に採用することは可能である。各部の具体的な構成は、上述した実施形態のみに限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々変形が可能である。
【符号の説明】
【００８０】
１…装置
１１…記憶部
１２…熱移動量算出部
１３…温度算出部
１４…伝熱式データ
１４ｅ…対流熱式データ
１７…比熱記憶部
１８…小領域熱量算出部
１９…小領域温度算出部
２０…断熱容器
２１…蓄熱材
Ａｒ１…収容用の内部空間
【図１】
2015187827.tif 000009
【図２】
2015187827.tif 000010
【図３Ａ】
2015187827.tif 000011
【図３Ｂ】
2015187827.tif 000012
【図４】
2015187827.tif 000013
【図５】
2015187827.tif 000014
【図６Ａ】
2015187827.tif 000015
【図６Ｂ】
2015187827.tif 000016
【図６Ｃ】
2015187827.tif 000017
【図７】
2015187827.tif 000018
【図８】
2015187827.tif 000019
【図９】
2015187827.tif 000020
【図１０Ａ】
2015187827.tif 000021
【図１０Ｂ】
2015187827.tif 000022
【図１１】
2015187827.tif 000023
【図１２Ａ】
2015187827.tif 000024
【図１２Ｂ】
2015187827.tif 000025
【図１３】
2015187827.tif 000026