Fully Ceramic Microencapsulated Fuels: Characteristics and Potential LWR Applications
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
In the last few decades, the nodal diffusion calculation accompanied with the assembly homogenization has been widely employed. The single assembly calculation with the reflective boundary condition (BC) would be carried out to determine the assembly-homogenized parameters. Because of the reflective BC, no axial and/or radial leakage rates have been taken into accounted for the lattice calculation. To represent the global influence of the reactor core on an assembly calculation, one mostly frequently uses a procedure called leakage model to approximately take into account the absorption due to the macroscopic structure of the physical flux in a heterogeneous situation. For the assembly homogenization, this absorption is introduced as the product of the leakage coefficient and the buckling obtained for an equivalent homogeneous medium using the well-known homogeneous B1 theory. Clearly, this is an imperfect representation of leakages in a heterogeneous medium. However, the errors introduced by the foregoing procedure are small for the whole? -assembly homogenization in a PWR. When it comes to the pin-by-pin calculation, it is not the case. The leakage rate? will be affected by the neutron streaming effect caused by the heterogeneity of the assembly. Especially when poison pins or empty channels are present, the streaming effect becomes important. The errors introduced appear particularly visible. The heterogeneous leakage model? which assumes that the neutron flux is variational and the leakage rate is not uniform over the lattice pitch is a more appropriate choice. When the heterogeneous leakage model is applied in the lattice calculation for the pin-cell homogenization, the pin-homogenized problem would be changed compared with the original assembly problem without any homogenization technique. So, the Cooperation of the heterogeneous leakage model with the pin-cell homogenization theory is unsettled. In this paper, both homogenous and heterogeneous leakage models are evaluated in pin-by-pin calculation. Their implements used in pin-cell homogenization of the lattice calculation are studied. The diffusion coefficient obtained from transport cross-section and the space-dependent leakage coefficient together with the two leakage models are analyzed for the practical reactor-core applications.
The Virtual Environment for Reactor Applications components included in this distribution include selected computational tools and supporting infrastructure that solve neutronics, thermal-hydraulics, fuel performance, and coupled neutronics-thermal hydraulics-fuel performance problems. The infrastructure components provide a simplified common user input capability and provide for the physics integration with data transfer and coupled-physics iterative solution algorithms. Neutronics analysis can be performed for 2D lattices, 2D core and 3D core problems for pressurized water reactor geometries that can be used to calculate criticality and fission rate distributions by pin for input fuel compositions. Two neutronics capabilities are provided in this distribution including MPACT (from the University of Michigan) that supports 2D lattices and Insilico (from Oak Ridge National Laboratory) that supports 2D lattices, 2D core planes, and 3D core geometries. Both codes include integrated cross section capabilities that provide problem-specific cross sections for the problems defined. A thermal-hydraulics capability is provided with CTF (an updated version of the COBRA-TF code developed by the Pennsylvania State University) that allows thermal-hydraulics analyses for single and multiple assemblies using the simplified VERA common input. This distribution also includes an initial coupled neutronics/thermal-hydraulics capability to allow calculations of single assemblies using Insilico or MPACT coupled with CTF. The VERA fuel rod performance subcomponent BISON-CASL calculates, on a 3D basis, fuel rod temperature, fuel rod internal pressure, free gas volume, clad integrity and fuel rod waterside diameter. These capabilities allow simulation of power cycling, fuel conditioning and deconditioning, high burnup performance, power uprate scoping studies, and accident performance. Specifically targeted at VERA’s core simulator functionality (VERA-CS), the existing VERA subchannel-neutronics capability based on CTF and Insilico has been extended to use the axisymmetric 2D R-Z fuel rod modeling functionality of BISON-CASL to produce a three-way coupled neutronics, thermal-hydraulics and fuel performance capability known as Tiamat. This distribution is intended to be used for testing and evaluation of the released components of VERA. Testing within CASL has focused specifically on Westinghouse four-loop reactor geometries and conditions with example problems included in the distribution.
The Rossendorf Research Reactor (RRR), or Rossendorfer Ringzonenreacktor, was an Argonaut type zero power research reactor located at the Rossendorf research facility in Dresden Germany. The RRR consists of an annular core surrounded by a gr¬aphite reflector. Inside the core is an additional removable internal graphite reflector. This internal reflector was replaced in December of 1972 with a rapid deployment lattice known as Schnelles Einsatz-Gitter (SEG). The insertion of the SEG experimental lattice leads to a fast-thermal coupled system with the annular RRR reactor being the thermal driver and the SEG being the fast insertion lattice. The RRR/SEG integral experimental setup was used to perform small sample reactivity measurements using the pile oscillator method. By loading the SEG lattice with specially designed unit cells, a specific adjoint spectra could be obtained (only SEG 6 had a special radial arrangement). With a nearly energy-independent adjoint spectrum, the slowing-down effect becomes a constant and essentially disappears, such that any chan¬¬¬¬ge in reactivity can be attributed to the absorption component of the sample material. SEG 4, 5, and 7 exhibit this characteristic. If the adjoint function is strongly energy-dependent, then the scattering component will dominate. In the case of both SEG 6 lattices, where the adjoint is also monotonously rising, the scattering effect will always be negative. This analysis aims to re-evaluate the central reactivity worths of different samples for the SEG 4, 5, 6, and 7 experiments using modern continuous energy monte-carlo.
Evaluation of thermal properties is a crucial factor for nuclear fuel performance. During reactor operation, the accumulation of fission products and irradiation-induced lattice defects are responsible for degradation in thermal conductivity. Consequently, it affects fuel temperature and fission gas release (FGR) among other Multiphysics processes important for economics and safety analysis. We analyze the implications of point defects (PD) accumulation described using a rate theory (RT) Model on lattice thermal conductivity of UO 2 . Here, we demonstrate that fission rate-dependent point defect concentrations have the largest impact on in-pile thermal conductivity in the periphery of light water reactor fuels below a temperature threshold governed by the migration barrier of defects. Our analysis provides a mechanistic description of this phenomena which current fuel performance codes treat empirically. The reduction of thermal conductivity in the low -temperature rim region acts as additional thermal resistance and leads to a temperature notably larger than suggested by Lucuta thermal conductivity correlation. These effects are anticipated to have notable impacts when fuels are exposed to accelerated radiation. The impact of such point defect-informed treatment of thermal conductivity on fuel performance is evaluated by a detailed analysis of fission gas behavior and its release. We consider several models capturing different stages of fission gas bubble evolution and fission gas release (FGR). Finally, a new fission rate-dependent correction to the Lucuta correlation is proposed. The results show a significant reduction in thermal conductivity at the fuels’ periphery and an increase in fuel centerline temperature specifically at low burnups. Ultimately a modified LC shows a higher FGR compared to the original LC, while the acceleration process results in a reduction in overall FGR.
The fully ceramic microencapsulated (FCM) fuel is based on tri-structural isotropic (TRISO) coated particles embedded in silicon carbide (SiC) matrix, which is among the accident tolerant fuel concept options. With the strength of high fission product retention, the FCM fuel provides potential superior safety performance. Compared to the conventional standard UO{sub 2} fuel, the deficiency of FCM fuel is its considerable lower fissile inventory, which results in a reduced reactor power or/and life. Supposing the FCM fuel is used as direct replacement of the fuel in commercial PWRs, the spectrum is softened due to the usage of SiC as matrix, and the moderator temperature coefficient is possible to be positive at the beginning of life, leading to high risk during operation. A neutronic feasibility study of FCM fuel utilized in commercial PWRs is performed in the paper. The fuel enrichment and under-moderated lattice are evaluated by analyzing the typical TRISO particle with UO{sub 2} kernels which is widely used in high temperature gas reactors and that with enlarged UN fuel kernels. It was found that under the mandatory constraint of low enriched uranium, the FCM fuel with UO{sub 2} kernels, although relatively mature, is only capable of small reactor designs, and the assembly lattice needs to be modified. The FCM fuel with UN kernels, when directly loading in the commercial PWRs, can yield comparable reactor power and life; however, a higher fuel enrichment and the fuel design and fabrication need to be optimized. (authors)
Zirconium nitride holds promise for improving nuclear reactor components, where thermal conductivity is important, highlighting the need for a detailed understanding of its phonons. However, the phonon dispersion derived from previous inelastic neutron scattering data is in disagreement with modern theoretical calculations, pointing to the necessity for further experimental validation. Here, in this study, we use inelastic neutron scattering experiments on polycrystalline ZrN at 6, 100, 300, 400, and 600 K, and density functional theory (DFT) calculations to study the momentum-resolved scattering and the phonon density of states. Both experimental and theoretical results reveal a phonon cutoff energy near 65 meV, in stark disagreement with previous single-crystal neutron scattering measurements indicating a cutoff near 75 meV. The 65 meV cutoff energy is, however, in reasonable agreement with previous measurements on ZrN 0.9 . Consistency between our theoretical and experimental results confirms that ZrN is a wide phonon band gap system and that DFT provides a reliable description of its phonons.
Thermal conductivity of neutron irradiated copper at low temperature noting lattice defects
The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.
The deterministic full core simulators require homogenized group constants covering the operating and transient conditions over the entire lifetime. Traditionally, the homogenized group constants are generated using lattice physics code over an assembly or block in the case of prismatic high temperature reactors (HTR). For the case of strong absorbers that causes strong local depressions on the flux profile require special techniques during homogenization over a large volume. Fuel blocks with burnable poisons or control rod blocks are example of such cases. Over past several decades, there have been a tremendous number of studies performed for improving the accuracy of full-core calculations through the homogenization procedure. However, those studies were mostly performed for light water reactor (LWR) analyses, thus, may not be directly applicable to advanced thermal reactors such as HTRs. This report presents the application of SuPer-Homogenization correction method to a hypothetical HTR core.
The MARVEL (Microreactor Applications Research Validation and Evaluation) microreactor utilizes natural circulation as core cooling mechanism and liquid metal as a primary coolant. Moreover, the reactor core has a pitch-to-diameter ratio of 1.056, which is considered a tight lattice configuration. Numerous studies have widely reported that Reynolds-Averaged Navier-Stokes (RANS) turbulence models inaccurately predict heat transfer in liquid metals and fail to capture flow pulsations that can occur within tight lattices, leading to further inaccuracies in simulation results. Therefore, evaluating the accuracy of RANS turbulence models in the thermal-hydraulic analysis of the MARVEL microreactor core is crucial for assessing reactor safety. In this study, a Large Eddy Simulation (LES) of the MARVEL microreactor core subchannel was conducted and compared with a RANS simulation to evaluate the accuracies and conservatism of the RANS model. The flow pulsation in a tight lattice predicted by LES enhanced the heat transfer, whereas the RANS model underpredicted it. Consequently, the RANS model predicted the peak cladding temperature higher than the LES model, but the discrepancy between the two approaches was not significant due to the good thermal characteristics of the liquid metal. It can be concluded that the steady-state RANS model is effective for the thermal analysis of liquid-metal-cooled MARVEL microreactor core and can provide conservative predictions from a safety analysis perspective.
The Virtual Environment for Reactor Applications components included in this distribution include selected computational tools and supporting infrastructure that solve neutronics, thermal-hydraulics, fuel performance, and coupled neutronics-thermal hydraulics problems. The infrastructure components provide a simplified common user input capability and provide for the physics integration with data transfer and coupled-physics iterative solution algorithms. Neutronics analysis can be performed for 2D lattices, 2D core and 3D core problems for pressurized water reactor geometries that can be used to calculate criticality and fission rate distributions by pin for input fuel compositions. MPACT uses the Method of Characteristics transport approach for 2D problems. For 3D problems, MPACT uses the 2D/1D method which uses 2D MOC in a radial plane and diffusion or SPn in the axial direction. MPACT includes integrated cross section capabilities that provide problem-specific cross sections generated using the subgroup methodology. The code can be executed both 2D and 3D problems in parallel to reduce overall run time. A thermal-hydraulics capability is provided with CTF (an updated version of COBRA-TF) that allows thermal-hydraulics analyses for single and multiple assemblies using the simplified VERA common input. This distribution also includes coupled neutronics/thermal-hydraulics capabilities to allow calculations using MPACT coupled with CTF. The VERA fuel rod performance component BISON calculates, on a 2D or 3D basis, fuel rod temperature, fuel rod internal pressure, free gas volume, clad integrity and fuel rod waterside diameter. These capabilities allow simulation of power cycling, fuel conditioning and deconditioning, high burnup performance, power uprate scoping studies, and accident performance. Input/Output capabilities include the VERA Common Input (VERAIn) script which converts the ASCII common input file to the intermediate XML used to drive all of the physics codes in the VERA Core Simulator (VERA-CS). VERA component codes either input the VERA XML format directly, or provide a preprocessor which can convert the XML into native input. VERAView is an interactive graphical interface for the visualization and engineering analyses of output data from VERA. The python-based software is easy to install and intuitive to use, and provides instantaneous 2D and 3D images, 1D plots, and alpha-numeric data from VERA multi-physics simulations. Testing within CASL has focused primarily on Westinghouse four-loop reactor geometries and conditions with example problems included in the distribution.
Solid Fueled Thorium Molten Salt Reactor (TMSR-SF1) is Chinese Fluoride-salt-cooled High Temperature Reactor (FHR) designed by the TMSR group. The construction of this test reactor is expected in the near years. Using TRISO particle based pebble fuel leads to various challenges in neutronics modeling, such as double heterogeneity effect, preserving exact number of particles, their random distributions in pebble, and etc. For the purpose of verifying the calculations, benchmark study of infinite lattice of pebble fuel within FLiBe coolant for the TMSR-SF1 design has been performed. Back-to-back comparisons have been made between TMSR and MIT researchers using different Monte-Carlo based codes. Very good agreements have been found and the results quantitatively evaluate the neutronic effect of various modeling approaches. In addition, sensitivity study has been carried out for the reference pebble-bed fuel lattice. (authors)
Accurate, high-resolution full core simulations of prismatic high temperature reactors for burnup and transient analysis remain a challenge. The cores are geometrically large, but some features are quite small (tri-structural isotropic (TRISO) particles, burnable poison pellets). The treatment of both neutron scattering in graphite and resonance capture are complex and not adequately captured using the methods traditionally used in High Temperature Gas-cooled Reactors (HTGR) and commonly used in Light Water Reactors (LWR). For transient analysis, temperature feedback is a function of the TRISO fuel form, but full core models cannot resolve phenomena at this scale without averaging over space and energy. For burnup calculations, a second level of heterogeneity must be resolved to capture the local effects of fuel and burnable poisons, while accurately propagating their effects through and between blocks. The HTGR Methods and Simulation group has been utilizing the INL-developed PHISICS/RELAP5-3D{sup C} (P/R), code suite for participation in the Organization for Economic Cooperation and Development (OECD) Modular High Temperature Gas Cooled Reactor (MHTGR) Benchmark, funded by the US Department of Energy (DOE) under the Advanced Reactor Technologies (ART) program. The MHTGR-350 benchmark consists of several lattice, steady-state and transient exercises, and the P/R code has been successfully applied to most of the cases defined. In addition to this code-to-code verification activity, it is planned that P/R will be validated using data from several Loss of Forced Cooling (LOFC) tests performed at the High Temperature Test Reactor (HTTR) in a joint effort with Japan Atomic Energy Agency (JAEA) under the Civil Nuclear Energy Working Group (CNWG) project. These verification and validation simulations will identify code and model discrepancies, sensitivities, and establish the current state-of-the-art in HTGR simulation. The Department of Astronautical, Electrical and Energy Engineering (DIAEE), 'Sapienza' University has been involved in the development of P/R since 2014 through various graduate and PhD student projects. In the framework of this collaboration the P/R optimization for application on the OECD MHTGR-350 and HTTR models will be carried out.
During reactor operation, nuclear fuels are subject to extreme temperature and irradiation conditions which can significantly degrade the fuel's thermal transport properties. The reduction in thermal conductivity of the fuel as a result of irradiation-induced lattice defects is arguably the most important fuel performance metric in regard to reactor efficiency and safety. Because thorium dioxide (ThO2) is suitable as a model system for more complex materials such as UO2 and its mixed oxides, we present a theoretical investigation of thermal conductivity reduction seen in defect-bearing thorium dioxide and compare directly to experimental measurements. Phonon-mediated thermal transport of the fuel is modeled by a solution to the Boltzmann transport equation (BTE) for phonons. A cluster dynamics (CD) model for lattice defect evolution during irradiation predicts defect densities which are used as input to the BTE for modeling phonon-defect scatterings. Phonon scatterings by lattice defects include those from point defects and vacancy clusters and interstitial clusters of various sizes. The CD model is benchmarked against structural defect characterization of irradiated thorium dioxide using electron microscopy. Thermal conductivity predicted by the BTE model is compared to measured values for irradiated thorium dioxide specimens below room temperature to isolate effects of phonon-defect scattering from intrinsic 3-phonon processes, which dominate at higher temperatures. The computed conductivity values are in partial agreement at temperatures close to room temperature while slight deviations are observed at the lowest measured temperatures, suggesting that implemented phonon-defect scattering cross-section expressions may not be adequate for low temperatures. The presented work provides a necessary investigation of the influence of irradiation induced defects on fuel performance and represents a first step toward a full characterization of phonon mediated thermal transport in irradiated materials with complex defect microstructure.
The Method of Characteristics (MOC) has seen wide interest in reactor physics because of its accuracy and efficiency in computing lattice physics problems. While most of its use has been in solving 2D problems, there has been recent interest in extending MOC to 3D in order to more accurately calculate 3D power distributions in LWRs. While the method is naturally extensible to 3D, it presents significant computational difficulties. OpenMOC, an open source MOC simulator, has been extended from 2D to 3D capabilities. Simulation results are presented in this study for the full core 3D simulation of the BEAVRS benchmark, adequately resolved in space and angle.
Many liquid metal plasma-facing component concepts for fusion reactors require a capillary porous material to wick the liquid metal to the surface of the component. In this work, samples of 3D-printed molybdenum lattices designed to have pore sizes ranging from 100 to 1000 µm were fabricated by two different suppliers. The dimensions of the lattice structures were characterized using an optical microscope, revealing that the lattice walls were almost systematically thicker than they were designed to be and in turn the pores are all smaller than they were designed to be. In many cases the pores were completely clogged. Several spectroscopic techniques were used to characterize the chemical composition of the materials, collectively indicating the presence of a thick layer of oxygen-rich material extending more than 1 µm below the surface. Compression testing revealed significant asymmetry in the mechanical behavior depending on which direction the force was applied relative to the print direction. Lastly, lithium wetting tests of several lattice samples were performed, and a sample with ~250-µm pores was found to soak up lithium quickly after initial wetting of the surface occurred. The implications of these results for future liquid metal plasma-facing component design are discussed.
The purpose of this paper is to present an analysis of the hafnium configuration of FLUOLE-2 program. FLUOLE-2 is a set of benchmark-type experiments dedicated to neutron attenuation analysis with the aim of improving the TRIPOLI-4{sup R} Monte Carlo code validation and industrial neutron physics scheme. The CEA developed this program, with the support of EDF, in order to be representative of 900 and 1450 MWe Pressurized Water Reactors operated in France. For that purpose, different stainless steel structures have been designed and appropriately positioned inside the EOLE facility located at Cadarache CEA center. EOLE is a pool type zero power reactor, composed of a cylindrical aluminum vessel with an over-structure of stainless steel, that is able to contain various types of core and related structures. The FLUOLE-2 core has been designed as a 29 x 29 pins square lattice of fuel rods. Different kind of dosimeters (cobalt, gold, tin, rhodium, indium, iron, nickel, titanium, aluminum, and vanadium) were irradiated inside and outside the core. (authors)