本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

改訂番号

改訂 28

提出年月日

平成 29 年 4 月 17 日

重大事故等対策の有効性評価に係るシビアアクシデント解析コードについて

平成 29 年 4 月
東北電力株式会社
東京電力ホールディングス株式会社
中部電力株式会社
中国電力株式会社

目次(Rev.08)
１．はじめに・・・・・・・・・・・・・・・・・・・・・・・・・・・1-1
２．有効性評価における物理現象の抽出・・・・・・・・・・・・・・・2-1
２．１ 炉心損傷防止・・・・・・・・・・・・・・・・・・・・・・2-7
２．２ 格納容器破損防止・・・・・・・・・・・・・・・・・・・・2-49
２．３ 運転停止中原子炉における燃料損傷防止・・・・・・・・・・2-70
３．抽出された物理現象の確認・・・・・・・・・・・・・・・・・・・3-1
３．１ BWR プラントシステムの階層構造分析と抽出された物理現象の対応確認・3-1
３．２ EURSAFE における物理現象と抽出された物理現象の対応確認・・・・ 3-2
４．適用候補とするコードについて・・・・・・・・・・・・・・・・・4-1
４．１ 適用候補コードの概要・・・・・・・・・・・・・・・・・・4-1
５．有効性評価に適用するコードの選定・・・・・・・・・・・・・・・5-1
５．１ 炉心損傷防止・・・・・・・・・・・・・・・・・・・・・・5-1
５．２ 格納容器破損防止・・・・・・・・・・・・・・・・・・・・5-5
５．３ 運転停止中原子炉における燃料損傷防止・・・・・・・・・・5-6
６．選定されたコードの有効性評価への適用性について・・・・・・・・6-1
添付資料１ 許認可解析と重大事故等対策の有効性評価の比較について
第１部 SAFER(Rev.15)
第２部 CHASTE (Rev.09)
第３部 REDY (Rev.14)
第４部 SCAT (Rev.16)
第５部 MAAP
本文 (Rev.13)
添付１ 高圧溶融物放出／格納容器雰囲気直接加熱の防止について (Rev.10)
添付２ 溶融炉心と冷却材の相互作用について (Rev.11)
添付３ 溶融炉心とコンクリートの相互作用について (Rev.15)
第６部 APEX(Rev.08)

今回ご提出箇所

Rev.15添付３

溶融炉心とコンクリートの相互作用について目次
1. まえがき ............................................................ 5-3-1
2. 現象の概要 .......................................................... 5-3-2
3. これまでの知見の整理................................................. 5-3-3
4. 評価モデルと不確かさに関する整理 ..................................... 5-3-6
5. 感度解析と評価 ..................................................... 5-3-29
6. まとめ ............................................................. 5-3-41
付録 .................................................................. 5-3-42
付録 1 ドライ条件での MCCI 実験 ....................................... 5-3-42
付録 2 注水を伴った MCCI 実験 ......................................... 5-3-49
付録 3 溶融物の水中落下時の粒子化実験と解析モデル .................... 5-3-74
付録 4 溶融物の拡がり実験 ............................................ 5-3-91
付録 5 粒子状ベッドの冷却性実験と解析モデル .......................... 5-3-105
参考文献 .............................................................. 5-3-107

１．まえがき
溶融炉心－コンクリート相互作用（MCCI）は，炉心溶融が進展し，溶融炉心が原子炉圧
力容器底部を溶融貫通し，格納容器下部のペデスタル床に落下した場合に発生する現象である。MCCI が継続すると，格納容器構造の侵食や非凝縮性ガス（水素等）発生による格納
容器過圧により格納容器破損に至る可能性があり，重要な格納容器破損モードと考えられてきており，これまで種々の実験や解析モデルの開発が行われてきている。これまでは，注水が無いドライ条件での MCCI 挙動や溶融炉心上に注水した場合（Top Flooding）の MCCI抑制効果を調べる研究が主に行われてきた。
MCCI の緩和対策として，溶融炉心落下後に注水する以外に，溶融炉心落下前にペデスタ
ルに事前水張りするマネジメントが効果的とされており[1]，本有効性評価ではこの対策を前提としている。事前水張りにより，溶融炉心が落下時に水中で粒子化され，溶融炉心が粒子状ベッドとして堆積し，デブリ冷却性の向上が期待されるためである(図 1-1 参照)。
本資料では，国内外で実施された実験等による知見を整理するとともに，解析モデルに関する不確かさの整理を行い，感度解析により有効性評価への影響を確認した。

原子炉圧力容器

冷却水

溶融炉心

粒子化

水

粒子状
ベッド

図1-1 ペデスタル事前水張りによるデブリ粒子化の概念

5-3-1

1

２．現象の概要
MCCI は，溶融炉心が原子炉圧力容器下部を溶融貫通して格納容器床面に落下し，床面のコンクリートと接触した場合に発生し，高温の溶融炉心によりコンクリートが加熱され熱分解し溶融侵食される現象である。侵食が継続した場合，格納容器バウンダリの破損や格納容器内支持構造の強度低下に至る可能性がある。また，コンクリートの熱分解により発生する水蒸気及び二酸化炭素が溶融炉心内を通過する間に未酸化金属成分と反応して水素や一酸化炭素等の非凝縮性の可燃性ガスが発生し，格納容器内の過圧や燃焼に至る可能性がある。このような要因により格納容器破損に至る可能性があり，MCCI は重要な格納容器
破損モードと考えられている。MCCI を停止するためには，落下した溶融物（コリウム）を冷却することにより，侵食限界に至る前にコリウムと接しているコンクリート温度を溶融
温度(約 1500 K)以下にすることが必要となる。
国内 BWR では，原子炉圧力容器破損前にペデスタルへ水張りを行い，下部ヘッドが破損
し溶融炉心が落下した際の溶融炉心の冷却を促進することにより，MCCI を緩和する対策を採っている。ペデスタルに落下した溶融炉心は，水プール中を落下する際に，一部は粒子化して水中にエントレインされ，残りは床面に落下して堆積し溶融プールを形成する。エ
ントレインされたデブリ粒子は，水中を浮遊しながら水へ膜沸騰による熱伝達が行われて冷却され，クエンチし，溶融プール上に堆積し，粒子状ベッドを形成する。
ペデスタル床に堆積した溶融炉心は，崩壊熱や化学反応熱により発熱しているが，水やコンクリートへの伝熱により徐々に冷却され，溶融炉心温度が固化温度を下回ると固化する。溶融炉心の冷却の過程では，中心に溶融プール(液相)，外周部にクラスト(固相)が形成される。溶融炉心上のプール水との伝熱は粒子状ベッドを介して行われるため，冷却効
果は，粒子状ベッドの冷却性に依存する。
ペデスタルの溶融炉心と接触しているコンクリートは，溶融炉心からの熱伝達により加熱され，その温度が融点を上回ると溶融し侵食される。この際にガス(水蒸気及び二酸化炭素)及びコンクリートスラグが発生し，溶融炉心に混入されて化学反応（未酸化 Zr による還元反応）により，反応熱とともに水素や一酸化炭素が発生する。
MCCI によるコンクリート侵食量は，溶融炉心から上部水プール，壁方向及び床方向のコ
ンクリートへの熱分配に支配される。国内 BWR の MCCI 対策としている事前水張りをしている場合には，溶融炉心上面からの除熱が支配的となるため，相対的にコンクリートへの熱移行は抑えられ，コンクリート侵食量は抑えられる。この場合には，上面熱流束が重要となる。また，万一事前水張りできず，溶融炉心上面からの除熱が十分でない場合は，コン
クリートへの熱移行量が増え，コンクリート侵食量は大きくなる。この場合には，コンクリートへの熱伝達が重要となる。

5-3-2

2

３．これまでの知見の整理
本章では，MCCI に関する実験の概要及びそこで得られた知見に関して整理を行う。
MCCI に関する実験としては，水による冷却を伴わない実験として米国アルゴンヌ国立研究所（ANL）で実施された ACE 実験[2][3]，米国サンディア国立研究所（SNL）で実施された
SURC 実験[4]（国際標準問題 ISP-24[5]）及び独カールスルーエ研究所（KfK）で実施された
BETA 実験[6]がある（付録１参照）。水による冷却を伴う実験（溶融物上に冷却水を注水した実験）としては，SNL で実施された SWISS 実験[7]及び WETCOR 実験[8]，米国電力研究所（EPRI）の主催で実施された MACE 実験[3][9][10]，我が国の原子力発電技術機構（NUPEC）により実施された COTELS 実験[1][11]，OECD プロジェクトとして ANL にて行われた MCCI 実験[10][12][13][41]がある（付録２参照）。
また，水中に溶融デブリを落下させ溶融物の粒子化を調べた実験としては，JRC-ISPRA にて実施した KROTOS 実験[14]や FARO 実験[14]，NUPEC にて実施した COTELS-FCI 実験[1][15]，ス
ウェーデン王立工科大学（KTH）で実施された DEFOR 実験[16]がある（付録３参照）。
ペデスタルでの溶融物の拡がり実験としては，水による冷却を伴わないドライ条件での実験として，国内 BWR を対象に BWR 産業界が実施した SPREAD 実験[19][20]や，EPR を対象とした実験が複数実施されており，ウェット条件での実験については実施例が少ないが，前述の SPREAD 実験や KTH の PULiMS 実験[21]，ANL にて不均質に堆積させたデブリベッドの拡がりを確認したセルフレベリング実験[22]が行われている。
各実験の詳細は付録資料にて示す。以下に，実機評価に関連する知見の概要を纏める。
・ドライ条件でのコンクリート侵食 （付録１参照）
溶融炉心に注水が行われない場合の侵食速度は，ACE 実験，SURC 実験及び BETA 実験で確認されている。また，MAAP コードによる実験解析が行われており，解析で得られた侵
食速度は実験結果と良く一致している（1 時間あたり 20cm 程度）[23][34]。
・ウェット条件（Top Flooding）でのコンクリート侵食（付録２参照）

SWISS 実験，WETCOR 実験，MACE 実験では，溶融デブリ上面へ注水をした場合（TopFlooding）のコンクリート侵食実験が行われたが，結果として，デブリ上面に安定なハードクラストが形成されてデブリ内への水の浸入を妨げ，コンクリート侵食が継続する結果が得られた。コンクリート侵食が停止しなかった主な原因は，デブリ模擬物上面に形成されたクラストが側壁（側壁に耐火物を採用した 1 次元侵食実験）あるいは電極と側壁の両者に固着し，クラストとその下のデブリ模擬物とが分離したことにより，デブリ
模擬物の効果的な除熱がなされなかったためであるとされている。
COTELS 実験では，溶融物の落下過程を含めて模擬した MCCI 実験が実施され，結果として，側壁侵食部に水が浸入し，コンクリート侵食が停止する結果が得られた。
実機スケールで安定クラストが形成されるか否かを解明するため，OECD/MCCI プロジェ

5-3-3

3

クトが行われ，クラスト強度，クラスト浸水，2 次元コンクリート侵食挙動等を分離効果的に調べる実験が実施された。その結果，実機スケールでは，安定クラストは形成されず，ハードクラストは割れて水が内部に浸入し，MCCI により発生したガスにより溶融物がクラストの割れ目から噴出する火山型のクラストが形成されるとの結果が得られている。また，CCI 実験結果から玄武岩系コンクリートに侵食の異方性が見られる結果が得られている。
これらの知見は，MCCI 評価のためにはデブリから水プールやコンクリートへの熱伝達の扱いが重要であることを示唆している。
・水張りによる溶融物の粒子化 （付録３参照）
溶融物ジェットが水中へ落下する場合の粒子化挙動については，FCI 実験等により確認されている。粒子化割合は，主に水深やジェット径に依存し，粒径（質量中央径）は，比較的大きく，実験条件（初期圧力，水深，コリウム落下速度，サブクール度）に対する依存性は低いと報告されている。各 UO2 混合物実験の平均的な粒子径は，溶融物量の多い FARO 実験では 2.6～4.8mm[14]，COTELS 実験では 6mm 程度[15]である。また，粒子化割合を評価する解析モデルが提案されている（Ricou-Spalding 相関式[25]や Saito の相関式[26]等）。
・溶融物の拡がり （付録４参照）
Mark-I 型格納容器のシェルアタックに関する研究[17][18][19][20]により，事前水張りが，溶融物の拡がりを抑制し，シェルアタック防止に有効であることが確認されるとともに，溶融物の拡がりを評価する解析コードが作成されている。溶融物の拡がりには，デブリ
落下流量，デブリ成分，デブリ過熱度が主に影響することが示されている。KTH では，水中での溶融物の拡がり挙動と固化性状を調べる PULiMS 実験[21]が実施され，溶融物の拡がる過程は，流体力学と固化の２つの競合プロセスに支配されるとして，流体力学には重力，慣性力，粘性力及び表面張力が影響し，固化には溶融物から周囲への熱伝達，崩壊
熱発生及び溶融物の相変化が影響するとして，実験と実機条件とを関連づけるスケーリ
ング則の提案[20][27][28]が行われている。また，独カールスルーエ研究センター（FZK）に
おいて KATS 実験[23][29][30]が実施され，溶融物の放出速度が比較的高い場合は，冷却材の有無によらず同様な拡がり挙動になることが示されている。溶融物の拡がり挙動の解析のために，種々の解析コードが提案されており，実験データを元に検証が行われている。
3 次元の拡がり評価も行われており，比較的広いペデスタルをもつ ABWR（ペデスタル半
径約 5.3m）の溶融物拡がり評価が実施され，床上に水がある場合でも，床全面に溶融物が拡がることが示されている[40]。
また，ANL でのセルフレベリング実験[22]により，水プール中に不均質に堆積させた粒子状ベッドにおいて，内部沸騰によって短時間でベッドの厚さが均一化されることが示さ

5-3-4

4

れている。
・粒子状ベッドの冷却性 （付録５参照）
粒子状ベッドの冷却性については，ドライアウト熱流束が種々の実験で確認されており，主に粒子径に依存するとの結果が得られている。また，ドライアウト熱流束を評価する解析モデルがいくつか提案されているが，Lipinski0-D モデル[24][31][32][33]が広く使わ
れている。

5-3-5

5

４．評価モデルと不確かさに関する整理本章では，MCCI 評価に使用している解析モデルと想定される不確かさについて整理する。
4.1

MCCI 評価モデル

MCCI 評価には，MCCI 現象を扱え，溶融炉心落下や注水によるデブリ冷却効果を評価でき
る MAAP コード[34]を使用している。以下に，MCCI 伝熱モデルの概要を示す（図 4-1 参照）。
コリウムの構成・ 均質に溶融したプールを形成すると仮定。
・ 上部，下部及び側面にクラスト層を仮定。
クラストの生成・ クラストのエネルギーバランス（溶融プールからの伝熱，水プール・コンクリートへの伝熱）によりクラスト厚さの変化率を計算。
コンクリートの侵食
・ コンクリート表面から侵食深さ方向に１次元熱伝導を解き温度分布を計算。
・ コンクリート溶融温度（入力値：1500K）以上で侵食開始。
・ 溶融炉心からの伝熱量と分解／溶融潜熱により侵食量を評価。
・ コンクリート分解による自由水・結合水・CO2の発生を考慮し，溶融コリウム中での化学熱力学平衡計算により金属との反応，H2及びCO発生を評価。反応熱は，溶融
炉心の崩壊熱に加算される。
伝熱モデル
・ 溶融炉心－コンクリート間の熱伝達； 溶融プールからクラスト層への対流熱伝達量とクラスト内崩壊熱の和で熱流束が与えられる。溶融プールからクラスト層への対流熱伝達は，粘性に影響されるため，固化割合に応じて補正されている。溶融炉
心から床方向及び側面方向のコンクリートへの熱流束を以下に示す。
=ℎ

−

ℎ =ℎ

1−

=ℎ

−

ℎ =ℎ

1−

床方向の熱流束
壁方向の熱流束

,

,

+
+

∙
∙

ここで，
， ： 床方向及び側面方向の熱流束[W/m2]ℎ ，ℎ ： 溶融プールからクラスト層への対流熱伝達係数[W/m2-K]ℎ ，ℎ ： 溶融プールが完全な液相の場合の対流熱伝達係数： 固化割合[－]： 固化効果項の指数[－]： 溶融プールの温度

[K]

5-3-6

6

[W/m2-K]

本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

,： デブリ融点： 体積発熱率，[K][W/m3]： 床面及び壁面のクラスト厚さ[m]

対流熱伝達係数ℎ 及びℎ は，BETA実験（2次元侵食実験）のベンチマーク解析[23][34]を元に，それぞれとが，固化効果項の指数 は，が推奨されている。
・ クラストの伝熱と厚さ； クラストのエネルギーバランス（溶融プールからの伝熱，水プール・コンクリートへの伝熱）によりクラスト厚さの変化率が計算されている。
側面方向及び床方向の溶融プールからの対流熱伝達係数が異なるため,下部クラス
トと側面クラストの厚さも異なる。内部発熱を考慮し，下部クラスト,側面クラスト,
上部クラストは放物線状の温度分布が仮定されている。クラスト成長の結果，溶融プールからの熱とクラスト内発生熱をコンクリートへ伝達できなくなると，クラス
ト内側で溶融が生じ，クラストは薄くなる。逆に,クラストが十分に薄く,溶融プー
ルからの熱とクラスト内発生熱以上を伝達できれば，凝固が生じ，クラストは厚くなる。
・ 上部クラスト－水プール間の熱伝達；

デブリが冠水されて表面温度が十分低下す

ると,上部クラストはクエンチ層として扱われ，以下の Kutateladze型の水平平板限界熱流束相関式が適用される。このクエンチ熱流束は,デブリ平均温度が水の飽和温度以上である場合に適用される。

q chf = Fchf h fg ρ g0.5 {σg (ρ l − ρ g )}

0.25

ここで，

q chf ：熱流束 [W/m2]
Fchf ：係数 [－]（ユーザー入力）
h fg ：蒸発潜熱 [J/kg]

ρ g ：蒸気の密度 [kg/m3]

ρ l ：水の密度 [kg/m3]
σ ：表面張力 [N/m]g：重力加速度[m/s2]

係数 Fchf として，SNLで実施されたSWISS実験において報告されている溶融物から水プールへの熱流束800kW/m2に相当する値である0.1がデフォルトとして設定されており，係数 Fchf を変更することにより粒子状ベッドから水プールへの熱流束を想定した解析が実施できる。係数 Fchf と熱流束 q chf の関係を図4-2に示す。
デブリ平均温度が水の飽和温度と等しくなると(デブリクエンチ)，水プールへの熱流束は，デブリの崩壊熱が全て水プールへ伝達されるとして計算され，上記限界

5-3-7

7

熱流束相関式は適用されない。デブリクエンチ 後は，デブリからコンクリートへのデブリからコンクリートへのしているコンクリート温度が侵食温度以
熱伝達は発生しなくなるため，デブリと 接しているコンクリート温度下となり，コンクリートの侵食は 停止する。

溶融プール

プール温度
プール融点

クラスト

(放物分布)

境界温度

コンクリート

図 4-1

MAAP コードの MCCI 伝熱モデル概要

コードのデブリ上面
上面の熱流束を与える
図 4-2 MAAP コードのデブリ
Fchf 入力 と限界熱流束の関係

5-3-8

8

4.2 MCCI評価における不確かさの整理
MCCI の過程毎に不確かさ要因を整理する。MCCI は以下の過程で段階的に進展する。
① 原子炉圧力容器下部ヘッド破損過程② 溶融物の落下・堆積過程③ MCCI 進行と注水によるデブリ冷却過程
各過程での物理現象及び解析モデルに関し，不確かさ要因を整理し，感度解析が必要なパラメータを確認する。
(1) 原子炉圧力容器下部ヘッド破損過程
この過程の不確かさは，下部ヘッドの破損面積や溶融炉心の流出量となる。それぞれについて，評価モデルと評価条件の考え方を示す。
下部ヘッド破損面積評価モデル： MAAP コードでは，制御棒駆動機構ハウジング溶接部の破損や下部ヘッド
クリープ破損等の種々の下部ヘッド破損モードがモデル化されており，最も早く破損条件に至った破損モードにより破損が生じる。破損面積は，溶融物流出に伴う破損口の径方向
侵食による拡大も考慮されている。
評価条件： MAAP 解析結果から，制御棒駆動機構ハウジング溶接部の破損による制御棒
駆動機構ハウジングの逸出が最も早く発生する。これにより約 15cm 径の破断口が下部ヘッ
ドに開き，下部プレナム内の溶融炉心がペデスタルに落下することになる。他の下部ヘッド貫通部が破損してもこれ以上の破損サイズにはならないことから，制御棒駆動機構ハウ
ジングの逸出が最も厳しい条件となっている。なお，制御棒駆動機構ハウジング溶接部が破損しても外部サポートが設置されている場合には，制御棒駆動機構ハウジングの逸出は起きにくいと考えられるため，この扱いは保守的である。
溶融炉心の落下量
評価モデル：

溶融炉心流出流量は駆動差圧として原子炉圧力容器内とペデスタルとの圧力差に水及び溶融炉心の水頭圧を考慮して計算されている。
評価条件： MAAP 解析結果から，下部ヘッド破損前には全炉心相当の溶融炉心が下部プ
レナム内に堆積した状態になっている。この状態に至る前には，以下の過程が起こる。即ち，溶融炉心が下部プレナムへ落下を開始すると，下部プレナム水により溶融炉心が冷却される。また，制御棒案内管内にも水が溜まっており，溶融炉心の冷却のヒートシンクになる（図 4-3 参照）。下部プレナムや制御棒案内管内の水が蒸発して喪失すると，制御棒案
内管の温度が上昇し溶融に至る。制御棒案内管は炉心重量を支持しているために，制御棒
案内管の溶融が起こると，炉心の支持機能が喪失し，全炉心が下部プレナムへ落下する。
その後，制御棒駆動機構ハウジング溶接部の破損が起こり，溶融炉心の放出が開始する。

5-3-9

9

以上のような過程の結果，下部ヘッド破損時に，下部プレナム内に堆積している全炉心相
当の溶融炉心が流出することになり，溶融炉心の落下流量としては，最も厳しい条件となっている。
以上に示したように，下部ヘッドの破損面積や溶融炉心の流出量について，いずれも最も厳しい条件が使用されているため，特に感度解析を行う必要はない。

図 4-3 炉心部から下部プレナム内への溶融炉心の移行状態を示す模式図5-3-1010本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

(2) 溶融物の落下・堆積過程
この過程では，下部ヘッドが破損し，溶融炉心が事前水張りされたペデスタルへ落下す
るため，溶融炉心が水中で粒子化することが想定され，粒子化されない溶融炉心はペデス
タル床に溶融プールとして堆積する。この過程では，溶融炉心の水中での粒子化に不確かさが想定される。
溶融炉心の粒子化評価モデル： MAAP コードでは，粒子化割合は，Ricou-Spalding モデル[25]で計算される。
Ricou-Spalding モデルは，エントレインメント量（粒子化量）を流入流体の速度（ジェット速度）と両流体の密度比に関連させたモデルであり，液液混合問題において広く利用されている相関式である。落下する溶融炉心は円柱ジェットで冷却水中に突入することを想定し，円柱ジェット外周部の侵食として粒子化割合を評価している。
円柱ジェットのエントレインメント速度 ment（m/s）は，次式で与えている。
1ρ  2
ment = Eo  w  udjρ  dj 

(1)

ここで， E o ：エントレインメント係数

u dj ：ジェット速度 [m/s]ρdj ：ジェット密度 [kg/m3]
ρ w ：水密度 [kg/m3]
水中でのジェット速度を一定と仮定し，水面から垂直方向に積分すると，プール底部（水深： ∆H pool ）におけるジェット直径が(2)式のように得られる。
1ρ  2d dj = d dj,o − 2E o  w  ∆H pool ρ dj (2)(2)式を用いて，溶融ジェット断面の減少分が粒子化割合であるため，次式のように計算される。

d dj,o − d dj2

Φ ent =

d dj,o2(3)2

ここで， ddj ：プール底部におけるジェット径 [m]d dj , 0 ：プール水面におけるジェット径 [m]Φ ent ：水中における粒子化割合 [-]以上の評価式において，デブリ粒子化割合に影響するパラメータは，エントレインメン
ト係数，溶融ジェット径及びプール水深であり，このうちモデルパラメータとしてエント
レインメント係数と溶融物ジェット径に不確かさが想定される。
エントレインメント係数については，MAAP では代表的な FCI の大規模実験である FARO 実
験に対するベンチマーク解析[34]によって範囲を設定しており，不確かさ範囲は，なっている。FARO 実験条件

[14]

～

は，水プールの水深は 0.87～2.05m，水プールのサブ

5-3-11

11

本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

クール度は 0（飽和）～124K，雰囲気圧力条件は 2～5.8MPa（高圧条件），0.2～0.5MPa（低圧条件）の範囲で行われており，一方，実機条件では，水深は 2～4m 程度，水プールサブ
クール度は飽和～90K 程度，雰囲気圧力は 0.2MPa[abs]程度であることを考慮すると，実験条件は有効性評価の特徴的な条件を包絡しており，ベンチマーク解析で検討されたエント
レインメント係数の範囲内で感度を確認すれば十分といえる。
評価条件： MAAP 解析ではエントレインメント係数のデフォルト値としてが使用されており，これは実験解析により確認された不確かさ範囲の代表値となっている。
溶融物ジェット径は，下部ヘッド破損口径により決まることから，前述のように大きい
値としている。これは，粒子化割合を小さめにする扱いをしていることになる。
因みに，Ricou-Spalding 相関式を使用して，実機で想定される水張り水深 2m, 溶融物
ジェット径 0.15m の条件の場合で 63%程度が粒子化される（図 4-4 参照）。[33]以上のことから，エントレインメント係数については，不確かさ範囲が設定されており，溶融炉心が水中に落下した際の粒子化割合とその過程での蒸気発生やデブリ粒子の酸化による水素発生に影響することから，不確かさ範囲内で感度解析を実施する。

図 4-4

Ricou-Spalding 相関式による粒子化割合のマップ[33]5-3-1212

(3) MCCI 進行と注水によるデブリ冷却過程
この過程では，溶融炉心の拡がり，デブリベッドの冷却性，コンクリート侵食に不確かさが想定される。
溶融炉心の拡がり
評価モデル： MAAP コードでは，落下した溶融炉心はペデスタル床上に均一に拡がると仮定している。
評価条件：

事前水張りされたペデスタルのプール水中に落下する溶融炉心は，前節の溶融物の落下・堆積過程で示したように，一部が粒子化した後に固化し，残りが溶融状態のまま床面に到達して床面上を拡がり，固化したデブリ粒子が床上の連続層の上に堆積して粒子状ベッドを形成するものと想定される（図 1-1 参照）。このような溶融炉心の拡がりにおいて，溶融物の拡がり距離と粒子状ベッドの堆積形状に不確かさが想定される。
まず，溶融物の拡がり距離の不確かさについて考察する。Mark-Ⅰ型格納容器のシェルア
タックに関する研究[17][18][19][20]によれば，水張りが無いと溶融炉心は格納容器シェル壁（拡がり距離 7m 程度）まで到達するが，水張りがあれば溶融炉心の拡がりは抑えられ，拡がり
距離は落下量等にも依存するが 5m 程度となる結果が得られている。有効性評価で想定している制御棒駆動機構ハウジングの逸出を想定すると，ABWR（ペデスタル半径約 5.3m）で床上に水がある場合でも，落下位置に依らず床全面に溶融物が拡がることが，3 次元溶融物拡がり解析により示されている[40]（付録４(5)a）参照）。また，PULiMS 実験では，水中での溶融物の拡がり挙動が観察され，高温の溶融物は，水張りがあったとしても落下直後は膜沸騰が維持された状態で拡がり，その後拡がりが停止することが確認されている。また，実験結果を元に，水中での溶融物の拡がり距離を求めるスケーリング則が提案されている。
このスケーリング則に，BWR の溶融炉心落下条件を適用して，水中での溶融物の拡がり距離を評価すると，約 18m となる（付録４(5)b）参照）。ペデスタル半径が，ABWR が約 5.3m，Mark-Ⅰ改良型格納容器が約 3.2m であることを考慮すると，水張りをしている場合でも溶融炉心は，落下位置に依らず床全面に拡がる可能性が高いと考えられる。これは，溶融物の拡がりに影響する因子のうち，主に溶融炉心落下流量が大きい（全炉心相当の溶融炉心が 1 分程度で落下することを想定）ことによるものと考えられる。溶融物の拡がり距離については，不確かさを考慮しても床全面に拡がるものと想定される。
次に，粒子状ベッドの堆積形状の不確かさについて考察する。落下直後の初期段階では，DEFOR 実験[16]で確認されたように不均一に堆積することが想定されるが，崩壊熱を発生する
デブリ粒子では，粒子状ベッド内の継続的沸騰による粒子の攪拌によるセルフレベリング
効果により，時間と共に堆積厚さが均一化することが想定される（付録４(4)参照）。但し，溶融炉心がペデスタルの中心からずれた位置に落下する場合を想定すると，セルフレベリ
ング効果により均一化に向かっていくが，偏った堆積形状で MCCI が進行する可能性も考えられる。堆積形状の相違により，溶融炉心から水やコンクリートへの伝熱面積が異なって

5-3-1313

くる。これにより，溶融炉心から水及びコンクリートへの熱配分が異なってくる。MAAP コー
ドでは，均一堆積形状（円柱）として扱うため，不均一な堆積形状は，正確には模擬できない。図 4-1 に示す溶融炉心からの伝熱のうち，溶融炉心から水への伝熱が，溶融炉心の冷却に支配的であり，MCCI の結果に大きく影響することを考慮して，均一堆積形状（円柱）と種々の不均一な堆積形状について，水への伝熱面積を比較した。1100MWe BWR-5 Mark-Ⅰ
改相当の比較結果を図 4-5 に示すが，MAAP コードで想定している均一堆積形状（円柱）が，最も水への伝熱面積が小さくなる。偏った堆積形状を代表する斜切円柱の場合でも，均一堆積形状（円柱）より水への伝熱面積は大きくなっている。円柱状の部分拡がりの場合でも，上面の伝熱面積は減少するが，側面の伝熱面積が加わることにより，均一堆積形状（円柱）より水への伝熱面積は大きくなっている。以上から，堆積形状の不確かさについては，
均一堆積形状（円柱）の扱いが，溶融炉心から水への伝熱を小さくし，結果的にコンクリー
トへの伝熱を大きくしているため，MCCI 評価の観点から保守的な扱いとなっていると考えられる。
以上の考察より，溶融炉心の拡がりについて，粒子状ベッドの堆積形状の不確かさが想定され，溶融炉心から水への伝熱に影響する可能性も考えられることから，拡がりを抑制した場合の感度解析の実施が必要と考えられる。ただし，溶融炉心の拡がりについては，
ペデスタルの形状や事前水張りの深さに依るところが大きいため，個別プラントにおけるこれらの状況を踏まえて感度解析等の取り扱いを行うことが適切と考えられる。

5-3-1414

への伝熱面積の比較
図 4-5 溶融炉心の堆積形状のバリエーションと 水への伝熱面積

5-3-1515

デブリから水への熱伝達
デブリベッドから水への熱流束は，溶融炉心上面の性状により大きく影響されるため，大きな不確かさが想定される。
評価モデル：

溶融炉心上面から水への熱流束は，上部クラスト－水プール間の熱伝達モデルである Kutateladze 型の水平平板限界熱流束相関式(Kutateladze 式)で与えられるが，係数 Fchf を調整することにより，溶融炉心上面の性状に応じた熱流束を設定することが可能である。これは，上面の性状に応じて水への熱流束が内部で計算されるモデルではなく，
Kutateladze 式による限界熱流束で与えるモデルとなっているためである。
評価条件：

事前水張り対策により溶融炉心の一部は粒子化して，デブリ上部に粒子状ベッドを形成し，下部には粒子化しなかった溶融炉心が溶融プールを形成すると想定される。これは，水中に溶融炉心を落下させた実験的知見[14][15][16]（付録３参照）から，安定ク
ラストは形成されず，粒子状ベッドが形成されていることから妥当と考えられる。なお，
SSWICS 実験では，コンクリートがデブリに含有されると上面熱流束が低下する結果が得られているが，これは火山型クラストが形成された場合の知見であり，粒子状ベッドに対する知見ではないと判断される。
デブリ上面に粒子状ベッドが形成されることを想定し，粒子状ベッドから水への熱流束の不確かさを評価する。ここでは，Lipinski 0-D モデル[24]を使用し，主なパラメータである粒子径，ポロシティ，堆積高さ，圧力を変えた場合の粒子状ベッドのドライアウト熱流束の不確かさを評価した。以下のように各パラメータを想定した。
・ 粒子径は，これまでの実験等の知見により 3～7mm と報告されているため[33]，この範囲を不確かさの範囲と想定する。
・ ポロシティは，概ね 0.3 以上が報告されているが[33]， 粒子の充填状態の幅を考慮して，0.26（面心立方格子，最稠密），0.32（体心立方格子），0.4（MAAP 標準値），
0.48（単純立方格子）の範囲を想定する（図 4-6 参照）。
・ 粒子状ベッドの堆積高さは，全炉心相当が落下した場合の堆積高さは 1m 以上となるが，破損口径の拡大を考慮しても 33%以上の粒子化が想定されるため，0.3m 以上とする。
・ 圧力は，MCCI が発生する時間では格納容器圧力は 0.4MPa[abs]以上となっているため（図 4-7a,図 4-7b 参照），0.4MPa[abs]以上とする。
Lipinski 0-D モデルによる評価結果を以下に示す。
・ 図 4-8a にドライアウト熱流束の圧力依存性を示す。上記のパラメータの範囲から，粒子径 3mm，ポロシティ 0.32 の熱流束特性が，現実的な範囲内で保守側の特性と想定され，そのときの熱流束は，大気圧状態で 800 kW/m2 程度となる。また，圧力が0.4MPa[abs]以上での下限値としては，ポロシティの最も小さい 0.26 のケースで 800
kW/m2 程度のドライアウト熱流束となっている。
・ 図 4-8b にドライアウト熱流束の粒子径依存性を示す。粒子径が 3mm 以上での下限値

5-3-1616

程度のドライアウトのドライアウト
としては，ポロシティの最も小さい 0.26 のケースで 800 kW/m2 程度
熱流束となっている。

・ 図 4-8c にドライアウト熱流束の 粒子ベッド高さ依存性を示す。粒子 ベッド高さが，以上ではドライアウト熱流束
熱流束はほとんどはほとんど変化がなく，下限値としては ，ポロシ
0.3m 以上ではドライアウトのドライアウト熱流束となっていとなっていティの最も小さい 0.26 のケースで 800 kW/m2 程度のドライアウト熱流束る。

以上の粒子状ベッドから水
ベッドから水への熱流束 の不確かさ評価結果より，有効性評価条件
有効性評価条件としてとして，大気圧状態で 800 kW/m2 で圧力依存性を考慮使用することは，妥当と判断考慮した上面熱流束を使用することはされる。これは，MAAP コードの Kutateladze 式の係数 Fchf をデフォルトの 0.1 に設定することに相当する。また，デブリ上面熱流束
デブリ上面熱流束の不確かさの下限値
かさの下限値は，800 kW/m 2 程度と考えら
粒子径やポロシティ等の不確 かさによる上面熱流束
よる上面熱流束の不確かさ の影響を確認すれるため，粒子径やポロシティ
るため，800 kW/m2 一定を仮定した感度解析 を実施する。

体心立方格子面心立方格子（最稠密 ）（ポロシティ 0.32）（ポロシティ
（ポロシティ 0.26）

図 4-6 粒子 の格子配列とポロシティ
格子配列とポロシティ

5-3-1717

格納容器限界圧力（854kPa[gage]（最高使用圧力 の 2 倍））

格納容器代替スプレイ作動による 圧力変化

原子炉圧力容器破損
プレナムへの 溶融炉心落下
下部プレナムへの溶融炉心落下

図 4-7a MCCI 評価時の格納容器圧力

原子炉圧力容器破損

図 4-7b MCCI 評価時の溶融炉心温度5-3-1818図 4-8a 圧力とドライアウト熱流束の関係（Lipinski 0-D モデル）

図 4-8b 粒子径とドライアウト熱流束の関係（Lipinski 0-D モデル）

5-3-1919図 4-8c 粒子状ベッド高さとドライアウト熱流束の関係（Lipinski 0-D モデル）

5-3-2020

コンクリート侵食
コンクリートの種類（玄武岩系や石灰岩系）により，侵食挙動やガス発生挙動が異なる。
また，玄武岩系コンクリートに関しては，壁方向と床方向の侵食量に相違が見られる知見が出ている。床方向に大きく侵食する結果（付録１(3) BETA 実験参照）と壁方向に大きく侵食する結果（付録２(5)OECD/MCCI 実験の CCI-1/3 実験参照）があり，この侵食の異方性は，未だ十分な理解が得られておらず，現状モデルでも扱われていない。また，実験ではドライ条件において侵食の異方性が確認されており（CCI-1/3 実験では，最初ドライ条件で侵食させている），実機のウェット条件でもこの侵食の異方性が発生するかは不明であるが，念のためモデルパラメータの感度解析により，この侵食異方性の影響を検討する。
評価モデル：・ コンクリート表面から深さ方向に１次元熱伝導を解き温度分布を計算している。
・ クラスト- コンクリート界面の温度がコンクリート融点以上である場合にコンク
リート侵食が計算される。
・ コリウムからの伝熱量と分解／溶融潜熱により侵食量を評価する。
・ コンクリート分解による自由水・結合水・CO2の発生を考慮し，溶融コリウム中での化学熱力学平衡計算による未酸化金属との反応によりH2及びCO発生が計算される。
反応熱は，溶融炉心の崩壊熱に加算される。
評価条件：
コンクリートの種類には，玄武岩系や石灰岩系があり，プラントに使用されているコン
クリートの種類により組成や物性値を使い分けている。表 4-1 にコンクリートの物性値，
表 4-2 にコンクリートの組成例を示しており，これらを MAAP の入力としている。融点は石灰岩系コンクリートの方が高いが，比熱は石灰岩系コンクリートの方が小さい。コンクリー
ト成分は骨材（aggregate）の種類に影響され，玄武岩系コンクリート（Basaltic）の成分
は SiO2 が主体であるが，石灰岩系コンクリート（Limestone）では，CaCO3 が主体である。
コンクリートの侵食モデルは，主に SURC 実験や ACE/MCCI 実験の 1 次元侵食の実験や，BETA 実験の 2 次元侵食のドライ実験に対して検証されている（付録１参照）。なお，コンク
リート侵食モデルは，コンクリート表面から侵食深さ方向に一様に侵食する扱いをしているのに対し，実験で得られた侵食深さは，一般的に侵食面で一様ではない。このため，コンクリート侵食モデルの検証では，実験の平均的な侵食深さと比較して侵食速度が良く一致していることが確認されている。
侵食面での侵食の不均一性と MAAP コードの一様侵食との関係を見るため，SURC-4 実験（1次元侵食）に対する MAAP コードのベンチマーク解析で得られた侵食深さに対して，実験で得られた侵食深さがどの程度の幅にあるかを確認した結果を図 4-9 に示すが，実験で得られた侵食深さは，MAAP コードによる予測侵食量に対して概ね±20%程度の範囲にある。同様に BETA-V5.1 実験（2 次元侵食）に対して確認した結果を同図に示すが，この場合も概ね±20%程度の範囲にある。従って，実験における侵食面での侵食深さは，MAAP コードによる予5-3-2121

測侵食量の概ね±20%程度の範囲にあると考えられる。
また，壁方向と床方向の 2 次元侵食を調べる OECD/MCCI プロジェクトの CCI 実験では，石灰岩系コンクリート（CCI-2 実験）と玄武岩系コンクリート（CCI-1 実験，CCI-3 実験）の侵食挙動が調べられている。
石灰岩系コンクリート（CCI-2 実験）では，壁方向と床方向に等方的に侵食している。一方，玄武岩系コンクリート（CCI-3 実験）では，壁方向と床方向の侵食の異方性が確認されており（付録２参照），玄武岩系コンクリートの壁方向と床方向の侵食速度の相違から，壁方向の熱流束が，床方向の熱流束の約 4 倍程度大きいと推定されている。なお，同様に玄武岩系コンクリートが使われた CCI-1 実験では，加熱条件が他の実験と異なり高く設定されていたため，代表性があるとは見なされず，侵食挙動の分析やベンチマーク解析には，
CCI-3 実験の結果の方が使われている。玄武岩系コンクリートの壁方向と床方向の侵食の異方性のメカニズムに関しては，明確な理解は得られていないが，CCI-3 実験に対する MAAP
コードのベンチマーク解析[35]において，壁方向への熱流束が床方向の 4 倍程度大きくなる
ように溶融プールからの対流熱伝達係数を設定することにより，実験の侵食挙動と良く一致する解析結果が得られている。この場合も，実験の侵食深さは，MAAP コードの予測侵食量に対して，上記の SURC-4 実験及び BETA-V5.1 実験のベンチマーク解析で設定した 20%以内の侵食深さに収まっている（図 4-10 参照）。
以上の MAAP コードのベンチマーク解析の分析から，コンクリート侵食の不確かさについて，以下のことが言える。
侵食面における侵食深さのばらつきの範囲は，MAAP コードの予測侵食量の 20%以内に入っている。
玄武岩系コンクリートの壁方向と床方向の侵食の異方性は，MAAP コードにおいて溶融プールからの対流熱伝達係数を変更することにより予測可能である。その場合も侵食のばらつきの範囲は，MAAP コードの予測侵食量の 20%以内に入っている。
以上のことから，コンクリート侵食の不確かさに関して，侵食面における侵食深さの不均一性については，概ね MAAP コードによる予測侵食量の 20%の範囲内であることがわかった。一方，玄武岩系コンクリートの壁方向と床方向の侵食の異方性については，実機における影響の確認のため感度解析が必要と考えられ，CCI-3 実験のベンチマーク解析を参考にして，溶融プールから壁方向と床方向の対流熱伝達係数を調整して感度解析を実施する。
壁方向の熱配分の最大ケースとして，床方向の 4 倍として感度解析を行い，コンクリート
侵食の壁方向と床方向の異方性の影響を確認する。
以上の各過程での，MCCI 現象の影響因子と感度解析パラメータの関係を図 4-11 に示し，
不確かさ要因，有効性評価の扱い及び感度解析の要否を表 4-3 に整理する。
また，溶融炉心とコンクリートの相互作用における実機で想定される現象と解析上の取り扱いの比較を図 4-12 に示す。

5-3-2222

本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

表 4-1 コンクリートの物性比較 (MAAP 入力)
項目
コンクリート組成
液相線温度 （K）
固相線温度 （K）
比熱 (J/kg-K)
溶融潜熱 (J/kg)

玄武岩系コンクリート
SiO2 が主体

石灰岩系コンクリート
CaCO3 が主体

表 4-2 代表的なコンクリートの組成例(重量%)成分
SiO2
TiO2,MnO,MgO
CaO
Na2O
K 2O
Fe2O3
Al2O3
Cr2O3
CO2
H2O（自由水，結合水）

玄武岩系コンクリート
54.84
7.21
8.82
1.80
5.39
6.26
8.32
0.00
1.50
5.86

5-3-2323[37]

石灰岩系コンクリート
3.60
5.80
45.40
0.0078
0.68
1.20
1.60
0.004
35.6985.94

本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

図 4-9

SURC-4 と BETA-V5.1 実験と MAAP4 コード予測の侵食量の比較[23]5-3-2424

本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

図 4-10 CCI-3 実験（玄武岩系コンクリート）の侵食深さ[12]5-3-25

25

6
2
3
5

26

MCCI 進行と注水によるデブリ冷却過程

溶融物の落下・堆積過程

下部ヘッド破損過程

MCCI に至る過程

図 4-11 MCCI における不確かさに関する流れ図（床/壁への熱分配の不確かさ）

溶融プール－クラスト間の熱伝達係数

コンクリートへの熱伝達
コンクリート種類（個別プラント条件）

溶融炉心の拡がり

上面熱流束係数（Fchf）

エントレインメント係数（全炉心相当を想定しており，
感度解析不要）（下部ヘッド貫通部の最大口径の破損を想定しており感度解析不要）

感度解析パラメータ

溶融炉心の拡がり

水への熱伝達デブリ上面性状

プール水深

溶融炉心ジェット径

溶融炉心の粒子化

溶融炉心落下量

下部ヘッド破損モード

MCCI 現象の影響因子

7
2
3
5

27
コンクリートの種類により，侵食挙動やガスの発生挙動が異なる。

コンクリート種類に応じて玄武岩系や石灰岩系を模擬可能。溶融プールからクラストへの対流熱伝達係数を変更することにより熱移行の異方性を評価可能。

解析上の取り扱い

図 4-12 溶融炉心とコンクリートの相互作用における実機で想定される現象と解析上の取り扱いの比較

水深が浅い場合は，溶融物
は完全には粒子化せず，床上に溶融炉心が拡がり，その上に粒子状ベッドが堆積する。水により拡がりが抑制され，不均一に堆積する。上面の性状により水への熱伝達が大きく影響される。

溶融物が水中を落下する
過程で粒子化する。この際デブリ粒子のクエンチにより水蒸気が発生し，デブ
リ粒子の酸化により水素が発生する。

下部ヘッド破損口を溶融
炉心が通過する際に破損口側面が溶かされ，破損口が拡大する。

実機で想定される現象

溶融炉心は，床上で均一に堆積し，均一に混合して外周部のクラスト層と中心部の溶融プールから構成される平板の発熱体として模擬。上面
熱流束は Kutateladze の水
平平板限界熱流束相関式で与えており，係数変更により
粒子状ベッドから水への熱流束の模擬が可能。

溶融物が水中を落下する過程での粒子化を模擬。粒子化
割合は Ricou-Spalding モデ
ルにより評価。この際デブリ粒子のクエンチによる水蒸気発生と，デブリ粒子の酸化による水素発生を模擬。

下部ヘッド破損口を溶融
炉心が通過する際の破損
口側面の溶融による破損
口拡大を模擬。

8
2
3
5

28

水による拡がり抑制

デブリ上面の性状

壁方向と床方向の熱分配

デブリから水への熱伝達

コンクリート種類
玄武岩系コンクリートの
侵食の異方性

溶融炉心ジェット径

粒子状ベッドの熱流束として，現実的な範囲内で保守側の想定をして，大気圧状態で800kW/m2 として圧力依存性を考慮
溶融プールが完全な液相の場合の対流熱伝達係数として，床方向クラストへは，壁方向クラストへはを使用

拡がりの知見から，全面に拡がることを想定し，均一堆積モデルを使用

破損口径に対応した径を考慮

床方向と壁方向の侵食の異方性の影響を見るため，CCI-3 実験で確認された熱分配比率を想定した感度解析を実施。

粒径やポロシティの不確かさを考慮して，下限の熱流束に対して感度解析を実施。

実験解析を元に，粒子化割合を少ない側に評価する値を使用

エントレインメント係数

溶融炉心の拡がり

溶融炉心の粒子化

解析モデルパラメータではないため，感度解析不要。
エントレインメント係数に不確かさがあるため，MAAP コードの不確かさ範囲内(FARO 実験解析に基づく)で感度解析を実施。
最も大きい径を想定して，粒子化割合を小
さくし，水に落下した際のデブリクエンチ
量を小さくして厳しい側の扱いをしているため，感度解析不要。
溶融炉心の拡がりを抑制した場合を想定
し，水張り深さ等の個別プラントの状況を踏まえた感度解析等の取り扱いを行うことが必要

ペデスタルの注水開始条件及び注水流量について，手順書規定に準じた操作を想定

プール水深

原子炉圧力容器内溶融進展

溶融炉心落下量

感度解析の要否
下部ヘッド貫通部の中で最も大口径の貫通部の破損を想定し，破損口径の拡大も考慮しているため，感度解析不要。
全炉心相当の最も厳しい落下量を想定しているため，感度解析不要。

破損部位
破損口侵食拡大

下部ヘッド破損モード

有効性評価の扱い
下部ヘッド貫通部の中で最も大口径の制御棒駆動機構ハウジングの瞬時破損を想定し，破損口径の拡大を考慮

MCCI 評価の不確かさに関する整理結果

全炉心相当を想定

不確かさ要因表 4-3

MCCI 現象への影響因子属しますので公開できません。

本資料のうち，枠囲みの内容は商業機密に

本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

５． 感度解析と評価
ペデスタル床上の溶融炉心の堆積高さが高い方が，溶融炉心の冷却性が悪くなり，MCCIの影響が厳しくなる。国内 BWR では，全炉心が落下した場合には，溶融炉心の堆積高さはMark-I 改や Mark-II で約 1m 程度であり，ABWR では約 0.5m 程度である。このため，ここでは，溶融炉心の堆積高さの高い BWR-5, Mark-I 改プラントを対象に感度解析を行う。前節の検討を踏まえ，BWR-5, Mark-I 改プラントの大 LOCA シーケンスを対象に抽出されたパラメータの感度解析を実施した。
また，ガスの発生の影響を見るために参考解析として石灰岩系コンクリートを想定した解析も実施した。

(1) エントレインメント係数の感度解析
エントレインメント係数は，ベースケースではを設定しているが，感度解析ケー

スでは，MAAP コードの FARO 実験解析を元に設定されている当該変数の不確かさ範囲のうち，最小値であると最大値であるの感度解析を実施する。感度解析ケースを表 5-1 に示す。解析結果を図 5-1-1a～図 5-1-2b 及び表 5-5 に示す。
エントレインメント係数が小さいと粒子化割合が減り，落下時のクエンチによる除熱量が減るため，床に堆積する溶融炉心の温度はベースケースより高めに推移し，コンクリートへの伝熱量も増加するため，侵食量は若干増えている。エントレインメント係数が大きいと粒子化割合が増えて逆のことが起こり，侵食量は若干減少している。
エントレインメント係数の不確かさ幅の範囲では，コンクリート侵食量への影響は小さい。

表 5-1 エントレインメント係数のベースケースと感度解析ケース
パラメータ①エントレイ
ンメント係数ベースケース（有効性評価）

感度解析ケース

(1)
(2)

（小）
（大）

5-3-2929

根拠(1)MAAP 不確かさ範囲*の最小値(2)MAAP 不確かさ範囲*の最大値*FARO 実験解析を元に設定

(2) 上面熱流束の感度解析
感度解析ケースを表 5-2 に示す。上面熱流束の不確かさを考慮した下限値として，
Lipinski 0-D モデルにより評価される粒子状ベッドのドライアウト熱流束において，想定される範囲でドライアウト熱流束を小さくする条件の値を選定することとし，圧力 0.4MPa[abs]，粒子径 3mm，及びポロシティ 0.26 を条件とした場合の熱流束である
800kW/m2（一定）を選定した。解析結果を図 5-2a～図 5-2b 及び表 5-5 に示す。
上面熱流束が小さいと水への熱伝達量が減少するため，溶融炉心の温度はベースケー
スより高めに推移し，コンクリートへの伝熱量も増加するため，侵食量は増えている。
上面熱流束の不確かさ幅の範囲では，コンクリート侵食量への影響は大きい。

表 5-2 上面熱流束係数のベースケースと感度解析ケース
パラメータ

ベースケース（有効性評価）②上面熱流束係数（Fchf）

Fchf=0.1 (default)
圧力依存性あり(大気圧状態で約800kW/m2)

感度解析ケース

800kW/m2 一定

5-3-3030

根拠
Lipinski 0-D モデルによる評 価 に お い て ， 圧 力0.4MPa[abs]，粒子径 3mm，
ポロシティ 0.26 のドライア
ウト熱流束本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

(3) 溶融プール－クラスト間の熱伝達係数の感度解析
感度解析ケースを表 5-3 に示す。熱伝達係数の不確かさ幅の中で，床方向への熱伝達を大きくしたケースと壁方向への熱伝達を大きくしたケースの解析結果を図 5-3-1a～
図 5-3-2b 及び表 5-5 に示す。
床方向への熱伝達を大きくすると，床方向へ侵食量は増えており，壁方向の侵食は若干減少している。壁方向への熱伝達を大きくすると，壁方向へ侵食量は増えており，床方向の侵食はあまり変わらない。
溶融プールの熱伝達係数の不確かさ幅の範囲では，コンクリート侵食量への影響は小さい。
表 5-3 溶融プール－クラスト間の熱伝達係数のベースケースと感度解析ケース
パラメータ③溶融プール－クラスト間の熱伝達係数ベースケース（有効性評価）
溶融プールが完全な液相の場合の対流熱伝達係数として;
床方向
壁方向

感度解析ケース(1) 床：
壁：(2) 床：
壁：

5-3-3131根拠(1)
床：壁の 1/4(CCI-3 実験反映)
壁：推奨入力値(2)
床：推奨入力値
壁：床の 1/4((1)の逆とした。)

(4) 石灰岩系コンクリートを想定した参考解析
解析条件を表 5-4 に示す。コンクリートの組成や物性は，表 4-1 及び表 4-2 に示す値を使用している。石灰岩系コンクリートを想定した場合の解析結果を図 5-4a～図 5-4d
及び表 5-5 に示す。石灰岩系コンクリートの場合の，床方向及び壁方向の侵食量は若干増加しているが，感度は小さい。また，格納容器内のガスモル分率もベースケースと同様な結果となる。このケースのようにコンクリート侵食量が小さい場合は，石灰岩系コンクリートの影響は小さい。
表 5-4 コンクリート種類を変えた参考解析条件
パラメータ

ベースケース（有効性評価）参考解析ケース根拠

コンクリート種類

玄武岩系

石灰岩系

ガス発生の影響を確認する。

5-3-3232

本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

表 5-5
ケース

有効性評価
（ベースケース）①-1 エントレインメン
ト係数 小①-2 エントレインメン
ト係数 大② 上面熱流束 下限③-1 溶融プール熱伝達
床方向大③-2 溶融プール熱伝達
壁方向大（参考解析）
石灰岩系コンクリート

MCCI 感度解析結果の一覧

パラメータ設定
・エントレインメント係数:
・上面熱流束: 約 800kW/m2 @ 1ata
(Fchf=0.1 デフォルト)
圧力依存性考慮・溶融プールからクラストへの熱伝達係数
床壁・コンクリート組成及び物性値：
玄武岩系
エントレインメント係数:
エントレインメント係数:
800kW/m2 一定
溶融プール熱伝達係数:
床壁
溶融プール熱伝達係数:
床壁
コンクリート組成及び物性値：
石灰岩系

5-3-3333

コンクリート侵食量

床面：
壁面：

2.4cm
2.3cm

床面： 3.6cm
壁面： 3.5cm
床面： 1.7cm
壁面： 1.7cm
床面： 22.5cm
壁面： 21.6cm
床面：
壁面：

2.4cm
2.2cm

床面：
壁面：

2.1cm
2.4cm

床面：
壁面：

2.8cm
2.8cm

本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。
「本製品（又はサービス）には，米国電力研究所（The Electric Power ResearchInstitute）の出資により電力産業用に開発された技術が取り入れられています。
」

図 5-1-1a コンクリート侵食量の比較（エントレインメント係数

図 5-1-1b 溶融炉心温度の比較（エントレインメント係数

5-3-3434））

本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

図 5-1-2a コンクリート侵食量の比較（エントレインメント係数

図 5-1-2b 溶融炉心温度の比較（エントレインメント係数

5-3-3535））

コンクリート侵食量の比較（上面熱流束 800kW/m2 一定 ）図 5-2a コンクリート侵食量

図 5-2b 溶融炉心温度の溶融炉心温度の比較（上面熱流束 800kW/m2 一定）

5-3-3636

図 5-3-1a コンクリート侵食量の比較（床方向対流熱伝達係数大）

図 5-3-1b 溶融炉心温度の比較（床方向対流熱伝達係数大）5-3-3737

図 5-3-2a コンクリート侵食量の比較（壁方向対流熱伝達係数大）

図 5-3-2b 溶融炉心温度の比較（壁方向対流熱伝達係数大）5-3-3838図 5-4a コンクリート侵食量の比較（石灰岩系コンクリート）

図 5-4b 溶融炉心温度の比較（石灰岩系コンクリート）

5-3-3939

図 5-4c ドライウェルガスモル分率（石灰岩系コンクリート）

図 5-4d サプレッションチェンバガスモル分率（石灰岩系コンクリート）

5-3-4040

６．まとめMCCI に関する種々の実験から得られた知見等に基づき MCCI の各過程における不確かさの分析を行い，不確かさがあり，評価結果に影響するモデルパラメータとして，エントレイ
ンメント係数，上面熱流束を抽出した。また，2 次元侵食実験結果で確認されている侵食の異方性については，未だ十分な解明が行われておらず，実機のウェット条件でも起こるのかは不明であるが，念のために異方性の影響をみるために，溶融プールから床方向と壁方向への熱伝達係数の感度解析を行った。これらのパラメータの感度解析以外に，コンクリー
ト種類を玄武岩系コンクリートから石灰岩系コンクリートに変更した場合の解析も参考に実施し，コンクリート侵食量への影響を確認した。
BWR-5, Mark-I 改プラントの有効性評価に使用しているノミナル条件を設定したベース
ケースでは，コンクリート侵食量は 2.4cm であった。
エントレインメント係数の感度解析結果より，水中での粒子化によるクエンチ効果に差が出て，床上に堆積した溶融炉心温度に影響し，コンクリート侵食量に影響するが，その影響は小さかった。
上面熱流束の感度解析結果より，上面熱流束は，床上に堆積した溶融炉心温度に大きく影響し，コンクリート侵食量に大きく影響する結果となった。
玄武岩系コンクリートに見られる床方向と壁方向のコンクリート侵食の異方性の影響を確認するため，溶融プールから床方向と壁方向への熱分配の異方性として扱い，溶融プー
ル－クラスト間の熱伝達係数を変えて，床方向への熱伝達が大きい場合と壁方向への熱伝達が大きい場合の感度を調べた。その結果，熱伝達の大きい方向の侵食量が大きくなり，
侵食の異方性が確認された。但し，ノミナル条件の場合のように侵食量が小さい場合は，侵食量に大きな感度がないことがわかった。
感度解析の結果，コンクリート侵食量に対して上面熱流束の感度が支配的であることが確認された。上面熱流束を想定される下限値とした場合でも，コンクリート侵食量は，22.5cm 程度に収まることが確認された。
また，実験で確認されている侵食面における侵食の不均一性については，MAAP モデルのベンチマーク解析結果から，実験の侵食面における侵食のばらつきが，MAAP コードによる予測侵食量の 20%の範囲内に収まっていることが確認され，上面熱流束の感度に比べて影響が小さいことが確認された。
また参考として，他の条件を同一として玄武岩系コンクリートから石灰岩系コンクリー
トに変更した解析も実施したところ，侵食量に大きな感度は無かった。また，ガス発生の影響も顕著には見られなかった。
最後に，MCCI に関係する現象は，複雑な多成分・多相熱伝達現象であり，現状でも知見が十分であるとは言えない。また事前水張り時の落下デブリの冷却性を直接調べた実験例が殆どないことから，今後も継続して検討を進め，新たな知見が出た場合には適宜評価に反映するように努めることが重要であると考えられる。

5-3-41

41

付録
付録 1 ドライ条件での MCCI 実験
デブリとして UO2 コリウムを用い，崩壊熱を模擬し，注水を行っていないドライ条件の実験として，米国サンディア国立研究所（SNL）の SURC 実験，米国アルゴンヌ国立研究所（ANL）の ACE 実験及び独カールスルーエ研究所（KfK）で実施された BETA 実験がある。以下にそれぞれの実験について説明する。これらの実験は，実機における MCCI 継続の影響の検討に参考となる。
(1) SURC 実験（米国 SNL）[4][23][34]MCCI 時の伝熱，化学反応，ガス及びエアロゾルの放出等を調べ，CORCON 等の解析コード
を検証・改良することを目的に SNL にて 4 回実験が実施された。そのうち実機と同様の組成の UO2 デブリを使用した実験は 2 回実施された（SURC-1,2。他の 2 回の実験は模擬デブリとして SUS304 が使用された）。実験装置の概要を図付 1-1 に示すが，直径 40cm のコンクリー
ト試験体（SURC-1 実験では石灰岩系（米国で使用されている），SURC-2 実験では玄武岩系のコンクリートを使用）を入れた MgO 製の円筒容器に約 250kg の模擬デブリ（組成：UO2 69％，
ZrO2 22%, Zr 9%）を誘導加熱して溶融させてコンクリートとの反応を継続させた。SURC-2のコンクリートの侵食深さの変化を図付 1-2 に示すが，侵食速度は 1 時間に 20cm 程度となっている。
また，SURC-4 は，コード比較のための国際標準問題（ISP－24）[5]に選定され，実験デー
タと種々の解析の比較が行われている。200kg のステンレス鋼と模擬 FP （Te: 0.5 kg，
La2O3: 1.17 kg， CeO2 ：1.23 kg， BaO：1.1 kg）が，コンクリート侵食が開始するまで加熱され，侵食開始後 14 分経ってから約 0.5 秒間に追加的に 20kg のジルコニウムが溶融物に添加されている。SURC-4 のコンクリートの侵食深さの変化を図付 1-3 に示すが，侵食速度は約 55 分の時点で 24.5～27.5cm で，MAAP 解析結果とよい一致を示している。なお，SURC-4
実験では外周部のコンクリート侵食深さが大きくなっており，MAAP 解析結果と若干相違が見られるが，これは実験では外周部から誘導加熱して溶融物を加熱しているため，外周部の溶融物への入熱量が大きくなることが原因と考えられる。従って，MAAP との相違は，実験固有の原因と考えられるため，ベンチマーク解析の観点からは問題とはならないと考えられる。

5-3-4242

図付 1-1

SURC 実験装置[23]図付 1-3

図付 1-2 SURC-2 実験の侵食深さ[23]SURC-4 実験の侵食深さと MAAP 解析との比較[23]5-3-4343

(2) ACE 実験（米国 ANL）[3][23][34]MCCI における熱水力学的及び化学的プロセスを検証し関連コードのデータベースを拡充することを目的に，国際プログラムとして ACE（Advanced Containment Experiments）計
画の一部として実施された。Phase-C で MCCI 時の FP エアロゾル放出の定量化の目的で実験が実施された。実験装置の概要を図付 1-4 に示す。実炉組成の約 300kg の模擬デブリを使用し，タングステン電極により直接通電により加熱し溶融させている。コンクリート侵食
はベースマット中にある熱電対によりモニターされている。実験マトリックスを表付 1-1に示すが，実験はコンクリートの種類やデブリの組成を替えて実施された。
L2 実験のコンクリートの侵食深さの MAAP 解析との比較を図付 1-5 に示すが，MAAP 解析
結果とよい一致を示している。

5-3-4444

表付 1-1

ACE/MCCI 実験マトリックス[23]出典：MAAP4 User’s Manual, EPRI

図付 1-4

ACE／MCCI 実験装置[3]図付 1-5

ACE-L2 実験と MAAP 解析の
侵食深さの比較[34]5-3-4545

(3) BETA 実験（独 KfK）[6][34]MCCI 解析コードの侵食挙動やエアロゾル生成挙動の検証用のデータベースを拡充することを目的として，独 KfK で実施された。実験装置の概要を図付 1-6 に示す。テルミット反応により金属溶融物(Fe,Cr,Ni,Zr)と酸化物溶融物(Al2O3,SiO2,CaO) を生成させて，玄武岩系コンクリート製のるつぼに落下させて，るつぼ内の溶融物は，誘導加熱により加熱している。ベースマット中にある熱電対により 2 次元のコンクリート侵食挙動が確認できるよ
うになっている。溶融物として Zr を多く含む溶融物を用いた実験シリーズの実験マトリックス(V5.1～V5.3)を表付 1-2 に示す。
実験後のコンクリート侵食状況を示す模式図を，図付 1-7 に示す。壁方向よりは，床方向の侵食量が大きくなっている。溶融物に含まれる Zr によりコンクリート成分の SiO2 が還元され，Si や SiO の生成が確認されている。V5.1 実験と V5.2 実験のコンクリートの侵食
深さの MAAP 解析との比較[34]を図付 1-8 及び図付 1-9 に示すが，MAAP 解析結果とよい一致を示している。この MAAP のベンチマーク解析により，2 次元侵食に関わるパラメータ設定の妥当性が判断されている。

図付 1-6BETA 実験装置 [34]5-3-4646

表付 1-2

図付 1-7

BETA 実験マトリックス[6]BETA 実験後のコンクリート侵食状況を示す模式図 [6]5-3-4747

本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

図付 1-8BETA-V5.1 実験と MAAP 解析の侵食深さの比較 [34]

図付 1-9BETA-V5.2 実験と MAAP 解析の侵食深さの比較 [34]5-3-4848

付録２ 注水を伴った MCCI 実験
注水を伴う実験（溶融物上に冷却水を注水した実験）としては，SNL で実施された SWISS
実験及び WETCOR 実験，ANL（EPRI）の主催で実施された MACE 実験及び OECD/MCCI 実験，
原子力発電技術機構（NUPEC）により実施された COTELS 実験がある。
以下の実験の知見を示す。
① SWISS 実験（米国サンディア研究所(SNL)）② WETCOR 実験（米国サンディア研究所(SNL)）③ MACE 実験（米国アルゴンヌ研究所(ANL)）④ COTELS 実験（日本 原子力発電技術機構（NUPEC））⑤ OECD/MCCI 実験（米国アルゴンヌ研究所(ANL)）(1)SWISS 実験（米国 SNL）[7]高周波加熱によりステンレスを溶融させ，溶融したステンレスとコンクリートとの反応中に注水し，溶融物の冷却性を調べる実験が 2 回行われた。実験結果 2 回とも同様の結果となっており，メルト上面に安定なクラストが形成されてメルト内部に冷却水が浸入しにくくなり，メルトによりコンクリートの侵食が継続した（図付 2-1 参照）。但し，本実験はメルトの発熱が実機で想定される値より１桁程度大きいことや l00％ステンレスであるため模擬崩壊熱に加え，水－金属反応の発熱が加わったことも一因と分析されている。また
コリウムより上方水プールへの熱流束は約 800kW/m2 程度であった（図付 2-2 参照）。

5-3-4949

図付 2-1

SWISS 実験の実験装置とコンクリート侵食深さの実験と解析（MAAP）の比較[7]5-3-5050図付 2-2

SWISS-2 実験の水プールへの熱流束[7](2) WETCOR 実験（米国 SNL）[8]Al2O3 と CaO の酸化物混合溶融物 34 kg を用いてコンクリートとの反応中に注水し溶融物の冷却性を調べる実験が行われた（図付 2-3 参照）。SWISS 実験と同様の理由によりコンクリートの侵食が継続した。コンクリート侵食速度は SWISS 実験と比べて約 l/3 に低下している。水プールへの熱流束は溶融時に 520kW/m2 程度，凝固時に 200kW/m2 程度であった。

5-3-5151

（実験装置の概要）（実験後の状態の模式図）（上部クラストの様子）（コンクリート侵食の時間変化）
図付 2-3

WETCOR 実験結果[8]5-3-5252

(3)MACE 実験（米国 ANL）[3][9][10][20]MACE 実験（Melt Attack and Coolability Experiment）は，米国電力研究所（EPRI）

が主催する格納容器内溶融炉心冷却性に関する国際協力計画である。この実験では，溶融炉心がコンクリート床に落下した後に，溶融炉心の上に注水する状況を模擬しており，大規模実験が含まれていることと，実機相当の UO2 及び ZrO2 を主成分とする溶融炉心模擬物を用いて，さらに，直接通電加熱により崩壊熱の発生を模擬していることが特徴である。実験装置の概要を図付 2-4 に示す。MACE 実験では規模の異なる M0，M1b，M3b，M4 の 4 回の実験が実施されており，その主要条件を表付 2-1 に示す。M0 実験のみコンクリートの側壁で，その他は MgO が用いられている。
M3b と M4 の実験結果を図付 2-5 及び図付 2-6 に示す。いずれのケースもコンクリート
侵食が継続する結果となっている。実験後の溶融炉心模擬物の固化状況を見ると，メルト上部に安定クラストが形成され側壁に固定化されており，クラスト下部に空洞が見られ，
クラストと溶融物とが分離した状態になっている。メルトから冷却水への熱流束を見ると，注水直後の短期間は，大きい熱流束となっているが，安定クラストが形成されクラストと溶融物が分離したと想定される状況では，熱流束は低下している。安定クラストの形成とメルトとの分離によりクラストの下にあるデブリの冷却が阻害される結果となったと分析されている。しかしながら，デブリの長期的な冷却に関与し得るメカニズムとして，クラストに生じる亀裂等への浸水及びコンクリート分解ガスによるクラスト開口部からのデブリの噴出が確認されている。

5-3-5353

表付 2-1

図付 2-4MACE 実験の主要条件[20]MACE 及び OECD/MCCI 実験装置[9]5-3-5454

（軸方向侵食データ）（メルトから水への熱流束）（実験後状態の模式）
図付 2-5MACE-M3b 実験結果[10]5-3-5555

（軸方向侵食データ）（メルトから水への熱流束）（実験後状態の模式）
図付 2-6

MACE-M4 実験結果[10]5-3-5656

(4) COTELS 実験（日本 NUPEC）[1][11]プール水中に実機組成の UO2 コリウムを落下させる実験（FCI）とコンクリート上に落下した UO2 コリウムに注水する実験(MCCI)を実施している。図付 2-7 にコンクリート上の UO2コリウムに注水する実験の実験装置を示す。約 60kg の実機組成の UO2 コリウムを電気炉で加熱溶融させ，玄武岩系のコンクリート床上に落下させ，その後高周波誘導加熱によりコンクリート上のコリウムを加熱して崩壊熱を模擬している。側壁と床面の両方の侵食を考慮しており，コリウムの高さと直径の比を実験パラメータとしている。注水後にコンクリー
ト温度は低下している。また，侵食深さは 2cm 程度であった。また，図付 2-8 に実験後に固化したコリウムの断面を観察した結果を示すが，上部は粒子状になっており，下部はイ
ンゴット状になっているのが観察されている。また，インゴット状の部分にはひび割れや流路が見られる。
COTELS 実験では，いずれもデブリによるコンクリート侵食が生じているが，先行の SWISS，
WETCOR，MACE とは異なり，コンクリート侵食は停止している。これは，以下の理由と分析されている。
・デブリの粒子化と塊状デブリ上への粒子化デブリの堆積・デブリとコンクリートとの境界における粗骨材を主体とする多孔質層の形成
・粗骨材層及び塊状デブリに形成された流路への浸水
溶融物から水プールへの熱流束は水プールへの熱流束は 100～650kW/m2 程度であったが，これらの熱流束は限界熱流束よりも低く，水がさらに高い除熱能力を有する可能性がある
ことが示唆されている。

図付 2-7

COTELS 実験装置（コンクリート上の UO2 コリウムに注水する実験）[11]5-3-5757

（コンクリート床温度変化）（実験後のコリウム断面の模式図）
図付 2-8COTELS C-5 実験の結果[11]5-3-5858

(5) OECD/MCCI 実験 (米国 ANL)[10][12][13][35][38][41]OECD/MCCI 実験は，MACE 実験や COTELS 実験の注水を伴う MCCI 実験で同定された個々のデブリ冷却メカニズムや火山型クラストの冷却性（図付 2-9 及び図付 2-10 参照）を定量的に解明することを目的として，以下の実験が実施された。
① クラスト浸水・強度実験 (SSWICS 実験)
SSWICS 実験は，クラスト内への浸水メカニズムと浸水による除熱量，さらにクラス
ト強度を調べるために実施された。実験装置の構造を図付 2-11 に示し，実験条件一覧を表付 2-2 に示す。テスト部内の非侵食性の床面上で，テルミット反応により UO2，ZrO2
及びコンクリート成分を主成分とするデブリ模擬溶融物を生成しその上に注水する。
崩壊熱とコンクリート分解気体の発生は模擬しない過渡クエンチ実験である。最大除熱量はクラスト上面におけるドライアウト熱流束により評価される。実験結果より，安定なクラストは形成されず，クラスト内に冷却材が浸入することにより，熱伝導律速以上の熱流束を期待することができるとしている。更に，SSWICS 実験で得られたドライアウト熱流束を実機の事故シーケンス解析において直接的に適用できるよう，
Lister/Epstein のドライアウト熱流束モデルへの適応が行われた。実験結果とこれに基づいた Lister/Epstein モデルによる予測の比較を図付 2-12 に示す。この実験結果
及びこれらに基づいた予測値によれば，コンクリートを全く含有しないデブリに対す
る上面水プールへのドライアウト熱流束は，約 400kW/m2 であり，コンクリート含有とともに熱流束は低下し，含有率 15%程度以上では約 100kW/m2 となっている。
② 溶融物噴出実験(MET 実験)
コンクリート分解気体がクラストの開口部を通過する際に，気体に随伴されてクラスト上に噴出するデブリ量の定量化を目的としている。クラスト浸水実験と同様に，テスト部内の非侵食性の床面上で，テルミット反応により UO2，ZrO2 及びコンクリート
成分を主成分とするデブリ模擬溶融物を生成しその上に注水する。クラスト浸水実験との違いは，本実験においては直接通電加熱により崩壊熱を模擬すること及び非侵食性床面に設けた多数の細孔からコンクリート分解気体を模擬する非凝縮性気体を供給している。テスト部で発生する水蒸気量の変化と実験後分析から得られるデブリ噴出量に基づいて，供給気体流量と噴出量の関係を評価している。
③ 長期 2 次元コンクリート侵食実験(CCI 実験)

デブリによる長期的な床部及び側壁の侵食挙動を明らかにすることを目的に実施された。実験装置のテスト部の構造を図付 2-13 に示し，実験条件一覧を表付 2-3 に示す。
テルミット反応により生成した UO2，ZrO2 及びコンクリート成分を主成分とするデブリ模擬溶融物を直接通電により加熱している。電極がある 2 面は非侵食性の側壁を採用し，残りの 2 面をコンクリートにより製作している。使用するコンクリートは実験
により種類が変えられている。注水なしの条件で，側壁のコンクリートの最大侵食深さが 30cm に達するまで侵食を継続した後に注水を行い，注水後，最大侵食深さが 35cm

5-3-5959

に到達した時点で実験を終了している。実験装置に荷重負荷装置を設置し，注水後に形成されるクラストに荷重を与えて，クラストの機械的強度に関するデータを得ている（図付 2-14 参照）。コンクリート種類を変えて 3 回実験が行われている。CCI 実験のコンクリート組成を表付 2-4 に示す。CCI 実験のデブリ-水間の熱流束を図付 2-15 に示
し，溶融物の平均温度を図付 2-16 に示す。また，実験後の侵食状態の模式図を図付 2-17に示す。石灰岩系コンクリートの場合，床方向と側壁方向の侵食は同程度となっているのに対し，玄武岩系コンクリートの場合，側壁方向の侵食が床方向より大きくなっており，侵食に異方性があることが示されている。表付 2-5 に示すように，侵食率か
ら側壁方向と床方向の熱流束が推定されている。CCI-2 実験（石灰岩系コンクリート）の場合，側壁方向と床方向は同程度と推定されているのに対し，CCI-3 実験（玄武岩系
コンクリート）の場合，側壁方向の熱流束は，床方向より 4 倍程度大きいと推定されている。
④ クラスト破損実験
本実験は上記①③の実験に付随して実施されており，①の SSWICS 実験で形成されたクラストの冷却後破損強度を調べる実験と，③の 2 次元侵食実験の間に高温状態で破損させる実験の 2 種の実験が行われている。
前者の実験では，クラスト浸水実験で形成された直径約 30cm のクラストに荷重を加えて，クラストの強度を測定している。低温クラスト破損実験装置の概要を図付 2-18に示す。クラスト上面の中央部に集中荷重を与え，クラスト破損時の荷重からクラス
トの機械的強度（破損時にクラストに発生する最大の応力）を評価している。クラス
ト上に堆積する冷却水層及び噴出デブリの荷重に起因してクラスト内に発生する応力と，本実験から得られるクラストの機械的強度を比較することにより，実機規模においては，クラストは破損し下部の溶融物と接触する可能性が高いことが示されている。
後者の 2 次元侵食実験では，注水前の高温状態で突き棒（lance）によりクラストが破損する荷重を測定し，注水後にも再度突き棒によりクラストが破損する荷重を測定している。これにより，現実に近い状態でのクエンチ前後のクラスト強度が確認され，クエンチ後は，クラスト強度が大きく低下することが確認された。
OECD/MCCI 実験結果から以下の知見が得られている。
・ガス噴出は，安定なクラストの形成を阻害し，熱伝達を促進する。
・クラストが冷却されると，亀裂が生じクラスト強度は大幅に低下する。
・クラストの亀裂により冷却水が内部に浸入し得る。
・溶融物の噴出により，多孔質の粒子状ベッドが形成される。
・コンクリートタイプにより侵食異方性がみられる。
これらの知見は，CORQUENCH コード等に反映される予定である。
また，CCI 実験について， MAAP コードによりベンチマーク解析が行われている[35]。

5-3-6060

CCI-2 実験（石灰岩系コンクリート）のベンチマーク解析では，固化効果項の影響を除いて対流熱伝達係数のサーベイが可能なように，対流熱伝達係数として床方向と側壁方向ともに 300W/m2-K，固化効果項の指数 n に 0 が設定されている。ここで 300W/m2-K は，実際のMCCI 時の対流熱伝達係数は，固化割合を考慮すると，完全液相状態の値よりも１桁程度小さい値（図付 2-19 において固化割合が 60%程度の値）になることを考慮して設定されている。また，CCI-2 実験では，床方向と側壁方向の熱流束が同程度と推定されていることから，
対流熱伝達係数も同じ値が設定されている。この設定を用いた解析により，コリウム温度の変化（図付 2-20 参照）は，実験データと良く合っている。また，床方向の侵食量は（図付 2-21 参照），実験では注水前には 25cm 程度侵食されているのに対し，MAAP 解析では 33cm
程度の侵食量となっている。側壁方向の侵食量は（図付 2-22 参照），実験では注水前には29cm 程度侵食されているのに対し，MAAP 解析では 33cm 程度の侵食量となっている。
CCI-3 実験（玄武岩系コンクリート）のベンチマーク解析では，対流熱伝達係数として床方向に 80W/m2-K，側壁方向に 300W/m2-K，n に 0 が設定されている。ここで，80W/m2-K は，
CCI-3 実験では，床方向の熱流束が側壁方向の 1/4 程度と推定されていることから，床方向の対流熱伝達係数を側壁方向の 1/4 程度として設定されている。この設定により，コリウ
ム温度の変化（図付 2-23 参照）は，より実験データに近くなっている。また，床方向の侵食量は（図付 2-24 参照），実験では注水前には 5cm 程度侵食されているのに対し，MAAP 解析では 5.7cm 程度の侵食量となっている。側壁方向の侵食量は（図付 2-25 参照），実験では注水前には 29cm 程度侵食されているのに対し，MAAP 解析では 27cm 程度の侵食量となっている。
CCI 実験を対象にした MAAP ベンチマーク解析により，対流熱伝達係数を適切に設定することにより侵食量が概ね合うことが確認されている。また，侵食の異方性を扱うためには，
床方向と側壁方向の対流熱伝達係数の比率の設定が重要であることが示されている。

5-3-61

61

図付 2-9 想定されているデブリ冷却メカニズム[38]

図付 2-10 想定されている火山型クラストの模式図[12]5-3-6262

表付 2-2

SSWICS 実験（デブリ冷却実験）条件一覧[12]図付 2-11 小規模浸水・クラスト強度実験（SSWICS）装置[12]5-3-6363

図付 2-12 SSWICS 実験データと Lister/Epstein モデルのドライアウト熱流束の比較[12]表付 2-3

CCI 実験（2 次元コンクリート侵食実験）条件一覧[12]5-3-6464

表付 2-4

表付 2-5

CCI 実験に使用されたコンクリートの組成 [12]CCI 実験の側壁方向／床方向の侵食率と熱流束の推定値 [12]5-3-6565

図付 2-13 CCI 実験のテスト部[12]図付 2-14 長期 2 次元侵食実験におけるクラスト破壊／強度測定装置[41]5-3-6666

図付 2-15 CCI 実験におけるデブリ-水間の熱流束 [12]図付 2-16 CCI 実験における溶融物平均温度 [12]5-3-6767

図付 2-17 CCI 実験における侵食状態の模式図 [12](a) CCI-1, (b) CCI-2, (c) CCI-3

5-3-6868

図付 2-18 クラスト浸水実験の実験後の低温クラスト強度測定装置[38]5-3-6969本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

図付 2-19 溶融プールからクラストへの対流熱伝達係数

図付 2-20 CCI-2 実験のコリウム温度の実験データと MAAP 解析結果の比較[35]（MAAP New 熱伝達係数（破線）：床方向 300W/m2-K，側壁方向 300W/m2-K，n=0）

5-3-7070

図付 2-21 CCI-2 実験の床方向侵食の実験データと MAAP 解析結果の比較[35]

図付 2-22 CCI-2 実験の側壁方向侵食の実験データと MAAP 解析結果の比較[35]5-3-71

71

図付 2-23 CCI-3 実験のコリウム温度の実験データと MAAP 解析結果の比較[35]（MAAP New 熱伝達係数（破線）：床方向 80W/m2-K，側壁方向 300W/m2-K，n=0）

図付 2-24 CCI-3 実験の床方向侵食の実験データと MAAP 解析結果の比較[35]5-3-7272

図付 2-25 CCI-3 実験の側壁向侵食の実験データと MAAP 解析結果の比較[35]5-3-7373

付録３ 溶融物の水中落下時の粒子化実験と解析モデル
溶融物が水中に落下した場合の粒子化現象については，主に溶融炉心と冷却材の相互作用（FCI）を調べるために行われてきた。ここでは，以下の実験についての知見を述べる。
① FARO 実験（イタリア JRC-ISPRA）② COTELS-FCI 実験（日本 NUPEC）③ DEFOR 実験（スウェーデン KTH）
また，MAAP コードでは，溶融炉心が圧力容器下部プレナム水等に落下する場合の粒子化割合を評価しており，Ricou-Spalding 相関式又は Saito の相関式が用いられている。ここではこれらの相関式の元になった実験についても説明する。
(1)

FARO 実験（イタリア JRC-ISPRA）[1][14]

主に圧力容器内で溶融物が炉心から下部プレナムの水プールに落下した場合の水蒸気爆
発発生の可能性を調べることを目的とした実験で，当初は高圧での実験を対象としてきたが，後に低圧での実験が３ケース実施された。高圧実験装置の構成は，FARO るつぼと実験
容器 TERMOS（内径；710mm，体積；1.5m3）である（図付 3-1 参照）。低圧実験装置では FATと呼ばれる実験容器の中に TERMOS と同一内径の円筒容器が置かれ，観測窓が設置された。
実験手順は原則的に高圧実験も低圧実験も同様で，FARO るつぼで生成された溶融コリウム
は一度リリースベッセルに保持し，その底部にあるフラップ又は弁を開放して溶融物を水プールに落下させる。高圧実験ではリリースベッセルは溶融物落下速度を調整するために加圧可能である。ほとんどの実験は酸化物コリウム（80wt％UO2＋20wt％ZrO2）で実施され，テスト L11 のみが金属 Zr を含むコリウム（77wt％UO2＋19wt％ZrO2＋4wt％Zr）が用いられた。また，多くの実験は高圧・飽和水条件で実施されたが，近年は低圧・サブクール水条
件でも３ケース実施されてきている。主な実験パラメータは，初期圧力（0.2～6MPa），コ
リウム質量（18～177kg），水深（87cm～2.05m），サブクール度（飽和～124K），コリウム組
成（１ケースのみ金属 Zr を 4.1wt%追加）である（表付 3-1 参照）。
実験からの主な知見は以下のとおりである。
・ 高圧・飽和水実験，低圧・サブクール実験の何れにおいても水蒸気爆発は発生していなかった。
・ いずれの実験においても，コリウムの一部が粒子化し，一部はパンケーキ状としてプール底部に堆積した。プール水中に設置した熱電対の測定結果より，一部のコリウムはジェット状でプール底面に衝突したと推定されている（図付 3-2 参照）。
・ 粒子化割合は，ほとんどの実験でプール水深に依存し，高圧飽和水実験の水深 1m 程度で落下コリウムの約 50%が粒子化したが，コリウムに少量の金属 Zr（4.1wt%）を添加した場合と低圧のサブクール水の実験では 100％近くが粒子化した（図付 3-3 参照）。
・ 粒子の質量中央径は比較的大きく，3.2～4.8mm の範囲であり，実験パラメータ（初期圧力，水深，コリウム落下速度，サブクール度）に依存していない。

5-3-7474

表付 3-1

図付 3-1

FARO 実験マトリックス[14]FARO 実験装置（L-27 実験）[14]

5-3-7575

図付 3-2

図付 3-3

FARO L-24 実験後のデブリの堆積状態（数値単位 mm）[14]

FARO/ L-33 実験後のデブリ取り出し段階毎のデブリ状態 [14]（低圧，サブクール水）

5-3-7676

(2)COTELS-FCI 実験（日本 NUPEC）[1][15]

本実験は，圧力容器外での FCI 挙動を調べるために実施された。UO2 混合物を溶融物として用い，その成分は炉内構造物のみならず，圧力容器下部ヘッド内の下部構造物も考慮して設定されている。図付 3-4 に COTELS-FCI 実験装置の概要を示す。実験装置は，UO2 混合物を溶融させる電気炉とその下に溶融物と水との相互作用を調べるメルトレシーバより構成
されている。電気炉は，グラファイト製るつぼと高周波誘導加熱コイルにより構成されている。グラファイト製るつぼは最高温度 3200K で，60kg の UO2 混合物を加熱溶融可能である。
るつぼ下部には，溶融物を落下させるための穴を設けるプラグ破壊装置が設置され，UO2 混合物が溶融後，瞬時に直径５cm の穴が作ることが可能である。メルトレシーバは，高さ約2.5ｍ，外径約 90cm で円筒状をしており，メルトレシーバの中には，高さが可変のメルト
キャッチャーが設置されており，メルトキャッチャーの高さを変化させることによりプール水深を変化させる実験を実施しており，雰囲気体積は変化させないように考慮されている。メルトキャッチャーの中心部は，実炉の格納容器床面を模擬するために厚さ３cm のコ
ンクリートが設置され，溶融コリウム／水／コンクリート間の相互作用を調べることが可能となっている。なお，コンクリートの成分は日本と同じ玄武岩系である。実験マトリッ
クスと主な結果を表付 3-2 に示す。
表付 3-2COTELS-FCI 実験マトリックス[1]5-3-7777

図付 3-4(3)

COTELS-FCI 実験装置[15]DEFOR 実験（スウェーデン KTH）[ 16]

水張り時のデブリ冷却性を調べる目的 で実施されている。溶融物の粒子化 について，水サブクール，水深，メルト成分，過熱度等 の影響を調べるために，溶融物 の落下実験が実施されている。実験装置を図付 3-5 に示示すが，上から，誘導加熱炉，メルト 容器，冷却タ

ンクより成り，冷却タンクはガラス
冷却タンクはガラス張りで 粒子化の観察が可能となっている
可能 となっている。溶融物はメ
容器下部のノズルからのノズルから重力落下
重力落下によりにより水中に落下する。DEFOR-E 実験の 実験マトリック
ルト容器下部
ケースの溶融物は CaO-B2O3（密度 2500kg/m3）の 酸化物が使用さスを表付 3-3 に示すが，6 ケースの溶融物れており，1 ケースはより重量
ケースはより重量の大きい WO3-CaO（密度 5000～8000kg/m3） が使用されている。

実験結果より，粒子化挙動，粒子化挙動は，水サブクールとサブクールと水深水深の影響が大きいとしている 。実験後のデブリ状態例を図付 3-6 に示すが，平均
平均ポロシティは
ポロシティは 0.55～0.7 程度と 大きい値が得られている。主な結果を表付 3-3 に併せて 示す。

5-3-7878

表付 3-3

DEFOR-E 実験マトリックス[16]図付 3-5

DEFOR-E 実験装置[16]5-3-7979

図付 3-6

DEFOR-E 実験後のデブリ状態[16]5-3-8080

(4) 溶融物落下時の水中での粒子化割合の解析モデル
水中落下時の溶融物の粒子化現象について，これまで種々の実験的研究が行われており，また粒子化割合については相関式も提案されている。ここでは，粒子化現象に関する相関式を用いて，溶融炉心ジェットのペデスタル水プールの条件に対する粒子化割合を評価した。相関式として，MAAP コードで使用されている以下の２つの式を用いて計算を行った。
Ricou-Spalding の相関式
Saito の相関式a) Ricou-Spalding の粒子化相関式[25][34]
Ricou-Spalding の相関式は，空気中へガスジェットを噴出させたエントレインメ
ント実験から得られたもので（実験装置は図付 3-7 参照），周囲の流体（空気）がガ
スジェットによってエントレインされ（巻き込まれ），ジェット質量流量が距離とともに増加することを相関づけた整理式である。式の検証例を図付 3-8 に示す。
これをメルトジェットのブレークアップに適用するため，周囲の流体（水）がメ
ルトジェットにエントレインされると考え，そのエントレインされた分だけメルト
ジェットが侵食され粒子化するとして，円柱ジェットのエントレインメント速度（m/s）を，次式で与えている（図付 3-9 参照）。
1ρ  2
ment = Eo  w  u djρ 
 dj 
ここで， E o ：エントレインメント係数

u dj ：ジェット速度 [m/s]ρ dj ：ジェット密度 [kg/m3]
ρ w ：水密度 [kg/m3]

5-3-81

81

図付 3-7

Ricou-Spalding のエントレイメント実験装置[25]

図付 3-8

Ricou-Spalding の相関式の検証例[25]5-3-8282

図付 3-9 溶融デブリのジェットブレークアップとクエンチの
溶融デブリのジェットブレークアップとクエンチの 概念[34]b) Saito の相関式[26]フレオン 11 あるいは液体窒素の高温水ジェットを吹き込む実験 により，ジェの中に高温水ジェットをットのブレークアップ長さに関する 相関式を導出している。実験装置 を図付 3-10 に示す。実験条件は，ノズル直径 5～ 40mm，ジェット速度 2.5～15m/s である 。実験結果
は，図付 3-11 に纏められ，それに 基づき下記相関式が導かれた。

2.1

ここで，

.

.

：フルード 数

L：ジェットのブレ ークアップ長さ [m]：ジェット 直径 [m]：ジェット 密度 [kg/m3]：冷却材密度 [kg/m3]：ジェット 速度 [m/s]：重力加速度 [m/s2]

5-3-8383

図付 3-10 Saito らの実験装置[33]

図付 3-11 Saito らの実験結果[33]5-3-8484本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

c) 相関式の適用性に関する検討[33]Ricou-Spalding の相関式は，ガスジェットの空気中でのガス巻き込み現象に基づき作成された相関式であるため，溶融炉心物質ジェットの水中での粒子化現象への適用性の確認のため，ALPHA-MJB 実験，FARO-TERMOS 実験を対象に比較計算された結果を示す。
合わせて Saito の相関式との比較も示す。Saito の相関式ではジェットのブレーク
アップ長さが与えられるので，ジェットの突入から消滅までジェット径が直線的に減少すると仮定して粒子化割合が求められている。
(a)ALPHA-MJB 実験
水中に鉛－ビスマスの溶融物を落下させた実験であり，密度や速度といった計
算条件は表付 3-4 の値が用いられている。実験結果と計算結果の比較を図付 3-12，
13 に示す。実験結果は，ジェットブレークアップ長さとして記載されている水深を，粒子化割合 100%として記す。
Ricou-Spalding 相関式における係数 E_0 の値は，MAAP 推奨値

を使用する

ことで，実験結果及び Saito の相関式とよく一致する結果が得られている。
表付 3-4
ケース
MJB01
MJB02

溶融物密度
(kg/m3)10453↑

ALPHA-MJB 実験の検証解析条件[33]水密度(kg/m3)
958
↑

ジェット径
(m)
0.03
0.021

（補足）溶融物密度は鉛－ビスマス 500K での値を使用

図付 3-12 ALPHA-MJB01 実験結果との比較[33]5-3-8585

ジェット速度
(m/s)
3.05
2.86

本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

図付 3-13 ALPHA-MJB02 実験結果との比較[33](b)FARO-TERMOS 実験
水中に実機燃料組成を模擬した溶融物を落下させた実験であり，密度や速度等の条件は表付 3-5 の値が用いられている。実験結果と計算結果の比較を図付 3-14
～19 に示す。実験結果は水プール水深に対する粒子化割合として示されている。
計算結果は，実験結果とよく合っているか粒子化割合を小さく評価している。な
お，Ricou-Spalding 相関式における係数 E_0 の値は
表付 3-5が使用されている。

FARO-TERMOS 実験の検証解析条件 [33]5-3-8686本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

図付 3-14 FARO-L06 実験結果との比較[33]

図付 3-15 FARO-L08 実験結果との比較[33]

図付 3-16 FARO-L11 実験結果との比較[33]5-3-8787本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

図付 3-17 FARO-L14 実験結果との比較[33]

図付 3-18 FARO-L19 実験結果との比較[33]

図付 3-19 FARO-L20 実験結果との比較[33]5-3-8888本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

d) 格納容器ペデスタル水張り水深に対する粒子化割合の評価[33]Ricou-Spalding の相関式及び Saito の相関式に基づき，大気圧飽和状態の水プールを想定し，水深とジェット直径及びジェット速度をパラメータとして粒子化割合を整理されたものを図付 3-20～22 に示す。

図付 3-20 Ricou-Spalding 相関式による粒子化割合のマップ[33]

図付 3-21 Saito の相関式による粒子化割合のマップ（Vj=5.0m/s）[33]

5-3-8989本資料のうち，枠囲みの内容は商業機密に属しますので公開できません。

図付 3-22 Saito の相関式による粒子化割合のマップ（Vj=10.0m/s）[33]

5-3-9090

付録４ 溶融物の拡がり実験
溶融物の拡がり挙動に関して，表付 4-1 に示す実験的研究が実施されている。この中で水が存在する場合の溶融物の拡がりに関する主要な知見について以下に整理する。
表付 4-1 溶融物拡がりに関する主な実験[27]実験名
BNL[18]SPREAD[19]CORINE[36]
VULCANO
KATS[23]
COMAS
ISPRA
S3E[27]
PULiMS[21]

実施国
米国
日本
フランス
フランス
ドイツ
ドイツ
EU JRAC
スウェーデン

拡がり形状
2D
1D, 2D
2D
2D
1D, 2D
1D, 2D2D
1D, 2D
2D

溶融物

その他
Pb
ドライ，ウエットステンレススチール
ドライ，ウエット
グリセロール他
ドライ，ウエット
コリウム
ドライ
テルミット
ドライ，ウエット
コリウム
ドライ
コリウム
ドライ
模擬物（ウッズメタル等） ドライ，ウエット
模擬物(B2O3-WO3 等)
ウエット(1) Mark-I 型格納容器のシェルアタックに関する研究[17][18][19][20]Mark-I 型格納容器のシェルアタックに関して ROAAM 手法による確率論的評価[17]が行われており，その中で水張り時のデブリの拡がりの評価が行われている。溶融物拡がり挙動解析には，MELTSPREAD コードが使用されている。溶融物の拡がり挙動には，溶融物落下流量，溶融物成分，溶融物過熱度が主に影響するとされている。
米国 BNL では，溶融物として鉛を使用した小規模な拡がり実験を，ドライ条件とウ
ェット条件で行っている[18]。溶融物の過熱度と水深が，デブリ拡がり挙動において重要であること，溶融物の拡がりは，水への熱伝達による溶融物の温度挙動と，溶融物の流体力学的挙動により制限されることが示されている。
国内 BWR 産業界では，SPREAD 実験[19][20]が実施されている。溶融物として約 2000K に加熱したステンレススチールを用い，溶融物落下流量，水深，溶融物過熱度等を変化させて，格納容器模擬試験体 (1 次元形状と 2 次元形状の 2 種類，2 次元形状は図付
4-1 を参照)に放出し，溶融物の拡がり挙動に関する基礎的なデータを収集した。実験
条件を表付 4-2 に示す。また，溶融物の拡がりを解析する SPREAD コードが開発されている。溶融物が重力で拡がると仮定し，熱源として崩壊熱及び化学反応を考慮し，床コンクリートや上面の水プールへの伝熱を考慮し，質量，運動量及びエネルギのバランスを解くことにより，溶融物の移動速度と溶融物厚さを求めている。
（図付 4-2 参照)
このモデルでは，流動限界固相率を超えた時点で流動が停止すると仮定されており，
実験データを再現する流動限界固相率として 55%が使用されている。2 次元体系に対する実験結果と解析結果の比較を図付 4-3 に示すが，よい一致が示されている。

5-3-9191

図付 4-1

溶融物拡がり実験の試験体（2 次元形状）[20]

図付 4-2SPREAD コードの概念図 [20]5-3-9292

図付 4-3

溶融物拡がり距離の実験値と解析値の比較 [20]
SPREAD 実験条件 [19]NozzleSlit Width
Floor
Diam.
(m)
Material(mm)
0.15
Concrete
70
↑
↑
↑
30
↑
↑
0.05
↑
↑
↑
↑
↑
↑
↑
↑
↑
↑
↑
↑
↑
↑
0.025
↑
↑
0.05
MgO
↑
Concrete
↑
↑
↑
↑
↑
0.025
↑
↑
0.05
70
↑
30
↑
↑
↑
↑
↑
↑
↑
↑
↑
↑
↑
0.025
↑
↑
0.05
↑
↑
↑
↑
↑
↑
↑
↑
0.025
↑
↑
0.05
↑
↑
↑
↑
↑

表付 4-2Test
No.

Melt Mass
(kg)

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25

95.0
7.5
24.0
43.0
19.0
66.0
115.0
66.0
94.2
60.2
93.4
67.0
93.1
68.1
63.6
116.0
106.2
108.8
61.9
64.0
69.1
65.0
60.0
65.0
63.8

Super
Heat
(K)
133
77
149
115
110
107
135
112
110
86
122
115
116
130
133
131
171
158
108
91
107
90
94
125
94

* O/C: Outside of Cylindrical Region

I/C: Inside of Cylindrical Region

5-3-93

93

Water
Depth
(mm)
0
0
0
0
0
40
0
0
0
0
40
40
40
0
0
0
0
0
0
10
40
70
40
40
Spray

Water
Added
Position *
O/C
O/C
O/C
O/C
O/C
O/C
O/C
O/C
O/C,I/CO/C

(2)EPR コアキャッチャの溶融物拡がりに関わる実験[23][29][30][36]
CEA/DRN/DTP で行われた CORINE 実験[36]では，低融点物質（グリセロール他）を模擬物質として使用して，水中での拡がり挙動を調べる実験が実施され，拡がり先端の移動速度や底部に形成されるクラストの影響が調べられた。
独カールスルーエ研究センター（FZK）で実施された KATS 実験[23][29][30]では，溶融物としてテルミット（Al2O3 約 150kg，Fe 約 150kg）が使用され，溶融物の放出速度や温度，拡がり形状（1D,2D），床の材質（コンクリート，セラミック，コーティング），水の有無をパラメータに溶融物の拡がり実験が行われている。実験装置を図付 4-4 及び図付 4-5 に示す。Al2O3 と Fe では密度が異なり成層化するため，溶融物の出口を 2 箇所設け，最初に Al2O3 が放出し，最後に Fe を放出することにより酸化物溶融物の拡がりと金属溶融物の拡がりを分けて実験が可能となっている。実験条件を表付 4-3 に示す。KATS-10 と KATS-11 の実験条件はほぼ同様であるが，KATS-10 の方は 1mm の水張りをしてあり，KATS-11 の方はドライ条件となっている。両者の拡がり結果を図付 4-6 に示すが，両ケースのように溶融物の放出速度が比較的高い場合は，冷却材の有無によらず同様な拡がり挙動になる結果となっている。

図付 4-4

KATS 実験装置の概要図 [29]5-3-9494

(1D)
図付 4-5(2D)
KATS 実験の 1D と 2D の拡がり形状の写真[30]5-3-9595

(3)表付 4-3

KATS 実験条件と拡がり距離（酸化物溶融物，1D 拡がり）[29]

図付 4-6

KATS 実験の水の有無の影響（酸化物溶融物，1D 拡がり）[29]

KTH の溶融物拡がり研究[21][27]

スウェーデン KTH では，溶融物拡がりの小規模実験として，S3E 実験[27]が行われている。
種々の溶融物と冷却材の組合せ（例ウッズメタル－水，塩（NaNO3-KNO3）－パラフィン）

で溶融物温度や溶融物流量等をパラメータに，溶融物の拡がり挙動に関するデータが取られている。溶融物が冷却される過程において，低密度の多孔質の層が形成される結果が得られており，溶融物の流れにより，この多孔質層が持ち上げられ，流量（慣性力）によって決定される距離まで拡がる挙動を示した。
また，水中での溶融物拡がり挙動を調べる PULiMS 実験[21]が実施されている。この実験では，浅い水プール（水深 20cm）へ溶融した Bi2O3-WO3 合金を流入させ，その拡がり挙動が観察された。実験装置を図付 4-7 に示し，実験条件を表付 4-4 に示す。水中へ流入した溶融

5-3-9696

物の拡がりの様子を図付 4-8 に示す。また，固化したデブリ層の性状を図付 4-9 に示す。
固化したデブリ層は３層を形成し，デブリ層下面とコンクリートは密着せずギャップが存在しており，水又は蒸気が存在していた可能性がある。３層の最下層は薄い 1.5-2mm のケーキ（クラスト）層で空隙は無く，中間層は割れ目や細長い孔が含まれている。最上部は不規則な構造で高多孔性である。

図付 4-7

PULiMS 実験装置の概要図[21]表付 4-4

PULiMS 実験条件 [21]5-3-9797

図付 4-8

PULiMS 実験 E4 ケースの溶融物拡がり挙動のスナップ写真 [21]図付 4-9(4)

PULiMS-E4 実験の固化デブリの性状写真 [21]

粒子状ベッドの拡がり[22]

ANL では，水プール中に不均質に堆積させたデブリベッドが，内部の沸騰により均一化す
る様子を確認するため，水プール中に UO2 ，SUS，Cu の 0.2～1mm 程度の粒子ベッドを非均一の厚さに堆積させ，誘導加熱により崩壊熱発生を模擬させた実験[22]を行っている。図付 4-10 に粒子ベッド厚さの均一化の概念図を示す。非均一の厚さに堆積した粒子ベッドは，

5-3-9898

誘導加熱により粒子ベッド内に沸騰が生じ，粒子が吹き上げられて再堆積する過程で厚さが均一化し，均一化に要した時間は 2～3 分程度であると報告されている。

図付 4-10 セルフレベリング実験におけるデブリベッドの均一化の概念図 [22](5)

溶融物の拡がり解析[20][27][28][39][40]

溶融物の拡がり挙動の解析のために，種々の解析コードが提案されている。主要なもの
を表付 4-5 に示す。これらの解析コードは，実験データを元に検証が行われている。また，
実験と実機条件とを関連づけるスケーリング則の提案も行われている。[20][27][28]
表付 4-5 溶融物拡がり挙動の主な解析コード[27]Code
CORFLOW *CROCOLAVA *MELTSPREAD *MECORASPLAV / SPREADSPREAD
THEMA *
Footnote

Country
Model characteristics
Germany
2D - mechanistic
FranceSimplified (ISPN)Japan
2D - Bingham fluidUSA
1D - complex physics and chemistry
Germany
Mechanistic (Ruhr Univ.)
Russia2D - hydrodynamicJapan
simplified
Francephenomenological (CEA/DRN)* Code receiving International recognition

5-3-9999

a)

3 次元溶融物拡がり解析例
溶融物の拡がり挙動を 3 次元で評価した例もある。図付 4-11 に示すとおり過酷事故

解析コード SAMPSON のデブリ拡がりモジュール(DSA) [39]を使用して，ABWR の格納容器
ペデスタル上におけるデブリ拡がり評価を実施している。[40]床面に 200mm の水位がある条件(Wet 床条件)で，制御棒駆動機構ハウジング 1 本の破損を想定して、全炉心相当の溶融炉心がペデスタル（半径約 5.3m）へ落下する場合の溶融炉心の床上の拡がり面積を評価している。床面に水が存在する場合でも，溶融炉心は落下位置に依らずペデ
スタルのほぼ全面に拡がり，最終的な堆積高さは 400～500mm 程度となる結果を得ている。

図付 4-11

ABWR を対象にした 3 次元溶融炉心拡がり評価結果例(SAMPSON コード) [40]b) スケーリング則を用いた溶融物拡がり相関式(3)に示した PULiMS 実験の結果を元に，水中での溶融物の拡がり距離を求めるスケ
ーリング則が提案されている[28]。これは，障害物の無い水平な床面に落下した溶融物の拡がり停止後の平均厚さや拡がり距離を求める相関式である。一般的に，溶融物の拡がり挙動は，溶融物の堆積高さに対応した重力による拡がりと，溶融物冷却に伴う粘性増加による流動性低下の競合する２つのプロセスにより支配される。（図付 4-12
参照）
スケーリング則では，競合する２つのプロセスの時間スケールとして，流体力学的
拡がり時間τconv を，拡がり停止後の平均厚さδv に達するまでの時間と定義し，固化時

5-3-100

100

間τsolid を，溶融物が冷却され固化するまでの時間と定義し，これらの時間の比を，拡/ " #$% ）と定義している。また，拡がり停止後のがりの無次元時間スケール T（

平均厚さδv と表面張力とのバランスで決まる厚さδcap の比を，無次元長さスケール L（& /& '( ）と定義している。

障害物の無い水平な床を想定し，落下した溶融物が軸対称に拡がると仮定すると，無次元長さスケール L は，無次元時間スケール T を用いて(1)式で表されるとしている。
添え字の i は 慣性力支配流れの値を示し，v は粘性力支配流れの値を示している。
)* ∙ , -/. ∙ /-/.

(1)

また，拡がり距離r 1 と拡がり時間 t の関係は，以下のように表されるとしている。
2$ 1

慣性力支配流れの場合

3$ ∙ 1 4/5

2 1

粘性力支配流れの場合(2)

3 ∙ 1 -/.

(3)

ここで，
& '( ：表面張力とバランスする厚さ（

267/ 8 ′）Cv：比例定数。PULiMS 実験より，Cv=1.30 が得られている。
N：粘性因子（τ :;<,<
τ :;<,>

β：拡がり係数

β$

∙

δ<

δ ?@

）

0.84

′D

" #$% ：デブリ固化特性時間(s) （

-/K

-/5

，β

& '( ∙

L@,M ∙∆OPQ@ RS∙TU）
8 ∙ V" RV" YV Z

′H I

0.67 G
Q@

X;

<∙

J

?@

Δ,"[( ：デブリ過熱度(K) ，*(,8 ：デブリ比熱(J/kg/K)\"[( ：デブリ上面の熱流束（放射及び対流熱伝達）(W/m2)\"% ：デブリ下面の熱流束(対流熱伝達) (W/m2)\ ：単位体積当たりの崩壊熱(W/m3)]^ ：デブリの溶融潜熱(J/kg)，η：流動停止までの有効潜熱割合(-)

σ：デブリ表面張力(N/m) ， 8 ：デブリ密度(kg/m3) ， _ ：水の密度(kg/m3)
g：重力加速度(m/s2) ，v：デブリの動粘性係数（m2/s）
g′：水中における重力加速度の補正 （3b b ：デブリ体積(m )（c -YdḾM

g ρ8 − ρ_ /ρ8 ）(m/s2)）， m：デブリ質量(kg)

G：エントレインされなかった連続体デブリの体積流量(m3/s)（
1 l# ：原子炉圧力容器からの溶融炉心放出時間(s)m8 ：デブリ落下時のエントレイン割合(-)

5-3-101

101

fg:g
hijk

）

′8：有効デブリ密度(kg/m3) （ ρ8 1 − ε8 ）， n8：デブリボイド分率(-)(1)式は，拡がり停止時のデブリ平均厚さが，無次元時間スケール（＝溶融物の拡がり時間／固化時間）の平方根に比例することを意味しており，溶融物の拡がりが遅いほど，固化するまでの時間が短いほど，無次元時間は大きくなり，結果として，拡がり停止時の平均厚さが大きく，拡がり面積は小さくなる傾向を示している。逆に，溶融物の拡がりが早いほど，固化までの時間が長いほど，無次元時間は小さくなり，平均厚さが薄くなり，拡がり距離は大きくなる。このスケーリング則について，PULiMS
実験との比較が行われており，図付 4-13 に，実験結果とスケーリング則による評価との比較を示すが，概ね一致している結果が得られている。
, 及び N は次式のように得られる。

上記のスケーリング則の式より，,/

G

g:g ∙pP:k>Xor∙β< ∙Z ?@ ∙L< ∙s

rZ ?@

-/4

J

g: ・

(4)式及び(5)式により
2

,

3 ∙

q/

,

-/.

./4

t
q/

(4)

u> v/I(5)u<, が求まり，次式で拡がり距離が得られる。

(6)

BWR の溶融物条件（表付 4-6 参照）を用いて，上記の PULiMS スケーリング則を適用
し，(6)式により溶融物の拡がり半径を評価した。図付 4-14 に評価結果を示す。制御
棒駆動機構ハウジング破損時の溶融物落下質量速度である約 8500 kg/s の場合，約 18mの拡がり半径を得る。溶融物質量を減らして落下質量速度を減らすと，拡がり距離は減少するが，BWR 条件の落下質量速度の 1/10 にしても，まだ，ペデスタル半径（Mark-I
改の場合約 3.2m）を超えた拡がり半径を得る。
従って，PULiMS スケーリング則を適用しても，BWR の溶融物条件では，溶融物はペデスタル床全面に拡がる結果となる。

5-3-102102

を示
図付 4-12 溶融炉心の拡がりに関 わる熱伝達及び流体力学的な力を示す模式図

図付 4-13 PULiMS 実験を対象にしたスケーリング 則による評価と実験
実験結果結果の比較[28]5-3-103103

本資料のうち，枠囲みの内容は商業機密
商業機密にしますので公開
公開できません
できません。
属しますので表付 4-6

図付 4-14

BWR 溶融物条件（BWR-5, Mark-I 改）

溶融物落下質量速度と拡 がり距離の関係（PULiMS スケーリング
スケーリング則適用
則適用）

5-3-104104

付録５

ベッドの冷却性実験と解析 モデル
粒子状ベッドの冷却性

粒子状
粒子状ベッドのベッドの冷却性については，高速増殖炉
高速増殖炉のの炉心損傷事故
炉心損傷事故の評価においてにおいて重要重要な評行われてきている[32]。また，価となることから，従来から多くの実験的及 び解析的研究が行われてきている
TMI-2 事故以降は，軽水炉を対象とした 研究も進められている。

デブリベッドデブリベッド冷却
冷却に関する研究の主要 な焦点の一つは，デブリベッドのデブリベッドの冷却限界が，崩るかどうかであり，判断基準 のパラメータはデブリベッドのドライアウト
パラメータはデブリベッドのドライアウト 熱流
壊熱を上回るかどうかであり考えられている
えられている。発熱する粒子状ベッドの 内を流下する冷却水と上昇 する蒸気の対向束と考二相流二相流がデブリベッドがデブリベッド内で形成される。粒子状デブリベッドの流動状態の模式図 を図付 5-1
局所的にデブリベッド内
内でドライアウトがでドライアウトが生じ ，デブリベッドに示す。崩壊熱が高いと，局所的にデブリベッド
熱流束がドライアウト
がドライアウト熱流束
熱流束となる。 粒子状デブリベの温度が急激に高くなる。この時の熱流束ッドのドライアウト熱流束は，種々の実験 が行われ，種々の相関式が提案
提案されている
されている。その中でも Lipinski 0-D モデル[24]（次頁参照 ）は，小さい粒径から大きい粒径 まで適用可能とされており，広く使われている。MAAP コード内でも圧力容器下部
でも圧力容器下部プレナムの
プレナムの粒子状
粒子状ベッモデルでは，ドライアウト
ドライアウト熱流束は，主ドの冷却の評価に使用されている。Lipinski 0 -D モデルでは，に粒子径，ポロシティ，堆積高さ，圧力 に依存する。

図付 5-1 粒子状デブリベッドの 流動状態の模式図[32]5-3-105105

Lipinski 0-D モデルについて（日本機械学会編「沸騰熱伝達と冷却」より抜粋）[32]

5-3-106106

参考文献:[1] （財）原子力発電技術機構（NUPEC）,「重要構造物安全評価（原子炉格納容器信頼性実証事業）に関する 総括報告書」 (2003).
[2] EPRI, Technical Foundation of Reactor Safety,Revision 1, Knowledge Base for
Resolving Severe Accident Issues, 1022186 (2010).
[3] B.R.Sehgal, et al., “ACE Project Phases C & D : ACE/MCCI and MACE Tests,”
Proceedings of the U.S. Nuclear Regulatory Commission, 19th Water Reactor
Safety Information Meeting, NUREG/CP-0119, Vol.2 (1991).
[4] E.R.Copus, at al.,“Sustained Uranium-Concrete Interactions (SURC)”, Reactor
Safety Research Semiannual Report, NUREG/CR-4805 (1986).
[5] G.A.Greene, R.A.Bari,”ISP-24 International Standard Problem No 24, SURC-4
ExperimentOnCore-ConcreteInteractions,

Final

Workshop

Summary

Report,”CSNI Report No 155, Volume 2 (1989).
[6] H.Alsmeyer, et al.,”BETA experiments on Zirconium Oxidation and Aerosol
Release during Melt-Concrete interaction,”Proceedings of the Second OECD(NEA)
CSNI Specialist

Meetingon

Molten Core Debris-Concrete Interactions,

NEA/CSNI/R(92)10, Karlsruhe, Germany (1992).
[7] R.E.Blose,etal.,“SWISS:

Sustained

HeatedMetallic

Melt/ConcreteInteractions With Overlying Water Pools,”NUREG/CR-4727 (1987).
[8] R.E.Blose, et al.,“Core-Concrete Interactions with Overlying Water Pools The WETCOR-1 Test,” NUREG/CR-5907 (1993).
[9] M.T.Farmer, et al., "Status of Large Scale MACE Core Coolability Experiments,"
OECD Workshop on Ex-Vessel Debris Coolability, Karlsruhe, Germany (1999).
[10]

M.T.Farmer, et al., ” Corium Coolability under

Ex-Vessel

Accident

Conditions for LWRs,” Nuc. Eng. and Technol., 41, 5 (2009).
[11]

H.Nagasaka, et al.,“COTELS Project (3): Ex-vessel Debris Cooling Tests,”

OECD Workshop on Ex-Vessel Debris Coolability, Karlsruhe, Germany (1999).
[12]

M.T.Farmer, et al.,“OECD MCCI Project Final Report,” OECD/MCCI-2005-TR06(2006).
[13]

M.T.Farmer, et al.,“OECD MCCI-2 Project Final Report,” OECD/MCCI-2010-TR07(2010).
[14]

D.Magallon,” Characteristics of corium debris bed generated in large-scale

fuel-coolant interaction experiments,”Nucl. Eng.Design, 236 1998-2009 (2006).
[15]

M. Kato, H. Nagasaka, et al.,“COTELS Fuel Coolant Interaction Tests under

Ex-Vessel Conditions,” JAERI-Conf 2000-015 (2000).
[16]

A. Karbojian, et al.,” A scoping study of debris bed formation in the DEFOR

5-3-107107

test facility,” Nucl. Eng. Design 239 1653–1659, (2009).
[17]

T.G.Theofanous, et al.,”The Probability of Liner Failure in a Mark-I

Containment,” NUREG/CR-5423,(1991).
[18]

G.A.Greene, et al., "Experimental Studies on Melt Spreading, Bubbling Heat

Transfer and Coolant Layer Boiling," Proceedings of 16th Water Reactor Safety
Information Meeting, NUREG/CP-0096, pp.341-358 (1988).
[19]

H.Suzuki, et al., ”Fundamental Experiment and Analysis for Melt Spreadingon Concrete Floor,” Proceedings of 2nd ASME/JSME Nuclear Engineering
Conference, Vol. 1, pp. 403-407 (1993).
[20]（社）日本原子力学会，シビアアクシデント熱流動現象評価 (2001)[21]

A.Konovalenko, et al.,”Experimental Results on Pouring and Underwater

Liquid Melt Spreading and Energetic Melt-Coolant Interaction,” Proceedings
of NUTHOS-9, N9P0303, Taiwan (2012).
[22]

J. D. Gabor, L. Baker, Jr., and J. C. Cassulo, (ANL), “Studies on Heat Removal

and

Bed

LevelingofInduction-heatedMaterials

Simulating

Fuel

Debris,”SAND76-9008 (1976).
[23]

Proceedings of the Second OECD(NEA) CSNI Specialist Meeting on Molten Core

Debris-Concrete Interactions, NEA/CSNI/R(92)10,Karlsruhe, Germany (1992).
[24]

R.J.Lipinski,”A

Model

for

BoilingandDryoutinPartic1e

Beds,“ NUREG/CR-2646, SAND82-0765 (1982).
[25]

F.P.Ricou, D.B.Spalding, “ Measurements of Entrainment by Axisymmetrical

Turbulent Jets,”Journal of Fluid Mechanics, Vol.11, pp.21-32 (1961).
[26]

M.Saito, K.Sato, S.Imahori,”Experimental Study on Penetration Behaviors of

Water Jet into Freon-11 and Liquid Nitrogen,” ANS Proc. of National Heat
Transfer Conference, Houston, Texas, USA, (1988).
[27]

T.N.Dinh, et al.,”Core Melt Spreading on a Reactor Containment Floor,”

Progress in Nuclear Energy，Vol. 36，No.
[28]

4，pp.405-468 (2000).

A.Konovalenko, P.Kudinov,”Development of Scaling Approach for Predictionof Terminal Spread Thickness of Melt Poured into a Pool of Water,” Proceedings
of NUTHOS-9, N9P0302, Taiwan (2012).
[29]

B.Eppinger, et al.，”KATS Experiments to Simulate Corium Spreading in the

EPR Core Catcher Concept,”FzK, Karlsruhe,Germany.
[30]

B.Eppinger, et al.，” Simulationsexperimente zum Ausbreitungsverhalten vonKernschmelzen: KATS-8 bis KATS-17,” FZKA 6589 (2001).
[31]

I.Lindholm, “A Review of Dryout Heat Fluxes and Coolability of Particle Beds,”

SKI Report 02:17 (2002).

5-3-108

108

[32]

日本機械学会編，「沸騰熱伝達と冷却: 第 8 章 デブリベッドの冷却」，日本工業

出版，平成元年[33]（財）原子力安全研究協会，「シビアアクシデント対策評価のための格納容器イベ

ントツリーに関する検討」，平成 13 年 7 月
付録 7-4 溶融炉心の粒子化割合，付録 7-5 粒子状物質ベッドの冷却性
付録 7-6 粒子状物質の粒子径[34]EPRI/FAI, MAAP4 Computer Code Manual (2007).

[35]

Q.Zhou, et al.,”Benchmark of MCCI Model in MAAP5.02 against OECD CCI

Experiment Series,” Proceedings of ICAPP 2014, ICAPP-14352, Charlotte,USA,(2014).
[36]

J.M.Veteau and R.Wittmaack.，"CORINE Experiments and Theoretical Modeling,"

Proceedings of FISA-95， Luxemburg EUR 16896 EN， pp.271-285 (1996).
[37]

R.K.Cole, et al.,”CORCON-Mod2: A Computer Program for Analysis of Molten

Core-ConcreteInteractions,”NUREG/CR-3920,SAND84-1246,

SandiaNational

Laboratories, Albuquerque, NM (1984).
[38]

独立行政法人原子力安全基盤機構，「平成１７年度 原子炉施設のアクシデントマ

ネージメントに係る知識ベースの整備に関する報告書 ＝環境への影響緩和（デブリ冷却）＝」，06 基シ報-0003，平成１８年７月[39]

M.Hidaka, et al.,”Verification for Flow Analysis Capability in the Modelof Three-Dimensional Natural Convection with Simultaneous Spreading, Melting
and Solidification for the Debris Coolability Analysis Module in the Severe
Accident Analysis Code ‘SAMPSON’, (II),”Journal of NUCLEAR SCIENCE andTECHNOLOGY, Vol. 39, No. 5, p. 520–530 (2002).
[40]

中島 他，SAMPSON コードによる ABWR 格納容器ペデスタル上の炉心デブリの 3

次元拡がり評価，日本原子力学会「2013 年秋の大会」,H12，2013 年 9 月[41]

M.T.Farmer, et al.,“OECD MCCI Project 2-D Core Concrete Interaction (CCI)

Tests: Final Report,” OECD/MCCI-2005-TR05 (2006).

5-3-109109