Search NASA⌕ Search

SEARCH · Search NASA

Results for “UNCERTAINTY PRINCIPLE”

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 199 records · Page 11

Performing k eff Validation of As-Loaded Criticality Safety Calculations Using UNF-ST&DARDS: Applicable Experiment Selection

The general method for performing validation of as loaded criticality safety calculations using UNF ST&DARDS is presented in a paper by Clarity, which includes a description of the UNF-ST&DARDS system. Proof-of-principle analyses were performed in the summer of 2019 for MPC-32 dual purpose canisters (DPCs) containing pressurized water reactor (PWR) fuel assemblies. Summaries of these results are presented in this and a companion paper for this conference. The current paper describes the TSUNAMI-IP calculations performed to select applicable experiments for validation of 11 MPC-32 DPCs. The companion paper discusses the TSUNAMI-3D calculations used to generate sensitivity data to support the experiment selections discussed here. Experiment selection is based on the sensitivity/uncertainty (S/U) methods used to validate criticality safety calculations of as-loaded DPCs containing pressurized water reactor (PWR) spent nuclear fuel (SNF). This process has been demonstrated and is summarized in this paper. The approach is similar to that used in NUREG/CR-7109, which provides an approach for validation of PWR burnup credit (BUC), including major and minor actinides and major fission products. The premise of S/U-based validation is that applicable experiments—those having a similar bias to a given application system—will have similar sensitivities for each isotope and reaction in the two systems. It is assumed that cross sections with larger uncertainties are more likely to contain data errors which contribute to the bias. The integral index c k thus propagates the system sensitivities with the nuclear covariance data to calculate a correlation coefficient representing the similarity of the two systems. In this work, a c k value of 0.8 or higher is interpreted as identifying an experiment with sufficient similarity for use in validation. This paper presents a brief summary of the characteristics of the 11 MPC-32 DPCs used in the proof of-principle analysis for as-loaded criticality safety calculation validation and an overview of the critical experiment suite with which each of these DPC models was compared. A summary and discussion of c k results is also presented, followed by conclusions and a discussion of future work to be performed for validation of UNF ST&DARDS as-loaded criticality safety calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dual-frequency (Ka-band and G-band) radar estimates of liquid water content profiles in shallow clouds

The profile of the liquid water content (LWC) in clouds provides fundamental information for understanding the internal structure of clouds, their radiative effects, propensity to precipitate, and degree of entrainment and mixing with the surrounding environment. In principle, differential absorption techniques based on coincident dual-frequency radar reflectivity observations have the potential to provide the LWC profile. Previous differential frequency radar reflectivity (DFR) efforts were challenged by the fact that the measurable differential attenuation for small quantities of LWC is usually comparable to the system measurement error. This typically renders the retrieval impractical, as the uncertainty can become many times greater than the retrieved value itself. Theoretically, this drawback can be mitigated following two interconnected approaches: (1) increasing the frequency separation between the dual-frequency radar system to measure greater differential attenuation and (2) increasing the radar operating frequency to reduce the instrument measurement random error. Our recently developed 239 GHz radar was deployed during the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) along with a variety of collocated remote sensing and in situ instruments. We have combined Ka-band (35 GHz) and G-band (239 GHz) observations to retrieve the LWC from more than 100 vertical profiles of shallow clouds with typical amounts of LWC smaller than 1 g m -3 . We theoretically and experimentally demonstrate that the Ka-band and G-band pair of frequencies offers at least a 65 % relative improvement in the LWC retrieval sensitivity compared to previous works reported in the literature using lower-frequency radars. This new technique provides a missing capability to determine the LWC in the challenging low liquid water path (LWP) range (< 200 g m -2 ) and suggests a way forward to characterize microphysical and dynamical processes more precisely in shallow clouds.

54 ENVIRONMENTAL SCIENCES↗

Suppression of antiferromagnetic order and strong ferromagnetic spin fluctuations in Ca(Co 1-x Ni x ) 2-y As 2 single crystals

CaCo 2–y As 2 is a unique itinerant system having strong magnetic frustration. Here, we report the effect of electron doping on the physical properties resulting from Ni substitutions for Co. The single crystals of Ca(Co 1–x Ni x ) 2–y As 2 were characterized by single-crystal x-ray diffraction, energy-dispersive x-ray spectroscopy, magnetization M versus temperature T, magnetic field H, time t, and heat capacity C p (H, T) measurements. The A-type antiferromagnetic (AFM) transition temperature T N = 52 K for x = 0 decreases to 22 K with only 3% Ni substitution and is completely suppressed for x > 0.16. For 0.11 ≤ x ≤ 0.52 strong ferromagnetic (FM) fluctuations develop as revealed by magnetic susceptibility χ(T) = M(T) / H measurements. For x = 0.11 and 0.16 competing AFM and FM interactions result in a reentrant spin-glass behavior below T N , as evidenced by the observations of thermomagnetic hysteresis and magnetic relaxation. Enhanced FM fluctuations are also found for the x = 0.21 and 0.31 crystals, where χ c increases significantly at low T. A large χ anisotropy in these compositions where χ c is up to a factor of two larger than χ ab suggests that the FM spin fluctuations are quasi-1D in nature. Weak ferromagnetic contributions to the magnetization are found at T = 2 K for x = 0.11 –0.31. Heat-capacity Cp(T) measurements reveal the presence of FM quantum spin fluctuations for 0.11 ≤ x ≤ 0.52, where a logarithmic T dependence of Cp(T) / T is observed at low T. Here, the suppression of AFM order by the development of strong FM fluctuations in Ca(Co 1–x Ni x ) 2–y As 2 crystals suggests the presence of a FM quantum-critical point at x ≈ 0.20. Our density-functional theory (DFT) calculations confirm that FM fluctuations are enhanced by Ni substitutions for Co in CaCo 2–y As 2 . The Sommerfeld electronic heat-capacity coefficient is enhanced for x = 0, 0.21, and 0.42 by about a factor of two compared to DFT calculations of the density of states (DOS) at the Fermi energy, suggesting an enhancement of the DOS from electron-phonon and/or electron-electron interactions. The crystals with x > 0.52 do not exhibit FM spin fluctuations or magnetic order at T ≥ 1.8 K, which was found from the DFT calculations to arise from a Stoner transition. Superconductivity is not observed above 1.8 K for any of the compositions. Neutron-diffraction studies of crystals with x = 0.11 and 0.16 in the crossover regime (0.1 ≲ x ≲ 0.2) show no evidence of A-type ordering as observed in the parent compound with x = 0. Furthermore, no other common magnetic structures, such as ferromagnetic (FM), helical stacking of in-plane FM layers, or in-plane AFM structure, are found with an ordered moment greater than the uncertainty of 0.05 μ B per transition-metal atom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Development of a health monitoring framework: Application to a supercritical pulverized coal-fired boiler

In this study, this work details the development of a physics-based equipment health monitoring framework for a supercritical boiler, using first-principles models to estimate the remaining useful life (RUL) of its components. The framework accounts for fatigue and creep life consumption, generating spatio-temporal variations in mechanical and thermal stress. Analysis of the stress profile throughout the boiler highlights the finishing superheater inlet steam header as a vulnerable location susceptible to damage from cycling operation. The framework also yields quantified uncertainty in the RUL projection for specific locations, accounting for uncertainties in material properties and boiler operation. Results indicate that operational uncertainties (e.g., seasonal variation and operational strategy) and material properties (e.g., rupture time coefficients, Young’s modulus, yield strength, and coefficient of thermal expansion) significantly impact the RUL of the finishing superheater inlet steam header. Additionally, case studies demonstrate the use of the health monitoring framework as a predictive tool for operational planning under uncertainty, including scenarios with and without updates on the operation of the boiler.

20 FOSSIL-FUELED POWER PLANTS↗

First-principles calculations of metal surfaces. I. Slab-consistent bulk reference for convergent surface properties

The first-principles computation of the surfaces of metals is typically accomplished through slab calculations of finite thickness. The extraction of a convergent surface formation energy from slab calculations is dependent upon defining an appropriate bulk reference energy. I describe a method for an independently computed, slab-consistent bulk reference that leads to convergent surface formation energies from slab calculations that also provides realistic uncertainties for the magnitude of unavoidable nonlinear divergence in the surface formation energy with slab thickness. In conclusion, the accuracy is demonstrated on relaxed, unreconstructed low-index aluminum surfaces with slabs with up to 35 layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Data-driven recovery of hidden physics in reduced order modeling of fluid flows

In this article, we introduce a modular hybrid analysis and modeling (HAM) approach to account for hidden physics in reduced order modeling (ROM) of parameterized systems relevant to fluid dynamics. The hybrid ROM framework is based on using first principles to model the known physics in conjunction with utilizing the data-driven machine learning tools to model the remaining residual that is hidden in data. This framework employs proper orthogonal decomposition as a compression tool to construct orthonormal bases and Galerkin projection (GP) as a model to build the dynamical core of the system. Our proposed methodology hence compensates structural or epistemic uncertainties in models and utilizes the observed data snapshots to compute true modal coefficients spanned by these bases. The GP model is then corrected at every time step with a data-driven rectification using a long short-term memory (LSTM) neural network architecture to incorporate hidden physics. A Grassmannian manifold approach is also adopted for interpolating basis functions to unseen parametric conditions. The control parameter governing the system's behavior is thus implicitly considered through true modal coefficients as input features to the LSTM network. The effectiveness of the HAM approach is then discussed through illustrative examples that are generated synthetically to take hidden physics into account. Furthermore, our approach thus provides insights addressing a fundamental limitation of the physics-based models when the governing equations are incomplete to represent underlying physical processes.

42 ENGINEERING↗

Complete collision data set for electrons scattering on molecular hydrogen and its isotopologues: III. Vibrational excitation via electronic excitation and radiative decay

We present state-resolved cross sections for vibrational excitation via electronic excitation followed by radiative decay (ERD), for electrons scattering on all bound vibrational levels of the ground electronic state ($X\hspace{1.5 mm} ^1Σ_g^+$) of molecular hydrogen and its isotopologues (H 2 , HD, HT, D 2 , DT and T 2 ). We consider excitation of the singlet $\textit{n}$ = 2–3 states (where n refers to the united-atoms-limit principle quantum number) and account for all possible decay pathways back to the bound vibrational levels of the ground electronic state ($X\hspace{1.5 mm} ^1Σ_g^+$) to produce an estimate for the ERD cross sections. A selection of the results are presented in graphical form and a full data set is provided as both numerical values and analytic fit functions in supplementary data files. The uncertainty in the cross sections is estimated to be 11%, except for scattering on the highest bound vibrational level of HD, HT, D 2 , and DT, where the uncertainty is 21%. The data can be downloaded from the MCCC database at mccc-db.org.

74 ATOMIC AND MOLECULAR PHYSICS↗

Probabilistic and maximum entropy modeling of chemical reaction systems: Characteristics and comparisons to mass action kinetic models

We demonstrate and characterize a first-principles approach to modeling the mass action dynamics of metabolism. Starting from a basic definition of entropy expressed as a multinomial probability density using Boltzmann probabilities with standard chemical potentials, we derive and compare the free energy dissipation and the entropy production rates. We express the relation between entropy production and the chemical master equation for modeling metabolism, which unifies chemical kinetics and chemical thermodynamics. Because prediction uncertainty with respect to parameter variability is frequently a concern with mass action models utilizing rate constants, we compare and contrast the maximum entropy model, which has its own set of rate parameters, to a population of standard mass action models in which the rate constants are randomly chosen. We show that a maximum entropy model is characterized by a high probability of free energy dissipation rate and likewise entropy production rate, relative to other models. We then characterize the variability of the maximum entropy model predictions with respect to uncertainties in parameters (standard free energies of formation) and with respect to ionic strengths typically found in a cell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interfacing electron and neutrino quasielastic scattering cross sections with the spectral function in GENIE

Progress in neutrino-nucleus cross section models is being driven by the need for highly accurate predictions for the neutrino oscillation community. These sophisticated models are being developed within a microscopic description of the nucleus with the goal of encompassing all reaction modes relevant for the accelerator neutrino program. The disconnect between these microscopic models and the event generators that will be used in the next generation of experiments represents a critical obstacle that must be overcome in order to precisely measure the neutrino oscillation parameters. To this end we have developed a hadron tensor interface for lepton-nucleus quasielastic (QE) scattering within the GENIE event generator as a proof of principle, with the broader goal of creating an efficient pipeline for incorporating advanced theoretical models in event generators. As a demonstration of this interface we have implemented the spectral function model into GENIE by connecting theorist provided fortran code through the hadron tensor interface. The spectral function model offers a more complete description of the nuclear ground state, as well as the ability to provide quantifiable theoretical uncertainties. Finally, we validate this implementation and compare its predictions against data and against QE models already available in GENIE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simple nonparametric reconstruction of parton distributions from limited Fourier information

Some calculations of parton distributions from first principles only give access to a limited range of Fourier modes of the function to reconstruct. We present a physically motivated procedure to regularize the inverse integral problem using a Gaussian process as a Bayesian prior. We propose to fix the hyperparameters of the prior in a meaningful physical fashion, offering a simple implementation, great numerical efficiency, and allowing us to understand and keep control easily of the uncertainty of the reconstruction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Inverse problem in the large momentum effective theory framework

One proposal to compute parton distributions from first principles is the large momentum effective theory (LaMET), which requires the Fourier transform of matrix elements computed nonperturbatively. Lattice quantum chromodynamics (QCD) provides calculations of these matrix elements over a finite range of Fourier harmonics that are often noisy or unreliable in the largest computed harmonics. It has been suggested that enforcing an exponential decay of the missing harmonics helps alleviate this issue. Using nonperturbative data, we show that the uncertainty introduced by this inverse problem in a realistic setup remains significant without very restrictive assumptions, and that the importance of the exact asymptotic behavior is minimal for values of 𝑥 where the framework is currently applicable. We show that the crux of the inverse problem lies in harmonics of the order of 𝜆 = 𝑧⁢𝑃 𝑧 ∼ 5–15, where the signal in the lattice data is often barely existent in current studies, and the asymptotic behavior is not firmly established. We stress the need for more sophisticated techniques to account for this inverse problem, whether in the LaMET or related frameworks like the short-distance factorization. We also address a misconception that, with available lattice methods, the LaMET framework allows a “direct” computation of the 𝑥-dependence, whereas the alternative short-distance factorization only gives access to moments or fits of the 𝑥-dependence.

Dutrieux, Hervé [Aix-Marseille Université, Marseil↗

A probabilistic graphical model foundation for enabling predictive digital twins at scale

A unifying mathematical formulation is needed to move from one-off digital twins built through custom implementations to robust digital twin implementations at scale. This work proposes a probabilistic graphical model as a formal mathematical representation of a digital twin and its associated physical asset. We create an abstraction of the asset–twin system as a set of coupled dynamical systems, evolving over time through their respective state spaces and interacting via observed data and control inputs. The formal definition of this coupled system as a probabilistic graphical model enables us to draw upon well-established theory and methods from Bayesian statistics, dynamical systems and control theory. The declarative and general nature of the proposed digital twin model make it rigorous yet flexible, enabling its application at scale in a diverse range of application areas. Here, we demonstrate how the model is instantiated to enable a structural digital twin of an unmanned aerial vehicle (UAV). The digital twin is calibrated using experimental data from a physical UAV asset. Its use in dynamic decision-making is then illustrated in a synthetic example where the UAV undergoes an in-flight damage event and the digital twin is dynamically updated using sensor data. The graphical model foundation ensures that the digital twin calibration and updating process is principled, unified and able to scale to an entire fleet of digital twins.

42 ENGINEERING↗

Proton stopping measurements at low velocity in warm dense carbon

Abstract Ion stopping in warm dense matter is a process of fundamental importance for the understanding of the properties of dense plasmas, the realization and the interpretation of experiments involving ion-beam-heated warm dense matter samples, and for inertial confinement fusion research. The theoretical description of the ion stopping power in warm dense matter is difficult notably due to electron coupling and degeneracy, and measurements are still largely missing. In particular, the low-velocity stopping range, that features the largest modelling uncertainties, remains virtually unexplored. Here, we report proton energy-loss measurements in warm dense plasma at unprecedented low projectile velocities. Our energy-loss data, combined with a precise target characterization based on plasma-emission measurements using two independent spectroscopy diagnostics, demonstrate a significant deviation of the stopping power from classical models in this regime. In particular, we show that our results are in closest agreement with recent first-principles simulations based on time-dependent density functional theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of the Li 6 asymptotic normalization constant onto α -induced reactions of astrophysical interest

Indirect methods have become the predominant approach in experimental nuclear astrophysics for studying several low-energy nuclear reactions occurring in stars, as direct measurements of many of these relevant reactions are rendered infeasible due to their low reaction probability. Such indirect methods, however, require theoretical input that in turn can have significant poorly quantified uncertainties, which can then be propagated to the reaction rates and have a large effect on our quantitative understanding of stellar evolution and nucleosynthesis processes. Here we present two such examples involving α-induced reactions, 13 C (α,n)⁢ 16 O and 12 C (α,γ)⁢ 16 O, for which the low-energy cross sections have been constrained with ( 6 Li,d) transfer data. In this Letter, we discuss how a first-principle calculation of 6 Li leads to a 21% reduction of the 12 C⁡(α,γ) ⁢ 16 O cross sections with respect to a previous estimation. This calculation further resolves the discrepancy between recent measurements of the 13 C (α,n)⁢ 16 O reaction and points to the need for improved theoretical formulations of nuclear reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Applications of emulation and Bayesian methods in heavy-ion physics

Abstract Heavy-ion collisions provide a window into the properties of many-body systems of deconfined quarks and gluons. Understanding the collective properties of quarks and gluons is possible by comparing models of heavy-ion collisions to measurements of the distribution of particles produced at the end of the collisions. These model-to-data comparisons are extremely challenging, however, because of the complexity of the models, the large amount of experimental data, and their uncertainties. Bayesian inference provides a rigorous statistical framework to constrain the properties of nuclear matter by systematically comparing models and measurements. This review covers model emulation and Bayesian methods as applied to model-to-data comparisons in heavy-ion collisions. Replacing the model outputs (observables) with Gaussian process emulators is key to the Bayesian approach currently used in the field, and both current uses of emulators and related recent developments are reviewed. The general principles of Bayesian inference are then discussed along with other Bayesian methods, followed by a systematic comparison of seven recent Bayesian analyses that studied quark-gluon plasma properties, such as the shear and bulk viscosities. The latter comparison is used to illustrate sources of differences in analyses, and what it can teach us for future studies.

Paquet, Jean-François (ORCID:0000000187368171)↗

The Reservoir Temperature Estimator (RTEst): A multicomponent geothermometry tool

The Reservoir Temperature Estimator (RTEst) is a multicomponent geothermometry tool for estimating reservoir geochemical parameters including reservoir temperature, CO 2 fugacity, mass of water lost or gained, and a reaction factor. It estimates these parameters and their associated uncertainties by minimizing an objective function that is the weighted sum of squares of the saturation indices of a user-selected set of minerals believed to be equilibrated with the reservoir fluid. RTEst accomplishes these estimates by combining the geochemical modeling capabilities of The Geochemist’s Workbench® with the optimization/parameter estimation resources of PEST®. An included interface aids the user in selecting plausible mineral phases to comprise the objective function and calculates their weighting factors. The working principles of RTEst are described and its efficacy is illustrated by presenting results of its application to various geothermal fields with known conditions. These examples show RTEst can account for the alteration of ascending reservoir fluid by mineral (calcite) re-equilibration with changes in temperature, reconstruct waters with CO 2 loss, correct for the deficit of water and other volatiles (CO 2,gas , H 2 S gas ) from boiling, and determine the amount of mixing of thermal and non-thermal waters. RTEst can use data with basis species below detection limit, missing, or unreliable either by assuming equilibrium with a controlling mineral (fixed-analyte method) or by treating the analyte concentration as an optimization parameter. The inverse of variance weighting method included in RTEst provides more representative results than either the normalization or unit weighting methods. Finally, the ability of RTEst to calculate reservoir temperatures, gas fugacity, and mixing fractions demonstrates its usefulness as a tool for evaluating geothermal systems.

15 GEOTHERMAL ENERGY↗

Enabling robust offline active learning for machine learning potentials using simple physics-based priors

Abstract Machine learning surrogate models for quantum mechanical simulations have enabled the field to efficiently and accurately study material and molecular systems. Developed models typically rely on a substantial amount of data to make reliable predictions of the potential energy landscape or careful active learning (AL) and uncertainty estimates. When starting with small datasets, convergence of AL approaches is a major outstanding challenge which has limited most demonstrations to online AL. In this work we demonstrate a Δ-machine learning (ML) approach that enables stable convergence in offline AL strategies by avoiding unphysical configurations with initial datasets as little as a single data point. We demonstrate our framework’s capabilities on a structural relaxation, transition state calculation, and molecular dynamics simulation, with the number of first principle calculations being cut down anywhere from 70%–90%. The approach is incorporated and developed alongside AMP torch , an open-source ML potential package, along with interactive Google Colab notebook examples.

Shuaibi, Muhammed↗