(57)【要約】【課題】対象物の殺菌の観点のみならず殺菌以外の観点でも加熱条件を評価し得るシミュレーション方法を提供する。【解決手段】シミュレーション方法は、対象物の加熱条件が設定されるステップＳ０２と、加熱条件に基づき、対象物における中心部の温度としての第１推定温度Ｔ１ｎを演算するステップＳ０３と、第１推定温度Ｔ１ｎに基づき、対象物における中心部の第１領域を囲む第２領域内の温度としての第２推定温度Ｔ２ｎを演算するステップＳ０４とを含む。【選択図】図１
2017215191.tif 000002
【特許請求の範囲】
【請求項１】
対象物の加熱条件が設定されるステップと、
前記加熱条件に基づき、前記対象物における中心部の温度としての第１推定温度を演算するステップと、
前記第１推定温度に基づき、前記対象物における前記中心部の第１領域を囲む第２領域内の温度としての第２推定温度を演算するステップとを含む
ことを特徴とするシミュレーション方法。
【請求項２】
前記第１推定温度、前記第２推定温度、前記対象物の体積に対する前記第１領域の体積の割合、及び前記対象物の体積に対する前記第２領域の体積の割合に基づき、前記対象物の温度履歴を評価するステップをさらに含む
請求項１に記載のシミュレーション方法。
【請求項３】
前記対象物の温度履歴を評価するステップでは、前記第１推定温度、前記第２推定温度、前記対象物の体積に対する前記第１領域の体積の割合、及び前記対象物の体積に対する前記第２領域の体積の割合に基づき第３推定温度を演算し、該第３推定温度により前記対象物の温度履歴を評価する
請求項２に記載のシミュレーション方法。
【請求項４】
前記加熱条件は、雰囲気温度と前記対象物の加熱時間とを含み、
前記第２推定温度は、前記対象物における前記第１推定温度と前記雰囲気温度との平均値である
請求項１乃至３のいずれか１項に記載のシミュレーション方法。
【請求項５】
演算装置に、加熱対象物の推定温度の演算を実行させるためのシミュレーションプログラムであって、
演算装置に、
対象物の加熱条件の設定を受け付けるステップと、
前記加熱条件に基づき、前記対象物における中心部の温度としての第１推定温度を演算するステップと、
前記第１推定温度に基づき、前記対象物における前記中心部の第１領域を囲む第２領域内の温度としての第２推定温度を演算するステップとを実行させる
シミュレーションプログラム。
【請求項６】
請求項５に記載のシミュレーションプログラムを内蔵した記憶媒体を含み、前記プログラムを演算装置が実行させるシミュレーション装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、加熱中の処理対象物の温度変化を計算するシミュレーション方法、シミュレーションプログラム、及びこのプログラムを内蔵した記憶媒体を含むシミュレーション装置に関する。
【背景技術】
【０００２】
従来、缶詰やレトルト食品といった包装食品を製造する際に、製造された包装食品の加熱殺菌が行われている。特定の加熱条件が食品の殺菌に適しているか否かは、一般に、温度と時間との関係で表される殺菌値であるＦ値により評価されている。特定の加熱条件下における食品の温度履歴が定められたＦ値を満たす場合には、加熱条件が食品の殺菌に適していると評価できる。例えば、レトルト食品のＦ値は、食品衛生法により、１２０．０℃４分相当以上とされている。
【０００３】
食品の温度履歴は、加熱中の食品の温度をセンサー等により実測することでも確認できるが、形状や大きさの異なる食品に対する温度履歴を実測により確認するためには、コスト及び時間が必要である。このようなコスト及び時間を低減するために、コンピュータを用いて加熱中の食品の推定温度を演算するシミュレーション方法が普及している（例えば、特許文献１参照）。
【０００４】
雰囲気からの伝熱により食品の温度が上昇することを前提とすると、食品全体のうちの中心部が、最も熱が伝わりにくく殺菌されにくい。従って、従来のシミュレーション方法を用いて加熱条件を評価する場合、食品の中心部の推定温度を指標として、加熱条件が評価されていた。なお、ここでいう中心部とは、物体の物理的中心ではなく、温度の上昇又は下降の最も遅れる部分を指す。
【０００５】
さらに、食品の加熱は、殺菌だけでなく食品に含まれるたんぱく質やビタミンの分解等にも影響する。また、食品以外の加熱についても同様であり、例えば、医薬品や医療機器の加熱は、医薬品に含まれる成分や医療機器を構成する材料等に悪影響を及ぼすおそれがある。そのため、処理対象物の加熱条件を、殺菌以外の観点でも評価することが求められている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第３０７１４１２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明は、斯かる実情に鑑み、対象物の殺菌の観点のみならず殺菌以外の観点でも加熱条件を評価し得るシミュレーション方法を提供するものである。
【課題を解決するための手段】
【０００８】
本発明にかかるシミュレーション方法は、対象物の加熱条件が設定されるステップと、前記加熱条件に基づき、前記対象物における中心部の温度としての第１推定温度を演算するステップと、前記第１推定温度に基づき、前記対象物における前記中心部の第１領域を囲む第２領域内の温度としての第２推定温度を演算するステップとを含むことを特徴とする。
【０００９】
上記シミュレーション方法により演算される第１、第２推定温度は、殺菌の観点や殺菌以外の観点で加熱条件を評価する際に用いられる。第１推定温度は、対象物の中心部の推定温度であるため、殺菌の観点で加熱条件を評価する際の指標に適している。一方、第２推定温度は、対象物における中心部の領域を囲む領域内の温度であるため、殺菌以外の観点で加熱条件を評価する際の指標として好ましい。このように、上記シミュレーション方法により、対象物の殺菌のみならず殺菌以外の観点でも加熱条件をより正確に評価することができる。
【００１０】
本発明にかかるシミュレーション方法の一態様として、前記第１推定温度、前記第２推定温度、前記対象物の体積に対する前記第１領域の体積の割合、及び前記対象物の体積に対する前記第２領域の体積の割合に基づき、前記対象物の温度履歴を評価するステップをさらに含んでもよい。
【００１１】
上記シミュレーション方法において、対象物の体積に対する第１、第２領域の体積の割合を考慮することで、より正確に（対象物の実物の温度に近似した値で）温度を求めることができ、加熱条件をより正確に評価することができる。
【００１２】
本発明にかかるシミュレーション方法の一態様として、前記対象物の温度履歴を評価するステップでは、前記第１推定温度、前記第２推定温度、前記対象物の体積に対する前記第１領域の体積の割合、及び前記対象物の体積に対する前記第２領域の体積の割合に基づき第３推定温度を演算し、該第３推定温度により前記対象物の温度履歴を評価してもよい。
【００１３】
上記シミュレーション方法において、第３推定温度という一つの温度を用いることで、第１、第２推定温度という複数の温度を用いる場合よりも、加熱条件をより簡単に評価することができる。
【００１４】
本発明にかかるシミュレーション方法の一態様として、前記加熱条件は、雰囲気温度と前記対象物の加熱時間とを含み、前記第２推定温度は、前記対象物における前記第１推定温度と前記雰囲気温度との平均値であってもよい。
【００１５】
上記シミュレーション方法において、第２推定温度は、第１推定温度と雰囲気温度との平均値であるため、熱収支等の複雑な関係式を用いて演算される場合よりも、より簡単に演算される。
【００１６】
本発明にかかるシミュレーションプログラムは、演算装置に、加熱対象物の推定温度の演算を実行させるためのシミュレーションプログラムであって、演算装置に、対象物の加熱条件の設定を受け付けるステップと、前記加熱条件に基づき、前記対象物における中心部の温度としての第１推定温度を演算するステップと、前記第１推定温度に基づき、前記対象物における前記中心部の第１領域を囲む第２領域内の温度としての第２推定温度を演算するステップとを実行させることを特徴とする。
【００１７】
第１推定温度は、対象物の中心部の推定温度であるため、殺菌の観点で加熱条件を評価する際の指標に適しており、第２推定温度は、対象物における中心部の領域を囲む領域の温度であるため、殺菌以外の観点で加熱条件を評価する際の指標として好ましい。このように、上記シミュレーションプログラムにより、殺菌のみならず殺菌以外の観点でも加熱条件を評価することができる。
【００１８】
本発明にかかるシミュレーション装置は、上記シミュレーションプログラムを内蔵した記憶媒体を含み、前記シミュレーションプログラムを演算装置が実行させることを特徴とする。
【００１９】
第１推定温度は、対象物の中心部の推定温度であるため、殺菌の観点で加熱条件を評価する際の指標に適しており、第２推定温度は、対象物における中心部の領域を囲む領域の温度であるため、殺菌以外の観点で加熱条件を評価する際の指標として好ましい。このように、上記シミュレーション装置により、殺菌のみならず殺菌以外の観点でも加熱条件を評価することができる
【発明の効果】
【００２０】
本発明によれば、対象物の殺菌の観点のみならず殺菌以外の観点でも加熱条件を評価し得るシミュレーション方法を提供することができる。
【図面の簡単な説明】
【００２１】
【図１】図１は、本発明の一実施形態にかかるシミュレーション方法のフローチャート図である。
【図２】図２は、本発明の一実施形態にかかるシミュレーション方法により推定温度を演算する対象物のモデル図である。
【図３】図３は、本発明の別の実施形態にかかるシミュレーション方法により推定温度を演算する対象物のモデル図である。
【発明を実施するための形態】
【００２２】
以下、本発明のシミュレーション方法について、添付図面を参酌して説明する。本実施形態のシミュレーション方法は、図１のフローチャート図に示すように、対象物の各種条件又は物性がシミュレーション装置に設定されるステップ（Ｓ０１）と、加熱条件がシミュレーション装置に設定されるステップ（Ｓ０２）と、シミュレーション装置が第１推定温度Ｔ１ｎを演算するステップ（Ｓ０３）と、シミュレーション装置が第２推定温度Ｔ２ｎを演算するステップ（Ｓ０４）と、シミュレーション装置が第３推定温度Ｔ３ｎを演算するステップ（Ｓ０５）とを含む。本実施形態のシミュレーション方法は、シミュレーション方法を実施できるプログラムを内蔵した記憶媒体と、プログラムを実行させる演算装置（ＣＰＵ）とを含むシミュレーション装置で実施される。シミュレーションの結果は、例えば、シミュレーション装置に設けられたディスプレイに表示されるとともに、メモリーカードに保存される。以下、本実施形態のシミュレーション方法に含まれる各ステップについて、順に説明する。なお、本実施形態では、対象物の中心部の温度を演算する方法として、ＡＴＳ法（Ａｍｂｉｅｎｔ  Ｔｅｍｐｅｒａｔｕｒｅ  Ｓｌｉｄｅ  ｍｅｔｈｏｄ）を採用した場合について説明する。
【００２３】
そこで、まず、ＡＴＳ法の概略について説明する。ＡＴＳ法では、雰囲気からの伝熱により対象物の温度が上昇することを前提とする。また、「雰囲気が対象物に与える熱量」と「対象物が雰囲気から受け取る熱量」とが一致することを前提として計算が行われる。ここで、単位時間をΔｔとし、対象物の中心点を中心点Ｐとし、雰囲気温度及び対象物の中心点温度をそれぞれＴｗｎ、Ｔｐｎとし（下添え字ｎは時間刻みのｎ番目を示す）、対象物の表面温度が雰囲気温度Ｔｗｎと一致することとし、対象物の内部の温度勾配を直線とし、対象物の表面から中心点までの距離をＬとし、対象物の熱伝導率をｋとし、対象物の表面積をＡとすると、単位時間当たりに「雰囲気が対象物に与える熱量」即ち下記式の左辺が得られる。
【００２４】
また、対象物の体積をＶとし、対象物の密度をρとし、対象物の比熱をｃｐとし、中心点温度Ｔｐｎは体積平均温度Ｔｐ＊ｎと一致するものとすると、単位時間当たりに「対象物が雰囲気から受け取る熱量」即ち下記式の右辺が得られる。
ｋＡ（Ｔｗｎ−１−Ｔｐｎ−１）Δｔ/Ｌ＝Ｖρｃｐ（Ｔｐ＊ｎ−Ｔｐ＊ｎ−１)
【００２５】
上記式においてｋ/（ρｃｐ)を熱拡散係数αとして整理すると、下記式が得られる。
Ｔｐ＊ｎ＝Ｔｐ＊ｎ−１＋αΔｔ（Ｔｗｎ−１−Ｔｐｎ−１）／Ｌ２
【００２６】
実際には、中心点温度Ｔｐｎは体積平均温度Ｔｐ＊ｎと異なるため、中心点温度Ｔｐｎと体積平均温度Ｔｐ＊ｎとのずれ率βを用いて、Ｔｐｎ＝Ｔｐ＊ｎβとすると、下記式が得られる。
Ｔｐｎ＝Ｔｐｎ−１＋αβΔｔ（Ｔｗｎ−１−Ｔｐｎ−１）／Ｌ２
【００２７】
さらに、αβΔｔ／Ｌ２を伝熱係数τとすると、下記式が得られる。下記式では、単位時間あたりに、対象物の表面温度Ｔｗｎ−１と中心点温度Ｔｐｎ−１との差分に伝熱係数τを乗じた分だけ、中心点温度が上昇することが表されている。
Ｔｐｎ＝Ｔｐｎ−１＋τ（Ｔｗｎ−１−Ｔｐｎ−１）
【００２８】
ところが、実際の対象物では表面と中心点とが離れているため、中心点の温度変化が表面の温度変化よりも遅れるが、この遅れは上記式に反映されていない。これに対して、上記式において、表面温度として、雰囲気温度Ｔｗｎ−１の代わりに、雰囲気温度Ｔｗｎ−１の時間変化をδ１だけ遅らせた仮想的な雰囲気温度Ｔｗ†ｎ−１を用いると、下記式が得られる。
Ｔｐｎ＝Ｔｐｎ−１＋τ（Ｔｗ†ｎ−１−Ｔｐｎ−１）
【００２９】
従って、ＡＴＳ法を採用する場合には、ステップＳ０１において、例えば、ユーザーの入力により、伝熱係数τ及び遅れ時間δ１がシミュレーション装置に設定される。伝熱係数τ及び遅れ時間δ１は、対象物を構成する材料や対象物の形状により異なる。ユーザーは、対象物の伝熱係数τ及び遅れ時間δ１のデータを所持している場合、このデータをシミュレーション装置に入力する。ユーザーは、対象物の伝熱係数τ及び遅れ時間δ１のデータを所持していない場合、伝熱係数τ及び遅れ時間δ１を求め、これをシミュレーション装置に入力する。伝熱係数τ及び遅れ時間δ１を求めるには、ユーザーが、適当な伝熱係数τ及び遅れ時間δ１を予測し、この予測値に基づきＡＴＳ法により第１推定温度Ｔ１ｎを演算し、この第１推定温度Ｔ１ｎと実測温度とを比較して、伝熱係数τ及び遅れ時間δ１の予測値を調整し、調整した伝熱係数τ及び遅れ時間δ１に基づき第１推定温度Ｔ１ｎを再度演算し、実測温度に近づくまでこれを繰り返すことにより適切な伝熱係数τ及び遅れ時間δ１を求めることができる。即ち、トライアルアンドエラーの手法を採用し得る。なお、対象物の実測温度は、例えば、サーミスタのような温度検出センサーを用いて測定される。
【００３０】
例えば、伝熱係数τが不明である場合、予測した伝熱係数τ’１に基づき第１推定温度Ｔ’１ｎを演算し、第１推定温度Ｔ’１ｎと実測温度との差を計算する。第１推定温度Ｔ’１ｎと実測温度との差が、所望の値以下である場合、伝熱係数τ’１をシミュレーション装置に入力する。
【００３１】
第１推定温度Ｔ’１ｎと実測温度との差が所望の値よりも大きい場合、伝熱係数τ’１に近い四つの伝熱係数τ’’１…を予測し、四つの伝熱係数τ’’１…各々に基づき第１推定温度Ｔ’’１ｎ…を演算し、四つの第１推定温度Ｔ’’１ｎ…と実測温度との差を計算する。
【００３２】
第１推定温度Ｔ’’１ｎ…と実測温度との差の最小値が、第１推定温度Ｔ’１ｎと実測温度との差よりも小さく、且つ、所望の値よりも小さい場合、第１推定温度Ｔ’’１ｎ…と実測温度との差の最小値の演算に用いた伝熱係数τ’’１をシミュレーション装置に入力する。
【００３３】
第１推定温度Ｔ’’１ｎ…と実測温度との差の最小値が、第１推定温度Ｔ’１ｎと実測温度との差よりも小さく、且つ、所望の値よりも大きい場合、第１推定温度Ｔ’’１ｎ…と実測温度との差の最小値の演算に用いた伝熱係数τ’’１に近い四つの伝熱係数τ’’’１…を予測し、四つの伝熱係数τ’’’１…各々に基づき第１推定温度Ｔ’’’１ｎ…を演算する等の一連の計算を繰り返す。
【００３４】
第１推定温度Ｔ’’１ｎ…と実測温度との差の最小値が、第１推定温度Ｔ’１ｎと実測温度との差よりも大きい場合、伝熱係数τ’１にさらに近似した四つの伝熱係数τ’’’１…を予測し、四つの伝熱係数τ’’’１…各々に基づき第１推定温度Ｔ’’’１ｎ…を演算する等の一連の計算を繰り返す、又は、伝熱係数τ’１をシミュレーション装置に入力する。
【００３５】
ステップＳ０２では、例えば、ユーザーの入力により、シミュレーション装置に対象物の加熱条件が設定される。加熱条件は、例えば、雰囲気温度Ｔｗｎ及び加熱時間である。
【００３６】
ステップＳ０３では、シミュレーション装置が、ステップＳ０１及びステップＳ０２において設定された各種条件に基づき、対象物における中心部の温度である第１推定温度Ｔ１ｎを演算する。ここでいう「対象物における中心部」として、対象物の中心点Ｐ又はその近傍を採用し得るが、本実施形態では、対象物の中心点Ｐを採用する。この場合、第１推定温度Ｔ１として、対象物の中心点温度Ｔｐｎが用いられる。なお、対象物の中心点Ｐとして、対象物の質量中心、即ち、対象物の重心を採用し得る。対象物の重心の位置は、公知の方法により求めることができる。
【００３７】
中心点温度Ｔｐｎは、上述したように、ＡＴＳ法で導出される下記式により演算される。
Ｔｐｎ＝Ｔｐｎ−１＋τ（Ｔｗ†ｎ−１−Ｔｐｎ−１）
【００３８】
ステップＳ０４では、シミュレーション装置が、第１推定温度Ｔ１ｎ及び雰囲気温度Ｔｗｎに基づき、前記第２領域内の温度としての第２推定温度Ｔ２ｎを演算する。本実施形態では、第２推定温度Ｔ２ｎは、中心点温度Ｔｐｎ（第１推定温度Ｔ１ｎに相当）及び雰囲気温度Ｔｗｎの単純平均値として表され、具体的には、下記式により演算される。
Ｔ２ｎ＝（Ｔｐｎ＋ＴＷｎ）／２
【００３９】
ステップＳ０５では、シミュレーション装置が、第１推定温度Ｔ１ｎ、第２推定温度Ｔ２ｎ、対象物１００の体積に対する第１領域ＡＲ１の体積ＶＡＲ１の割合、及び対象物１００の体積に対する第２領域ＡＲ２の体積ＶＡＲ２の割合に基づいて、第３推定温度Ｔ３ｎを演算する。具体的には、第３推定温度Ｔ３ｎは、第１領域ＡＲ１と第２領域ＡＲ２との体積平均温度Ｔｐ＊ｎとして演算される。第１推定温度Ｔ１ｎとしては、上述したように、中心点温度Ｔｐｎが用いられる。第２推定温度Ｔ２ｎとしては、上述したように、中心点温度Ｔｐｎ及び雰囲気温度Ｔｗｎの単純平均値が用いられる。これにより、体積平均温度Ｔｐ＊ｎは、下記式により、中心点温度Ｔｐｎ及び雰囲気温度Ｔｗｎに基づき演算される。
Ｔｐ＊ｎ＝０．２９６Ｔｐｎ＋０．７０４（Ｔｗｎ＋Ｔｐｎ）／２
【００４０】
対象物の形状としては、立方体形状、直方体形状、その他の形状が考えられるが、本実施形態では、対象物の形状が、図２に示すような立方体形状であると仮定する。対象物１００は、対象物１００の中心部の領域である第１領域ＡＲ１と、第１領域ＡＲ１を囲む第２領域ＡＲ２とを含む。第１領域ＡＲ１が中心部の領域であれば、対象物をどのように二分割してもよいが、本実施形態では、対象物１００を中心部分と表層部分とに二分割し、中心部分を第１領域ＡＲ１とし、表層部分を第２領域ＡＲ２としている。
【００４１】
上記式の導出は、以下のように行う。対象物１００における一辺の長さは、図２に示すように、何れもＬ１とする。これにより、対象物１００の体積ＶはＬ１３と表され、対象物１００の表面積Ａは６Ｌ１２と表される。ここで、第１領域ＡＲ１の表面から第２領域ＡＲ２の表面までの距離をＬ２とすると、距離Ｌ２は対象物１００における表層厚み（Ｖ／Ａ）と一致するため、距離Ｌ２は下記のように演算される。
Ｌ２＝Ｖ／Ａ＝Ｌ１／６
【００４２】
この場合、第１領域ＡＲ１の体積ＶＡＲ１は、下記のように演算される。
ＶＡＲ１＝（Ｌ１−２・Ｌ２）３＝（２／３）３Ｌ１３
【００４３】
一方、体積平均温度Ｔｐ＊ｎは、中心点温度Ｔｐｎに、対象物１００の体積Ｖに対する第１領域ＡＲ１の体積ＶＡＲ１の比率を乗じたものと、中心点温度Ｔｐｎ及び雰囲気温度Ｔｗｎの単純平均値に、対象物１００の体積Ｖに対する第２領域ＡＲ２の体積ＶＡＲ２の比率を乗じたものとを足し合わせることで得られる。具体的には、体積平均温度Ｔｐ＊ｎは、下記式により演算される。
Ｔｐ＊ｎ＝Ｔｐｎ・ＶＡＲ１／Ｖ＋（（Ｔｗｎ＋Ｔｐｎ）／２）・（Ｖ−ＶＡＲ１）／Ｖ
上記式に、ＶＡＲ１＝（２／３）３Ｌ１３、及び、Ｖ＝Ｌ１３を代入することで、体積平均温度Ｔｐ＊ｎを演算する下記式が得られる。
Ｔｐ＊ｎ＝（２／３）３Ｔｐｎ＋（１−（２／３）３）・（Ｔｗｎ＋Ｔｐｎ）／２
上記式における数値を、小数点第４位で四捨五入すると下記式が得られる。
Ｔｐ＊ｎ＝０．２９６Ｔｐｎ＋０．７０４（Ｔｗｎ＋Ｔｐｎ）／２
【００４４】
このように、本実施形態に係るシミュレーション方法を実施することで、第１、第２、第３推定温度を演算することができる。演算された第１、第２、第３推定温度は、例えば、対象物がレトルト食品である場合、この食品の加熱条件を、殺菌の観点等の複数の観点で評価する場合の指標として用いられる。以下、殺菌の観点及びたんぱく質の分解の観点という２つの観点で、第１、第２、第３推定温度を指標として、食品の加熱条件を評価する場合について説明する。
【００４５】
特定の加熱条件が食品の殺菌の観点で適切か否かは、一般に、加熱時の食品の温度履歴が、殺菌評価積算値であるＦ値を満たすか否かで評価される。また、殺菌の観点では、最も殺菌されにくい食品の中心部の温度を評価することが好ましい。そのため、殺菌の観点で加熱条件を評価するには、第１推定温度Ｔ１ｎから第１推定温度履歴Ｔ１を求めて、この第１推定温度履歴Ｔ１が定められたＦ値に相当するか否かを評価する。
【００４６】
特定の加熱条件がたんぱく質の分解の観点で適切か否かは、例えば、加熱時の食品の温度履歴が、温度と時間との関係で表されるたんぱく質分解評価積算値であるＣ値を満たすか否かで評価される。また、たんぱく質の分解の観点では、加熱中の食品全体の温度を評価することが好ましい。そのため、たんぱく質の分解の観点で加熱条件を評価するには、第３推定温度Ｔ３ｎから第３推定温度履歴Ｔ３を求めて、この第３推定温度履歴Ｔ３が定められたＣ値を満たすか否かを評価する。
【００４７】
このように、殺菌の観点で加熱条件を評価する際には、第１推定温度Ｔ１ｎが指標として用いられ、たんぱく質の分解の観点で加熱条件を評価する際には、第３推定温度Ｔ３ｎが指標として用いられる。以下、本実施形態の効果をまとめて説明する。
【００４８】
本実施形態のシミュレーション方法により演算される第１推定温度Ｔ１ｎは、上述のように、殺菌の観点で加熱条件を評価する際の指標に適している。一方、第２推定温度Ｔ２ｎは、対象物における中心部の領域を囲む領域内の温度であるため、殺菌以外の観点で加熱条件を評価する際の指標として好ましい。このように、本実施形態のシミュレーション方法により、対象物の殺菌のみならず殺菌以外の観点でもより正確に加熱条件を評価することができる。
【００４９】
本実施形態のシミュレーション方法において、対象物１００の体積Ｖに対する第１、第２領域ＡＲ１、ＡＲ２の体積ＶＡＲ１、ＶＡＲ２の割合を考慮することで、より正確に（対象物１００の温度に近似した値で）温度を求めることができ、加熱条件をより正確に評価することができる。
【００５０】
第３推定温度Ｔ３ｎという一つの温度を用いることで、第１、第２推定温度Ｔ１ｎ、Ｔ２ｎという複数の温度を用いる場合よりも、加熱条件をより簡単に評価することができる。
【００５１】
第２推定温度Ｔ２ｎは、第１推定温度Ｔ１ｎと雰囲気温度Ｔｗｎとの平均値として算出されることにより、熱収支等の複雑な関係式等を用いて演算される場合よりも、より簡単に演算される。
【００５２】
なお、本発明のシミュレーション方法は、上記実施形態の方法に限定されるものではなく、本発明の要旨を逸脱しない範囲で種々の変更が可能である。
【００５３】
上記実施形態では、対象物の形状が立方体形状であると仮定したが、これに限定されない。例えば、対象物の形状が直方体形状、円柱形状や球形状等その他の形状である仮定してもよい。図３に示すように、対象物２００の形状が直方体形状であると仮定すると、以下のように、中心点温度Ｔｐｎ及び雰囲気温度Ｔｗｎに基づく体積平均温度Ｔｐ＊ｎ、即ち、第３推定温度Ｔ３ｎを演算することができる。
【００５４】
対象物２００が、中心部の領域である第１領域ＡＲ１と、第１領域ＡＲ１を囲む第２領域ＡＲ２とを含むものとする。第２領域ＡＲ２は、対象物２００における表層部分とする。対象物２００において、底面の短辺の長さをｗａとし、底面の長辺の長さをｗｂとし、高さをＨとする。この場合、対象物２００の体積ＶはＶ＝ｗａ・ｗｂ・Ｈで表される。また、対象物２００の表面積Ａは、２（ｗａ・ｗｂ＋ｗａ・Ｈ＋ｗｂ・Ｈ）で表される。一方、第１領域ＡＲ１の表面から第２領域ＡＲ２の表面までの距離Ｌ２とすると、第１領域ＡＲ１の体積ＶＡＲ１は、下記式で表される。
ＶＡＲ１＝（ｗａ−２Ｌ２）・（ｗｂ−２Ｌ２）・（Ｈ−２Ｌ２）
【００５５】
距離Ｌ２としては表層厚み（Ｖ／Ａ）を採用できるため、距離Ｌ２は下記式により演算される。
Ｌ２＝Ｖ／Ａ＝ｗａ・ｗｂ・Ｈ／（ｗａ・ｗｂ＋ｗａ・Ｈ＋ｗｂ・Ｈ）
【００５６】
さらに、体積平均温度Ｔｐ＊ｎは下記式で表される。
Ｔｐ＊ｎ＝Ｔｐｎ・ＶＡＲ１／Ｖ＋（（Ｔｗｎ＋Ｔｐｎ）／２）・（Ｖ−ＶＡＲ１）／Ｖ
【００５７】
上記式に、体積Ｖ及び体積ＶＡＲ１の数値を代入することで、中心点温度Ｔｐｎ及び雰囲気温度Ｔｗｎに基づき、体積平均温度Ｔｐ＊ｎが演算される。
【００５８】
このように、体積平均温度Ｔｐ＊ｎ、即ち、第３推定温度Ｔ３ｎを求めることで、対象物２００の形状により適した推定温度が演算される。
【００５９】
また、上記実施形態では、第２推定温度Ｔ２ｎとして、中心点温度Ｔｐｎと雰囲気温度Ｔｗｎとの単純平均値を用いていたが、これに限定されない。例えば、第２推定温度Ｔ２ｎとして、別の平均値を用いてもよい。具体的には、第２推定温度Ｔ２ｎとして、第２領域内の温度勾配を考慮した平均値、第２領域内の形状や比熱を考慮した平均値、第２領域を任意に二分割しこれらの体積を考慮した平均値等を用いてもよい。第２推定温度Ｔ２ｎとして、このような平均値を用いることで、対象物の実物の温度により近似した温度を得ることができる。
【００６０】
上記実施形態では、第３推定温度Ｔ３ｎを演算し、この第３推定温度Ｔ３ｎから求めた第３推定温度履歴Ｔ３を指標として、殺菌以外の観点で加熱条件を評価していたが、これに限定されない。例えば、第３推定温度Ｔ３ｎを演算せず、第１推定温度Ｔ１ｎ及び第２推定温度Ｔ２ｎを指標として、殺菌以外の観点で加熱条件を評価してもよい。第１推定温度Ｔ１ｎ及び第２推定温度Ｔ２ｎを指標とすれば、第３推定温度を指標する場合よりも、第１、第２領域の温度履歴をそれぞれ個別に評価できるため、対象物の実物の温度履歴に近似したより正確な加熱条件の評価が可能となる。
【００６１】
上記実施形態では、対象物を二分割して第３推定温度を演算していたが、これに限定されない。例えば、対象物を三以上の複数の領域に分割し、それぞれの領域における推定温度を演算してもよい。また、各推定温度から求めた推定温度履歴を指標として加熱条件を評価してもよい。三以上の複数の領域の各々における推定温度履歴を指標とすることで、対象物の温度履歴により適した加熱条件の評価が可能となる。
【００６２】
上記実施形態では、第１推定温度をＡＴＳ法により演算したが、これに限定されない。例えば、Ｂａｌｌの数式法等の別の方法によって、第１推定温度を演算してもよい。
【００６３】
上記実施形態のシミュレーション装置は、対象物を加熱する装置と別の装置であってもよいし、対象物を加熱する装置に組み込まれていてもよい。
【００６４】
上記実施形態では、対象物はレトルト食品であったが、これに限定されない。例えば、缶詰等の他の包装食品や食品以外のもの、例えば、医薬品、注射器などの医療機器等であってもよい。
【００６５】
上記実施形態では、対象物は食品であり、対象物を加熱する際の加熱条件を、殺菌とたんぱく質の分解との２つの観点で評価したが、これに限定されない。例えば、対象物が食品である場合、殺菌以外の評価の観点として、ビタミンの分解、酵素の分解、テクスチャーや色の劣化などが挙げられる。また、対象物が医薬品や医療機器等である場合、殺菌以外の評価の観点として、医薬品に含まれる成分や医療機器を構成する材料の変質が挙げられる。
【産業上の利用可能性】
【００６６】
本発明のシミュレーション方法は、レトルト食品、缶詰、医薬品、医療機器等の加熱による温度変化を計算する際に利用することができる。
【符号の説明】
【００６７】
１００…対象物、ＡＲ１…第１領域、ＡＲ２…第２領域、Ｐ…中心点
【図１】
2017215191.tif 000003
【図２】
2017215191.tif 000004
【図３】
2017215191.tif 000005