First results and next step of experimental investigation on melt spreading on top of the shield plate penetrated by a forest of CRGTs and IGTs�
Di Fang
Royal Institute of Technology (KTH)
Division of Nuclear Science and Engineering (NSE)
Half-time webinar for APRI-12, June 4, 2025
Ex-vessel melt behavior
2
Illustration of corium injection into the drywell
Lessons from Fukushima accident
3
[1] Farmer, M. T. (2018). The MELTSPREAD Code for Modeling of Ex-Vessel Core Debris Spreading Behavior (ANL--18/30, 1483992; p. ANL--18/30, 1483992)
[2] Yamashita, T., Sato, T., Madokoro, H., & Nagae, Y. (2022). BWR lower head penetration failure test focusing on eutectic melting. Annals of Nuclear Energy, 173, 109129.
Reactor vessel head
Horizontal Beams
Control rod drive housings (185)
Hanger rods
Horizontal support bars, grid plates and grid clamp
Instrument tube TYP. (55)
Railings
Platform I-beams (4)
Wheels
Doorway
Sump
Water level
CRD flanges
Illustration of Below Vessel Structure in a Mark I containment [1]
Illustration of Below Vessel Structure in a Mark I System [2]
Pedestal wall
Vessel support skirt
4
Depiction of corium interaction with
below vessel structure in a Mark I containment
CRD housing brackets now (Unit 3)
2017.7.24 Unit 3 PCV internal investigation (Preliminary report of July 22 investigation). 2017. [video] TEPCO: TEPCO
Lessons from Fukushima accident (contd.)
CRD housing brackets�before the accident (Unit 3)
5
Concerned problems in the Nordic BWR
Simplified into
BWR 75 – reactor containment [1]
[1] Nilsson, L. (2006). Development of an Input Model to MELCOR 1.8.5 for Oskarshamn 3 BWR.
Approach for melt-structure interaction investigation
6
Nordic BWR
DEFOR-S
Simulation
MPS or MELTSPREAD
Experiment
MELCOR
FDNPP
COSIN
Melt spreading
Structure melting
Melt distribution
Dynamic characteristics
Debris final configuration
Heat transfer characteristics
Extended prediction, guideline
Boundary condition, database
Scaling analysis
focused phenomenon
Simplified structures design
Prototypical prediction
Test separate effect
Test integral effect
Simplified scaling analysis
7
Jet releasing
Melt spreading
Jet impingement
Jet along structures
BWR
Bernoulli's equation
Dinh’s spreading model
Model of liquid metal flow in bundles
Jet impingement model
Vj , Hfree
DS, G, Hs, Dr
T, L
Hs, Lpicth, Dd
Radiation shield
Simplified scaling analysis
8
Parameters | Prototypical | Scaling ratios | Test | Selection |
Pressure (bar) | 2 - 5 | - | 1 | 1 |
Jet diameter from outlet (mm) | 70 -223 | 1:6 | 13 - 37 | 20 |
Melt mass (kg) | 64600 | 1:529 | 122 | 20 |
Releasing time of metallic melt (s) | 181 | 1:6 | 30 | 22 |
Velocity at horizontal substrate (m/s) | 4.26 -8.31 | 1:2.4 | 1.74 - 3.39 | 2.5 |
Structure height (m) | 0.475 | 1:6 | 0.08 | 0.08 |
Structure diameter or CRD housing (mm) | 190 | 1:6 | 33 | 33 |
Distance between Radiation protection and damper (m) | 4.5 | 1:6 | 0.75 | 0.75 |
Damper height (mm) | 122 | 1:6 | 20 | 20 |
Damper Diameter (mm) | 265 | 1:6 | 44 | 44 |
Pitch diameter ratio of structures | 1.58 | 1:1 | 1.58 | 1.58 |
Jet and structure diameter ratios | 0.37 – 1.23 | 1:1 | 0.37 – 1.23 | 0.6,1 |
Test facility - COSIN
9
Induction heating furnace for melt preparation
Funnel
Spreading plate (quartz)
1000 x 1000 mm
Structures
Nozzle
TCs | Purpose |
Thermocouples (1.5 mm above substrate) | Boundary layer temperature on the top surface of substrate |
Thermocouples (5 mm above substrate) | Melt top surface temperature |
Thermocouples insides structures (1.5 mm above substrate) | Heat transfer through structures |
Thermocouples at different height | Melt temperature at different height |
Instrumentation details
Test matrix-melt spreading with structures (COSIN)
10
Experiment | Simulant | Spreading plate | Melt mass (kg) | Melt superheat (℃) | Velocity (m/s) | Free fall height (mm) | Jet diameter (mm) | With structures | Structure diameter (mm) | Structure material |
COSIN-S1 | Zn | Quartz | 20 | 102 | 2 | 200 | 10 | No | - | - |
COSIN-S2 | Zn | 76 | 2 | 200 | 20 | No | - | - | ||
COSIN-S3 | Zn | Quartz | 80 | 2 | 200 | 20 | Yes | 33 | SS | |
COSIN-S4 | Zn | SS | 20 | 80 | 2 | 200 | 20 | Yes | 33 | SS |
COSIN-S5 | Zn | Quartz | 80 | 2.6 | 400 | 20 | Yes | 33 | SS | |
COSIN-S6 | Zn | 80 | 2 | 200 | 30 | Yes | 33 | SS | ||
COSIN-S7 | Bi2O3-WO3 | 80 | 2 | 200 | 20 | Yes | 33 | SS | ||
COSIN-S8 | Zn | 80 | 2 | 200 | 20 | Yes | 33 | Zn or Sn |
The new test: COSIN-S3
11
COSIN-S3: Visualization (compared with COSIN-S2)
12
COSIN-S2
COSIN-S3
Tsup=76 ℃; Hfree=200 mm
Without structures
Tsup=80 ℃; Hfree=200 mm
With structures
0.2X
0.2X
COSIN-S3: Final configuration (compared with COSIN-S2)
13
Melt spreading profile
COSIN-S2
Melt thickness
Cavity formation
COSIN-S3
Conclusion
14
Future work
15