Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Benchmark”

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.

At least 1,675 records · Page 93

An Overview of the State-of-the-Art Reactor Consequence Uncertainty Assessment Accident Progression Insights

The U.S. Nuclear Regulatory Commission (NRC) with Sandia National Laboratories (Sandia) have completed three uncertainty analyses (UAs) as part of the State-of-the-Art Reactor Consequence Analyses (SOARCA) program. The SOARCA UAs included an integrated evaluation of uncertainty in accident progression, radiological release, and offsite health consequence projections. The UA for Peach Bottom, a boiling-water reactor (BWR) with a Mark I containment located in the State of Pennsylvania, analyzed the unmitigated long-term station blackout SOARCA scenario. The UA for Sequoyah, a 4-loop Westinghouse pressurized-water reactor (PWR) located in the State of Tennessee, analyzed the unmitigated short-term station blackout SOARCA scenario, with a focus on issues unique to the ice condenser containment and the potential for early containment failure due to hydrogen deflagration. The UA for Surry, a 3-loop Westinghouse PWR with a sub-atmospheric large dry containment located in the State of Virginia, analyzed the unmitigated short-term station blackout SOARCA scenario including the potential for thermally-induced steam-generator tube rupture. These three UAs are currently documented in three NUREG/CR reports. This report provides input to planned NRC documentation on the insights and findings from the SOARCA UA program. The purpose of the summary report is to provide a useful reference for regulatory applications that require the evaluation of offsite consequence risk from beyond design basis event severe accidents. This report focuses on the accident progression and source term insights developed from the MELCOR analyses. MELCOR is the NRC's best-estimate, severe accident computer code used in the SOARCA UAs. In anticipation of the SOARCA UA insights work, NRC and Sandia benchmarked the response of the Peach Bottom model to selected reference calculations from the Peach Bottom SOARCA UA. Peach Bottom was the first SOARCA UA performed and was completed in 2015 using the MELCOR 1.8.6 code. The PWR SOARCA UAs evolved the original methodology and utilized the updated MELCOR 2.2 computer code. The Peach Bottom model has been systematically updated for other NRC research efforts and has been updated to MELCOR 2.2. computer code. The findings from the new reference calculations using the updated model with the MELCOR 2.2 code are also integrated into the report. A second objective is an assessment of the applicability of the results to the other nuclear reactors in the U.S. As the key findings are reviewed, judgments are presented on the applicability of the results to other U.S. nuclear power plants. An important objective of the SOARCA program relied on high- fidelity plant-specific modeling. However, the nature of the insights and conclusions allowed judgements to be made on the applicability of the various insights to the same general classification of plant (i.e., BWR or PWR) or the entire fleet of plants. Finally, the results from the SOARCA UA accident progression calculations contain a wealth of information not previously documented in the NUREG/CRs. This report includes new but related information that can be used to benchmark past or support future regulatory decisions related to severe accidents. The new work includes a benchmark of the NUREG-1465 licensing source term definitions, the variability of key accident progression events and timing to radionuclide release, and an improved understanding of the timing and source terms from consequential steam generator tube ruptures. iii ACKNOWLEDGEMENTS The Sandia authors gratefully acknowledge the significant technical and programmatic contributions from the NRC SOARCA team which are reflected throughout the report. Dr. Tina Ghosh has been involved throughout the SOARCA UAs, providing the primary managerial and technical oversight. The long lists of NRC and Sandia contributors from the SOARCA UAs are cited in the three NUREG/CRs and are also gratefully acknowledged by the small team of authors compiling the results of their efforts. Significant technical contributions, advice, and reviews were provided by Dr. Hossein Esmaili, Dr. Alfred Hathaway, and Dr. Edward Fuller (retired) of the NRC. Dr. Randal Gauntt (retired), Mr. Patrick Mattie, Mr. Joseph Jones (retired), and Dr. Doug Osborn from Sandia are recognized as the SOARCA UA managers guiding the past efforts. There is a comparable list of project managers at the NRC including Ms. Patricia Santiago, Dr. Salman Haq, and Mr. Jon Barr. Sadly, we have lost Mr. Charlie Tinkler and Mr. Robert Prato, who were important contributors to the original SOARCA project. Finally, Mr. Kyle Ross and Mr. Mark Leonard have also retired but were significant technical contributors. Mr. Kyle Ross was the technical lead on all three SOARCA UAs and the original pressurized water reactor SOARCA study. Mr. Leonard was the technical lead on the original boiling water reactor SOARCA study and a key contributor to the first Peach Bottom SOARCA UA. iv

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗

Evaluation and improvement of the parameterization of aerosol hygroscopicity in global climate models using in-situ surface measurements (Final Report)

Aerosols are tiny particles suspended on the atmosphere that can interact with incoming solar radiation and affect the Earth radiative budget. They do so by scattering and absorbing solar radiation, and these properties vary depending on the aerosol size and chemical composition. Moreover, by taking up water from the surrounding air, hygroscopic aerosol particles will grow in size and change their chemical composition, thus modifying their optical properties (scattering and absorption) and their final impact on radiative forcing calculations. An accurate knowledge of aerosol hygroscopicity is crucial for estimating the net radiative impact of aerosols. We took a three-pronged approach to improve our understanding of aerosol hygroscopicity and how it is implemented in Earth system models. In the first part of our project, we developed a benchmark dataset from existing aerosol hygroscopic growth measurements made by tandem nephelometer humidogram systems. We analyzed, using a standardized methodology, observations from 26 in-situ stations around the globe. Measurement data was collected from multiple data providers, reviewed and harmonized to create a consistent dataset of the scattering enhancement factor due to aerosol water uptake. This dataset is archived in several publicly available databases for use by interested researchers. In the second part of the project, we used the benchmark hygroscopicity dataset to perform a global study on aerosol hygroscopicity and aerosol optical properties. Measurements show a global picture of scattering enhancement with larger values for Arctic and marine sites and lower for urban and desert sites. We assessed the RH dependence of aerosol radiative forcing and showed that the overall effect of aerosol hygroscopicity on DARF is an increase in the absolute forcing effect (negative sign) by a factor of up to 4 compared to dry conditions (RH<40%). Finally, we explored using aerosol single scattering albedo (SSA) and scattering Angstrom exponent (SAE) as possible proxies for aerosol hygroscopicity. SSA showed more promise as a surrogate for the scattering enhancement factor than SAE, but neither was ideal. In the third part of the project we evaluated the output of ten Earth system models (ESMs) against the benchmark hygroscopicity dataset. ESMs utilize various schemes for aerosol hygroscopicity which had not been previously tested against observations on a global scale. We found that ESMs currently overestimate scattering enhancement due to hygroscopic growth. Model parameterizations of hygroscopicity and model chemistry are two main factors driving the observed diversity in hygroscopicity simulations among the models. In addition, our study makes several suggestions for modelers, including improving the parameterizations of organic and sea salt aerosol hygroscopicity. Future hygroscopicity model evaluation experiments should include the model data related to particles size which was not available for our study.

54 ENVIRONMENTAL SCIENCES↗

Algal Biomass Conversion to Fuels via Combined Algae Processing (CAP): 2020 State of Technology and Future Research

The annual State of Technology (SOT) assessment is an essential activity for platform research conducted under the Bioenergy Technologies Office. It allows for the impact of research progress to be quantified in terms of economic improvements in the overall biofuel production process for a particular biomass processing pathway, whether based on terrestrial or algal biomass feedstocks. As such, initial benchmarks can be established for currently demonstrated performance, and progress can be tracked towards out-year goals to ultimately demonstrate economically viable biofuel technologies. NREL's algae state of technology benchmarking efforts focus both on front-end algal biomass production and separately on back-end conversion to fuels through NREL's "combined algae processing" (CAP) pathway. The production model is based on outdoor long-term cultivation data, enabled by comprehensive algal biomass production trials conducted under the Development of Integrated Screening, Cultivar Optimization, and Verification Research consortium (DISCOVR) efforts, driven by data furnished by Arizona State University (ASU) at the Arizona Center for Algae Technology and Innovation (AzCATI) testbed site. The CAP model is based on experimental efforts conducted under NREL research and development projects. This report focuses on back-end conversion of algal biomass through the CAP pathway, highlighting the 2020 updates to minimum fuel selling price (MFSP). Notably, this analysis marks the first time that polyurethane (PU) is included as a value-added coproduct in the SOT case. Relative to prior 2019 benchmarks, the 2020 CAP SOT achieved further MFSP reductions of over $4/gallon gasoline equivalent, reflective of (a) the inclusion of PU co-production, (b) reduced SOT biomass production costs, and (c) higher pretreatment sugar yield data.

09 BIOMASS FUELS↗

Biochemical Conversion of Lignocellulosic Biomass to Hydrocarbon Fuels and Products: 2020 State of Technology and Future Research

The annual State of Technology (SOT) assessment is an essential activity for biochemical platform research. It allows the impact of research progress to be quantified in terms of economic improvements in the overall cellulosic biofuel production process for a particular conversion pathway. As such, initial benchmarks can be established for currently demonstrated performance and progress can be tracked towards out-year goals to ultimately demonstrate cost-competitive cellulosic biofuel technology. The purpose of this report is to benchmark the latest experimental developments across a number of potential bioconversion pathways as quantified by modeled minimum fuel selling prices (MFSPs), as a measure of current status relative to those final targets. For this state of technology, TEA models were run for two separate biological conversion pathways to fuels, based on available data for integrated biomass deconstruction and hydrolysate processing; namely carboxylic acids (primarily butyric acid) and diols (2,3-butanediol [BDO]), reflecting NREL's recently-published 2018 biochemical design report focused on those two pathways. The models were run across three scenarios for lignin utilization, namely combustion, conversion to coproducts based on "base case" performance with biomass hydrolysate, and conversion to coproducts based on "high" performance demonstrated with model lignin monomer components. Relative to prior 2019 benchmarks, the 2020 SOT achieved further MFSP reductions of approximately $1/gallon gasoline equivalent, primarily attributed to improvements in pretreatment that led to both cost reductions and yield improvements.

09 BIOMASS FUELS↗

Algal Biomass Production via Open Pond Algae Farm Cultivation: 2020 State of Technology and Future Research

The annual State of Technology (SOT) assessment is an essential activity for platform research conducted under the Bioenergy Technologies Office. It allows for the impact of research progress to be quantified in terms of economic improvements in the overall biofuel production process for a particular biomass processing pathway, whether based on terrestrial or algal biomass feedstocks. As such, initial benchmarks can be established for currently demonstrated performance, and progress can be tracked towards out-year goals to ultimately demonstrate economically viable biofuel technologies. NREL's algae state of technology benchmarking efforts focus both on front-end algal biomass production and separately on back-end conversion to fuels through NREL's "combined algae processing" (CAP) pathway. The production model is based on outdoor long-term cultivation data, enabled by comprehensive algal biomass production trials conducted under Development of Integrated Screening, Cultivar Optimization, and Verification Research (DISCOVR) consortium efforts, driven by data furnished by Arizona State University (ASU) at the Arizona Center for Algae Technology and Innovation (AzCATI) testbed site. The CAP model is based on experimental efforts conducted under NREL research and development projects. This report focuses on front-end algal biomass production, documenting the pertinent algal biomass cultivation parameters that were input to the NREL open pond algae farm model based on the latest DISCOVR cultivation performance data. Relative to prior 2019 benchmarks, the 2020 algae farm SOT achieved further reductions in modeled minimum biomass selling price of approximately $70/ton (roughly 10%), primarily attributable to a 16% demonstrated improvement in annual average cultivation productivity rates.

09 BIOMASS FUELS↗

Transport in Stochastic Media with Random Chord Length Distributions

Thermal radiation transport computations in binary Markovian random mixtures rely almost exclusively on the Levermore-Pomraning (LP) model which is obtained by applying a heuristic closure to the ensemble averaged random medium transport equation. The validity of this model has been extensively tested by comparing numerical results over a broad parameter range (material types and mixing parameters) against benchmark solutions in planar geometry. The conditions under which the LP-model provides useful results and when it breaks down are now well established, but work to date has been largely restricted to homogeneous mixing statistics, i.e., the mean chord lengths of both materials are taken to be spatially constant. In recent work, this limitation was relaxed by allowing the mean chord lengths and, in a consistent fashion, the volume fractions in the LP-model to vary spatially and in a follow-up investigation benchmark solutions were obtained by ensemble averaging results over material realizations sampled from a nonhomogeneous Poisson process (NHPP). Numerical experiments in rod geometry with specifically linear and quadratic spatial dependence of chord lengths showed that the material averaged radiation intensities vary nonmonotonically with depth into the medium, in stark contrast to solutions obtained assuming uniform chord lengths. Moreover, depending on the local optical thickness and strength of scattering, the LP-model results showed locally more nuanced deviations from the benchmark solutions than was the case with constant chord lengths. These limited numerical investigations highlight the nontrivial qualitative and quantitative consequences of nonhomogeneous mixing statistics, in particular that closure approximations may not be uniformly valid or invalid over the problem domain.

42 ENGINEERING↗

Closing the Loop on Automotive Carbon Fiber Prepreg Manufacturing Scrap

The project demonstrated how to “close the loop” on carbon fiber by integrating industrial carbon fiber scrap into new functional components in an automotive lightweighting application. The project serves as a validation of discontinuous recycled carbon fiber in a commercial context, while generating comprehensive material data throughout the production chain. To this end, the project exhibited increasing complexity as material evaluation progressed from benchtop to commercial scale through full-scale part production, with key material properties thoroughly characterized throughout the process. Of particular focus was the form of the fiber that was fed into compounding, as recycled fiber has historically been problematic to feed at commercial-scale. Carbon fiber is energy-intensive to manufacture, so reuse of existing fiber material can reduce costs and increase sustainability. Additionally, by integrating recovered short fiber into a thermoplastic, regrind processes can be used to provide feedstock for later generations of product. While regrind plastics are not “infinitely recyclable”, reusing the manufacturing scrap over several generations of products can greatly increase material sustainability and lower the fractional embodied energy of each successive product. As such, this project supports the IACMI technical goals of (1) 25% lower carbon fiber-reinforced polymer (CFRP) cost, (2) 50% reduction in CFRP embodied energy, and (3) 80% composite recyclability into useful products. The initial stage of the project involved down selecting surface treatment (sizing) chemistries. Sizing evaluations were performed on Vartega’s chemically recycled intermediate modulus fiber along with standard modulus dry scrap which was oversized with sizing provided by Michelman. More dramatic improvements from sizing were found on the standard rather than the intermediate modulus fibers. The strength of the chemically recycled individual fibers were evaluated by Michelman and ORNL through single fiber testing and found to be comparable to similarly evaluated virgin fibers. UDRI’s mechanical testing on injection molded test specimens identified similar mechanical properties and fiber distribution relative to benchmark specimens. Additional surface chemistry tests and visualizations were performed by the Colorado School of Mines to confirm close conformance between the benchmark and recycled-fiber specimens. As the mechanical test results exceeded the 80% threshold established as the go/no-go(GNG), the project scale was increased to use commercial-scale equipment that would both better characterize the manufacturing utility of the target product format and allow qualitative assessment of a complex commercial part. An upscaled compounding evaluation was performed with a 27 mm twin-screw compounding extruding using oversized standard modulus fibers that were formatted to improve bulk solids transfer. The project team anticipated that milestone mechanical benchmarks could be achieved given the favorable performance of the sized standard modulus material identified in the initial micro-compounding trials. While the mechanical performance did meet the milestone target for that phase of the project, mechanical properties for this standard modulus-based compound were still less than those of the Ford specification. To compare the performance, the project team oversized intermediate modulus dry fibers and compounded them with the project resin at BASF using a 40 mm compounder. Test specimen mechanical performance exceeded the targets laid out in both the project milestone and the Ford specification. A series of prototype parts were successfully molded, albeit with instances of short shot components due to the high thermal conductivity of the carbon fiber compared to glass fiber for which the prototype tool was designed for. The project demonstrated that recycled carbon fiber is a viable option in fiber reinforced compound, providing greatly increased strength and modulus for applications that require them. The “agglomerated” format that facilitated effective bulk solids transfer of recovered fiber showed promise for industrial application.

36 MATERIALS SCIENCE↗

Evaluation of China Experimental Fast Reactor Start-up Tests (Final Report, Revision 1)

This report documents the results and observations on the Coordinated Research Project (CRP) of the International Atomic Energy Agency (IAEA) on “Neutronics Benchmark of CEFR Start-Up Tests.” The China Experimental Fast Reactor (CEFR) is a 65MWt sodium-cooled fast reactor with highly enriched uranium oxide fuel. The reactor achieved first criticality in 2010, and a series of start-up tests were conducted to measure various reactor physics parameters. In 2018, the IAEA launched the CRP for validation and qualification of member states' computation capabilities in the field of fast reactor simulation by utilizing the measured data in the CEFR start-up test. Twenty-nine international organizations from eighteen member countries, including Argonne National Laboratory, have participated in the CRP. The CEFR start-up tests offer a rare opportunity to validate U.S. nuclear engineering software because it is a well-specified benchmark with corresponding measurements for a recently built fast reactor starting up with a known fuel composition (fresh fuel). Previous fast reactor validation benchmarks frequently involve reactors that have already been started up and contain irradiated fuel that is difficult to characterize with high certainty.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

FY20 Proxy App Suite Release: Report for ECP Proxy App Project Milestone ADCD-504-10

Version 4.0 of the ECP Proxy App Suite is practically unchanged from the previous release. The current set of proxies has proven useful for many aspects of benchmarking and co-design and we see little reason to alter the suite. Although there have been few changes to the ECP suite, the team has been hard at work in other areas. In the area of Machine Learning (ML) we have now created a separate proxy suite dedicated to this scientific applications of ML. The suite includes: miniGAN (Generative Adversarial Networks), miniRL (Reinforcement Learning), CRADL (inline inference), Cosmoflow-Benchmark (Convolutional Neural Network), and MLPerf-DeepCam (Climate Segmentation Benchmark). Section 3 contains more information about these proxies as well as the principles that are guiding the development of the suite. We have surveyed available proxies for several application domains including Computational Fluid Dynamics, Quantum Chemistry, Quantum Computing Simulation, Molecular Dynamics, Monte Carlo Transport, and Density Functional Theory to identify gaps in proxy coverage. Several new proxy apps are either already available or will be released soon to fill these gaps. Section 4 provides full details. Finally, section 5 reports on our continued collaboration with the ECP Continuous Integration (CI) effort to use proxy apps to help identify problems and roadblocks to cross-lab CI. We have also assisted the El Capitan Center of Excellence (COE) to stand up a CI system that will be used to test releases of HPE and AMD software stacks. We hope that the COE effort can serve as a model for ECP by showing how the Proxy App Team can work with AD and ST teams identify critical features, kernels, patterns, etc. and incorporate them into a CI system that will help ensure that Frontier and Aurora will provide those needed capabilities.

97 MATHEMATICS AND COMPUTING↗

Biochemical Conversion of Lignocellulosic Biomass to Hydrocarbon Fuels and Products (2021 State of Technology and Future Research)

The annual State of Technology (SOT) assessment is an essential activity for biochemical platform research. It allows the impact of research progress to be quantified in terms of economic improvements in the overall cellulosic biofuel production process for a particular conversion pathway. As such, initial benchmarks can be established for currently demonstrated performance and progress can be tracked towards out-year goals to ultimately demonstrate cost-competitive cellulosic biofuel technology. The purpose of this report is to benchmark the latest experimental developments across a number of potential bioconversion pathways as quantified by modeled minimum fuel selling prices (MFSPs), as a measure of current status relative to those final targets. For this state of technology, TEA models were run for two separate biological conversion pathways to fuels, based on available data for integrated biomass deconstruction and hydrolysate processing; namely carboxylic acids (primarily butyric acid) and diols (2,3-butanediol [BDO]), reflecting NREL's recently-published 2018 biochemical design report focused on those two pathways. The models were run across three scenarios for lignin utilization, namely combustion, conversion to coproducts based on "base case" performance with biomass hydrolysate, and conversion to coproducts based on "high" performance demonstrated with model lignin monomer components. A key improvement reflected in the 2021 SOT is centered around making use of the latest lignin conversion data, which over the past year focused primarily on production of ß-ketoadipate (BKA) as a more optimal molecule compared to the closely-related adipic acid coproduct of prior recent focus, both in terms of superior product properties and biology, as well as reduced processing complexity (reducing two steps for sequential production of muconate followed by hydrogenation to adipic acid down to a single step for direct production of BKA). This update translated to a roughly 17% increase in mass yield of final coproduct output at a nearly four-fold increase in fermentation productivity on lignin monomers relative to prior 2020 SOT benchmarks for muconic/adipic acid production.

09 BIOMASS FUELS↗

Continued Verification of MOOSE Structural Mechanics Tools for Modeling Core Bowing Phenomena in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY22, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21. Several benchmark verification examples were selected from the IAEA verification and validation report. These examples involve clusters of ducts representative of a sector of a hexagonal reactor core which bow into each other and cause contact and load pad elevations, as well as single ducts subjected to irradiation fields undergoing swelling and subsequent bowing. The MOOSE-based results were compared to both IAEA benchmark participants’ results, analytic equations as available, and NUBOW-3D, a beam model code developed by Argonne National Laboratory. In every case, the MOOSE results agreed with other simulations results, providing additional verification basis of the tools for this particular physics application.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CERBERUS: CED-2 Final Design Report

The goal of the Critical Experiment Reflected By copper to bEtteR Understand Scattering [CERBERUS] is to design a critical experiment that maximizes sensitivities to copper (Cu) reactions, particularly in the intermediate energy region (0.625 eV – 100 keV) as well as the 100 – 600 keV energy region. Despite the number of experiments evaluated in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook1 , there is still a lack of benchmarks sensitive to neutrons in the intermediate energy region. The ZEUS series was designed specifically to provide data in the intermediate energy range. The ZEUS series used Cu reflectors to allow for the construction of an intermediate energy system that could be constructed within a reasonable size. While the reflection provided by Cu reduces the size of the system, it creates a system that is very sensitive to the angular scattering in Cu. A neutron scattering off Cu at some angles reflects back into to the system, while at other angles, it is lost entirely. Of the ICSBEP benchmarks sensitive to neutrons in the intermediate energy region, very few are also sensitive to Cu in that region. Improving Cu nuclear data is important outside of the ZEUS series, because it is present in many bronze and aluminum alloys, which are used in various nuclear operations such as motors, wires, and some containers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Subcriticality Measurements for a Changing Concentration Uranyl Nitrate Solution Tank by Californium Source-Driven Noise Analysis

Subcritical californium source-driven noise analysis (CSDNA) measurements were performed at the Oak Ridge Critical Experiments Facility in 1983 using a stainless-steel tank with an inside diameter of 76.2 cm and a height of 91.4 cm but only filled with a uranyl (93.16 wt.% 235 U) nitrate solution to a height of 76.2 cm for all measurements. The free acid content of the uranyl nitrate solution was less than 0.01 N to minimize hazards in handling. The total available highly enriched uranium for these measurements was 5,296 g. The concentration of the solution was varied from the highest density of 14.71 grams of uranium per liter (g U/L) initially in 16 steps, lowering the concentration to 0.3492 g U/L, and finally using a tank filled with only water. The tank had a Plexiglas lid to minimize evaporation. The Cf source was contained in a re-entrant Lexan tube which could be in the solution parallel to the axis of the cylindrical tank. The detectors were 3 He proportional counters in the solution in shrink-fit tubing to isolate the counters from the fissile solution and adjacent to the outside of the tank. In some cases, the detectors were scintillators adjacent to the outside of the tank. In other cases, the source was external to the tank. Some data presented in this report are from notes and are not in the logbook. The prompt neutron decay constant was obtained from fitting the date from the CSDNA measurements. The purpose of this report is to document the experimental information for the measurements performed so that later, researchers could perform the required uncertainty and calculational analyses and documentation to use these data for an International Nuclear Criticality Safety Benchmark Program (ICSBEP) or a Nuclear Energy Agency benchmark. The prompt neutron decay measurements could be the basis of as International Reactor Physic Benchmark Program. The data from these measurements are available from the Records Management Services Department of Oak Ridge National Laboratory, and the logbook is also available from the ICSBEP at Idaho National Laboratory. Preparation of the present report is part of a larger cooperative effort between Idaho National Laboratory (INL) and Oak Ridge National Laboratory (ORNL) to document more than 15 undocumented critical and subcritical experiments enumerated in ORNL/TM-2019/18 and performed by ORNL at ORCEF and other USDOE critical experiments facilities using more than 500 operational days of critical facility time.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

FY22 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY22. As an extension of the effort from FY21, this year’s work focuses primarily on the two-phase modeling of the NSTF using RELAP5-3D, particularly with the inclusion of the cavity model. The results from simulations were used to compare against experimental data for benchmarking purposes of the RELAP5 deck. Additionally, RELAP5 was used as a predictive tool to guide planned test operations and identify expected system behaviors. In the first part of this report, details are provided of the cavity omitted model where heat flux is applied directly as a boundary condition to the risers. The general trend predicted by the RELAP5 model matches that from the experimental data when a single-phase natural circulation flow is first established, followed by an oscillatory two-phase period and finally a stable two-phase flow. However, the onset of oscillations is predicted early by the model due to the smaller thermal mixing region in the tank. However, by expanding the simulated thermal mixing region in the tank, the onset of oscillations predicted by the model is able to match that from the experiment. These oscillations where studied in depth and are deduced to be flashing-induced instability. The model was then modified to simulate an accident scenario case where a representative heat load based on the full-scale Framatome’s 625 MW t SC-HTGR was applied directly to the riser channels. The simulated initial and boundary conditions were identical to those performed experimentally, facilitating direct comparisons between the predicted and experimental results. It was determined that the results showed some discrepancies remain, likely due to the overprediction of vapor generation rate by the computer model. In the second part of this report, the cavity model is re-introduced where it is observed that the RELAP5 prediction is now able to capture the major trends of the observed flow commonly observed during two-phase conditions. However, the onset of oscillations is once again predicted early by the model, possibly caused by the underprediction of heat loss from the heater and cavity. This is likely due to the omission of support structures in the cavity that can act as additional pathways for heat to escape to the environment. To overcome the underprediction of heat loss, part of the insulation surrounding the cavity side panels and the back of the heaters are removed to allow heat to escape directly to the environment, which then improves the RELAP5 prediction. Parametric studies are also performed to investigate the effects of heater power, tank inventory level, and tank gas space pressure on flow behaviors, also in direct comparison to conditions tested experimentally. User option-61 in the RELAP5-3D input deck, which changes the heat transfer coefficient correlations used for calculating the vapor generation, is also investigated where it is found that by enabling the option, the overall duration of oscillations is increased and matches that from the experiment better. The RELAP5 model is further benchmarked with a header inlet- throttling case where it is observed that the prediction from the model fails to capture some major features observed in the experiment. By using a modified loss coefficient curve for the valve, the accuracy of the prediction is improved where most of the major features observed in the experiment are predicted by the model. Lastly, the model is benchmarked with an inventory depletion scenario where it is observed that despite the modeling limitation of RELAP5, the prediction shows good agreement with the experimental data where major trends and features are captured by the model. Future work will see continued development of the current RELAP5-3D input deck of the NSTF to both improve the accuracy of the model’s predictive capability and continuing serving the experimental program. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

How can a diverse set of integral and semi-integral measurements inform identification of discrepant nuclear data? [Slides]

Importance of nuclear data to bias changed significantly when including diverse measurement sets. LLNL pulsed sphere measurements provided information for 9 Be nuclear data above 2 MeV and subcritical benchmarks provided information for 240 Pu nuclear data between 100 keV and 10 MeV. Pulsed sphere leakage spectra and neutron noise observables are differently sensitive to nuclear data compared to critical benchmarks. Leakage spectra are sensitive to nuclear data above 5 MeV. Count rate and Feynman Y are more sensitive to nuclear data. A change in bias can help evaluators identify discrepant nuclear data. In the future, EUCLID plans to perform RAFIEKI analysis with additional benchmark and measurement sets.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SCALE 6.2.4 Validation: Nuclear Criticality Safety

The computational bias of criticality safety computer codes must be established through the validation of the codes to critical experiments. A large collection of suitable experiments has been vetted by the International Criticality Safety Benchmark Evaluation Project (ICSBEP) and made available in the International Handbook of Evaluated Criticality Safety Benchmark Experiments (ICSBEP Handbook). More than 600 cases from this handbook have been prepared and reviewed within the Verified, Archived Library of Inputs and Data (VALID), which is maintained by the Reactor and Nuclear Systems Division at Oak Ridge National Laboratory. The performance of the KENO V.a and KENO-VI Monte Carlo codes within the SCALE 6.2.4 code system is assessed using the VALID models of benchmark experiments. A range of nuclear cross section libraries based on Evaluated Nuclear Data File (ENDF)/B-VII.1 in both multigroup (MG) and continuous energy (CE) formats is considered. The critical experiments available to validate the KENO V.a code cover 15 broad categories of systems. These systems use a range of fissile materials, including a range of uranium enrichments, various plutonium isotopic vectors, and some mixed uranium/plutonium oxides. The physical forms of the fissile material also vary and are represented as metal, solutions, or arrays of rods or plates in a water moderator. The neutron energy spectra of the systems also vary and cover fast, intermediate, mixed, and thermal spectra. Over 550 of the total cases use the KENO V.a code for the four nuclear data libraries considered in this report.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Review of SCALE Validations Applicable to Spent Nuclear Fuel Shielding Calculations

This report presents a review of publicly available documents providing validations of SCALE code system capabilities for depletion and shielding calculations. The validation studies primarily used information available in the Spent Fuel Isotopic Composition Database, the Shielding Integral Benchmark Archive and Database, and the International Criticality Safety Benchmark Evaluation Project Handbook. Validation results based on radiochemical assay data and shielding benchmark experiments relevant to spent nuclear fuel shielding applications are summarized in this report. The information summarized in this report supports pressurized water reactor and boiling water reactor spent fuel storage and transportation safety analyses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗