Search NASA⌕ Search

SEARCH · Search NASA

Results for “Probability estimates”

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 361 records · Page 20

Learning Many-Body Hamiltonians with Heisenberg-Limited Scaling

Learning a many-body Hamiltonian from its dynamics is a fundamental problem in physics. Here, in this Letter, we propose the first algorithm to achieve the Heisenberg limit for learning an interacting N-qubit local Hamiltonian. After a total evolution time of $\mathscr{O}$⁡(ε –1 ), the proposed algorithm can efficiently estimate any parameter in the N-qubit Hamiltonian to ε error with high probability. Our algorithm uses ideas from quantum simulation to decouple the unknown N-qubit Hamiltonian H into noninteracting patches and learns H using a quantum-enhanced divide-and-conquer approach. The proposed algorithm is robust against state preparation and measurement error, does not require eigenstates or thermal states, and only uses polylog⁡(ε –1 ) experiments. In contrast, the best existing algorithms require $\mathscr{O}$⁡⁡(ε –2 ) experiments and total evolution time. We prove a matching lower bound to establish the asymptotic optimality of our algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Induced seismicity and geothermal energy production in the Salton Sea Geothermal Field, California

We analyze the relationship between geothermal energy production and seismic hazards in the Salton Sea Geothermal Field (SSGF) between 1972 and 2022. A clear increase in seismic activity accompanies geothermal energy production and is greatest to the east of the Brawley fault, where the amount of injection exceeds the amount of production. We estimate that, whereas there was a 2% chance of a M6.07 earthquake nucleating inside the SSGF within a fifty-year period prior to energy production (pre-production), there is an equal probability that a M6.70 earthquake nucleates there in a fifty-year period during which historically average energy production conditions persist (syn-production). Similarly, we estimate a 2% chance of a M7.15 and M6.92 occurring in the broader Brawley Seismic Zone (excluding the SSGF) within a fifty-year period during the pre- and syn-production eras, respectively. We find that linear regression models fail to reliably forecast the background seismicity rate as a function of fluid production and injection.

15 GEOTHERMAL ENERGY↗

Quantifying Uncertainty in PV Energy Estimates Final Report

Uncertainty in PV energy estimates is "one of the most critical areas of lack of understanding" according to independent engineers, financiers, PV model developers, and other industry stakeholders. The primary problem is a lack of rigorous, transparent, widely accepted methods for quantifying uncertainty in energy production estimates. Uncertainty in energy production estimates arises from variability of the solar resource, inexact PV performance models and their parameters, and system reliability considerations. Uncertainty in annual energy production is frequently calculated for larger projects in order to quantify financial risk. Key statistics for energy, such as the P-values "P50" and "P90" (the annual energy values that are exceeded in future years with 50\% and 90\% probability, respectively) are used by financing institutions to calculate the repayment risk for the project. The current methods to estimate these statistics are typically proprietary, specialized, and involve significant post-processing of commercial performance model results. This black-box approach leads to inconsistent P-value estimates from different parties, which reduces investors' confidence in the results. Since the financial community bases its risk assessment on these estimates, reduced confidence increases perceived project risk, and consequently financing costs. The goal of this project was to establish a set of best practices for quantifying uncertainty in energy production estimates, including identifying what sources of uncertainty must be considered with clear definitions and metrics, determining which sources are the biggest drivers of uncertainty, and providing a computationally efficient framework for combining different sources of uncertainty that is flexible enough to accommodate substitutions of data or methods when better information is available. We engaged a wide set of stakeholders to ensure industry endorsement and adoption, and leveraged complementary projects investigating individual sources of uncertainty in great detail, as well as others' work that started down this path.

ENERGY PLANNING, POLICY, AND ECONOMY,SOLAR ENERGY↗

Uncertainty-aware Continuous Implicit Neural Representations for Remote Sensing Object Counting

Many existing object counting methods rely on density map estimation (DME) of the discrete grid representation by decoding extracted image semantic features from designed convolutional neural networks (CNNs). Relying on discrete density maps not only leads to information loss dependent on the original image resolution, but also has a scalability issue when analyzing high-resolution images with cubically increasing memory complexity. Furthermore, none of the existing methods can offer reliable uncertainty quantification (UQ) for the derived count estimates. To overcome these limitations, we design UNcertainty-aware, hypernetwork-based Implicit neural representations for Counting (UNIC) to assign probabilities and the corresponding counting confidence over continuous spatial coordinates. We derive a sampling-based Bayesian counting loss function and develop the corresponding model training algorithm. UNIC outperforms existing methods on the Remote Sensing Object Counting (RSOC) dataset with reliable UQ and improved interpretability of the derived count estimates. Our code is available at https://github.com/SiyuanXu-tamu/UNIC.

97 MATHEMATICS AND COMPUTING↗

Surrogate Modeling For Efficiently Accurately and Conservatively Estimating Measures of Risk

We present a surrogate modeling framework for conservatively estimating measures of risk from limited realizations of an expensive physical experiment or computational simulation. We adopt a probabilistic description of risk that assigns probabilities to consequences associated with an event and use risk measures, which combine objective evidence with the subjective values of decision makers, to quantify anticipated outcomes. Given a set of samples, we construct a surrogate model that produces estimates of risk measures that are always greater than their empirical estimates obtained from the training data. These surrogate models not only limit over-confidence in reliability and safety assessments, but produce estimates of risk measures that converge much faster to the true value than purely sample-based estimates. We first detail the construction of conservative surrogate models that can be tailored to the specific risk preferences of the stakeholder and then present an approach, based upon stochastic orders, for constructing surrogate models that are conservative with respect to families of risk measures. The surrogate models introduce a bias that allows them to conservatively estimate the target risk measures. We provide theoretical results that show that this bias decays at the same rate as the L 2 error in the surrogate model. Our numerical examples confirm that risk-aware surrogate models do indeed over-estimate the target risk measures while converging at the expected rate.

97 MATHEMATICS AND COMPUTING↗

Formal Experimentation and Analysis of Handheld RFID Readers as a Tool for Nuclear Material Accounting

Commercial, off-the-shelf Radio Frequency Identification (RFID) systems have been successfully deployed for inventory tracking in numerous industries, but their viability in the tracking of complex environments containing nuclear material is less understood. Of primary interest in this setting is a RFID tracking system for nuclear material accounting which can reliably identify as many tags as possible with high accuracy, while also reducing an operator’s exposure to radiation. In this work, we develop a formal statistical approach to identify relevant handheld RFID reader settings which optimize tagging performance. To achieve this goal, we design a full factorial split-plot experiment for a static shelf configuration scene with 50 randomly placed nuclear material containers affixed with RFID tags. We use Bayesian inference to fit a second-order response surface model which expresses the probability of a successful match for each container as a function of the experimental factors. Such effects are allowed to vary across individual containers and the containers’ population in its entirety to estimate overall effects. Uncertainties of estimates and predictions are quantified via their corresponding posterior distributions. Following extensive model checking and validation, the fitted model is used to identify experimental factors which maximize matching probabilities at both the container-level and for the full shelf configuration scene. We also analyze sensitivity of performance to relevant factors.

25 ENERGY STORAGE↗

Bottomonium suppression and elliptic flow using Heavy Quarkonium Quantum Dynamics

We introduce a framework called Heavy Quarkonium Quantum Dynamics (HQQD) which can be used to compute the dynamical suppression of heavy quarkonia propagating in the quark-gluon plasma using real-time in-medium quantum evolution. Using HQQD we compute large sets of real-time solutions to the Schrödinger equation using a realistic in-medium complex-valued potential. We sample 2 million quarkonia wave packet trajectories and evolve them through the QGP using HQQD to obtain their survival probabilities. The computation is performed using three different HQQD model parameter sets in order to estimate our systematic uncertainty. After taking into account final state feed down we compare our results to existing experimental data for the suppression and elliptic flow of bottomonium states and find that HQQD predictions are good agreement with available data for R AA as a function of N part and p T collected at \( \sqrt{s_{\mathrm{NN}}} \) = 5 . 02 TeV. In the case of v 2 for the various states, we find that the path-length dependence of ϒ(1 s ) suppression results in quite small v 2 for ϒ(1 s ). Our prediction for the integrated elliptic flow for ϒ(1 s ) in the 10 – 90% centrality class, which now includes an estimate of the systematic error, is v 2 [ϒ(1 s )] = 0 . 003 ± 0 . 0007 ± \( {}_{0.0013}^{0.0006} \) . We also find that, due to their increased suppression, excited bottomonium states have a larger elliptic flow. Based on this observation we make predictions for v 2 [ϒ(2 s )] and v 2 [ϒ(3 s )] as a function of centrality and transverse momentum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fluorine Limits and Impacts in High-Level Waste Glass Compositions

The impact of elevated fluorine (F) content on Hanford high-level waste (HLW) glasses has not previously been studied in detail. This effort represents the first systematic study to determine what F concentration limits should be used for the design of alkali-borosilicate-based Hanford Waste Treatment and Immobilization Plant (WTP) HLW glasses, and to document the technical basis for that limit. If alkali borosilicate glass made from Hanford HLW can accommodate a large amount of F, the large capital costs for complex sludge washing facilities may be avoided, as would much of the operational costs and negative schedule impacts associated with handling the large volumes of water required to dissolve these salts. In order to determine a limit for F in likely HLW glass compositions, an evaluation was conducted on glasses with F ≤ 0.90 mass% from previous nuclear waste glass studies. The collected dataset contains 239 glasses (232 HLW glasses and 7 LAW glasses) including 109 glasses with 0.9 ≤ F mass% ≤ 2.5, 116 with 2.5 < F mass% ≤ 8.0, and 14 with F mass% ≥ 8 (max. F mass% = 17.42). The collected composition and property data were analyzed to determine the basis for the F tolerance, i.e. the maximum F concentration that can be processed without potential issues. Fluorine volatility, product consistency test (PCT) response, liquidus temperature (T L ), glass melt viscosity, and crystallinity have been investigated. No limits for F concentration can be made based on F volatility, T L , or glass melt viscosity, because the data show that high F in glasses do not indicate, with high probability, being restricted by those property constrains. However, crystallinity and PCT response were used to estimate the F tolerance. The results show that glasses with high F (≥ 0.90 mass%) are more likely to form large fractions of F-containing crystal phases which may increase PCT responses, i.e. decrease the glass durability. Based on the results of crystallinity and PCT data, the F tolerance of 4.5 mass% is estimated. There is no evidence of other glass components, such as calcium oxides and alkali metal oxides have combined impacts with F on the glass properties. Overall, the available high-F glass data is limited, especially in the designed HLW glass composition regions. Future work on formulation and testing of HLW glasses with F ≥ 0.9 mass% will close the data gaps and expand operational flexibility with respect to the fluoride tolerances. Volatility of F from melters and corrosion of materials in contact with glass melts are important for processing of high-F wastes; yet no test data are currently available. It is recommended tests be conducted to address these two potential issues.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Spatial patterns in occupancy and density of larval lampreys in freshwater habitats restored to a Stage 0 condition

Abstract We examined occupancy and density of larval lampreys ( Entosphenus tridentatus and Lampetra spp.) in two rivers in Oregon (USA) restored to a Stage 0 condition 1–5 years prior, using a multiscale occupancy model and a zero‐inflated Poisson mixture model. We sampled lampreys using backpack electrofishing in randomly distributed, paired, 1‐m 2 quadrats and recorded environmental data. Probabilities of occupancy and density were higher when water velocity was low, the substrate was noncompacted, and sediment was dominated by fines (<4 mm). At mean water depth (0.34 m) and velocity (0.09 m/s), estimated densities in occupied quadrats were 4.8 lampreys/m 2 (95%: 3.4–6.9) when the substrate was compacted, and fines were not dominant, and 21.1 lampreys/m 2 (95%: 17.7–25.3) when the substrate was noncompacted and fines were dominant. Probabilities of detecting occupancy in a 1‐m 2 quadrat sampled by backpack electrofishing were 0.76 (95%: 0.64–0.87) when captured after visual observation and 0.80 (95%: 0.71–0.88) with blind sweeps (i.e., constantly moving the net regardless of observation). The probability of capturing a single lamprey in a quadrat sampled by blind sweeps was 0.32 (95%: 0.27–0.37). Sampling in paired 1‐m 2 quadrats facilitated concurrent examination of patterns in occupancy and density while accounting for capture probability, which could aid temporal monitoring of restored habitats. To the best of our knowledge, this is the first study to document occupancy and estimate densities of larval lampreys in habitats that underwent valley floor restoration to Stage 0. We observed both lamprey genera within 5 years of restoration. Aquatic restoration that increases low‐velocity, noncompacted, fine sediment habitats could benefit lampreys.

Harris, Julianne E.↗

Variational control forces for enhanced sampling of nonequilibrium molecular dynamics simulations

We introduce a variational algorithm to estimate the likelihood of a rare event within a nonequilibrium molecular dynamics simulation through the evaluation of an optimal control force. Optimization of a control force within a chosen basis is made possible by explicit forms for the gradients of a cost function in terms of the susceptibility of driven trajectories to changes in variational parameters. We consider probabilities of time-integrated dynamical observables as characterized by their large deviation functions and find that in many cases, the variational estimate is quantitatively accurate. Additionally, we provide expressions to exactly correct the variational estimate that can be evaluated directly. We benchmark this algorithm against the numerically exact solution of a model of a driven particle in a periodic potential, where the control force can be represented with a complete basis. Additionally, we then demonstrate the utility of the algorithm in a model of repulsive particles on a line, which undergo a dynamical phase transition, resulting in singular changes to the form of the optimal control force. In both systems, we find fast convergence and are able to evaluate large deviation functions with significant increases in statistical efficiency over alternative Monte Carlo approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reduced Order Modeling conditioned on monitored features for response and error bounds estimation in engineered systems

Reduced Order Models (ROMs) form essential tools across engineering domains by virtue of their function as surrogates for computationally intensive digital twinning simulators. Although purely data-driven methods are available for ROM construction, schemes that allow to retain a portion of the physics tend to enhance the interpretability and generalization of ROMs. However, physics-based techniques can adversely scale when dealing with nonlinear systems that feature parametric dependencies. This study introduces a generative physics-based ROM that is suited for nonlinear systems with parametric dependencies and is additionally able to provide numerical error bounds associated with the respective estimates. A main contribution of this work is the conditioning of these parametric ROMs to features that can be derived from monitoring measurements, feasibly in an online fashion. This is contrary to most existing ROM schemes, which remain restricted to the prescription of the physics-based, and usually a priori unknown, system parameters. Our work utilizes conditional Variational Autoencoders to continuously map the required reduction bases to a feature vector extracted from limited output measurements, while additionally allowing for a probabilistic assessment of the ROM-estimated Quantities of Interest. An auxiliary task using a neural network-based parametrization of suitable probability distributions is introduced to re-establish the link with physical model parameters. We verify the proposed scheme on a series of simulated case studies incorporating effects of geometric and material nonlinearity under parametric dependencies related to system properties and input load characteristics.

Conditional VAEs↗

Quantifying concordance of correlated cosmological data sets

We develop estimators of agreement and disagreement between correlated cosmological data sets. These account for data correlations when computing the significance of both tensions and excess confirmation while remaining statistically optimal. We discuss and thoroughly characterize different approaches commenting on the ones that have the best behavior in practical applications. We complement the calculation of their statistical distribution within the Gaussian model with one estimator that takes non-Gaussianities fully into account. To illustrate the use of our techniques, we apply these estimators to supernovae measurements of the distance-redshift relation, absolutely calibrated by the local distance ladder. The suite of best estimators that we discuss finds results that are in excellent agreement between estimators and find no indications of significant internal inconsistencies in this data set above the 1% probability threshold. Furthermore, this shows the robustness of local determinations of the Hubble constant to features in the distance-redshift relation.

79 ASTRONOMY AND ASTROPHYSICS↗

Adapting CLUTCH methodology to multigroup TSUNAMI-3D for eigenvalue sensitivity calculations

The sensitivity of the eigenvalue to uncertainties in nuclear data and its evaluation are important for nuclear criticality safety. TSUNAMI-3D sequences within the SCALE code system offer several options to the user community for calculating eigenvalue sensitivity coefficients with multigroup (MG) and continuous energy (CE) 3D transport capabilities. TSUNAMI-3D sequences implement the adjoint-based perturbation theory with MG KENO code, the Contributon Linked eigenvalue sensitivity/Uncertainty estimation via Track length importance CHaracterization (CLUTCH) method with CE KENO code, and the Iterated Fission Probability (IFP) method with CE KENO and Shift codes. Each method has benefits and limitations depending on the problem that is run. The work presented here aims to adapt the CLUTCH method, which enables the Contributon method's mesh-free, memory-efficient approach for calculating adjoint-weighted tallies for sensitivity calculations, to the MG TSUNAMI-3D sequence. This application would eliminate the explicit adjoint KENO calculation, as well as the memory-consuming mesh flux moment tallies required by the conventional MG TSUNAMI-3D. Smaller memory footprints in the CLUTCH methodology and relatively shorter runtimes in MG KENO transport can make MG TSUNAMI-3D a viable method for some complex problems. Moreover, this adaptation allows MG sensitivity calculations with Shift, ORNL's next-generation high-performance Monte Carlo transport code, which currently does not offer any sensitivity capabilities with MG particle transport simulations. Initial implementation of the new MG TSUNAMI-3D sequence and its preliminary results with a selected critical benchmark experiment in the Verified, Archived Library of Inputs and Data (VALID) are presented in this study.

KENO↗

Impact of Increased Latent Generations on Sensitivity Calculations with SCALE [Slides]

Depending on the system, energy, and/or nuclide, altering the number of latent generations can significantly impact generated sensitivity coefficients. In general, as the number of latent generations increases, the accuracy of the generated sensitivity value compared with the DP value increases and there is a significant increase in the uncertainty with generated sensitivity values. IFP generated sensitivities appear to be more stable than those with CLUTCH. Complex systems, even with IFP and increased latent generations can still produce poor results (i.e., MPC- 32). This reiterates the importance of performing DPs to confirm sensitivities. IFP-Shift and CLUTCH can take advantage of parallel computing abilities. Work continues in this area as additional parameters and calculational methods are examined to provide analysts insights to successfully generating sensitivity coefficients for confirmatory analyses and validation efforts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measuring Young Stars in Space and Time. II. The Pre-main-sequence Stellar Content of N44

The Hubble Space Telescope survey Measuring Young Stars in Space and Time (MYSST) entails some of the deepest photometric observations of extragalactic star formation, capturing even the lowest-mass stars of the active star-forming complex N44 in the Large Magellanic Cloud. We employ the new MYSST stellar catalog to identify and characterize the content of young pre-main-sequence (PMS) stars across N44 and analyze the PMS clustering structure. To distinguish PMS stars from more evolved line of sight contaminants, a non-trivial task due to several effects that alter photometry, we utilize a machine-learning classification approach. This consists of training a support vector machine (SVM) and a random forest (RF) on a carefully selected subset of the MYSST data and categorize all observed stars as PMS or non-PMS. Combining SVM and RF predictions to retrieve the most robust set of PMS sources, we find ∼26,700 candidates with a PMS probability above 95% across N44. Employing a clustering approach based on a nearest neighbor surface density estimate, we identify 16 prominent PMS structures at 1σ significance above the mean density with sub-clusters persisting up to and beyond 3σ significance. The most active star-forming center, located at the western edge of N44's bubble, is a subcluster with an effective radius of ∼5.6 pc entailing more than 1100 PMS candidates. Furthermore, we confirm that almost all identified clusters coincide with known H ii regions and are close to or harbor massive young O stars or YSOs previously discovered by MUSE and Spitzer observations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The anomaly of the CMB power with the latest Planck data

Abstract The lack of power anomaly is an unexpected feature observed at large angular scales in the maps of Cosmic Microwave Background (CMB) produced by the COBE, WMAP andPlancksatellites. This signature, which consists in a missing of power with respect to that predicted by the ΛCDM model, might hint at a new cosmological phase before the standard inflationary era.The main point of this paper is taking into account the latestPlanckpolarisation data to investigate how the CMB polarisation improves the understanding of this feature. With this aim, we apply to the latestPlanckdata, both PR3 (2018) and PR4 (2020) releases, a new class of estimators capable of evaluating this anomaly by considering temperature and polarisation data both separately and in a jointly way. This is the first time that the PR4 dataset has been used to study this anomaly. To critically evaluate this feature, taking into account the residuals of known systematic effects present in thePlanckdatasets, we analyse the cleaned CMB maps using different combinations of sky masks, harmonic range and binning on the CMB multipoles.Our analysis shows that the estimator based only on temperature data confirms the presence of a lack of power with a lower-tail-probability (LTP), depending on the component separation method, ≤ 0.33% and ≤ 1.76% for PR3 and PR4, respectively. To our knowledge, the LTP≤ 0.33% for the PR3 dataset is the lowest one present in the literature obtained fromPlanck2018 data, considering thePlanckconfidence mask. We find significant differences between these two datasets when polarisation is taken into account most likely due to a different level of systematics. Especially, the analysis with PR3 data, unlike that with PR4, seems to point towards a lack of power at large scales also for polarisation.Moreover, we also show that for the PR3 dataset the inclusion of the subdominant polarisation information provides estimates that are less likely accepted in a ΛCDM cosmological model than the only-temperature analysis over the entire harmonic-range considered. In particular, at ℓ max = 26, we found that no simulation has a value as low as the data for all the pipelines.

Astronomy & Astrophysics↗

Research Trends and Applications of PMUs

This work is a survey of current trends in applications of PMUs. PMUs have the potential to solve major problems in the areas of power system estimation, protection, and stability. A variety of methods are being used for these purposes, including statistical techniques, mathematical transformations, probability, and AI. The results produced by the techniques reviewed in this work are promising, but there is work to be performed in the context of implementation and standardization. As the smart grid initiative continues to advance, the number of intelligent devices monitoring the power grid continues to increase. PMUs are at the center of this initiative, and as a result, each year more PMUs are deployed across the grid. Since their introduction, myriad solutions based on PMU-technology have been suggested. The high sampling rates and synchronized measurements provided by PMUs are expected to drive significant advancements across multiple fields, such as the protection, estimation, and control of the power grid. This work offers a review of contemporary research trends and applications of PMU technology. Most solutions presented in this work were published in the last five years, and techniques showing potential for significant impact are highlighted in greater detail. Being a relatively new technology, there are several issues that must be addressed before PMU-based solutions can be successfully implemented. This survey found that key areas where improvements are needed include the establishment of PMU-observability, data processing algorithms, the handling of heterogeneous sampling rates, and the minimization of the investment in infrastructure for PMU communication. Solutions based on Bayesian estimation, as well as those having a distributed architectures, show great promise. The material presented in this document is tailored to both new researchers entering this field and experienced researchers wishing to become acquainted with emerging trends.

42 ENGINEERING↗