Search NASA⌕ Search

SEARCH · Search NASA

Results for “Separate Effect Test”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Sodium-Cooled Fast Reactor Reference Plant Model

This report details the progress of Idaho National Laboratory (INL) in creating a reference plant multiphysics model for the Advanced Burner Test Reactor (ABTR). This model was developed under Task 13 of the U.S. Nuclear Regulatory Commission project “Development and Modeling Support for Advanced Non-Light Water Reactors,” and is an extension of the reference plant model developed in Task 4b, which was improved upon in the following ways. (1) The discrete ordinates method was used in lieu of the super-homogenization (SPH)-corrected diffusion approximation in order to better capture the anisotropic scattering contribution and the neutron leakage change due to thermal expansion. (2) The novel neutronic spatial discretization approach, termed the ring-heterogeneous (RH) approximation, was conceptualized and introduced to capture the differential expansion of the materials in the core. This new technique proved capable of preserving fission rates and maintaining the eigenvalue within 2.5% and 266 pcm with 9 neutron energy groups, respectively. Separating the different materials in the core enables the differential expansion of materials to be explicitly accounted for, eliminating the need for problem-specific cross-section functionalization techniques. (3) The SAM model for the core and system thermal-hydraulics analysis was updated to include 61 channels instead of just four representative ones. This enables users to obtain improved spatial resolution for sodium temperature and density scalar fields. (4) All the mesh files were created via the Multiphysics Object-Oriented Simulation Environment (MOOSE) Reactor module, eliminating all reliance on external tools for mesh creation. (5) Finally, the fuel axial expansion now leverages the HT9 and UPuZr material properties that have been validated against experimental data. The reference plant model was used to perform a full-core unprotected loss of flow (ULOF) transient calculation, including neutronics, thermal and mechanical feedback mechanisms. Future work will be devoted to further enhancements of the model. Potential improvements to the model include the addition of the control rod driveline expansion feedback and the upgrading of the support plate model so as to explicitly include 3D effects. Additionally, a Nuclear Energy Advanced Modeling and Simulation funded parallel effort has completely automated the creation of the ring-heterogeneous (RH) mesh from the fully heterogeneous (FH) geometry, thus maximizing user friendliness for the sodium fast reactor sodium-cooled fast reactor (SFR) workflow and will be incorporated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

The Use of Silver Chloride Injection in Remediation of Iodine-129 by In Situ Capture as Silver Iodide at the F-Area Seepage Basin - 20225

The Savannah River Site (SRS) produced tritium, plutonium, and special nuclear materials for national defense, medicine, and the space programs. As part of operations, the F-Area Seepage Basins operated until 1988 for the disposition of deionized acidic wastewater from the F Separations Facility. The wastewater contained dilute nitric acid and low concentrations of non-radioactive metals, and radionuclides, with the major isotopes being Cs-137, Sr-90, U-235, U-238, Pu-239, Tc-99, I-129, and tritium. The seepage basins were closed in 1988 and backfilled and capped by 1991. The groundwater emanating from beneath the closed and capped seepage basins is acidic and contains elevated levels of both chemical and radiological contaminants. Releases from the groundwater plumes sourced from the F-Area Seepage Basins have impacted the water quality of Fourmile Branch, which is a small tributary to the Savannah River, a regional water source. A large pump-and-treat system was constructed in 1997 and operated until 2003 in an attempt to capture the releases to Fourmile Branch. The system in F Area and a similar system in H Area were expensive (∼$1.3 M/month) to operate and produced large quantities of radioactive waste, with concentrations of I- 129 too high to be disposed of at the SRS. In 2004, SRS replaced pump-and-treat with a funnel and gate system that along with operation of a base injection system at the gates reduces the flux of contaminants to the wetlands adjacent to Fourmile Branch. The alkaline solution injected into the aquifer neutralizes the acidic plume and immobilizes many of the cationic constituents. However, base injection is not effective in managing the release of iodine-129, an anionic contaminant. To address iodine-129 SRS and the Savannah River National Laboratory developed an in situ technology that uses ultra-fine ground silver chloride (AgCl) as an injectable capture medium for the sequestration of iodine-129. The AgCl amendment has a very small particle size and is designed to be injected into the contaminated aquifer to capture iodine-129. Dissolved iodine-129 forms a stable and highly insoluble solid (silver iodide) upon contact with AgCl. Laboratory studies, a field scale pilot test (2009), and three deployments (2011, 2015 and 2019) of AgCl have been successfully performed at the F-Area Seepage Basins. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

COMPASS-FME Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) Experiment Level 1 Sensor Data v1-2

This is the version 1-2 Level 1 (L1) data release for COMPASS-FME environmental sensors located at our Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) experimental site. This manipulative, ecosystem-scale TEMPEST experiment addresses the potential for freshwater and estuarine-water disturbance events to alter tree function, species composition, and ecosystem processes in a deciduous coastal forest in MD, USA. The experiment uses a large-unit (2000 m2), un-replicated experimental design, with three 50 m × 40 m plots serving as control, freshwater, and estuarine-water treatments.L1 data are close to raw, but are units-transformed and have out-of-instrument-bounds and out-of-service flags added. Duplicates and missing data are removed but otherwise these data are not filtered, and have not been subject to any additional algorithmic or human QA/QC. Any scientific analyses of L1 data should be performed with care. **This dataset will be updated quarterly with new data for the duration of the project**This dataset includes:- An overall dataset README file that describes the current version, gives citation and contact information, etc.- Site- and year-specific folders, each holding up to 12 CSV (comma separated value) data files for each site and plot in that year.- Metadata files within each site-year folder provide full information on data units, expected ranges, contact information, detailed flood times, as well as a general description of the site.- Environmental sensor types that appear in the data files include weather (ClimaVUE50, CS, RM Young, and LI instruments in the graphs below); soil conditions (TEROS12); soil redox state (Redox); groundwater variables (AquaTROLL200 and AquaTROLL600); open water sondes (Exo); tree sap velocity (Sapflow); and system voltage and state (Datalogger). Data are normally logged every 15 minutes.Please see v1-2 TEMPEST L1 Sensor Package Quick Start.pdf for detailed information on data package structure, temporal coverage, and versioning.The TEMPEST flood events occurred on the following dates. They lasted for ~10 hours each day and delivered ~80,000 gallons to each plot; many data streams are available at 1 or 5 minute frequency during these periods.* Tests: Aug 25 (fresh plot) and Sep 9 (salt plot), 2021* TEMPEST 1: June 22, 2022* TEMPEST 2: June 6-7, 2023* TEMPEST 3: June 11-13, 2024

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

COMPASS-FME Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) Experiment Level 1 Sensor Data v2-1

This is the version 2-1 Level 1 (L1) data release for COMPASS-FME environmental sensors located at our Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) experimental site. This manipulative, ecosystem-scale TEMPEST experiment addresses the potential for freshwater and estuarine-water disturbance events to alter tree function, species composition, and ecosystem processes in a deciduous coastal forest in MD, USA. The experiment uses a large-unit (2000 m2), un-replicated experimental design, with three 50 m × 40 m plots serving as control, freshwater, and estuarine-water treatments. L1 data are close to raw, but are units-transformed and have out-of-instrument-bounds, out-of-service, and outlier flags added. Duplicates and missing data are removed but otherwise these data are not filtered, and have not been subject to any additional algorithmic or human QA/QC. Any scientific analyses of L1 data should be performed with care. **This dataset will be updated quarterly with new data for the duration of the project** This dataset includes: - An overall dataset README file that describes the current version, gives citation and contact information, etc. - Site- and year-specific folders, each holding variable-specific CSV (comma separated value) data files for each site and plot in that year. - Metadata files within each site-year folder provide full information on data units, expected ranges, contact information, detailed flood times, as well as a general description of the site. - Environmental sensor types that appear in the data files include weather (ClimaVUE50, CS, RM Young, and LI instruments in the graphs below); soil conditions (TEROS12); soil redox state (Redox); groundwater variables (AquaTROLL200 and AquaTROLL600); open water sondes (Exo); tree sap velocity (Sapflow); and system voltage and state (Datalogger). Data are normally logged every 15 minutes. Please see v2-1 TEMPEST L1 Sensor Package Quick Start.pdf for detailed information on data package structure, temporal coverage, and versioning. The TEMPEST flood events occurred on the following dates. They lasted for ~10 hours each day and delivered ~80,000 gallons to each plot; many data streams are available at 1 or 5 minute frequency during these periods. * Tests: Aug 25 (fresh plot) and Sep 9 (salt plot), 2021 * TEMPEST 1: June 22, 2022 * TEMPEST 2: June 6-7, 2023 * TEMPEST 3: June 11-13, 2024 This dataset was updated 2026-03-12: (i) data now go through 2025-12-31 (previous end was 2025-06-30) and (ii) dataset and file names updated to “…v2-1” (previously was “v2-0”).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Production of Carbon Nanomaterials and Sorbents from Domestic U.S. Coal (Final Report)

The main goal of this project was to produce high-value carbon nanomaterials and carbon sorbents from domestic coal resources in a cost-effective manner. Four types of domestic coal samples were processed through a combination of deashing, devolatilization, oxidation, reduction, and activation treatments to produce graphene oxide (GO), reduced graphene oxide (RGO), and activated carbon (AC). The precursors and developed materials were extensively characterized by various methods to investigate the impact of the coal feedstock type on the yield and quality of each product. A commercial graphite-based GO sample was included in the experimental work as the baseline material for comparison with the coal-based materials developed in this work. The performed work also included a technoeconomic analysis and cost estimation for a plant processing 20 tons/day coal, a market evaluation for the graphene materials, and a technology gap analysis. A simple process by concentrated nitric acid oxidation is used to oxidize coal precursors for the production of GO. Fine oxidized coal particles that were separated from larger oxidized coal particles had significantly higher oxygen contents and were identified as coal-based GO samples. Coarse oxidized coal particles were used as precursors for production of AC. Based on the Raman spectroscopy results, coal-based GO samples exhibited G and D bands similar to those of graphite-based GO samples. X-ray photoelectron spectroscopy revealed that coal-based GO samples had surface oxygen contents of ~26-35% that were higher than the oxygen contents of graphite-based GO samples. Larger particles of oxidized coal samples were activated under different conditions to produce high surface area functionalized AC. Prepared materials had surface areas exceeding 1,500 m 2 /g and pore volumes more than 1 cm 3 /g with different pore size distributions. Reduced graphene oxide samples were prepared by thermal reduction of both coal-based and graphite-based GO samples at 170-2800 ºC. Coal-based and graphite-based RGO samples exhibited similar carbon contents, Raman spectra, and XRD profiles. Heat treatment above 1500 ºC shifted RGO to synthetic graphite. Among anthracite, bituminous, subbituminous, and lignite coals tested, anthracite was the best precursor to produce carbon nanomaterials exhibiting properties similar to those of graphite-based materials. Anthracite-based carbon nanomaterials also had the highest production yields. Technoeconomic analysis estimated the production cost of GO and RGO for anthracite-based samples at about 2,600 and 4,200 $/ton, respectively, which is about two orders of magnitude lower than the current estimated price of graphite-based materials. Several gaps to further develop the proposed technology were identified and discussed that include process and equipment optimization, process integration, need for additional bench- and pilot-scale experiments, and other items to reduce the scale up risk. Market analysis reports suggested a compound annual growth rate of 40% for graphene and related materials. Short-term applications include composites, inks, and coatings. However, energy storage applications appear to be the dominant potential long-term applications. Different applications of coal-based GO and RGO need to be explored and clear metrics and standards for each application need to be developed.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Update on Full-Scale Retrieval and Transfer System Testing for Pneumatic Conveyance of High-Level Calcine Waste - 20241

The Idaho National Laboratory processed spent nuclear fuel from 1953 to 1994. The process used solvent extraction to separate uranium, resulting in a liquid waste. A calcination process was used to convert the liquid waste to a high-level radioactive solid powder, referred to as calcine. The calcine is stored in stainless steel bins contained within concrete vaults called the Calcined Solids Storage Facilities. There are a total of seven CSSFs. The Calcine Retrieval Project is tasked with retrieving the calcine from the bins of CSSF 1 and pneumatically transferring it to the bins of CSSF 6. CSSF 1 was constructed to be a final calcine repository with no forethought of retrieval. It represents the greatest technical hurdle to calcine extraction of all the CSSFs. The challenges inherent with accessing the calcine, safely extracting it, and securely transporting it are copious and complex. Once the calcine is removed from CSSF1, that facility may be closed, reducing the footprint of waste storage at the INL. The same process for removing Calcine from CSSF1 may then be utilized to remove the calcine from the remaining CSSFs for final treatment and transport to a long-term storage facility. To overcome these inherent challenges with accessing CSSF 1 and extracting the calcine in a timely manner, a full-scale integrated mockup has been constructed. The mockup includes a full-size model of a concentric group of bins from CSSF 1 and a storage vessel representing CSSF 6. All essential elements of the transport design are represented in the Mockup, including motive and makeup air supply (air compressor), a replica of the CSSF6 cyclone, pre-filter, HEPA filter, vacuum pump, and retrieval control devices. Construction on the Mockup was completed in 2018 and testing began in February 2019. It has already proven very beneficial to the CRP. Experience in running the Mockup have aided in the development of retrieval and transport methodologies and the selection of process parameter values. It has also uncovered some unforeseen technological challenges and will continue to be valuable to the testing team as these hurdles are overcome and the process is proven effective. Operations personnel have also benefitted from its use as a training tool and will continue to use the Mockup as a simulator to reproduce problems encountered during extraction and transfer of Calcine from CSSF 1. The knowledge gained during this process will also be of continued benefit as future work is done to remove various types of calcine from the remaining CSSFs. The methodology of safe and secure transfer of radioactive material may be of benefit to others throughout the U.S. Department of Energy complex. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Active Aerodynamic Load Control for Wind Turbines

The goal of this project was to develop and demonstrate an advanced dielectric barrier discharge (DBD) plasma actuator technology. We set out to demonstrate the efficacy and impact of the new actuator technology as a key component of an active load control system for wind turbines. A DBD plasma actuator consists of a thin layer of dielectric material separating a pair of offset electrodes. One of those electrodes is embedded between the dielectric layer and a non-conductive substrate, while the other is exposed to air. When driven by an appropriate high voltage waveform, the device ionizes the air near the surface of the dielectric and adjacent to the exposed electrode. Collisions between the ions in the plasma and neutral air molecules result in a wall-jet – a region of induced air velocity that can be used to modify the flow around a lifting body. Use of the device near the trailing edge of a wind turbine blade, designed in such a way as to amplify the effects of the flow-modifying device, can result in large changes to the global forces experienced by the blade. Because of the fast response time of the device, it can allow the turbine to react in real time to changes in the wind associated with turbulence, wind shear, gusts, and the like, when paired with appropriate sensors and control algorithms. The goal of this project was to increase the capacity of the device to induce velocity on its surface at the levels required by large, utility-scale wind turbines. The first technical goal was to increase the induced velocity from the current industry-best of about 4 m/s to 10 m/s by modifying the electrical waveform used to drive the device and by introducing a semi-conductive surface coating to control electrostatic charge build-up on the surface. The second goal was to use the device to modify the lift on a representative airfoil in a wind tunnel, demonstrating a reduction in lift coefficient of 0.2 or better. The third goal was to produce an actuator-induced velocity of 20 m/s. Alongside our partners at the University of Texas at Dallas, we also applied modern, advanced design methods to optimize the impact of the technology on the design of wind turbines. We also had planned to install a segmented, active load control system on a test turbine to demonstrate the ability to reduce unsteady aerodynamic forces on the turbine associated with changes in the wind. The highest induced velocity achieved during this research was 11 m/s. However, practical design constraints limited the change in lift coefficient to about 0.12 for a representative airfoil in the wind tunnel at a Reynolds number of 400,000. The primary conclusion is that the plasma actuator control authority remains insufficient for practical purposes when extrapolated to Reynolds numbers over 1 million. The potential impact of the active lift control concept was evaluated through detailed simulations for three wind turbine sizes: a 3.4MW onshore turbine, a 10MW offshore turbine, and a 15MW offshore turbine. A feedback control system was designed for each turbine within two scenarios: one where the active lift control is used as a retrofit capability on the baseline design, and the other where the designers were permitted to “upscale” the turbines in order increase annual energy production. The levelized cost of energy was then evaluated for the range of turbine sizes and design configuration. It was found that the LCOE reduction associated with active lift control fell in the range of 0.7% to 7.2%, with the highest reduction associated with upscaling the 3.4MW turbine.

17 WIND ENERGY↗
DOE OSTI2020

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Neutron irradiation & thermomechanical experiment (NITE) - design

For the reliable long-term operation of fusion power plants, it is crucial to understand and predict the lifetime of materials in use. These materials include all structural and functional materials utilized at the first wall, blanket, magnets, and shielding. The key challenge is, that the harsh environment including high heat fluxes, high thermal stress and stress cycling, neutron irradiation, and sputtering on such materials should not be viewed separately. Currently, the synergistic loads cannot be evaluated experimentally because of the lack of adequate facilities. The purpose of that work is to design a synergetic Neutron Irradiation and Thermomechanical Experiment (NITE) for fusion materials. This design will leverage the existing Advanced-Test-Reactor (ATR), a fission reactor at the Idaho National Laboratory. We also acknowledge that with existing fission reactors the exact fusion condition cannot be created, and the limitations are critically discussed. The combination of neutron irradiation with a high heat flux is the focus. This is realized with an irradiation capsule design that includes a TRISO fueled region inside the capsule to enable a steady-state heat flux on one side of the specimen. In conclusion, the experimental design modeling showed that steady-state heat fluxes of 2.4 MW/m 2 with a thermal gradient of above 250°C can be achieved in a 5 mm thick specimen.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Physics Division Strategic Plan Fiscal Years 2020-2024

The vision of the Physics Division (PHY) at Argonne National Laboratory is to continue enhancing its role as a world-leading institution in basic nuclear physics research and its applications. Key to this vision is for PHY to continue to safely and effectively operate and evolve the capabilities of the Argonne Tandem Linac Accelerator System (ATLAS) facility to best serve its users. ATLAS is the Department of Energy (DOE) accelerator facility for low-energy nuclear physics research. The research carried out by PHY covers many themes in contemporary science but can be distilled into five areas of focus. These themes are interconnected, with continuously evolving synergies between the various groups and facilities in PHY and the broader Laboratory. This five-year strategic plan serves to illustrate our current capabilities and new directions related to these five areas of focus (ATLAS also develops an independent strategic plan): Accelerator research and design. The goal of this theme is to design, fabricate, test, and implement novel accelerator systems, with a focus on high-intensity ion and electron systems. These R&D activities have led to enhancements in the ATLAS accelerator system. Among others, the division’s accelerator systems are in use or planned for use at Fermilab, the Advanced Photon Source, and in a broad variety of applications at the Facility for Rare Isotope Beams (FRIB). Atom trapping and fundamental symmetries. The goal of this theme is to explore and exploit the uses of advanced laser cooling and trapping techniques to manipulate atoms. There are three main areas of focus. First is in the application of the atom trap trace analysis (ATTA) technique for age determination of groundwater and ice by using radio-krypton dating. The next two trapping-based programs involve tests of fundamental symmetries in nature: the cooling and trapping of radium-225 with the aim of determining limits on an observation of its electric dipole moment and precision measurements of the beta decay properties of helium-6 to set limits on the tensor coupling constant. This is complemented by measurements of similar properties in lithium-8 and boron-8. Nuclear astrophysics. The goal of this theme is to enhance our understanding of how elements are created in the universe via explosive nucleosynthesis and how stars evolve. To meet the challenges of this theme, many of the capabilities of ATLAS are being enhanced, including the development of new beams through a new in-flight separator (RAISOR) and the anticipated neutron-generator upgrade of the Californium Rare Isotope Breeder Upgrade facility (nuCARIBU). These are coupled to state-of-the-art instruments such as the Canadian Penning Trap, Helical Orbit Spectrometer, Gammasphere, GRETINA, Multi-Sampling Ionization Chamber, the Fragment Mass Analyzer, and Argonne Gas-filled Fragment Analyzer, and a new low-background experimental area for decay studies. Nuclear structure. The goal of this theme is to understand the structure of nuclei, both stable and radioactive, in terms of single-particle properties, their shapes, and the dynamics governing reactions between them. These include questions such as what are the limits of nuclear stability and what are the properties of super heavy nuclei. As with the nuclear astrophysics theme, the capabilities of ATLAS, guided by the ATLAS user community, are continuously being enhanced to this end. The Division has strategic initiatives to play a leading role in the development of instrumentation and research programs at the FRIB. Quantum chromodynamics (QCD) and hadron physics. The goal of this theme is to lead major research programs focused on revealing the quark and gluon structure of protons, neutrons, nuclei, and short-lived mesons and baryons. These involve major programs at Jefferson Lab, Fermilab, and smaller facilities. They are complemented by theoretical endeavors centered on the question of how hadrons and their properties emerge from QCD. Out of these activities arise strategic initiatives to play a major role in the forthcoming Electron-Ion Collider (EIC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Acquisition and Implementation of a Comprehensive Environmental Permits Linking Tool at Savannah River Site - 20234

Historically, Savannah River Nuclear Solutions LLC (SRNS) tracked environmental regulatory commitments and the requirements from hundreds of permits at the Savannah River Site (SRS) using several separate methods, making integrated compliance assurance cumbersome and labor-intensive. When SRNS experienced an increase in environmental issues in 2017, SRNS management and U.S. Department of Energy - Savannah River (DOE-SR) management decided a single, proactive approach was needed to capture environmental permit information (including regulations, Consent Orders, DOE Orders, and any other state or federally issued statement of requirements), track the tasks necessary to ensure compliance with these requirements, and thereby mitigate the risk of noncompliance. SRNS developed a list of mandatory objectives that the tool must meet to function as a Comprehensive Environmental Permits Linking Tool (CEPLT). A key requirement was the ability to map Site permits to their governed locations and display the associated requirements at the compliance point (e.g., outfall, stack, waste unit, etc.). Several options included modifying existing onsite resources, building a custom onsite solution, purchasing an off-the-shelf solution, and contracting an offsite developer to build a custom solution. SRNS concluded an off-the-shelf solution with configuration and customization options would provide the flexibility to fit the unique needs of Savannah River Site (SRS) while taking advantage of industry-tested software and providing a reduced deployment timeline. SRNS chose Gensuite{sup R} a, a cloud-based solution that offers numerous a la carte applications in the environmental, health, and safety arenas, as the best candidate. SRNS selected three (3) integrated applications (Compliance Calendar, Permit Manager, and Mapper) to function as the CEPLT. The Compliance Calendar module allows for creation and tracking of regulatory commitment tasks assigned to responsible environmental professionals. Permit Manager organizes permits and other requirement documents, linking the commitments in each to Compliance Calendar tasks and/or implementing procedures. Mapper provides GIS capability for mapping the data from the other two modules to their physical onsite locations. The CEPLT was configured to allow for other SRS Site Tenants to eventually utilize the applications. Additionally, DOE-SR uses the CEPLT to provide an overview of contractor environmental compliance activities and to organize DOE-specific documents and tasks. The CEPLT fulfilled the requirement of meeting current compliance needs as well as providing the ability to grow as new organizations are incorporated and functionality is expanded. SRNS now uses the CEPLT to more effectively manage permit requirements, improve knowledge transfer, and increase the Site's overall protection of the environment, the Site worker, and the public. Possible future uses of the system include integration of mobile applications for timely communication of potential non-compliant conditions and DOE complex-wide deployment allowing for enhanced DOE site and Head Quarters oversight. Implementation of CEPLT will result in cost savings, both in terms of dollars and man-hours. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Data from: "Warming and the dependence of limber pine (Pinus flexilis) establishment on summer soil moisture within and above its current elevation range"

This data package contains data that were used for analysis in “Warming and the dependence of limber pine (Pinus flexilis) establishment on summer soil moisture within and above its current elevation range”, by Moyes et al. 2013. All data collection and field research were completed on Niwot Ridge, Colorado, USA.This data package contains nine comma-separated-values (.csv) files, one text (.txt) file, and two zipped seedling folders that were used for leaf area analysis. One zipped folder contains 468 .jpg photographs of seedlings, and the second contains 468 corresponding Image J-processed .jpg images that include silhouette leaf area values. .csv and .txt files can be opened using any compatible simple text-editor software such as TextEdit (Mac) and Notepad (Windows); .csv’s can also be opened using R and Microsoft Excel. Image files can be opened using Preview (Mac) and Photos (Windows). In addition, there are a total of 31 Microsoft Excel files: three .xlsx files, and 28 raw Li-Cor output .xls files. This data user’s guide is available in .pdf format, and can be opened using Adobe Acrobat Reader, or any other compatible file viewing software. Geospatial data showing field site locations are also included in the archive for use and reference. There are two geospatial formats in this archive: ESRI shapefiles (.shp) and keyhole markup-language (.kml) files. Both file types contain bounding box information, with the former being polygons, and the latter containing corner coordinates for each site. ESRI shapefiles can be opened using any geospatial software compatible with the file type (such as ESRI’s ArcGIS suite and QGIS), and .kml files are compatible with Google Earth and Google Maps.--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------Continued changes in climate are projected to alter the geographic distributions of plant species, in part by affecting where individuals can establish from seed. We tested the hypothesis that warming promotes uphill redistribution of subalpine tree populations by reducing cold limitation at high elevation and enhancing drought stress at low elevation. We seeded limber pine (Pinus flexilis) into plots with combinations of infrared heating and water addition treatments, at sites positioned in lower subalpine forest, the treeline ecotone, and alpine tundra. In 2010, first-year seedlings were assessed for physiological performance and survival over the snow-free growing season. Seedlings emerged in midsummer, about 5–8 weeks after snowmelt. Low temperature was not observed to limit seedling photosynthesis or respiration between emergence and October, and thus experimental warming did not appear to reduce cold limitation at high elevation. Instead, gas exchange and water potential from all sites indicated a prevailing effect of summer moisture stress on photosynthesis and carbon balance. Infrared heaters raised soil growing degree days (base 5 °C, p < 0.001) and August–September mean soil temperature (p < 0.001). Despite marked differences in vegetation cover and meteorological conditions across sites, volumetric soil moisture content (θ) at 5–10 cm below 0.16 and 0.08 m^3 m^(−3) consistently corresponded with moderate and severe indications of drought stress in midday stem water potential, stomatal conductance, photosynthesis, and respiration. Seedling survival was greater in watered plots than in heated plots (p = 0.01), and negatively related to soil growing degree days and duration of exposure to θ < 0.08 m^3 m^(−3) in a stepwise linear regression model (p < 0.0001). We concluded that seasonal moisture stress and high soil surface temperature imposed a strong limitation to limber pine seedling establishment across a broad elevation gradient, including at treeline, and that these limitations are likely to be enhanced by further climate warming.--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------Maintenance log:Amended file name error in Data User's Guide May 19 2022

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗
DOE OSTI2020

Trace Chemical Detection Using Intercalated MXenes as a Signal Enhancing Substrate in Optical Probes

MXenes are 2D materials composed of layered transition metal nitrides or carbides. These materials are synthesized by HF exfoliation from MAX phases (Ti{sub 3}AlC{sub 2}). The 2D nanomaterial was synthesized by the removal of the 'A' element, resulting in a Mxene product (Ti{sub 3}C{sub 2}). MXenes have the general formula M{sub n+1}X{sub n}T{sub x}, where M is an early transition metal, X is Carbon and/or Nitrogen, and T accounts for surface terminated functional groups such as Fluoride, hydroxyl, and oxygen. These materials have very unique properties, similar to graphene, that allows them to be applied in a variety of trace detection techniques including surface-enhanced Raman spectroscopy (SERS). MXenes have also been demonstrated to selectively uptake uranyl ion, UO{sub 2}{sup 2+}. If this property can be combined with SERS or fluorescence detection, it may be possible to use MXenes as the basis for an alternative method to kinetic phosphorescence analysis (KPA) for trace uranyl measurements. Objectives: To confirm that MXene Nano materials are suitable substrates for SERS and sensor development by enhancing Raman signaling. To determine if certain MXene preparation methods yield materials that are more suitable for trace sensing methods. To determine uranyl uptake properties of these MXene materials and test them for analytical signals. Sample Preparation: Preparation of Ti{sub 3}C{sub 2}MXene (at FSU). MXenes were prepared by etching Al from Ti{sub 3}AlC{sub 2} (MAX phase)material. Two etching techniques yield different MXene products: LiF/HCl: Milder reaction, larger MXene flakes. HF: Harsher reaction, smaller flakes, larger layer separation. Products washed to remove etchant, vacuum filtered, and dried. Dried MXene flakes are air-stable. Film preparation for sensor testing (at SRNL): Suspend powder in diH{sub 2}O, purge with Ar, sonicate for 30 min. Centrifuge and collect supernate with suspended particles. Observed LiF-etched Mxene yielded a higher density of particles and darker collected solution. Drop-cast (4 ml) supernate onto slides and dried with Ar. For Rhodamine B (RhB) testing, drop-cast 4 ml drops onto Mxene spots and dried with Ar. Scanning Electron Microscopy conditions: 10 kV Beam energy, high vacuum; Working distance of 8 mm; beam penetration depth appx. 4 microns, beam spot size appx. 2 nanometers. Results: Detection of aluminum correlates with bright spots on image. Presence of aluminum shows that LiF/HCl etching was less thorough than HF etching. Trace Cl detection in LiF images suggests incomplete rinsing. HF has smaller feature size, more layer structure, and increased homogeneity, consistent with expectations. Macroscopic Raman spectroscopy measurements: 532 nm excitation, ∼50 mW with a ∼100 micron spot size (InPhotonics RPB probe). Kaiser Optical Holospec f/1.8 spectrometer with cooled (-60 deg.C) Andor iDus OE420 CCD. LiF 1x supernate showed good signal for trace measurements of Rhodamine B. HF and 1/4x LiF supernates showed little Mxene or Rhodamine B signal. Low deposition densities led to excess background signal from glass slides. For LiF film, response is linear with Rhodamine B concentration over range tested. Will retest with Raman microscope (∼1 micron spot size) to characterize SERS of more dilute LiF and HF etched Mxenes. Conclusions: The LiF etched material was more suitable for macroscopic SERS measurements because it was more concentrated, resulting in a thicker film than the HF etched Mxene and diluted LiF sample. However, the other materials may give greater SERS enhancements, which we hope to determine from measurements with the Raman microscope. From characterization with SEM we concluded that the HF etched Mxene is more uniform/homogenous and has smaller particle size than the LiF etched Mxene. There is still aluminum present in both samples indicating that etching wasn't complete, but the removal of the aluminum was more efficient in the HF method. Path Forward: Observe SERS with Raman microscopy, to obtain better signals for the more diluted samples and be able to compare enhancement effects for the different MXenes. Characterize uranyl sorption into MXene films and test Raman and fluorescence signals. Revisit the etching conditions to improve removal of aluminum. FSU and SRNL will continue to collaborate to create and characterize different Mxene materials and test their usefulness for sensor applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗
DOE OSTI2024

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Initial Assessment of TeF 6 Adsorption

Research and development that supports the management of off-gases from nuclear fuel reprocessing has historically been focused on the off-gas streams that arise from aqueous reprocessing technology. However, as Gen-IV reactor development pathways move toward deployment, alternative spent nuclear fuel (SNF) processing and disposition pathways have been considered more actively. This work is focused upon aspects of fluoride volatility (FV) processing. The versatility of this method for deployment against multiple types of spent fuel encourages the continued advancement of both the primary separations processes and the secondary processes, including waste treatment, material control and accountability, and engineering designs. Recent work noted that TeF 6 , the most highly volatile fluorinated compound produced during FV processing, did not have a clear abatement technology recommended in the literature. Thus, this work performed scoping tests on activated alumina and copper shot to assess whether they could be used to remove TeF 6 from gas streams that bear F 2 . Previous work on this topic was not well described in the literature and was not performed with excess F 2 in the stream as would be typical of spent fuel processing via FV. To support an understanding of the concentration of TeF 6 contacting the adsorbent beds, a series of preliminary testing identified the TeF 6 production rate and the equilibrium concentration of TeF 6 in the gas stream contacting the adsorbent. Excess F 2 was determined to not affect the ability of activated alumina to remove TeF 6 from the gas stream quickly and completely. A determination of whether excess F 2 affected the distribution depth of Te in the sorbent bed is still pending analysis of the used sorbent. Literature suggests that when activated alumina is near saturation, TeF 6 could migrate from the sorbent bed. Future testing should investigate this possibility. Copper metal did not adsorb TeF 6 in the presence of F 2 across the sorbent temperature range of 50 to 335°C. F 2 was fully removed by the copper bed. Although preliminary thermodynamics would indicate adsorption to be energetically favorable, other factors that impact adsorption (e.g., slow kinetics, excess fluorine on the copper surface, an unfavorable transition state) are likely preventing the adsorption of TeF 6 at an easily measurable rate. The work also investigated TeF 6 production from multiple forms of Te in the temperature range of 100 to 250°C. No previous study had assessed the initial reaction rates for this process. A carefully designed study allowed determination of the activation energy for the production of TeF 6 from Te metal. The substantial amount of data collected during this study merits analysis beyond what is described here. A more in-depth kinetic analysis will be pursued. Additional analytical results will provide the ability to benchmark adsorption coefficients for TeF 6 , understand the distribution of TeF 6 within the alumina bed, and better understand the effect of F 2 partial pressure on Te fluorination. The data from this report, as supplemented by these additional analyses, will be submitted to a peer-reviewed journal.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗