LANL Contribution to NRAP-Open-IAM
NRAP-Open-IAM is an open-source software product that enables quantification of containment effectiveness and leakage risk at storage sites in the context of system uncertainties and variability.
Chen, Bailian↗
SEARCH · Engineering Papers
Results for “uncertainty Quantification”
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.
Compaction of crushed salt for safe containment – a summary of the KOMPASS projects
Abstract. For the underground disposal of high-level nuclear waste in rock salt formations, the safety concept includes the backfilling of open cavities with crushed salt. For the prognosis of the sealing function of the backfill for the safe containment of the radioactive waste, it is crucial to have a comprehensive process understanding of the crushed-salt compaction behavior. The crushed-salt compaction process is influenced by internal properties (e.g., grain size, mineralogy, and moisture content) and boundary conditions (e.g., temperature, stress state, and compaction rate) and, therefore, involves several coupled thermal–hydro–mechanical (THM) processes (Hansen et al., 2014; Kröhn et al., 2017). With the paradigm shift from the limited release of radionuclides to safe containment due to the German Repository Site Selection Act passed in 2017, the importance of crushed salt as geotechnical barrier has increased, with a focus on the evolution of its hydraulic properties. Based on the knowledge gaps in the current process understanding, the “Compaction of crushed salt for safe containment” (KOMPASS) projects were initiated to improve the scientific basis behind using crushed salt for the long-term isolation of high-level nuclear waste within rock salt repositories. The efforts to improve the prediction of crushed-salt compaction begun during the first phase of the KOMPASS projects (Czaikowski et al., 2020) and were followed up in a second phase ending in June 2023. The primary achievements of the projects are as follows (Czaikowski et al., 2020; Friedenberg et al., 2022): specification of the KOMPASS reference material, an easily available and reproducible synthetic crushed-salt material, for generic investigations; development of pre-compaction methods and successful production of samples in the short term and under in situ loading conditions; formulation of an extended laboratory program addressing the isolated investigation of known relevant factors influencing the compaction behavior of crushed salt (Düsterloh et al., 2022); execution of long-term compaction tests addressing isotropic and deviatoric load changes, temperature, and compaction state; construction of a backfill body using the KOMPASS reference material in the Sondershausen mine through collaboration with the SAVER (Entwicklung eines salzgrusbasierten Versatzkonzepts unter der Option Rückholbarkeit) project (Schaarschmidt and Friedenberg, 2022); advancement of the tools for microstructure investigation methods (Svensson and Laurich, 2022); generation (first stages) of a microphysical process list combining literature research with our own findings; benchmarking of long-term compaction test for model development and optimization of various existing models as well as the development of new models; application of a virtual demonstrator (2D model representing a backfilled drift in rock salt) for the visualization of developments and the quantification of the models (Rabbel, 2022). In summary, the KOMPASS projects contributed to the reduction of uncertainties and the strengthening of the safety case for using crushed salt within rock salt repositories.
Friedenberg, Larissa↗
Quantifying uncertainty in analysis of shockless dynamic compression experiments on platinum. II. Bayesian model calibration
Dynamic shockless compression experiments provide the ability to explore material behavior at extreme pressures but relatively low temperatures. Typically, the data from these types of experiments are interpreted through an analytic method called Lagrangian analysis. Here, in this work, alternative analysis methods are explored using modern statistical methods. Specifically, Bayesian model calibration is applied to a new set of platinum data shocklessly compressed to 570 GPa. Several platinum equation-of-state models are evaluated, including traditional parametric forms as well as a novel non-parametric model concept. The results are compared to those in Paper I obtained by inverse Lagrangian analysis. The comparisons suggest that Bayesian calibration is not only a viable framework for precise quantification of the compression path, but also reveals insights pertaining to trade-offs surrounding model form selection, sensitivities of the relevant experimental uncertainties, and assumptions and limitations within Lagrangian analysis. The non-parametric model method, in particular, is found to give precise unbiased results and is expected to be useful over a wide range of applications. The calibration results in estimates of the platinum principal isentrope over the full range of experimental pressures to a standard error of 1.6%, which extends the results from Paper I while maintaining the high precision required for the platinum pressure standard.
Brown, Justin Lee↗
An adaptive synchronous extraction (ASE) method for estimating intensity and footprint of surface urban heat islands: A case study of 254 North American cities
The urban heat island (UHI) effect has attracted great attention due to its potential impacts on rapidly growing urban areas. Using remotely sensed estimates of land surface temperature (LST), a large number of studies have focused on the surface UHI (SUHI) effect, which can be characterized by its two fundamental properties: intensity and footprint. The SUHI intensity reflects the LST difference between the urban area and the background reference area (BRA), and the SUHI footprint indicates the spatial extent influenced by the heat island. Currently, numerous methods have been developed to estimate the SUHI intensity and footprint, but are still greatly challenged by three main issues. Namely, the discrepancy in BRA selection criterion brings great uncertainty to the estimated SUHI intensity, the estimation of SUHI footprint is largely constrained by the predefined models, and the quantification of SUHI effect is potentially influenced by several confounding factors. Here, we proposed an adaptive synchronous extraction (ASE) method, which is capable of adaptively selecting the most optimal BRA while removing the influence of confounding factors, and achieving synchronous estimation of SUHI intensity and footprint. We applied the ASE method to 254 North American cities and conducted an in-depth comparative analysis to discuss its applicability and benefits. The main results include: (1) The ASE method avoids the limitations of existing methods in BRA selection and model presetting, and shows resilience to parameter variations. This makes the ASE method highly applicable to quantify the SUHI intensity and footprint in cities with various thermal characteristics. (2) The ASE method can better highlight the spatial, seasonal and day-night contrasts in the estimated SUHI intensity. This superiority is particularly evident when comparing it to methods based on the equal-area buffer or the simplified urban-extent algorithm. (3) Confounding factors pose non-negligible impacts on the quantification of the SUHI effect. Typically, ignoring the influence of topographic relief or missing LST data can lead to an overall overestimation of the SUHI intensity, while not removing surrounding urban areas will cause some underestimation of the SUHI intensity. In conclusion, overall, the proposed ASE method provides a new generalizable tool for quantifying the SUHI effect, which has great potentials for future studies and urban climate assessments.
54 ENVIRONMENTAL SCIENCES↗
Albedo susceptibility of northeastern Pacific stratocumulus: the role of covarying meteorological conditions
Abstract. Quantification of the radiative adjustment of marine low clouds to aerosol perturbations, regionally and globally, remains the largest source of uncertainty in assessing current and future climate. One of the important steps towards quantifying the role of aerosol in modifying cloud radiative properties is to quantify the susceptibility of cloud albedo and liquid water path (LWP) to perturbations in cloud droplet number concentration (Nd). We use 10 years of spaceborne observations from the polar-orbiting Aqua satellite to quantify the albedo susceptibility of marine low clouds to Nd perturbations over the northeast (NE) Pacific stratocumulus (Sc) region. Mutual information analysis reveals a dominating control of cloud state (e.g., LWP and Nd) on low-cloud albedo susceptibility, relative to the meteorological states that drive these cloud states. Through a LWP–Nd space decomposition of albedo susceptibilities, we show clear separation among susceptibility regimes (brightening or darkening), consistent with previously established mechanisms through which aerosol modulates cloud properties. These regimes include (i) thin non-precipitating clouds (LWP < 55 g m−2) that exhibit brightening (occurring 37 % of the time), corresponding to the Twomey effect; (ii) thicker non-precipitating clouds, corresponding to entrainment-driven negative LWP adjustments that manifest as a darkening regime (36 % of the time); and (iii) another brightening regime (22 % of the time) consisting of mostly precipitating clouds, corresponding to precipitation-suppression LWP positive adjustments. Overall, we find an annual-mean regional low-cloud brightening potential of 20.8±2.68 W m−2 ln(Nd)−1, despite an overall negative LWP adjustment for non-precipitating marine stratocumulus, owing to the high occurrence of the Twomey–brightening regime. Over the NE Pacific, clear seasonal covariabilities among meteorological factors related to the large-scale circulation are found to play an important role in grouping conditions favorable for each susceptibility regime. When considering the covarying meteorological conditions, our results indicate that for the northeastern Pacific stratocumulus, clouds that exhibit the strongest brightening potential occur most frequently within shallow marine boundary layers over a cool ocean surface with a stable atmosphere and a dry free troposphere above. Clouds that exhibit a darkening potential associated with negative LWP adjustments occur most frequently within deep marine boundary layers in which the atmospheric instability and the ocean surface are not strong and warm enough to produce frequent precipitation. Cloud brightening associated with warm-rain suppression is found to preferably occur either under unstable atmospheric conditions or humid free-tropospheric conditions that co-occur with a warm ocean surface.
54 ENVIRONMENTAL SCIENCES↗
Quantification of human contribution to soil moisture-based terrestrial aridity
Abstract Current knowledge of the spatiotemporal patterns of changes in soil moisture-based terrestrial aridity has considerable uncertainty. Using Standardized Soil Moisture Index (SSI) calculated from multi-source merged data sets, we find widespread drying in the global midlatitudes, and wetting in the northern subtropics and in spring between 45°N–65°N, during 1971–2016. Formal detection and attribution analysis shows that human forcings, especially greenhouse gases, contribute significantly to the changes in 0–10 cm SSI during August–November, and 0–100 cm during September–April. We further develop and apply an emergent constraint method on the future SSI’s signal-to-noise (S/N) ratios and trends under the Shared Socioeconomic Pathway 5-8.5. The results show continued significant presence of human forcings and more rapid drying in 0–10 cm than 0–100 cm. Our findings highlight the predominant human contributions to spatiotemporally heterogenous terrestrial aridification, providing a basis for drought and flood risk management.
54 ENVIRONMENTAL SCIENCES↗
Assessing the nexus between groundwater and solar energy plants in a desert basin with a dual-model approach under uncertainty
Globally, many solar power plants and other types of renewable energy are being located in water-scarce regions. Many projects rely on groundwater resources whose sustainability is uncertain. In the Chuckwalla Basin in California, quantification of recharge and trans-valley underflow is needed to estimate the impacts of solar project withdrawals on the water table. However, such estimates are highly challenging due to data scarcity, heterogeneous soils and long residence times. Conventional assessment employs isolated groundwater models configured with crude and uniform estimates of recharge. Here, we employ a data-constrained surface subsurface processes model, PAWS+CLM, to provide an ensemble of recharges and underflows with perturbed parameters. Then, the Parameter Estimation (PEST) package is used to calibrate MODFLOW aquifer conductivity and filter out implausible recharges. The novel dual-model approach, potentially applicable in other arid regions, can effectively assimilate groundwater head observations, reject unrealistic parameters, and narrow the range of estimated drawdowns. Simulated recharge concentrates along alluvial fans at the mountain fronts and ephemeral washes where run-off water infiltrates. If an evenly distributed recharge was assumed, it resulted in under-estimated drawdown and larger uncertainty bounds. The withdrawals are approaching total inflow, suggesting the system will be nearing, if not exceeding, its sustainable groundwater production capacity, and a boom of such projects will not be sustainable. Especially, the cost/benefit of pumped-storage projects is called into question as the initial-fill phase depletes entire area’s recharge. Our study highlights the stress on groundwater resources of solar development, and that the speed of groundwater recovery does not indicate sustainability. Main point 1: A novel dual model approach, involving an integrated surface/subsurface model and a groundwater parameter-estimation model, was able to better constrain the model. Main point 2: The groundwater system may be nearing, if not exceeding, its sustainable groundwater production capacity and the speed of recovery is not indicative of sustainability. Main point 3: Results from using conventionally-assumed uniform recharge distort calibrated K fields and impacts assessment
14 SOLAR ENERGY↗
Impact of Wildfires on Solar Generation, Reserves and Energy Prices
Wildfire seasons in the Western U.S. become more prolonged and intense in recent years bringing significant variability and uncertainty of solar generation. It greatly challenges the bulk power system and electricity market operation managed by California Independent System Operator (CAISO) as California leads both the wildfire records and the solar power integration. This study presents a screening-level analysis of the impact of wildfires on solar generation, operating reserve and energy prices applying historical real-world wildfire and market operation data. To the best of the authors knowledge, it is a first-of-its-kind study and will lay the foundation for market impact quantification and wildfire mitigation strategies design based on projected wildfire activities in future years.
electricity price↗
Uranium Holdup Monitoring with Compton Imaging as Function of Depth and Mass
Material Control and Accountability (MC&A) programs track masses of materials flowing through production processes, where some material remains in process equipment as holdup. Quantification of these masses using gamma emissions typically involves usage of a well characterized radiation detector and assumptions of geometric configuration to constrain detection efficiencies, where the Generalized Geometry Holdup (GGH) method is often used. Uncertainties on the order of ±50% are observed, where the difficulties arise based on the asymmetric and gradient distributions of material that do not fall within the bounds of GGH. Gamma imagers provide a possible means to improve this geometric knowledge. This work describes the process used in creating a series of mass standards to aid the development of imaging algorithms for quantification.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Assessing arsenic species in foods using regularized linear regression of the arsenic K-edge X-ray absorption near edge structure
Solvent-free quantifications of arsenic species in foods are attained by applying lasso regression to the analysis of X-ray absorption spectra and assigning uncertainties by bootstrapping.
Chemistry↗
LA17-V-InverseRadTransportProblem-PD2Kb. Final Project Report
In the inverse radiation transport problem, a radiation signal emanating from an unknown object is used to reconstruct geometric and material parameters of the object , e.g., the mass of special nuclear material (SNM), the isotopic composition of the SNM (e.g., the enrichment of 235U), the thickness of shielding, etc. The inverse radiation transport problem is one of the most fundamental and important in the emergency response mission, and has been addressed in previous NA-22 and NA-42/84 investments by using derivative-free optimization methods. Our project was aimed at (i) consolidating previous work by applying polynomial chaos expansion (PCE) for reducing significantly the number of forward model simulations required by derivative-free optimization methods; (ii) developing a general-purpose framework and stand-alone software for solving the inverse radiation transport problem by applying a new “predictive modeling” methodology (predictive modeling of coupled multi-physics systems, PM_CMPS) which enables the fusing of experimental and computational information; and (iii) developing several benchmarks (for unscattered gamma rays) for demonstrating the quantification of non-Gaussian features (e.g., skewness and kurtosis, which characterize asymmetries and long tails) of response distributions by deterministically and efficiently propagating parameter uncertainties (covariances and higher-order moments, where available) using high-order sensitivities computed by the recently developed comprehensive adjoint sensitivity analysis methodology (C-ASAM). With regard to (i), application of PCE, the most significant accomplishment of this project was an implementation of PCE into INVERSE that significantly reduced computation time for the DiffeRential Evolution Adaptive Metropolis (DREAM) and Differential Evolution (DE) methods. For inverse transport problems with measurements of neutron leakage, neutron multiplication, and neutron-induced gamma-rays, the solution times for the DREAM and DE methods were reduced by factors of 30-300, reducing run times from hours to minutes while maintaining the same accuracy as the standard (non-PCE) DREAM and DE. Speedups were further improved by exploring sparse grid quadrature sets using the Toolkit for Adaptive Stochastic Modeling and Non-Intrusive ApproximatioN (Tasmanian) toolset. An initial exploration of sparse grids indicated that further speedups factors of 5-10 are obtainable. The implementation of PCE in INVERSE is not in a form that users can easily access. With regard to (ii), solving the inverse problem using a predictive modeling paradigm, the most significant accomplishment of this project were the derivation and implementation of adjoint-based second derivatives in neutron and photon transport theory. There are two codes that compute these derivatives. SENSPG computes them for uncollided gamma rays and SENSMG computes them for neutron leakage. SENSMG is open-source code available on github. Other significant accomplishments included a stand-alone code, MULTI-PRED, that implements the predictive modeling paradigm, along with a detailed manual. The code was implemented in the existing INVERSE code for inverse problems, but a head-to-head comparison of inverse methods was not achieved. We also derived and implemented the capability to compute sensitivities for multiplicity measurements into the SENSMG code, and this capability was used in a data assimilation application, but it was not used for solving inverse problems. With regard to (iii), developing benchmarks to demonstrate the importance of high-order sensitivity and uncertainty analysis, the most significant accomplishment of this project were two test problems that showed how important second-order sensitivities are.
07 ISOTOPE AND RADIATION SOURCES↗
Quantification of regional net CO 2 flux errors in the Orbiting Carbon Observatory-2 (OCO-2) v10 model intercomparison project (MIP) ensemble using airborne measurements
Inverse model intercomparison projects (MIPs) provide a chance to assess the uncertainties in inversion estimates arising from various sources. However, accurately quantifying ensemble CO 2 flux errors remains challenging and often relies on the ensemble spread. This study proposes a method for quantifying the errors in regional net surface–atmosphere CO 2 flux estimates from models taken from the Orbiting Carbon Observatory-2 (OCO-2) v10 MIP by using independent airborne CO 2 measurements for the period 2015–2017. We first calculate the root mean square error (RMSE) between the ensemble mean of posterior CO 2 concentrations and airborne observations and then isolate the CO 2 concentration errors caused solely by the ensemble mean of posterior net fluxes by subtracting the observation, representation, and transport errors from seven regions. Our analysis reveals that the flux errors projected onto CO 2 space account for 55 %–85 % of the regional average RMSE over the 3 years, ranging from 0.88 to 1.91 ppm. In five regions, the error estimates based on observations exceed those computed from the ensemble spread of posterior fluxes by a factor of 1.33–1.93, implying an underestimation of the actual flux errors, while their magnitudes are comparable in two regions. The adjoint sensitivity analysis identifies that the underestimation of flux errors is prominent where the magnitudes of fossil fuel emissions exceed those of terrestrial-biosphere fluxes by a factor of 3–31 over the 3 years. This suggests the presence of systematic biases in the inversion estimates associated with errors in the prescribed fossil fuel emissions common to all models. Our study emphasizes the value of airborne measurements for quantifying regional errors in ensemble net CO 2 flux estimates.
54 ENVIRONMENTAL SCIENCES↗
Open World Dempster-Shafer Theory/The Transferable Belief Model with Intervals: A Practitioner's Guide to DST and TBM
Dempster-Shafer theory (DST) is a mathematical framework that allows for uncertainty or ignorance to be quantified and included when making predictions from evidence. This is in contrast to Bayesian theory, which does not allow for any quantification of ignorance. The framework is described in great detail in [7]. DST is particularly useful for problems where the inclusion of additional evidence (for example, data from another sensor) could lead to a different conclusion. Thus, it is a useful data fusion method, especially in applications not suited to maximum likelihood or maximum a posteriori estimations due to limited samples or incomplete prior knowledge.
97 MATHEMATICS AND COMPUTING↗
Modeling biological nitrogen fixation in global natural terrestrial ecosystems
Biological nitrogen fixation plays an important role in the global nitrogen cycle. However, the fixation rate has been usually measured or estimated at a particular observational site. To quantify the fixation amount at the global scale, process-based models are needed. This study develops a biological nitrogen fixation model to quantitatively estimate the nitrogen fixation rate by plants in a natural environment. The revised nitrogen module better simulates the nitrogen cycle in comparison with our previous model that has not considered the fixation effects. The new model estimates that tropical forests have the highest fixation rate among all ecosystem types, which decreases from the Equator to the polar region. The estimated nitrogen fixation in global terrestrial ecosystems is 61.5 Tg N yr-1 with a range of 19.8–107.9 Tg N yr-1 in the 1990s. Our estimates are relatively low compared to some early estimates using empirical approaches but comparable to more recent estimates that involve more detailed processes in their modeling. Furthermore, the contribution of nitrogen made by biological nitrogen fixation depends on ecosystem type and climatic conditions. This study highlights that there are relatively large effects of biological nitrogen fixation on ecosystem nitrogen cycling. and the large uncertainty of the estimation calls for more comprehensive understanding of biological nitrogen fixation. More direct observational data for different ecosystems are in need to improve future quantification of fixation and its impacts.
54 ENVIRONMENTAL SCIENCES↗
Restoration Strategy for Active Distribution Systems Considering Endogenous Uncertainty in Cold Load Pickup
Cold load pickup (CLPU) phenomenon is identified as the persistent power inrush upon a sudden load pickup after an outage. Under the active distribution system (ADS) paradigm, where distributed energy resources (DERs) are extensively installed, the decreased outage duration can induce a strong interdependence between CLPU pattern and load pickup decisions. In this paper, we propose a novel modelling technique to tractably capture the decision-dependent uncertainty (DDU) inherent in the CLPU process. Subsequently, a two-stage stochastic decision-dependent service restoration (SDDSR) model is constructed, where first stage searches for the optimal switching sequences to decide step-wise network topology, and the second stage optimizes the detailed generation schedule of DERs as well as the energization of switchable loads. Further, to tackle the computational burdens introduced by mixed-integer recourse, the progressive hedging algorithm (PHA) is utilized to decompose the original model into scenario-wise subproblems that can be solved in parallel. The numerical test on modified IEEE 123-node test feeders has verified the efficiency of our proposed SDDSR model and provided fresh insights into the monetary and secure values of DDU quantification.
24 POWER TRANSMISSION AND DISTRIBUTION↗
PVEMC: isolating the flavor-dependent EMC effect using parity-violating inelastic scattering in SoLID
In order to better understand the EMC effect, we propose a clean and precise measurement of the flavor dependence of the EMC effect using parity-violating deep inelastic scattering on a 48 Ca target. This measurement will provide an extremely sensitive test for flavor dependence in the modification of nuclear parton distribution functions (PDFs) for neutron-rich nuclei. A measurement of the flavor dependence will provide new and vital information and help to explain nucleon modification at the quark level. In addition to helping understand the origin of the EMC effect, a flavor-dependent nuclear pdf modification could significantly impact a range of processes, including neutrino-nucleus scattering, nuclear Drell-Yan processes, and e-A observables at the Electron-Ion Collider. The parity-violating asymmetry A PV from 48 Ca using an 11 GeV beam at 80 μA will be measured using the SoLID detector, proposed for a series of measurements in Hall A at Jefferson Lab. In 68 days of data taking, we will reach 0.7–1.3% statistical precision for 0.2 < x < 0.7 with 0.6–0.7 systematic uncertainties. Here, the goal is to make the first direct measurement of the flavor dependence of the EMC effect. The precision of the measurement will allow for quantification of the flavor-dependent effects, greatly improving our ability to differentiate between models of the EMC effect and constraining the u- and d-quark contributions in neutron rich nuclei.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Quantification of 238U Holdup using Iterative, Point-Cloud-Based Compton Imaging
U holdup has impacts on several aspects of operation: Worker dose, Criticality safety, Safeguards, Outage planning. Generalized Geometry Holdup (GGH) currently estimates U mass within a high uncertainty band (± 50%). This method seeks to improve this uncertainty via quantitative Compton imaging.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
An improved method for uranium separation from urine for rapid bioassay by inductively coupled plasma mass spectrometry
Bioassay monitoring is an important capability for any facility that works with significant amounts of radioactive material. The speed at which these analyses can be made is critical in the event of a potential exposure. Here, we present a new method that facilitates the rapid separation of uranium from a urine matrix for assay and isotopic determination by inductively coupled plasma mass spectrometry. This method is validated with archived urine samples that were previously analyzed using established, more time consuming methods. The sources of uncertainty for uranium assay and isotopic measurements are modelled using the Monte Carlo method to better understand the precision and accuracy of this method at the limit of quantification.
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗