SEARCH · Engineering Papers

Results for “probability distribution”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Adapting Reinforcement Learning Treatment Policies Using Limited Data to Personalize Critical Care

We report reinforcement learning (RL) demonstrates promise for developing effective treatment policies in critical care settings. However, existing RL methods often require large and comprehensive patient data sets and do not readily lend themselves to settings in which certain patient subpopulations are severely underrepresented. In this study, we develop a new method, noisy Bayesian policy updates (NBPU), for selecting high-performing reinforcement learning–based treatment policies for underrepresented patient subpopulations using limited observations. Our method uses variational inference to learn a probability distribution over treatment policies based on a reference patient subpopulation for which sufficient data are available. It then exploits limited data from an underrepresented patient subpopulation to update this probability distribution and adapts its recommendations to this subpopulation. We demonstrate our method’s utility on a data set of ICU patients receiving intravenous blood anticoagulant medication. Our results show that NBPU outperforms state-of-the-art methods in terms of both selecting effective treatment policies for patients with nontypical clinical characteristics and predicting the corresponding policies’ performance for these patients.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Deriving the Landauer Principle From the Quantum Shannon Entropy

We derive an expression to determine the equilibrium probability distribution of a quantum state in contact with a noisy thermal environment that formally separates contributions from quantum and classical forms of probabilistic uncertainty. A statistical mechanical interpretation of this probability distribution enables us to derive an expression for the minimum free energy costs for arbitrary (reversible or irreversible) quantum state changes. In conclusion, based on this derivation, we demonstrate that–in contrast to classical systems–the free energy required to erase or reset a qubit depends sensitively on both the fidelity of the target state and on the physical properties of the environment, such as the number of quantum bath states, due primarily to the entropic effects of system-bath entanglement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Atmospheric Drivers of Wind Turbine Blade Leading Edge Erosion: Review and Recommendations for Future Research

Leading edge erosion (LEE) of wind turbine blades causes decreased aerodynamic performance leading to lower power production and revenue and increased operations and maintenance costs. LEE is caused primarily by materials stresses when hydrometeors (rain and hail) impact on rotating blades. The kinetic energy transferred by these impacts is a function of the precipitation intensity, droplet size distributions (DSD), hydrometeor phase and the wind turbine rotational speed which in turn depends on the wind speed at hub-height. Hence, there is a need to better understand the hydrometeor properties and the joint probability distributions of precipitation and wind speeds at prospective and operating wind farms in order to quantify the potential for LEE and the financial efficacy of LEE mitigation measures. However, there are relatively few observational datasets of hydrometeor DSD available for such locations. Here, we analyze six observational datasets from spatially dispersed locations and compare them with existing literature and assumed DSD used in laboratory experiments of material fatigue. We show that the so-called Best DSD being recommended for use in whirling arm experiments does not represent the observational data. Neither does the Marshall Palmer approximation. We also use these data to derive and compare joint probability distributions of drivers of LEE; precipitation intensity (and phase) and wind speed. We further review and summarize observational metrologies for hydrometeor DSD, provide information regarding measurement uncertainty in the parameters of critical importance to kinetic energy transfer and closure of data sets from different instruments. A series of recommendations are made about research needed to evolve towards the required fidelity for a priori estimates of LEE potential.

17 WIND ENERGY↗

DOE OSTI2023

Perturbative readout-error mitigation for near-term quantum computers

Readout errors on near-term quantum computers can introduce significant error to the empirical probability distribution sampled from the output of a quantum circuit. These errors can be mitigated by classical postprocessing given the access of an experimental response matrix that describes the error associated with the measurement of each computational basis state. However, the resources required to characterize a complete response matrix and to compute the corrected probability distribution scale exponentially with the number of qubits, n . In this work, we modify standard matrix inversion techniques using perturbative approximations with significantly reduced complexity and bounded error when the likelihood of high-order bit-flip events is strongly suppressed. Given a characteristic error rate q , we discuss a method to recover the probability of the all-zeros bit string p 0 by sampling only a small subspace of the response matrix before inverting readout error, resulting in a relative speedup of poly [ 2 n / ( n w ) ] , which we motivate using a simplified error model for which the approximation incurs only O ( q w ) error for some integer w . We then provide a generalized technique to efficiently recover full output distributions with O ( q w ) error in the perturbative limit. These approximate techniques for readout-error correction may greatly accelerate near-term quantum computing applications.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Magnetic tunnel junction random number generators applied to dynamically tuned probability trees driven by spin orbit torque

Abstract Perpendicular magnetic tunnel junction (pMTJ)-based true-random number generators (RNGs) can consume orders of magnitude less energy per bit than CMOS pseudo-RNGs. Here, we numerically investigate with a macrospin Landau–Lifshitz-Gilbert equation solver the use of pMTJs driven by spin–orbit torque to directly sample numbers from arbitrary probability distributions with the help of a tunable probability tree. The tree operates by dynamically biasing sequences of pMTJ relaxation events, called ‘coinflips’, via an additional applied spin-transfer-torque current. Specifically, using a single, ideal pMTJ device we successfully draw integer samples on the interval [0, 255] from an exponential distribution based on p -value distribution analysis. In order to investigate device-to-device variations, the thermal stability of the pMTJs are varied based on manufactured device data. It is found that while repeatedly using a varied device inhibits ability to recover the probability distribution, the device variations average out when considering the entire set of devices as a ‘bucket’ to agnostically draw random numbers from. Further, it is noted that the device variations most significantly impact the highest level of the probability tree, with diminishing errors at lower levels. The devices are then used to draw both uniformly and exponentially distributed numbers for the Monte Carlo computation of a problem from particle transport, showing excellent data fit with the analytical solution. Finally, the devices are benchmarked against CMOS and memristor RNGs, showing faster bit generation and significantly lower energy use.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2022

A Power Outage Data Informed Resilience Assessment Framework

Catastrophic impacts to power systems due to disruptive events have increased significantly during the last decade. These events highlight the need to develop approaches to assess the resilience of power systems against extreme events. However, the availability of data that capture power system performance during and after disruptive events is scarce. This paper proposes an assessment framework to evaluate the performance aspects of the grid system during extreme outage events using the Environment for Analysis of Geo-Located Energy Information (EAGLE-I) data. EAGLE-I includes information related to the number of impacted customers, duration, and location of power outages in the United States. Statistical analyses were conducted to extract resilient-based outage data and derive probability distribution functions of their impact and recovery characteristics. A list of extreme events is identified based on few predetermined threshold values. Metrics from other power outage assessments were used to measure the characteristics of each event, including impact rate and duration, recovery rate and duration, and impact level. A probability distribution function is obtained for each metric. The obtained results provide a representation of national grid performance during extreme events, which can be applied as a framework to evaluate various resilience enhancement techniques.

EAGLE-I↗

DOE OSTI2026

Wasserstein normalized autoencoder for anomaly detection

A novel anomaly detection algorithm is presented. The Wasserstein normalized autoencoder (WNAE) is a normalized probabilistic model that minimizes the Wasserstein distance between the learned probability distribution—a Boltzmann distribution where the energy is the reconstruction error of the autoencoder (AE)—and the distribution of the training data. This algorithm has been developed and applied to the identification of semivisible jets—conical sprays of visible standard model (SM) particles and invisible dark matter states—with the CMS experiment at the CERN LHC. Trained on jets of particles from simulated SM processes, the WNAE is shown to learn the probability distribution of the input data in a fully unsupervised fashion, such that it effectively identifies new physics jets as anomalies. The model exhibits stable, convergent training and recovers strong classification performance for a wide range of signals against the selected background process, for which a standard AE fails because of outlier reconstruction. In addition, the model improves upon standard normalized autoencoders while remaining fully agnostic to the signal. The WNAE directly tackles the problem of outlier reconstruction, a common failure mode of autoencoders in anomaly detection tasks.

Hayrapetyan, Aram [Yerevan Phys. Inst.]↗

DOE OSTI2025

Continuing Development of the Nuclear Data Processing Code AMPX [Poster]

The ENDF/B-VIII.1 evaluation library has seen a great growth in the thermal neutron scattering sub-library. The SCALE code system has traditionally approached CE transport by assuming that the CE library on disk represented the fully expanded cumulative probability distributions, conditional on exiting angle and marginal on exiting energy. While this is a complete description of the data, it comes at the potential cost of large amounts of on-disk storage. This approach was strained by several TSL files in ENDF/B-VIII.1, such as graphite, which contained data for a large number of Bragg edges. In the fully expanded probability distributions, this was found to be a disproportionately large fraction of the SCALE CE library.

GNDS↗

DOE OSTI2022

Typological representation of the offshore oceanographic environment along the Alaskan North Slope

Erosion and flooding impacts to Arctic coastal environments are intensifying with nearshore oceanographic conditions acting as a key environmental driver. Robust and comprehensive assessment of the nearshore oceanographic conditions require knowledge of the following boundary conditions: incident wave energy, water level, incident wind energy, ocean temperature and salinity, bathymetry, and shoreline orientation. The number of offshore oceanographic boundary conditions can be large, requiring a significant computational investment to reproduce nearshore conditions. This present study develops location-independent typologies to reduce the number of boundary conditions needed to assess nearshore oceanographic environments in both a Historical (2007–2019) and Future (2020–2040) timespan along the Alaskan North Slope. We used WAVEWATCH III® and Delft3D Flexible Mesh model output from six oceanographic sites located along a constant ~50 m bathymetric line spanning the Chukchi to Beaufort Seas. K-means clustering was applied to the energy-weighted joint-probability distribution of significant wave height (H s ) and peak period (T p ). Distributions of wave and wind direction, wind speed, and water level associated with location-independent centroids were assigned single values to describe a reduced order, typological rendition of offshore oceanographic conditions. Reanalysis data (e.g., ASRv2, ERA5, and GOFS) grounded the historical simulations while projected conditions were obtained from downscaled GFDL-CM3 forced under RCP8.5 conditions. Location-dependence for each site is established through the occurrence joint-probability distribution in the form of unique scaling factors representing the fraction of time that the typology would occupy over a representative year. As anticipated, these typologies show increasingly energetic ocean conditions in the future. They also enable computationally efficient simulation of the nearshore oceanographic environment along the North Slope of Alaska for better characterization of coastal processes (e.g., erosion, flooding, or sediment transport).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Efficient First-Order Algorithms for Large-Scale, Non-Smooth Maximum Entropy Models with Application to Wildfire Science

Maximum entropy (MaxEnt) models are a class of statistical models that use the maximum entropy principle to estimate probability distributions from data. Due to the size of modern data sets, MaxEnt models need efficient optimization algorithms to scale well for big data applications. State-of-the-art algorithms for MaxEnt models, however, were not originally designed to handle big data sets; these algorithms either rely on technical devices that may yield unreliable numerical results, scale poorly, or require smoothness assumptions that many practical MaxEnt models lack. In this paper, we present novel optimization algorithms that overcome the shortcomings of state-of-the-art algorithms for training large-scale, non-smooth MaxEnt models. Our proposed first-order algorithms leverage the Kullback–Leibler divergence to train large-scale and non-smooth MaxEnt models efficiently. For MaxEnt models with discrete probability distribution of n elements built from samples, each containing m features, the stepsize parameter estimation and iterations in our algorithms scale on the order of O(mn) operations and can be trivially parallelized. Moreover, the strong ℓ1 convexity of the Kullback–Leibler divergence allows for larger stepsize parameters, thereby speeding up the convergence rate of our algorithms. To illustrate the efficiency of our novel algorithms, we consider the problem of estimating probabilities of fire occurrences as a function of ecological features in the Western US MTBS-Interagency wildfire data set. Our numerical results show that our algorithms outperform the state of the art by one order of magnitude and yield results that agree with physical models of wildfire occurrence and previous statistical analyses of wildfire drivers.

Physics↗

DOE OSTI2023

Biophotonic rogue waves in red blood cell suspensions

Rogue waves are ubiquitous in nature, appearing in a variety of physical systems ranging from acoustics, microwave cavities, optical fibers, and resonators to plasmas, superfluids, and Bose–Einstein condensates. Unlike nonlinear solitary waves, rogue waves are extreme events that can occur even without nonlinearity by, for example, spontaneous synchronization of waves with different spatial frequencies in a linear system. Here, we report the observation of rogue-wave-like events in human red blood cell (RBC) suspensions under weak light illumination, characterized by an abnormal L-shaped probability distribution. Such biophotonic extreme events arise mostly due to the constructive interference of Mie-scattered waves from the suspended RBCs, whose biconcave shape and mutable orientation give rise to a time-dependent random phase modulation to an incident laser beam. We numerically simulate the beam propagation through the colloidal suspensions with added disorder in both spatial and temporal domains to mimic random scattering due to Brownian motion. In addition, at high power levels, nonlinear beam self-focusing is also observed, leading to a dual-exponential probability distribution associated with the formation of multiple soliton-like spots. Such rogue wave events should also exist in environments with cells of other species such as swimming bacteria, and understanding of their underlying physics may lead to unexpected biophotonic applications.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Dynamics of heavy quarks in strongly coupled $\mathcal{N}$ = 4 SYM plasma

We calculate the probability distribution P(k) for a heavy quark with velocity v propagating through strongly coupled N = 4 SYM plasma in the ’t Hooft limit (N c → ∞, λ = g 2 N c → ∞) at a temperature T to acquire a momentum k due to interactions with the plasma. This distribution encodes the well-known drag coefficient η D and the transverse and longitudinal momentum diffusion coefficients κ T and κ L . The jet quenching parameter $\hat{q}$ can be extracted from P(k) for v = 1. Going beyond these known Gaussian characteristics of P(k), our calculation determines all of the higher order and mixed moments to leading order in 1/$\sqrt{λ}$ for the first time. These non-Gaussian features of P(k) include qualitatively novel correlations between longitudinal energy loss and transverse momentum broadening at nonzero v. We show that all higher moments scale characteristically with an effective temperature of the boosted plasma in the heavy quark rest frame, and we demonstrate that these non-Gaussian characteristics can be sizable in magnitude and even dominant in physically relevant situations. We use these results to derive a Kolmogorov equation for the evolution of the probability distribution for the total momentum of a heavy quark that propagates through strongly coupled plasma. This evolution equation accounts for all higher order correlations between transverse momentum broadening and longitudinal energy loss, which we have calculated from first principles. It reduces to a Fokker-Planck equation when truncated to only include the effects of η D , κ T and κ L . Remarkably, while heavy quarks do not reach kinetic equilibrium with the plasma if evolved with this Fokker-Planck equation, by showing that the Boltzmann distribution is a static solution of the all-order Kolmogorov equation that we have derived we demonstrate that heavy quarks do reach kinetic equilibrium if evolved with this equation. Our results thus provide a dynamically complete framework for understanding the thermalization of a heavy quark that may be initially far from equilibrium in the strongly coupled N = 4 SYM plasma — as well as new insight into heavy quark transport and equilibration in quark-gluon plasma.

Holography and Hydrodynamics↗

DOE OSTI2023

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Using Signal Clustering Similarity for Detecting CAN Masquerade Attacks

The computer code assumes that time series representing the physical signals of the vehicle have been extracted from the CAN bus. The main input of the computer code is the multivariate time series representation of the signals in the CAN bus. The computer code cluster these time series using agglomerative hierarchical clustering from benign and attack datasets. Based on this, it generates probability distributions from the similarity of the obtained clusters based in each scenario---benign and attack---using the CluSim method (https://github.com/Hoosier-Clusters/clusim). Finally, it compares how a new data collection compares with the previous distribution to provide and probability score for an intrusion.

Moriano, Pablo↗

DOE OSTI2023

Approximate Boltzmann distributions in quantum approximate optimization

Approaches to compute or estimate the output probability distributions from the quantum approximate optimization algorithm (QAOA) are needed to assess the likelihood it will obtain a quantum computational advantage. We analyze output from QAOA circuits solving 7200 random MaxCut instances, with $n$ = 14–23 qubits and depth parameter $p$ ≤ 12 and find that the average basis state probabilities follow approximate Boltzmann distributions: The average probabilities scale exponentially with their energy (cut value), with a peak at the optimal solution. Furthermore, we describe the rate of exponential scaling or effective temperature in terms of a series with a leading-order term $T$ ~ $C$ min /$n$ $\sqrt{p}$, with $C$ min the optimal solution energy. Using this scaling, we generate approximate output distributions with up to 38 qubits and find these give accurate accounts of important performance metrics in cases we can simulate exactly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Simulated effect of edge plasma density parameters on lower hybrid wave scattering in EAST

The incorporation of lower hybrid (LH) wave spectrum broadening in the poloidal wavenumber (⁠k θ ) space at the last close field surface (LCFS) is reported to lead to better agreement of the modeled LH wave current profile with the experimental results [Baek et al., Nucl. Fusion 61, 106034 (2021)]. To further understand its underlying mechanism and find the possible influence factors, effects of wave scattering caused by drift-wave type density fluctuation on the probability distribution of the LH wave polar refractive index (⁠N θ ) at the LCFS are studied under density parameters in the scrape-off-layer. According to a scattering model [P. T. Bonoli and E. Ott, Phys Fluids 25(2), 359–375 (1982)], scattering probability and scattering angle distribution are two main factors that determine the degree of spectral broadening. Studies presented here show that the total scattering probability increases first and then decreases as the wave propagates toward a smaller normalized radius of poloidal magnetic flux (⁠ρ⁠). The degree of spectral broadening is found to depend on the density and density fluctuation together by changing the intensity and a proportion of the geometrical optics approximation term and the E×B drift term in the scattering model. Additionally, the fluctuation correlation length can significantly modify the probability distribution of N θ at the LCFS, which is found to significantly change the LH wave current profile.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

A Eulerian three-fluid flow framework for the study of fuel dispersal behavior under loss-of-coolant accident conditions

As the nuclear industry explores new fuel designs to accommodate increased burnup, studying fuel behavior during loss-of-coolant accidents is essential for ensuring the safe operation of light-water reactors. Here, this study presents a computational model to simulate the intricate three-phase flow dynamics of fuel dispersal following cladding breach, where solid fuel particles and carrier fission gases interact with the surrounding liquid or gas medium. The model utilizes a Eulerian three-fluid framework to depict the phase interaction within interpenetrating continua, treating fragmented fuel as granular material. Closure relationships for bulk-flow properties of the solid phase and detailed interfacial exchange terms in high-temperature, high-speed gas-solid-liquid flow scenarios are proposed. The model was validated using adiabatic experimental results in the literature, focusing on characterizing fuel dispersal behavior during loss-of-coolant accidents conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched probability density functions reported in experimental studies. Key highlights of this study include the theoretical description of a coupled three-phase flow with phase change, interpretation of the high-pressure boundary conditions and insights into the transient behavior of gas-solid-liquid phase dynamics during a simulated high-speed dispersal event.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Asymptotic inconsistency of the cumulative algorithm for laser-induced damage probability analysis

The “cumulative algorithm” is a data analysis method that has been proposed to provide an objective, nonparametric determination of laser-induced damage probability as a function of fluence from experimental data that contain both damaged sites and undamaged sites (i.e., 1-on-1 or S-on-1 testing protocols). In this work, the limitations of this approach are explored by considering the asymptotic limit of a large number of test sites. It is shown that the cumulative algorithm does not converge to the true probability distribution and significantly underestimates the damage probability near the damage onset. Here, based on the results of this work, the cumulative algorithm is not recommended for accurate estimation of damage probability.

Computational methods↗