(57)【要約】      （修正有）【課題】従来の冷却性評価モデルと同等の予測性能を維持しつつ、過渡解析にかかる時間を短縮する。【解決手段】高速炉デブリベッド冷却の解析装置は、高速炉において、原子炉内に設けられたコアキャッチャーに堆積したデブリベッドを分割した第１サイズの複数のメッシュのそれぞれについて、デブリベッドの未沸騰領域における等価熱伝導度と、デブリベッドの沸騰領域における等価熱伝導度よりも大きい模擬的な等価熱伝導度とに基づいて各時間における非定常熱伝導方程式を解くことにより、デブリベッドの未沸騰領域の温度及び沸騰領域の範囲を過渡解析する第１解析部と、特定の時間において解析された沸騰領域の範囲に対応するデブリベッドを第１サイズよりも小さい第２サイズの複数のメッシュに分割し、デブリベッドの等価熱伝導度に基づく定常解析を行うことにより沸騰領域の温度及び水率の分布を解析する第２解析部と、を備える。【選択図】図４
2019211340.tif 000002
【特許請求の範囲】
【請求項１】
コンピュータが実行する、
高速炉において、原子炉内に設けられたコアキャッチャーに堆積した粒子状の燃料デブリであるデブリベッドを分割した第１サイズの複数のメッシュのそれぞれについて、デブリベッドの未沸騰領域における等価熱伝導度と、デブリベッドの沸騰領域における等価熱伝導度よりも大きい、前記沸騰領域における模擬的な等価熱伝導度とに基づいて各時間における非定常熱伝導方程式を解くことにより、前記デブリベッドの未沸騰領域の温度及び沸騰領域の範囲を過渡解析する第１解析ステップと、
前記第１解析ステップにより解析された各時間の沸騰領域の範囲のうち、特定の時間における沸騰領域の範囲に対応する前記デブリベッドを前記第１サイズよりも小さい第２サイズの複数のメッシュに分割し、当該第２サイズの複数のメッシュのそれぞれについて、前記デブリベッドの沸騰領域における等価熱伝導度に基づく定常解析を行うことにより、前記特定の時間における前記沸騰領域の温度及び水率の分布を解析する第２解析ステップと、
を備える高速炉デブリベッド冷却の解析方法。
【請求項２】
前記沸騰領域における模擬的な等価熱伝導度は、前記デブリベッドの沸騰領域の開始温度から所定範囲内において、最低の熱伝導度であり、前記所定範囲外の温度において、前記デブリベッドの等価熱伝導度よりも高い熱伝導度である、
請求項１に記載の高速炉デブリベッド冷却の解析方法。
【請求項３】
コンピュータを、
高速炉において、原子炉内に設けられたコアキャッチャーに堆積した粒子状の燃料デブリであるデブリベッドを分割した第１サイズの複数のメッシュのそれぞれについて、デブリベッドの未沸騰領域における等価熱伝導度と、デブリベッドの沸騰領域における等価熱伝導度よりも大きい、前記沸騰領域における模擬的な等価熱伝導度とに基づいて各時間における非定常熱伝導方程式を解くことにより、前記デブリベッドの未沸騰領域の温度及び沸騰領域の範囲を過渡解析する第１解析部、及び、
前記第１解析部により解析された各時間の沸騰領域の範囲のうち、特定の時間における沸騰領域の範囲に対応する前記デブリベッドを前記第１サイズよりも小さい第２サイズの複数のメッシュに分割し、当該第２サイズの複数のメッシュのそれぞれについて、前記デブリベッドの沸騰領域における等価熱伝導度に基づく定常解析を行うことにより、前記特定の時間における前記沸騰領域の温度及び水率の分布を解析する第２解析部、
として機能させる高速炉デブリベッド冷却の解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、高速炉デブリベッド冷却の解析方法及び解析プログラムに関する。
【背景技術】
【０００２】
次世代ナトリウム冷却高速炉においては、炉心損傷事象の炉容器内終息を目標とした設計がなされている。図１は、炉心損傷事象に伴い発生した燃料デブリが炉容器下部プレナムに設置されたコアキャッチャー上に堆積した状態を説明する図である。炉心損傷事象に伴い発生した溶融燃料は、炉心領域外に排出され、冷却材との相互作用により分散・微粒化する。そして、溶融燃料は、図１に示すように、最終的に下部プレナムに設置されたコアキャッチャー１０１にベッド状に堆積したデブリベッド１０２となる。デブリベッドは、主に数百μｍの粒子で形成され、崩壊熱により発熱する。炉内終息の達成のためには、デブリベッドを長期間安定に冷却する必要があり、冷却性評価モデルの検討が重要となる。
【０００３】
冷却性評価モデルとしては、下記式（１）に示すデブリベッドの等価熱伝導度ｋｅｑを用いた非定常熱伝導方程式を解く冷却性評価モデルが知られている（例えば、非特許文献１参照）。
【数１】
2019211340.tif 000003
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Nakamura, H. et al., “A Method to Calculate Boiling and Dryout in Debris Bed Using a Heat Conduction Code”, in Proceedings of International Topical Meeting on Fast Reactor Safety, Knoxville, Tennessee, U.S.A..1985
【発明の概要】
【発明が解決しようとする課題】
【０００５】
図２は、デブリベッドの沸騰領域における等価熱伝導度を示す図である。図２に示すように、デブリベッドの沸騰領域における等価熱伝導度は、数℃で大きく値が変動する。したがって、式（１）に示す非定常熱伝導方程式を解く冷却性評価モデルにおいて、沸騰領域の正確な温度及び水率分布を予測するためには、非常に小さなメッシュサイズ及びタイムステップを必要とする。
【０００６】
図３は、従来の過渡解析手法を説明するための図である。従来の過渡解析の実施手法では、過渡解析において沸騰領域の正確な温度及び水率分布を予測するため、非常に小さなメッシュサイズ及びタイムステップを用いて計算を行うため、過渡解析の実施に膨大な時間を要していた。
【０００７】
そこで、本発明はこれらの点に鑑みてなされたものであり、従来の冷却性評価モデルと同等の予測性能を維持しつつ、過渡解析にかかる時間を短縮することができる高速炉デブリベッド冷却の解析方法及び解析プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明者らは、デブリベッドにおける沸騰領域が定常的な挙動を示すことに着目し、デブリベッドにおける沸騰領域における等価熱伝導度を大きな値に単純化して非定常熱伝導方程式を解くことにより過渡解析を実施し、過渡解析では直接予測できなくなった沸騰領域のデブリベッドの温度及び水率分布を必要に応じて定常解析から予測することを見出した。
【０００９】
本発明の第１の態様に係る高速炉デブリベッド冷却の解析方法は、コンピュータが実行する、高速炉において、原子炉内に設けられたコアキャッチャーに堆積した粒子状の燃料デブリであるデブリベッドを分割した第１サイズの複数のメッシュのそれぞれについて、デブリベッドの未沸騰領域における等価熱伝導度と、デブリベッドの沸騰領域における等価熱伝導度よりも大きい、前記沸騰領域における模擬的な等価熱伝導度とに基づいて各時間における非定常熱伝導方程式を解くことにより、前記デブリベッドの未沸騰領域の温度及び沸騰領域の範囲を過渡解析する第１解析ステップと、前記第１解析ステップにより解析された各時間の沸騰領域の範囲のうち、特定の時間における沸騰領域の範囲に対応する前記デブリベッドを前記第１サイズよりも小さい第２サイズの複数のメッシュに分割し、当該第２サイズの複数のメッシュのそれぞれについて、前記デブリベッドの沸騰領域における等価熱伝導度に基づく定常解析を行うことにより、前記特定の時間における前記沸騰領域の温度及び水率の分布を解析する第２解析ステップと、を備える。
【００１０】
前記沸騰領域における模擬的な等価熱伝導度は、前記デブリベッドの沸騰領域の開始温度から所定範囲内において、最低の熱伝導度であり、前記所定範囲外の温度において、前記デブリベッドの等価熱伝導度よりも高い熱伝導度であってもよい。
【００１１】
本発明の第２の態様に係る高速炉デブリベッド冷却の解析プログラムは、コンピュータを、高速炉において、原子炉内に設けられたコアキャッチャーに堆積した粒子状の燃料デブリであるデブリベッドを分割した第１サイズの複数のメッシュのそれぞれについて、デブリベッドの未沸騰領域における等価熱伝導度と、デブリベッドの沸騰領域における等価熱伝導度よりも大きい、前記沸騰領域における模擬的な等価熱伝導度とに基づいて各時間における非定常熱伝導方程式を解くことにより、前記デブリベッドの未沸騰領域の温度及び沸騰領域の範囲を過渡解析する第１解析部、及び、前記第１解析部により解析された各時間の沸騰領域の範囲のうち、特定の時間における沸騰領域の範囲に対応する前記デブリベッドを前記第１サイズよりも小さい第２サイズの複数のメッシュに分割し、当該第２サイズの複数のメッシュのそれぞれについて、前記デブリベッドの沸騰領域における等価熱伝導度に基づく定常解析を行うことにより、前記特定の時間における前記沸騰領域の温度及び水率の分布を解析する第２解析部、として機能させる。
【発明の効果】
【００１２】
本発明によれば、従来の冷却性評価モデルと同等の予測性能を維持しつつ、過渡解析にかかる時間を短縮することができるという効果を奏する。
【図面の簡単な説明】
【００１３】
【図１】炉心損傷事象に伴い発生した燃料デブリが炉容器下部プレナムに設置されたコアキャッチャー上に堆積した状態を説明する図である。
【図２】デブリベッドの沸騰領域における等価熱伝導度を示す図である。
【図３】従来の過渡解析手法を説明するための図である。
【図４】本実施形態に係る高速炉デブリベッド冷却の解析手法を説明するための図である。
【図５】本実施形態に係る高速炉デブリベッド冷却の解析装置の構成を示す図である。
【図６】本実施形態に係る第１解析部が非定常熱伝導方程式を解く場合に用いるデブリベッドの沸騰領域における模擬的な等価熱伝導度を示す図である。
【図７】パラメータであるデブリベッド高さに対する最高温度及び最低水率のピーク値を示す図である。
【図８】本感度解析においてドライアウトに至らなかった最大のデブリベッド高さにおける最高温度及び最低水率の過渡変化を示す図である。
【発明を実施するための形態】
【００１４】
［高速炉デブリベッド冷却の解析装置の概要］
図１に示すように、原子炉容器１００の下部プレナムには、炉心損傷事象に伴い発生した燃料デブリを受け止めるコアキャッチャー１０１が設けられている。炉心損傷事象に伴い発生した溶融燃料は、炉心領域外に排出され、冷却材との相互作用により分散・微粒化する。そして、溶融燃料は、最終的にはコアキャッチャー１０１にベッド状に堆積したデブリベッド１０２となる。
【００１５】
図４は、本実施形態に係る高速炉デブリベッド冷却の解析手法を説明するための図である。本実施形態に係る高速炉デブリベッド冷却の解析装置１は、コアキャッチャー１０１に堆積したデブリベッド１０２の温度及び水率を解析する装置である。高速炉デブリベッド冷却の解析装置１は、例えば、コンピュータである。以下の説明において、高速炉デブリベッド冷却の解析装置１を、単に解析装置１という。
【００１６】
本実施形態において、解析装置１は、図４に示す第１の解析として、従来の冷却性評価モデルにおいて用いられているメッシュサイズよりも大きい第１サイズでデブリベッドを分割し、第１サイズの複数のメッシュを生成する。そして、解析装置１は、デブリベッドの未沸騰領域における等価熱伝導度と、従来の冷却性評価モデルにおいて用いられる沸騰領域における等価熱伝導度の値に比べて大きな値に単純化した沸騰領域における模擬的な等価熱伝導度とを用いて非定常熱伝導方程式を解くことにより、各時間におけるデブリベッドの未沸騰領域の温度及び沸騰領域の範囲を過渡解析する。図４に示すグラフにおいて、縦軸はデブリベッドの高さを示しており、横軸は、デブリベッドの温度及び水率を示している。ここで、未沸騰領域は、デブリベッドの沸点未満の温度領域である。また、沸騰領域は、デブリベッドの沸点以上の温度領域であり、デブリベッドがドライアウトとなる温度領域も含むものとする。
【００１７】
続いて、解析装置１は、図４に示す第２の解析として、第１の解析により解析された各時間のそれぞれの沸騰領域の範囲のうち、特定の時間に対応する沸騰領域の範囲のデブリベッドを、第１サイズよりも小さい第２サイズにより複数のメッシュに分割する。そして、解析装置１は、第２サイズの複数のメッシュのそれぞれについて、従来の冷却性評価モデルにおいて用いられるデブリベッドの等価熱伝導度に基づく定常解析を行うことにより沸騰領域の温度及び水率の分布を解析する。
【００１８】
このように、解析装置１は、第１の解析において、従来の冷却性評価モデルのメッシュサイズに比べて大きなメッシュサイズを用いて、従来のタイムステップよりも大きなタイムステップで高速に解析を行うとともに、第２の解析において、各時間のうち特定の時間に対応する沸騰領域の範囲の温度及び水率分布を正確に予測する。これにより、解析装置１は、従来の冷却性評価モデルと同等の予測性能を維持しつつ、過渡解析にかかる時間を短縮することができる。
【００１９】
［解析装置１の構成例］
続いて、解析装置１の構成について説明する。図５は、本実施形態に係る解析装置１の構成を示す図である。
解析装置１は、入力部１１と、表示部１２と、通信部１３と、記憶部１４と、制御部１５とを備える。
【００２０】
入力部１１は、例えば、マウス又はキーボード等により構成されており、解析装置１のユーザから操作入力を受け付ける。
表示部１２は、例えば、液晶ディスプレイや有機ＥＬ（Electro-Luminescence）ディスプレイ等により構成される。表示部１２は、制御部１５の制御に応じて各種情報を表示する。
通信部１３は、例えば、解析装置１が他の装置と通信するための通信インタフェースである。
【００２１】
記憶部１４は、例えば、ＲＯＭ（Read Only Memory）及びＲＡＭ（Random Access Memory）等である。記憶部１４は、解析装置１を機能させるための各種プログラムを記憶する。例えば、記憶部１４は、解析装置１の制御部１５を、後述する第１解析部１５１、第２解析部１５２、及び出力部１５３として機能させる高速炉デブリベッド冷却の解析プログラムを記憶する。
【００２２】
制御部１５は、例えばＣＰＵ（Central Processing Unit）である。制御部１５は、記憶部１４に記憶されている高速炉デブリベッド冷却の解析プログラムを実行することにより、第１解析部１５１、第２解析部１５２、及び出力部１５３として機能する。
【００２３】
第１解析部１５１は、高速炉において、原子炉内に設けられたコアキャッチャー上にデブリベッドが堆積している場合を想定し、コアキャッチャーに堆積した粒子状の燃料デブリであるデブリベッドを分割することにより、第１サイズの複数のメッシュを生成する。
【００２４】
本実施形態では、第１解析部１５１は、コアキャッチャー上の一点に堆積しているデブリベッドを、デブリベッドの高さ方向に第１サイズにより分割することにより、第１サイズの複数のメッシュを生成する。なお、本実施形態では、第１解析部１５１は、デブリベッドを、デブリベッドの高さ方向に第１サイズにより分割することとしたが、これに限らず、デブリベッドを、高さ方向、水平方向（幅方向、奥行方向）のうち、一以上の方向に分割してもよい。例えば、第１解析部１５１は、デブリベッドを高さ方向及び幅方向に分割してもよいし、高さ方向、幅方向及び奥行き方向に分割してもよい。
【００２５】
第１解析部１５１は、生成した第１サイズの複数のメッシュのそれぞれについて、デブリベッドの未沸騰領域における等価熱伝導度と、デブリベッドの沸騰領域における等価熱伝導度よりも大きい、デブリベッドの沸騰領域における模擬的な等価熱伝導度とに基づいて各時間における非定常熱伝導方程式を解くことにより、デブリベッドの未沸騰領域の温度及び沸騰領域の範囲を過渡解析する。
【００２６】
ここで、非定常熱伝導方程式は、以下の式（２）に示されるものとする。式（２）中の（ρＣｐ）Ｂはデブリベッドの熱容量、ＴＢはデブリベッドの温度、ｔは時間、Ｑはデブリベッドの発熱密度、ｋｅｑは等価熱伝導度、ｚはデブリベッドの高さ方向の位置を示している。
【数２】
2019211340.tif 000004
【００２７】
図６は、第１解析部１５１が非定常熱伝導方程式を解く場合に用いるデブリベッドの沸騰領域における模擬的な等価熱伝導度を示す図である。図６における実線及び一点鎖線は模擬的な等価熱伝導度を示しており、破線は、図２に示したデブリベッドの等価熱伝導度を示している。図６に示すように、模擬的な等価熱伝導度は、デブリベッドの沸騰領域の開始温度から所定範囲内を示す範囲Ｒ１内において、最低の熱伝導度である。
【００２８】
また、模擬的な等価熱伝導度は、デブリベッドの沸騰領域の開始温度から所定範囲外の温度を示す範囲Ｒ２内の温度において、図６に破線で示すデブリベッドの等価熱伝導度よりも高い熱伝導度である。デブリベッドの沸騰領域の開始温度から所定範囲外の温度における模擬的な等価熱伝導度は、例えば、同温度範囲におけるデブリベッドの等価熱伝導度のピーク値に対して数倍以上高い値を有している。また、模擬的な等価熱電導度は、図６に一点鎖線で示すように、最低の熱伝導度から、図２に示すデブリベッドの等価熱伝導度よりも高い熱伝導度に移行するまでの間に極端に不連続な値とならないように定義される。
【００２９】
図６に示すように、沸騰領域の広い範囲における模擬的な等価熱伝導度を、デブリベッドの等価熱伝導度に比べて大きい値とすることにより、メッシュサイズを大きくしてもメッシュ間で熱が十分に伝わるという前提で解析することができる。これにより、第１解析部１５１は、デブリベッドの温度及び水率の過渡解析を高速に行うことができる。
【００３０】
第１解析部１５１が各メッシュについて温度を算出した結果、沸騰点よりも低い温度となるメッシュに対応する位置のデブリベッドは未沸騰領域に属し、沸騰点以上の温度となるメッシュに対応する位置のデブリベッドは沸騰領域に属する。また、未沸騰領域の水率は１となり、沸騰領域の水率は１以下となる。ここで、第１解析部１５１は、図６に示す模擬的な等価熱伝導度を用いてデブリベッドの温度及び水率を算出するので、沸騰領域の温度及び水率の精度が低くなってしまう。
【００３１】
沸騰領域において、沸騰点からドライアウト（等価熱伝導度が再びほぼゼロとなる温度）までの温度差は数℃しかなく、その間の等価熱伝導度の平均値は非常に大きいことから、沸騰領域は、過渡においても定常的な挙動を示すと考えられる。そこで、解析装置１は、第２解析部１５２により、第１解析部１５１が解析した沸騰領域の範囲に対して、図２に示すデブリベッドの等価熱伝導度に基づく定常解析を行い、沸騰領域の温度及び水率の分布を解析する。
【００３２】
具体的には、まず、第２解析部１５２は、第１解析部１５１により解析された各時間の沸騰領域の範囲のうち、特定の時間における沸騰領域の範囲に対応するデブリベッドを第１サイズよりも小さい第２サイズの複数のメッシュに分割する。本実施形態では、第１解析部１５１がデブリベッドを高さ方向に分割していることから、第２解析部１５２も、沸騰領域の範囲に対応するデブリベッドを高さ方向に分割する。これにより、第２解析部１５２は、第２サイズの複数のメッシュを生成する。なお、本実施形態では、第２解析部１５２は、デブリベッドを、デブリベッドの高さ方向に第２サイズにより分割することとしたが、これに限らず、デブリベッドを、高さ方向、水平方向（幅方向、奥行方向）のうち、一以上の方向に分割してもよい。例えば、第２解析部１５２は、デブリベッドを高さ方向及び幅方向に分割してもよいし、高さ方向、幅方向及び奥行き方向に分割してもよい。
【００３３】
ここで、特定の時間は、解析装置１が、解析結果をユーザが閲覧可能な情報として出力するタイミングに対応する時間である。例えば、第１解析部１５１が第１時間間隔ごとに解析を行うのに対し、第２解析部１５２は、第１時間よりも長い第２時間間隔ごとに解析を行う。
【００３４】
続いて、第２解析部１５２は、第２サイズの複数のメッシュのそれぞれについて、デブリベッドの沸騰領域における等価熱伝導度に基づく定常解析を行うことにより、特定の時間における沸騰領域の温度及び水率の分布を解析する。上述したとおり、沸騰領域は、過渡においても定常的な挙動を示すと考えられることから、式（２）の左辺のｄＴ／ｄｔがほぼ０となる。第２解析部１５２は、式（２）の左辺の値を０として式（２）を解くことにより、第２サイズの複数のメッシュのそれぞれに対応する温度及び水率を算出する。
【００３５】
出力部１５３は、第２解析部１５２が解析した結果を出力する。例えば、出力部１５３は、第２解析部１５２が解析した結果を示す情報を記憶部１４に記憶させたり、表示部１２に表示させたりする。
【００３６】
［従来手法との比較結果］
続いて、典型的なナトリウム冷却高速炉のデブリベッド条件を用いてデブリベッド高さをパラメータとしたときの、従来の冷却性評価モデルに係る解析手法を用いた解析結果と、本発明に係る解析手法を用いた解析結果とを比較した結果について説明する。
【００３７】
下記の表１は、典型的なナトリウム冷却高速炉のデブリベッド条件を示す。
【表１】
2019211340.tif 000005
【００３８】
以下の式（３）は、従来の解析手法を用いた解析結果と、本発明に係る解析手法を用いた解析結果とを比較する際に用いた、沸騰領域における非定常熱伝導方程式である。
【数３】
2019211340.tif 000006
【００３９】
（ρＣｐ）Ｂはデブリベッドの熱容量、ＴＢはデブリベッドの温度、ｔは時間、εはデブリベッドの空隙、ρｖは蒸気密度、Ｈｌｖは気化潜熱、ｓは水率を示している。ｓはＴＢの関数である。Ｑはデブリベッドの発熱密度、ＫＢは気液移流を除く等価熱伝導度、ｚはデブリベッドの高さ方向の位置、Ｋｌｖは気液移流の等価熱伝導度を示している。
【００４０】
したがって、式（３）の左辺の第１項は顕熱に係わる慣性項，左辺の第２項は蒸発潜熱に係わる慣性項であり、右辺の第１項はデブリベッドの発熱密度、第２項は気液移流を除く伝熱項、第３項は気液移流に伴う伝熱項となる。
【００４１】
解析を行うにあたり、デブリベッドの周囲冷却材温度は、簡易的に４００℃一定とした。図７は、パラメータであるデブリベッド高さに対する最高温度及び最低水率のピーク値を示す図である。図８は、本感度解析においてドライアウトに至らなかった最大のデブリベッド高さ（０．１６［ｍ］）における最高温度及び最低水率の過渡変化を示す図である。図７及び図８に示すように、本発明に係る解析手法による解析結果と、従来の解析手法による解析結果とはよく一致しており、本手法に係る解析手法の妥当性が確認できた。また、従来の解析手法では、数時間から半日程度の時間を要していた計算が、同じ計算機環境にて数分から数十分程度まで短縮でき、計算時間が約１／５０に短縮できた。
【００４２】
［本実施形態における効果］
以上のとおり、本実施形態に係る解析装置１は、原子炉内に設けられたコアキャッチャーに堆積した粒子状の燃料デブリであるデブリベッドを分割した第１サイズの複数のメッシュのそれぞれについて、デブリベッドの未沸騰領域における等価熱伝導度と、デブリベッドの沸騰領域における等価熱伝導度よりも大きい、デブリベッドの沸騰領域における模擬的な等価熱伝導度とに基づいて各時間における非定常熱伝導方程式を解くことにより、デブリベッドの未沸騰領域の温度及び沸騰領域の範囲を過渡解析する。そして、解析装置１は、解析された各時間の沸騰領域の範囲のうち、特定の時間における沸騰領域の範囲に対応するデブリベッドを第１サイズよりも小さい第２サイズの複数のメッシュに分割し、当該第２サイズの複数のメッシュのそれぞれについて、デブリベッドの沸騰領域における等価熱伝導度に基づく定常解析を行うことにより、特定の時間における沸騰領域の温度及び水率の分布を解析する。
【００４３】
このように、解析装置１は、第１の解析において、従来の冷却性評価モデルのメッシュサイズに比べて大きなメッシュサイズを用いて、従来のタイムステップよりも大きなタイムステップで高速に解析を行うとともに、第２の解析において、各時間のうち特定の時間に対応する沸騰領域の範囲の温度及び水率分布を正確に予測する。これにより、解析装置１は、従来の冷却性評価モデルと同等の予測性能を維持しつつ、過渡解析にかかる時間を短縮することができる。
【００４４】
以上、本発明を実施の形態を用いて説明したが、本発明の技術的範囲は上記実施の形態に記載の範囲には限定されず、その要旨の範囲内で種々の変形及び変更が可能である。例えば、上述の実施形態では、高速炉デブリベッドの冷却性の解析を行う例について説明したがこれに限らない。例えば、本発明に係る第１の解析及び第２の解析を、高速炉デブリベッドとは異なる、微小粒子で形成された充填ベッドの冷却性の評価解析に適用してもよい。
【００４５】
この場合、第１解析部１５１は、微小粒子で形成された充填ベッドを分割した第１サイズの複数のメッシュのそれぞれについて、充填ベッドの未沸騰領域における等価熱伝導度と、充填ベッドの沸騰領域における等価熱伝導度よりも大きい、沸騰領域における模擬的な等価熱伝導度とに基づいて各時間における非定常熱伝導方程式を解くことにより、充填ベッドの未沸騰領域の温度及び沸騰領域の範囲を過渡解析する。また、第２解析部１５２は、第１解析部１５１が解析した各時間の沸騰領域の範囲のうち、特定の時間における沸騰領域の範囲に対応する充填ベッドを第１サイズよりも小さい第２サイズの複数のメッシュに分割し、当該第２サイズの複数のメッシュのそれぞれについて、充填ベッドの沸騰領域における等価熱伝導度に基づく定常解析を行うことにより、特定の時間における沸騰領域の温度及び水率の分布を解析する。この場合における、沸騰領域における模擬的な等価熱伝導度は、例えば、充填ベッドの沸騰領域の開始温度から所定範囲内において最低の熱伝導度であり、所定範囲外の温度において充填ベッドの等価熱伝導度よりも高い熱伝導度である。また、この場合において、解析装置１は、充填ベッドを高さ方向に分割し、解析を行うものとするが、これに限らず、高さ方向、水平方向（幅方向、奥行方向）のうち、一以上の方向に分割して解析を行ってもよい。例えば、解析装置１は、充填ベッドを高さ方向及び幅方向に分割して二次元解析を行ったり、高さ方向、幅方向及び奥行き方向に分割して三次元解析を行ったりしてもよい。
【００４６】
また、装置の分散・統合の具体的な実施の形態は、以上の実施の形態に限られず、その全部又は一部について、任意の単位で機能的又は物理的に分散・統合して構成することができる。また、複数の実施の形態の任意の組み合わせによって生じる新たな実施の形態も、本発明の実施の形態に含まれる。組み合わせによって生じる新たな実施の形態の効果は、もとの実施の形態の効果を合わせ持つ。
【符号の説明】
【００４７】
１・・・高速炉デブリベッド冷却の解析装置、１１・・・入力部、１２・・・表示部、１３・・・通信部、１４・・・記憶部、１５・・・制御部、１５１・・・第１解析部、１５２・・・第２解析部、１５３・・・出力部、１００・・・原子炉容器、１０１・・・コアキャッチャー、１０２・・・デブリベッド
【図１】
2019211340.tif 000007
【図２】
2019211340.tif 000008
【図３】
2019211340.tif 000009
【図４】
2019211340.tif 000010
【図５】
2019211340.tif 000011
【図６】
2019211340.tif 000012
【図７】
2019211340.tif 000013
【図８】
2019211340.tif 000014