Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model development and 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 451 records · Page 25

Prediction of PbLi fluid flow and temperature field in a thermal convection loop for qualification of fusion materials

Two modeling tools, a thermohydraulics code and a computational model in COMSOL Multiphysics, have been developed, tested and then applied to the analysis of fluid flow and heat transfer in a thermal convection loop (TCL). Such a device has recently been used to experimentally evaluate corrosion compatibility of APMT (Fe–21Cr–5Al–3Mo) steel with high-temperature molten eutectic lead-lithium (PbLi) alloy at Oak Ridge National Lab, TN, USA. Here, the 1-D thermohydraulics code allows for rapid calculations of the loop parameters as a function of the applied heat flux. The supplemental COMSOL finite-element computations provide detailed 3-D velocity and temperature field data but are much more time consuming. Both modelling tools demonstrate an excellent agreement in the computed circulation velocity as well as maximum and minimum temperatures. The performed TCL analysis focuses on flow development in the “hot” and “cold” legs, formation of Dean vortices in corner regions, blocking effect of the immersed samples, and radiative and convective cooling effects under the experiment-relevant conditions for Prandtl number Pr = 0.015, Grashof numbers Gr ~10 7 , and Reynolds numbers Re ~10 4 .

36 MATERIALS SCIENCE↗

Comprehensive Cradle to Grave Life Cycle Analysis of On-Road Vehicles in the United States based on GREET

To properly compare and contrast the environmental performance of one vehicle technology against another, it is necessary to consider their production, operation, and end-of-life fates. Since 1995, Argonne’s GREET® life cycle analysis model (Greenhouse gases, Regulated Emissions, and Energy use in Technologies) has been annually updated to model and refine the latest developments in fuels and materials production, as well as vehicle operational and composition characteristics. Updated cradle-to-grave life cycle analysis results from the model’s latest release are described for a wide variety of fuel and powertrain options for U.S. light-duty and medium/heavy-duty vehicles. Light-duty vehicles include a passenger car, sports utility vehicle (SUV), and pick-up truck, while medium/heavy-duty vehicles include a Class 6 pickup-and-delivery truck, Class 8 day-cab (regional) truck, and Class 8 sleeper-cab (long-haul) truck. Powertrain coverage includes internal combustion (spark ignition and compression ignition) engines, hybrid electric, plug-in hybrid, full battery electric, and fuel cell vehicles powered by conventional and low carbon energy sources. The results offer insights into the current state of these technologies, as well as a projection of the likely environmental implications of future fuel and vehicle advancements through a time-series evaluation of life cycle greenhouse gas emissions.

Kelly, Jarod C.↗

ROM-Based Surrogate Systems Modeling of EBR-II

We report the System Analysis Module (SAM), developed and maintained by Argonne National Laboratory, is designed to provide whole-plant transient safety analysis capabilities for a number of advanced non-light water reactors, including sodium-cooled fast reactor (SFR), lead-cooled fast reactor (LFR), and molten salt reactor (MSR)/fluoride-salt-cooled high-temperature reactor (FHR) designs. SAM is primarily constructed as a systems-level analysis tool, with the potential to incorporate reduced order models from three-dimensional computational fluid dynamics (CFD) simulations to improve characterization of complex, multidimensional physics. It is recognized that the computational expense associated with CFD can be intractable for various engineering analyses, such as uncertainty quantification, inference, and design optimization. This paper explores the reducibility of a SAM model using recent advances in randomized linear algebra techniques, which attempt to find recurring patterns in the various realizations generated by a model after randomly perturbing all its input parameters. The reduction is described in terms of fewer degrees of freedom (DOFs), referred to as the active DOFs, for the model variables such as input model parameters and model responses. The results indicate that there is significant room for additional reduction that may be leveraged for additional computational gains when employing SAM for engineering-intensive analyses that require repeated model executions. Different from physics-based reduction approaches, the proposed approach allows one to estimate upper bounds on the reduction errors, which are rigorously developed in this work. Finally, different methods for surrogate model construction, such as regression and neural network-based training, are employed to correlate the input and output active DOFs, which are related back to the original variables using matrix-based linear transformations.

42 ENGINEERING↗

A Comprehensive Chemistry Evaluation and Diagnostics Package for E3SM – ChemDyg Version 1.1.0

The Chemistry Evaluation and Diagnostics Package (ChemDyg) is an open-source tool designed for the Energy Exascale Earth System Model (E3SM) developed by the U.S. Department of Energy. ChemDyg facilitates routine evaluation, tailored development, and in-depth analysis of atmospheric chemistry through its modular architecture, allowing users to compare model outputs with observational data. Version 1.1.0 introduces a robust set of diagnostic capabilities, including climatology, time evolution of key tracers, diurnal and annual cycle analyses, and extensive budget diagnostics. These features help identify model discrepancies and enhance the representation of atmospheric chemistry in E3SM. Each self-contained diagnostic set includes dedicated scripts and documentation for ease of use. The interactive HTML output improves data accessibility, accelerating chemistry model development. Additionally, ChemDyg's flexible framework allows for customization, enabling users to create unique diagnostic sets for specific scientific contributions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Practical Insights on Applying Simulation-Based Control Methods in Experimental Studies

Advanced nuclear reactors are crucial to the future of energy both in the United States and around the globe. In contrast to the current operating fleet, they are characterized as being deployable in remote locations and able to operate in semi-autonomous or autonomous fashion. This leap forward necessitates a new reactor control paradigm. Because advanced nuclear reactors are still under development in the United States, the creation of new control methods to achieve autonomous operations has been based on systems modeling and simulation. However, an important factor in successfully deploying these new control methods is the ability to seamlessly transition from simulation environments to real-world settings. Control methods tested in both simulation and experimental settings need to be investigated in the context of advanced reactor applications. This work developed a series of simple controllers for Idaho National Laboratory (INL)’s Microreactor Applications Research Validation and Evaluation (MARVEL) microreactor operating in load-following scenarios. These controllers were tested in both simulation and experimental settings, and a comparative performance analysis was performed. The simulation tests leveraged the Control and Optimization Modular Modeling Application for Nuclear Deployment (COMMAND) software developed in a previous stage of the current effort, along with the MARVEL Reactor Excursion and Leak Analysis Program (RELAP5-3D) and Monte Carlo N-Particle (MCNP) models. The experimental tests leveraged the COMMAND software, MARVEL models, and the U.S. Department of Energy Microreactor Program’s Microreactor Automated Control System (MACS). MACS was developed to serve as a control method testbed. It was customized to mirror the MARVEL microreactor, and COMMAND enabled MACS to emulate the physics of MARVEL. The load-following controller was developed using the simulation platform, with efforts to emulate real systems by introducing actuator saturation and noise. These factors were incrementally accounted for in the controller design. After finalizing the controller design, it was implemented with the experimental setup. The experimental conditions tested included an initial test under conditions similar to the final simulation test, and two additional scenarios. The first scenario introduced additional actuator saturation to account for equipment aging over time, which was unknown to the controller. The second scenario introduced sensor delay, a phenomenon anticipated with the use of remote operations or wireless communication in advanced reactors. These tests revealed several notable differences. While the controller performed well in simulation, it exhibited several limitations when transitioning to hardware. The main challenges involved maintaining the steady-state target power, as evidenced by larger error values between the true reactor power and setpoint power, as well as persistent oscillations in controlled reactor power. These issues could lead to unacceptable transient conditions in real reactor testing. Introducing actuator aging and stochastic delays in the experimental setup significantly impacted controller performance, resulting in increased overshoot and undershoot, and exacerbated error and oscillations previously mentioned. These findings underscore the importance of experimental testbeds for testing and validating control methods, as controllers developed using only theory and/or simulation may perform unexpectedly when applied to actual hardware. This research emphasizes the need for an experimental testbed for achieving such validation.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Real-Time Drilling Optimization System for Improved Overall Rate of Penetration and Reduced Cost Per Foot in Geothermal Drilling

The key to success in geothermal drilling is economic feasibility, and a major cost in the development of geothermal resources is the actual drilling of the wells. In this project, a real-time drilling optimization system for geothermal drilling was developed. The system couples three individual components while drilling. The first component is a drill stem vibration analysis model, the second is Mechanical Specific Energy (MSE) analyses, and the third is a detailed PDC Rate of Penetration (ROP) drill bit model for optimum RPM and WOB combinations. The benefit of the coupled system is that the range of WOB and RPM could be selected to avoid drill stem vibrations. Secondly, MSE is used as an efficiency measure and the detailed PDC drill bit model ensures the drill bit does not endure temperatures that exceed the temperature at which the PDC cutters experience accelerated wear. The new detailed PDC bit model is based on rock/bit interaction that physically tracks the PDC cutter wear flats as the bit drills ahead giving the capability to calculate the temperature being generated underneath the worn cutters to better advise on operational parameters to avoid accelerated cutter wear and failure and to ensure that operational parameters are applied so that overall ROP is maximized. By combining the drill stem vibrations and the detailed PDC bit cutter wear and “safe” non-accelerated cutter wear temperature and optimum ranges of operating parameters, it results in higher ROP and lower cost drilling. Single cutter PDC testing performed in different lithologies at Sandia was utilized to verify the PDC cutter forces and depth of cut for new and worn cutters. Based on single cutter PDC temperature modeling, verification using single cutter data from the testing done by National Oilwell Varco (NOV) was performed. Sandia’s Hard-Rock Drilling Facility (HRDF) was utilized to test different drill bit configurations with different cutter designs and wear status with different induced modes of vibration to obtain the critical bit RPM/WOB ranges resulting in ineffective drilling and low ROP. The collected test data were further used to verify and calibrate the full hole PDC ROP model that was developed based on single cutter interaction data. A full coupled drill stem vibration model was formulated and verified with geothermal field data from the Chocolate Mountain Aerial Gunnery Range (CMAGR). A graphical user interface (GUI) was developed using Tkinter library in the computer programming language Python, which integrates all the developed models in one system. The developed system consists mainly of the PDC ROP model, PDC bit wear model, PDC cutter temperature model, Mechanical Specific Energy (MSE) model, and drillstring vibration model integrated into one system. The developed system can be used for both, post well analysis and real-time optimization using different criteria such as ROP maximization or MSE minimization. The software uses Differential Evolution Algorithm (DEA) to find optimum values for operational parameters based on last foot drilled while avoiding the drillstring vibration and cutter temperature critical operating parameters.

15 GEOTHERMAL ENERGY↗

From multivariate to functional data analysis: Fundamentals, recent developments, and emerging areas

Functional data analysis (FDA), which is a branch of statistics on modeling infinite dimensional random vectors resided in functional spaces, has become a major research area for Journal of Multivariate Analysis. We review some fundamental concepts of FDA, their origins and connections from multivariate analysis, and some of its recent developments, including multi-level functional data analysis, high-dimensional functional regression, and dependent functional data analysis. Here, we also discuss the impact of these new methodology developments on genetics, plant science, wearable device data analysis, image data analysis, and business analytics. Two real data examples are provided to motivate our discussions.

97 MATHEMATICS AND COMPUTING↗

Load Composition Analysis in Support of the NERC Load Modeling Task Force. 2019-2020 Field Test of the Composite Load Model

In 2015, NERC’s reliability standards were revised to require the use of dynamic load models in transmission planning studies. To comply with the standards, planners must use load models that explicitly represent the dynamic behavior of the different constituents of load at each load bus within their transmission planning models. The most important of these constituents are motor-driven and power electronics-based loads. Collectively, these representations are known as composite load models. In anticipation of the compliance date for the new standards, NERC’s Load Modeling Task Force (LMTF), in 2019, initiated a field test of composite load models involving the regional reliability planning entities. In support of the field test, DOE and BPA researchers developed region-specific composite load models that could be assigned to each non-industrial load bus in the planning models for each of the North American interconnections. Separate models were developed for each hour of a summer peak day, a winter peak day, and a spring light-load day. This report is the technical documentation for the load composition analysis that was conducted to develop these non-industrial composite load models.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Impact of system parameters and geospatial variables on the reliability of residential systems with PV and energy storage

A reliable power supply is the foundation of modern society, enabling technologies used to function within a society. Residential systems are places where the end users directly consume power, enabling technologies to sustain life. With the emergence of behind-the-meter resources, the end-users have some control over power supply reliability. The intermittency and variability of these resources impact residential system reliability. In this work, we study the reliability performance of a grid-supplemented residential system with behind-the-meter Distributed Energy Resources(DER) subject to various system parameters and geospatial variables. We propose a multistate reliability model for the behind-the-meter microinverter-based Photo Voltaic (PV) system and integrated inverter-based energy storage (ES) system. A sequential Monte Carlo method is then presented to evaluate the reliability indices for the residential system with behind-the-meter DERs as the main supply and the grid as the backup. The sequential Monte Carlo method is used to analyze the reliability performance of the residential model at the top 100 populous counties of the United States, where actual load and solar irradiance data at the counties is used. In the analysis, the sensitivity of indices to system parameters such as DER size and the sensitivity of the indices to climate zone and gross horizontal irradiation (GHI), which affects the load conditions and the PV output of the residential system, is performed. The analysis results show that the system parameters and the geospatial variables significantly impact the residential system’s reliability. The insights from this analysis will be of immense value to the distribution system planners to provide zone-specific guidelines for DER system sizing and toward the evolution of utility business models. Finally, the methodology developed can be used to extend the analysis to other locations.

14 SOLAR ENERGY↗

CO 2 zonal injection rate allocation and plume extent evaluation through wellbore temperature analysis

Temperature analysis during a pause in injection operations, known as warmback analysis, has been used in the petroleum industry for evaluating the injection conformance and estimating the location of the flooded front in applications, such as waterflooding oil reservoirs. Here in this work, methods are introduced to extend the application of temperature warmback analysis to estimate the zonal CO 2 injection rate and zonal CO 2 plume extent during geologic CO 2 storage in a saline aquifer. First, novel analytical solutions are developed to model transient temperature in the aquifer during the injection and subsequent shut-in periods considering two-phase flow (gaseous CO 2 and aqueous brine) conditions in the aquifer. The solution involves a discretization of the aquifer into regions; the energy and mass conservation equations for the regions are solved simultaneously considering appropriate boundary conditions at the interfaces. Two solutions techniques are presented: multi-region and three-region solutions. Inverse models are developed accordingly to evaluate the injection profile and estimate the extent of the plume front in the reservoir during the injection period. The multi-region solution results in an inversion approach that requires regression analysis. However, the three-region formulation results in a simple graphical technique for inverse modeling. The analytical solutions are validated against a thermally coupled reservoir simulation tool using different synthetic cases for CO 2 injection in deep saline aquifers. The results of the developed solutions provide a good match with numerical results during forward and inverse modeling.

54 ENVIRONMENTAL SCIENCES↗

Building an Efficient Cluster Cosmology Software Package for Modeling Cluster Counts and Lensing

We introduce a software suite developed for galaxy cluster cosmological analysis with the Dark Energy Survey Data. Cosmological analyses based on galaxy cluster number counts and weak-lensing measurements need efficient software infrastructure to explore an increasingly large parameter space, and account for various cosmological and astrophysical effects. Our software package is designed to model the cluster observables in a wide-field optical survey, including galaxy cluster counts, their averaged weak-lensing masses, or the cluster's averaged weak-lensing radial signals. To ensure maximum efficiency, this software package is developed in C++ in the CosmoSIS software framework, making use of the CUBA integration library. We also implement a testing and validation scheme to ensure the quality of the package. We demonstrate the effectiveness of this development by applying the software to the Dark Energy Survey Year 1 galaxy cluster cosmological data sets, and acquired cosmological constraints that are consistent with the fiducial Dark Energy Survey analysis.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A time-domain electromechanical co-simulation framework for power system transient analysis with retainment of user defined models

This paper proposes a time-domain electromechanical co-simulation framework to accomplish concurrent running of multiple simulators in order to preserve user-defined models and allow maximum simulation accuracy for power system transient analysis. Here, a middleware interface is developed to manage the interaction and synchronization of the multiple simulators via boundary model. Consequently, it eliminates the requirement of model conversion of each subsystem and thus ensure the user-defined models in each specific software are retained for entire system situational awareness. Preliminary testing results on 140-bus Northeastern Power Coordinating Council (NPCC) model validate the effectiveness and accuracy of the proposed framework for PSS/E and PSLF co-simulation.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Pumped-Storage Hydropower using Abandoned Underground Mines (PSH-AUM) as an Innovative Energy Storage Technology for Fossil-Integrated Systems

Pumped-storage hydropower (PSH) provides around 95% of all utility-scale energy storage in the U.S. and globally. Additional deployment of PSH is hampered by excessively long permitting and commissioning requirements and is constrained to locations for which natural topography provides suitable elevation relief between the upper and lower reservoirs (the ΔH challenge). The purpose of this research was to evaluate Pumped-Storage Hydropower using Abandoned Underground Mines (PSH-AUM) as a means to solve the ΔH challenge and initiate the commercialization pathway for a promising new energy storage technology. Four primary tasks were conducted: Screening and ranking of candidate sites for project development; multiphase reservoir modeling to evaluate mine performance; techno-economic analysis and preliminary designs for PSH-AUM systems integrated with fossil-fuel power units; and stakeholder engagement to identify pathways to commercialization of this new technology. Key results were achieved in each of the four primary tasks. Candidate site screening determined that nearly 10,000 underground mines were spatially locatable, of which more than 100 sites appear suitable for integration with existing fossil power plants that are expected to remain in longer-term operation. Mine reservoir models were developed using PNNL’s STOMP simulator and parameterized using candidate site data to evaluate interactions with the surrounding groundwater system and confirm the potential for some mines to accommodate inflows and outflows on the order of 100 m3/s over 8-hour durations (1.6 GWh system) without excessive aqueous pressures. Techno-economic analysis resulted in a project cost optimization scheme to identify key sensitivities and the development of preliminary designs that can minimize overall costs of deployment. Finally, stakeholder engagement with industry, state government, and local economic development leaders confirmed the viability of PSH-AUM as a promising new technology. Our Phase I project results suggest that PSH-AUM technology has a domestic market potential on the order of $100+ Billion with ample space for technology development to commercialization within the next decade.

13 HYDRO ENERGY↗

Automated ICRF heating surrogate modeling via machine learning

This work introduces automated machine learning workflows that address critical bottlenecks in surrogate model development for Ion Cyclotron Range of Frequencies (ICRF) heating applications. The automated framework includes data analysis tools that transform raw datasets into actionable insights in seconds, replacing weeks of manual exploratory effort and ensuring consistent, reproducible dataset characterization. By integrating advanced hyperparameter optimization (HPO) methods including Bayesian optimization via BoTorch and Tree-structured Parzen Estimators (TPE), the framework significantly reduces model development time from weeks to hours, decreasing computational cost and required expertise, while enabling high-accuracy surrogate models. Compared to traditional hyperparameter scanning (HPS) techniques such as methodical, randomized, and grid searches, HPO methods achieve superior convergence and predictive performance, even when compared to already well-tuned reference models. On NSTX High Harmonic Fast Wave (HHFW) heating datasets, both Random Forest Regressor (RFR) and neural network surrogates demonstrate improved accuracy, achieving R 2 values beyond 0.97 and 0.98, respectively. The results show that while HPO gains are modest for robust architectures like RFR, they become essential for more sensitive models such as neural networks, highlighting the trade-offs across optimization strategies. Through automated workflows that eliminate manual hyperparameter tuning and require minimal ML expertise, this work enables widespread adoption of high-fidelity surrogate models across the fusion community for real-time plasma control, uncertainty quantification, rapid experimental scenario development, and integrated system optimization.

Sanchez-Villar, Alvaro [Princeton Plasma Physics L↗

Advanced Measurements for Resilient Integration of Inverter-Based Resources: PROGRESS MATRIX Final Report

As nearly every aspect of the electric power grid undergoes rapid change, measurement technologies that support grid operation and planning must evolve as well. The rapid large-scale deployment of inverter-based resources (IBRs) vital to achieving the nation’s clean energy goals has in some cases led to negative impacts on the reliability and security of the bulk power system (BPS). Advanced power system measurements, including synchronized phasor and waveform measurements, are key to making IBR integration secure and reliable. To this end, the Department of Energy (DOE) initiated the PROGRESS MATRIX project to develop advanced measurement capabilities and analytics that will accelerate adoption of IBRs while improving the reliability and resilience of the BPS. This report discusses the outcomes of the project, which was a joint effort between the Pacific Northwest National Laboratory (PNNL), Oak Ridge National Laboratory (ORNL), the National Renewable Energy Laboratory (NREL), and Lawrence Berkeley National Laboratory (LBNL). In the project’s first year, PNNL, NREL, and ORNL partnered with the Bonneville Power Administration (BPA), the Western Area Power Administration (WAPA), and Kauai Island Utility Cooperative (KIUC) to understand their existing measurement capabilities and the gaps limiting deployment of IBR-focused measurement systems and analytics. The other primary activity in the first year was deployment of GridSweep instruments, which provide unprecedented precision in waveform measurement while probing distribution systems. The instruments were deployed at Dominion Energy and the University of California, Riverside. In the project’s second year, the input from partner utilities and collected measurements were used to advance measurement capabilities. Twelve analytical methods spanning disturbance analysis, power plant evaluation, feeder evaluation, and modeling were developed. Two software tools were developed, one to analyze GridSweep measurements and another to automatically evaluate the control performance of power plants connected to the BPS. Testbeds at ORNL and NREL were augmented to better enable studies of IBR integration. The project culminated in demonstrations of these analytical methods, software tools, and testbeds, both in the field and in the laboratory. This report discusses these various accomplishments and documents the significant progress in developing advanced measurement capabilities to support the secure, reliable, and accelerated adoption of IBRs in the BPS.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Effect of Anoxic Iron Corrosion on WIPP Brine Geochemistry FY23 Final Report (U)

A 280-day study was completed to evaluate the effect of zero-valent iron (Fe 0 ) on the Waste Isolation Pilot Plant (WIPP) brine geochemistry under anticipated reducing conditions. Hydrogen (H 2 ) gas is expected to be present in the repository after closure due to the anoxic corrosion of a vast quantity of iron contained in the waste forms disposed at WIPP; therefore, a background argon atmosphere containing H 2 was chosen for this study. WIPP groundwater brine pH and E h will impact the mobility and fate of plutonium within the repository. Modeling and laboratory results for Castile WIPP brine indicate that equilibrium fa values relative to the standard hydrogen electrode (SHE) are 40 mV more reducing (i.e., more negative) than those for Salado WIPP brine (-480 mV vs. -440 mV, respectively) because of the higher pH of the Castile brine (pH 9 .3 for Castile vs. pH 8.8 for Salado). The E h and pH data were corrected for the effects of high ionic strength. The experimental results for both brines are consistent with thermodynamic predictions using OLI Systems' Mixed Solvent Electrolyte chemical equilibrium model. The measured and corrected pH and E h data from this study are provided in Table ES-I and Table ES-2, respectively. The experimental study, with four test conditions in triplicate, was performed in a dual glovebox with a nominally 3 vol.% H 2 in argon atmosphere (target H 2 range: 3 ± I vol.%). Simulants containing MgO only ( experimental control) and MgO+Fe 0 (WIPP base case) were prepared for both the Salado and Castile brines. MgO was included in all simulants to account for the use of bulk magnesium oxide in the WIPP repository. Fe 0 was included in some simulants to incorporate the effects of the anoxic corrosion of iron and in-situ hydrogen generation in the study. The brine compositions were developed by Sandia National Laboratory (SNL; Xiong, 2008) and have been used in previous WIPP evaluations. The test method (agitation, etc.) is partially based on ASTM D3987-12. Twelve rounds of periodic measurements of pH and E h were performed over the course of the study. Chemical analysis results for liquids and solids (ICP-MS, ICP-ES, IC Anion, TIC, SEM-EDX) are consistent with the pH, E h , and thermodynamic modeling results. This study included the following conditions that deviate from anticipated post-closure conditions following brine intrusion, but were selected to facilitate bench-scale testing to validate modeling of pH and E h for the post-closure WIP P repository: an anoxic glove box atmosphere containing ≤ 4 vol. % H 2 vs. substantially higher H 2 gas concentrations assumed in the WIPP Performance Assessment (PA); a significantly higher liquid-to-solid test ratio compared to the much lower phase ratio anticipated in the WIP P repository; agitation of the simulant bottles to maximize mass transfer; and finally the use of Fe 0 reagents having a much greater surface area than expected in the WIP P repository. Non-representative conditions were chosen for various reasons such as: to provide bounding conservative results, to provide a margin of safety for testing, or to facilitate simulant sub-sampling and analysis. In a parallel effort, aqueous electrolyte thermodynamic models were developed for the synthetic Salado and Castile brines to inform the experimental design, facilitate laboratory data interpretation, and allow extension of evaluations beyond the parameters tested. Thermodynamic modeling simulations including the MgO and Fe 0 additives that are directly relevant to the experimental measurements (e.g., pH calibration curve, ORP corrections) are included in this report. The measured fa of the simulants was close to the OLI model predictions for both brines and was largely controlled by the background H 2 partial pressure in the vapor phase as well as H 2 generated in situ in the aqueous phase by the Fe 0 corrosion. The H 2 gas-phase concentration tested and thermodynamically evaluated was much lower than is assumed in the WIPP PA; however, H 2 (g) concentrations significantly below this level are still predicted to result in very reducing conditions. In conclusion: • The experimental results are consistent with thermodynamic model predictions for fa, pH, and the effects of high ionic strength. • Evidence to date suggests that the H2 concentration in the glovebox atmosphere ultimately determined the final E h values of the simulants and resulted in highly reducing conditions. As a result, little difference was observed between the control simulants containing only MgO and the WIPP base-case simulants that contained MgO and Fe 0 . • This test methodology is recommended for future studies evaluating WIPP repository conditions. The methodology includes: (1) background H 2 in argon with agitation ( or could alternatively include in-situ-generated H 2 in sealed bottles); (2) carefully measured and corrected ORP data ( with much effort focused on allowing the probes to fully stabilize); and (3) ionic-strength-corrected pH data. Other best practices, such as simulant sparging/handling, ORP probe replacement, etc., should also be considered. • The coupling of experimental studies and thermodynamic modeling is also highly recommended because these methods inform and direct one another leading to greater confidence in and understanding of the results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Risk Importance Ranking of Fire Data Parameters to Enhance Fire PRA Model Realism

Fire is historically and analytically a significant contributor to nuclear power plant risk. The level of fire risk and the methods, tools and data for modeling this risk is highly debated by experts. One area of debate is the input data used in fire modeling and how to deal with this data’s high uncertainty. This report outlines initial work performed for determining the key parameters causing this uncertainty and how it propagates into nuclear power plant models. This research paves the way for the development of methods to reduce fire data uncertainty used in modeling. The Nuclear Regulatory Commission has mandated that nuclear power plants perform fire risk modeling. However, there are several issues with the current risk modeling implementation that affect the results. Approved modeling methods can be overly conservative and often do not match plant experience. Also, the data used in the modeling can have high uncertainties and is influenced by expert judgement. To evaluate input data uncertainty, researchers performed an initial review of several fire experiments done at Sandia National Laboratories. Uncertainties for fire data can come from many sources, such as experiment design constraints, environmental conditions, or other plant-specific aspects. There are many different significant and insignificant parameters driving the uncertainty. Additionally, the uncertainty of the different input data used in the fire modeling could have a significant or insignificant effect on the entire plant risk. A four-step methodology was developed to perform Integrated Probabilistic Risk Assessment Importance Ranking. A demonstration case using these steps was set up and three of the four steps were completed in fiscal year (FY) 2019 and the fourth step done FY 2020. These steps are: 1. The qualitative analysis of potential sources was conducted with the following items identified for the demonstration. • Maximum heat release rate • Time to maximum heat release rate • Duration of max heat release rate • Time to decay • Thermal conductivity of concrete • Specific heat of concrete • Density of concrete • Cable jacket thickness 2. A quantitative characterization of dominant sources of uncertainty was performed. A list of distributions and determined values of the dominant sources is shown in Appendix A. 3. A quantitative screening of the potential sources of uncertainty using Morris Elementary Effects Analysis was completed. An experimental model using the physics-based fire modeling tool Fire Dynamics Simulator was developed and coupled with the Risk Analysis Virtual Environment. The Morris analysis identified at least two parameters that can be eliminated as significant contributors (specific heat of concrete and cable jacket thickness). 4. Global importance measure (Global IM) analysis to generate a comprehensive ranking based on their influence on the plant risk. In this research, a moment-independent Global IM is used since it can address (a) uncertainty in the input parameters of the fire model, (b) uncertainty in the risk outputs, and (c) non-linearity and interactions among input parameters in the fire model, more accurately than the correlation-based and variance-based global methods. The observations from the research showed that, depending on the initial and boundary conditions of the fire scenarios, fire-induced damage could have a very small probability and could be dominated by the tail of the uncertainty distribution; hence, the accuracy of the correlation-based and variance-based methods is questionable. The moment-independent Global IM analysis in this research provides a better understanding of how experimental uncertainty data affects industry’s plant models and where improvements in that data will have the largest benefit for improving fire modeling accuracy in causing core damage. Among the five unscreened parameters obtained from the Morris EE analysis, the Global IM analysis results for the case study indicated that max heat release rate and fire location are the most important parameters. The report also outlines benefits of using a unified computational platform that integrates the underlying simulations (e.g., a fire progression model), quantitative screening (using the Morris EE method), and the Global IM analysis. A unified platform can (i) facilitate the ranking of input parameters considering multiple key fire scenarios simultaneously, rather than considering one scenario at a time, (ii) contribute to more explicit and accurate treatment of dependencies at multiple levels of Fire PRA, (iii) facilitate the sampling-based uncertainty quantification for Fire PRA, and (iv) help generating both “industry-wide” and "plant-specific" ranking of uncertainty sources in Fire PRA. Future research should be done to include additional parameters such as detection/suppression or cable fire spread. Adding a PRA software such as SAPHIRE to the RAVEN platform would help with plant model integration and improve treatment of fire-induced dependency. The I-PRA risk importance ranking methodology offered in this report can provide valuable information for efficiently (a) enhancing the realism of Fire PRA for existing plants and (b) supporting the development of Dynamic Fire PRA for advanced reactors and new plants.

97 MATHEMATICS AND COMPUTING↗

Coupled SAM/Griffin Model of a Reference Fluoride-Salt-Cooled High-Temperature Reactor for Multi-Physics Simulations

A multi-physics coupled simulation model of a reference pebble bed fluoride-salt-cooled high-temperature reactor (PB-FHR) has been developed with SAM and Griffin computer codes for transient safety analysis. The reference problem of a prototypical reactor design serves as the foundation for the U.S. NRC (Nuclear Regulatory Commission) to verify the adequacy of computer codes and evaluation models for specific reactor types. In this work, the previously developed SAM model for PB-FHR primary loop has been updated for the coupled simulation. The updated SAM PB-FHR model includes a 2-D axial symmetric core region and external core components in 0-D/1-D. In addition to the primary loop, a detailed model of the RCCS (reactor cavity cooling system) is added. The 2-D and 1-D domains are tightly coupled using the single-solve approach developed in SAM. In the pebble bed region, the SAM multiscale explicit pebble model is applied to calculate the pebble and TRISO fuel kernel temperatures. The Griffin model used in this work is based on a model developed at Idaho National Laboratory in collaboration with the U.S. NRC. The Griffin neutronics model and SAM thermal hydraulics model is coupled with the Comprehensive Reactor Analysis Bundle (CRAB or alternately BlueCRAB) application. Both steady-state and transient scenarios are simulated to demonstrate the model's suitability for multi-physics simulations of PB-FHR transients.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗