Search NASA⌕ Search

SEARCH · Search NASA

Results for “fuel cycle length”

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.

225 records · Page 13

Calibration of the Diffusivity Predictions of Centipede Using Approximate Bayesian Computation and Applications in Nyx (Engineering Scale) and Xolotl-MARMOT (Meso-Scale) Simulations

Fission gas evolution and release in UO 2 nuclear fuel are important fuel performance metrics and occur in several distinct stages: 1) nucleation, growth and resolution of intra-granular bubbles, 2) diffusion to grain boundaries and 3) nucleation and growth of bubbles at grain boundaries, which eventually form a connected network (percolation) enabling release of gas from grain boundaries through connections to triple junctions, grain edges or free surfaces. The NE-SciDAC project is developing several computational tools to model this problem, which are connected in a hierarchical multi-scale framework. The information transfer in the multi-scale framework is a critical step that, in addition to best-estimates, should include uncertainty quantification. Despite taking a first-principles multi-scale approach, there is a need to perform parameter calibration to ensure consistency with available experimental data. In the present study, uncertainty quantification (UQ) and parameter calibration is demonstrated for one of the lower length scale codes in the multi-scale framework (Centipede) and then the results, including instances of the propagated uncertainties, are used in other codes within the framework, specifically Nyx and Xolotl-MARMOT. We calibrated the model parameters in Centipede, a computer code used to predict diffusivities of uranium (U) and xenon (Xe) in the context of the simulation of fission gas in uranium oxide (UO 2 ) nuclear fuel. The Centipede code depends on 183 parameters, all of which are subject to uncertainty. The three data sets used in our calibration effort are taken from the literature. This data is available as a set of measurements, including measurement errors. Our goal is to calibrate a statistical model that predicts both the value of the measurement and the uncertainty associated with the measurement. We perform a Bayesian calibration of the model parameters using a dedicated approximate Bayesian computation (ABC) likelihood function. To avoid excessive computational costs, we replace the expensive Centipede simulation code by a higher-order surrogate model, constructed using only the 9 most important parameters. These important parameters are identified by a preliminary global sensitivity analysis (GSA) study. Among the important parameters are T0 (the temperature at which UO 2 is perfectly stoichiometric) and Hf_pO2 (the temperature dependence of the oxygen (O) partial pressure) that should be considered as operating conditions to be estimated along with the other parameters. We consider two different cases: one where we define one set of these operating conditions for all data sets, and one where we define distinct operating condition parameters for each data set. The Xe diffusivities predicted by the latter case show distinct features that could not be observed in the former. Next, we use the diffusivity predictions by Centipede as input to Nyx, a reduced order fuel performance code focused on gas behavior alone, in order to estimate quantities associated with inter-granular bubble formation at conditions specified by the experiments. Finally, the diffusivities obtained from the calibrated Centipede runs were used in coupled Xolotl-MARMOT simulations of intra- and inter-granular gas evolution. The results are compared to simulations using the baseline diffusivities from Turnbull et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Scoping Study for Fast Flux Testing in the Advanced Test Reactor

The value of fast spectrum reactors remains prominent in the nuclear technology portfolio. The performance of these reactors can be maximized with advancements in nuclear fuel technologies, but development of these technologies is currently held back by lack of fast spectrum test reactors available to the United States. Spectral modification of experiment positions in the thermal spectrum Advanced Test Reactor (ATR) has long been used to support fast reactor fuel development, but these methods have not been progressed to their full potential. This study investigated the use of concentric rings of aluminum-clad fuel plates in ATR flux traps and thermal neutron absorbing filters to increase fast neutron flux on test specimens. This concept was termed the Boosted Energy Advanced Spectrum Test (BEAST). This approach will enable irradiation of advanced fuel designs in prototypic-length fuel pins and representative flux environment to support post irradiation exams, enable transient testing, and produce the type of data that will permit lead test assembly irradiations in true Sodium Fast Reactors (SFRs) when they become available. Neutronic predictions were performed to investigate BEAST design options and thermal hydraulic models were produced to ensure feasibility of BEAST. Two versions were considered based on the geometric limitations of ATR’s small and large flux traps. The small version was found to be preferable due to slightly higher fast flux and fast-to-thermal neutron ratio. Perhaps more influentially, the small flux trap option was also preferred to avoid conflict with ongoing very high temperature reactor fuel irradiation programs in ATR’s large northeast flux trap. The small flux trap option provided less than half the test volume of the large version, but still had adequate volume for seven SFR pins in cross section which could be stacked two-high in ATR’s 1.2m long core to accommodate up to 14 EBR-II size pins. The preference for the small flux trap configuration should be revisited if additional collaborative test programs emerge with the need to irradiate a significant volume of additional specimens. Calculations were performed regarding a lithium deuteride ring to convert thermal neutrons into 14 MeV fusion neutrons. At the time this report was written these calculations were partially complete and it remains to be seen whether the concept would be worth including in BEAST. Given the preference for the small flux trap option, which does not afford enough room for the 14 MeV ring, it was concluded to defer future work on the lithium deuteride ring. This decision could be revisited if fusion material research programs emerge for collaborative testing in BEAST. A cadmium-lined specimen holder design was found to be adequate in filtering thermal neutrons and preferred over other neutron absorbers based on past experience with cadmium baskets. It was acknowledged that cadmium-bearing hardware would become depleted and need to be replaced occasionally, which appeared feasible from a mechanical design perspective. Neutronic studies investigated different enrichment levels in the booster fuel using uranium-molybdenum alloy dispersion fuel which has performed well in past ATR irradiations. Both options were able to drive fuel pins to SFR-like fission heating rates. The high enriched booster fuel option outperformed the low enriched option by ~20% on key metrics including fast flux and fast-to-thermal ratio, but the low enriched option was favored in order to broaden options for potential fuel suppliers. The preferred BEAST design options including cadmium filter with low enriched booster fuel in the small flux trap configuration was predicted to achieve 6.2E14 n/cm2sec fast flux (>0.1 MeV) with a fast-to-thermal ratio of 44.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Codisposal Waste Package Loading Options for DOE SNF and HLW PPT

The U.S. Department of Energy (DOE) is responsible for managing spent nuclear fuel (SNF) that is currently in, or will in the future come into, its possession. DOE must continue to safely store that SNF, transport it to an interim storage site or a repository, and dispose of it. These fuels come from a wide range of reactor types with various cladding materials and enrichments. Many of these reactors, now decommissioned, had unique design features, such as core configuration, fuel element and assembly geometry, moderator and coolant materials, operational characteristics, and neutron spatial and spectral properties, resulting in a large diversity of reactor and fuel designs. Because of the wide variety and conditions of SNF, a robust DOE Standard Canister was proposed that would confine radionuclides and preclude moderator. The DOE Standard Canister had four variations: 3.05-meter (10-foot) or 4.57-meter (15-foot) length, and 45.7-cm (18-inch) or 61.0-cm (24-inch) diameter. For ultimate disposal in the Yucca Mountain Repository, these canisters were to be grouped with 61.0-cm (24-inch) diameter high level waste (HLW) canisters in a 2.13-meter (84-inch) diameter co-disposal waste package. The smaller 45.7-cm (18-inch) diameter DOE Standard Canister could be placed in the middle of five HLW canisters. The larger 61.0-cm (24-inch) diameter DOE Standard Canister would take the place of one of the five HLW canisters on the outer ring in the co-disposal waste package. No DOE Standard Canisters have been loaded. A preliminary evaluation has estimated the number of elements of a fuel type that can fit into the different sizes of the DOE Standard Canister, but no definitive loading configuration has been selected. Changing the loading configuration could impact the number of loadable DOE Standard Canisters and the number of co-disposal waste packages needed for eventual disposition. This paper conveys the ranges of DOE Standard Canisters and HLW canisters that may be produced under certain conditions. It also examines the differences in the estimated number of co-disposal waste packages produced for eventual disposal when using different loading strategies in the DOE Standard Canister for Advanced Test Reactor (ATR), Peach Bottom, and High Flux Isotope Reactor (HFIR) SNF. Changing the loading configurations of ATR, Peach Bottom, and HFIR SNF slightly impacted the number of co-disposal waste packages that may be needed for ultimate disposal. The change in loading configuration was more impactful when a different canister was used, as opposed to varying the number of elements that could fit inside the same size canister. In one case, a reduction of co-disposal waste packages could be achieved by allowing mixing of short HLW canisters with long DOE Standard Canisters. The main conclusion from this analysis is that the ratio between HLW canisters and DOE Standard Canisters will drive the total number of co-disposal waste packages. If too many HLW canisters (i.e., more than five times the number of 18-inch DOE standard canisters) or DOE Standard Canisters are produced, some co-disposal waste packages may not have all positions filled. A co-disposal waste package may be filled with all HLW with no DOE Standard Canister, or a co-disposal waste package could be filled with a single DOE Standard Canister. A ratio that does not closely align to optimum could allow for the design of waste packages that hold just HLW canisters or just DOE SNF canisters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Develop Accurate Techniques for Passive SiC Temperature Monitoring of Miniature Samples for Cross-Cutting Applications

Passive thermometry is critically important because most fuels and materials irradiation experiments are not instrumented, and it is necessary to understand the irradiation temperature to properly interpret any post-irradiation examination data, including evolving properties and/or microstructures. The standard passive thermometry approach uses continuous dilatometry to evaluate changes in the instantaneous coefficient of thermal expansion during post-irradiation thermal annealing. This approach has limitations in terms of sample size (minimum length requirements) and the maximum irradiation temperature that can be accurately determined, which is limited by the reduced swelling (and therefore recovery) following higher temperature irradiation and limitations on the furnaces used with push-rod dilatometers. This work evaluates two new proposed techniques for post-irradiation evaluation of passive SiC temperature monitors: differential scanning calorimetry (DSC) and Raman spectroscopy. DSC is an extremely sensitive technique that can be used for any specimen geometry and is capable of higher temperature operation. Raman spectroscopy is similar in that it is a surface technique capable of examining extremely small samples (submillimeter), can be used with a heated stage up to 1,500°C (planned for future work), and is capable of mapping local irradiation temperatures throughout a sample. Existing SiC samples that were previously irradiated over a wide range of temperatures were cut into multiple pieces to allow for annealing studies using multiple different techniques: dilatometry, DSC, and Raman spectroscopy. This approach mitigates the concern that samples analyzed using one technique may have a slightly different irradiation history than those analyzed using a different technique. Recovery was clearly observed during annealing using both DSC and dilatometry. In some cases, a direct comparison could not be made due to some of the DSC runs accidentally including material from multiple specimens and issues with using an alternative DSC sample holder for the highest temperature annealing studies. Nevertheless, one trend was clear: the DSC runs resulted in higher irradiation temperatures compared to those of the dilatometry runs. Part of this could be attributed to the higher temperature ramp rates used during the DSC runs, which are often preferred to reduce noise in the measurements. By comparison, dilatometry has previously been shown to produce better data at lower ramp rates. Future work should further investigate the ideal ramp rate for both techniques to produce consistent results. Additional work should evaluate the best holder material to use for DSC runs exceeding 1,000°C to provide reliable data while preventing interactions between SiC and the holder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Develop Accurate Techniques for Passive SiC Temperature Monitoring of Miniature Samples for Cross-Cutting Applications

Passive thermometry is critically important because most fuels and materials irradiation experiments are not instrumented, and it is necessary to understand the irradiation temperature to properly interpret any post-irradiation examination data, including evolving properties and/or microstructures. The standard passive thermometry approach uses continuous dilatometry to evaluate changes in the instantaneous coefficient of thermal expansion during post-irradiation thermal annealing. This approach has limitations in terms of sample size (minimum length requirements) and the maximum irradiation temperature that can be accurately determined, which is limited by the reduced swelling (and therefore recovery) following higher temperature irradiation and limitations on the furnaces used with push-rod dilatometers. This work evaluates two new proposed techniques for post-irradiation evaluation of passive SiC temperature monitors: differential scanning calorimetry (DSC) and Raman spectroscopy. DSC is an extremely sensitive technique that can be used for any specimen geometry and is capable of higher temperature operation. Raman spectroscopy is similar in that it is a surface technique capable of examining extremely small samples (submillimeter), can be used with a heated stage up to 1,500°C (planned for future work), and is capable of mapping local irradiation temperatures throughout a sample. Existing SiC samples that were previously irradiated over a wide range of temperatures were cut into multiple pieces to allow for annealing studies using multiple different techniques: dilatometry, DSC, and Raman spectroscopy. This approach mitigates the concern that samples analyzed using one technique may have a slightly different irradiation history than those analyzed using a different technique. Recovery was clearly observed during annealing using both DSC and dilatometry. In some cases, a direct comparison could not be made due to some of the DSC runs accidentally including material from multiple specimens and issues with using an alternative DSC sample holder for the highest temperature annealing studies. Nevertheless, one trend was clear: the DSC runs resulted in higher irradiation temperatures compared to those of the dilatometry runs. Part of this could be attributed to the higher temperature ramp rates used during the DSC runs, which are often preferred to reduce noise in the measurements. By comparison, dilatometry has previously been shown to produce better data at lower ramp rates. Future work should further investigate the ideal ramp rate for both techniques to produce consistent results. Additional work should evaluate the best holder material to use for DSC runs exceeding 1,000°C to provide reliable data while preventing interactions between SiC and the holder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Conceptual Designs for Irradiation Creep Testing of SiC in HFIR

Understanding irradiation creep of nuclear fuel cladding is important to properly size the initial fuel-cladding gap and understand when pellet-cladding contact is expected to occur due to a combination of fuel swelling and cladding creep-down. Irradiation creep also plays a role in relaxing stresses that develop in-pile. Silicon carbide fiber–reinforced silicon carbide matrix (SiC/SiC) composites are the leading long-term accident-tolerant fuel cladding concept for light-water reactors (LWRs). Although some limited data are available regarding irradiation creep of the individual constituents (fibers, matrix), data regarding irradiation creep of SiC/SiC composites are currently insufficient. Additional data regarding irradiation creep compliance and the rupture lifetime (combination of creep and slow crack growth) are needed to understand material limitations. This work describes the design and development of two irradiation vehicles that are being pursued for testing SiC/SiC concepts in the High Flux Isotope Reactor (HFIR). The first is a passive experiment, referred to as the PRECISE experiment, that leverages the constant coolant pressure of HFIR to compress a metallic bellows and provide a well-characterized load to drive creep in a SiC/SiC dog bone specimen. The total creep strain would be quantified post-irradiation by measuring dimensional changes of the specimen length as well as local dimensional changes within the gauge region. Non-stressed specimens would also be irradiated under the same conditions to provide an indication of dimensional changes due to radiation-induced swelling in the absence of creep. A second, more complex experiment, referred to as the INSITE experiment, is being designed in parallel that would use pneumatics to pressurize a metal bellows and linear variable differential transformers (LVDTs) to measure the specimen displacement in situ during irradiation. Such an experiment would provide significantly more data regarding the evolution of the creep compliance as a function of dose and applied stress within a single experiment but would require significantly more development time and cost to execute. The primary concern with the INSITE experiment is the accuracy, reliability, and expected lifetime of the LVDTs during irradiation at elevated temperatures. Efforts are being made to adjust the experiment design and operating procedure to limit LVDT temperatures and mitigate or otherwise compensate for uncertainties due to factors such as temperature fluctuations, creep in the surrounding structural materials, and drift of the LVDTs. This work describes the experiment designs, thermal and structural analysis that were performed to ensure that the desired temperature and stress conditions can be achieved, some initial sensitivity analyses to predict the evolution of the radiation-induced specimen displacements, and potential sources of uncertainty in the measurements. Out-of-pile testing is being performed in parallel to confirm that the test trains achieve the expected stress states in the specimens and do not result in prohibitive stress concentrators (e.g., in the grip regions) that might risk pre-mature failure. The PRECISE experiments are proceeding toward fabrication and assembly with HFIR insertion planned during fiscal year 2026. The INSITE experiment is progressing toward out-of-pile demonstrations, which will provide more conclusive evidence regarding the feasibility of executing these tests in HFIR or whether alternative displacement monitoring techniques may need to be considered.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Organic Iodide Sorption from Dilute Gas Streams

Reprocessing used nuclear fuel releases volatile radionuclides, including 129 iodine (I), into the off-gas of a processing plant. Volatile radioiodine could be present in several forms, depending on the chemistry of the process used and the off-gas stream. Inorganic I 2 is expected to be the predominant I species in the dissolver off-gas (DOG), with minor organic iodides present. The bulk of the I is expected to volatilize into the DOG in parts-per-million-level (ppm) concentrations. In contrast, in the vessel off-gas (VOG), most of the volatile I is expected to be found as organic iodides, such as CH 3 I, C 4 H 9 I, and C 12 H 25 I. These species are expected to be present in parts-per-billion-level (ppb) concentrations but require abatement, even at low expected concentrations, to meet regulatory emissions limits in the United States. Historically, studies of I abatement by Ag-functionalized sorbents have focused on inorganic I in the DOG, but in the last few years, more research attention has been given to organic iodides, especially longer chain species, such as C 4 H 9 I, and C 12 H 25 I. This report has three main goals: (1) to present new data generated at Oak Ridge National Laboratory (ORNL) in FY21 on the sorption behavior of organic iodides on AgZ, (2) to summarize and synthesize organic iodide data produced by ORNL and Idaho National Laboratory (INL) over the last 4 years to answer questions on organic iodides behavior outlined in the 2018 joint test plan (Jubin et al. 2018), and (3) to propose a VOG abatement system design that can provide the capture efficiencies required to meet I emission limits. The 2021 ORNL experimental campaign tested the effects of organic iodide speciation and concentration in the off-gas, superficial velocity of the off-gas, and effects of aging on AgZ sorbent capacity. These studies found that the sorption rate of organic iodides by AgZ depends on the hydrocarbon chain length and the concentration in the off-gas. Higher molecular weight organic iodides adsorb to AgZ more slowly than I. At a concentration of 50 ppm concentration in the off-gas, CH 3 I loads 8% slower, C 4 H 9 I loads 20% slower, and C 12 H 25 I loads 40% slower than I. The lowest concentration loading rates calculated were in 5 ppm organic iodide gas streams in which AgZ gained on average 0.14 mg I/g sorbent/hour in the bench scale test system. Thus, longer sorbent beds might be needed to accommodate slower loading rates onto AgZ in lower concentration gas streams. Although sorption rate varies as a function of hydrocarbon chain length, the saturation concentration of the sorbent for these I-bearing species does not vary. Aging AgZ in a humid air stream for 9 months drops the overall sorbent capacity by ~35% for CH 3 I, ~50% for C 4 H 9 I, and ~40% for C 12 H 25 I. This results in a saturation capacity between 35 and 70 mg I/g sorbent for the aged AgZ. In conjunction with recent data produced by INL, these data are used to estimate the mass transfer zone (MTZ) and decontamination factor (DF) for sorbent beds of AgZ. Sorption tests performed with iodide gas concentrations of about 1 ppm and higher at a superficial gas velocity of 10 m/min, indicate that MTZ depths for these conditions tend to range between about 8-20 cm. Tests performed at lower concentrations between 50-90 ppb and at gas superficial velocities of 1, 10, and 20 m/min indicate that the MTZ depth increases with increasing superficial gas velocity. The 20 m/min test indicates that the MTZ for those conditions was at least 22 cm, and doubling the superficial gas velocity from 10 to 20 m/min could roughly double or triple the MTZ depth. Doubling and tripling the bounding MTZ depth of 20 cm for the body of MTZ estimates made at with 10 m/min superficial velocity would extend the MTZ for a superficial velocity of 20 m/min to 40-60 cm. This bounding limit applies to all of the organic iodides that have been tested. These results also indicate that the sorption rate-limiting step is not sensitive to the superficial gas velocity; otherwise the MTZ depth would not have increased approximately in proportion to the increase in the gas superficial velocity. This further suggests that the rate limiting step is not associated with mass transfer of the sorbate to the sorbent surface, or mass transfer of the reaction byproducts from the sorbent surface, but is associated with sorption or chemical reactions on the sorbent surface or in sorbent pores. Deep bed testing at INL has established DFs of >2,000 for I, CH 3 I, and C 4 H 9 I under a range of conditions (Soelberg et al. 2021, Bruffey et al. 2019). DF does not seem to be affected by the concentration of the organic iodide in the gas stream over the range of 1 to 50 ppm. Thus, if the MTZ is accommodated in sorbent bed design for the DOG and VOG, then regulatory DFs will be met. To meet the third objective outlined in this report, these experimental data were used to update an engineering evaluation of the VOG first completed in 2016. The updated VOG design can be found in an accompanying document (Welty et al., 2021; INL- LTD-21-64587). This report finds that the VOG will decrease in both size and complexity, relative to previous designs, and will still meet regulatory requirements for all iodine forms.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development and Validation of SAM Multi-dimensional Flow Model for Thermal Mixing and Stratification Modeling

Thermal mixing and stratification in large pools or enclosures are very important phenomena that are critical to nuclear reactor safety. Because of the wide ranges of time and length scales associated with such phenomena, accurate modeling and simulation of thermal mixing and stratification remain as the key unresolved, challenging problems for reactor transient analyses. In traditional system analysis codes, simplified zero-dimensional (0-D) models are widely used for their high numerical efficiency, but they generally suffer from very limited prediction accuracies or range of applicability. Like traditional system analysis codes, the current version of the SAM code has implemented such simplified 0-D and one-dimensional mixing models. On the other hand, high-resolution Computational Fluid Dynamics (CFD) tools are often used to model complex thermal mixing and stratification phenomena. They are, however, generally numerically expensive, and they require large amounts of computational resources. It is therefore desirable to implement advanced and efficient thermal mixing and stratification modeling capabilities embedded in a system analysis code. This approach will improve the accuracy of reactor safety analyses when thermal mixing and stratification are involved, and also avoid using the large computational resources needed for high-resolution CFD analysis. Currently, with the support of the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation program, an effort has been launched to develop and implement a multi-dimensional flow model in the system analysis code SAM, and demonstrate its applications to model thermal mixing and stratification phenomena in large enclosures. The main outcomes of this research and development activity are summarized in this report, which presents an attempt to include a built-in advanced multi-dimensional flow model in a system analysis code with the focus on overcoming the simulation challenges of thermal mixing and stratification phenomena. In this report, we start with the introduction of existing SAM code capabilities to simulate thermal mixing and stratification phenomena, which is followed by a short summary of the multi-dimensional model implemented in the SAM code, including both the physical model and the Finite Element Method code implementation. In this study, two options were implemented in the SAM code to model turbulent flows: a relatively simple built-in turbulence model and an interface to accept externally computed turbulent viscosities (e.g., from a high- resolution CFD simulation). Code validation studies on this newly added capability were then carried out to compare SAM simulation results with experimental data from the SUPERCAVNA facility, which was designed to study the complex flow recirculation and thermal stratification phenomena relevant to sodium fast reactor designs. In this study, one transient and two steady-state test cases were used for code validation. Different approaches have been used to model the complex turbulence flow fields in the SUPERCAVNA facility. A highly simplified zero-equation turbulence model was first used, but it was determined that it is too simple to capture the complex turbulence flow fields in these test cases. Subsequently, the code validation continued with the use of turbulent viscosity data from high-resolution STAR-CCM+ CFD simulations to improve the accuracy of the results. Using this approach, the SAM simulation results showed very good agreement with both the SUPERCAVNA experimental data and STAR-CCM+ simulation results. In this report, we demonstrate the development, implementation, and successful validation of a multi-dimensional flow model in the SAM code, which aims to improve the simulation accuracy for complex thermal mixing and stratification phenomena. Lessons have also been learned, including that in cases where the flow fields are not well predicted by the zero-equation model, the iinput of turbulent viscosities from an external source can enhance the overall predictive capabilities needed to accurately capture complex thermal-hydraulic phenomena. Therefore, future research will be needed to further improve the code’s capabilities, such as by developing a more efficient and robust approach to capture the turbulence effects in the SAM code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗