Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Best Estimate Plus Uncertainty”

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.

MatLib-1.2.2: Nuclear Material Properties Library

The U.S. Nuclear Regulatory Commission (NRC) uses the computer code Fuel Analysis under Steady-state and Transients (FAST) to model steady-state and transient fuel behavior to support regulatory decisions. FAST relies on a material properties library (MatLib) that contains the thermal and mechanical properties of the nuclear materials and coolants of interest to support the U.S. commercial nuclear industry. MatLib contains properties for a variety of nuclear fuels, cladding and other structural materials, gases, and coolants. In this document, material property correlations for the materials contained within MatLib are presented and discussed. When available, comparisons are made between the material property correlations and available data. Additionally, uncertainties are quantified on the material properties, which is then used by the NRC to support uncertainty quantification for best-estimate plus uncertainty safety evaluation reviews. This document describes MatLib-1.2.2, updated from MatLib-1.2.1 to include additional properties for metallic fuel. It is one of a series of documents on FAST; the other documents detail the models used by FAST as well as its integral assessment to experiments and commercial data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Hierarchical Bayesian modeling for Inverse Uncertainty Quantification of system thermal-hydraulics code using critical flow experimental data

The best estimate plus uncertainty methodology in nuclear system thermal-hydraulic studies necessitates a comprehensive understanding of uncertainties in system code predictions. The forward uncertainty quantification (UQ) process involves the propagation of input uncertainties through the computational models to obtain uncertainties in the outputs. To this end, achieving an accurate estimation of input uncertainties is important, which is the focus of inverse UQ (IUQ). Traditionally, research in Bayesian IUQ within the nuclear engineering domain has largely relied on single-level Bayesian inference. While being effective for relatively small datasets, this approach encounters limitations for cases with large datasets. The use of a single-level model may prove inefficient, as the resultant posterior distributions can significantly differ when distinct subsets of data are employed. To address this issue, we employ an hierarchical Bayesian model for IUQ. Furthermore, this approach involves organizing observations into different groups based on the test conditions, thereby accommodating varying calibration parameters across these distinct groups. In this study, we developed and implemented a hierarchical Bayesian IUQ method to consider the grouping effect of critical flow measurement data from various geometries. Comparing the outcomes of IUQ under different selections of test data using hierarchical Bayesian IUQ against those obtained from single-level Bayesian IUQ, the forward propagation of hierarchical Bayesian IUQ results demonstrates a notably improved agreement with the experimental data.

42 ENGINEERING

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion

Development and application of two-step uncertainty propagation and sensitivity analysis methodology for fast reactor safety analysis

Uncertainty quantification (UQ) in nuclear reactors for transients is directly linked with safety assessment through the cross-sections uncertainties, provided as a covariance matrix, which are propagated through the reactor system to output of interest pertaining to reactor safety, such as peak temperatures in fuel/clad/coolant. Using a two-step approach, uncertainties are first quantified and propagated from basic input variables (such as reaction cross-sections) to intermediate quantities (such as reactivity feedback coefficients) through lattice level calculations. Uncertainties of intermediate quantities (from the first step) are then propagated through the system transient calculations, in the second step, to obtain uncertainties on reactor safety output parameters of interest. The scope of this work consists of Uncertainty Quantification & Propagation of nuclear data uncertainties that are highly correlated through unprotected transient overpower and unprotected loss of flow to assess their impact on core safety parameters. This two-step approach in the presence of covariance renders the sensitivity analysis very challenging. In fact, usually the sensitivity analysis is restricted to each step, which limits its application since the sensitivities between the system output quantities and the basic input variables are difficult to obtain. Here, in this work, we address this issue by proposing a simple, general methodology to combine the sensitivity indices obtained in each step by assuming the model behavior being linear. For the first step Generalized Perturbation theory based indices are used while in the second step the recently studied Johnson indices. The uncertainty quantification and sensitivity methodologies discussed here are demonstrated on a generic LFR design which is based on the 500 MWth demonstration Lead-cooled fast reactor (DLFR) using oxide fuel, developed by Westinghouse Electric Company (WEC).

42 - ENGINEERING

Robust calibration of non-perturbative models with history matching

We apply, for the first time, Bayes Linear Emulation and History Matching to the calibration of non-perturbative models in Monte Carlo event generators. In contrast to the usual approach of "Monte Carlo tuning", History Matching does not result in best-fit plus ellipsoidal parameter uncertainty estimates but instead identifies all parameter space regions that are consistent with data. This approach leads to a systematic and robust quantification of parametric uncertainties in the models, especially in those challenging cases where different, possibly disjoint, regions of parameter space deliver similar results, which are usually not properly treated with current methodology. We highlight the power of this method with the hadronisation models available through Sherpa : the built-in cluster fragmentation Ahadic and string fragmentation through an interface to Pythia .

Iskauskas, Andrew [Durham U. (main)] (ORCID:000000

SeaWiFS Technical Report Series: SeaWiFS Calibration Topics - Volume 40

For Earth-observing satellite instruments, it was standard to consider each instrument band to have a spectral response that is infinitely narrow, i.e., to have a response from a single wavelength. The SeaWiFS bands, however, have nominal spectral bandwidths of 20 and 40 nm. These bandwidths effect the SeaWiFS measurements on orbit. The effects are also linked to the manner in which the instrument was calibrated and to the spectral shape of the radiance that SeaWiFS views. The spectral shape of that radiance will not be well known on orbit. In this technical memorandum, two source spectra are examined. The first is a 12,000 K Planck function, and the second is based on the modeling results of H. Gordon at the University of Miami. By comparing these spectra, the best available corrections to the SeaWiFS measurements for source spectral shape, plus estimates of the uncertainties in these corrections, can be tabulated.

Barnes, Robert A.

All-Sky Observational Evidence for An Inverse Correlation Between Dust Temperature and Emissivity Spectral Index

We show that a one-component variable-emissivity-spectral-index model (the free- model) provides more physically motivated estimates of dust temperature at the Galactic polar caps than one- or two-component fixed-emissivity-spectral-index models (fixed- models) for interstellar dust thermal emission at far-infrared and millimeter wavelengths. For the comparison we have fit all-sky one-component dust models with fixed or variable emissivity spectral index to a new and improved version of the 210-channel dust spectra from the COBE-FIRAS, the 100-240 micrometer maps from the COBE-DIRBE and the 94 GHz dust map from the WMAP. The best model, the free-alpha model, is well constrained by data at 60-3000 GHz over 86 per cent of the total sky area. It predicts dust temperature (T(sub dust)) to be 13.7-22.7 (plus or minus 1.3) K, the emissivity spectral index (alpha) to be 1.2-3.1 (plus or minus 0.3) and the optical depth (tau) to range 0.6-46 x 10(exp -5) with a 23 per cent uncertainty. Using these estimates, we present all-sky evidence for an inverse correlation between the emissivity spectral index and dust temperature, which fits the relation alpha = 1/(delta + omega (raised dot) T(sub dust) with delta = -.0.510 plus or minus 0.011 and omega = 0.059 plus or minus 0.001. This best model will be useful to cosmic microwave background experiments for removing foreground dust contamination and it can serve as an all-sky extended-frequency reference for future higher resolution dust models.

Liang, Z.

Atmospheric opacity in the Schumann-Runge bands and the aeronomic dissociation of water vapor

Knowledge of the aeronomic production of odd hydrogen in the dissociation of water vapor is limited by uncertainties in the penetration of solar irradiance in the Schumann-Runge bands of O2 and by incomplete information concerning the products of photolysis at Lyman alpha. Consideration of all error sources involved in computing the H2O dissociation rate in the wavelength region 175-200 nm leads to an estimated uncertainty of plus or minus 35% at an altitude of 90 km for an overhead sun. The uncertainty increases with decreasing altitude such that the true dissociation rate at 60 km for an overhead sun lies between 0.45 and 1.55 times the results computed using the best input parameters currently available. Calculations of the H2O dissociation rate by Lyman alpha should include the variation in O2 opacity across the solar line width. Neglect of this can lead to errors as large as 50% at altitudes where the process is the major source of odd hydrogen.

Frederick, J. E.

Effective radiative forcing and adjustments in CMIP6 models

The effective radiative forcing, which includes the instantaneous forcing plus adjustments from the atmosphere and surface, has emerged as the key metric of evaluating human and natural influence on the climate. We evaluate effective radiative forcing and adjustments in 17 contemporary climate models that are participating in the Coupled Model Intercomparison Project (CMIP6) and have contributed to the Radiative Forcing Model Intercomparison Project (RFMIP). Present-day (2014) global-mean anthropogenic forcing relative to pre-industrial (1850) levels from climate models stands at 2.00 (±0.23) W/sq. m, comprised of 1.81 (±0.09) W/sq. m from CO2, 1.08 (± 0.21) W/sq. m from other well-mixed greenhouse gases, −1.01 (± 0.23) W/sq. m from aerosols and −0.09 (±0.13) W/sq. m from land use change. Quoted uncertainties are 1 standard deviation across model best estimates, and 90 % confidence in the reported forcings, due to internal variability, is typically within 0.1 W/sq. m. The majority of the remaining 0.21 W/sq. m is likely to be from ozone. In most cases, the largest contributors to the spread in effective radiative forcing (ERF) is from the instantaneous radiative forcing (IRF) and from cloud responses, particularly aerosol–cloud interactions to aerosol forcing. As determined in previous studies, cancellation of tropospheric and surface adjustments means that the stratospherically adjusted radiative forcing is approximately equal to ERF for greenhouse gas forcing but not for aerosols, and consequentially, not for the anthropogenic total. The spread of aerosol forcing ranges from −0.63 to −1.37 W/sq. m, exhibiting a less negative mean and narrower range compared to 10 CMIP5 models. The spread in 4×CO2 forcing has also narrowed in CMIP6 compared to 13 CMIP5 models. Aerosol forcing is uncorrelated with climate sensitivity. Therefore, there is no evidence to suggest that the increasing spread in climate sensitivity in CMIP6 models, particularly related to high-sensitivity models, is a consequence of a stronger negative present-day aerosol forcing and little evidence that modelling groups are systematically tuning climate sensitivity or aerosol forcing to recreate observed historical warming.

effective radiative forcing

Updating PV and Battery Bill Savings Calculations for Net Billing: New Best Practices for Input Data and Uncertainty

Jurisdictions are increasingly adopting compensation structures for distributed PV and PV-battery systems that price exported energy lower than energy consumed onsite (net billing rates). Standard methods for calculating the bill savings from PV and PV-battery systems were developed for net metering structures, and applying these same methods to net billing structures (such as using Typical Meteorological Year weather with actual year load) introduces bias errors that underestimate PV-battery system bill savings by between 1.5\% and 9\%, depending on the utility rate. We assess the magnitude of these errors and compare them to other sources of uncertainty when estimating the bill savings from PV and PV-battery systems under more complex utility rates.

14 SOLAR ENERGY