炉心溶融挙動を予測する新しい数値シミュレーションコードの開発

デブリの詳細な組成分布の推定に光が見えた

2018/03/23　国立研究開発法人日本原子力研究開発機構

【発表のポイント】

炉心溶融事故における炉内構造物、核燃料の溶融・蓄積、これら核燃料の再臨界の可能性などの予測は難しく、数値シミュレーション手法の開発が切望されている。
原子炉過酷事故時で溶融した炉内構造物などのより実態に近い溶融蓄積挙動を予測するため、数値シミュレーションコード「JUPITER」を新たに開発した。
JUPITERと核計算手法を組み合わせ、原子炉過酷事故時に懸念される再臨界の可能性が極めて低いことを推定した。
今後、JUPITERの計算条件の最適化を図るなどの改良を進め、廃炉作業や過酷事故対策の高度化などを通じて、軽水炉の安全性向上への貢献が期待できる。

国立研究開発法人日本原子力研究開発機構（理事長：児玉敏雄、以下「原子力機構」という。）の山下晋研究員は、今般、炉心が溶融した際に溶け落ちる構造物や核燃料などの物質の動きについて、より実態に近い溶融蓄積挙動を予測することができる数値シミュレーションコード「JUPITER」（ジュピター）を独自に開発しました。

既存手法で必要であった、溶融進展のシナリオの指定や構造物の形状などの簡略化を排除して、材料の熱力学的挙動や流動挙動を含めて溶融物の移動や蓄積の過程を厳密に計算することにより、非常に複雑な分布で堆積する溶融蓄積を可視化しました。また、JUPITERと核計算手法とを組み合わせた計算では、仮想的な初期条件の下では原子炉過酷事故時に懸念される再臨界¹⁾の可能性が極めて小さくなることを推定しました。

今後、JUPITERに対しては、溶融燃料挙動を詳細に検討できるように、大規模な挙動に影響を及ぼす部分的な挙動評価から炉心全体まで拡張した条件でのシミュレーションを可能にするなどの改良を進めます。これにより、東京電力ホールディングス(株)福島第一原子力発電所の廃止措置に向けた燃料デブリの取り出しなどで重要となる溶融物の詳細な堆積状況といった炉内状況把握への貢献も期待できます。

JUPITERに関しては、学術誌「Nuclear Engineering and Design」に掲載されました。また、本研究成果は、「日本原子力学会誌」に掲載される予定です。

【研究開発の背景】

原子力発電所の過酷事故では、冷却材の喪失により高温になった核燃料と周囲の構造物が溶融し、その混合溶融物が原子炉格納容器内に落下、蓄積することが想定されます。しかし、その混合溶融物の核燃料物質（核分裂生成物^2）を含む）と構造物の混合状態や組成分布、混合溶融物からの核分裂生成物の放出、さらに核燃料の再臨界の可能性の予測は難しく、原子炉過酷事故時評価における課題となっています。

また炉心溶融事故では、原子炉停止から全電源停止までの時間差などで生じる炉内温度分布の違いといったわずかな状況の差で、溶融物の流れ方や堆積状況などに違いがあることが想定されます。このため、数値シミュレーションにより、考えられる多様な状況での混合溶融物の組成分布などを予測することが求められています。一方、現在使われている過酷事故の数値シミュレーション手法には、米国スリーマイルアイランド事故³⁾やその後の少数の実験で得られた知識が取り入れられていますが、このような手法では計算時間短縮のため、近似やモデルの簡素化に加え、原子炉内で起こる変化を予め決めておいたシナリオに従ってシミュレーションしており、精度良く炉内状況を推定することが困難でした。

【研究開発の内容】

原子力機構では、より精緻な（近似手法を用いることなく厳密に混合溶融物を取り扱う）数値シミュレーションコード「JUPITER」を開発しました。JUPITERは、決められたシナリオに従ってシミュレーションするのではなく、物理的/熱力学的な基本方程式のみを用いて現象の進行をシミュレーションする点で、従来の予測手法と大きく異なります。JUPITERは、高精度な計算を行うことで数値シミュレーションによる予測に対する信頼性を向上させるのみならず、時空間的に大きな範囲の複雑な現象を数値シミュレーションするため、世界最速クラスのスーパーコンピュータにも適用可能です。

本手法により、従来の手法では推定することのできなかった溶融物の詳細な挙動や蓄積分布を予測することが可能になります。また、この予測結果は過酷事故時の炉心冷却や効率的な廃止措置方法の検討への貢献が期待できます。さらに、福島第一原子力発電所の廃炉作業に対しては、原子炉格納容器内の燃料デブリの組成、核分裂生成物の化学的形態⁴⁾や付着状況、燃料デブリの再臨界の可能性などをより詳細に推定し、安全かつ合理的な廃炉作業の実施に貢献することが期待できます。

【研究成果の概要】

図１に示す体系および条件設定のもとに、JUPITERによる混合溶融物の蓄積数値シミュレーションを行い、３つの項目についてJUPITERの適用可能性を確認しました。

原子炉格納容器下部への混合溶融物の蓄積挙動とその組成
- 図２は、混合溶融物がペデスタル⁵⁾底部へ落下した後の成分と相（固相または液相）分布です。図より、サンプピット⁶⁾内部では、重い成分である溶融燃料成分（黄色）が先に流入して下に蓄積し、その後、軽い成分である構造材分（オレンジ色）が流入して溶融燃料成分の上に積もっている様子が分かります。
- サンプピットの壁面と燃料成分が接している領域は、冷却され凝固している（茶色領域）ことが確認できます。さらに、部分的に凝固することから、混合溶融物は内部で複雑な分布となることが分かりました。
今後は、原子力機構や他の研究機関で行われている実験との比較を行い、より信頼性の高い数値シミュレーションコードに発展させていく予定です。
溶融物中における核分裂生成物の残留推定
- 廃炉作業時には、放射線を出す核分裂生成物の分布やソースターム⁷⁾の情報が重要となります。1.で得られたペデスタル底部の凝固物の温度条件（約2700℃と仮定）から、代表的な放射性核分裂生成物であるセシウム⁸⁾とストロンチウム⁹⁾などが含まれる燃料デブリからの放出挙動を推定しました。その結果、揮発性が高いセシウムは溶融物が凝固する前にほぼ全て気体として放出され、凝固物の中に存在しない一方で、揮発性が低いストロンチウムは凝固物中にある可能性が高いことが分かりました。
- JUPITERでは、空間的な温度分布とその経時変化を得られることから、これまでの数値シミュレーション手法に比べて詳細な核分裂生成物の分布を得ることができます。これは、蓄積した核分裂生成物の安全な取り出し手順を検討する際に必要な知見となります。
再臨界の可能性凝固物の臨界の可能性を検討するために、計算された凝固物の組成分布から、連続エネルギーモンテカルロコードMVP¹⁰⁾を用いて、その再臨界の可能性を調べました。その結果、ペデスタル内部に冷却水がない条件では臨界には至らないことを確認しました。さらに、核燃料成分内に水が存在すると仮定した上で、さらに核燃料を含む物質が空間的に無限に存在する極端に再臨界になりやすい非現実的な条件を設定した場合を除き、再臨界には至らないことが明らかになりました。

図１　数値シミュレーション体系と仮定

図２　ペデスタル底部における溶融物成分分布

【今後の展開】

今後は、より現実的な数値シミュレーションのために、溶融物成分または凝固物成分間での化学反応を考慮できるよう、JUPITERの改良を進めます。さらに、数値シミュレーションに用いる構造物の形状について、より詳細な形状を用いることや、想定される様々なパラメータ（溶融物の流入量、流入速度、温度など）を考慮した数値シミュレーションおよびスーパーコンピュータ「京」や原子力機構の大型計算機を用いた超大規模の数値シミュレーションを進め、想定されるより詳細な移行挙動および分布状況を予測結果として蓄積していく予定です。JUPITERが目指す詳細な凝固物分布の予測が可能になれば、廃炉作業への貢献が期待できます。