(57)【要約】【課題】対象物の実際の温度変化に近似した推定温度を演算できる加熱シミュレーション方法を提供する。【解決手段】加熱シミュレーション方法では、設定された加熱条件下における加熱対象物の温度変化を演算により求める。当該加熱シミュレーション方法は、加熱対象物の物性値として第１物性値と第２物性値とを設定し、第１物性値から前記第２物性値へと加熱対象物の物性値が変化する温度を変換温度として設定し、加熱対象物の演算された温度が変換温度未満の場合には、加熱対象物の物性値として第１物性値を用いて加熱対象物の温度変化を演算し、加熱対象物の演算された温度が変換温度以上の場合には、加熱対象物の物性値として第２物性値を用いて加熱対象物の温度変化を演算する。【選択図】図１
2017212913.tif 000002
【特許請求の範囲】
【請求項１】
設定された加熱条件下における加熱対象物の温度変化を演算により求める加熱シミュレーション方法であって、
前記加熱対象物の物性値として第１物性値と第２物性値とを設定し、
前記第１物性値から前記第２物性値へと前記加熱対象物の物性値が変化する温度を変換温度として設定し、
前記加熱対象物の演算された温度が前記変換温度未満の場合には、前記加熱対象物の物性値として前記第１物性値を用いて前記加熱対象物の温度変化を演算し、
前記加熱対象物の演算された温度が前記変換温度以上の場合には、前記加熱対象物の物性値として前記第２物性値を用いて前記加熱対象物の温度変化を演算する
加熱シミュレーション方法。
【請求項２】
前記演算中の前記加熱対象物の温度が前記変換温度以上となった場合、その後、演算された温度に関わりなく前記加熱対象物の物性値として前記第２物性値を用いて前記加熱対象物の温度変化を演算する
請求項１に記載の加熱シミュレーション方法。
【請求項３】
演算装置に、設定された加熱条件下における加熱対象物の温度変化の演算を実行させるための加熱シミュレーションプログラムであって、
演算装置に、
前記加熱対象物の物性値として第１物性値と第２物性値との設定を受け付けるステップと、
前記第１物性値から前記第２物性値へと前記加熱対象物の物性値が変化する温度を変換温度として設定を受け付けるステップと、
前記加熱対象物の演算された温度が前記変換温度以上であるか否かを判定するステップと、
前記加熱対象物の演算された温度が前記変換温度未満の場合には、前記加熱対象物の物性値として前記第１物性値を用いて前記加熱対象物の温度変化を演算するステップと、
前記加熱対象物の演算された温度が前記変換温度以上の場合には、前記加熱対象物の物性値として前記第２物性値を用いて前記加熱対象物の温度変化を演算するステップと
を実行させる加熱シミュレーションプログラム。
【請求項４】
請求項３に記載の加熱シミュレーションプログラムを内蔵した記憶媒体を含み、前記加熱シミュレーションプログラムを演算装置が実行させる加熱シミュレーション装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、加熱中の処理対象物の温度変化を計算する加熱シミュレーション方法、加熱シミュレーションプログラム、及びこのプログラムを内蔵した記憶媒体を含む加熱シミュレーション装置に関する。
【背景技術】
【０００２】
従来、缶詰やレトルト食品といった包装食品を製造する際に、製造された包装食品の加熱殺菌が行われている。特定の加熱条件が食品の殺菌に適しているか否かは、一般に、温度と時間との関係で表される殺菌値であるＦ値により評価されている。特定の加熱条件下における食品の温度履歴が定められたＦ値を満たす場合には、加熱条件が食品の殺菌に適していると評価できる。例えば、レトルト食品のＦ値は、食品衛生法により、１２０．０℃４分相当以上とされている。
【０００３】
食品の温度履歴は、加熱中の食品の温度をセンサー等により実測することでも確認できるが、形状や大きさの異なる食品に対する温度履歴を実測により確認するためには、コスト及び時間が必要である。このようなコスト及び時間を低減するために、コンピュータを用いて加熱中の食品の推定温度を演算するシミュレーション方法が普及している（例えば、特許文献１参照）。食品の推定温度は、例えば、ＡＴＳ法（Ａｍｂｉｅｎｔ  Ｔｅｍｐｅｒａｔｕｒｅ  Ｓｌｉｄｅ  ｍｅｔｈｏｄ）により導出される演算式により演算される（例えば、非特許文献１参照）。具体的には、食品の加熱条件及び食品の物性値（例えば、伝熱係数等）がシミュレーション装置に設定され、この加熱条件及び物性値に基づき、シミュレーション装置が、演算式により推定温度を演算する。なお、食品以外のもの、例えば、医薬品のような処理対象物に対する加熱時の推定温度も、同様のシミュレーションにより演算される。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第３０７１４１２号公報
【非特許文献】
【０００５】
【非特許文献１】向井勇、外１名、「温度履歴曲線の相似関係によるＡＴＳ法の理論的課題の解明」、日本食品工学会誌、ｖｏｌ．１６，Ｎｏ．３、ｐｐ２０９−２１７、Ｓｅｐ．２０１５
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上述のように、シミュレーションによって対象物の推定温度を演算するためには、まず、対象物の物性値を設定する必要があるが、対象物の物性値が不明である場合、仮設定した物性値を用いてシミュレーションを行い、トライアルアンドエラーの手法を用いて、物性値を求める。具体的には、対象物の物性値を予測し、この物性値の予測値を用いてシミュレーションを行って推定温度を演算し、演算した推定温度と実測温度とを比較し、これらの差異が縮まるように物性値の予測値を修正し、修正後の物性値を用いて再度シミュレーションを行って求めた推定温度と対象物の実測温度とを比較し、これらの差異がさらに縮まるように物性値を修正するという一連の手法を、推定温度と実測温度との差異が十分に縮まるまで繰り返す。この差異が十分に小さくなった際の物性値が、対象物の推定温度を演算するための物性値として用いられる。しかしながら、対象物によっては、物性値の修正を繰り返しても、推定温度と実測温度との差異が十分に縮まらないことがあった。
【０００７】
本発明は、斯かる実情に鑑み、対象物の実際の温度変化に近似した推定温度を演算できるシミュレーション方法を提供するものである。
【課題を解決するための手段】
【０００８】
発明者が上記課題について鋭意検討したところ、上記のような対象物においては、加熱によって、当該対象物がその特有の温度になると物性が変化しており、物性の変化を考慮したシミュレーションを行うことにより、シミュレーションにより求めた対象物の推定温度と実測温度とを近似させうることを見出した。
【０００９】
本発明にかかる加熱シミュレーション方法は、設定された加熱条件下における加熱対象物の温度変化を演算により求める加熱シミュレーション方法であって、前記加熱対象物の物性値として第１物性値と第２物性値とを設定し、前記第１物性値から前記第２物性値へと前記加熱対象物の物性値が変化する温度を変換温度として設定し、前記加熱対象物の演算された温度が前記変換温度未満の場合には、前記加熱対象物の物性値として前記第１物性値を用いて前記加熱対象物の温度変化を演算し、前記加熱対象物の演算された温度が前記変換温度以上の場合には、前記加熱対象物の物性値として前記第２物性値を用いて前記加熱対象物の温度変化を演算することを特徴とする。
【００１０】
上記加熱シミュレーション方法では、対象物の推定温度が変換温度未満の場合に第１物性値に基づき温度変化が演算され、対象物の推定温度が変換温度以上の場合に第２物性値に基づき温度変化が演算される。従って、上記加熱シミュレーション方法では、対象物の実際の温度変化により近似した推定温度が演算される。
【００１１】
本発明にかかる加熱シミュレーション方法の一態様として、前記演算中の前記加熱対象物の温度が前記変換温度以上となった場合、その後、演算された温度に関わりなく前記加熱対象物の物性値として前記第２物性値を用いて前記加熱対象物の温度変化を演算してもよい。
【００１２】
上記加熱シミュレーション方法では、対象物の推定温度が変換温度に到達したことがあると、第２物性値に基づき、対象物の温度変化が演算される。そのため、この加熱シミュレーション方法は、対象物の物性が一旦変化すれば、さらなる加熱や冷却を行っても、対象物の物性は元に戻らない場合に適している。従って、上記加熱シミュレーション方法では、例えば、材料に液卵を含む卵豆腐等を対象物としても、対象物の実際の温度変化により近似した推定温度が演算される。
【００１３】
本発明にかかる加熱シミュレーションプログラムは、演算装置に、設定された加熱条件下における加熱対象物の温度変化の演算を実行させるための加熱シミュレーションプログラムであって、演算装置に、前記加熱対象物の物性値として第１物性値と第２物性値との設定を受け付けるステップと、前記第１物性値から前記第２物性値へと前記加熱対象物の物性値が変化する温度を変換温度として設定を受け付けるステップと、前記加熱対象物の演算された温度が前記変換温度以上であるか否かを判定するステップと、前記加熱対象物の演算された温度が前記変換温度未満の場合には、前記加熱対象物の物性値として前記第１物性値を用いて前記加熱対象物の温度変化を演算するステップと、前記加熱対象物の演算された温度が前記変換温度以上の場合には、前記加熱対象物の物性値として前記第２物性値を用いて前記加熱対象物の温度変化を演算するステップとを実行させることを特徴とする。
【００１４】
上記プログラムでは、対象物の推定温度が変換温度未満の場合に第１物性値に基づき温度変化が演算され、対象物の推定温度が変換温度以上の場合に第２物性値に基づき温度変化が演算される。従って、上記プログラムでは、対象物の実際の温度変化により近似した推定温度が演算される。
【００１５】
本発明にかかる加熱シミュレーション装置は、加熱シミュレーションプログラムを内蔵した記憶媒体を含み、前記加熱シミュレーションプログラムを演算装置が実行させることを特徴とする。
【００１６】
上記装置では、対象物の推定温度が変換温度未満の場合に第１物性値に基づき温度変化が演算され、対象物の推定温度が変換温度以上の場合に第２物性値に基づき温度変化が演算される。従って、上記装置では、対象物の実際の温度変化により近似した推定温度が演算される。
【発明の効果】
【００１７】
本発明のシミュレーション方法によれば、対象物の実際の温度変化に近似した推定温度を演算できる。
【図面の簡単な説明】
【００１８】
【図１】図１は、本発明の一実施形態にかかる推定温度を演算するフローチャート図である。
【図２】図２は、本発明の一実施形態にかかる対象物の実測温度履歴及び推定温度履歴を示すグラフである。
【発明を実施するための形態】
【００１９】
以下、本発明の加熱シミュレーション方法について、添付図面を参酌して説明する。本実施形態のシミュレーション方法は、シミュレーション方法を実行できるプログラムを内蔵した記憶媒体と、プログラムを実行させる演算装置（ＣＰＵ）とを含むシミュレーション装置で実施される。シミュレーションの結果は、例えば、シミュレーション装置に設けられたディスプレイに表示されるとともに、メモリーカードに保存される。
【００２０】
本発明の加熱シミュレーション方法は、設定された加熱条件下における加熱対象物の温度変化を演算により求めるものである。本実施形態の加熱シミュレーション方法は、加熱により物性が変化する対象物、特に、加熱により比熱が変化するような対象物の推定温度を演算するために適している。そのような対象物としては、例えば、卵豆腐が挙げられる。この場合、加熱開始時点における液卵を含む材料液は、加熱により卵豆腐に変化する。このように、対象物の状態が加熱により液体状からゲル状に変化することで、対象物の比熱も変化し、対象物の推定温度を演算するための物性値も変化する。なお、対象物が液卵のように加熱により変性するたんぱく質を含む場合、加熱により対象物の物性が一旦変化すると、冷却しても元の物性に戻ることは無い。
【００２１】
以上を踏まえて、本発明の加熱シミュレーション方法は、加熱対象物の物性値として第１物性値と第２物性値とを設定し、第１物性値から第２物性値へと加熱対象物の物性値が変化する温度を変換温度として設定し、加熱対象物の演算された温度が変換温度未満の場合には、加熱対象物の物性値として前記第１物性値を用いて加熱対象物の温度変化を演算し、加熱対象物の演算された温度が変換温度以上の場合には、加熱対象物の物性値として第２物性値を用いて加熱対象物の温度変化を演算する。以下、本実施形態の加熱シミュレーション方法の具体的なステップについて説明する。
【００２２】
本実施形態の加熱シミュレーション方法（以下、シミュレーション方法と称する）は、図１のフローチャート図に示すように、対象物の第１、第２物性値、及び変換温度Ｔｃが加熱シミュレーション装置（以下、シミュレーション装置と称する）に設定されるステップ（Ｓ０１）と、加熱条件がシミュレーション装置に設定され、且つ、シミュレーション装置が状態フラグＳに０をセットし、ｎに１をセットするステップ（Ｓ０２）と、シミュレーション装置が、加熱条件から雰囲気温度Ｔｗｎを読み込むステップ（Ｓ０３）と、シミュレーション装置が「状態フラグＳが１か否か」を判定するステップ（Ｓ０４）と、シミュレーション装置が「推定温度Ｔｎ−１が変換温度Ｔｃ以上か否か」を判定するステップ（Ｓ０５）と、シミュレーション装置が第１物性値に基づき推定温度Ｔｎを演算するステップ（Ｓ０６）と、シミュレーション装置が状態フラグＳに１をセットするステップ（Ｓ０７）と、シミュレーション装置が第２物性値に基づき推定温度Ｔｎを演算するステップ（Ｓ０８）と、シミュレーション装置がｎに１を加算するステップ（Ｓ０９）と、シミュレーション装置が「Δｔ×ｎが加熱時間以上か否か」を判定するステップ（Ｓ１０）とを含む。状態フラグＳは、対象物の推定温度が変換温度に到達したことがあることを示すフラグである。Δｔは単位時間を示し、ｎは時間刻みのｎ番目を示す。なお、本実施形態では、対象物全体の推定温度は均一であると仮定する。
【００２３】
以下、本実施形態のシミュレーション方法に含まれる各ステップについて、順に説明する。本実施形態では、対象物の推定温度を演算する方法として、対象物の中心部の温度を演算する方法であるＡＴＳ法（Ａｍｂｉｅｎｔ  Ｔｅｍｐｅｒａｔｕｒｅ  Ｓｌｉｄｅ  ｍｅｔｈｏｄ）を採用した場合について説明する。なお、ここでいう中心部とは、物体の物理的中心ではなく、温度の上昇又は下降の最も遅れる部分を指す。
【００２４】
そこで、まず、ＡＴＳ法の概略について説明する。ＡＴＳ法では、雰囲気からの伝熱により対象物の温度が上昇することを前提とする。また、「雰囲気が対象物に与える熱量」と「対象物が雰囲気から受け取る熱量」とが一致することを前提として計算が行われる。ここで、上述したように単位時間をΔｔとし、対象物の中心点を中心点Ｐとし、雰囲気温度及び対象物の中心点温度をそれぞれＴｗｎ、Ｔｐｎとし（下添え字ｎは、上述のｎと同様に時間刻みのｎ番目を示す）、対象物の表面温度が雰囲気温度Ｔｗｎと一致することとし、対象物の内部の温度勾配を直線とし、対象物の表面から中心点までの距離をＬとし、対象物の熱伝導率をｋとし、対象物の表面積をＡとすると、単位時間当たりに「雰囲気が対象物に与える熱量」即ち下記式の左辺が得られる。
【００２５】
また、対象物の体積をＶとし、対象物の密度をρとし、対象物の比熱をｃｐとし、中心点温度Ｔｐｎは体積平均温度Ｔｐ＊ｎと一致するものとすると、単位時間当たりに「対象物が雰囲気から受け取る熱量」即ち下記式の右辺が得られる。
ｋＡ（Ｔｗｎ−１−Ｔｐｎ−１）Δｔ/Ｌ＝Ｖρｃｐ（Ｔｐ＊ｎ−Ｔｐ＊ｎ−１)
【００２６】
上記式においてｋ/（ρｃｐ)を熱拡散係数αとして整理すると、下記式が得られる。
Ｔｐ＊ｎ＝Ｔｐ＊ｎ−１＋αΔｔ（Ｔｗｎ−１−Ｔｐｎ−１）／Ｌ２
【００２７】
実際には、中心点温度Ｔｐｎは体積平均温度Ｔｐ＊ｎと異なるため、中心点温度Ｔｐｎと体積平均温度Ｔｐ＊ｎとのずれ率βを用いて、Ｔｐｎ＝Ｔｐ＊ｎβとすると、下記式が得られる。
Ｔｐｎ＝Ｔｐｎ−１＋αβΔｔ（Ｔｗｎ−１−Ｔｐｎ−１）／Ｌ２
【００２８】
さらに、αβΔｔ／Ｌ２を伝熱係数τとすると、下記式が得られる。下記式では、単位時間あたりに、対象物の表面温度Ｔｗｎ−１と中心点温度Ｔｐｎ−１との差分に伝熱係数τを乗じた分だけ、中心点温度が上昇することが表されている。
Ｔｐｎ＝Ｔｐｎ−１＋τ（Ｔｗｎ−１−Ｔｐｎ−１）
【００２９】
ところが、実際の対象物では表面と中心点とが離れているため、中心点の温度変化が表面の温度変化よりも遅れるが、この遅れは上記式に反映されていない。これに対して、上記式において、表面温度として、雰囲気温度Ｔｗｎ−１の代わりに、雰囲気温度Ｔｗｎ−１の時間変化をδだけ遅らせた仮想的な雰囲気温度（ｔ−δ）Ｔｗｎ−１を用いると、下記式が得られる。
Ｔｐｎ＝Ｔｐｎ−１＋τ１（（ｔ−δ）Ｔｗｎ−１−Ｔｐｎ−１）
【００３０】
従って、ＡＴＳ法を採用する場合には、ステップＳ０１において、物性値として伝熱係数τ及び遅れ時間δを設定する。具体的には、ステップＳ０１では、例えば、ユーザーの入力により、第１物性値（第１伝熱係数τ１及び第１遅れ時間δ１）、第２物性値（第２伝熱係数τ２及び第２遅れ時間δ２）、及び変換温度Ｔｃがシミュレーション装置に設定される。伝熱係数τ及び遅れ時間δ１は、対象物を構成する材料、対象物の形状、対象物の状態により定められる。変換温度Ｔｃは、第１、第２物性値を変換させる温度である。ユーザーは、第１物性値（第１伝熱係数τ１及び第１遅れ時間δ１）、第２物性値（第２伝熱係数τ２及び第２遅れ時間δ２）、及び変換温度Ｔｃのデータを所持していれば、このデータをシミュレーション装置に入力する。ユーザーは、第１物性値（第１伝熱係数τ１及び第１遅れ時間δ１）、第２物性値（第２伝熱係数τ２及び第２遅れ時間δ２）、及び変換温度Ｔｃのデータを所持していなければ、これらを求めてシミュレーション装置に入力する。
【００３１】
第１物性値、第２物性値、及び、変換温度Ｔｃを求める際には、まず、変換温度Ｔｃを求めた後、対象物の推定温度が変換温度Ｔｃ未満の場合に対応する第１物性値を求め、さらに対象物の推定温度が変換温度Ｔｃ以上の場合に対応する第２物性値を求める。以下、変換温度Ｔｃを求める方法、及び、第１、第２物性値を求める方法について、順に説明する。
【００３２】
変換温度Ｔｃは、実測温度と、物性値の予測値（伝熱係数τ及び遅れ時間δの予測値）に基づき演算した推定温度とを比較して求められる。これについて、図２のグラフを用いて説明する。このグラフは、対象物の加熱条件に対する対象物の温度変化を示すグラフである。グラフの横軸は時間（ｓｅｃ）を示し、グラフの縦軸は温度（℃）を示す。一点鎖線は、加熱条件として設定される雰囲気温度Ｔｗｎを示す。本実施形態では、加熱時の雰囲気温度Ｔｗｎは一定であり、冷却時の雰囲気温度Ｔｗｎも一定である。実線は、対象物の実測温度をプロットして得られる実測温度履歴を示す。点線は、物性値の予測値（伝熱係数τ及び遅れ時間δの予測値）に基づき演算された推定温度をプロットして得られる推定温度履歴を示す。変換温度Ｔｃは、このグラフのように、加熱時の実測温度履歴が上に凸の形状を採り、推定温度履歴が実測温度よりも下側に位置する場合、推定温度と実測温度との差異が最大となる時刻ｔｃにおける実測温度として求められる。ここでいう上側はグラフにおける高温度側、下側はグラフにおける低温度側を言う。
【００３３】
なお、物性値の予測値によっては、推定温度履歴が実測温度履歴よりも上側に位置することがある。この場合、推定温度履歴が実測温度履歴よりも下側に位置するように、物性値の予測値を修正し、修正後の物性値の予測値に基づき推定温度を演算し、修正後の推定温度履歴と実測温度履歴との差異により変換温度Ｔｃを求める。
【００３４】
また、対象物によっては、加熱時の実測温度履歴が下に凸の形状を採ることもある。この場合、推定温度履歴が実測温度履歴よりも上側に位置するように、物性値の予測値を選択し、この物性値の予測値に基づき推定温度を演算すればよい。この場合も、変換温度Ｔｃは、推定温度と実測温度との差異が最大となる時刻ｔｃにおける実測温度として求められる。
【００３５】
第１物性値（第１伝熱係数τ１及び第１遅れ時間δ１）を求める際には、トライアルアンドエラー、即ち、まず、第１物性値について、ユーザーが、適当な伝熱係数τ及び遅れ時間δを予測し、この予測値に基づきＡＴＳ法により推定温度Ｔｎを演算し、この推定温度Ｔｎと実測温度とを比較して、伝熱係数τ及び遅れ時間δ１の予測値を修正し、修正した伝熱係数τ及び遅れ時間δに基づき推定温度Ｔｎを再度演算し、推定温度Ｔｎが実測温度に近づくまでこれを繰り返すことにより適切な第１物性値（第１伝熱係数τ１及び第１遅れ時間δ１）を求めることができる。なお、対象物の実測温度は、例えば、サーミスタのような温度検出センサーを用いて測定される。
【００３６】
以下、具体的な求め方について説明する。例えば、伝熱係数τが不明である場合、予測した伝熱係数τ’１に基づき推定温度Ｔ’ｎを演算し、推定温度Ｔ’ｎと実測温度との差異を計算する。推定温度Ｔ’ｎと実測温度との差異が所望の値以下である場合、伝熱係数τ’１を第１伝熱係数τ１として採用する。
【００３７】
推定温度Ｔ’ｎと実測温度との差異が所望の値よりも大きい場合、伝熱係数τ’１に近い四つの伝熱係数τ’’１…を設定し、四つの伝熱係数τ’’１…各々に基づき推定温度Ｔ’’ｎ…を演算し、四つの推定温度Ｔ’’ｎ…と実測温度との差異を計算する。
【００３８】
推定温度Ｔ’’ｎ…と実測温度との差異の最小値が、推定温度Ｔ’ｎと実測温度との差異よりも小さく、且つ、所望の値よりも小さい場合、推定温度Ｔ’’ｎ…と実測温度との差異の最小値の演算に用いた伝熱係数τ’’１を第１伝熱係数τ１として採用する。
【００３９】
推定温度Ｔ’’ｎ…と実測温度との差異の最小値が、推定温度Ｔ’ｎと実測温度との差異よりも小さく、且つ、所望の値よりも大きい場合、推定温度Ｔ’’ｎ…と実測温度との差異の最小値の演算に用いた伝熱係数τ’’１に近い四つの伝熱係数τ’’’１…を再度設定し、四つの伝熱係数τ’’’１…各々に基づき推定温度Ｔ’’’ｎ…を演算する等の一連の計算を繰り返す。
【００４０】
このように、一連の計算を繰り返すトライアルアンドエラーの手法を実行することで伝熱係数τ’１を求めることができる。
【００４１】
第２物性値についても、上述したトライアルアンドエラーの手法により求める。
【００４２】
ステップＳ０２では、例えば、ユーザーの入力により、シミュレーション装置に対象物の加熱条件が設定される。加熱条件は、例えば、雰囲気温度Ｔｗｎ、加熱時間である。また、シミュレーション装置は、状態フラグＳに０をセットし、ｎに１をセットする。
【００４３】
ステップＳ０３では、シミュレーション装置が、加熱条件から雰囲気温度Ｔｗｎを読み込み、ステップＳ０６又はステップＳ０８で用いる雰囲気温度Ｔｗｎをセットする。例えば、雰囲気温度Ｔｗｎが雰囲気温度Ｔｗｎ−１と異なる場合、ステップＳ０６又はステップＳ０８で用いる雰囲気温度は、雰囲気温度Ｔｗｎに変更される。
【００４４】
ステップＳ０４において、シミュレーション装置が、「状態フラグＳが１か否か」を判定した後、ステップＳ０５では、シミュレーション装置が「推定温度Ｔｎ−１が変換温度Ｔｃ以上であるか否か」を判定する。ステップＳ０４では、対象物の推定温度が変換温度に到達したことがあるか否かを判定している。上述のように、状態フラグＳは、対象物の推定温度が変換温度に到達したことがあることを示すため、ステップＳ０５でＹｅｓと判定されたときに１がセットされる。ステップＳ０５では、推定温度Ｔｎ−１が、変換温度Ｔｃ未満であるか、変換温度Ｔｃ以上であるかを判定している。
【００４５】
ステップＳ０４においてＮｏであって、ステップＳ０５においてＮｏ、即ち、対象物の推定温度が変換温度に到達したことが無く、且つ、推定温度Ｔｎ−１が変換温度Ｔｃ未満であるとき、ステップＳ０６では、シミュレーション装置が、ＡＴＳ法により導出される下記式により、第１物性値（第１伝熱係数τ１及び第１遅れ時間δ１）に基づき、推定温度Ｔｎを演算する。
Ｔｐｎ＝Ｔｐｎ−１＋τ１（（ｔ−δ１）Ｔｗｎ−１−Ｔｐｎ−１）
【００４６】
ステップＳ０４においてＮｏであって、ステップＳ０５においてＹｅｓ、即ち、対象物の推定温度が変換温度に到達したことが無く、且つ、推定温度Ｔｎ−１が変換温度Ｔｃ以上であるとき、ステップＳ０７で状態フラグＳに１がセットされた後、ステップＳ０８では、シミュレーション装置が、ＡＴＳ法により導出される下記式により、第２物性値（第２伝熱係数τ２及び第２遅れ時間δ２）に基づき、推定温度Ｔｎを演算する。
Ｔｐｎ＝Ｔｐｎ−１＋τ２（（ｔ−δ２）Ｔｗｎ−１−Ｔｐｎ−１）
【００４７】
ステップＳ０４においてＹｅｓ、即ち、対象物の推定温度が変換温度Ｔｃに到達したことがあるとき、ステップＳ０８では、シミュレーション装置が、ＡＴＳ法により導出される下記式により、第２物性値（第２伝熱係数τ２及び第２遅れ時間δ２）に基づき、推定温度Ｔｎを演算する。
Ｔｐｎ＝Ｔｐｎ−１＋τ２（（ｔ−δ２）Ｔｗｎ−１−Ｔｐｎ−１）
【００４８】
即ち、演算中の対象物の温度が変換温度Ｔｃ以上となった場合、その後、演算された温度に関わりなく対象物の物性値として第２物性値（第２伝熱係数τ２及び第２遅れ時間δ２）を用いて対象物の温度変化を演算する。
【００４９】
ステップＳ０９においてシミュレーション装置がｎに１を加算した後、ステップＳ１０では、シミュレーション装置が「Δｔ×ｎが加熱時間以上か否か」を判定する。
【００５０】
ステップＳ１０においてＮｏ、即ち、Δｔ×ｎが加熱時間未満であるとき、シミュレーション装置は、ステップＳ０３〜ステップＳ０９を繰り返す。
【００５１】
ステップＳ１０においてＹｅｓ、即ち、Δｔ×ｎが加熱時間以上であるとき、シミュレーション装置は、シミュレーション方法の実施を終了する。
【００５２】
このように、本実施形態に係るシミュレーション方法を実施することで、対象物の実際の温度変化に近似した対象物の推定温度を演算することができる。この推定温度は、例えば、食品の加熱条件を殺菌の観点で評価する場合の指標として用いられる。以下、殺菌の観点で、推定温度を指標とした食品の加熱条件の評価について説明する。
【００５３】
特定の加熱条件が食品の殺菌の観点で適切か否かは、一般に、加熱時の食品の温度履歴が、殺菌評価積算値であるＦ値を満たすか否かで評価される。また、殺菌の観点では、雰囲気からの熱が伝わりにくく、最も殺菌されにくい食品の中心部の温度を評価することが好ましい。そのため、殺菌の観点で加熱条件を評価するには、食品の中心部において推定温度Ｔｎから推定温度履歴Ｔを求めて、この推定温度履歴Ｔが定められたＦ値に相当するか否かを評価する。以下、本実施形態の効果をまとめて説明する。
【００５４】
本実施形態のシミュレーション方法では、対象物の推定温度が変換温度Ｔｃ未満である場合には、第１物性値（第１伝熱係数τ１及び第１遅れ時間δ１）に基づき推定温度Ｔｎが演算され、対象物の推定温度が変換温度Ｔｃ以上である場合には、第２物性値（第２伝熱係数τ２及び第２遅れ時間δ２）に基づき推定温度が演算される。従って、上記シミュレーション方法では、対象物の実際の温度変化により近似した推定温度が演算される。
【００５５】
本実施形態のシミュレーション方法は、加熱条件が設定されるステップＳ０２と、推定温度と変換温度Ｔｃとの判定をするステップＳ０５との間に、対象物の推定温度が変換温度に到達したことがあるか否かを判定するステップＳ０４を含む。また、対象物の推定温度が変換温度Ｔｃに到達したことがある（状態フラグＳが１である）とき、ステップＳ０５〜Ｓ０８を実行する代わりに、第２物性値に基づき、対象物における推定温度Ｔｎを演算する。即ち、演算中の対象物の温度が変換温度Ｔｃ以上となった場合、その後、演算された温度に関わりなく、対象物の物性値として第２物性値を用いて対象物の温度変化を演算する。このシミュレーション方法は、対象物の物性が一旦変化すれば、さらなる加熱や冷却を行っても、対象物の物性は元に戻らない場合に適している。従って、上記シミュレーション方法では、例えば、材料に液卵を含む卵豆腐等を対象物としても、対象物の実際の温度変化により近似した推定温度が演算される。
【００５６】
なお、本発明にかかるシミュレーション方法は、上記実施形態の構成に限定されるものではなく、本発明の要旨を逸脱しない範囲で種々の変更が可能である。
【００５７】
上記実施形態では、変換温度が一つであり、第１、第２物性値により推定温度を演算していたが、これに限定されない。変換温度を二つ以上とし、第１、第２、第３物性値、…という三以上の複数の物性値により推定温度を演算してもよい。
【００５８】
上記実施形態では、第１、第２物性値として、伝熱係数τ及び遅れ時間δの両方を用いていたが、これに限定されない。例えば、変換温度の前後で、伝熱係数τが変化するものの遅れ時間δが変化しない場合、第１、第２物性値として、伝熱係数τのみを用いてもよい。
【００５９】
上記実施形態では、変換温度は、実測温度と推定温度との温度差が最大となる温度としていたが、これに限定されない。例えば、変換温度は、実測温度に対する実測温度と推定温度との差異の比率の平均値が最大となる温度としてもよい。即ち、差異の比率は、以下の数式により求めることができる。なお、（実測温度−推定温度）が負の値となる場合は、絶対値を用いる。
｛（実測温度−推定温度）／実測温度｝×１００
【００６０】
上記実施形態では、一つの基準により変換温度を求めていたが、これに限定されない。例えば、実測温度と推定温度との差異が最大となる第１変換温度と、及び実測温度に対する実測温度と推定温度との差異の比率の平均値が最大となる第２変換温度をそれぞれ求めて、第１、第２変換温度に基づき推定温度をそれぞれ演算し、第１、第２変換温度のうち演算された推定温度がより実測温度と近い変換温度を用いてもよい。
【００６１】
上記実施形態では、第１物性値、第２物性値、及び、変換温度Ｔｃを求める際に、変換温度Ｔｃを求めた後、変換温度Ｔｃ未満の温度に対応する第１物性値を求めて、変換温度Ｔｃ以上の温度に対応する第２物性値を求めた。しかしながら、これらの求める順序は上述の順序と異なってもよいし、第１、第２物性値を同時に求めてもよい。
【００６２】
上記実施形態では、加熱時の雰囲気温度Ｔｗｎは一定であり、冷却時の雰囲気温度Ｔｗｎも一定であったが、これに限定されない。加熱時及び冷却時の雰囲気温度をそれぞれ複数種類ずつ設定してもよい。
【００６３】
上記実施形態では、対象物は卵豆腐であったが、これに限定されない。例えば、レトルト食品、缶詰等の他の包装食品や食品以外のもの、例えば、医薬品等であってもよい。
【００６４】
上記実施形態では、対象物の加熱過程において変換温度が存在していたが、これに限定されない。例えば、対象物の冷却過程において変換温度が存在することも考えられる。このような対象物としては、冷却により状態が変化するゼリー等が挙げられる。ゼリーのような対象物については、シミュレーション方法において、対象物の推定温度が変換温度以上であるか否かを判定し、それぞれの場合において異なる物性値により推定温度を演算すればよい。
【００６５】
上記実施形態では、対象物の物性が一旦変化すれば、さらなる加熱や冷却を行っても、対象物の物性は変化しないこととしたが、これに限定されない。例えば、対象物が、ゼリーである場合、加熱開始時にゲル状であり、加熱中に液体状に変化し、冷却後にゲル状に戻る。この場合、シミュレーション方法では、図１のフローチャート図における、状態フラグの判定ステップＳ０４を実行せず、推定温度が変換温度Ｔｃに到達した場合でも、推定温度と変換温度Ｔｃとの判定により、第１物性値又は第２物性値に基づき、推定温度Ｔｎを演算すればよい。従って、この場合には、演算された対象物の推定温度が変換温度Ｔｃよりも低くなった場合、第１物性値に基づき推定温度Ｔｎが演算される。
【００６６】
上記実施形態では、加熱過程及び冷却過程において、対象物の状態が同じであるため、同一の第２物性値に基づき推定温度を演算したが、これに限定されない。例えば、対象物の状態が同じであっても、加熱過程及び冷却過程における対象物の実際の温度変化が異なる場合には、これに適するよう、加熱過程及び冷却過程において異なる物性値に基づき推定温度を演算してもよい。
【００６７】
上記実施形態では、推定温度をＡＴＳ法により演算したが、これに限定されない。例えば、Ｂａｌｌの数式法等の別の方法によって、推定温度を演算してもよい。
【００６８】
上記実施形態のシミュレーション装置は、対象物を加熱する装置と別の装置であってもよいし、対象物を加熱する装置に組み込まれていてもよい。
【産業上の利用可能性】
【００６９】
本発明のシミュレーション方法は、レトルト食品、缶詰、医薬品等の加熱による温度変化を計算する場合に利用することができる。
【符号の説明】
【００７０】
Ｔｃ…変換温度、Ｔｎ…推定温度、Ｔｗｎ…雰囲気温度
【図１】
2017212913.tif 000003
【図２】
2017212913.tif 000004