Engineering Papers⌕ Search

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.

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

At least 685 records · Page 38

Solving Stochastic Inverse Problems for Property–Structure Linkages Using Data-Consistent Inversion and Machine Learning

Determining process–structure–property linkages is one of the key objectives in material science, and uncertainty quantification plays a critical role in understanding both process–structure and structure–property linkages. In this work, we seek to learn a distribution of microstructure parameters that are consistent in the sense that the forward propagation of this distribution through a crystal plasticity finite element model matches a target distribution on materials properties. This stochastic inversion formulation infers a distribution of acceptable/consistent microstructures, as opposed to a deterministic solution, which expands the range of feasible designs in a probabilistic manner. Furthermore, to solve this stochastic inverse problem, we employ a recently developed uncertainty quantification framework based on push-forward probability measures, which combines techniques from measure theory and Bayes’ rule to define a unique and numerically stable solution. This approach requires making an initial prediction using an initial guess for the distribution on model inputs and solving a stochastic forward problem. To reduce the computational burden in solving both stochastic forward and stochastic inverse problems, we combine this approach with a machine learning Bayesian regression model based on Gaussian processes and demonstrate the proposed methodology on two representative case studies in structure–property linkages.

36 MATERIALS SCIENCE↗

Inequality in the availability of residential air conditioning across 115 US metropolitan areas

Continued climate change is increasing the frequency, severity, and duration of populations’ high temperature exposures. Indoor cooling is a key adaptation, especially in urban areas, where heat extremes are intensified—the urban heat island effect (UHI)—making residential air conditioning (AC) availability critical to protecting human health. In the United States, the differences in residential AC prevalence from one metropolitan area to another is well understood, but its intra-urban variation is poorly characterized, obscuring neighborhood-scale variability in populations’ heat vulnerability and adaptive capacity. We address this gap by constructing empirically derived probabilities of residential AC for 45,995 census tracts across 115 metropolitan areas. Within cities, AC is unequally distributed, with census tracts in the urban “core” exhibiting systematically lower prevalence than their suburban counterparts. Moreover, this disparity correlates strongly with multiple indicators of social vulnerability and summer daytime surface UHI intensity, highlighting the challenges that vulnerable urban populations face in adapting to climate-change driven heat stress amplification.

54 ENVIRONMENTAL SCIENCES↗

Decoding ionic conductivity and reordering in cation-disordered pyrochlores

The ordered structure A 2 B 2 O 6 O’ in pyrochlores engenders twin rows of inequivalent anion sublattices each centred on alternating cations. While it is known that cation antisite disorder augments the ionic conductivity by several orders of magnitude, the local cation environment around the anions and the dynamic anion reordering during the cation disordering are not well-elucidated. Using atomistic simulations on Gd 2 Zr 2 O 7 , we first show that the anions engage in concerted hops to the neighbouring tetrahedral sites mostly along with the $\langle$1 0 0$\rangle$ direction while completely avoiding the octahedral sites. While the initially vacant 8 a sites start accommodating oxygen ions with increasing cation disorder, they show noticeable reluctance even at significant levels of disorder. We have also tracked both the distribution of available oxygen sites following random cation disorder, which is dependent only on cation disordering, and the probability of occupation of these sites. Interestingly, the available oxygen sites show a non-monotonic dependence on the number of B ions in the nearest neighbouring shell while the occupation probability of all the available oxygen sites increases monotonically. A tetrahedral oxygen site thus has a better probability of being occupied when it has a greater number of second neighbour B ions. This article is part of the Theo Murphy meeting issue 'Understanding fast-ion conduction in solid electrolytes'.

Science & Technology - Other Topics↗

Photonuclear tomography in ultraperipheral heavy-ion collisions

We present a theoretical investigation of photonuclear tomography as a novel technique for probing the internal structure of nuclei. In this approach, ultraperipheral heavy-ion collisions (UPCs) serve as a source of intense fluxes of virtual photons, which induce coherent production of vector mesons. By analyzing the probabilities and cross sections of these photon-induced processes, we propose a methodology for reconstructing the spatial distribution of nucleons within the nucleus. Our framework provides a systematic way to access information on the nuclear geometry probed in UPCs, offering new opportunities for studies of nuclear structure using particle production as a probe. Numerical calculations for selected examples illustrate the feasibility and potential of this method.

particle & resonance production↗

Cybersecurity Value-at-Risk Framework: Preprint

As more variable renewable energy sources are added to the grid, the critical role of hydropower as a reliable baseline and firming resource is rapidly growing. However, the U.S hydropower fleet is not fully prepared to face modern issues such as cybersecurity threats. Hydropower accounts for 37% of U.S. utility-scale renewable electricity but is challenged by diverse infrastructure and legacy devices that predate modern security practices. While new cybersecurity solutions cannot simply be added to current hydropower generation and operation technologies, custom cybersecurity assessments can reveal system-specific threats and risk probabilities and identify mitigating enhancements.

cybersecurity valuation↗

Cumulative short-lived photofission product yields for nuclear forensic application

Experiments have long been carried out with linear accelerator (linac) produced bremsstrahlung X-rays to study high-energy photon-induced fission, known as “photofission”. One benefit of photofission, when compared with neutron-induced fission, is the ability to investigate fundamental physics of fission at excitation energies lower than what is possible with neutrons [1]. Absorption of an incident neutron results in a compound nucleus with an excitation energy at least equal to the neutron binding energy; thus, the lowest energy region of the fission barrier cannot be studied. Fission product yield distributions are dependent on the incident particle energy and the target nucleus mass. They are generally represented by a double-humped curve, with low and high mass “peaks” and an intermediate mass “valley”. As incident particle energy increases, the excitation energy of the nucleus increases, thereby increasing the overall yield of isotopes in the valley with respect to the two maxima or peaks of the fission product curve [3]. Photofission is most probable for incident photons around 14 MeV due to the giant dipole resonance. Induced photofission around this energy results in a fission product yield distribution that resembles 14 MeV DT neutron fission more than thermal or fast (500 keV) neutron fission. In post-detonation nuclear forensics, the peak-to-valley ratio of the fission product yield curve could be used for characterizing device type by providing an indication of the neutron energy spectrum that had induced fission (e.g., 500 keV vs. 14 MeV). The goal of this work is to improve nuclear data for isotopes of interest in nuclear forensics and special nuclear material (SNM) detection, specifically the fission product yield data via photofission experiments for the isotopes of 238U and 232Th. Photofission has been proposed as a less costly and more flexible production method for the valley isotopes that are representative of DT neutron fission product yield distributions, which would benefit nuclear forensics exercises in testing for analysis techniques. Natural uranium and thorium targets were irradiated with bremsstrahlung photons (nominal endpoint energies of 8, 14, and 20 MeV), generated via impingement of high-energy electrons on a tungsten electron-to-photon radiator. Gamma-ray spectroscopy was used to observe and identify the decay of the resulting fission products. The photofission yields of 102 fission products were measured for the set of the 238U and 232Th targets at the three bremsstrahlung X-ray endpoint energies 8, 14, and 20 MeV. These fragments included masses ranging from A=84 to A=144, with half-lives as short as 1.07 seconds.

07 ISOTOPE AND RADIATION SOURCES↗

Spin inhibition in γ -decay probabilities for states above S n in Sm and Dy nuclei

When a compound nucleus is formed at an excitation energy above the neutron-separation threshold, it is assumed that its de-excitation will proceed via neutron emission. However, if the excited nucleus is in a high-spin state, but does not have enough excitation energy to conserve angular momentum by either photon emission after neutron emission or relative angular momentum carried off by the neutron, the nucleus will de-excite via γ -emission instead. This effect, the spin inhibition, provides an insight into the structure of the compound nuclei and can aid the understanding of the distribution of the populated spin states in a compound nuclear reaction. In this work, the effects of spin inhibition on the γ-decay probabilities from states in 146,147Sm and 160Dy are presented. For high-spin states above the Sn, spin inhibition is able to suppress neutron emission, and de-excitation via γ -ray emission is observed for states up to 3 MeV above the Sn .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An Investigation of Time Distributions for Task Primitives to Support the HUNTER Dynamic Human Reliability Analysis

As an effort to support data collection for dynamic human reliability analysis (HRA), this study investigates time distributions for task primitives defined in the Goals, Operators, Methods, and Selection rules (GOMS)–Human Reliability Analysis (HRA) method. GOMS-HRA was developed to provide cognition-based time and human error probability information for dynamic HRA calculation in the Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER) framework. HUNTER is a framework to support the dynamic modelling of human error in conjunction with other modelling tools. In this paper, we investigate time distributions using experimental data collected from the Simplified Human Error Experimental Program (SHEEP) study, which suggests an HRA data collection framework to complement full-scope simulator research as well as collect input data for dynamic HRA using simplified simulators such as the Rancor Microworld Simulator. In this study, time required for GOMS-HRA task primitives to satisfy thirteen statistical distributions is investigated. Then, the time distributions from student operators and professional operators are compared and discussed. As a result, this study identified several time distributions on five GOMS-HRA task primitives at a statistically significant level. According to analyses to date, a greater number of significant time distributions was found in abnormal or emergency operating procedures rather than standard operating procedures. In the future, it is expected that the result of this study can provide objective reference on elapsed time data for task primitives as well as help to realistically simulate scenarios within dynamic HRA.

99 GENERAL AND MISCELLANEOUS↗

Direct numerical simulations of three-component Rayleigh–Taylor mixing and an improved model for multicomponent reacting mixtures

We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan ( J. Fluid Mech. , vol. 974, 2023, A35) and Jacobs & Dalziel ( J. Fluid Mech. , vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.

42 ENGINEERING↗

A North Atlantic synthetic tropical cyclone track, intensity, and rainfall dataset

Tropical Cyclones (TCs) cause significant socio-economic damages to the US and Caribbean coastal regions annually, making it important to understand TC risk at the local-to-regional scales. However, the short length of the observed record and the substantial computational expense associated with high-resolution climate models make it difficult to assess TC risk using either approach. To overcome these challenges, we developed a database of synthetic TCs using the Risk Analysis Framework for Tropical Cyclones (RAFT). The database includes 40,000 synthetic TC tracks, along-track intensities and storm-induced precipitation. TC tracks generated in RAFT are in reasonable agreement with the observed spatial distribution of TC tracks and basin-scale TC statistics. Specifically, along the coast, spatial variations in TC crossing probability and extreme winds upon landfall are well-reproduced by RAFT with R-squared values of 0.81 and 0.73, respectively. In summary, the synthetic TC database constructed with RAFT provides a reasonable pathway for the robust assessment of North Atlantic TC wind and rainfall risks.

54 ENVIRONMENTAL SCIENCES↗

Quantum Frequency Combs with Path Identity for Quantum Remote Sensing

Quantum sensing promises to revolutionize sensing applications by employing quantum states of light or matter as sensing probes. Photons are the clear choice as quantum probes for remote sensing because they can travel to and interact with a distant target. Existing schemes are mainly based on the quantum illumination framework, which requires quantum memory to store a single photon of an initially entangled pair until its twin reflects off a target and returns for final correlation measurements. Existing demonstrations are limited to tabletop experiments, and expanding the sensing range faces various roadblocks, including long-time quantum storage and photon loss and noise when transmitting quantum signals over long distances. We propose a novel quantum sensing framework that addresses these challenges using quantum frequency combs with path identity for remote sensing of signatures (“qCOMBPASS”). The combination of one key quantum phenomenon and two quantum resources—namely, quantum-induced coherence by path identity, quantum frequency combs, and two-mode squeezed light—allows for quantum remote sensing without requiring quantum memory. The proposed scheme is akin to a quantum radar based on entangled frequency-comb pairs that uses path identity to detect, range, or sense a remote target of interest by measuring pulses of one comb in the pair that never traveled to the target but that contains target information “teleported” by quantum-induced coherence by path identity from the other comb in the pair that traveled to the target but is not detected. We develop the basic qCOMBPASS theory, analyze the properties of the qCOMBPASS transceiver, and introduce the qCOMBPASS equation—a quantum analog of the well-known LIDAR equation in classical remote sensing. We also describe an experimental scheme to demonstrate the concept using two-mode squeezed quantum combs. qCOMBPASS can strongly impact various applications in remote quantum sensing, imaging, metrology, and communications. These applications include detection and ranging of low-reflectivity objects, measurement of small displacements of a remote target with precision beyond the standard quantum limit (SQL), standoff hyperspectral quantum imaging, discreet surveillance from space with low detection probability (detect without being detected), very-long-baseline interferometry, quantum Doppler sensing, quantum clock synchronization, and networks of distributed quantum sensors. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DNTTD (Distributed Non-Negative Tensor Train Decomposition)

The era of exascale computing opens new venues for innovations and discoveries in many scientific, engineering, and commercial fields. However, with the exa flops also come the extra-large high-dimensional data generated by high performance computing. High-dimensional data is presented as multidimensional arrays, aka tensors. The presence of latent (not directly observable) structures in the tensor allows a unique representation and compression of the data by classical tensor factorization techniques. However, the classical tensor methods are not always stable or they can be exponential in their memory requirements, which makes them not suitable for high-dimensional tensors. Tensor train (TT) is a state-of-the-art tensor network introduced for factorization of high-dimensional tensors. TT transforms the initial high-dimensional tensor in a network of three-dimensional tensors that requires only a linear storage. Many real-world data, such as, density, temperature, population, probability, etc., are non-negative and for an easy interpretation, the algorithms preserving non-negativity are preferred. Here, we introduce a distributed non-negative tensor-train and demonstrate its scalability and the compression on synthetic and real world big datasets.

Bhattarai, Manish↗

Cybersecurity Value-at-Risk Framework

As more variable renewable energy sources are added to the grid, the role of hydropower as a reliable baseline and firming resource is growing more critical. However, the U.S hydropower fleet is not fully prepared to face modern issues such as cybersecurity threats. Hydropower accounts for 37% of U.S. utility-scale renewable electricity but is challenged by diverse infrastructure and legacy devices that predate modern security practices. While new cybersecurity solutions cannot simply be added to current hydropower generation and operation technologies, custom cybersecurity assessments can reveal system-specific threats and risk probabilities and identify mitigating enhancements.

cybersecurity valuation methodology↗

Load balancing for multi-beam additive manufacturing systems

As powder bed fusion (PBF) additive manufacturing (AM) becomes a more mature field, system configurations are gradually moving away from the classic single heat source, layer-by-layer system configurations towards unconventional system configurations that offer higher throughput. Higher throughput systems allow PBF systems to be considered for a larger variety of industrial applications. However, the inclusion of multiple heat sources, or beams, also increases the complexity of the control schemes needed. For multi-beam systems with overlapping fields of view, the distribution of workload, or load balancing, across these beams directly affects the total print time for a build. Additionally, the probability of any beam failing in a multi-beam system increases with the number of beams. While manual methods of load balancing and dealing with beam failures are reasonable for current generation multi-beam systems, as system configurations become more complex, manual methods will become prohibitively inefficient. Here, this paper introduces two different ways to load balance multi-beam systems of various configuration types, regardless of their complexity, which are highly performant. A consequence of this performance is the enablement of on-the-fly load balancing in the event a beam fails, thus improving system robustness.

36 MATERIALS SCIENCE↗

Computational fluid dynamics modeling and analysis of silica nanoparticle synthesis in a flame spray pyrolysis reactor

Flame Spray Pyrolysis (FSP) is a method for large-scale production of nanoparticles and nanoscale powders employed in a wide range of industrial applications. Particle size and morphology are complex functions of the physicochemical phenomena occurring in the FSP reactor. An extensive study of FSP-related phenomena can be utilized to develop effective strategies for achieving desired particle size/morphology and scaling up the overall yield of an FSP system. In this work, a computational fluid dynamics (CFD) model of an FSP reactor is developed to simulate the coupling of key phenomena involved in the particle synthesis process: liquid spray breakup and evaporation, mixing, combustion, and particle formation/growth of silica nanoparticles. Herein, the particle sizes and their distributions from the CFD simulations are validated against experimental data. Subsequently, the simulations are utilized to investigate the impact of process parameters on the resultant flame dynamics and particle growth. Firstly, the CFD results show that the particle sizes are strongly correlated with the precursor concentration in the solvent. At lower precursor concentrations, the spread of the distribution is relatively insensitive to the value of the concentration. At higher concentrations, the spread is higher as the collision probability between particles is higher. Secondly, increasing the pilot flow rate increases the length of the pilot flames impacting the local ignition location of the spray flame. Lastly, it is shown that the dispersion gas flow rate strongly influences the spray flame shape. This shape can be used for control of particle growth as it helps determine the regions of high temperature and the residence time of the particles in the high temperature region enabling the design and process optimization of the FSP reactor.

42 ENGINEERING↗

Statistical analysis of non-Maxwellian electron distribution functions measured with angularly resolved Thomson scattering

Angularly resolved Thomson scattering is a novel extension of Thomson scattering, enabling the measurement of the electron velocity distribution function over many orders of magnitude. Here, details of the theoretical basis of the technique and the instrument designed for this measurement are described. Angularly resolved Thomson-scattering data from several experiments are shown with descriptions of the corresponding distribution functions. A reduced model describing the distribution function is given and used to perform a Monte Carlo analysis of the uncertainty in the measurements. The electron density and temperature were determined to a precision of 12% and 21%, respectively, on average, while all other parameters defining the distribution function were generally determined to better than 20%. It was found that these uncertainties were primarily due to limited signal to noise and instrumental effects. Furthermore, measurements with this level of precision were sufficient to distinguish between Maxwellian and non-Maxwellian distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimizing temperature distributions for training neural quantum states using parallel tempering

Parametrized artificial neural networks (ANNs) can be very expressive ansatzes for variational algorithms, reaching state-of-the-art energies on many quantum many-body Hamiltonians. Nevertheless, the training of the ANN can be slow and stymied by the presence of local minima in the parameter landscape. One approach to mitigate this issue is to use parallel tempering methods, and in this work, we focus on the role played by the temperature distribution of the parallel tempering replicas. Using an adaptive method that adjusts the temperatures in order to equate the exchange probability between neighboring replicas, we show that this temperature optimization can significantly increase the success rate of the variational algorithm with negligible computational cost by eliminating bottlenecks in the replicas' random walk. Furthermore, we demonstrate this using two different neural networks, a restricted Boltzmann machine and a feedforward network, which we use to study a toy problem based on a permutation invariant Hamiltonian with a pernicious local minimum and the 𝐽 1 −𝐽 2 model on a rectangular lattice.

Neural network simulations↗

Fast HARDI Uncertainty Quantification and Visualization with Spherical Sampling

In this paper, we study uncertainty quantification and visualization of orientation distribution functions (ODF), which corresponds to the diffusion profile of high angular resolution diffusion imaging (HARDI) data. The shape inclusion probability (SIP) function is the state‐of‐the‐art method for capturing the uncertainty of ODF ensembles. The current method of computing the SIP function with a volumetric basis exhibits high computational and memory costs, which can be a bottleneck to integrating uncertainty into HARDI visualization techniques and tools. We propose a novel spherical sampling framework for faster computation of the SIP function with lower memory usage and increased accuracy. In particular, we propose direct extraction of SIP isosurfaces, which represent confidence intervals indicating spatial uncertainty of HARDI glyphs, by performing spherical sampling of ODFs. Our spherical sampling approach requires much less sampling than the state‐of‐the‐art volume sampling method, thus providing significantly enhanced performance, scalability, and the ability to perform implicit ray tracing. Our experiments demonstrate that the SIP isosurfaces extracted with our spherical sampling approach can achieve up to 8164× speedup, 37282× memory reduction, and 50.2% less SIP isosurface error compared to the classical volume sampling approach. We demonstrate the efficacy of our methods through experiments on synthetic and human‐brain HARDI datasets.

97 MATHEMATICS AND COMPUTING↗