シビアアクシデント進展解析コード
MELCORについて
原子力学会核燃料部会「溶融事故における核燃料関連の課題検討」ワーキンググループ
平成24年2月20日
日本原子力研究開発機構 安全研究センター
リスク評価・防災研究グループ
丸山 結
Page 1

シビアアクシデント時の現象とその評価の特徴¾ 現象が多岐に亘り複雑（マルチスケール・マルチフィジィックス）¾ 現象の解明が不十分であること等に起因して、一般的に不確かさが大きい¾ 不確かさが大きいことを考慮して過度に保守的な評価を行うと、非現実的な結果を導く可能性がある¾ 不確かさの低減を図ることが重要で、そのためには安全上重要な現象や不確かさを支配する因子を的確に同定し、効率的に問題を解決することが重要¾ 大規模実験や総合効果実験（システム挙動の評価）では計測等に困難を伴う場合が多い¾ 必要に応じ、理論や機構論的な解析、模擬物質を使った個別効果実験を活用¾ 工学的な判断を要する場合が少なからずある理論
小・中規模
個別効果実験を中心とした実験研究機構論的解析手法の構築

SA進展解析コード（システムコード）整備

3要素の効果的なカップリングと研究資源の適正な配分が重要

PSA
決定論的評価根拠のある工学的判断
etc.

Page 2

シビアアクシデント進展解析コードの役割¾ 個別プラントの設計やAM操作の特徴を反映した
SA事象進展の総合的評価
事象発生のタイミングや放射性物質放出に対する障壁への負荷、ソースタームの評価
不確かさや重要因子の評価¾ SA時個別現象に係る実験や詳細解析コードに対する初期・境界条件の提供
z AMやSA対応設計の有効性評価
z 格納容器内ソースタームの更新（立地評価）
z レベル2PSA （事故進展解析、ソースターム評価）
z 実効性の高い防災計画の立案
z 研究ニーズや優先順位の検討
Page 3

シビアアクシデント進展解析コードの位置付け
個別現象に係る実験• 実験ニーズ• 優先順位• 実験条件• 検証• 機構論的モデル構築

個別現象詳細解析コード• 詳細解析コード
開発のニーズ• 優先順位• 解析条件• 技術情報（簡易化モデルの組み込み等）•検証•モデル構築SA進展解析コード（MELCOR・THALES・MAAP等）

Page 4

シビアアクシデント進展解析コードへの要求¾シビアアクシデント時に想定される現象を網羅していること¾個別モデルやシステム挙動の十分な検証がなされていること¾最新の知見が反映されていること¾それが容易にできること¾数値解析の頑健性が高いこと¾（ユーザーが多いこと）¾…
Page 5

シビアアクシデント時の諸現象（PWR）
冷却機能喪失失敗

格納容器

Phase-I AM（SAの発生防止）

SAの発生
失敗

爆発的な現象による短時間のエネルギー発生
A 早期格納容器破損

Phase-II AM（SAの影響緩和）

崩壊熱やZr酸化による継続的なエネルギー発生
B 後期格納容器破損

圧力容器破損格納容器リーク・破損
格納容器バイパス
環境へのFP放出

環境へのFP放出

燃料からのFP放出
原子炉冷却系内FP移行
格納容器内FP挙動

格納容器
リーク・破損格納容器
バイパス 環境へのFP放出

圧力容器
圧力容器内炉心溶融進展
圧力容器内
溶融炉心冷却

高圧融体放出格納容器直接加熱（DCH）

水素燃焼／爆燃／爆轟
圧力容器破損

溶融炉心−冷却材相互作用（FCI）（水蒸気爆発含む）
圧力容器外溶融炉心冷却

溶融炉心−コンクリート相互作用（MCCI）
Page 6

シビアアクシデント時の諸現象（BWR）
フランジ部破損
環境へのFP放出

燃料からのFP放出・
原子炉冷却系内FP移行

格納容器

圧力容器

水素燃焼／爆燃／爆轟

格納容器内FP挙動
炉容器内炉心溶融進展
環境へのFP放出

炉容器内溶融炉心冷却

貫通部破損

炉容器破損

高圧融体放出格納容器直接加熱(DCH)

溶融炉心／冷却材相互作用 (水蒸気爆発を含む)
FPのプール
スクラビングによる除去

炉容器外溶融炉心冷却

溶融炉心/コンクリート相互作用(MCCI)Page 7放射性物質の移行に係わる諸現象
粒子の成長(凝集)

層流／乱流
再浮遊

粒子状FPの沈着挙動

拡散

重力沈降
気液平衡

スクラビング

気液平衡気泡捕集

凝縮
気液平衡
拡散泳動

スプレイ
液滴

さえぎり

拡散FP
放出源

フィルタ
除去

環境

慣性衝突

スプレイによるFP除去

重力沈降

気体

再蒸発

乱流

FP
放出源

熱泳動
拡散泳動

液体

化学吸収

拡散

拡散

S/CFPの沈着挙動

漏洩

環境

溶解

Page 8

MELCORコードの概要¾ 米国原子力規制委員会／サンディア国立研究所
において開発されているシビアアクシデント進展解
析コード¾ 炉心溶融進展やソースタームに係わる主要な現象を網羅¾ 原子炉の主要な機器や系統（ECCS、格納容器スプレイ等）を考慮可能¾ ユーザーが多いSA解析コードの一つ¾ 解析コードの開発にユーザーの使用経験（検証を含む）を反映

Page 9MELCORコードの現象に係わるパッケージ群¾¾¾¾¾¾¾¾¾¾¾¾¾¾¾

BH：BWR Bottom Head→下部プレナム・下部ヘッドの熱応答等
BUR：Burn→水素・一酸化炭素の燃焼・爆燃
CAV：Cavity→溶融炉心／コンクリート相互作用（MCCI）
CND：Condenser→アイソレーションコンデンサー（IC）・受動的格納容器冷却系（PCCS）
COR：Core→炉心・炉内構造物の熱応答・リロケーション等
CVH/CVT/FL ：Control Volume Hydrodynamics/Control Volume
Thermodynamics/Flow Path →コントロールボリュームの熱水力・熱流動等
DCH：Decay Heat→崩壊熱
FCL：Fan Cooler→ファンクーラー内質量・熱移行
FDI：Fuel Dispersal Interactions→下部ヘッド破損時の低圧・高圧融体放出
HS：Heat Structure→熱構造体の熱伝導
MP：Material Properties→物性
NCG/H2O ：Non-Condensable Gas/Water →非凝縮性気体・蒸気表
PAR：Passive Autocatalytic Hydrogen Recombiner→触媒式水素再結合器
RN：Radionuclide→放射性物質の放出・移行挙動
SPR：Containment Spray→スプレイ液滴／格納容器雰囲気間の熱伝達

Page 10MPパッケージ¾ 組み込まれている燃料関連物質の種類
UO2、U、ジルカロイ、ZrO2、ステンレス鋼、ステンレス鋼の酸化物、制御材（B4C及びAg-In-Cd）¾ 組み込まれている物性*比エンタルピ、比熱、密度、融点、溶融潜熱¾ その他の物質について入力データで定義可能*主要参考文献(1) D. L. Hagrman, G. A. Reymann, and R. E. Mason, MATPRO VERSION 11 (Revision 1) A
Handbook of Materials Properties for Use in the Analysis of Light Water Reactor Fuel Rod
Behavior, NUREG/CR-0497 and TREE-1280 Rev. 1, EG&G Idaho, Inc.,Idaho Falls, ID (February1980).
(2) D. L. Hagrman, Materials Properties for Severe Core Damage Analysis, EGG-CDD-5801, EG&G,
Idaho Falls, ID (May 1982).
(3) F. P. Incropera and D. P. Dewitt, Fundamentals of Heat Transfer, John Wiley and Sons, New York(1981).
(4) M. M. Wakil, Nuclear Heat Transport, International Textbook Co., New York (1971).
(5) CRC Handbook of Chemistry and Physics, 63rd ed., CRC Press, Boca Raton, FL (1982).
Page 11CORパッケージにおけるノード分割¾ 炉心及び下部プレナムを軸対称の複数セルに分割¾ 炉心に関連する主な入力パラメータ セルの径及び高さ
 径及び軸方向の出力分布
 セル内構成物質の量
 燃料ペレット及び被覆管の半径
 ペレットと被覆管間のギャップ幅
 燃料棒のピッチ
 炉型及び制御材の種類
 チャンネルボックスの肉厚（BWR）
Page 12

酸化モデル¾ ジルカロイ、ステンレス鋼、B4Cの酸化を考慮
反応速度：
ステンレス鋼／水蒸気[3]Zr／水蒸気[1]B4C／水蒸気（簡易）[4]Zr／酸素[2] 水蒸気あるいは酸素が極めて少ない場合における気相内の物質伝達による酸化速度の制限を考慮（熱伝達とのアナロジー）
 熱化学平衡（Gibbs自由エネルギー最小化手法）に基づいたB4C酸化モデルも利用可能[1] V. F. Urbanic and T. R. Heidrich, "High-Temperature Oxidation of Zircaloy-2 and Zircaloy-4 in Steam," J. Nuc. Matls., 75,
pp. 251-261 (1978).
[2] A. S. Benjamin, D. J. McCloskey, D.A. Powers, and S. A. Dupree, Spent Fuel Heatup Following Loss of Water During
Storage, SAND77-1371, NUREG/CR-0649, Sandia National Laboratories, Albuquerque, NM, March 1979.
[3] J. F. White et al., "Fifth Annual Report--High Temperature Material Programs, Part A," GEMP-400A (February 1966).
[4] L. J. Ott, C. F. Weber, and C. R. Hyman, "Station Blackout Calculations for Browns Ferry," Proceedings of the Thirteenth
Water Reactor Safety Research Information Meeting, Gaithersburg, MD (October 1985).
Page 13

物質相互作用モデル¾ 共晶反応 ジルカロイ／インコネル：1400K[1] ジルカロイ／ステンレス鋼：1400K[1] B4C／ステンレス鋼：1520K[2]¾ 溶解（Dissolution）
溶解速度[1][3]

溶融ジルカロイ中へのUO2及びZrO2の溶解に係わる実験をベース[1] P. Hofmann et al., “Reactor Core Materials Interactions at Very High Temperatures,” Nuclear Technology, 87, pp. 146-186, August 1989.
[2] W. Hering and K. Muller, “Modelling of Eutectic Interactions in KESS-III (Module EUTECT),” International CORA Workshop 1992,
Karlsruhe, FRG, October 5-8, 1992.
[3] G. Eriksson, "Thermodynamic Studies of High Temperature Equilibria: XII. SOLGASMIX, a computer program for calculation ofequilibrium compositions in multiphase systems," Chemica Scripta, 8, 1975, pp. 100-103.

Page 14

RNパッケージ
放射性物質のグループ化

Page 15

RNパッケージの放射性物質放出モデル¾ 燃料からの放射性物質放出に関して以下の3つのモデルを導入
CORSOR[1]f：fraction/min
CORSOR-M[1]CORSOR-Booth[2]→燃料内拡散と気相内物質伝達（抵抗）を考慮

燃料から放出された後の化学反応（放射性物質／制御材、放射性物質間等）は限定的
TMI-2事故から得られた知見等に基づいて、燃料から放出された後のIはCsI、残りのCsはCsOHの化学形を取ると仮定するのが一般的[1] M. R. Kuhlman, D. J. Lehmicke, and R. O. Meyer, CORSOR User’s Manual, BMI-2122, NUREG/CR-4173 (March 1985).
[2] M. Ramamurthi and M. R. Kuhlman, Final Report on Refinement of CORSOR—An Empirical In-Vessel Fission Product Release Model,
Battelle Memorial Institute (October 31, 1990).

Page 16まとめ¾ MELCORコードは、燃料の溶融や燃料からの放射性物質の放出を含めたシビアアクシデント時の主要な現象を網羅¾ 放射性物質の化学的挙動（放出時や放出後の化学形）について検討の余地あり¾ 福島第一原子力発電所のシビアアクシデント解析に適用（米国、JNES、JAEA等）

Page 17