Fuel Performance Modeling of High Burnup Mixed Oxide Fuel for Hard Spectrum LWRs
by Yanin Sukjai
M.S., Nuclear Science and Engineering (2014) Massachusetts Institute of Technology
B.Eng., Mechanical Engineering (2001) King Mongkut’s University of Technology North Bangkok
Submitted to the Department of Nuclear Science and Engineering in Partial Fulfillment of the Requirements for the Degree of
DOCTOR OF PHILOSOPHY IN NUCLEAR SCIENCE AND ENGINEERING at the MASSACHUSETTS INSTITUTE OF TECHNOLOGY
February 2018
© 2018 Massachusetts Institute of Technology. All rights reserved
Signature of Author: ______Yanin Sukjai Department of Nuclear Science and Engineering September 28, 2017
Certified by: ______Koroush Shirvan, Ph.D. Assistant Professor of Nuclear Science and Engineering Thesis Supervisor
Certified by: ______Ronald G. Ballinger, Ph.D. Professor of Nuclear Science and Engineering Professor of Materials Science and Engineering Thesis Reader
Accepted by: ______Ju Li, Ph.D. Battelle Energy Alliance Professor of Nuclear Science and Engineering Professor of Materials Science and Engineering Chair, Department Committee on Graduate Students
1
This Page Left Intentionally Blank
2
Fuel Performance Modeling of High Burnup Mixed Oxide Fuel for Hard Spectrum LWRs
by
Yanin Sukjai
Submitted to the Department of Nuclear Science and Engineering on September 28, 2017 in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy in Nuclear Science and Engineering
ABSTRACT
According to the future of the nuclear fuel cycle study at MIT, a reactor with a conversion ratio around one can achieve desired objectives in the long-term sustainability of uranium and reduction of transuranic wastes. This finding relaxes the need for sodium fast reactors (SFR) in a closed-loop nuclear fuel cycle and enables high-conversion light water reactors (HC-LWR) to be used as an alternative. HC-LWRs have two major advantages over SFRs. First, apart from the reactor core, the remaining reactor system can be based on existing LWR technology. Second, extensive operating experience and a proven record of high reliability of LWRs would ease licensing and commercialization processes. Therefore, operating HC-LWRs instead of SFRs may be more economically and technically viable with lower capital and development cost for the near term. This type of reactor is being developed by Hitachi Ltd. under the name of resource- renewable boiling water reactor (RBWR).
This study focuses on RBWR-TB2, transuranic burning version of RBWR. To demonstrate that the RBWR-TB2 can operate safely within design constraints and regulatory limits, the thermo- mechanical behavior of this reactor has been analyzed through fuel performance modeling.
Due to its unique design characteristics, several physical phenomena at high temperature and high burnup typically ignored in most LWR fuel performance codes can potentially become active under RBWR’s operating conditions. These phenomena involve migration of fuel constituents and fission products, the evolution of O/M ratio with burnup, high burnup structure (HBS) formation, accelerated corrosion, hot pressing, gaseous fuel swelling, hydride precipitation and hydrogen migration in the cladding. Semi-empirical models describing porosity and cesium migration behaviors have been replaced with mechanistic models. All of these phenomena have been successfully implemented in a modified version of FRAPCON-3.5 known as FRAPCON-3.5 EP where EP stands for enhanced performance.
The fuel performance comparison between RBWR-TB2 and ABWR fuel rods suggest that because of high axial peaking factors and relatively flat power history, fuel temperature is
3 significantly higher in fissile zones of the RBWR-TB2 leading to various undesirable effects such as excessive fission gas release and cladding deformation. Local fuel burnup in fissile zones of RBWR-TB2 is multiple times higher than that of ABWR leading to excessive fuel swelling, accelerated cladding oxidation, and PCMI at fissile-blanket interfaces. Even if the RBWR-TB2 has to operate under such demanding conditions with a small margin to fuel melting, a steady- state fuel performance analysis still shows that this reactor can operate safely with an acceptable thermo-mechanical performance.
In the future optimization of RBWR-TB2 performance, several fuel design strategies are recommended based on a series of sensitivity studies. The sensitivity study on key design parameters indicates that using annular fuel geometry and more hypostoichiometric fuel (lower O/M ratio) could reduce fuel temperature at high burnup. For better resistance to cladding corrosion and PCMI, it is recommended to increase cladding thickness and decrease fuel density.
Thesis Supervisor: Koroush Shirvan Title: Assistant Professor, Department of Nuclear Science and Engineering
Thesis Reader: Ronald G. Ballinger Title: Professor, Department of Nuclear Science and Engineering and Department of Materials Science and Engineering
4
Acknowledgement
I would like to express my most profound gratitude to my former and late thesis advisor, Professor Mujid S. Kazimi, for his continuous guidance, support, and encouragement during the course of my academic journey. Without his patience and understanding during my time at MIT, this thesis would not be possible. Not only did he serve as my former thesis advisor, but he was also my role model professionally. His ability to analyze complicated technical issues and clearly convey the essential information for the audiences is one of the most important characteristics of a successful researcher. I feel honored to have the opportunity to work with him. May your soul rest in peace; you will always be remembered, Professor Kazimi.
I owe a great debt to my current thesis supervisor, Professor Koroush Shirvan, for his continuous support and help throughout my thesis. This thesis would have been greatly delayed had he not immediately taken me into his supervision and research group. I was deeply impressed and influenced by his diligence, positive attitude, attention to detail, and unending desire for learning. He assisted me in every step of my research: literature review, model simplification, code diagnostics, data analysis, and troubleshooting. More than an advisor, he is a real friend who helped and encouraged me at difficult moments.
I am equally grateful to Professor Ronald Ballinger, my thesis reader, for taking me in as his student immediately after the passing of Professor Kazimi. The discussions we had during group meetings have always been of great value, along with his comments, and suggestions on the thesis. I also would like to thank Professor Ju Li for having accepted my request and served as a thesis committee member. As a world-renowned expert in materials science and engineering, he has broadened my vision and provided insights from a different point of view in this work.
My thanks are also extended to Dr. Aydin Karahan who had previously developed FRAPCON- EP to the point at which this work had started.
I wish to thank NSE colleagues and Thai friends at MIT for their hospitality, friendship, and advice throughout my study. My special thanks go to my fellow classmate and Thai friend, Mr. Ittinop Dumnernchanvanit, who has always support me socially and academically since the first
5 day I landed at the Logan airport.
I would like to thank my parents, Mr. Ekkachai Sukjai and Mrs. Nichapa Sukjai, my brother, Mr. Thanakorn Sukjai, my wife, Ms. Lalita Urasuk, my son, Ratchanin Sukjai, and my daughters, Pimonnart Sukjai and Natacha Sukjai, for their unconditional love and support. Thanks for taking care of me and being a source of knowledge, morale, and inspiration throughout my entire life.
Finally, I would like to express my deep appreciation to the Royal Thai Government, the King Mongkut's University of Technology Thonburi (KMUTT), Hitachi Limited of Japan, and the Center for Advanced Nuclear Energy Systems (CANES) at MIT for their financial support for my graduate study.
6
Table of Contents Abstract ...... 3 Acknowledgement ...... 5 Table of Contents ...... 7 List of Figures ...... 10 List of Tables ...... 21 Nomenclature ...... 22 Chapter 1 ...... 23 1.1 Thesis objective ...... 23 1.2 Motivation ...... 23 1.2.1 Transuranic waste incineration in high-conversion LWRs ...... 23 1.3 Scope of work...... 30 1.4 Nuclear fuel performance modeling ...... 32 1.4.1 FRAPCON-3 fuel performance code ...... 34 1.4.2 FRAPCON-EP ...... 37 1.5 Thesis organization ...... 39 Chapter 2 ...... 41 2.1 History of HC-LWR development ...... 41 2.2 General description and design characteristics of RBWR-TB2 ...... 54 2.3 Material challenges associated with RBWR-TB2 ...... 60 Chapter 3 ...... 64 3.1 Parameters affecting thermal conductivity ...... 65 3.2 Thermal conductivity correlations for mixed oxide fuels ...... 68 3.3 Comparison of thermal conductivity correlations ...... 86 3.4 Benchmarking with experimental data ...... 96
3.5 Effect of PuO2 content on MOX thermal conductivity ...... 106 Chapter 4 ...... 117 4.1 Porosity migration and central void formation ...... 118 4.2 Plutonium migration ...... 125 4.3 Oxygen-to-metal ratio variation with burnup ...... 133 4.4 Oxygen-to-metal ratio variation with temperature ...... 144
7
4.5 Cesium migration and formation of Joint Oxyde-Gaine (JOG) ...... 151 4.6 High-burnup structure and RIM porosity ...... 171 4.7 Hot-pressing ...... 175 4.8 Fuel swelling from gaseous fission products ...... 177 4.8.1 Intragranular gas swelling ...... 179 4.8.2 Intergranular gas swelling ...... 181 4.9 Accelerated corrosion at high burnup ...... 187 4.10 Hydrogen migration and hydride precipitation in cladding ...... 195 Chapter 5 ...... 203 5.1 Modification of FRAPCON-3.5 EP for fast reactor conditions ...... 203 5.2 Validation of plutonium and porosity migration models ...... 205 5.3 Validation of cesium migration model in-pile reactor experiments ...... 216 5.4 Validation of cesium migration model with out-of-pile experiments ...... 236 5.5 Sensitivity study of cesium migration model on cesium fuel swelling rate ...... 248 5.6 Validation of hydrogen migration model ...... 262 Chapter 6 ...... 269 6.1 Reactor condition and fuel rod geometry ...... 269 6.2 Radial power profile and fast neutron flux ...... 270 6.3 Power history and axial peak factor ...... 276 6.4 Results of fuel performance simulation ...... 279 6.4.1 Rod-average and local fuel burnup ...... 279 6.4.2 Average fuel temperature ...... 281 6.4.3 Fuel centerline temperature ...... 282 6.4.4 Plenum pressure and fission gas release ...... 287 6.4.5 Cladding corrosion ...... 290 6.4.6 Structural radial gap and interfacial pressure ...... 292 6.4.7 Cladding hoop stress and strain ...... 294 6.4.8 Porosity migration and central void formation ...... 296 6.4.9 Plutonium redistribution ...... 298 6.4.10 Oxygen-to-metal ratio radial redistribution ...... 300 6.4.11 Cesium migration and JOG formation ...... 302
8
6.4.12 Hydrogen redistribution and precipitation ...... 306 6.5 Parametric study on important fuel design parameters ...... 309 6.5.1 Initial fuel-clad gap thickness ...... 310 6.5.2 Fuel density ...... 317 6.5.3 Oxygen-to-metal ratio ...... 326 6.5.4 Helium pressure ...... 333 6.5.5 Central void diameter ...... 340 6.5.6 Cladding thickness ...... 351 6.6 Code-to-code comparison between FRAPCON-3.5 and FRAPCON-3.5 EP ...... 361 6.6.1 RBWR-TB2 design ...... 362 6.6.2 ABWR design ...... 376 Chapter 7 ...... 390 7.1 Summary ...... 390 7.1.1 Evaluation of thermal conductivity correlation for mixed oxide fuel ...... 390 7.1.2 Development of physical phenomena at high temperature and high burnup ...... 391 7.1.3 Validation of FRAPCON-3.5 EP with experimental data ...... 394 7.1.4 Fuel performance modeling of RBWR-TB2 ...... 396 7.2 Conclusions ...... 399 7.3 Recommendations for Future Work ...... 401 7.3.1 Full-core fuel performance analysis for RBWR-TB2 ...... 401 7.3.2 Fuel performance modeling of RBWR-TB2 in transient conditions ...... 403 7.3.3 Out-of-pile experiments for fuel constituent migration under temperature gradient… ...... 403 7.3.4 Experimental evaluation of MOX thermal conductivity at high burnup ...... 404 7.3.5 In-pile experiment for accelerated corrosion at high neutron fluence ...... 405 7.3.6 Theoretical and experimental study of O/M ratio evolution with burnup ...... 406 References ...... 408 Appendix A: Material properties of HT-9 and SS-304 Stainless Steel ...... 432 Appendix B: Sample input files ...... 458
9
List of Figures Figure 1: Once-through fuel cycle [3]...... 25 Figure 2: Closed-loop fuel cycle [3]...... 26 Figure 3: Natural uranium requirement over the course of simulation [5]...... 28 Figure 4: Total amount of TRU in the system [5]...... 28 Figure 5: High level waste in geological repository [5]...... 29 Figure 6: Levelized cot of electricity [5]...... 30 Figure 7: Complexity of nuclear fuel performance modeling [6]...... 34 Figure 8: Simplified FRAPCON-3 flowchart [7]...... 36 Figure 9: Classification of high-conversion LWR [11]...... 42 Figure 10: Moderator-to-fuel ratio and breeding ratio comparison [36]...... 50 Figure 11: Rate of formation and consumption of TRUs [37]...... 53 Figure 12: (a) Reactor pressure vessel of RBWR (b) Horizontal cross-section of RBWR reactor core [42]...... 55 Figure 13: Axial and hexagonal configuration of RBWR-TB2 fuel bundle [42]...... 56 Figure 14: Horizontal configuration of RBWR-TB2 [42]...... 56 Figure 15: Axial LHGR of the RBWR-TB2 as a function of core height [39]...... 57 Figure 16: Local fuel burnup of the RBWR-TB2 as a function of axial node and time step...... 58 Figure 17: Axial void fraction of RBWR-TB2 as a function of relative core height [42]...... 59 Figure 18: Comparison of normalized neutron spectra for the RBWR-AC, RBWR-Th, SFR and ABWR [43]...... 60 Figure 19: Axial peaking factor vs. core height of ABWR and RBWR-TB2...... 61 Figure 20: Effect of fuel swelling and thermal stress [44]...... 62 Figure 21: Consequence of pellet-cladding interaction (PCI) [44]...... 62 Figure 22: Parameters affecting thermal conductivity [56]...... 66
Figure 23: Electron density of state of UO2 (a), ThO2 (b), and PuO2 (c). Cross-hatched bands indicate occupied level [59]...... 68 Figure 24: Thermal conductivity of MOX at 0 MWd/kgHM and O/M = 2.0...... 90 Figure 25: Thermal conductivity of MOX at 100 MWd/kgHM and O/M = 2.0...... 90 Figure 26: Thermal conductivity of MOX at 0 MWd/kgHM and O/M = 1.95...... 91 Figure 27: Thermal conductivity of MOX at 100 MWd/kgHM and O/M = 1.95...... 92 Figure 28: Thermal conductivity of MOX at 0 MWd/kgHM and O/M = 2.0...... 93 Figure 29: Thermal conductivity of MOX at 100 MWd/kgHM and O/M = 2.0...... 94 Figure 30: Thermal conductivity of MOX at 0 MWd/kgHM and O/M = 1.95...... 95 Figure 31: Thermal conductivity of MOX at 100 MWd/kgHM and O/M = 1.95...... 96
Figure 32: Measured vs. calculated thermal conductivity of MOX and UO2 using Duriez- modified NFI and Duriez-Lucuta correlations...... 99
Figure 33: Measured vs. calculated thermal conductivity of MOX and UO2 using Inoue-modified NFI and Inoue-Lucuta correlations...... 100
10
Figure 34: Measured vs. calculated thermal conductivity of MOX and UO2 using Baron, Halden, and Amaya correlations...... 101
Figure 35: Measured vs. calculated thermal conductivity of MOX and UO2 using Duriez- modified NFI and Duriez-Lucuta correlations...... 103
Figure 36: Measured vs. calculated thermal conductivity of MOX and UO2 using Inoue-modified NFI and Inoue-Lucuta correlations...... 104
Figure 37: Measured vs. calculated thermal conductivity of MOX and UO2 using Baron, Halden, and Amaya correlations...... 105
Figure 38: Thermal conductivity of fresh MOX as a function of PuO2 content [78]...... 107 Figure 39: Thermal conductivity correlation recommended by Nichenko and comparison with experimental result (100% TD) [65]...... 110
Figure 40: Effect of PuO2 on thermal conductivity by Matsumoto et al. [69]...... 110 Figure 41: The variation of thermal conductivity of MOX as a function of uranium composition [70]...... 111
Figure 42: Thermal conductivity as a function of temperature at varying PuO2 weight fraction [78]...... 113 Figure 43: Comparison of Duriez-modified NFI correlation with measured thermal conductivity
of MOX at 5 and 30 wt% PuO2 from Gibby [78]...... 113 Figure 44: Calculated thermal conductivity of MOX as a function of PuO2 weight fraction from empirical correlation proposed by Gibby [78]...... 114
Figure 45: Multiplying factor to the Duriez-modified NFI correlation with reference PuO2 weight fraction at 30 wt%...... 115 Figure 46: Cross section of mixed oxide fuel rod irradiated at 56 kW/m to 25 MWd/kgHM [126]...... 119 Figure 47: Schematic of restructured regions [126]...... 120 Figure 48: Comparison of fuel temperature calculated from three-region model and measured porosity [127]...... 122 Figure 49: Electron probe microanalysis (EPMA) false color X-ray maps showing radial distribution of U, Pu, and Am. The concentration increases in the following order: green, yellow and red [140]...... 127 Figure 50: Fuel restructuring and radial distribution of Pu and Am after 10 minutes and 24 hours irradiation at high LHGR conditions [140]...... 128
Figure 51: The oxygen potential at four different radial positions at high burnup UO2 fuel [141]...... 135 Figure 52: Burnup dependence of the oxygen potential at 750 oC for different irradiated oxide fuel [144]...... 136 Figure 53: Oxygen potential measurement of irradiated Phoenix fuel of initial composition of
(U0.8Pu0.2)U1.98 [146]...... 136 Figure 54: A unit volume of fresh and irradiated mixed oxide fuel as a constant mass system [126]...... 137
11
Figure 55: partial molar enthalpy and entropy of oxygen in mixed oxide fuel [126]...... 142 Figure 56: Axial 137Cs activity profiles in fuel rod C-10 and C-19 [159]...... 152 Figure 57: The axial distribution of Cs and I under simulated temperature gradients [173]...... 153 Figure 58: 137Cs and 134Cs radial distribution profile at 22.9 and 48.26 MWd/kgHM [170]. .... 154 Figure 59: 137Cs and 134Cs radial distribution profile at higher fuel temperature [170]...... 155 Figure 60: Localized cladding strain as a result of cesium migration [161]...... 155 Figure 61: X-ray mapping of the JOG compounds at peak power node at high burnup. The gap region is filled with chemical compounds of mainly Cs, Mo and O [163]...... 156 Figure 62: Gamma scan of 137Cs of solid and annular pellets [164]...... 157 Figure 63: EPMA mapping images of JOG showing complex compounds of Cs and Mo [172]...... 159 Figure 64: Axial distribution of 137Cs intensity along the fission column [164]...... 160 Figures 65: Ceramographs of transverse section from bottom to top of fuel column [164]...... 160 Figure 66: EPMA mapping of the gap region of position B [164]...... 162 Figure 67: EPMA mapping of the gap region of position C [164]...... 162 Figure 68: As-fabricated grain (left) vs. high burnup structure (right) [49]...... 173 Figure 69: Porosity evolution as a function of local burnup in HBS region...... 175
Figure 70: (a) Ceramograph of gas bubbles in irradiated UO2 [196] and (b) schematic representation of intra-granular and inter-granular gas bubbles [197]...... 179 Figure 71: A schematic representation of a TKD grain of equal side [199]...... 182 Figure 72: Schematic illustration of connected tunnel of grain edge bubbles [200]...... 183 Figure 73: Two-stage oxidation process of zirconium alloys [201]...... 188
Figure 74: Possible phases of zirconium oxide (ZrO2) in LWR conditions [202]...... 188 Figure 75: Schematic of zirconium oxide formation showing barrier and porous layers [201]. 189 Figure 76: Schematic of 3-stage oxidation model proposed by Zhou et al. [203]...... 191 Figure 77: Oxide thickness, hydrogen content and rod growth as a function of fast fluence [205]...... 192 Figure 78: Acceleration in oxide layer thickness at high burnup in Zircaloy-2 [208]...... 193 Figure 79: Acceleration in HPUF at high burnup in Zircaloy-2 [203]...... 193 Figure 80: Hydrogen concentration as predicted by a default correlation in FRAPCON-3.5. ... 195 Figure 81: Non-uniform distribution of zirconium hydride in a typical LWR cladding [212]. .. 196 Figure 82: Radial distribution of hydride in Zircaloy cladding [213] [214]...... 196 Figure 83: Azimuthal distribution of hydride in Zircaloy cladding [213]...... 197 Figure 84: Hydride distribution at the inter-pellet gap cladding compared to mid-pellet cladding [215]...... 198 Figure 85: Power history of B11 experiment [227]...... 207 Figure 86: Axial peaking factor of B11 experiment...... 207 Figure 87: Central void diameter as a function of axial height for Am-1-1-1 rod...... 208 Figure 88: Central void diameter as a function of axial height for Am-1-2-1 rod...... 209 Figure 89: Central void diameter as a function of axial height for Am-1-2-2 rod...... 209
12
Figure 90: Plutonium concentration as a function of fuel radius for Am-1-1-1 rod...... 210 Figure 91: Plutonium concentration as a function of fuel radius for Am-1-2-1 rod...... 210 Figure 92: Plutonium concentration as a function of fuel radius for Am-1-2-2 rod...... 211 Figure 93: Power history of B14 experiment [227]...... 212 Figure 94: Axial peaking factor of B14 experiment [124]...... 213 Figure 95: Central void diameter as a function of axial height for PTM0001 rod...... 214 Figure 96: Central void diameter as a function of axial height for PTM0002 rod...... 214 Figure 97: Central void diameter as a function of axial height for PTM0003 rod...... 215 Figure 98: Central void diameter as a function of axial height for PTM0010 rod...... 215 Figure 99: Plutonium concentration as a function of fuel radius for PTM0001 rod...... 216 Figure 100: Plutonium concentration as a function of fuel radius for PTM0010 rod...... 216 Figure 101: Power history of ACO-1 fuel rod...... 219 Figure 102: Axial peaking factor of ACO-1 fuel rod...... 219 Figure 103: Relative 137Cs concentration as a function of axial height of ACO-1 fuel rod...... 220 Figure 104: Cladding strain at EOL as a function of axial height of ACO-1 fuel rod...... 221 Figure 105: Power history of fuel rods in ACO-3 experiment [166]...... 222 Figure 106: Peak cladding midwall temperature of fuel rods in ACO-3 experiment [166]...... 223 Figure 107: Axial peaking factor of ACO-3 experiment [166]...... 223 Figure 108: Relative 137Cs concentration as a function of axial height of 150073 fuel rod...... 225 Figure 109: Relative 137Cs concentration as a function of axial height of 150080 fuel rod...... 225 Figure 110: Relative 137Cs concentration as a function of axial height of 150088 fuel rod...... 225 Figure 111: Relative 137Cs concentration as a function of axial height of 150094 fuel rod...... 226 Figure 112: Cladding strain at EOL as a function of axial height of 150073 fuel rod...... 227 Figure 113: Cladding strain at EOL as a function of axial height of 150080 fuel rod...... 227 Figure 114: Cladding strain at EOL as a function of axial height of 150088 fuel rod...... 228 Figure 115: Cladding strain at EOL as a function of axial height of 150094 fuel rod...... 228 Figure 116: Histories of maximum LHGR and maximum cladding temperature of fuel rods in C3M experiment [165]...... 230 Figure 117: Axial profiles of fast neutron fluence and life-averaged cladding midwall temperature of C3M experiment [165]...... 231 Figure 118: Relative 137Cs concentration as a function of axial height of G305 fuel rod...... 232 Figure 119: Relative 137Cs concentration as a function of axial height of G339 fuel rod...... 232 Figure 120: Relative 137Cs concentration as a function of axial height of G357 fuel rod...... 233 Figure 121: Cladding strain at EOL as a function of axial height of G305 fuel rod...... 234 Figure 122: Cladding strain at EOL as a function of axial height of G339 fuel rod...... 235 Figure 123: Cladding strain at EOL as a function of axial height of G357 fuel rod...... 235 Figure 124: Schematic representation of experimental device for Experiment #1 [170]...... 237 Figure 125: Schematic representation of experimental device for Experiment #2 [170]...... 238 Figure 126: Schematic representation of stainless steel test capsule showing the location where the additives (Cs, I, Te, and Mo) were introduced at hot end...... 239
13
o Figure 127: Axial cesium distribution of Peehs’s Experiment #1, Tmax = 1000 C...... 242 o Figure 128: Axial cesium distribution of Peehs’s Experiment #1, Tmax = 1200 C...... 242 o Figure 129: Axial cesium distribution of Peehs’s Experiment #1, Tmax = 1320 C...... 243 o Figure 130: Axial cesium distribution of Peehs’s Experiment #1, Tmax = 1400 C...... 243 Figure 131: Location of Cs-137 peaks as function of time of Peehs’s Experiment #2...... 245 Figure 132: Axial distribution of Cs-137 activity before and after thermal treatment ...... 246 Figure 133: Comparison of cesium diffusion coefficients used during simulation with correlations from literature ...... 248 Figure 134: Relative 137Cs concentration as a function of axial height of 150073 fuel rod...... 251 Figure 135: Relative 137Cs concentration as a function of axial height of 150080 fuel rod...... 251 Figure 136: Relative 137Cs concentration as a function of axial height of 150088 fuel rod...... 252 Figure 137: Relative 137Cs concentration as a function of axial height of 150094 fuel rod...... 252 Figure 138: Average fuel temperature of high and low swelling models as a function of time. 254 Figure 139: Average fuel temperature at EOL of high and low swelling models as a function of axial node...... 254 Figure 140: Centerline fuel temperature of high and low swelling models as a function of time...... 255 Figure 141: Centerline fuel temperature at EOL of high and low swelling models as a function of axial node...... 255 Figure 142: Plenum pressure of high and low swelling models as a function of time...... 256 Figure 143: Fission gas release of high and low swelling models as a function of time...... 257 Figure 144: Minimum structural radial gap of high and low swelling models as a function of time...... 258 Figure 145: Structural radial gap at EOL of high and low swelling models as a function of axial node...... 258 Figure 146: Interfacial pressure of high and low swelling models as a function of time...... 259 Figure 147: Cladding hoop stress of high and low swelling models as a function of time...... 260 Figure 148: Cladding hoop stress at EOL of high and low swelling as a function of axial node...... 260 Figure 149: Cladding hoop strain of high and low swelling models as a function of time...... 261 Figure 150: Cladding hoop strain at EOL of high and low swelling models as a function of axial node...... 262 Figure 151: Axial distribution of hydrogen concentration in Sawatzky’s experiment #1...... 264 Figure 152: Axial distribution of hydrogen concentration in Sawatzky’s experiment #2...... 265 Figure 153: Power history of fuel rod 1079 of Gravelines nuclear power station [242]...... 267 Figure 154: Axial peaking factor used in Lacroix [244] and this study...... 267 Figure 155: Total hydrogen concentration as a function of clad radius...... 268 Figure 156: Normalized radial peaking factors for RBWR-TB2 and ABWR...... 272 Figure 157: Cross-section view of RBWR-TB2 assemblies in SERPENT...... 273 Figure 158: Fast neutron flux of RBWR-TB2 at 0, 40, and 80 MWd/kgHM...... 274
14
Figure 159: Specific power of RBWR-TB2 at 0, 40, and 80 MWd/kgHM...... 275 Figure 160: Specific fast flux of RBWR-TB2 at 0, 40, and 80 MWd/kgHM...... 276 Figure 161: Power history of RBWR-TB2 and ABWR...... 278 Figure 162: Axial peaking factor of RBWR-TB2 and ABWR...... 278 Figure 163: Rod-average fuel burnup of RBWR-TB2 and ABWR as a function of time...... 280 Figure 164: Local fuel burnup of RBWR-TB2 and ABWR at EOL...... 280 Figure 165: Average fuel temperature of RBWR-TB2 and ABWR as a function of time...... 281 Figure 166: Average fuel temperature of RBWR-TB2 and ABWR as a function of relative axial length...... 282 Figure 167: Centerline fuel temperature of RBWR-TB2 and ABWR as a function of time...... 283 Figure 168: Centerline fuel temperature of RBWR-TB2 and ABWR as a function of relative axial length...... 284 Figure 169: Axial variation of fuel centerline and fuel melting temperature at EOL of RBWR- TB2...... 285 Figure 170: Axial variation of fuel melting temperature and plutonium content at EOL of RBWR-TB2...... 286 Figure 171: Reduction of fuel melting temperature of RBWR-TB2 with plutonium content at BOL and EOL...... 287 Figure 172: Plenum pressure of RBWR-TB2 and ABWR as a function of time...... 288 Figure 173: Fission gas release of RBWR-TB2 and ABWR as a function of time...... 290 Figure 174: Oxide layer thickness of RBWR-TB2 and ABWR as a function of time...... 291 Figure 175: Average hydrogen concentration of RBWR-TB2 and ABWR as a function of time...... 292 Figure 176: Structural radial gap of RBWR-TB2 and ABWR as a function of time...... 293 Figure 177: Interfacial pressure of RBWR-TB2 and ABWR as a function of time...... 294 Figure 178: Cladding hoop stress of RBWR-TB2 and ABWR at peak axial node as a function of time...... 295 Figure 179: Cladding hoop strain of RBWR-TB2 and ABWR at peak axial node as a function of time...... 296 Figure 180: Central void diameter of RBWR-TB2 along axial nodes...... 297 Figure 181: Central void diameter of ABWR along axial nodes...... 298 Figure 182: Radial distribution of plutonium of RBWR-TB2 at peak axial node...... 299 Figure 183: Radial distribution of plutonium of ABWR at peak axial node...... 300 Figure 184: Radial distribution of O/M ratio of RBWR-TB2 at peak axial node...... 301 Figure 185: Radial distribution of O/M ratio of ABWR at peak axial node...... 302 Figure 186: Axial distribution of cesium of RBWR-TB2...... 303 Figure 187: Axial distribution of cesium of ABWR...... 304 Figure 188: Radial distribution of cesium of RBWR-TB2...... 305 Figure 189: Radial distribution of cesium of ABWR...... 305
15
Figure 190: Axial distribution of hydride precipitation of RBWR-TB2 at cladding outer surface...... 307 Figure 191: Axial distribution of hydride precipitation of ABWR at cladding outer surface. ... 307 Figure 192: Radial distribution of hydrogen solute of RBWR-TB2...... 308 Figure 193: Radial distribution of hydrogen solute in ABWR...... 309 Figure 194: Average fuel temperature of RBWR-TB2 at 55 and 110 μm as a function of time. 311 Figure 195: Centerline fuel temperature of RBWR-TB2 at 55 and 110 μm as a function of time...... 311 Figure 196: Average fuel temperature of RBWR-TB2 at 55 and 100 μm at EOL as a function of axial node...... 312 Figure 197: Centerline fuel temperature of RBWR-TB2 at 55 and 100 μm at EOL as a function of axial node...... 313 Figure 198: Fission gas release of RBWR-TB2 at 55 and 110 μm as a function of time...... 314 Figure 199: Plenum pressure of RBWR-TB2 at 55 and 110 μm as a function of time...... 314 Figure 200: Interfacial pressure of RBWR-TB2 at 55 and 110 μm at peak axial node as a function of time...... 315 Figure 201: Cladding hoop stress of RBWR-TB2 at 55 and 110 μm at peak axial node as a function of time...... 316 Figure 202: Cladding hoop strain of RBWR-TB2 at 55 and 110 μm at peak axial node as a function of time...... 317 Figure 203: Rod average burnup of RBWR-TB2 at fuel density of 80-95 %TD as a function of time...... 318 Figure 204: Local burnup of RBWR-TB2 at fuel density of 80-95 %TD as a function of axial node...... 319 Figure 205: Average fuel temperature of RBWR-TB2 at fuel density of 80- 95 %TD as a function of time...... 320 Figure 206: Centerline fuel temperature of RBWR-TB2 at fuel density of 80-95 %TD as a function of time...... 321 Figure 207: FGR of RBWR-TB2 at fuel density of 80-95 %TD as a function of time...... 322 Figure 208: Plenum pressure of RBWR-TB2 at fuel density of 80-95 %TD as a function of time...... 322 Figure 209: Total void volume of RBWR-TB2 at fuel density of 80-95 %TD as a function of time...... 323 Figure 210: Interfacial pressure of RBWR-TB2 at fuel density of 80-95 %TD as a function of time...... 324 Figure 211: Cladding hoop stress of RBWR-TB2 at fuel density of 80-95 %TD as a function of time...... 325 Figure 212: Cladding hoop strain of RBWR-TB2 at fuel density of 80-95 %TD as a function of time...... 326
16
Figure 213: Average fuel temperature of RBWR-TB2 at O/M ratios of 1.90-2.00 as a function of time...... 327 Figure 214: Centerline fuel temperature of RBWR-TB2 at O/M ratios of 1.90-2.00 as a function of time...... 328 Figure 215: Fission gas release of RBWR-TB2 at O/M ratios of 1.90-2.00 as a function of time...... 329 Figure 216: Plenum pressure of RBWR-TB2 at O/M ratios of 1.90-2.00 as a function of time. 330 Figure 217: Interfacial pressure of RBWR-TB2 at O/M ratios of 1.90-2.00 as a function of time...... 331 Figure 218: Cladding hoop stress of RBWR-TB2 at O/M ratios of 1.90-2.00 as a function of time...... 332 Figure 219: Cladding hoop strain of RBWR-TB2 at O/M ratios of 1.90-2.00 as a function of time...... 332 Figure 220: Average fuel temperature of RBWR-TB2 at helium pressure of 1.0-4.0 MPa as a function of time...... 334 Figure 221: Centerline fuel temperature of RBWR-TB2 at helium pressure of 1.0-4.0 MPa as a function of time...... 335 Figure 222: Fission gas release of RBWR-TB2 at helium pressure of 1.0-4.0 MPa as a function of time...... 336 Figure 223: Plenum pressure of RBWR-TB2 at helium pressure of 1.0-4.0 MPa as a function of time...... 337 Figure 224: Interfacial pressure of RBWR-TB2 at helium pressure of 1.0-4.0 MPa as a function of time...... 338 Figure 225: Cladding hoop stress of RBWR-TB2 at helium pressure of 1.90-2.00 as a function of time...... 339 Figure 226: Cladding hoop strain of RBWR-TB2 at helium pressure of 1.90-2.00 as a function of time...... 340 Figure 227: Rod average burnup of RBWR-TB2 at central void diameter of 0.0-2.0 mm as a function of time...... 341 Figure 228: Local burnup of RBWR-TB2 at central void diameter of 0.0-2.0 mm as a function of axial node...... 342 Figure 229: Average fuel temperature of RBWR-TB2 at central void diameter of 0.0-2.0 mm as a function of time...... 343 Figure 230: Average fuel temperature at EOL of RBWR-TB2 at central void diameter of 0.0-2.0 mm as a function of axial node...... 344 Figure 231: Centerline fuel temperature of RBWR-TB2 at central void diameter of 0.0-2.0 mm as a function of time...... 345 Figure 232: Fission gas release of RBWR-TB2 at central void diameter of 0.0-2.0 mm as a function of time...... 346
17
Figure 233: Plenum pressure of RBWR-TB2 at central void diameter of 0.0-2.0 mm as a function of time...... 347 Figure 234: Void volume of RBWR-TB2 at central void diameter of 0.0-2.0 mm as a function of time...... 348 Figure 235: Interfacial pressure of RBWR-TB2 at central void diameter of 0.0-2.0 mm as a function of time...... 349 Figure 236: Cladding hoop stress of RBWR-TB2 at central void diameter of 0.0-2.0 mm as a function of time...... 350 Figure 237: Cladding hoops strain of RBWR-TB2 at central void diameter of 0.0-2.0 mm as a function of time...... 351 Figure 238: Average fuel temperature of RBWR-TB2 at cladding thickness of 0.6-0.9 mm as a function of time...... 353 Figure 239: Centerline fuel temperature of RBWR-TB2 at cladding thickness of 0.6-0.9 mm as a function of time...... 354 Figure 240: Fission gas release of RBWR-TB2 at cladding thickness of 0.6-0.9 mm as a function of time...... 355 Figure 241: Plenum pressure of RBWR-TB2 at cladding thickness of 0.6-0.9 mm as a function of time...... 356 Figure 242: Oxide layer thickness of RBWR-TB2 at cladding thickness of 0.6-0.9 mm as a function of time...... 357 Figure 243: Clad-average hydrogen concentration of RBWR-TB2 at cladding thickness of 0.6- 0.9 mm as a function of time...... 358 Figure 244: Interfacial pressure of RBWR-TB2 at cladding thickness of 0.6-0.9 mm as a function of time...... 359 Figure 245: Cladding hoop stress of RBWR-TB2 at cladding thickness of 0.6-0.9 mm as a function of time...... 360 Figure 246: Cladding hoop strain of RBWR-TB2 at cladding thickness of 0.6-0.9 mm as a function of time...... 361 Figure 247: Average fuel temperature of RBWR-TB2 as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 364 Figure 248: Average fuel temperature of RBWR-TB2 at EOL as calculated by FRAPCON-3.5 and 3.5 EP as a function of axial node...... 365 Figure 249: Centerline fuel temperature of RBWR-TB2 as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 366 Figure 250: Centerline fuel temperature of RBWR-TB2 at EOL as calculated by FRAPCON-3.5 and 3.5 EP as a function of axial node...... 367 Figure 251: Fission gas release of RBWR-TB2 as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 368 Figure 252: Plenum pressure of RBWR-TB2 as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 369
18
Figure 253: Oxide layer thickness of RBWR-TB2 as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 370 Figure 254: Hydrogen concentration of RBWR-TB2 as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 371 Figure 255: Structural radial gap of RBWR-TB2 as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 372 Figure 256: Interfacial pressure of RBWR-TB2 as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 374 Figure 257: Cladding hoop stress of RBWR-TB2 as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 375 Figure 258: Cladding hoop strain of RBWR-TB2 as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 376 Figure 259: Power history of an average fuel rod in ABWR...... 377 Figure 260: Average fuel temperature of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 378 Figure 261: Centerline fuel temperature of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 380 Figure 262: Fission gas release of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 381 Figure 263: Plenum pressure of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 382 Figure 264: Total void volume of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 382 Figure 265: Fuel axial extension of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 383 Figure 266: Oxide layer thickness of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 384 Figure 267: Hydrogen concentration of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 384 Figure 268: Structural radial gap of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 385 Figure 269: Gap conductance of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 386 Figure 270: Interfacial pressure of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 387 Figure 271: Cladding hoop stress of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 388 Figure 272: Cladding hoop strain of ABWR as calculated by FRAPCON-3.5 and 3.5 EP as a function of time...... 389
19
Figure 273: Sample core map of RedTail predictions for maximum plenum pressure for (a) Zirconium-based cladding and (b) SiC-based cladding [54]...... 402 Figure A.1: Poisson’s ratio as a function of temperature Zircaloy-2, SS-304 and HT-9...... 433 Figure A.2: Specific heat as a function of temperature Zircaloy-2, SS-304 and HT-9...... 435 Figure A.3: Thermal conductivity as a function of temperature Zircaloy-2, SS-304 and HT-9. 436 Figure A.4: Thermal expansion as a function of temperature Zircaloy-2, SS-304 and HT-9. ... 438 Figure A.5: Void swelling as a function of neutron fluence for SS-304...... 440 Figure A.6: Void swelling as a function of neutron fluence for HT-9...... 440 Figure A.7: Irradiation growth as a function of neutron fluence for Zircaloy-2, SS-304 and HT-9...... 441 Figure A.8: Cladding hardness as a function of temperature for Zircaloy-2, SS-304 and HT-9. 442 Figure A.9: Emissivity as a function of temperature for Zircaloy-2, SS-304 and HT-9...... 443 Figure A.10: Elastic modulus as a function of temperature for Zircaloy-2, SS-304 and HT-9. . 446 Figure A.11: Shear modulus as a function of temperature for Zircaloy-2, SS-304 and HT-9. ... 446 Figure A.12: Oxide growth rate constant as a function of temperature for SS-304...... 451 Figure A.13: Oxide thickness as a function of time for HT-9...... 451 Figure A.14: Oxide thermal conductivity as a function of temperature for Zircaloy-2, SS-304 and HT-9...... 452 Figure A.15: Yield and ultimate strength as a function of temperature for SS-304...... 454 Figure A.16: Yield and ultimate strength as a function of temperature for HT-9...... 455 Figure A.17: Multiplying factor to account for irradiation hardening effect of stainless steel. . 455
20
List of Tables Table 1: Spent fuel inventories in spent fuel pools and dry-cask storage as of the end of 2007 [1]...... 24 Table 2: Plant specification and core design [42]...... 54 Table 3: Experimental data of MOX thermal conductivity...... 97
Table 4: Experimental data of UO2 thermal conductivity at high burnup...... 98 Table 5: Root mean square and standard deviation of kcalculated-kmeasured of thermal conductivity correlations...... 102
Table 6: Root mean square and standard deviation of kcalculated-kmeasured of thermal conductivity correlations after exclusion of unirradiated data...... 105 Table 7: Probable chemical and physical states of fission products in near-stoichiometric mixed oxide fuel [126]...... 134 Table 8: Elemental fission product yield in a fast neutron spectrum [126]...... 139 Table 9: Fuel rod characteristics of JOYO B11 experiments [227]...... 206 Table 10: Fuel rod characteristics of JOYO B14 experiments [227]...... 211 Table 11: Design parameters for ACO-1 and ACO-3 experiments [228]...... 217 Table 12: Irradiation condition of fuel rods in ACO-3 experiments [166]...... 221 Table 13: Design parameter of C3M fuel rods [165]...... 229 Table 14: Diffusion coefficient and heat of transport of cesium used in simulation ...... 246 Table 15: Description of Sawatzky’s experiments [241]...... 263 Table 16: Fuel rod design and geometry of fuel rod 1079 of Gravelines nuclear power station [244]...... 266 Table 17: Fuel rod design and geometry of RBWR-TB2 and ABWR...... 270 Table 18: Key design parameters of RBWR-TB2 fuel rods...... 309 Table 19: Comparison of fuel geometry at different cladding thickness...... 352 Table 20: A summary of key design parameters, their effects on fuel performance, and recommended values for further revisions...... 399 Table 21: A comparison of fuel design characteristics of current and future RBWR-TB2...... 401 Table A.1: Density, melting point and heat of fusion of Zircaloy, SS-304, and HT-9...... 432
21
Nomenclature
ABWR advanced boiling water reactor MPa megapascal(s)
BOL beginning-of-life MWd megawatt-day(s)
EOL end-of-life MWth megawatt(s)-thermal
EPRI Electric Power Research Institute NRC Nuclear Regulatory Commission
FBR fast breeder reactor O/M oxygen-to-metal ratio
FGR fission gas release PCMI pellet cladding mechanical interaction FFTF fast flux testing facility PNNL Pacific Northwest National HBS high burnup structure Laboratory HC-LWR high-conversion light water PuO reactor 2 plutonium dioxide
HPUF hydrogen pickup fraction PWR pressurized water reactor
JAEA Japan Atomic Energy Agency RBWR resource-renewable boiling water reactor JOG Joint Oxyde-Gaine RIA reactivity insertion accident kgU kilogram(s) of uranium RPV reactor pressure vessel kgHM kilogram(s) of heavy metals SBA station blackout accident LHGR linear heat generation rate SFR sodium fast reactor LOCA loss of coolant accident T.D. theoretical density LWR light water reactor TRU transuranic elements MOL middle-of-life U-235 uranium-235 MOX mixed oxide U.S. United States of America MIT Massachusetts Institute of uranium dioxide Technology UO2
22
Chapter 1
Introduction
1.1 Thesis objective
The objective of this thesis is to propose and analyze advanced fuel designs for the applications of plutonium and transuranic (TRU) waste incineration in light water reactors (LWR). With TRU burning and multi-recycling capability, this advanced water-cooled reactor is designed for use in a closed-loop nuclear fuel cycle to alleviate the accumulation of spent nuclear fuel from the current fleet of LWRs. This thesis addresses key design challenges of fuel element associated with TRU burning through fuel performance modeling. Thermal, chemical and mechanical behaviors of fuel rods are investigated to ensure an acceptable performance and reliability throughout reactor operation.
1.2 Motivation
1.2.1 Transuranic waste incineration in high-conversion LWRs
Although nuclear power is an extensive and reliable source of carbon-free energy, it is associated with various challenges that it needs to overcome. One of the most critical challenges for nuclear power sustainability is long-term radioactive waste management from depleted reactor fuels. Currently, global production of spent fuel from light water reactors (LWR) is approximately 10,500 metric tons of heavy metal per year, with roughly 8,500 tons of heavy metal going into long-term storage and about 2,000 tons of heavy metal being reprocessed. By the end of 2009, the global inventories of spent fuel were about 240,000 metric tons, mostly stored on-site [1]. Based on this figure in 2009 and the current global production rate of LWR spent fuel, the projected global inventories of spent fuel will have likely increased to about 313,500 metric tons by the end of 2016. For a typical LWR spent fuel with an average burnup of 50 MWd/kgHM, the composition is made of about 93.4% uranium (~0.8% U-235), 5.2% fission products, 1.2% plutonium and 0.2% minor actinides (neptunium, americium, and curium) [1]. During the first hundred years, the decay heat and radioactivity are dominated by fission products. After that, the
23
source of total radioactivity primarily comes from transuranic elements which will last for several hundred thousands of years until the ingestion toxicity of spent fuel become less than that of natural uranium [1]. With current waste disposal method, the time span of the nuclear waste toxicity can persist far beyond the lifespan of human civilization. The top 10 countries holding spent fuel are shown in Table 1. The United States so far has been the largest holder and producer of spent fuel. By the end of 2010, the total U.S. amount of spent fuel was 64,500 tons including 15,350 tons in dry casks [1]. This number has been increased to nearly 70,000 metric tons according to data published by the U.S. Energy Information Administration (EIA) in 2015 [2].
Table 1: Spent fuel inventories in spent fuel pools and dry-cask storage as of the end of 2007 [1]. Country Spent fuel inventory Spent fuel policy (tons of heavy metal) USA 61,000 Direct disposal Canada 38,400 Direct disposal Japan 19,000 Reprocessing France 13,500 Reprocessing Russia 13,000 Direct disposal and reprocessing South Korea 10,900 Storage, disposal undecided Germany 5,850 Direct disposal United Kingdom 5,850 Reprocessing but future unclear Sweden 5,400 Direct disposal Finland 1,600 Direct disposal
Direct disposal is a final step in a once-through open fuel cycle (OTC) where spent nuclear fuel is discarded as waste in geological repositories without extracting any fissile materials left in the spent fuel. Figure 1 illustrates the concept of the once-through fuel cycle. Although not very effective in term of resource utilization, the once-through fuel cycle offers the most convenient way to handle spent nuclear fuel. To date, this is the most economic fuel cycle to operate nuclear power plants because the fuel cost represents only a small fraction of total cost of electricity generation. Most of the world’s reactors operate based on this cycle [3].
24
Figure 1: Once-through fuel cycle [3].
Nevertheless, the once-through fuel cycle has some disadvantages in terms of fuel sustainability and waste disposal. The current knowledge of uranium resources indicates that there are some 3.5 million tons of uranium exploitable at below $80 per kg [4]. The current estimate of all expected uranium resources is four times more. At the current consumption rate of fresh uranium ore of about 67,000 tons per year, this figure will continue to increase as nuclear power expands [4] . It is expected that by the end of 2035 the consumption will be in the range of 98,000 and 136,000 tons of uranium per year [4]. Although the uranium resource is plentiful at least for the first half of the 21th century, long-term uranium scarcity in the next 200 years is inevitable. Assuming no additional discovery of large uranium resources, the 14 million tons of uranium ore would be depleted by 2226 at the current consumption rate of 67,000 tons per year. However, if the consumption rate continues to increase, the year of uranium depletion should occur earlier: 2159 and 2119 for the consumption rate of 98,000 and 136,000 tons, respectively. Another issue for the once-through fuel cycle is the disposal of spent nuclear fuel. In this case, the capacity of long-term repositories may not be sufficient. In fact, suitable sites for geological repository may not be available in certain countries.
It is understandable that the land area and geological characteristics of each country highly influence the decision on disposal policy of spent nuclear fuel. It can be seen that, for countries with large area like the United States, Russia and Canada or having low population density like Sweden and Finland, direct disposal of spent fuel is more favorable. However, for countries with limited land area and high population density, reprocessing seems to be unavoidable.
Reprocessing refers to a partially closed-looped nuclear fuel cycle where spent nuclear fuel from LWRs is reprocessed to extract the remaining fissile materials for further use. In this cycle, plutonium, transuranic elements (TRUs), and long-live fission products are separated by various
25 chemical processes. The recovered fissile contents are then mixed with depleted uranium, a leftover from uranium enrichment process, passed through fuel fabrication and returned to LWRs for power production. Figure 2 illustrates the closed-loop fuel cycle. Ideally, it is possible to sustain reactor operation without external supply of fissile material because, in fast breeder reactors, they can breed more fuel than they consume. However, a fully closed fuel cycle is not currently achievable because of technology and economic constraints [3].
Figure 2: Closed-loop fuel cycle [3].
Since the economics of nuclear power is dictated by the capital cost of nuclear reactors and their associated infrastructure, the type of nuclear reactor plays a major role in the fuel cycle option [3]. As of now, reprocessing is not a viable option mainly because the capital cost of fast breeder reactors (FBRs) is high and uncertain. Furthermore, to resolve technical and reliability issues of FBR technology, significant investment in research and development is necessary. However, if the capital cost of FBRs could be reduced with a less expensive reactor type and if limited operational experience of FBRs could be replaced with more mature technology then the technical and economic challenges associated with fully closed fuel cycle could be overcome. It will then make the closed fuel cycle more appealing and economically competitive for deployment.
26
When referring to an FBR, a sodium-cooled fast breeder reactor (SFR) is typically chosen as a representative of this technology due to historical reasons, and past experience of building and operating several experimental FBRs. According to a fuel cycle study at MIT, it is suggested that a reactor with conversion ratio (the ratio of fissile material at discharge to fissile material at loading) near unity may be more desirable because it relaxes the technical limitations of using SFR in a closed-loop fuel cycle while maintaining similar objectives in term of long-term sustainability of uranium and reduction in spent fuel mass and repository needs [3].
Technical limitations that hinder the widespread adoption of SFR technology include the fact that metallic sodium as reactor coolant reacts violently with water/steam which is used as working fluid in a power conversion cycle of SFRs. Additionally, sodium coolant can become activated from neutron irradiation during operation. In addition, liquid sodium is opaque, making visual maintenance somewhat difficult and complicated. Several engineering workarounds are required to overcome these issues such as the inclusion of a secondary sodium loop, advanced sealing technology, and high-precision measurement and handling.
Studies also indicate that LWRs can be designed to be self-sustaining with conversion ratio slightly above 1.0; no external fissile material is needed at equilibrium cycle [3] [5]. High- conversion LWRs, LWRs with conversion ratio near unity, have several advantages as an alternative to SFRs. First, apart from the reactor core, the remaining reactor system can be based on existing LWR technology. Second, extensive operating experience and proven record of high reliability of LWRs would ease the licensing and commercialization processes. From that, it can be inferred that the capital and operating costs would be similar to existing LWRs. Consequently, operating high-conversion LWRs instead of SFRs in a close-looped fuel cycle may be more economically and technically viable with higher reliability but lower capital and development cost [3].
To compare the performance of high-conversion LWRs and FBRs in a closed-loop fuel cycle, a fuel cycle simulation has been performed using CAFCA (Code for Advanced Fuel Cycle Analysis) [5] and the results have shown relatively comparable performance in terms of natural uranium requirement, total TRU inventory, repository storage capacity as shown in Figures 3, 4, and 5, respectively. In this study, a reduced-moderation boiling water reactor (RBWR) was used
27 as a representative of high-conversion LWRs with conversion ratio around one while a sodium fast reactor was used as a representative of fast breeder reactor with a similar conversion ratio. It can be clearly seen that fuel recycling and reprocessing in a closed-loop fuel cycle can conserve natural uranium resource and significantly reduce TRU inventory and repository capacity.
6.E+06 OTC Reference Scenario 5.E+06 RBWR Reference Scenario 4.E+06 FBR Reference Scenario 3.E+06
2.E+06
1.E+06
Natural Uranium Requirement [MT] Natural Uranium Requirement 0.E+00 2010 2030 2050 2070 2090 2110 Time [year] Figure 3: Natural uranium requirement over the course of simulation [5].
14000 OTC Reference Scenario 12000 RBWR Reference Scenario 10000
8000 FBR Reference Scenario 6000
4000
TRU in the System [MTHM] TRU 2000
0 2010 2030 2050 2070 2090 2110 Time [year] Figure 4: Total amount of TRU in the system [5].
28
450000 OTC Reference Scenario 400000 RBWR Reference Scenario 350000 FBR Reference Scenario 300000 250000 200000 [MTHM] 150000 100000 50000 HLW in Repository -HLW equivalent YM 0 2010 2030 2050 2070 2090 2110 Time [year] Figure 5: High level waste in geological repository [5].
However, any process that involves fuel recycling and reprocessing will increase the overall cost of electricity. As previously mentioned, the once-throughout fuel cycle currently is the most economically favorable fuel cycle. With a lower capital cost of higher-conversion LWRs, the negative impact on levelized cost of electricity can be mitigated. As shown in Figure 6, the increase in electricity cost of an RBWR was around 5% higher than OTC while that of FBR could be twice as much as RBWR cost increase.
29
94 OTC Reference Scenario 92 RBWR Reference Scenario FBR Reference Scenario 90
88
86 LCoE [$/MWhr] 84
82
80 2010 2020 2030 2040 2050 2060 2070 2080 2090 2100 2110 Time [year] Figure 6: Levelized cot of electricity [5].
Although the RBWR was designed based on the proven technology of the advanced boiling water reactor (ABWR), the unique characteristics of RBWR necessitate further research. From a fuel performance prospective, the technical challenges associated with RBWR are: (1) higher axial peaking factor, (2) higher local burnup, (3) harder neutron spectrum and (4) higher temperature gradient at the seed-blanket interfaces. Additional information on RBWR designs and challenges are described in Chapter 2.
1.3 Scope of work
This thesis focuses on the application of high-conversion LWRs for TRU incineration specifically a transuranic burning version of the reduced-moderation boiling water reactor (RBWR-TB2)1. To demonstrate its performance and safety during operation, thermal, chemical and mechanical behaviors of RBWR-TB2 fuel rods are analyzed through fuel performance modeling. The ultimate goal is to ensure that the fuel rod can survive throughout operation without fuel melting or cladding failure.
1 TB2 stands for transuranic burner 2.
30
RBWR-TB2 fuel is designed to operate at higher linear heat generation rate (LHGR) and higher local burnup beyond typical operating conditions of LWRs; therefore, several physical phenomena generally experienced in fast reactors which operate at high temperature and fuel burnup are expected to occur in RBWR-TB2 fuel rods.
A steady-state fuel performance modeling code, FRAPCON, is used as a primary tool in this study. However, FRAPCON was originally designed to include only physical phenomena relevant to LWRs. Thus, physical models at high temperature and high burnup are not included in the code. In this work, the code required modifications to better represent them in RBWR- TB2 conditions. Furthermore, important physical models such as fuel restructuring, plutonium redistribution, and cesium migration have been validated with experimental data from sodium fast reactors.
For model validation purposes, the code needed to be modified to better reflect experimental conditions in sodium fast reactors. By default, zirconium-based alloys are the only cladding options available in FRAPCON. Since fast reactors use stainless steel as cladding; therefore, the material properties of stainless steels needed to be added. This work includes two types of stainless steels: HT-9 and SS-304 as representative of ferritic/martensitic steels and austenitic steels, respectively. In addition, the coolant properties needed to be changed from water to sodium.
To illustrate the differences in thermo-mechanical behavior of RBWR-TB2 and a typical BWR, fuel performance comparison between RBWR-TB2 and ABWR are performed. Key performance indicators necessary for licensing such as maximum fuel temperature, oxide layer thickness, and cladding stress and strain are compared and discussed. The results from the simulation can also provide critical feedback to the neutronic design team. This work also includes sensitivity analysis on key design parameters so that the results can be used as guidance for fuel design changes.
31
1.4 Nuclear fuel performance modeling
The reactor design process typically involves various calculations: neutronic, thermal hydraulic, fuel performance, and transient/accident analyses. It is an iterative and circular operation as it can start at any process in the cycle but eventually all design constraints have to be met.
Neutronic simulation predicts the behavior of neutrons, fission reactions, and heat generation rate. The scale of calculation can range from a single fuel rod to an entire reactor core. Then, thermal hydraulic simulation describes the behavior of cooling water during operation. It can cover only a single flow channel around a fuel rod or multiple channel of a core. In LWRs, water serves as both coolant and neutron moderator; therefore, the changes in cooling water condition will have direct impact on the neutronic behavior of a reactor core. Perturbations from steady- state conditions such as startup, shutdown, or power fluctuation are covered in transient analysis. Accident analysis addresses the behavior of the reactor and its auxiliary system during accidents such as reactivity insertion accident (RIA), loss of coolant accident (LOCA), or station blackout accident (SBA).
Fuel performance simulation looks into the physical characteristics of the fuel operation; it predicts thermal, mechanical and chemical behavior of a fuel rod under extreme conditions in a reactor core and estimates whether the fuel rod is going to survive throughout operation without fuel melting or cladding failure. It is a multi-disciplinary field as it involves several aspects of sciences: nuclear physics, material science, chemistry, and mechanical engineering. The scale of problem can range from the atomistic scale involving the dynamics of material properties under irradiation up to engineering scales such as the effect of pressure, temperature, and stress on fuel behavior. One of the most fundamental assumptions in fuel performance modeling is that each fuel rod does not directly interact with other rods. Therefore, fuel performance is modeled and analyzed separately for each fuel rod. Typically, fuel performance modeling looks at the peak rod condition i.e. the fuel rod that is exposed to highest burnup or higher LHGR throughout the fuel cycle. If the peak rod can survive; it is assumed that the remaining rods can also survive.
In addition, nuclear regulators usually establish several constraints based on predictions from fuel performance modeling. For example, fuel temperature has to be below melting point.
32
Zircaloy cladding temperature must not exceed 1473 K (1200 oC) to prevent autocatalytic reaction between zirconium and water. Maximum fuel burnup may not exceed 62 MWd/kgHM to prevent cladding failure during transient conditions. End-of-life cladding strain may not exceed 1% to allow coolable geometry of a fuel rod. Therefore, fuel performance modeling is an important part of the reactor licensing process.
The behavior of fuel rods depends on various parameters, and some of them are highly interrelated. Thermal parameters such as the temperature distribution within the fuel rod are not only a function of thermal conductivity and heat transfer coefficient across material boundaries, they also depend on mechanical and chemical effects such as the size of fuel-cladding gap and thickness of cladding oxide layer. Likewise, mechanical parameters such as the stress and strain depend on thermal and chemical effects such as thermal expansion and material property change. Chemical parameters such as fission gas release and oxidation rate varies according to mechanical and thermal states of the materials. All of these inter-correlations in nuclear fuel performance modeling are illustrated in Figure 7. Parameters in red boxes are designated as input while ones in blue boxes are boundary conditions. Cladding and fuel parameters are placed in purple and green boxes, respectively.
33
Figure 7: Complexity of nuclear fuel performance modeling [6].
1.4.1 FRAPCON-3 fuel performance code
FRAPCON is a computer code that calculates the steady-state behavior of light-water reactor fuel rods as a function of burnup. The code was originally developed by the Pacific Northwest National Laboratory (PNNL) for use by the U.S. Nuclear Regulatory Commission (NRC) to calculate thermal, mechanical and material evolution of the fuel and the cladding of a single fuel rod as a function of time and burnup based on initial core conditions and power history up to maximum rod-average burnup of 62 MWd/kgHM [7].
The phenomena modeled in the code include (1) heat conduction through the fuel and cladding to the coolant; (2) cladding elastic and plastic deformations; (3) fuel-cladding mechanical interactions; (4) fission gas release from the fuel and rod internal pressure; and (5) cladding oxidation [7]. The code is designed to simulate the behavior of a single fuel rod under the slowly
34
changing conditions during in-core performance (typically called steady state, although the power derived from the rod can vary during its residence in the core). By definition, steady-state conditions imply that power and boundary conditions changes must be sufficiently slow for a quasi “steady-state” to exist during a portion of the irradiation. . This includes situations such as long periods at constant power and slow power ramps that are typical of normal power reactor operation.
FRAPCON uses a single-channel coolant enthalpy rise model to calculate the axial distribution of the bulk coolant temperature. It uses a finite difference heat conduction model to calculate the temperature distribution within a fuel pellet. Variable mesh spacing is also implemented to accommodate the power peaking at the pellet edge that occurs in high-burnup fuel. The code can calculate the variation with time of all significant fuel rod parameters, including fuel and cladding temperatures, cladding hoop strain, cladding oxidation, fuel irradiation swelling, fuel densification, fission gas release, and rod internal gas pressure. The code calls for material properties from the MATPRO material properties subroutines [7].
FRAPCON can be used to simulate light water and heavy water reactor fuels. Available materials for the fuel pellet, gas gap and cladding include uranium dioxide (UO2) and mixed
oxide pellet ceramic, integrated burnable absorber fuel, zirconium diborate coated UO2, pellet material mixed with gadolinium, zirconium-based alloys cladding which comprises: Zircaloy-2, Zircaloy-4, ZIRLO and M5 [7].
FRAPCON solves the equations iteratively by calculating the interrelated effects of fuel and cladding temperature, rod internal gas pressure, fuel and cladding deformation, release of fission product gases, fuel swelling and densification, cladding thermal expansion and irradiation- induced growth, and cladding corrosion as functions of time and linear power. The calculation procedure of FRAPCON is illustrated in a simplified flowchart as shown in Figure 8. The calculation begins by processing input data. Then, the initial fuel rod state is determined through a self-initialization calculation. Time is advanced according to the input-specified time-step size, a steady-state solution is performed, and the new fuel rod state is determined. The new fuel rod state provides the initial state conditions for the next time step. The calculations are cycled in this manner for the number of time steps as specified by users. The response at each time step is
35
determined by repeated cycling through two nested loops of iterative calculations until the fuel- cladding gap temperature difference and internal gas pressure converge [7].
Figure 8: Simplified FRAPCON-3 flowchart [7].
FRAPCON was chosen as a primary tool for this work because its source code is publicly available so that it is possible to modify the code with additional physical models or material properties. Given its creditability as an independent audit tool in NRC’s reviews of industry fuel performance codes in the United States, FRAPCON should provide reasonably accurate results within the limitations of the code.
36
Given the fact that FRAPCON was originally developed for conventional LWR fuel designs, several physical phenomena at high burnup and high temperature regime are not modeled in the code. For example, given a regulatory burnup limit of 62 MWd/kgHM, a high burnup structure (HBS) phenomena that usually occurs beyond this burnup limit are not included. Fuel constituent transport models under temperature gradient are not covered. This is understandable because in typical LWR fuel rods where fuel temperatures are well below its melting point, these phenomena are negligible. Nevertheless, FRAPCON was determined to serve as a reliable platform for advanced fuel design in this work because the code has been extensively benchmarked and validated with experimental data [7]. After additional models that represent physical phenomena at high temperature and high burnup are included, it can credibly be used to accurately analyze the behavior of RBWR fuel rods.
1.4.2 FRAPCON-EP
FRAPCON-EP is a modified version of FRAPCON-3.3 where EP stands for enhanced performance [8] [9] [10]. It was developed to improve prediction capability of FRAPCON-3.3 beyond a regulatory burnup limit of 62 MWd/kgHM and to address physical phenomena that become important only at high temperature and high burnup. Similar to other versions, FRAPCON-3.3 was originally developed to simulate fuel rod behaviors in typical LWR conditions which operate at relatively low temperature i.e. fuel temperature is less than half of melting point. Because of this limitation, the code did not explicitly model migration behavior of fuel constituents such as uranium, plutonium, oxygen, porosity, and fission products which normally occur at high temperature under steep temperature gradients. In addition, fuel swelling from gaseous fission products was assumed negligible in FRAPCON-3.3 which may be the case under typical LWR operating conditions. Given a burnup limit of 62 MWd/HM, high burnup structure (HBS) formation which usually occurs beyond this burnup limit are not covered.
FRAPCON-EP introduces several physical models that represent time-dependent behavior of fuel constituents under steep temperature gradient and high temperature. From a uniform radial distribution of plutonium, oxygen and fission products, FRAPCN-EP employs mechanistic models based on Fick’s law of diffusion and thermal migration (Soret’s effect) to describe the evolution of these species under the effect of temperature. A change in radial distribution of
37
plutonium will directly impact the radial power profile of a fuel pellet and fuel thermal conductivity. Oxygen migration affects the oxygen-to-metal ratio (O/M) and fuel thermal conductivity locally. As a representative of volatile fission products, cesium migration was modeled in FRAPCON-EP because of its high mobility and tendency to cause localized fuel swelling. Fuel porosity migration was also modeled in FRAPCON-EP. Under high temperature gradient, as-fabricated porosity migrates to higher temperature region, concentrates at fuel center and eventually forms a central void. The formation of central void changes fuel geometry from solid pellet to an annular pellet which affects temperature distribution in the fuel.
For physical phenomena at high burnup, FRAPCON-EP includes HBS formation, O/M variation with burnup and acceleration in corrosion and hydriding. High burnup structure is a change in fuel microstructure from a dense crystal structure with grain size of around 10 microns to a porous structure with a significantly smaller grain size (0.01-0.1 microns). Usually occurs at fuel periphery where temperature is below 1273 K (1000 oC), the occurrence of HBS increase fuel porosity and degrade fuel thermal conductivity at high burnup. As a result of fission reaction which releases two oxygen atoms in oxide fuel and contamination of fission products in fuel matrix, O/M ratio increases with burnup. This phenomenon had been captured in FRAPCON-EP.
In zirconium alloy cladding, oxidation resistance comes from the existence of secondary phase particles (SPP) which are intermetallic precipitate of zirconium and alloying elements. However, these particles will become amorphous or completely dissolved in the fuel cladding matrix as free metallic particles at high burnup. Once a complete dissolution of SPP takes place, acceleration in corrosion and hydriding of cladding can occur. FRAPCON-EP uses neutron flux and irradiation time to estimate threshold neutron fluence in which SPPs are completely dissolved. After that, a higher oxidation rate is assumed which will result in higher concentration of hydrogen in the cladding.
Although, various physical phenomena at high temperature and high burnup have been modeled and described in FRAPCON-EP, there was still significant room for improvement. For example, several physical models such as porosity migration and cesium migration were based on empirical models. Fundamentally, the overall behaviors of these models are heavily depends on arbitrary empirical constants and their applicability is often limited to a specific set of empirical
38
data on which they were developed. To improve the generality of these models, a more mechanistic and general approach was used to describe these phenomena. Furthermore, a base version of FRAPCON-EP needed to be upgraded to a more recent version of FRAPCON so that the code can gain benefits from enhanced capability, updated material property correlations, and improved fuel behavior models.
1.5 Thesis organization
This thesis consists of seven chapters and two appendices covering background information, material property development, physical model implementation, experimental validation and sensitivity analysis on key design parameters.
Chapter 1 describes the main objectives and motivation of using advanced LWR for plutonium and transuranic waste incineration. A brief description of a fuel performance code FRAPCON and its variants are also given.
Chapter 2 describes a brief history of high-conversion LWR designs as plutonium breeder and transuranic waste burner including Hitachi RBWR designs. With emphasis on transuranic waste incineration, specific design of RBWR-TB2 is described in details about its characteristics and material challenges in term of fuel performance.
Chapter 3 presents a literature review of various thermal conductivity correlations for MOX fuels developed over the years. A benchmarking of these correlations is performed to highlight the effects of dependent variables on thermal conductivity values. Comparison between correlations and experimental data are performed in order to identify the most appropriate one for RBWR-TB2 analyses.
Chapter 4 presents physical phenomena relevant to RBWR-TB2 conditions which operate at high temperature and high burnup. Physical phenomena at high temperature such as fuel constituent migration and gaseous fuel swelling are covered. Degradation of material properties at high burnup including corrosion, high burnup structure formation are also discussed in this chapter.
39
Chapter 5 presents the validation effort of some of the physical models implemented in FRAPCON-3.5 EP. To compare the code with experimental data from fast reactors, additional code modification to better reflect fast reactor conditions such as sodium coolant and stainless steel cladding is discussed.
Chapter 6 includes results of fuel performance modeling of RBWR-TB2. Several key performance indicators in fuel modeling such as fuel temperature, oxide layer thickness, cladding stresses and strains are presented and discussed. Sensitivity study of fuel rod behavior relative to key design parameters including a different type of cladding material are included in this chapter.
Chapter 7 summarizes the observations and conclusions drawn from the work performed for the thesis. Further improvements and opportunities for future research are also discussed.
40
Chapter 2
High-conversion Light Water Reactors
2.1 History of HC-LWR development
The High-Conversion LWR (HC-LWR) is a collective term used to describe a type of advanced LWR design with a conversion ratio higher than currently operating LWRs. From a fuel cycle prospective, LWR serves in three different roles: fuel consumer, power producer, and waste producer. On the other hand, the HC-LWR was designed to produce power at the same level as an LWR but consume less fuel and generate less waste with multi-recycling or Transuranic Waste (TRU) burning capability. Being able to recycle its spent nuclear fuel indefinitely translates to a significant reduction in natural uranium ore and amount of waste to be stored in geological repositories, thereby, increasing fuel resource sustainability and at the same time reducing environmental burden from nuclear waste.
The HC-LWR can be categorized based on its fuel composition: U-Pu or U-Th cycles. Then, it can be further categorized by its coolant characteristics: pressurized water or boiling water. In addition, it is also possible to use heavy water (D2O) instead of light water (H2O) as coolant to attain high conversion ratio and multi-recycle capability. However, since it does not exactly fall under the definition of the HC-LWR design, heavy water cooled reactors are not included in this review. A chronological list of some of HC-LWR designs is shown in Figure 9 [11].
41
High-Conversion Light Water Reactor
Fuel composition Uranium-Plutonium Uranium-Thorium
Boiling Pressurized Boiling Pressurized Coolant Water Water Water Water
RMWR (Takeda 1995) B&W PWR HCBWR LWBR (Edlund 1975) (Downar 2001) (Connors 1979)
BARS (Toshiba 2001) KWU APWR RMPWR Seed Blanket (Lindley 2014) (Broeders 1985) (Nunez-Carrera 2008) FLWR (JAEA 2009) HGLWBR (Radkowsky 1988) RBWR-Th (Ganda, 2012) RBWR (Hitachi 2009) RCVS (Framatome 1988)
Big Mac (Ronen 1998)
RMR-PWR (JAERI 2003)
Figure 9: Classification of high-conversion LWR [11].
The development of the HC-LWR began roughly a decade after fast breeder reactors during the 1970s. At that time, nuclear power was rapidly expanding and there was an ongoing concern over scarcity of natural uranium resource. By the 1980s, advances in geological exploration and uranium mining technology proved that uranium resource was not going to be a problem at least until the end of 21st century although the long-term issue of spent nuclear fuel from LWRs still persist [12]. Fast breeder reactors—nuclear reactors that can produce more fuel than they consume—were specifically designed as a countermeasure to this problem [13]. These reactors need to be cooled by a non-moderating coolant such as sodium, helium, or lead. With limited experience with such designs and successful deployment of light water as coolant in the nuclear navy program [14], there were some interests to pursue a light water reactor design with higher
42
conversion ratio as an alternative development route [15]. HC-LWR achieves multi-recycling or TRU burning capability by reducing neutron moderation in the core; thereby, shifting the neutron spectrum toward higher energy. Even though various HC-LWR designs have been proposed over the years, most of them share this characteristic to some extent. Depending on the nature of coolant, different design strategies are used to harden the neutron spectrum.
High-Conversion Pressurized Water Reactor
In pressurized water systems, reduction in neutron moderation is generally achieved by displacing water volume with fuel volume. Edlund was the first to propose the idea of retrofitting an existing PWR with new fuel assembly design capable of achieving higher conversion ratio [16]. In this design, the fuel assemblies are arranged in hexagonal lattice so that they can pack more fuel per unit area than a standard square lattice. Fuel rod pitch and assembly pitch are reduced; thus displacing water volume even further. As a result, the conversion ratio is increased from 0.5 to 0.9. A Babcock and Wilcox PWR was used as a template for this reactor concept.
This reactor is designed to use uranium-plutonium mixed oxide (UO2-PuO2 or MOX) as fuel. A negative void reactivity coefficient was achieved based on neutronic calculations [17].
Broeders and Donne proposed a high conversion advanced PWR (APWR) based on a German Kraftwerk Union (KWU) PWR using MOX fuel in a hexagonal tight lattice geometry. In a homogeneous core design, fuel rod geometry is identical; only fissile plutonium content changes in the radial direction to smoothen the radial power distribution. In a heterogeneous core design, the reactor core is divided into two radial zones: seed and blanket. Seed area is defined by an
inner core area with MOX fuel whereas the blanket area is an outer core area with depleted UO2 [18]. With radial seed-blanket designs, the conversion ratio of KWU APWR increased from 0.90 in a homogeneous core design to 0.96 in the heterogeneous core design and maintained a negative void reactivity coefficient from increased neutron leakage from seed to blanket zones.
However, PWR designs with a tight lattice and simple heterogeneous core can only achieve conversion ratio slightly below 1.0 at best. This is because the presence of single-phase liquid water in PWRs makes it difficult to increase conversion ratios greater than 1.0 unless some
43
advanced features such as a spectral shift mechanism, multi-stage operation with rapid reprocessing or complex core heterogeneity are used.
Radkowsky and Shayer achieved a higher conversion ratio of 1.08 in their High Gain Light Water Breeder Reactor (HGLWBR) design [19]. HGLWBR uses multi-stage and multi-core operation with rapid fuel reprocessing to achieve the high conversion ratio under the condition of PWR. HGLWBR operates in two stages of operation using two reactor cores: pre-breeder and breeder. In the pre-breeding stage, spent fuel from an LWR is fed to pre-breeder to maximize Pu- 241 production. Pu-241 is needed to achieve breeding capability in HGLWBR because the η value (number of neutrons emitted per neutron absorbed) is higher for Pu-241 than Pu-239 for neutron energies below 0.1 MeV. Given short half-life of 14.4 years for Pu-241, the fuel reprocessing and fabrication time has to be reduced to 3 months between discharging from pre- breeders cores and loading into breeder cores. To further increase breeding potential, high- density metallic Pu-Zr fuel is used in seed region whereas depleted UO2 is used in blanket region.
Spectral shift mechanism was used in RCVS reactor to promote fuel breeding where RCVS stands for Réacteur Convertible à Variation de Spectre [20]. The core was designed with an inner seed region using either enriched UO2 or MOX fuel, surrounded by an axial and radial blanket
region of depleted UO2 and stainless steel reflector. Fuel assemblies were arranged in a hexagonal lattice with capability of changing neutron spectrum with burnup. Essentially, they are fuel assemblies with water holes and removable fuel rods inside. Initially, they are filled with
depleted UO2 to make use of excess reactivity for plutonium breeding during the beginning of cycle (BOC). As excess reactivity decreases with burnup, these removable fuel rods are gradually pulled out of the core for further reprocessing. The remaining water holes are then filled with water thus allowing moderator-to-fuel ratio to increase and neutron spectrum to shift toward thermal spectrum for better fissioning of fissile materials. This design achieved a conversion ratio of 0.95 with negative void reactivity coefficient.
An example of micro-heterogeneity configuration was proposed in the so-called “Big Mac” design which is a high conversion PWR that stacks alternating layers between MOX (seed) and
natural UO2 (blanket) within a fuel rod [21]. The core dimension is similar to typical PWRs – 3.6
44
m height and 3.37 m diameter. Each layer is 3 cm thick. Because of neutron leakage from seed to blanket, this arrangement promotes plutonium breeding in blanket layer and increases conversion ratio up to 0.92 with negative void reactivity coefficient.
The RMWR-PWR is a Japanese design of a high-conversion PWR as a variant in the RMWR series which mostly focuses on high-conversion BWRs [22]. This design achieves conversion ratio around 1.0 and negative void reactivity coefficient with the use of zirconium hydride
(ZrH1.7) rods to soften the neutron spectrum in the blanket zone [23]. The RMWR-PWR employs both micro-heterogeneity at both fuel rod and assembly level to promote fuel breeding. A fuel rod is stacked with axially alternating layers between MOX and depleted UO2. The fuel
assembly comprises of inner seed fuel rods (MOX) and outer blanket fuel rods (depleted UO2). The fuel assemblies are arranged in a hexagonal tight lattice to reduce neutron moderation. Stainless steel was used instead of zirconium because of the high linear heat rate and the allowance of a reduction in cladding thickness. Subsequent design iterations of the RMWR-PWR eliminated zirconium hydride rods in blanket zone through optimization of various parameters including the length of seed and blanket in a fuel rod and number of rings in seed and blanket in an assembly [24]. Cladding material was also changed from stainless steel to Zircaloy to improve neutron economy in the core.
The RMWR-PWR was recently examined by Andrews et al. [25] for further improvement of the conversion ratio. The goal of this study was to increase the conversion ratio beyond 1.0, thus, achieving fuel breeding capability in a pressurized water environment. A fuel assembly of RMWR-PWR as proposed by Shelley et al. [24] was used as a reference design. With the use of higher density fuel form such as uranium-plutonium mononitride, (U,Pu)N, and a replacement of water rods with voided rods, it was possible to further reduce the moderator-to-fuel ratio of the assembly and increase the conversion ratio from 1.0 to around 1.03 at a single batch discharge burnup of 35 MWd/kgHM. This work also investigated some safety parameters including thermal hydraulics and reactivity coefficients. Although, the moderator density coefficients in cases of nitride fuels were found to be positive, the coupled power coefficients were more negative and should ensure the safety of the reactor in off-normal conditions. The author also pointed out that the radially reflective boundary conditions which are typically used in assembly-
45 level neutronic calculations might produce somewhat more conservative results than a full-core analysis which considered neutron leakage in radial directions.
It is worth to mention the Shippingport Light Water Breeder Reactor (LWBR) as a proven demonstration of breeding capability in a pressurized light water environment [15]. It is a high- conversion PWR design based on uranium-plutonium cycle in the first and second cores and uranium-thorium fuel cycle in the third and final core. Closely-packed hexagonal lattice and seed-blanket assembly design again were used to increase the conversion ratio. In the first two cores, fuel composition was made of highly enriched uranium oxide (93% U-235) as fissile material in the seed zone and natural uranium oxide as fertile material in the blanket zone. The last core was made of fissile uranium oxide (U-233) as seed and natural thorium oxide as blanket. It was the only high-conversion LWR that progressed through conceptual design, construction, operation and decommission – a full life cycle. Fuel breeding capability was proven in the final core of the LWBR with UO2-ThO2 fuel; approximately 1.4% more fissile materials were found in the core after it had been operating for 5 years during the last stage of its operation from 1977 to 1982 [26].
Achieving a conversion ratio greater than 1.0 is less difficult in uranium-thorium fuel than uranium plutonium, because of the neutronic characteristic of U-233 as a fissile nuclide. Naturally, the η value of U-233 is greater than 2 in a wider energy range than that of other fissile nuclides, making it possible for fuel breeding even in a thermal spectrum.
A more recent design of high-conversion PWR featuring thorium fuel was proposed by Lindley and Parks under the name of reduced-moderation PWR (RMPWR) [27] [28]. RMPWR was designed for TRU incineration instead of fuel breeding or fuel self-sustainability through multi- recycling of Th-TRU fuel in a pressurized water environment. One unique feature of this design is the range of fuel isotopic composition which is considerably more complex than U-TRU or U- Pu fuel. In oxide form, the fuel composition of RMPWR comprises Th-TRU fuel where TRUs come from LWR spent fuel and Th-U fuel where various isotopes of uranium (U-233, U-234, and U-236) are generated through neutron absorption in thorium. RMPWR was designed to be compatible with current generation of PWR such as EPR or AP1000 [28]; therefore, the fuel assembly is arranged in square lattice with similar fuel rod pitch as typical PWR. Neutron
46
moderation is reduced by increasing fuel rod diameter so that the moderator-to-fuel ratio was reduced from 1.98 in a reference PWR design to 1.09 in the RMPWR. Negative void reactivity coefficient is achieved through multiple zoning of fuel composition. In heterogeneous fuel assembly design, Th-U fuel rods are placed in the central region and Th-TRU fuel rods in the assembly periphery. In heterogeneous core design, fuel assemblies of entirely Th-U and Th-TRU rods are arranged in a perfect checkerboard lattice. Depending on fuel design, the TRU incineration rate of RMPWR ranges from approximately 170-190 kg per year per 1 GWth reactor. Axial heterogeneity is not incorporated in the RMPWR. Unlike breeder or break-even reactors where the use of thorium fuel increases the fuel breeding potential because of the η
value of U-233, for TRU burning application, ThO2-based fuel offers little advantage over UO2- based fuel because of complexity in fuel composition [29].
High-Conversion Boiling Water Reactor
In terms of fuel breeding, the boiling water reactor (BWR) has some advantages over pressurized water reactors (PWR) because it operates with the presence of boiling and two-phase coolant. Naturally, a two-phase mixture of liquid water and steam has lower density than single-phase liquid water. Therefore, the neutron moderation somewhat diminishes toward the top of reactor core. In traditional BWRs, it is necessary to include water rods – empty fuel rods without end caps – to serve as water flow channels and allow more liquid water to reach through the top without boiling and improve neutron moderation around the vicinity. In-core boiling and lower neutron moderation makes it relatively easier to increase conversion ratio in a BWR. Unlike a PWR which relies only on water volume reduction, in BWRs, both water volume and water density can be manipulated to achieve the objectives of hardened neutron spectrum and increasing conversion ratio. The concept of a high-conversion BWR received particular attention in Japan because it is more suitable to the country’s situation where BWRs are a majority of its reactor fleet.
This review includes three high-conversion BWR designs from Toshiba, Japan Atomic Energy Agency (JAEA) and Hitachi. Although, they share many similarities in designs features such as a shorter core, closely-packed fuel lattice, and high void fraction in the core, some differences in
47
fuel composition, geometry, and functionality do exist. Most of these designs are based on uranium-plutonium cycle thus the fuel composition is made of either enriched uranium or MOX.
Toshiba proposed a high-conversion BWR and named it BWR with Advanced Recycle System (BARS) [30]. The design was based on an existing design of advanced boiling water reactor (ABWR) – the latest generation of currently operating BWR – to reduce capital and development cost. Unlike any other high-conversion LWR, BARS uses a tightly-packed triangular fuel lattice arranged in a square fuel assembly. Core height was reduced from a full-length at 3.7 m in the ABWR design to 1.6 m in a normal fuel assembly and 0.8 m in a partial fuel assembly. BARS features a seed-blanket design with full-length and half-length fuel assemblies. The partial-length fuel assembly serves as a neutron-streaming channel in the upper half of the core to increase neutron leakage thus achieving a negative void reactivity coefficient through the cycle. The average void fraction of BARS is about 0.6 leading to a conversion ratio of 1.04. Another technique to improve the conversion ratio is to reduce cladding thickness so that fuel volume is increased but sufficient cooling is preserved. However, it comes with some penalty in term of neutron economy because stainless steel cladding is required given reduced strength because of reduction in thickness. To flatten power peaking, the fuel assembly of BARS contained different zones of plutonium weight fraction: 11% in low enrichment zone and 14-17% in high enrichment zone. The latest design effort has included thermal hydraulic analyses [31] [32] [33].
JAEA applied the concept of multi-role reactors in their high-conversion BWR design called Flexible Light Water Reactor (FLWR). It is also designed based on the ABWR with tight hexagonal lattice and seed-and-blanket MOX fuel design and increased average void fraction in the core. FLWR cores can be configured into different roles in a fuel cycle: TRU burner or plutonium breeder. The TRU burning version of FLWR is called high-conversion FLWR (HC- FLWR) and plutonium breeder version is called Reduced-Moderation Water Reactor (RMWR). In HC-FLWR, the conversion ratio is 0.84 with capability to use processed fuel that contains minor actinides [34]. HC-LWR fuel rods are segregated into three zones: lower blanket zone at the bottom, seed zone at the middle and another blanket zone at the top.
RMWR is a fuel breeder version of FLWR with a conversion ratio of 1.04. In the RMWR, the moderator-to-fuel ratio is reduced even further than the HC-FLWR by reducing fuel rod gap and
48
increasing fuel rod diameter. In addition, RMWR fuel rods have five zones. From the bottom to the top, there are lower blanket, lower seed, middle blanket, upper seed and upper blanket zones. These two reactor cores share some similarities in size and shape of fuel assembly so they can be converted into one another. This feature allows some flexibility in the fuel cycle if future needs arise. To maintain necessary neutron leakage and negative void reactivity coefficient, the core height of these reactors are significantly shorter than original length of ABWR: 0.95 m in HC- FLWR and 1.255 m in RMWR.
Hitachi also designed a series of high-conversion BWRs under the development name of Resource-Renewable Boiling Water Reactor (RBWR). This reactor concept was first proposed by Takeda et al. in 1995 [35] and is still under active development to date [36]. Similar to previous high-conversion BWRs described above, the reactor core of the ABWR was used as a design template and then several techniques were incorporated to reduce moderator-to-fuel ratio and harden the neutron spectrum. Hitachi also conceptualized RBWRs as multi-role reactors in a multi-stage fuel cycle; therefore, several design variants of RBWR were proposed. Classified by their roles in the fuel cycle, there are 5 different variants of RBWR to date: RBWR-T3, RBWR- AC, RBWR-Th, RBWR-TB, and RBWR-TB2. Each core design has different geometry, composition and purposes but they share the similar goal of being a better alternative to Sodium Fast Reactor (SFR), whether in the role of initiator, breeder reactor or TRU burner. As an initiator, RBWR-T3 uses enriched UO2 to operate and is designed to maximize plutonium production left in the spent fuel. The RBWR-AC is designed to be a self-sustaining reactor, where it breeds fuel as much as it consumes during its lifetime. The fuel of RBWR-AC is composed of TRUs extracted from spent nuclear fuels in depleted uranium matrix. At equilibrium, the fuel composition at the end of cycle will be approximately the same as at the beginning of cycle. On the other hand, RBWR-TB uses a slightly different geometry and composition to maximize the fission of TRUs [11].
Regardless of the design purposes, all RBWRs share the same technical challenge in that H/HM ratio (the ratio between moderator and fuel) has to be decreased to reduce the moderation in the reactor core so that they can breed fuels or incinerate TRU wastes via high-energy neutrons. Several design techniques have been used to harden the neutron spectrum and promote fuel breeding: (1) reducing water density by increased boiling, (2) reducing water volume with tight
49 lattice in a hexagonal geometry, and (3) increasing neutron absorption in a blanket region with axial heterogeneity by alternating blanket and seed regions in the fuel rod [11]. Figure 10 compares moderator-to-fuel ratio and breeding ratio of conventional BWR, RBWR-AC, RBWR- TB and RBWR-TB2 [36]. It can be seen that as the moderator-to-fuel ratio decreases; breeding ratio increases accordingly. With breeding ratio close to 1.0, RBWR-AC and RBWR-TB can sustain a continuous cycle without the need for external fissile material. Having lower breeding ratio by design, RBWR-TB2 requires external fissile TRUs from LWR spent fuels.
Figure 10: Moderator-to-fuel ratio and breeding ratio comparison [36].
The RBWR-T3 serves as an initiator in a closed-loop nuclear fuel cycle with the purpose of maximizing plutonium production from natural uranium without multi-recycling capability [11]. The spent fuel of the RBWR-T3 can be fed to RBWR-AC where multi-recycling of its own fuel is possible. Given its objective of once-through plutonium breeding, the neutron spectrum of RBWR-T3 is more thermalized than any other versions of the RBWR. The fuel composition of the RBWR-T3 is enriched UO2 in a hexagonal lattice. To allow more neutron moderation, the fuel rod pitch and inter-assembly pitch are wider than in other RBWR designs. The average void
50
fraction is similar to current ABWRs. In addition, water rods are included in the center of RBWR-T3 fuel assembly to increase neutron moderation in the upper region of the core.
The RBWR-AC is a self-sustaining version of the RBWR [36]. In other words, the RBWR-AC can perpetually operate with its own spent fuel except for the first fuel loading which requires an external source of fissile materials. The first batch of RBWR-AC fuel can be loaded with TRUs extracted from spent fuel from a LWR or RBWR-T3 to initiate an indefinite fuel recycling. Several design techniques are used to the harden neutron spectrum in the RBWR-AC in order to achieve self-sustainability or conversion ratio slightly above 1.0 and, at the same time, maintain a negative void reactivity coefficient. The RBWR-AC fuel assembly is arranged in a tightly- packed hexagonal lattice with small pitch-to-diameter ratio and narrow inter-assembly gap. Core flow rate is reduced by a factor of 2.5 when compared to that of the original ABWR; this makes core-averaged void fraction increase because of additional boiling. RBWR-AC fuel rods have 5 axial zones. From the bottom to the top, there are lower blanket, lower fissile, middle blanket, upper fissile and upper blanket zones. Seed-and-blanket design in conjunction with a much short core in RBWR-AC results in higher neutron leakage from the fissile zone which helps mitigate reactivity insertion upon core voiding. Fuel compositions in the blanket zone are made of
depleted uranium in the form of oxide (UO2) whereas in fissile zones they are made of mixed oxides of uranium, plutonium and minor actinides. Multiple plutonium enrichment zones reduce radial power peaking within a fuel assembly.
The RBWR-Th is a companion design of RBWR-AC with different fuel composition and axial seed-blanket zoning [28] [37] [38] [39]. As the name suggests, the RBWR-Th uses ThO2 instead of depleted UO2 in both fissile and blanket zones. It was designed to be a fuel-self-sustaining reactor which operates on the U-Th fuel cycle with U-233 as fissile material. Common design characteristics for high-conversion BWRs are also found in RBWR-Th: tight hexagonal lattice, high void fraction and high exit steam quality. However, the core height of RBWR-Th is not as short as any other designs. In fact, the length of the fuel rod is similar to the original ABWR at around 3.8 m whereas it is around 1.4 m in RBWR-AC design. In the RBWR-Th, there are 3 axial zonings: lower blanket, middle fissile zone and upper blanket. Middle blanket zone was eliminated; the upper and lower fissile zones were combined and elongated into a middle fissile zone. The RBWR-Th also achieves fuel sustainability with conversion ratio slightly above 1.0,
51
meaning that only externally supplied natural ThO2 is needed to feed the reactor. Reactivity coefficients for fuel temperature, coolant void, and power are all negative in this design.
The RBWR-TB is a TRU burning version of the RBWR [36]. The primary goal of this reactor is to completely eliminate the total amount of TRU through multi-recycling. Initiation of the RBWR-TB cycle can be made by loading TRU materials from an RBWR-AC. After that, it can consume its own spent fuel. Unlike the RBWR-AC which was designed to operate indefinitely, RBWR-TB was designed to operate for a limited amount of time until TRUs are completely incinerated. Performance of the RBWR-TB is defined by fission efficiency instead of a conversion ratio where fission efficiency is defined as the net decrease in TRUs divided by the total amount of actinides at the end of cycle. Currently, the RBWR-TB is able to attain fission efficiency of 51% meaning that at discharge it can burn half of TRUs loaded at the beginning of cycle. Therefore, the amount of TRUs can be reduced from 100% down to less than 1% after 8 cycles of operation in the RBWR-TB. Given that the total fuel residence time in the RBWR is approximately 4 years; it would take less than 40 years to completely incinerate the TRUs left from the fuel cycle. This may be viewed as a counter-productive measure if we are to sustain nuclear power generation; however, in the case of a nuclear phase-out scenario, RBWR-TB would significantly reduce the future burden of safeguarding TRU materials which would remain radioactive for hundreds of thousands of years. The core design of the RBWR-TB differs from the RBWR-AC in that it has a shorter core height, smaller fuel rod diameter, smaller rod pitch, and larger inter-assembly gap. The height of the fissile and blanket zones are also different from RBWR-AC. The lower blanket zone is eliminated. All of these adjustments are required to enable the RBWR-TB to recycle of TRUs multiple times under a different in neutron spectrum.
The RBWR-TB2 is a modified version of RBWR-TB proposed by the Electric Power Research Institute (EPRI) [40]. Unlike RBWR-TB which was specifically designed for a nuclear phase-out scenario, the purpose of RBWR-TB2 is to reduce of the amount of TRUs from spent fuel of existing LWRs and to sustain nuclear power generation. The design of the RBWR-TB2 is somewhat similar to that of RBWR-TB except that it requires additional fissile plutonium enrichment to compensate for the reactivity penalty of LWR spent fuels. Fission efficiency of RBWR-TB2 is currently rated at 45% whereas TRU production efficiency is 22% in ABWRs. Fission efficiency is defined as a net decrease in TRUs divided by total amount of fissioned
52
actinides through the total fuel residence time in the core. The TRU production efficiency is an inverse quantity of the fission efficiency; it is a net increase in TRU divided by total amount of fissioned actinide at the end of cycle. Essentially, these two values indicate that the amount of TRUs produced from two units of ABWR would be suppressed by one unit RBWR-TB2.
Figure 11 compares the formation rate of TRUs (both fissile and non-fissile isotopes) in metric tons per year per reactor of conventional BWRs with that of RBWR-TB and RBWR-TB2. For conventional BWRs, these values are positive. For RBWR-TB and RBWR-TB2, negative values in both fissile and non-fissile TRU isotopes means there is a net decrease in their quantity, every operation cycle [41].
Figure 11: Rate of formation and consumption of TRUs [41].
For the application of TRU incineration in advanced LWRs, this thesis specifically focuses on RBWR-TB2 design which recently received more attention than any other variants of RBWR [42] [43] [44] [45] [46]. More details about its design, characteristics and challenges in term of fuel performance will be discussed in the next section.
53
2.2 General description and design characteristics of RBWR-TB2
To reduce capital and development cost, the RBWR was designed with a commonality concept in that, except the reactor core, all other components in the power plant such as steam separator, steam dryer, steam turbine, condenser, feedwater pumps are the same as currently operating ABWRs. Even the reactor core of the RBWR was designed to fit within reactor pressure vessels of existing ABWRs. The common plant specifications of RBWR and ABWR are listed in Table 2 [42]. At the plant scale, both reactors are designed to deliver the same amount of power in the same reactor pressure vessel (RPV) at the same operating pressure.
Table 2: Plant specification and core design [42]. Reactor RBWR ABWR Thermal Power (MW) 3,926 3,926 Electrical Power (MW) 1,356 1,356 RPV Diameter 7.1 7.1 Core pressure 7.2 7.2 Number of fuel bundles 720 872 Fuel lattice type Hexagonal Square Lattice pitch (mm) 199 155 Number of control rods 223 205 Control rod type Y-type Cross-shape
Figure 12(a) shows the reactor pressure vessel of the RBWR [42]. A core configuration in horizontal cross-section view is shown in Figure 12(b) [42]. It can be seen that the RBWR reactor core is composed of 720 fuel bundles, and 223 control rods which is common to all design variants of RBWR.
54
(a) (b) Figure 12: (a) Reactor pressure vessel of RBWR (b) Horizontal cross-section of RBWR reactor core [42].
RBWR core is designed as a parfait core concept as the axial configuration of RBWR fuel assembly employs axial segregation between seed and blanket regions to enhance fuel breeding and/or promote transmutation reactions of TRUs. As an example of this concept, an RBWR-TB2 fuel bundle is shown in Figure 13 [42]. In this design, an internal blanket of depleted uranium is placed between lower and upper fissile zones. Then the upper and lower blankets are attached the upper and lower fissile zones, respectively. This parfait design is necessary to keep the void reactivity coefficient negative when the axial power changes due to the change of void fraction during transients. Furthermore, neutron absorber zones are attached above and below the fuel zones to increase negativity of void reactivity coefficient. The upper neutron absorber rod is a
sealed tube filled with B4C pellets while the lower neutron absorber zone is filled with B4C pellets at the bottom of the fuel rod. So, in each fuel rod, there are six zones: (1) lower neutron absorber, (2) lower fissile, (3) internal blanket, (4) upper fissile, (5) upper blanket, and (6) plenum.
55
Figure 13: Axial and hexagonal configuration of RBWR-TB2 fuel bundle [42].
Figure 14 shows a horizontal cross-section view of the fissile zones of the RBWR-TB2 fuel bundle [42]. For this specific design, the bundle-average plutonium content as expressed in term of plutonium weight fraction is approximately 80% and 70% in the lower and upper fissile zones, respectively.
Figure 14: Horizontal configuration of RBWR-TB2 [42].
56
Because of its unique pattern of axial segregation between seed and blanket regions, considerable variation in linear heat generation rate (LHGR) between each zone exists. The axial power distribution in the RBWR-TB2 exhibits multiple peak spots and a very steep gradient at the fissile-blanket interfaces as shown in Figure 15 [47]. It can be seen that the LHGR in the upper fissile zone could be as high as 35 kW/m at the middle-of-life and stay above 25 kW/m throughout the cycle. This is quite a significant departure from typical operating condition of LWRs which is around 15-25 kW/m [48].
35 BOL MOL 30 EOL
25
20
15 Local LHGR (kW/m) Local LHGR 10
5
0 0 20406080100120 Axial Length (cm) Figure 15: Axial LHGR of the RBWR-TB2 as a function of core height [47].
Because of these peaking regions, local burnup in fissile zones are extremely high. It can be seen from Figure 16 that local fuel burnup at end-of-life are on the order of 120 and 160 MWd/kgHM in the lower and upper fissile zones, respectively. This ultra-high level of burnup has never been achieved before in LWRs [49] [50] and could pose significant challenges in term of material property degradation. Double peaking axial power distribution profiles and ultra-high local fuel burnup are one of the most important design departures of RBWRs from BWRs.
57
180 BOL 160 MOL EOL 140
120
100
80
60 Local Burnup (MWd/kgHM) 40
20
0 0 20 40 60 80 100 120 Time (Days) Figure 16: Local fuel burnup of the RBWR-TB2 as a function of axial node and time step.
As result of water volume reduction due to the tight hexagonal fuel lattice and increased boiling by reduction of coolant mass flow rate, core-averaged void fraction of the RBWR-TB2 is 56% whereas that of an ABWR is 36%. Relatively speaking, this is a 55% increase from current operating condition of an ABWR as shown in Figure 17 [51]. The void fraction gradient in RBWR-TB2 also differs from that of the ABWR due to axial heterogeneity configuration. It can be noticed that the void fraction rises significantly faster in the lower fissile zone of RBWR-TB2. Then it tends to decrease in the internal blanket zone. Upon reaching the upper fissile zone, the void fraction rises again to reach 80% void fraction at core exit. On the other hand, the void fraction in the ABWR steadily increases from core inlet to exit. Although the core flow rate of RBWR-TB2 is reduced by a factor of 2.5 when compared to ABWR, with significantly higher core exit quality and void fraction, steam flow rate to turbine can be maintained at a similar level of an ABWR.
58
Figure 17: Axial void fraction of RBWR-TB2 as a function of relative core height [51].
A comparison of the neutron spectra of the RBWR-AC, RBWR-Th, ABWR as a representative of LWR and an SFR as a representative of fast breeder reactor is given in Figure 18 [52]. The neutron spectra of the RBWR-TB and RBWR-TB2 are not included here because they are not publically available yet. However, they are expected to be quite similar to that of RBWR-AC and RBWR-Th. As can be seen from the figure, the ABWR shows a typical neutron spectrum for a water-cooled reactor—double peaks in neutron flux at thermal (~0.1 eV) and fast (~1 MeV) energy regions. The neutron spectrum of the SFR also exhibits a usual behavior of fast reactors—a mountain-like spectrum with a single peak at around 0.5-1 MeV with higher neutron flux intensity than LWRs. With several design techniques to reduce moderation, it can be clearly seen that neutron thermalization and thermal flux peaks no longer exist in RBWRs. In fact, their neutron spectra look a lot like that of fast reactors. Furthermore, both types of RBWRs have a significantly higher portion of neutron at energy above 1 MeV and in epithermal energy than the SFR making it possible to breed sufficient fissile materials to achieve a conversion ratio of 1.0 and can perpetually recycle plutonium and minor actinides. With this hardened neutron spectrum, it is also possible to burn TRU wastes in RBWR-TB2 design.
59
Figure 18: Comparison of normalized neutron spectra for the RBWR-AC, RBWR-Th, SFR and ABWR [52].
2.3 Material challenges associated with RBWR-TB2
Although the RBWR-TB2 was designed based on the proven technology of the ABWR, its unique characteristics necessitate further considerations especially in terms of thermo- mechanical behaviors throughout fuel cycle. As can be clearly seen from Figure 19, the axial peaking factor of RBWR-TB2 is extremely high and discontinuous especially at fuel-blanket interfaces in order to keep core thermal output the same as with the ABWR with much shorter fuel length (around 1 m in RBWR-TB2 and 3.7 m in ABWR). In fact, the active fuel zone of RBWR-TB2 is even shorter than 1 m because it has to accommodate the blanket region which roughly takes approximately half of the total fuel length. As a result, the axial peaking factor of the RBWR-TB2 peaks at around 2.5-2.6 in the upper fissile zone at the beginning-of-life (BOL) and middle-of-life (MOL). These peaks tend to decrease as burnup increases and the blanket zone takes part in power generation toward the end-of-life (EOL). On the contrary, the axial peaking factor of the ABWR is considerably smoother with lower-magnitude peaks at around 1.5-1.6—a typical values of LWRs. Significant LHGR variation at the fuel-blanket interfaces would directly result in a very steep temperature gradient and it would likely cause additional thermal stress in the cladding due to differential thermal expansion.
60
4 ABWR BOL 3.5 ABWR MOL ABWR EOL RBWR-TB2 BOL 3 RBWR-TB2 MOL RBWR-TB2 EOL 2.5
2
1.5
1
0.5
0 0 0.5 1 1.5 2 2.5 3 Axial peaking factor Figure 19: Axial peaking factor vs. core height of ABWR and RBWR-TB2.
When compared to the ABWR, RBWR-TB2 fuel was designed to operate at higher peak burnup, higher linear heat generation rate (LHGR), and higher fast neutron fluence. Operating at higher LHGR may result in significant increase in fuel temperature and create several undesirable effects such as increased fuel swelling, fuel thermal expansion and fission gas release. High local burnup in the fissile zones could impact the degradation of material properties such as fuel density, fuel thermal conductivity which is heavily affected by fission product contamination, oxygen-to-metal (O/M) ratio, porosity, and microstructural changes. Therefore, several physical phenomena at high temperature and high burnup regime such as fuel restructuring, fuel constituent redistribution, and fission product migration are expected to occur in RBWR-TB2.
Figure 20 schematically shows the effect of excessive fuel swelling and thermal stresses on fuel and cladding deformation. From a uniform cylindrical shape, thermal stress causes a fuel pellet to crack and form an hourglass shape. During operation, fuel swelling will cause the fuel and cladding to be in direct contact when the gap between fuel and cladding closes at high burnup. At
61 this stage, the cladding will deform into a shape reflecting that of the fuel pellet making the cladding to look like a bamboo stalk. Direct contact between fuel and cladding adversely affects cladding integrity because it increases mechanical load and the risk of cladding failures from chemical interactions with corrosive fission products such as iodine and cadmium which will eventually lead to stress-corrosion cracking [53]. The failure from hard contact between fuel and cladding is termed pellet-cladding interaction (PCI). An example of cladding failure from stress corrosion cracking from hard contact and missing pellet fragment is shown in Figure 21.
Figure 20: Effect of fuel swelling and thermal stress [53].
Figure 21: Consequence of pellet-cladding interaction (PCI) [53].
In addition, the fast neutron flux of RBWR-TB2 is higher than a conventional BWR because of harder neutron spectrum; therefore, at the same burnup, irradiation damage in structural
62
components such as fuel cladding, control rods, channel box from fast neutrons is expected to be higher. Particularly for fuel cladding, higher fast neutron fluence may lead to acceleration of corrosion and hydriding in zirconium alloy cladding after a complete dissolution of secondary phase particles (SPPs) after neutron fluence threshold is exceeded [8] [9] [54]. The presence of zirconium hydride severely degrades mechanical strength of metallic zirconium cladding and results in a loss of ductility. As water flow through the core and is exposed to neutron and gamma irradiation, a wide variety of radiolysis products are generated. The most commons
species are H2 and H2O2. H2O2 can then decompose to O2 and H2O and remain dissolved in the
water. The presence of oxidants like H2O2 and dissolved O2 in the water increases the corrosion potential of the cladding. Naturally, free hydrogen radicals produced from radiolysis can recombine with the dissolved oxygen, become water molecules, and buffer the corrosion potential. However, this dissolved hydrogen can easily be lost once the water turns into steam
near top of the core as the H2 concentration in the water drops dramatically after boiling occurs. With increased boiling, the water chemistry of RBWR-TB2 is expected to be somewhat different
than that of the ABWR. In this case, the concentration of dissolved H2 in the core of RBWR-TB2 should be lower from higher core exit quality and core-averaged void fraction. In addition, the
effect of radiolysis, the concentration of dissolved O2, and the corrosion potential of water in the core of RBWR-TB2 are expected to be higher than in ABWR because of an increase in fast neutron flux.
Therefore, to accurately model the behavior of fuel rods under the conditions of RBWR-TB2, these design characteristics and physical phenomena at high temperature and high burnup have to be taken in account when performing thermomechanical analysis of fuel behavior. All of these phenomena relevant to RBWR-TB2 are modeled and discussed in Chapter 4.
63
Chapter 3
Evaluation of Thermal Conductivity Correlation for Mixed Oxide Fuels
Thermal conductivity is one of the most important material properties for fuel performance modeling because it dictates the fuel temperature distribution within fuel pellets. Fuel temperature relates directly or indirectly to all life-limiting parameters such as fission gas release, fuel melting, fuel swelling, and pellet-cladding interactions (PCI). Therefore, accurate prediction of thermal conductivity is a key element to achieve reliable simulation results for the thermo-mechanical behavior of fuel rods.
The objective of this work is to find an appropriate correlation for MOX fuel for use in the RBWR-TB2 which likely deviates from the applicability of default correlations in FRAPCON- 3.5 in terms of fuel burnup and plutonium weight fraction effects.
FRAPCON-3.5 currently has three options for thermal conductivity correlations: (1) Modified
NFI for UO2 and (U,Gd)O2, (2) Duriez-Modified NFI for MOX and (3) Halden for UO2,
(U,Gd)O2 and MOX. The Halden correlation is provided as an option while the other two are default models. These correlations are recommended up to a rod-average burnup of 62 MWd/kgHM which corresponds to a peak local burnup of around 150 MWd/kgHM at the pellet edge in LWRs. However, RBWR-TB2 fuel rods are designed to reach a higher rod-average burnup of roughly 70 MWd/kgHM with a peak local burnup of around 160 MWd/kgHM in fissile zones. Therefore, the applicability of these correlations needs to be re-evaluated using more recent experimental data.
For the MOX fuel option in FRAPCON-3.5, the Duriez-Modified NFI correlation is applicable for plutonium content from 3 to 15% weight fraction. According to Duriez [55], the effect of plutonium content on thermal conductivity in this range is negligible. In other words, thermal
conductivities of MOX at 3 to 15 wt% PuO2 are the same in Duriez-Modified NFI correlation.
However, the RBWR-TB2 was designed to use (U-TRU)O2 mixed oxide fuel which comprises
64
20% depleted UO2, 70% PuO2, and 10% Minor Actinides (NpO2, AmO2, CmO2, etc.) by weight. The effect of higher weight fraction of plutonium with addition of minor actinides should thus be carefully examined.
In this section, the brief history of MOX thermal conductivity correlations during the last three decades is presented. Emphasis is given on the selection and evaluation of default correlations from earlier to recent versions of FRAPCON. Then, we will compare a default correlation for MOX in FRAPCON-3.5, Duriez-Modified NFI, with other alternatives from the open literature.
Generally, development of thermal conductivity correlations can be categorized by type of samples: fresh and irradiated fuels and by type of measurements: in-pile fuel temperature and out-of-pile thermal diffusivity. Irradiated fuels can be measured both in-pile and out-of-pile but fresh fuels have to be measured out-of-pile to avoid irradiation effects. Direct out-of-pile measurement tends to be more reliable because important parameters can be controlled individually whereas reproducibility for in-pile measurement cannot be fully ensured. From out- of-pile measurement, thermal conductivity is derived from thermal diffusivity, density and heat capacity which are well-defined in the laboratory environment. The evaluation of thermal conductivity from in-pile temperature measurement is less accurate because of its integral nature in the reactor environment. To calculate thermal conductivity from in-reactor centerline temperature measurement, one has to solve a radial heat transfer equation from fuel pellet center to cladding outer surface which involves some uncertainties from estimation of fuel-cladding gap conductance.
3.1 Parameters affecting thermal conductivity
Fuel thermal conductivity depends on various parameters acting at different length scales. Most of them dynamically change during in-pile irradiation. At the macroscopic scale, temperature and burnup are two dominant factors because they generally represent the state of the fuel. High temperature represents the state of disorder in the crystal lattice and triggers another heat transfer mechanism through electron conduction. Fuel burnup represents the state of impurities and microstructural damage of the fuel. It can be viewed as a collective term which implicitly includes various effects from irradiation such as fission product contamination and
65 microstructural damage. Microscopic parameters such as stoichiometry (oxygen-to-metal ratio), and additive content (Pu, Gd, Cr) also play a critical role in determining thermal conductivity.
For non-conducting ceramics such as UO2 and MOX, heat is primarily transferred though atomic vibration. In this case, any microscopic perturbation to perfect crystal lattice would contribute to thermal conductivity degradation. Figure 22 illustrates parameters affecting thermal conductivity at microscopic scales where it can be further subdivided into mesoscopic and atomic level contributions.
Figure 22: Parameters affecting thermal conductivity [56].
Mesoscopic-level parameters involve imperfections in the microstructure of the oxides such as porosity, fission product precipitates, and intra-granular and inter-granular bubbles of fission gases. For fission products that are insoluble in the fuel matrix, they exist in the form of oxide precipitates, metallic inclusions, and inert gas atoms. Their impacts on thermal conductivity can be estimated by composite material correlations such as the Maxwell-Eucken or Loeb formulae [57] [58].
Atomic-level parameters are point defects including vacancies, interstitials, and substitutions of impurity atoms in the lattice such as Pu, Gd and fission products. Both point defects and substitutional atoms induce static displacements of host atoms from their equilibrium lattice sites. They serve as phonon scattering centers due to the differences in inter-atomic bonding
66 potential, ionic radii, and atomic mass, thereby limiting the phonon mean free path and lowering the rate of heat transport. Their effects are usually modeled in a form of k = 1/(A+BT) where k is thermal conductivity, A is a collective term representing thermal resistance due to phonon scattering by point defects and substitutional atoms and BT represents intrinsic lattice thermal resistance from phonon-phonon scattering. Stoichiometry or oxygen-to-metal ratio (O/M) has a major impact on overall thermal conductivity because it represents the state of oxygen vacancies in hypo-stoichiometric fuel (O/M < 2.0) or oxygen interstitials in hyper-stoichiometric fuel where (O/M > 2.0.).
During operation, fuel chemical composition, lattice structure and microstructure change markedly under intense neutron irradiation and high temperature. Because of their complexity and inter-correlation, it is difficult to isolate and model their effects individually. In practice, a number of parameters in thermal conductivity correlations are often reduced by retaining only the most influential parameters in the correlations while the effect of other parameters are implicitly taken into account by adjusting the models to experimental results.
In general, thermal conductivity of oxide fuels can be expressed as
1 k= +Ce (1) A+BT where the first term 1/(A+BT) represents lattice conduction by phonons and the second term CeDT represents high-temperature conduction through electron pair mobility. At temperatures below 1600 K, electron conduction is negligible. However, it becomes an important mode of heat transfer above 2000 K. Electronic contribution to thermal conductivity at high temperature is an important feature of UO2 because it does not exist in ThO2 and PuO2. This is because the presence of 5f valence electron at -2.3 eV in UO2 helps narrowing the gap between valence and conduction bands [59]. As shown in Figure 23, both ThO2 and PuO2 have larger valence- conduction band gaps at -6 eV [59]. This large energy barrier makes it very difficult for valence electrons in ThO2 and PuO2 to jump over the band gaps to the conduction bands and transfer the heat even at high temperature. From this information, ThO2 and PuO2 are classified as pure electronic insulators while UO2 can be viewed as semi-conductor at high temperature.
67
Figure 23: Electron density of state of UO2 (a), ThO2 (b), and PuO2 (c) [59].
Fuel burnup introduces several atomic perturbations such as soluble fission products, non- stoichiometry, interstitials, vacancies, and substitutions of foreign atoms. In general, thermal conductivity decreases with increasing burnup. To take into account the effect of burnup into the phonon conductivity term 1/(A+BT), there are two approaches: (1) modifying the empirical constants A and BT and (2) formulating multiplying factors based on fuel burnup. In the first approach, linear or quadratic equations may be used to capture burnup dependency in the phonon-defect (A) and phonon-phonon (BT) terms. For example, from constant values of fresh fuel, they may be changed to a burnup-dependent equation as A0+A1*bu and (B0+B1*bu)T where bu is fuel burnup, T is fuel temperature and A0, A1, B0, and B1 are new empirical constants derived from irradiated data. In the second approach, multiplying factors are formulated based on irradiated data without modifying original empirical constants derived from unirradiated data.
3.2 Thermal conductivity correlations for mixed oxide fuels
Due to complexity and inter-correlation of various parameters, especially burnup and material defects, thermal conductivity correlations for nuclear fuels are often derived from experiments by which certain mathematical functions are formulated with some empirical parameters to
68 match experimental data. Although there have been some efforts to evaluate thermal conductivity from computational atomistic simulation, these simulation results are only applicable to fresh fuel [60] [61] [62] [63] [64] [65] [68] [67] [70] [69] [70] [71]. They did not extend their evaluation beyond zero burnup potentially because of the complexity and uncertainty in representing atomic models of nuclear fuels at high burnup.
For oxide fuels, a number of thermal conductivity measurements have been carried out since the 1960s and since then a number of empirical correlations based on original measurements or literature review have been proposed [55] [72] [73] [74] [75] [76] [77]. In this section, the focus is given to experimental programs and thermal conductivity correlations developed for MOX fuel. Direct experimental programs for MOX thermal conductivity are somewhat underrepresented when compared with UO2 [56]. Before describing thermal conductivity correlations, it is worth to mention experimental results published over the years on thermal diffusivities for both fresh and irradiated samples. Experiments related to fresh samples will be covered first followed by irradiated ones. It is also beneficial to mention the difference between
MOX for light water reactors (LWR) and fast reactors (FBR). For LWR MOX, the PuO2 weight fraction is typically below 15% whereas it will be greater than 20% for FBR MOX. This is because of differences in fission cross sections in thermal and fast spectrum.
For unirradiated samples, Gibby evaluated thermal conductivity of fresh samples of UO2, PuO2 and MOX from 5 to 30 wt% PuO2 from laser-flash thermal diffusivity measurements [78]. The study observed a systematic reduction in thermal conductivity as a result of plutonium addition into the UO2 matrix. Using a simplified theory of lattice defect thermal resistance in dielectric solids, the model predicts that the thermal conductivity of MOX is lowest at 70 wt% PuO2.
Fukushima et al. [79] studied the effect of fission product addition (Nd and Eu) to MOX thermal conductivity. In this work, thermal diffusivities of near-stoichiometric MOX at 20 wt% of PuO2 with small additions of Nd and Eu up to 10 mol% were measured. It was found that the addition Nd and Eu gradually decreased MOX thermal conductivity. Thermal resistivity derived from lattice defect models agreed well with measurements. It was also found that the effect of lattice strain is more important than the effect of mass difference.
69
However, Schmidt et al. [80] [81] and Beauyy [82] [83] published contradicting results to that of
Gibby. They observed that the inclusion of Pu into UO2 matrix reduces overall thermal conductivity but not in a continuous way. In the range of 0 to 22 wt% PuO2, Beauyy noticed two local maximum at around 3 and 15 wt% PuO2 while Schmidt et al. found a local maximum at 15 wt% PuO2.
Bonnerot [84] conducted thermal diffusivity measurement of 18 values of PuO2 weight fraction from 0 to 100 wt%. It was found that the addition of PuO2 has moderate effect of MOX thermal conductivity. For stoichiometric MOX at 20 wt% of PuO2, the thermal conductivity was found to be 6% lower than stoichiometric UO2.
Duriez et al. [55] carried out a number of measurements on fresh MOX samples with Pu concentration from 3 to 15% and found that MOX thermal conductivity was significantly lower than that of UO2. It was also observed that the effect of PuO2 in the range of 3 to 15 %wt was negligible. Thermal diffusivity measurement of higher PuO2 concentration at 21.4 wt% was also conducted and the results were found to be significantly lower than samples in the range of 3-15 wt% PuO2. The study concluded that there was a significant difference in thermal conductivity for FBR and LWR MOX due to PuO2 concentration and microstructure caused by the addition of PuO2 into UO2 matrix.
Morimoto et al. [85] [86]measured thermal conductivity of (U,Pu,Am)O2 at 30 wt% of PuO2 and
0.7-3 wt% of AmO2 from 900 to 1773 K. The studied found that the presence of Am slightly decreased the thermal conductivity of the mixtures. The experimental results fitted well with a classical phonon transport model of 1/(A+BT) up to about 1500K. It was observed that the coefficient A increased linearly with Am content but small variation for the coefficient B.
Sengupta et al. [87] investigated the thermal conductivity of stoichiometric MOX at 44 wt%
PuO2—the highest PuO2 concentration in open literature. As expected, thermal conductivity of
MOX at 44 wt% PuO2 was found to be noticeably lower than MOX at 30 wt% PuO2. o Reportedly, at 1273 K (1000 C), thermal conductivity of MOX containing 44 wt% PuO2 was found to be 1.803 W/m/K while that of 30% PuO2 was 2.326 W/m/K. Similar to Gibby [78] and
70
Martin [88], this study confirmed the effect of PuO2 content on thermal conductivity reduction of MOX.
Morimoto [89] performed thermal diffusivity measurement of hypo-stoichiometric and stoichiometric MOX fuels at temperature from 990 to 2190 K using the laser flash method. This study intended to correlate the effect of high temperature to oxygen-to-metal ratio of the samples. They observed a reduction in O/M for stoichiometric samples above 1800 K while they found no significant changes in O/M for samples with O/M less than 1.95.
Staicu [90] investigated the effect of heterogeneity in fresh MOX fuel by comparing samples manufactured differently with different microstructures. For homogenous MOX fuel, samples came from SBR (Short Binderless Route) route with 4.8, 5.6 and 11.1 wt% PuO2 and a sol-gel
MOX with 7.8 wt% PuO2. The heterogeneous fuels were MIMAS (Micronized Master Blend)
MOX with 7.0 and 9.0 wt% PuO2. It was observed that the thermal conductivity of homogenous
MOX fuel was close to that of UO2 and the weight fraction of PuO2 does not have significant impacts on thermal conductivity. However, thermal conductivity of heterogeneous MOX fuel was found to be significantly lower than homogenous MOX and UO2. Given a similar PuO2 content of both homogenous and heterogeneous samples, it was pointed out that the main cause of thermal conductivity difference was due to the stoichiometry effects and this difference tends to disappear at high burnup.
Prieru et al. [91] measured various thermal properties including thermal conductivities of fresh samples of Np-MOX and Am-MOX which are mixed oxides of uranium and plutonium with small additional neptunium and americium. The uranium and plutonium weight fraction were 74% and 22%, respectively. The addition of Np was in the range of 0.5-2wt% and 0.35-2 wt% for Am. It was found that Np-MOX has a slightly higher thermal conductivity than Am-MOX.
However, both samples had lower thermal conductivity than UO2 but comparable, within 10% measurement uncertainty, to thermal conductivity of LWR MOX and FBR MOX from literature sources.
71
For irradiated samples, Yamamoto et al. [92] measured thermal conductivities of FBR MOX with 28.8 wt% at 0.5 MWd/kgHM and 17.7 wt% Pu at 8, 10 and 35 MWd/kgHM but surprisingly did not clearly observe thermal conductivity degradation due to burnup.
A series of experiments have been carried out at the Halden reactor to investigate the effect of burnup on thermal properties of both UO2 and MOX [93] [94] [95]. High burnup was achieved by retrieving fuel rods from commercial BWRs, refabricating them into instrumented fuel rods equipped with centerline thermocouples, pressure gauge, cladding elongation detector, and finally re-irradiating these test rods in the Halden boiling water reactor (HBWR). Reported burnup of these samples for MOX and UO2 test rods were 84 and 72 MWd/kgHM, respectively [96]. Fresh MOX fuel rods were also tested at the Halden reactor [97]. With on-line instrumentation capability, it was possible to continuously monitor the evolution of fuel centerline temperature and internal rod pressure since the beginning of irradiation which corresponds to fuel burnup from 0-32 MWd/kgHM. However, due to the nature of in-pile reactor experiments, fuel thermal conductivities were not directly measured. Instead, the researchers had to rely on various in-situ measurements such as centerline temperature, internal rod pressure, fuel dimension to formulate thermal conductivity correlations which introduced many uncertainties in both measurements and theoretical models. Alternatively, these experimental results were used to benchmark other correlations developed from direct out-of-pile thermal diffusivity measurements [96] [97] [98].
Cuzzo et al. [99] reported thermal diffusivity of irradiated UO2 and MOX samples with fuel burnup in the range of 31-36 MWd/kgHM. Thermal conductivity was inferred from laser-flash diffusivity measurement. It was found that, at similar burnup, the thermal conductivity of UO2 and MOX was in the same range. One thing to note about this finding was that the PuO2 content in MOX samples before irradiation was relatively small (3.7 wt%). So, the effect of plutonium depletion in MOX and plutonium buildup in UO2 may contribute to this behavior.
Staicu et al. [100] measured thermal diffusivity of LWR MOX at 23, 42, 44 and 47 MWd/kgHM from MOX samples fabricated from different methods i.e. SBR, MIMAS, and OCOM. SBR (Short Binderless Route) yields a homogenous distribution of plutonium in the matrix whereas MIMAS (Micronized Master Blend) and OCOM (Optimized Co-milling Method) provide
72 heterogeneous microstructure of the MOX samples. The results suggested that there thermal conductivities were not significantly different between homogeneous and heterogeneous MOX.
Nakae et al. reported another experiment at the Halden reactor on high burnup fuel rods up to 74
MW/kgHM [101]. Fuel samples were LWR MOX (8.4 wt% PuO2), and UO2 (8% enrichment).
Both heterogeneous (MIMAS) and homogenous (SBR) MOX and UO2 fuel rods were taken from commercial PWRs, re-fabricated and re-instrumented into instrumented fuel assemblies (IFA). Thermal conductivities were indirectly inferred from plots of centerline temperatures vs. linear heat generation rates (LHGR). In this case, higher thermal conductivity corresponded with to lower fuel temperature given the same amount of LHGR increase. It was unexpectedly found that the thermal conductivity of UO2 was less than that of MOX at 80 MWd/kgHM. Various speculations were given to explain this experimental finding such as non-stoichiometry, cracking, porosity and microstructure. However, no definitive causes of these counter-intuitive findings were confirmed and more direct thermal conductivity or diffusivity measurements were needed to confirm these results.
There are a sizable number of correlations for fresh MOX that have been proposed over the years [102] [55] [72] [75]; most of them were developed from out-of-pile thermal diffusivity measurements or literature reviews of published experimental data. In this thesis, the work of Philipponneau et al., Duriez et al., Inoue et al. and Baron et al. are presented.
Philipponneau [102]proposed a correlation based on a literature review of thermal diffusivity measurements of FBR MOX at approximately 20 wt% PuO2. Both stoichiometric and non- stoichiometric fuels were included in the formulation. In this work, the fuel composition was written in short form as U0.8Pu0.2U2-x where x represents a deviation from stoichiometry. In this case, when x is equal to zero, it means that the fuel is at stoichiometric composition and the subscripted numbers represent mole fraction of each constituent. For fresh MOX, the correlation takes into account the effect of temperature, non-stoichiometry, and porosity. However, Philipponneau suggested that the effect of plutonium content in the range of 15-30 wt% in MOX fuels should be ignored. The Philipponneau’s correlation is given by:
1 = + 76.38 × 10 (2) 1.528√ + 0.00931 − 0.1055 + 2.885 × 10
73
where k95 = thermal conductivity of MOX at 95% TD (W/m/K)
T = temperature in range of 500 K < T < 3000 K
x = deviation from stoichiometry = |2.0 – O/M|