Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Accident Analysis”

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 505 records · Page 28

Radiation Shielding Analysis for the PIP-II Linac at Fermilab

The Proton Improvement Plan-II (PIP-II) [1] has been developed at Fermilab to provide powerful proton beams to the laboratory’s experiments. An 800-MeV superconducting linear accelerator—a centerpiece of the project—is currently under construction in Batavia, Illinois (USA). After completion, the superconducting linac will be the starting point for the 1.2 MW (Phase 1) and 2.4-MW (Phase 2) proton beam that is needed for the Long-Baseline Neutrino Facility (LBNF) at Fermilab [2]. Due to unavoidable loss of a fraction of the beam in the accelerator components, a certain level of radiation will be generated in the accelerator tunnel both during normal operation and at accidents. This work deals with radiation shielding design for the accelerator facility.

Rakhno, Igor↗

Carbon Dioxide Utilization for Plasma Nanosynthesis of Carbon

This DOE Phase I SBIR project has developed a glow discharge-based plasma process to convert carbon dioxide to sold carbon and water. By adding natural gas, the additional enthalpy can be used to reform the output carbon dioxide to value added carbon and EPA-compliance-grade water, thus creating a process that has no gaseous products at all, such as: CO 2 + CH 4 → 2C + 2H 2 O, which is simply the exothermic reverse of the syngas reaction. During the Phase I work, several milestones were accomplished successfully. a. First, a synthesis reactor was designed, assembled, and successfully operated using the combined simultaneous operation of glow discharge and microwaves. DC discharge and microwaves can also be operated individually, allowing multiple configurations to be used This is one of the first reactors of its kind, and perhaps the first ever. b. Second, the research team successfully produced carbon nanomaterials using carbon dioxide as the source of carbon via CO 2 + 2H 2 → C + 2H 2 O. This is certainly one of the first processes to use plasma to convert carbon dioxide to nanocarbon, and perhaps the first ever demonstration. c. Third, it was also shown that carbon monoxide and carbon dioxide can be used with methane to reduce the synthesis temperature of ASI graphene decorated carbon nanotubes (GDCNT). Hitherto GDCNT were only fabricated above 1200 °C. GDCNT was synthesized at 780 °C. CO 2 + C → 2CO; 2 CO + CH 4 → 3 C + 2H 2 O, and also, CH 4 + CO 2 → 2 C + 2H 2 O. d. Fourth, closed cycle synthesis of nanomaterial was demonstrated with over 80% conversion of carbon dioxide to carbon and water was achieved, on the basis of pressure drop during glow discharge. e. Fifth, plasma operation was achieved via glow discharge and microwaves at 500 Torr, demonstrating a path forward to atmospheric pressure operation f. Sixth, using carbon monoxide and carbon dioxide in an electrothermal furnace, we produced 350 grams per hour of high-quality nanotubes in a single furnace, compared to 200 grams per hour for the commercial process. Controlling the enthalpy prevents self-extinguishing hot gas pyrolysis. Thus, carbon dioxide enhances the properties and reduces the production cost. This is a near term advance that can go into trial production in the near term. g. Seventh, economic analysis suggests an economic plan forward, producing carbon nanotubes for less than $50 per pound. The use of hydrogen and carbon dioxide as the feedstocks is acceptable from a cost basis, and in fact hydrogen is less expensive than carbon monoxide. Niche applications are available in the near term, though quantities relevant to sequestration remain a far-term goal. h. Eighth, though it was not the research team’s intention to create a new process for producing carbon monoxide, it was realized that almost by accident we demonstrated plasma reforming of carbon dioxide via a plasma version of the reverse Boudouard reaction. Accordingly, we have contacted a major gas manufacturer (Matheson Gas) to determine whether there is interest in creating an industrial version of the process. i. Ninth, a long-range strategy is identified in which a progressively more carbon-conscious world will seek to create solid carbon as a means of sequestering carbon dioxide, impeding its return to the atmosphere. In the long term the price must drop by orders of magnitude with larger production, and assuming success, solid carbon could be the preferred form for sequestration rather than underground storage. If carbon can be made as cheaply as dirt, or nearly so, then it can serve a useful purpose as structural material for housing, soil amendments in agriculture and fillers to modify low areas.

03 NATURAL GAS↗

Coupling Shared E-scooters and Public Transit: A Spatial and Temporal Analysis

The integration of shared e-scooters with public transit is a promising solution for urban mobility's first/last-mile challenge. This study explores spatiotemporal factors influencing this integration, using 35-day e-scooter trip data from Chicago. Employing a random-effect negative binomial approach, we modeled the frequency of e-scooter trips to access/egress to/from bus stops and train stations. Results indicate that weather conditions, design features like intersection density, and multimodal network density significantly influence usage. The transit system characteristics such as service frequency have a positive effect on the integration of e-scooters and trains while a similar effect for bus and e-scooter integration was not significant. Furthermore, safety-related variables such as accident and crime rates as well as demographic characteristics were also revealed to be significant factors in our study. These findings offer vital insights to urban planners and policymakers for infrastructure, safety enhancements, and interventions to encourage efficient e-scooter-public transit integration.

Chicago↗

Characterization of Fuel-to-Coolant Heat Transfer During Reactivity-Initiated Accidents Using Tightly Coupled Thermal Hydraulics and Fuel Thermomechanics

The reactivity-initiated accident (RIA) is a complex scenario with several tightly interacting physical phenomena. Accurately predicting fuel behavior during these transients is difficult due to limitations in the modeling of fuel-to-coolant heat transfer. Common approaches to simulate RIAs involve standalone calculations using either a fuel performance code or a thermal-hydraulic code. The complex interdependencies of thermal-hydraulic and fuel mechanical behavior suggest that a tight coupling between these codes may provide more accurate predictions of fuel-to-coolant heat transfer and cladding mechanical response. Here, RELAP5-3D and BISON are coupled in this paper to simulate RIAs, and a sensitivity analysis is performed to rank key thermal properties and two-phase heat transfer parameters relevant for fuel-to-coolant heat transfer and cladding failure mechanisms in UO 2 –Zircaloy-4 systems. Gas gap conductance, film boiling heat transfer uncertainty, pulse width, fuel-specific heat capacity, and cladding-specific heat capacity were identified as important parameters. Variations in figures of merit resulting from changes to pulse width and the material thermal properties indicate that time-dependent heat transfer rates are significant for safety-relevant mechanical parameters due to the time dependence of cladding ductility and pellet-cladding mechanical interaction loading. The results suggest that the thermal-hydraulic factors have a nonnegligible influence on the thermomechanical solution and vice versa. Tight coupling of both sets of physics is recommended to improve prediction of fuel behavior during RIAs. Highlights include the following: 1. The RELAP5-3D thermal-hydraulic code and the BISON fuel performance code are tightly coupled for simulation of RIA transients with energy depositions at the Zircaloy-4 cladding failure threshold. 2. Departure from nucleate boiling occurred for all simulated cases. Due to the ductility of fresh fuel, substantial ballooning occurred in most cases. 3. Gas gap conductance, fuel-specific heat capacity, cladding-specific heat capacity, transient pulse width, and film boiling heat transfer were the dominant thermal factors impacting the safety figures of merit at energy depositions.

Critical Heat Flux (CHF)↗

Diagnosing nuclear power plant pipe wall thinning due to flow accelerated corrosion using a passive, thermal non-destructive evaluation method: Feasibility assessment via numerical experiments

Flow accelerated corrosion (FAC) in nuclear power plant pipes is one of the leading causes of accidents, fatalities, damage and outages. Current FAC identification methods employ expensive sensing technology and are “active” methods, where the response of the piping system to an externally-generated thermal, mechanical or optical excitation must be measured. As a result, these techniques require a disruptive and time-consuming setup. Here we propose a method that utilizes pipe surface temperature measurements to passively monitor for FAC-induced pipe wall thinning without the need for expensive equipment or post-installation setup time. This diagnostic method utilizes a simulation data-driven diagnostic model to estimate the amount of thickness reduction in a pipe based on changes in measured steady-state pipe temperatures. In order to reduce the computational burden of generating large, simulation-based datasets, the behavior of the insulation of the pipe was modeled using a suitably calibrated heat transfer boundary parameter. Additionally, global sensitivity analysis was performed to determine system parameter(s), such as the temperature of water flowing inside the pipe, which significantly affect the steady state pipe wall temperature and could cause errors in diagnosis. Two diagnostic models, one using only the change in steady-state temperature as an indicator for FAC-induced pipe wall thinning and the other using water temperature as an additional diagnostic model input were evaluated for their ability to estimate thickness reductions in a pipe using simulated pipe wall temperature data. For the numerical experiments conducted in this work, both models estimated wall thickness with errors within 0.5 mm, indicating that the proposed technique can potentially be used as a low-cost, first-pass method for FAC monitoring.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Expansion of Hazards and Probabilistic Risk Assessments of a Light-Water Reactor Coupled with Electrolysis Hydrogen Production Plants

This report builds upon the body of work sponsored by the Department of Energy (DOE) Light-Water Reactor Sustainability (LWRS) Flexible Power Operation and Generation (FPOG) program that presented generic probabilistic risk assessments (PRAs) for the addition of a heat extraction system (HES) to light-water reactors to support the co-location of a high temperature hydrogen electrolysis facility (HTEF). Probabilistic and deterministic hazards assessments and risk analyses are leveraged throughout this report. Several improvements and new analyses are included in this report. First, higher amounts of detail in the specifications of the generic HTEFs are used to produce scaled results for a 100, 500, and 1000 MW nominal hydrogen production facility. An additional hazard assessment of 1000 kg of hydrogen storage is performed. The facility hazards and footprint are assessed to determine the safe distance required for placement near the nuclear power plant (NPP). Second, specific designs for corresponding HESs for the different levels of support required by the HTEFs are analyzed in the PRA model. Third, a hazards analysis of the specified HTEFs leads not only to effects of the quantified risk assessment for the NPP, but also qualitative hazards assessment for the community. Finally, a seismic analysis and a high winds analysis have each been added to the PRA. The results investigate the applicability of the potential licensing approaches which do not require a full United States (U.S.) Nuclear Regulatory Commission (NRC) licensing review. The PRAs are generic and include listed assumptions. The HTEF design built for this project has further eliminated many conservative assumptions from the prior PRAs in this series. The PRA results indicate that the 10 CFR 50.59 licensing approach is justified due to the minimal increase in initiating event frequencies for all design basis accidents, with none exceeding 7.7%. The PRA results for core damage frequency and large early release frequency support the use of NRC Regulation Guide 1.174 as further risk information that supports a change without a full licensing amendment review. The hazard analyses and PRA confirm the need for engineered blast barriers of storage tanks and the common production header leaving the HTEF. The hazards analyses and PRA also confirm with high confidence that using the assumptions of design in this report that the safety case for licensing an HES addition and an HTEF sited with its unprotected high-pressure stage components 187 meters from the NPP’s transmission towers (the most fragile structure, system, and component) is strong.

08 HYDROGEN↗

Development and Integration of a Stochastic Clad Damage Propagation Model into PRONGHORN-SC Subchannel Analysis Code

The failure of fuel pins in nuclear reactors is intrinsically stochastic. Typically, a combination of variation in manufacturing that affects the material characteristics and the fuel assembly dimensions, variation in operating conditions, such as local power, coolant flow rate, and irradiation induced changes in material properties lead to a large uncertainty in failure margin of the fuel pins. Failure, therefore, may occur in exceptional pins with adverse combinations of these variations. Upon a metal fuel pin (U-Pu-Zr/HT9) failure, depressurization of the fuel pin takes place by release of fission gas, liquid sodium bond, and potentially solid fuel particles or molten/eutectic fuel droplets through the hole in cladding. The effect of a fission gas jet on neighbor fuel pins and possible propagation of a clad damage during normal operation was studied experimentally in 1970s and it was found that the post-failure fission gas jet insulates the jet impingement area of the target fuel pin surface and could increase the target pin’s surface temperature by as much as 100 – 200 K during the failed pin depressurization. It was concluded that the effect should not lead to fuel pin failure propagation during normal operation. In accident scenarios of sodium and lead fast reactors such as Unprotected Loss-Of-Flow (ULOF) or Unprotected Transient Over Power (UTOP), the fuel pins can be subjected to higher clad temperatures and fuel pin pressures or fuel clad mechanical/chemical interaction where thermal creep margin becomes significantly lower compared to the normal operation conditions. Therefore, possible stochastic failure and the post-failure fission gas/fuel jet impingement could be critical in order to predict fuel pin failure propagation. Pin depressurization due to fission gas release may degrade the heat transfer by formation of a gas blanket on a neighboring pin surface, which is a local phenomenon, and by causing coolant flow deceleration and starvation, which could affect a surrounding region as well. Furthermore, the potential presence of solid fuel particles or molten fuel at the time of clad failure could boost post-failure jet induced degradation even further. The present study models the U-Pu-Zr/HT9 metal fuel pin failure and stochastic clad damage propagation by biased sampling based on a Cumulative Damage Fraction (CDF) type clad failure criterion and the normal distribution of fuel failure probability density as a function of logarithm of Cumulative Damage Fraction. In addition, the effect of post-failure fission gas jet on heat transfer degradation is modeled for the target pins. This model is called stochastic Clad Damage Propagation (CDAP). The CDAP model is now fully integrated into developmental version of PRONGHORN-SC subchannel analysis code, allowing for modeling local failures and its propagation potential. Section 2 describes the components of the CDAP models. Section 3 describes the model implementation to PRONGHORN-SC and input specifications. Section 4 describes the CDAP model validation coupled to PRONGHORN-SC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Advancing Multi-Hazard Risk and Safety Considerations for Aging Nuclear Facilities (Revision 1)

This project will demonstrate a multi-hazard time-dependent probabilistic risk assessment (PRA) approach for nuclear facilities considering aging-related deterioration of structures. A generic pressurized water reactor (PWR) reactor subjected to seismic mainshock-aftershock sequences considering the aging of the containment structure will be used as a case study to demonstrate the multi-hazard PRA approach. Using advanced modeling and simulation, seismic mainshock-aftershock fragility functions will be simulated for the containment structure considering aging effects. A multi-hazard PRA model for a generic PWR reactor will be built to quantify the multi-hazard core damage frequency (CDF) and large early release frequency (LERF) with explicit time-dependent modeling of event sequences. To date, the cascading impacts of multi-hazards are not adequately accounted for in the PRA models for nuclear facilities. In addition, for both initial and periodic evaluation of facilities to withstand natural phenomena hazards (NPH), deterioration of the SSCs due to aging and other effects is not adequately considered. By advancing multi-hazard considerations accounting for aging effects, this project will contribute to improved understanding of the safety of aging nuclear facilities. Project outcomes such as the multi-hazard CDF and LERF will also allow facility owners to optimize upgrade and retrofit protocols. This project is divided into two components: (1) advanced modeling and simulations; and (2) multi-hazard PRA. For the first component, Idaho National Laboratory’s (INL) Multi-hazard Analysis for STOchastic time-DOmaiN phenomena (MASTODON) and BlackBear codes will be used to simulate the seismic response and damage of a containment structure under cascading mainshock-aftershock sequences with aging effects. Uncertainties related to the seismic inputs, material parameters, and environmental factors such as temperature and humidity will be identified and propagated to the fragility functions. These fragility functions, which consider aging effects, will be time dependent. The second component will take the fragility information from the first component to build a multi-hazard time-dependent PRA model starting from a generic PWR model to evaluate the multi-hazard CDF and LERF and characterize the associated consequences. For this component, OpenPRA Web Application and SAPHIRE code will be used. Events such as the Fukushima Daiichi accident have highlighted the importance of considering the cascading impacts of multi-hazards for PRA. Moreover, many of the reactors in the current nuclear fleet in the United States (US) are already operating well beyond their initially planned design life, and applications for further extensions to operating licenses are being considered. Therefore, considering multi-hazard effects and aging deterioration in the NPH risk assessment process will contribute to the safety of both existing and future nuclear facilities. Outcomes of this project will thus directly benefit the standards DOE-STD-1020-16 and DOE-HDBK-1224-18.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modeling Molten Salt Spreading and Heat Transfer using MELTSPREAD – Model Development Updates

This report summarizes spreading and heat transfer results calculated using MELTSPREAD for a scenario of molten FLiNaK spilling onto a flat stainless steel substrate that was inspired by the “maximum credible accident” scenario of the Molten Salt Reactor Experiment (MSRE). MELTSPREAD was developed at Argonne to model the one-dimensional flowing and freezing of molten corium and was applied to model molten salts for the first time during FY21. The molten salt MELTSPREAD model that was developed during FY21 was updated as part of this work using a corrected heat of fusion value of FLiNaK and to accommodate a larger salt spill volume. The model was run with and without the inclusion of contributions from decay heat and a sensitivity analysis of initial spill conditions was performed to determine the importance of those factors on model outcome. The results provide insight into the expected spreading and heat transfer behavior of simulated and irradiated fuel salt that has been spilled and indicate aspects of the model to be revised in future versions of the MELTSPREAD model for applications to molten salt.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Multicriterion Benefit Evaluation of Deploying New Battery Technology with Increased Capacity at a Generic Nuclear Power

Nuclear power plant (NPP) safety improvements are routinely made by plant licensees and regulators. Examples of plant improvements include accident-tolerant fuel, diverse and flexible coping strategies, passive cooling systems, and increased battery capacity. A combined use of these plant improvements could lead to plant designs with enhanced resilience, allowing NPPs to better cope with both internal and external hazards and keep the plant operating safely, efficiently, and economically. This paper focuses on increased battery capacity and evaluates the potential costs and benefits of deploying batteries with increased capacities at a generic boiling water reactor (BWR) NPP. A multicriterion benefit evaluation methodology is used for the cost-benefit analysis. Ten alternatives for extending battery capacity are developed, including nine alternatives to provide additional direct-current power and one alternative to provide additional alternating-current power. Potential benefits of reducing plant risk are quantified through incorporating the alternatives into loss-of-offsite-power scenarios of the generic BWR probabilistic risk assessment model. Potential costs of implementing the alternatives are qualitatively discussed and ranked. The alternatives are then compared based on their impacts on plant risk and economics.

99 GENERAL AND MISCELLANEOUS↗

Multicriterion Benefit Evaluation of Deploying New Battery Technology with Increased Capacity at a Generic Nuclear Power Plant

This presentation is prepared for the DOE LWRS RISA Pathway Stakeholder Engagement Meeting on October 13-14, 2021. Nuclear power plant (NPP) safety improvements are routinely made by plant licensees and regulators. Examples of plant improvements include accident-tolerant fuel, diverse and flexible coping strategies, passive cooling systems, and increased battery capacity. A combined use of these plant improvements could lead to plant designs with enhanced resilience, allowing NPPs to better cope with both internal and external hazards and keep the plant operating safely, efficiently, and economically. This presentation focuses on increased battery capacity and evaluates the potential costs and benefits of deploying batteries with increased capacities at a generic boiling water reactor (BWR) NPP. A multicriterion benefit evaluation methodology is used for the cost-benefit analysis. Ten alternatives for extending battery capacity are developed, including eight alternatives to provide additional direct-current power and one alternative to provide additional alternating-current power. Potential benefits of reducing plant risk are quantified through incorporating the alternatives into loss-of-offsite-power scenarios of the generic BWR probabilistic risk assessment model. Potential costs of implementing the alternatives are qualitatively discussed and ranked. The alternatives are then compared based on their impacts on plant risk and economics.

99 GENERAL AND MISCELLANEOUS↗

Evaluating transient fission gas release in high burnup light water reactor fuel during loss of coolant accident conditions via new capabilities

In this work, the role of transient fission gas release (tFGR) in the cladding burst behavior of high burnup fuel during a loss-of-coolant accident (LOCA) in commercial light water reactors was further investigated via the use of a new apparatus. During the LOCA-related temperature ramps of high burnup fuel, the release of fission gases exceeds the steady-state release observed under normal operating conditions. An enhancement was made to the Oak Ridge National Laboratory Severe Accident Test Station (SATS) to probe the various factors influencing tFGR. Experiments were performed on commercially irradiated, zirconium-clad uranium dioxide fuel, and this paper details the design of the experimental setup, the initial test results, and the subsequent post-test analyses. Notably, the first test on high burnup fuel demonstrated a LOCA-relevant tFGR of 5.3% from an unpressurized fuel segment. The ultimate tFGR was 10.7% for beyond LOCA conditions. A follow-up test on similar fuel revealed a tFGR of 12.6% under comparable conditions. Microstructural analysis and an analysis of the released gas provide some insight regarding the source of tFGR in the fuel. Finally, a grain boundary bubble model may aid in the interpretation of the results and offer a guide for future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Porous Flow Modeling of Axial Gas Redistribution in Fragmented LWR Fuel Rods using MOOSE

Understanding how gas axially redistributes within fragmented fuel pellets is crucial for predicting the behavior of Light Water Reactor (LWR) fuel rods, particularly during transient and accidental scenarios. The time scale of this phenomenon plays a fundamental role in determining the progression and hazard of a Loss Of Coolant Accident (LOCA), especially when high burn-up fuel in a severe state of fragmentation is involved. Here, this study presents a Computational Fluid Dynamics (CFD) model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) to predict the time-scale of plenum depressurization in Light-Water Reactor (LWR) fuel rods driven by axial gas transport through fragmented pellets. The model examines the effects of incorporating non-linearities in the friction term by comparing the results with experimental data. The experiment employed surrogate fuel rods containing pellets subjected to mechanical and/or thermal loadings to simulate various severity of cracking, and aimed at studying the influence of fuel conditions on axial gas redistribution. The results of this analysis indicate that under certain flow regime conditions - determined by the value of an equivalent Reynolds number - accounting for the non-linear friction term in Navier-Stokes equations guarantees better predictions for the time-scale of plenum depressurization. Also, the model enabled the simulation of the pressure decay by assigning distinct permeability values to each pellet instead of a single uniform value. Multiple simulations were run across all possible pellet position combinations, having each pellet assigned with values of permeability extracted from the experimental data. This allows to quantify the impact of the considering various non-uniform distributions of permeability on the dynamics of axial gas redistribution. The present work findings enhance the understanding of axial gas transport, and provide valuable insights for the integration of a model for predicting the axial gas redistribution during a LOCA scenario into the BISON fuel performance code.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evaluation of Detritiation Strategies for Mitigating Tritium Releases

Tritium confinement is performed using different barriers to minimize releases to the environment. Primary confinement is usually performed by process piping and components, secondary confinement typically by inerted gloveboxes connected to a continuously operating tritium stripper system. Some facility designs employ tertiary tritium confinement by using process room confinement with a tritium stripper system activated after an accident to further mitigate tritium releases to the environment. Secondary confinement of some tritium systems such as radiation-hardened or shielded-cell enclosures are not easily achieved using designs for inerted glovebox. These enclosures can be challenging to seal for secondary tritium confinement due to heating, ventilation, and air conditioning (HVAC) needs for temperature control, as well as feedthroughs for equipment and instrumentation, are potential pathways for tritium leaks or loses. Additional tritium leak paths can be created in secondary confinement enclosures due to a Design Basis Accident (DBA) seismic event. These new, larger tritium leak rates from the DBA event are usually not part of the stripper system design basis so all tritium released to the secondary confinement enclosure is assumed released to the environment. Detritiation strategies for shielded enclosures and/or gloveboxes with significant leak rates are necessary to minimize off-site radiological dose consequences. To address this, the relationship between the confinement system leak rate and stripper system (recirculating) flow rate was evaluated to meet prescribed maximum allowable environmental tritium emissions from the event. The simple analysis described assumes a constant leak rate from the confinement system (i.e. shielded enclosure or glovebox) while a recirculating stripper system strips tritium from the system– a competition between tritium release and tritium recovery. At a high level, the resulting expression summarizes the relationship between system leak rate (F L ) and stripper flow rate (F S ) relative to allowable tritium release (Q A ) and initial tritium release (Q 0 ): $\frac{F_L}{F_s}$ ≈ $\frac{Q_A}{Q_0}$. This report derives the relationship between these parameters, examines the impact of finite stripping times followed by purging of the stripped volume. The analysis provides example results for up to 30 gram tritium releases: the maximum tritium inventory of a Hazard Category III (Department of Energy) Nuclear Facility. The detritiation model presented represents a high-level analysis of tritium recovery versus losses for “leaky” confinement enclosures such as shielded cells after large tritium releases. Key parameters for doing the analyses are the system leak rate, stripper flow rate, initial tritium release, and allowable tritium release values. The analyses show a proportional decrease in releases with reduction in initial tritium release but a non-linear increase in tritium recovery with increased stripper flow rate or reduced leak-rate. The analyses also apply to glovebox confinement systems which could have significant increases in leak rates after a DBA. The analysis recognizes but does not include tritium absorption followed by re-emission from the walls of the confinement volume – an analysis which could be pursued in future studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Assessment and validation of NEAMS tools for high-fidelity multiphysics transient modeling of microreactors: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The NEAMS Multiphysics Applications team aims at providing assessment of code useability and functionality for microreactor design and analyses, together with demonstration of their capabilities to properly capture the steady-state and time-dependent behavior of different microreactor concepts. In FY-24, significant progress was achieved in improving multi-physics models of several microreactors systems: HP-MR, GC-MR and KRUSTY. These efforts focused on solving more complex multiphysics problems enabled by enhanced tools capability, verifying and validating results obtained, providing feedback to developers for suggested improvements, and sharing these models to facilitate user training. A series of new multiphysics transients were completed on the HP-MR (using Griffin/BISON/Sockeye) with core startup transient, control drum inadvertent rotation accident, and hydrogen leakage from hydride moderator (also including SWIFT). On the GC-MR, a new full-core model was developed and analyzed through a series of new multiphysics (Griffin/BISON/SAM) transients to simulate moderator leakage (also including SWIFT), flow blockage and coolant depressurization. Additional and updated TRISO failure analyses were completed on the HP-MR unit-cell and GC-MR assembly models leveraging improved TRISO modeling capabilities. The amount of SiC failure following accidental transients at end-of-life was null. However, GC-MR assembly TRISO analysis highlighted Pd penetration rate can be problematic and may require design changes on the studied microreactor concept. The neutronics discrepancies observed on the KRUSTY model in previous years were resolved using hybrid set of Monte Carlo/Deterministic cross-sections. The multiphysics (Griffin neutronics / BISON thermal-mechanics) 15₵ insertion transient simulation displayed good agreement when comparing with experimental data. Initial modeling of the 30 ₵ reactivity insertion also displays promising results. Such close agreement provides important validation data that can be leveraged by the NEAMS program and by microreactor vendors to support licensing of their technology. Finally, important experience was gathered with the NEAMS tools leading to several user feedback shared with tools developers, especially with regards to MOOSE mesh generator and Griffin. This project led to many publications demonstrating modeling capabilities, and to three models shared on the Virtual Test Bed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

AUTOMATIC GENERATION OF EVENT TREES AND FAULT TREES: A MODEL-BASED APPROACH

In the past few decades, increasing complexity in modern engineering systems has been driven by the integration of a large number of components and by the fact that the system operations involve many disciplines (e.g., thermal-hydraulics, plant operations, cyber-security). Current safety/reliability modeling approaches to such systems are labor intensive, difficult to learn, and rely heavily on simplistic Boolean logic to depict failure propagation and accident progression. While these methods serve well for simple systems (i.e., linear causal systems with limited small inter- and intra-system interactions), their results are difficult to verify when modeling complex systems (typically performed through the extensive use of modeling assumptions). The development of new methods is addressed to meet these challenges through a model-based system engineering (MBSE) lens. Under MBSE philosophy, every aspect of the system (form or function) is represented by a model that completely characterizes its architecture or behavior. MBSE approach greatly improves the management of design, analysis and verification of complex systems. An integration of Dynamic Probabilistic Risk Assessment (DPRA) methods with MBSE models is proposed to perform safety/reliability analyses of engineering systems. In particular, MBSE representation of the system (performed using Systems Modeling Language [SysML]) is coupled with DPRA methods to automatically generate event trees and fault trees.

97 - MATHEMATICS AND COMPUTING↗

Submersible Multistage Centrifugal Pump for Versatile Test Reactor Cartridge Test Loop

Submersible multistage centrifugal pumps are ideal for pumping in narrow confined spaces and achieving necessary head pressures and flow rates. Once a diameter is determined then manipulation of the number of stages and motor speed are all that are required to meet desired flow conditions. The Versatile Test Reactor (VTR) closed loop cartridge systems will need forced convection cooling independent of the main reactor. A multistage centrifugal pump can meet the necessary flow rates and pumping pressures while minimizing space taken. The pump considered for this work was based off a deep well submersible pump, a variation of a multistage centrifugal pump. We experimented with two pump sizes, 5 cm (2 inch) and 7.5 cm (3 inch) diameters. These diameters were chosen to fit into the inner diameter of standard 5 and 7.5 (2 and 3 inch) Schedule 40 pipe, respectively. This made the design for the test loop both simpler and less expensive as the need for an engineered pump housing was eliminated. Initial test cartridge planning indicated space for only a 5 cm (2 inch) diameter pump, though early testing of this size showed the need for an abnormally high-speed and high-power motor. Fine tuning of the cartridge design allowed a pump size increase to 7.5 cm (3 inches), which was the pump size most extensively tested in this work. The test loop is composed of various sizes of PVC and aluminum piping components in a loop configuration. The pump is driven by a Pittman 250 W (1/3 horsepower) electric motor with maximum speed of 3,450 RPM. Testing consisted of running the pump at a constant motor speed while varying a control valve to restrict flow through the loop, with differential pressure and flow rate recorded. This was done for one and two stage configurations for the 5 cm (2 inch) diameter impeller design and one, two, and three stage configurations for the 7.5 cm (3 inch diameter) impeller design, respectively. Due to pumping power requirements, two and three stage 7.5 cm (3 inch) diameter impeller testing at higher flowrates lowered the motor speed substantially. In regions where motor speed could not be maintained constant, the data were discarded. The test loop was also reconfigured to allow for the pump to be tested for pressure drop in a stalled or inoperable (0 RPM) flow condition. Demonstration of adequate natural convection cooling of the test cartridge fuel type is necessary under accident conditions, and this will depend upon the flow resistance through the impeller assembly when the pump is not operating. Thus, accurate knowledge of the effective impeller assembly loss coefficient is important for safety evaluations. The test loop was modified to provide water inlet and outlets on either side of the pump impeller stack, and a metered flow of lab water was provided in order to measure the pressure drop across the cartridges as a function of flowrate. Data from the pump head curve testing developed as part of this work and supported by analysis using pump head affinity laws indicates that a three stage 7.5 cm (3 inch) pump impeller design will meet target requirements for coolant flow within the VTR cartridge sodium cartridge at full power conditions [1] with margin; this corresponds to a flowrate of 45 l/min (12 gpm) at a pressure drop of 6.1 m (20 feet) of water head. The results of the pressure loss measurements across the impeller assembly when the pump is stationary (i.e., at 0 RPM) indicate that the pressure loss coefficient is 0.921 for a two impeller stack configuration; this value is calculated based on the flow velocity through the minimum available flow area within a single stage of the impeller which corresponds to 1.4 cm2.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

FY23 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 FY23. As a continuation of the modeling effort from FY22, this year’s work focuses on improving the RELAP5-3D model developed previously for two-phase flow simulations. The RELAP5-3D model is updated to more accurately capture the heat loss experienced by the facility. The updated model is compared against experimental data for benchmarking purposes of the RELAP5-3D input model. By correctly accounting for heat loss, the updated RELAP5-3D model can now predict the two-phase baseline case more accurately. The onset and the duration of instability are captured well by the model. Furthermore, analyses are performed to better understand the instability mechanism experienced by the flow where the expansion of the boiling boundary in the chimney is studied in details and the fundamental frequencies of the oscillations are obtained. The updated RELAP5-3D model is further compared against four fault conditions, namely the reduction of riser header inlet flow area, depletion of system inventory, blocked riser channels, and static boiling scenario. For each fault condition, minor modifications and tuning are performed to improve the predictions of the model. The purpose of the analyses is to investigate the capability of RELAP5-3D in predicting complex two-phase flows in possible accident scenarios in actual Reactor Cavity Cooling System (RCCS). Overall, the model is able to capture the behaviors and trends of these fault conditions relatively well. Some discrepancies remain between the experimental data and the predictions, many of which are likely due to the differences in the predicted and experimental vapor generation rate. Future work will focus on the continued development of the current RELAP5-3D input model of the NSTF to both improve the accuracy of the model’s predictive capability and continue supporting the experimental program needs. 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↗