SEARCH · Search NASA
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.
Type 347H stainless steel (347H SS), used in commercial concentrating solar power (CSP) thermal energy storage to store solar-salt at a temperature of 565°C, has been reported in the literature to be susceptible to stress-relaxation cracking (SRC). The welded heat-affected zone (HAZ) and fusion zone (FZ) of 347H SS, particularly in thick sections, are known to be susceptible to failure during post-weld heat treatment (reheat cracking). SRC could also occur after months or years under an elevated-temperature service environment. Two conditions must be present for failure to occur in the HAZ and/or FZ: 1) a modified or sensitized microstructure and 2) sufficiently high tensile residual stresses present at the elevated service temperature. The overarching goal of this project is to recommend SRC mitigation protocols to avoid susceptibility to fail through SRC at temperatures relevant for CSP. Post weld heat treatment conditions and alternative alloys are investigated as potential mitigation solutions to SRC. We used Gleeble thermomechanical simulation tests and finite element (FE) models to understand the susceptibility of 347H SS to SRC as a function of temperature, stress, and microstructure. We started with a literature review of weldability issues with 347H SS and techniques to mitigate SRC in 347H SS weldments. Next, we used Gleeble experiments to determine reheat cracking susceptibility in the simulated HAZ of 347SS weldments and an alternative alloy, 316L SS (NUCL 167 SPH) with boron added. We also performed Gleeble experiments to compare the reheat cracking susceptibility of 347H cross welded with two different fillers: E347, which is used in some commercial CSP TES tanks, and E16.8.2. We validated the experimental results with FE models of the residual stresses. We found that although the 316L (NUCL 167 SPH) is less susceptible to reheat cracking, it is slightly weaker than 347H and ASME BP&V codes limit its use to a service condition of 565°C. We also found that welds using E16.8.2 as the weld filler are less susceptible to failure than those using E347, likely due to the higher creep ductility of E16.8.2, and that it may be used as an alternative filler for repair welding of 347H welds or as the primary choice of filler for newly developed weld joints. We also found that post weld heat treatment could be a viable solution for mitigating stress in E347-347H SS thick, constrained welds, like those found in CSP tanks, and propose several options for mitigating SRC in existing and future TES tanks.
Nearly 20 years since the origin of the infamous 'Kord Smith Challenge', and after 10 years and hundreds of millions of dollars invested by the US Department of Energy into the development of the Virtual Environment for Reactor Applications (VERA), the capability to accurately simulate a modern pressurized water reactor (PWR) fuel cycle at the fuel rod level with high-fidelity Monte Carlo (MC) stochastic transport methods still evades us. The deterministic methods in VERA make approximations that allow it to be successfully benchmarked against measured data from hundreds of PWR fuel cycles, but none of these benchmarks provided information at the local fuel rod level or allowed quantification of these approximations' impact for modern PWR reload designs. Using MPACT to establish the estimated isotopic distributions and thermal hydraulic conditions in a recent and challenging reactor design, a 2D quarter-core benchmark problem was established and solved consistently with both MPACT, the 51-energy group deterministic transport solver in VERA, and Shift, the continuous-energy MC transport solver in VERA. Impacts of energy groups, ray spacing, and scattering treatment were evaluated. Though MPACT has been previously shown to perform well compared to MC methods for initial Cycle 1 cores with only fresh fuel, this is the first benchmark of MPACT for a challenging modern reload core with depleted fuel and burnable absorbers. This work establishes additional credibility for the methods in MPACT and highlights the efficiency of these methods relative to those in use by MC-based tools. (authors)
During the first phase of the Consortium for Advanced Simulation of Light Water Reactors (CASL) program, the Virtual Environment for Reactor Applications (VERA) was developed with a focus on capabilities for high-fidelity, multiphysics simulation of pressurized water reactors (PWRs). During this development effort, a set of progression problems was created ranging from smaller pin cell calculations to larger 3D full-core calculations. These progression problems helped to guide the development of the software and served as benchmarks against which to test VERA. Since 2019, efforts have been made to extend the capabilities of VERA to model boiling water reactors (BWRs). Because BWR simulations come with many unique challenges, a set of BWR progression problems was developed to aid in this new effort. The BWR progression problems range from 2D lattice calculations to 3D mini-core problems, and reference neutronic solutions were computed using continuous-energy Monte Carlo codes. MPACT, one of the neutronics code in VERA, was benchmarked using the BWR progression problems. The code is capable of computing solutions to all problems. The eigenvalues computed by MPACT agree well with the Monte Carlo reference solutions. (authors)
The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling.
VERA core simulator capabilities to predict spent fuel decay heat have been explored. The comparison with the Serpent2 and NEWT/TRITON sequence from the SCALE package, performed on the NPP Krško fuel test cases, has enabled independent verification of the successful implementation of the VERA depletion/decay calculation sequence applied on the extended nuclide set. Since the VERA is a 3-D core simulator, the focus in the rest of the paper was on the investigation of the 3-D effects and non-linearity of the soluble boron concentration averaging process, which are difficult to assess with most of available computational tools. A difference in the decay heat prediction between the 3-D and 2-D model approaches was analyzed. In addition, the impact of the Inconel grids on the decay heat was determined. Averaging the concentration of soluble boron over the fuel irradiation period has been found to be a viable approach, since it results in slight decay heat overprediction, which is considered conservative for most applications. Finally, demonstration on the Watts Bar Unit 1 cycle 1 depletion has shown formidable VERA capabilities to accurately predict spent fuel decay heat.
The U.S. Department of Energy (DOE) Office of Integrated Waste Management is planning for the eventual transportation, storage, and disposal of spent nuclear fuel (SNF) and high-level radioactive waste (HLW) from nuclear power plant and DOE sites. The Stakeholder Tool for Assessing Radioactive Transportation (START) is a web-based, geospatial decision-support tool developed for evaluating routing options and other aspects of transporting SNF and HLW, covering rail, truck, barge, and intermodal infrastructure and operations in the continental United States. The verification and validation (V&V) process is intended to independently assess START to provide confidence in the ability of START to accurately provide intended results. The V&V process checks the START tool using a variety of methods, ranging from independent hand calculations to comparison of START performance and results to those of other codes. The V&V activity was conducted independently from the START development team with opportunities to provide feedback and collaborate throughout the process. The V&V analyzed attributes of transportation routes produced by START, including route distance and both population and population density captured within buffer zones around routes. Population in the buffer zone, population density in the buffer zone, and route distance were all identified as crucial outputs of the START code and were subject to V&V tasks. Some of the improvements identified through the V&V process were standardizing the underlying population data in START, changing the projection of the population raster data, and changes to the methodology used for population density to improve its applicability for expected users. This collaboration also led to suggested improvements to some of the underlying shape file segments within START. A formal series of test routes went through the V&V process from START version 3.2.1 following the implementation of the recommended improvements to START. The V&V of the buffer zone population and distance reported in START were achieved by utilizing a custom workflow developed in QGIS. This platform was chosen because it is fully independent of START. The buffer zone population and route distances in START showed excellent agreement to independent V&V test results. Over 200 route test cases were run; in all cases, the percent difference in the population within the buffer zone was less than +/- 5%, and the majority of cases were below +/- 1%. These cases were also duplicated in ArcMap, which is Environmental Systems Research Institute’s (ESRI) GIS desktop application. The ArcMap results also show good agreement; however, they cannot be considered independent of START because START also uses ESRI software. The V&V tasks are currently addressing population density and estimated radiation dose in START as well as implementing further automation of the V&V process to allow larger suites of test routes to be run more easily and quickly by less experienced GIS users. These improvements to the V&V process will make it a valuable tool for the development team and the stakeholders so that each new iteration of START can easily be checked for quality and consistency. This paper and presentation will describe the V&V process, findings and implications, and ongoing V&V activities.
Modeling and simulation (M&S) methods are able to predict uncertainties in experimental parameters (e.g., power and fission density) during irradiation. A shortfall exists in predicting how sensitive some of the parameters will behave during the experimental process. Sensitivity and Uncertainty Quantification (SUQ) is critical in support of qualification and licensing reactor fuels. The application of a method to quantify the uncertainty in these experiments is critical to the prediction of their performance. In this work, we propose the use of a polynomial chaos expansion (PCE) method to quantify the sensitive parameters in these simulations and, in an extension, their experimental surrogates. We propose to perform M&S using PCE uncertainty quantification on a previously irradiated fueled experiment in order to provide a validation case for Griffin and expand its use as a verification and validation (V&V) tool for experiments with a neutronics component. Griffin is an advanced, deterministic neutronics analysis code built using the MOOSE (multiphysics object-oriented simulation environment) framework which can provide state-of-the-art neutronic analysis on M&S of experiments. We will use the stochastic tools module (STM) in MOOSE to provide PCE uncertainty quantification on the proposed experimental setup. Idaho National Laboratory (INL) does not yet have an in-house developed code with V&V approval for experiments performed on-site; this work would provide a necessary addition of support for experiments performed at INL. The Nuclear Regulatory Commission (NRC) has explicitly requested uncertainties in calculated values such as fuel power and burnup, and the development of this capability would benefit the relationship between INL and the NRC.
Aeroelastic modeling (AM) is the primary methodology for structural and performance assessment of any wind turbine whereby providing an understanding of the impact of design parameters on its loading and power response before witnessing it in the field. Despite these advantages, the use of AM in the Distributed Wind Technology (DWT) sector is limited, especially within the less established manufacturers. This project represents an in-depth assessment of the status of AM and its role within the Standards for the DWT industry. The study gathered input and feedback from a large number of national and international stakeholders, reviewed technical strengths and weaknesses of the current edition of the design standards, analyzed recent industry workshops' and meetings' minutes, collected publicly available AM templates, and provided an evaluation of the existing AM codes. The study achieved several goals including providing strategies for the load assessment categorization of turbines based on rotor swept area and archetype, and guidance for AM verification and validation (V&V), which includes discussions of measurement requirements and a sample test-plan useful for future V&V campaigns and design standard development. This document summarizes the different tasks conducted in the course of the project and highlights the steps required to improve the AM adoption based on a multifaceted approach that encompasses: 1) augmenting AM software capabilities, 2) publishing AM best-practices and design-basis, 3) creating new model templates, 4) providing guidance for V&V of codes and specific turbine models leveraging field testing best-practice, and 5) addressing weaknesses in the current standards. Many of the future objectives identified in this study could leverage NREL's upcoming testing campaigns of three modern distributed wind turbines. Recommendations within this study will advance the value and the ease-of-use of AM, thereby allowing the industry to better capitalize this underutilized tool resulting in a more efficient design process, an easier path to certification, and overall better and more distributed reliable wind turbine products.
As part of a cooperative research and development agreement between the US department of Energy's Oak Ridge National Laboratory and Holtec International subsidiary SMR LLC, core simulator models of the SMR-160, a light water-cooled and -moderated, natural circulation, small modular reactor, were developed and compared. The model results indicate pin-by-pin power histories, critical boron concentration, and other key performance indicators useful for developing multi-cycle core loading patterns. One model uses industry standard nodal methods with a simple thermal-hydraulics solution, and the other uses an advanced deterministic neutronics solver with a coupled subchannel calculation. The goals of the comparison are to provide the reactor designer with a detailed and coupled core physics reference solution, inform the reactor designer of nodal model best practices based on the code-to-code comparison, and set up a workflow for performing subsequent detailed core simulator calculations. The results of this work will be used in the development of coupled neutronic, thermal hydraulic, and fuel performance models for advanced analysis of the SMR-160. (authors)
Effects of temperature and flow gradients in Westinghouse designed three-loop Pressurized Water Reactor (PWR) cold legs, the piping between the main coolant pump and the reactor vessel, were evaluated using Computational Fluid Dynamics (CFD) code STAR CCM+ and coupled neutronic and thermal-hydraulic (T/H) code system VERA. In the parametric study, several symmetric and asymmetric temperature gradients that were significantly larger than those observed from plant measurements were applied to the cold leg inlets for comparison with the base case without any temperature gradient. A CFD model using the STAR-CCM+ code was developed for a portion of the RCS region between the Reactor Coolant Pump (RCP) and the core inlet based on previously validated modeling approach. The CFD simulation results were processed for the temperature and flow rate distributions at the core inlet as input to the VERA calculations. The VERA code system consists of COBRA-TF (CTF) for thermal-hydraulics, MPACT for reactor physics and neutron transport, and ORIGEN for isotopic depletion. The VERA model was for depletion calculations of a high-burnup loading pattern with the reactor core in pin-by-pin and subchannel resolution. The results of the study indicate that the postulated temperature gradients within the PWR cold legs do not result in any significant changes in the core inlet temperature distributions and the core power distributions during the reactor operation. (authors)
Recently US and European nuclear data libraries have been released, namely ENDF/B-VIII.0 and JEFF-3.3 libraries, which are the result of years of evaluation and validation work in both communities. Consequently, efforts have been made to validate the evaluations in calculations of integral experiments (critical benchmarks, ..etc). Nevertheless, less stringent testing and validation efforts were performed on burnup applications before releasing the libraries. The presented work focusses on the testing of these recent nuclear data libraries for burnup calculations on two benchmarks at the pin and at the assembly level. Monte-Carlo depletion calculations were performed using the VESTA 2.2 code. The K{sub ∞} results between nuclear data libraries are compared. A strong k{sub ∞} bias is observed with burnup using both JEFF-3.3 and JEFF-4T0 compared to all other libraries, and especially ENDF/B-VIII.0, consisting in a strong k{sub ∞} over-estimation at low burnup and a high under-estimation at high burnup. JEFF-3.3 {sup 235}U and {sup 239}Pu evaluations mainly explain this result, as well as fission yields. New {sup 235}U and {sup 239}Pu evaluations were proposed for JEFF-4T0, but they do not address totally the bias issue, even if JEFF-4T0 {sup 239}Pu allows slightly reducing the k{sub ∞} underestimation at high burnup.
With growing renewables penetration, there is increased interest in flexible power operation for nuclear reactors. For multi-unit SMRs, in particular the NuScale SMR, which is an integral pressurized water reactor, there are opportunities to optimize flexible power operation across multiple units to limit the degradation of structural and control components. Here, we focus on degradation of in-core components, specifically the control rods and reactor pressure vessel. To perform these studies a high-fidelity, multi-cycle representation of the NuScale SMR is required, with a detailed representation of the structural and control components. To this end, the NuScale SMR has been modelled using the Virtual Environment for Reactor Applications (VERA) software. The entire transition to equilibrium is simulated, from Cycle 1 through to the equilibrium cycle. The equilibrium cycle model shows a good agreement with the NuScale design certification application (DCA) results, with differences attributable to a combination of using public domain data for the present study, and methodological differences. K-effective, power distributions, reactivity coefficients, and boron letdown curves are compared and all found to closely match. This shows that the VERA model is suitable for further studies. (authors)
During the last two years, the Virtual Environment for Reactor Applications (VERA) has been extended to simulate boiling water reactors (BWRs). The thermal hydraulic effects present in BWRs are far more complex than those in pressurized water reactors (PWRs). Therefore, the runtime is significantly increase and convergence behavior worse compared to PWRs. Most of the additional runtime is spent during the first few coupled iterations, when the power shape is still rapidly evolving, dramatically affecting the thermal hydraulics (TH). To alleviate the increased computational expense, a multiscale TH coupling approach was developed in VERA in which a highly efficient, simplified TH model is solved for several coupled iterations until the power shape is partially converged. The simplified solution is used to inform the assembly-wise flow distribution in the high-fidelity TH solver CTF, reducing the amount of work required to properly balance the pressure drop in each channel and reducing runtime for the couple calculation. Because the final TH calculations are performed with CTF, there is ultimately no impact on the accuracy of the converged solution. This paper presents the details of this multiscale TH coupling approach, along with results for a 4x4 array of GE-14 fuel bundles and whole-core coupled simulations of the Hatch core. These two cases show that the multiscale approach dramatically reduces the runtime of coupled BWR simulations. (authors)
Explore the source record for details and available documents.
Manuscript to submit to the SPIE Optics + Photonics 2022 conference, August 21–25, 2022, San Diego, California, for publication in the SPIE Digital Library (Proc. SPIE) (https://spie.org/conferences-and-exhibitions/optics-and-photonics?SSO=1).
Reactor licensing centers on the protection of the public and environment from the inadvertent release of radioactive material. Therefore, mechanistic source term analysis, or the realistic evaluation of radionuclide transport from the source to the environment for specific transient scenarios, is vital to reactor licensing efforts. Developed to resolve a gap in mechanistic source term modeling capabilities for sodium fast reactors (SFRs), the Simplified Radionuclide Transport (SRT) code created by Argonne National Laboratory (Argonne) is now utilized by advanced reactor vendors, universities, and research institutions to support a multitude of design, licensing, and research efforts. Recently, SRT version 2.1 was released, which includes improvements to code models and verification and validation (V&V) suite to support the SRT user community. The following work provides an overview of the improvements made as part of the release of SRT version 2.1. This effort is supported by the U.S. Department of Energy Office of Nuclear Energy (DOE:NE) Advanced Reactor Technologies (ART) Fast Reactor Program (FRP), as part of the program’s support of national laboratory design and safety analysis computer codes utilized by the fast reactor industry. The expansion of SRT code capabilities and improvements to code V&V associated with version 2.1 are in response to user requests and lessons learned from recent source term analyses performed by Argonne and advanced reactor vendors. They also align with the evolving role of SRT, from research and development tool to software utilized for reactor licensing calculations. The report is structured in alignment with the code improvements, as outlined in Figure 1-1. Section 2 provides background information on SRT, including its history, capabilities, and utilization. Section 3 details new code capabilities as part of version 2.1, while Section 4 focuses on the expansion of the code’s V&V suite. Lastly, Section 5 provides a summary and discussion of next steps.
A fully implicit particle-in-cell method for handling the v ∥ -formalism of electromagnetic gyrokinetics has been implemented in XGC. By choosing the v ∥ -formalism, here we avoid introducing the nonphysical skin terms in Ampère's law, which are responsible for the well-known “cancellation problem” in the p ∥ -formalism. The v ∥ -formalism, however, is known to suffer from a numerical instability when explicit time integration schemes are used due to the appearance of a time derivative in the particle equations of motion from the inductive component of the electric field. Here, using the conventional δf scheme, we demonstrate that our implicitly discretized algorithm can provide numerically stable simulation results with accurate dispersive properties. We verify the algorithm using a test case for shear Alfvén wave propagation in addition to a case demonstrating the ion temperature gradient-kinetic ballooning mode (ITG-KBM) transition. The ITG-KBM transition case is compared to results obtained from other δf gyrokinetic codes/schemes, whose verification has already been archived in the literature.