Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conditional 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.

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

At least 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Grid Topology Discovery Algorithm Evaluation of Suitability for Utility Deployment (CRADA 606 Final Report)

This work presents the results of a field-informed demonstration aimed at evaluating the practical suitability of a topology discovery algorithm for utility environments. We demonstrated an algorithm that uses a graph-theory-informed state estimation approach for model selection. In collaboration with Survalent and Peninsula Light Co., the algorithm was applied to real feeder models and field measurements from supervisory control and data acquisition (SCADA) and advanced metering infrastructure (AMI) systems to identify the operational topology of a power distribution system. The demonstration assessed the algorithm’s performance under realistic data conditions, including sparse and noisy measurements, and examined its ability to identify the most likely network configurations. The results confirmed that the approach can effectively narrow down feasible topologies, providing operators with improved situational awareness of network status. Key lessons learned emphasize the need for systematic data validation and strategic sensor placement to enhance observability. These insights inform future deployment strategies and guide refinements for broader adoption in utility operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Fundamental chemical physics revealed by scattering reactive open-shell atoms from surfaces

The dynamics of reactive atoms at surfaces are centrally important to areas such as heterogeneous catalysis, corrosion, materials degradation in extreme environments, and plasma etching. Remarkably detailed understanding of dynamical processes at surfaces has been extracted from scattering molecules and inert atoms under well-defined conditions. However, traditional techniques for generating beams of reactive atoms often result in impure mixtures, broad energy distributions, and poorly defined contributions of metastable electronically excited atoms. In this perspective article, we review the state-of-the-art in reactive atom surface scattering with a focus on experiments performed under controlled conditions on well-defined surfaces. We highlight a new technique for controlled state-to-state scattering of polyelectronic atoms from surfaces, based on vacuum ultraviolet photolysis and state-selective ion imaging. The new capabilities provide an avenue for research into the underexplored area of excited state and spin selective chemical dynamics at surfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Elucidating the phase transformations and grain growth behavior of O3-type sodium-ion layered oxide cathode materials during high temperature synthesis

Understanding the formation mechanism of layered oxide cathodes via solid-state synthesis is imperative to achieving controllability over their materials properties and electrochemical behaviors. In this work, we investigate the phase and microstructure evolution during the synthesis of NaNi 1/3 Fe 1/3 Mn 1/3 O 2 , a model sodium-ion layered oxide cathode, using a combination of imaging, diffraction, and spectroscopic techniques. We unravel the synthetic mechanistic pathways involved in the high-temperature calcination reaction, as well as elaborate the synthesis-microstructure-performance relationship of this material. The formation of the final layered oxide phase involves a gradual transformation through a sodiated oxyhydroxide intermediate. During the reaction, the precursor dehydration reaction dominates at 250–550 °C, while the major sodiation reaction occurs at 550–850 °C. Alongside multiple stages of phase transformations, the final grain structure formation occurs through the continuous growth of the (003) and (104) facets. During the reaction, Mn acts as the charge-compensating element and exhibits depth-dependent characteristics. When the sodiation reaction dominates over dehydration, the reaction intermediates undergo gradual electronic structure changes with increasing temperature, as indicated by the spectral features of TM3d-O2p hybrid states. Calcination duration is also a critical parameter governing the microstructure, surface reactivity, phase fraction distribution and electrochemical performance of the material. The optimal calcination duration was determined to be 18 hours at 850 °C under the conditions evaluated here. Calcination beyond this duration was found to be detrimental to electrochemical performance due to Na and O loss and heterogeneous sodium distribution throughout the particles. Our work sheds light on the complex crystallographic-chemical-microstructural evolution of sodium ion layered oxide cathodes and provides insight into precisely tuning material properties which are intimately linked to battery performances.

25 ENERGY STORAGE↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Learning earthquake ground motions via conditional generative modeling

Predicting high-fidelity ground motions for future earthquakes is crucial for seismic hazard assessment and infrastructure resilience. Conventional empirical simulations suffer from sparse sensor distribution and geographically localized earthquake locations, while physics-based methods are computationally intensive and require accurate representations of Earth structures and earthquake sources. We propose an artificial intelligence (AI) spectrogram generator, Conditional Generative Modeling for Ground Motion (CGM-GM). CGM-GM leverages earthquake magnitudes and geographic coordinates of earthquakes and sensors as inputs, when postprocessed with phase information, capturing spatially continuous Fourier amplitude spectra (FAS) as well as properties such as P and S arrivals, and waveform durations, without explicit physics constraints. This is achieved through a probabilistic autoencoder that extracts latent distributions in the time-frequency domain and variational sequential models for prior and posterior distributions. We evaluate the performance of CGM-GM using small-magnitude earthquake records from the San Francisco Bay Area, a region with high seismic risks. Here, we report that CGM-GM demonstrates potential for complementing physics-based simulations and non-ergodic empirical ground motion models, as well as shows promise in seismology and beyond.

geophysics↗
DOE OSTI2025

Return on Investment and Sustainability of HVDC Links: Role of Diagnostics, Condition Monitoring, and Material Innovations

HVDC cable systems are becoming an upscaled technical option, compared to AC, because of various factors, including easier interconnections, lower losses, and longer transmission distances. In addition, renewables providing direct DC energy, electrified transportation, and aerospace where DC can be favored because of higher carried specific power all point in the direction of broad future usage of HV and MV DC links. However, contrary to AC, there is little return from on-field installation as regards long-term cable reliability and aging processes. This gap must be covered by intensive research, and contributing to this research is the purpose of this paper. The focus is on key points for HVDC (and MVDC) cable reliability and sustainability, from design modeling able to account for voltage transients and extrinsic aging (such as that caused by partial discharges) to the impact of aging on insulation conductivity (which rules the electric field distribution, thus aging rate). Also, recyclable and nanostructured materials, as well as health conditions, are considered. It is shown how cable design can account for accelerated aging due to voltage transients, as well as for aging-time dependence of conductivity, and how design can be free of extrinsic aging caused by PDs. Algorithms for health condition evaluations, which have additional value in a relatively new technology such as HVDC polymeric cables, are applied to insulation system aging under partial discharges, showing how they can provide an indication of insulation degradation globally or locally (weak spots) and of possible maintenance times. All of this can effectively contribute to reducing the risk of major cable breakdown and damage under operation, which would significantly affect the return on investment (ROI).

Montanari, Gian Carlo (ORCID:0000000320258693)↗
DOE OSTI2024

Updimensioning strategy derived from synthetic equiaxed grain structures for approximating 3D grain size distributions from 2D visualizations with 1D parameters

We generated synthetic equiaxed grain structures using computer graphics software to explore the relationship between various grain size determination methods and true three-dimensional (3D) grain diameters. Mirroring grain measurement techniques, the synthetic 3D grain structures are imaged as 2D micrographs which are measured to yield 1D grain size parameters. Synthetic grain structures provide data at a mass scale and permit exploration of both polished and fractured surface micrographs, revealing one-to-one correspondence between exposed 2D grain cross-sections and individual 3D grains. Analysis of this correspondence yielded a procedure to approximate 3D equiaxed grain size and volume distributions based on the mode of the 2D fractograph grain size distribution. The 3D approximation procedure is shown to be less susceptible to different imaging conditions that affect small, undiscernible grains compared to the standard planimetric and linear intercept methods, which by design also tend to underestimate the 3D grain diameter. The procedure requires larger sample sizes to lower variance and a deeper analysis which could become more practical with machine learning (ML) models for grain boundary segmentation, which synthetic grain structures can help train. This work lays the foundation for analyzing other grain distributions such as columnar and composite grains in similar depth.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Supersaturation, Nucleation, and Phase Separation of Mesoscopic Systems

Supersaturation, nucleation, and phase separation are ubiquitous phenomena of great interest in both science and industry. However, a unified, quantitative understanding of these phenomena has yet to be achieved for mesoscopic systems. Here, we present a set of general equations that determine the monomer saturation degree, the size distribution, and the free energy of mesoscopic systems, as well as their phase-transition conditions. These equations reveal that, under supersaturation, the largest cluster size (LCS) is an important state variable; the supersaturation degree decreases with the LCS, approaching unity in the macroscopic limit. We identify the critical supersaturation, at which the nuclei undergo the phase transition to form large crystals. Below this critical supersaturation, the nucleus size distribution is either a unimodal function or a monotonically decreasing function of size, depending on the system and temperature. We also predict the most probable nucleus size and the direction of spontaneous changes of the LCS. Our theory provides a unified, quantitative explanation of the nucleus-size-distribution across six different systems, including nanoparticles and biological condensates. This work serves as a general theoretical framework useful for understanding and designing nucleation and phase transitions of mesoscopic systems.

Kang, Jingyu↗
DOE OSTI2024

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion↗
DOE OSTI2025

Quasiclassical sampling and Wigner sampling of initial vibrational coordinates and momenta for polyatomic molecules in Monte Carlo molecular dynamics simulations

In a quasiclassical trajectory simulation, the vibrational modes are initialised with quantised vibrational energies, but vibrational phases are sampled by Monte Carlo. This requires an algorithm to assign coordinates and momenta to the various atoms. In this work, we present two methods for implementing this for nonrotating polyatomic molecules, namely, fixed-energy vibrational-state-selected initial conditions and thermal initial conditions. We also present a method for initiating classical trajectories with a ground-state Wigner distribution. These vibrational treatments are sufficient to initialise trajectories for unimolecular processes, and we also show how they can be applied to simulate bimolecular collision processes. The treatments of unimolecular and bimolecular collision processes are available in two Python codes called wigner_state_selected.py and bimolecular_collision.py, respectively, which will generate initial condition files that are recognisable by the SHARC and SHARC-MN computer programs for dynamics calculations. Both codes are available as standalone programs, as well as being included in SHARC-MN, and they will be included in future versions of SHARC. Here, the methods implemented in these codes are mostly also available in the ANT computer program, and those that are not available in ANT will be incorporated in future versions of ANT.

Wigner distribution↗
DOE OSTI2024

Detecting Anomalies for Fire Prevention in Distribution Systems: Challenges and Analytical Techniques

Electric utilities in California have historically been linked to up to 10% of wildfires. To mitigate this risk, Southern California Edison has invested significantly in wildfire prevention strategies, including undergrounding cables and enhancing equipment inspections. This article explores a novel approach to fire prevention by detecting anomalies in the distribution system that may indicate potential fire hazards. The focus is on identifying arcing conditions through high-resolution point-on-wave (POW) measurements. Arcing, a precursor to fires, is challenging to detect due to its subtle transients and complex system topology. The article discusses the use of advanced signal processing and machine learning techniques, such as spectral correlation function and discrete wavelet transform, to extract features from POW data and accurately identify arcing events. The study demonstrates a high accuracy rate in detecting arcing, paving the way for improved fire prevention measures in electric distribution systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Estimating the Contributions to Human Error Probability from the Convolution of the Distribution of Time Available and Time Required

As part of their duties, Human Reliability Analysis must often evaluate if crews in nuclear power plants (NPPs) can complete tasks associated with a human-failure event within time limits. For example, the time required in NPP scenarios is determined by systematic and structured walkthroughs, feasibility studies, recorded times from training exercises, and interviews with experienced operators and experts. Typically, a point estimate is derived for the estimate (mean, maximum, or 95th percentile of time required). Using point-estimate values can mask the risk associated with variability among crews, plant conditions and set-up, environmental conditions, and other impact factors under which these actions are executed. While point estimates for time required and time available have served the industry well, without considering the uncertainty they could lead to biased understanding about the risk. The Integrated Human Event Analysis System - General Methodology (IDHEAS-G) model (developed by the US Nuclear Regulatory Commission, NRC) for human error probability calculates human error probability by summing two probabilities: insufficient time and cognitive error. As such, the model takes a more holistic approach by considering the full distributions for time required and time available to calculate the human error probability because the time available to complete the task is insufficient. In this study, we expand on the work of the NRC and discuss methods for estimating these time considerations. For example, for the time required, the impact of Performance Influencing Factors (PIFs) on the distribution was divided into impacts that are aleatory in nature, such as crew-to-crew variability, and those that are epistemic (i.e., the PIFs). Starting with the factors that introduce aleatory uncertainty, a first-order distribution was developed from a large set of time required (i.e., NPP task completion times) data for the range of operator actions that occur in the NPP control room under simulated accident conditions. The first-order distribution can then be adjusted to account for epistemic uncertainty using research associated with the impact of applicable PIFs on the time required. We also develop guidance for analysts to address the probability distributions for the time available. The guidance we developed on how to estimate time required and time available distributions is based on the identification of pertinent research and data, data analyses, and expert knowledge elicitation.

human error probability, human performance, time e↗
DOE OSTI2024

Leveraging Extremal Dependence to Better Characterize the 2021 Pacific Northwest Heatwave

Abstract In late June, 2021, a devastating heatwave affected the US Pacific Northwest and western Canada, breaking numerous all-time temperature records by large margins and directly causing hundreds of fatalities. The observed 2021 daily maximum temperature across much of the U.S. Pacific Northwest exceeded upper bound estimates obtained from single-station temperature records even after accounting for anthropogenic climate change, meaning that the event could not have been predicted under standard univariate extreme value analysis assumptions. In this work, we utilize a flexible spatial extremes model that considers all stations across the Pacific Northwest domain and accounts for the fact that many stations simultaneously experience extreme temperatures. Our analysis incorporates the effects of anthropogenic forcing and natural climate variability in order to better characterize time-varying changes in the distribution of daily temperature extremes. We show that greenhouse gas forcing, drought conditions and large-scale atmospheric modes of variability all have significant impact on summertime maximum temperatures in this region. Our model represents a significant improvement over corresponding single-station analysis, and our posterior medians of the upper bounds are able to anticipate more than 96% of the observed 2021 high station temperatures after properly accounting for extremal dependence. Supplementary materials accompanying this paper appear online.

Zhang, Likun (ORCID:0000000154584556)↗
DOE OSTI2024

A potential model study of the nucleon's charge and mass radius

We study the charge and mass distributions within a nucleon and compute the associated squared radii based on a potential model approach. Different constituent quark configurations such as Δ, Y, and quark-diquark are considered and compared, with model parameters calibrated by experimental measurements of the proton and neutron charge radius. The results suggest that while the charge radius is dictated by quark dynamics, the mass radius is strongly influenced by nonperturbative QCD contributions to a nucleon's mass that are not sensitive to the constituent quarks. As a result, the mass radius could become substantially different from the charge radius. The obtained nucleon mass distributions of different configurations are further used for simulations of the initial conditions in heavy ion collisions. Furthermore, the computed eccentricities ε 2 and ε 3 are found to demonstrate a considerable sensitivity to the input nucleon profiles, especially to the mass radius in the peripheral region as well as for systems with fewer participants.

Gallimore, Daniel [Indiana University, Bloomington↗
DOE OSTI2025

Unavoidable multilevel biaxial symmetry breaking in chiral hybrid liquid crystals

Chiral nematic or cholesteric liquid crystals (LCs) are fluid mesophases with long-ranged orientational order featuring a quasilayered periodicity imparted by a helical director configuration but lacking long-range positional order. Doping molecular cholesteric LCs with strongly anisotropic uniaxial colloidal particles adds another level of complexity because of the interplay between weak surface-anchoring boundary conditions and bulk-based elastic distortions near the particle-LC interface. Using cylindrical colloidal disks and rods with different geometric shapes and surface conditions, we demonstrate that these colloidal inclusions generically exhibit biaxial orientational probability distributions which may impart anomalously strong local biaxiality onto the hybrid cholesteric LC structure. Unlike nonchiral hybrid molecular-colloidal LCs, where biaxial order emerges only at critical colloid volume fractions exceeding some uniaxial-biaxial transition value, the orientational probability of the colloidal inclusions immersed in chiral nematic hosts is unambiguously biaxial even at infinite dilution. We demonstrate that the colloids induce local biaxial perturbations within the molecular orientational order of the LC host medium which strongly enhances the weak but native biaxial order of chiral nematic LC induced by the chiral symmetry breaking of the director field. With the help of analytical modeling and computer simulations based on the Landau–de Gennes free energy of the host LC around the colloids, we rationalize the observed multilevel biaxial order and conclude that it is not only unavoidable but also strongly enhanced compared to both achiral hybrid LCs and purely molecular cholesteric LCs.

biaxial liquid crystals↗
DOE OSTI2024

Precise measurements of W - and Z -boson transverse momentum spectra with the ATLAS detector using pp collisions at $\sqrt{s} = 5.02$ TeV and 13 TeV

This paper describes measurements of the transverse momentum spectra of W and Z bosons produced in proton–proton collisions at centre-of-mass energies of $\sqrt{s}$ = 5.02 TeV and $\sqrt{s}$ = 13 TeV with the ATLAS experiment at the Large Hadron Collider. Measurements are performed in the electron and muon channels, W → $\ell$$v$ and Z → $\ell$$\ell$ ($\ell$ = e or μ), and for W events further separated by charge. The data were collected in 2017 and 2018, in dedicated runs with reduced instantaneous luminosity, and correspond to 255 and 338 pb -1 at $\sqrt{s}$ = 5.02 TeV and 13 TeV, respectively. These conditions optimise the reconstruction of the W-boson transverse momentum. The distributions observed in the electron and muon channels are unfolded, combined, and compared to QCD calculations based on parton shower Monte Carlo event generators and analytical resummation. The description of the transverse momentum distributions by Monte Carlo event generators is imperfect and shows significant differences largely common to W - , W + and Z production. The agreement is better at $\sqrt{s}$ = 5.02 TeV, especially for predictions that were tuned to Z production data at $\sqrt{s}$ = 7 TeV. Higher-order, resummed predictions based on DYTurbo generally match the data best across the spectra. Distribution ratios are also presented and test the understanding of differences between the production processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Software For Advanced Large-scale Analysis Of Magnetic Confinement For Numerical Design, Engineering & Research (salamander)

As magnetic confinement fusion energy gains traction internationally to enable abundant energy production, designing components for fusion systems is a pressing challenge. During the planned lifetime of a fusion device, components evolve in extreme environments and must withstand large, repeated thermal loads and bombardment by 14 MeV neutrons, plasma ions, and neutral particles (deuterium, tritium, and helium), corrosive conditions, etc. All these physical processes take place simultaneously, interact in intricate ways, and impose important constraints that can affect performance. Experimental data is rare and costly to obtain, making design particularly challenging. Predictive computational frameworks must be an integral part of an accelerated and cost-effective design process by modeling fusion system performance in simulated environments. To better understand component degradation and operational impacts on their performance, the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER) is designed as an open-source, fully integrated, multiphysics, multiscale, NQA-1 compliant framework facilitating 3D, high-fidelity fusion system modeling. To that end, SALAMANDER is a MOOSE-based framework, and therefore leverages MOOSE upstream libraries such as PETSc and libMesh to deliver sophisticated finite element, finite volume, and nonlinear solver technology for fusion energy simulations. SALAMANDER couples MOOSE physics module capabilities—such as thermal hydraulics, heat conduction, Navier-Stokes, and thermomechanics—with tritium transport via TMAP8, neutronics via Cardinal, and nascent particle-in-cell capabilities. Direct simulation Monte Carlo methods will be used to address neutral transport near the walls. By coupling all these physics in an integrated application, SALAMANDER will enable high-fidelity modeling of irradiation levels and plasma exposure conditions of plasma facing components and their impact on heat and tritium distributions, as well as the resulting mechanical constraints experienced by the plasma facing components and performance of blanket systems. Furthermore, SALAMANDER will be particularly suited for engineering studies thanks to the stochastic tool module readily available in MOOSE, allowing for extended uncertainty quantification and risk analysis studies. It is also able to use computer-aided design (CAD) meshes to model complex geometries, which is indispensable for fusion systems. SALAMANDER therefore supports design, safety, engineering, and research projects for magnetic confinement fusion systems

Simon, Pierre-Clement [Idaho National Laboratory (↗
DOE OSTI2025

Thermal Analysis of the Bus Bars for the 800MeV Injection Era (Fermilab Summer School 2023)

My internship at IARC has been focused on the thermal analysis of the bus bars of the Booster. The bus bars are the conductors that bring the bias current from the bias supply to the ferrite tuners in each RF cavity of the synchrotron. Since they are not superconductors, during normal operations they will heat up due to the Joule effect. The goal of my work is to build a model of the bus bars assembly that describes effectively the experimental data and that can be used to predict the temperature increase and distribution among the components of the assembly, in present and future operational conditions.

43 PARTICLE ACCELERATORS↗