Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “empirical validation”

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 145 records · Page 8

Co-design optimization of combined heat and power-based microgrids

With the emergent need for clean and reliable energy resources, hybrid energy systems, such as the microgrid, are widely adopted in the United States. A microgrid can consist of various distributed energy resources, for instance, combined heat and power (CHP) systems. Here, the CHP module is a distributed cogeneration technology that produces electricity and recaptures heat generated as a by-product. It is an energy-efficient technology converting heat that would otherwise be wasted to valuable thermal energy. For an optimal system configuration, this study develops a novel co-design optimization framework for CHP-based cogeneration microgrids. The framework provides the stakeholder with a method to optimize investments and attain resilient operations. The proposed co-design framework has a mixed integer programming (MIP) model that outputs decisions for both plant designs and operating controls. The microgrid considered in this study contains six components: the CHP, boiler, heat recovery unit, thermal storage system, power storage system, and photovoltaic plant. After solving the MIP model, the optimal design parameters of each component can be found to minimize the total installation cost of all components in the microgrid. Furthermore, the online costs from energy production, operation, maintenance, machine startup, and disruption-induced unsatisfied loads are minimized by solving the optimal control decisions for operations. Case studies based on designing a CHP-based microgrid with empirical data are conducted. Moreover, we consider both nominal and disruptive operational scenarios to validate the performance of the proposed co-design framework in terms of a cost-effective, resilient system.

42 ENGINEERING↗

Enhancing high-fidelity neural network potentials through low-fidelity sampling

The efficacy of neural network potentials (NNPs) critically depends on the quality of the configurational datasets used for training. Prior research using empirical potentials has shown that well-selected liquid–solid transitional configurations of a metallic system can be translated to other metallic systems. This study demonstrates that such validated configurations can be relabeled using density functional theory (DFT) calculations, thereby enhancing the development of high-fidelity NNPs. Training strategies and sampling approaches are efficiently assessed using empirical potentials and subsequently relabeled via DFT in a highly parallelized fashion for high-fidelity NNP training. Our results reveal that relying solely on energy and force for NNP training is inadequate to prevent overfitting, highlighting the necessity of incorporating stress terms into the loss functions. To optimize training involving force and stress terms, we propose employing transfer learning to fine-tune the weights, ensuring that the potential surface is smooth for these quantities composed of energy derivatives. This approach markedly improves the accuracy of elastic constants derived from simulations in both empirical potential-based NNPs and relabeled DFT-based NNPs. Overall, this study offers significant insights into leveraging empirical potentials to expedite the development of reliable and robust NNPs at the DFT level.

97 MATHEMATICS AND COMPUTING↗

Demonstrate new plasticity models for doped UO 2 that capture dislocation mechanisms

In light water reactors, fuel vendors are investigating the use of dopants to modify the properties of UO 2 pellets, with the goal of improving pellet-cladding mechanical interactions during operation. Dopants are expected to ‘soften’ the pellets; that is, the doped pellets have higher plastic deformation than conventional UO 2 . This leads to a reduction in the severity of mechanical pellet-cladding interactions, helping to reduce the hoop strain on the cladding. By minimizing the strain exerted by the pellet on the cladding, it is anticipated that cladding performance under accident conditions can be enhanced (i.e., lowering the risk of burst during a LOCA). Dopants such as chromium (Cr) promote grain growth during pellet fabrication, leading to larger grains; therefore, understanding the link between chemistry, microstructure and mechanical deformation (enhanced creep rates) behavior of UO 2 is critical to helping operators further substantiate the benefits of doping UO 2 . Historically, the nuclear energy industry has relied on empirical models to make assessments of performance. Compared to empirical models, mechanistic physics-based models provide benefits, such as, fewer data points for validation and better extrapolation where experimental data is scarce or non-existent. In this report, Bayesian inference techniques have been applied to a previously developed lower length-scale-informed diffusional creep model. The objective is to i) infer lower-length-scale parameter distributions from available experiment and then ii) determine the uncertainties in the measurable quantity (in this case creep rates) after propagating the inferred lower length scale parameter uncertainties. The approach requires many evaluations of the model, which becomes computationally insurmountable; therefore, a neural-network model is trained to data obtained by sampling the full model over the most important parameters. This neural-network is then used in the Bayesian inference approach to determine probability distributions in the parameter values that represent the uncertainty in the model given what is known from the experiments (posterior). A significant reduction compared to conservative initial (prior) uncertainties is achieved through inference against the experimental data, demonstrating the efficacy of this approach. Furthermore, by accounting for uncertainties in the experimental conditions and sample non-stoichiometry, it is possible to resolve apparent discrepancies in experimental measurements within a self-consistent grain boundary (Coble) creep model that is sensitive to chemistry. This work has been written up and submitted to Nuclear Technology for a special issue on accelerated fuel qualification (AFQ). This uncertainty quantification (UQ) work not only improves the diffusional model, while accounting for uncertainty, but also establishes a framework which can readily be applied to the mechanistic models of dislocation deformation developed in this study. The most likely values from the Bayesian analysis are incorporated into our UO 2 diffusional creep model and a lower length scale-informed irradiation UO 2 creep mechanistic model to generate a dataset. This dataset has been provided to our INL collaborators for training an artificial neural network surrogate model, which will be implemented in the BISON fuel performance code to assess how the results differ from those currently obtained using a fully empirical model and that of using the nominal (uncalibrated) atomic scale parameters in our mechanistic model. Plastic deformation (creep and glide) in UO 2 is a complex phenomenon, governed by multiple underlying processes such as local defect concentrations, applied stresses, and microstructural characteristics. Consequently, there is a need for a meso-scale model with polycrystalline resolution capable of extrapolating to large grain sizes applicable to doped UO 2 , where data is limited and the model can help bridge the knowledge gap. By integrating atomistic data into the polycrystal LApx code, it becomes possible to predict dislocation climb and glide plasticity that simple analytical models cannot accurately represent. The application of atomic-scale data within LApx demonstrated the importance of climb and glide mechanisms in reproducing high-stress UO 2 behavior. Behaviors such as this are crucial to capture and implement in BISON, as parts of the fuel pellet can reach temperatures where glide can occur before pellet cracking. This model which captures dislocation based mechanisms for UO 2 is then used to stand up the doped model accounting for larger grain sizes. It was found that larger grain sizes can lead to enhanced deformation rates in the glide regime, and therefore can help with the pellet cladding mechanical interaction. Therefore if the fuel pellet reaches conditions (stress/temperature) where glide is active, the enhanced creep rates for larger grains in the glide regime (doped UO 2 ) can help with pellet cladding mechanical interactions. Plastic deformation in UO 2 involves multiple mechanisms, including diffusional creep, dislocation climb, and glide. This milestone contains two parts: (1) UQ of a pre-existing lower length scale informed mechanistic diffusional creep model, and (2) development of a new LApx based model for dislocation-mediated creep mechanisms in UO 2 , with application to large-grain doped UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Separate Effects Validation for Subcooled Boiling in CTF

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which is focused on developing high-fidelity and high-resolution modeling and simulation (M&S) tools for modeling of light water reactors (LWRs). This hub, the Consortium for Advanced Simulation of LWRs (CASL), has developed an LWR simulation tool called the Virtual Environment for Reactor Applications (VERA). The multi-physics capability of VERA is achieved through the coupling of single-physics codes, including CTF (the CASL version of Coolant Boiling in Rod Arrays— Three Field (COBRA-TF)), Michigan Parallel Characteristics Transport (MPACT), BISON, and Materials Performance and Optimization (MPO) Advanced Model for Boron Analysis (MAMBA). As part of its M&S efforts, CASL has identified various challenge problems, including Crud Induced Power Shift (CIPS), Crud-Induced Localized Corrosion (CILC), Pellet-Cladding Interaction (PCI), and Departure from Nucleate Boiling (DNB). This work addresses CASL milestone L2:VVI.P19.03, which focuses on uncertainty quantification of crud, which is relevant to both CIPS and CILC. This is achieved through an analysis and separate effects validation of the thermal hydraulic phenomenon known as subcooled boiling. As part of this work, various sources of experimental data are examined and compared to different options for empirical modeling of subcooled boiling. Through this analysis, a complete understanding of the underlying models and their implementation details are understood. A subset of these data are incorporated into a separate effects validation study of CTF. The Westinghouse Advanced Loop Tester (WALT) and Rohsenow experiments are modeled, and it is shown that the newly-implemented Gorenflo correlation is more accurate than the existing Chen and Thom correlations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A high‐order discontinuous Galerkin approach for physics‐based thermospheric modeling

Abstract The accurate prediction of aerodynamic drag on satellites orbiting in the upper atmosphere is critical to the operational success of modern space technologies, such as satellite‐based communication or navigation systems, which have become increasingly popular in the last few years due to the deployment of constellations of satellites in low‐Earth orbit. As a result, physics‐based models of the ionosphere and thermosphere have emerged as a necessary tool for the prediction of atmospheric outputs under highly variable space weather conditions. This paper proposes a high‐fidelity approach for physics‐based space weather modeling based on the solution of the Navier–Stokes equations using a high‐order discontinuous Galerkin method, combined with a matrix‐free strategy suitable for high‐performance computing on GPU architectures. The approach consists of a thermospheric model that describes a chemically frozen neutral atmosphere in nonhydrostatic equilibrium driven by the external excitation of the Sun. A novel set of variables is considered to treat the low densities present in the upper atmosphere and to accommodate the wide range of scales present in the problem. At the same time, and unlike most existing approaches, radial and angular directions are treated in a nonsegregated approach. The study presents a set of numerical examples that demonstrate the accuracy of the approximation and validate the current approach against observational data along a satellite orbit, including estimates of established empirical and physics‐based models of the ionosphere‐thermosphere system. Finally, a one‐dimensional radial derivation of the physics‐based model is presented and utilized for conducting a parametric study of the main thermal quantities under various solar conditions.

Engineering↗

Quantifying Epistemic Uncertainty in Binary Classification via Accuracy Gain

ABSTRACT Recently, a surge of interest has been given to quantifying epistemic uncertainty (EU), the reducible portion of uncertainty due to lack of data. We propose a novel EU estimator in the binary classification setting, as the posterior expected value of the empirical gain in accuracy between the current prediction and the optimal prediction. In order to validate the performance of our EU estimator, we introduce an experimental procedure where we take an existing dataset, remove a set of points, and compare the estimated EU with the observed change in accuracy. Through real and simulated data experiments, we demonstrate the effectiveness of our proposed EU estimator.

97 MATHEMATICS AND COMPUTING↗

Evaluation of semi-empirical soot models for nonpremixed flames with increased stoichiometric mixture fraction and strain

Due to the size and complexity of industrial-scale systems and computing limitations, semi-empirical soot models are often employed in CFD simulations rather than computationally-expensive detailed models. Many of these models were developed for specific applications with unique characteristic timescales and/or validated only under fuel-air combustion conditions. Hence, their use in different contexts, such as oxy-combustion, could lead to inaccurate predictions. In this study, twelve semi-empirical models (1-step or 2-step) are evaluated on their ability to respond to changes in stoichiometric mixture fraction (Z st ) and strain in a series of ethylene counterflow flames which span the experimental sooting-to-non-sooting (yellow to blue) transition. A unique approach of plotting soot formation rate in normalized, local equivalence ratio space is introduced to aid in the analysis. Results show that no existing model is able to predict a blue flame when Z st is increased beyond the experimentally-measured sooting limit. Many models give the counter-intuitive result of increasing peak and/or integrated soot volume fraction (svf) as Z st increased, contrary to experimental observations. Two models are able to successfully predict a blue flame upon increase of strain to the sooting limit value. For two-step models, the strong dependence of growth rate on surface area results in high sensitivity of svf to changing flame boundary conditions, but with mixed accuracy. There is still a significant need for a robust semi-empirical model which can accurately predict soot fraction in systems with oxygen-enrichment or variable strain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Imposing correct jellium response is key to predict the density response by orbital-free DFT

Orbital-free density functional theory constitutes a computationally highly effective tool for modeling electronic structures of systems ranging from room-temperature materials to warm dense matter. Its accuracy critically depends on the employed kinetic energy (KE) density functional, which has to be supplied as an external input. In this work we consider several nonlocal and Laplacian-level KE functionals and use an external harmonic perturbation to compute the static density response at T = 0 K in the linear and beyond-linear response regimes. We test for the satisfaction of exact conditions in the limit of uniform densities and for how approximate KE functionals reproduce the density response of realistic materials (e.g., Al and Si) against the Kohn-Sham DFT reference, which employs the exact KE. The results illustrate that several functionals violate exact conditions in the uniform electron gas (UEG) limit. We find a strong correlation between the accuracy of the KE functionals in the UEG limit and in the strongly inhomogeneous case. This empirically demonstrates the importance of imposing the limit of UEG response for uniform densities and validates the use of the Lindhard function in the formulation of kernels for nonlocal functionals. This conclusion is substantiated by additional calculations for bulk aluminum (Al) with a face-centered cubic (fcc) lattice and silicon (Si) with an fcc lattice, body-centered cubic (bcc) lattice, and semiconducting crystal diamond state. The analysis of fcc Al, and fcc as well as bcc Si data follows closely the conclusions drawn for the UEG, allowing us to extend our conclusions to realistic systems that are subject to density inhomogeneities induced by ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First Simultaneous Measurement of Differential Muon-Neutrino Charged-Current Cross Sections on Argon for Final States with and without Protons Using MicroBooNE Data

We report the first double-differential neutrino-argon cross section measurement made simultaneously for final states with and without protons for the inclusive muon neutrino charged-current interaction channel. The proton kinematics of this channel are further explored with a differential cross section measurement as a function of the leading proton’s kinetic energy that extends across the detection threshold. These measurements use data collected with the MicroBooNE detector from 6.4 × 10 20 protons on target from the Fermilab booster neutrino beam with a mean neutrino energy of ∼ 0.8 GeV . Extensive data-driven model validation utilizing the conditional constraint formalism is employed. This motivates enlarging the uncertainties with an empirical reweighting approach to minimize the possibility of extracting biased cross section results. The extracted nominal flux-averaged cross sections are compared to widely used event generator predictions revealing severe mismodeling of final states without protons for muon neutrino charged-current interactions, possibly from insufficient treatment of final state interactions. These measurements provide a wealth of new information useful for improving event generators which will enhance the sensitivity of precision measurements in neutrino experiments. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Validation Metrics for Fixed Effects and Mixed-Effects Calibration

The modern scientific process often involves the development of a predictive computational model. To improve its accuracy, a computational model can be calibrated to a set of experimental data. A variety of validation metrics can be used to quantify this process. Some of these metrics have direct physical interpretations and a history of use, while others, especially those for probabilistic data, are more difficult to interpret. In this work, a variety of validation metrics are used to quantify the accuracy of different calibration methods. Frequentist and Bayesian perspectives are used with both fixed effects and mixed-effects statistical models. Through a quantitative comparison of the resulting distributions, the most accurate calibration method can be selected. Two examples are included which compare the results of various validation metrics for different calibration methods. It is quantitatively shown that, in the presence of significant laboratory biases, a fixed effects calibration is significantly less accurate than a mixed-effects calibration. This is because the mixed-effects statistical model better characterizes the underlying parameter distributions than the fixed effects model. The results suggest that validation metrics can be used to select the most accurate calibration model for a particular empirical model with corresponding experimental data.

97 MATHEMATICS AND COMPUTING↗

EMPIRE-PIC: A Performance Portable Unstructured Particle-in-Cell Code

In this study we introduce EMPIRE-PIC, a finite element method particle-in-cell (FEM-PIC) application developed at Sandia National Laboratories. The code has been developed in C++ using the Trilinos library and the Kokkos Performance Portability Framework to enable running on multiple modern compute architectures while only requiring maintenance of a single codebase. EMPIRE-PIC is capable of solving both electrostatic and electromagnetic problems in two- and three-dimensions to second-order accuracy in space and time. In this paper we validate the code against three benchmark problems — a simple electron orbit, an electrostatic Langmuir wave, and a transverse electromagnetic wave propagating through a plasma. We demonstrate the performance of EMPIRE-PIC on four different architectures: Intel Haswell CPUs, Intel's Xeon Phi Knights Landing, ARM Thunder-X2 CPUs, and NVIDIA Tesla V100 GPUs attached to IBM POWER9 processors. This analysis demonstrates scalability of the code up to more than two thousand GPUs, and greater than one hundred thousand CPUs.

97 MATHEMATICS AND COMPUTING↗

Validation of MFUEL Metal Fuel Performance Models of SAS4A/SASSYS-1

The fuel characterization models of SAS4A/SASSYS-1 (SAS) have recently been extended to include a new U-Pu-Zr metal fuel model, MFUEL. MFUEL is equipped with mechanistic physics based models to predict the pre-transient characterization and transient response of metal fuel, with emphasis on fuel melting, cladding failure, and the metal fuel’s impact on core reactivity. The MFUEL model will be available in the full version of SAS4A/SASSYS-1 5.7, which is scheduled to be released in June 2023. Fast reactor fuel pins that operated in EBR-II and FFTF with low smear density U-Zr and U-Pu-Zr metal fuels and irradiation resistant ferritic-martensitic cladding showed significant advantages in achieving high burnups and assuring inherent safety characteristics during anticipated transients, design basis events, and beyond design basis events. For safety analysis, a fuel performance model must be able to predict (1) Fuel pin mechanics and compositional and dimensional changes, (2) Clad failure, and (3) Fuel pin thermal resistance. Achieving these high level goals accurately is strongly related to the model performance of individual physical processes taking place within a fuel pin during its lifetime. Metal fuels typically operate above the mid-point temperature of melting during normal, as well as off-normal, conditions. At these elevated temperatures, the availability of thermal activation provides a driving force for various diffusional processes leading to complex phase transformations, micro-structure evolution, significant amounts of fuel swelling, interconnected porosity formation, excessive amounts of fission gas release, and fuel clad chemical interactions. Clad failure in fast reactors primarily occurs as a result of creep rupture augmented by clad wastage formation. The reaction is driven by thermal creep induced dislocation motion, grain boundary cavity nucleation, growth and breakup of grain boundaries. The high level complexity and limited available data requires introducing physics-based modeling approaches to gain extrapolation ability and sensitivity with respect to various conditions. The objective of this report is to perform validation of MFUEL using the experimental data for (1) Normal operation EBR-II fuel behavior, (2) Normal operation PHENIX fuel behavior, (3) HT9 Pressure tube ramp-and-hold creep rupture tests, (4) Whole Pin Furnace (WPF) creep strain, creep rupture and eutectic tests, (5) Fuel Behavior Test Apparatus (FBTA) eutectic tests, and (6) TREAT M5-7 OverPower tests up to clad failure. Section-2 includes a brief description of the MFUEL models. A detailed description of the MFUEL physics-based, semi-empirical models will be presented in the SAS V 5.7 theory manual. Section-3, 4, and 5 describes the validation effort for the pre-transient irradiation, furnace transients, and TREAT M-Series transients, respectively.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Finding Liebig’s law of the minimum

Abstract Liebig’s law of the minimum (LLM) is often used to interpret empirical biological growth data and model multiple substrates co‐limited growth. However, its mechanistic foundation is rarely discussed, even though its validity has been questioned since its introduction in the 1820s. Here we first show that LLM is a crude approximation of the law of mass action, the state of art theory of biochemical reactions, and the LLM model is less accurate than two other approximations of the law of mass action: the synthesizing unit model and the additive model. We corroborate this conclusion using empirical data sets of algae and plants grown under two co‐limiting substrates. Based on our analysis, we show that when growth is modeled directly as a function of substrate uptake, the LLM model improperly restricts the organism to be of fixed elemental stoichiometry, making it incapable of consistently resolving biological adaptation, ecological evolution, and community assembly. When growth is modeled as a function of the cellular nutrient quota, the LLM model may obtain good results at the risk of incorrect model parameters as compared to those inferred from the more accurate synthesizing unit model. However, biogeochemical models that implement these three formulations are needed to evaluate which formulation is acceptably accurate and their impacts on predicted long‐term ecosystem dynamics. In particular, studies are needed that explore the extent to which parameter calibration can rescue model performance when the mechanistic representation of a biogeochemical process is known to be deficient.

54 ENVIRONMENTAL SCIENCES↗

Discovering the 3D Structure and Dynamics of the Sun-Interstellar Medium System on a Global Scale Los Alamos LDRD Report

We discover the 3D plasma structure of the Sun-interstellar medium interaction region, called the heliosheath. This is enabled by (1) the completion of a full solar cycle of observations of hydrogen energetic neutral atoms (ENAs) emanating from this region and detected by the Los Alamos-led IBEX-Hi ENA imager on NASA’s Interstellar Boundary Explorer (IBEX) mission, and (2) the development of tools to fully exploit all data from the IBEX-Hi imager, leading to a significant increase in the effective instrument sensitivity. This effort prepares us for optimizing the design and data exploitation of the IMAP-Hi imager, the Los Alamos-led next generation imager that will fly on the Interstellar Mapping and Acceleration Probe (IMAP) mission, currently under production. This project realized its goals as evidenced by two main outcomes: (1) The publication of the first-ever empirical 3D map of the interstellar boundary, and (2) the development and release of a new validated ENA data set that has tripled the amount of data available for scientific analysis of the outer heliosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Perspective on Traditional and Data Driven Electrochemical Modeling and Analysis

To understand the behavior of electrochemical systems, we need to reduce the dimensionality of the measured current-voltage-time (I-V-t) data by fitting models, thus enabling us to analyze and compare the governing physics. Traditionally, the process for this is an 'expert first' approach: defining the model and its explicit assumptions based on inductive reasoning or empirical observation, fitting small portions of the I-V-t data where assumptions are most valid or carefully designing experiments to enforce key assumptions, and then interpreting the model parameters. However, modern data-driven methods enable a new paradigm: a 'data first' approach, where the latent behaviors governing the system's measured response are identified directly using machine-learning models that optimize both model structure and parameters from the I-V-t data, guaranteeing that the learned model explains as much of the observed system response as possible. After model identification, the model can then be interrogated by an expert to connect observed behaviors with underlying physics. This talk will review several different types of electrochemical analysis (electrochemical impedance, differential voltage-capacity, electrochemical kinetics) and compare the traditional and data-driven methods for analyzing the data.

42 ENGINEERING↗

Machine-Learning Assisted Identification of Battery Life Models

Predictive battery life models are commonly utilized to extrapolate degradation trends observed during accelerated aging tests for simulation of degradation in real-world applications. Thus, fitting accelerated aging data as accurately as possible and with low uncertainty is crucial for making believable projections of battery lifetime, but it is challenging to identify algebraic expressions that accurately fit multivariate degradation trends. A review of models published in literature reveal some common expressions for fitting calendar aging data, which is only dependent on temperature and state-of-charge, but no consistency across many models for fitting cycle aging data, indicating the need for a statistically rigorous data driven approach for developing empirical models. This talk will describe a machine-learning assisted method for identification of predictive battery life models utilizing bilevel optimization and symbolic regression. Bilevel optimization with cross-validation is used to statistically determine cell- and stress-dependent model parameters, while symbolic regression identifies both linear and multiplicative candidate expressions to predict stress-dependent degradation rates by selecting low-order subsets of features from a generated feature library. Because model expressions are identified empirically, it is crucial to ensure resulting models behave according to physical expectations, so the stability of models for interpolation or extrapolation is interrogated qualitatively through simulation and quantitatively through cross-validation and uncertainty quantification via bootstrap resampling. This model identification approach substantially improves upon models identified purely using expert judgement in terms of both accuracy and uncertainty. Model simulation and validation is then conducted by deriving a state-equation form of the predictive model, enabling simulation of battery aging under dynamic stresses. This enables validation of the predictive battery model on lab-based tests with varying conditions or on drive-cycle or application-cycle testing protocols. Parameter uncertainty can be carried forward into model simulation, giving lifetime estimates and confidence windows for cell- or system-level lifetime. The financial impact of battery model uncertainty can be estimated by incorporating uncertainty into a technoeconomic model.

battery↗

Verification, Validation, and Uncertainty Quantification in Thermal Hydraulics, Freeman Scholar Lecture (2019)

Engineering problems are generally solved by analytical models or computer codes. These models, in addition to conservation equations, also include many empirical relationships and approximate numerical methods. Each of these components contributes to the uncertainty in the prediction. A systematic approach to judge the applicability of the code to the intended application is needed. It starts from verification of implementation of formulation in the code, identification of important phenomena, finding relevant tests with quantified uncertainty for these phenomena, and validation of the code by comparing predictions with the relevant test data. The relevant tests must address phenomena as expected in the intended application. In case of small size or limited condition tests, the scaling analyses are needed to assess the relevancy of the tests. Finally, a statement of uncertainty in the prediction is needed. Systematic approaches are described to aggregate uncertainties from different components of the code for intended application. Here, verification, validation, and uncertainty quantifications (VVUQs) are briefly described.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SAM Code Validation on Frictional Pressure Drop through Pebble Beds

The System Analysis Module (SAM) is currently under development at Argonne National Laboratory as a modern system-level modeling and simulation tool for safety analyses of advanced non-light water reactors. This report presents a recent effort to validate the capability of SAM to predict the frictional pressure drop through pebble beds. Selected experimental data were used for code validation, including data from test facilities at Texas A&M University, Missouri University of Science and Technology, and North-West University of South Africa. SAM implements three empirical correlations to predict frictional pressure drop: the classical Ergun correlation; the KTA correlation, which is widely used in high-temperature gas-cooled reactor applications; and the Eisfeld and Schnitzlein correlation, which explicitly considers wall effect. Code validation was performed using all three correlations. For all selected experimental data, the KTA correlation shows the best performance and agrees very well with experimental measurement; the Eisfeld and Schnitzlein correlation, explicitly considering wall effects, shows acceptable accuracy, while there is no evidence that it is better than the KTA correlation; the Ergun correlation, however, over-predicts frictional pressure drop for most selected data points.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗