Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Probability density function”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Probabilistic structural analysis methods for improving Space Shuttle engine reliability

Probabilistic structural analysis methods are particularly useful in the design and analysis of critical structural components and systems that operate in very severe and uncertain environments. These methods have recently found application in space propulsion systems to improve the structural reliability of Space Shuttle Main Engine (SSME) components. A computer program, NESSUS, based on a deterministic finite-element program and a method of probabilistic analysis (fast probability integration) provides probabilistic structural analysis for selected SSME components. While computationally efficient, it considers both correlated and nonnormal random variables as well as an implicit functional relationship between independent and dependent variables. The program is used to determine the response of a nickel-based superalloy SSME turbopump blade. Results include blade tip displacement statistics due to the variability in blade thickness, modulus of elasticity, Poisson's ratio or density. Modulus of elasticity significantly contributed to blade tip variability while Poisson's ratio did not. Thus, a rational method for choosing parameters to be modeled as random is provided.

Boyce, L.↗
DOE OSTI2024

Density estimation via measure transport: Outlook for applications in the biological sciences

Abstract One among several advantages of measure transport methods is that they allow or a unified framework for processing and analysis of data distributed according to a wide class of probability measures. Within this context, we present results from computational studies aimed at assessing the potential of measure transport techniques, specifically, the use of triangular transport maps, as part of a workflow intended to support research in the biological sciences. Scenarios characterized by the availability of limited amount of sample data, which are common in domains such as radiation biology, are of particular interest. We find that when estimating a distribution density function given limited amount of sample data, adaptive transport maps are advantageous. In particular, statistics gathered from computing series of adaptive transport maps, trained on a series of randomly chosen subsets of the set of available data samples, leads to uncovering information hidden in the data. As a result, in the radiation biology application considered here, this approach provides a tool for generating hypotheses about gene relationships and their dynamics under radiation exposure.

gene expression data↗
DOE OSTI2020

DFT Plus DMFT Study of Spin-Charge-Lattice Coupling in Covalnt LaCoO{siub 3}.

We study energetics and the nature of both homogeneous and mixed spin (MS) states in LaCoO3 incorporating structural changes of the crystal volume expansion and the Co-O bond disproportionation (BD) during the spin-state transition using the density functional theory plus dynamical mean field theory (DFT + DMFT) method. DFT + DMFT predicts that energetics of both excited spin states are almost the same while DFT + U calculations of the same structures energetically favor the MS states and produce various metastable solutions whose energetics depend sensitively on final spin states. Within DFT + DMFT, the homogeneous spin state in the expanded crystal volume shows the multiconfigurational nature with non-negligible occupancy probabilities of both high spin (HS) and low spin (LS) states along with d(6) and d(7) charge configurations indicating the dynamically fluctuating nature of spin and charge states due to the Co-O covalency. The nature of the MS state under the BD structure reveals that Co sites with the long Co-O bonds develop a Mott insulating state and favor HS with a d(6) configuration, while more covalent Co sites with the short Co-O bonds occupy more LS states with a d(7) configuration and behave as a band insulator, as a result, charge ordering is induced in the BD structure from the spin-state ordering. We also find that both energetics and electronic structure sensitively depend on the Co-O covalency effect, which can be tuned by changing the double counting potential and the resulting d occupancy (N-d), and N-d close to 6.7 is consistent with the nature of the spin-state transition. Our results show that structural changes during the spin-state transition can play an important role in understanding energetics and electronic structure of LaCoO3.

Park, Hyowon↗
DOE OSTI2023

The DESI One-Percent Survey: Constructing Galaxy–Halo Connections for ELGs and LRGs Using Auto and Cross Correlations

In the current Dark Energy Spectroscopic Instrument (DESI) survey, emission line galaxies (ELGs) and luminous red galaxies (LRGs) are essential for mapping the dark matter distribution at z ~ 1. We measure the auto and cross correlation functions of ELGs and LRGs at 0.8 < z ≤ 1.0 from the DESI One-Percent survey. Following Gao et al., we construct the galaxy–halo connections for ELGs and LRGs simultaneously. With the stellar–halo mass relation for the whole galaxy population (i.e., normal galaxies), LRGs can be selected directly by stellar mass, while ELGs can also be selected randomly based on the observed number density of each stellar mass, once the probability P sat of a satellite galaxy becoming an ELG is determined. We demonstrate that the observed small scale clustering prefers a halo mass-dependent P sat model rather than a constant. With this model, we can well reproduce the auto correlations of LRGs and the cross correlations between LRGs and ELGs at r p > 0.1 Mpc h –1 . We can also reproduce the auto correlations of ELGs at r p > 0.3 Mpc h –1 (s > 1 Mpc h –1 ) in real (redshift) space. Although our model has only seven parameters, we show that it can be extended to higher redshifts and reproduces the observed auto correlations of ELGs in the whole range of 0.8 < z ≤ 1.6, which enables us to generate a lightcone ELG mock for DESI. With the above model, we further derive halo occupation distributions for ELGs, which can be used to produce ELG mocks in coarse simulations without resolving subhalos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Database-wide hazard modelling of the onset of DIII-D tearing modes with field features

The rate of onset (hazard) of tearing modes is modelled probabilistically using statistical learning algorithms. Axisymmetric energy-density equilibrium fields are taken as raw high-dimensional input features which are reduced with principal component analysis. Signal processing of non-axisymmetric magnetics fluctuation array data provides the target information from which to learn. Model selection, visualization and calibration assessment procedures are detailed. Here, the analysis is deployed at large scale across the DIII-D tokamak database. Standard model selection criteria suggest that the energy-density post-processed feature is a better choice for modelling the onset rate compared to the non-processed equilibrium reconstruction solution. Two example applications of the learned rate function are demonstrated: (i) proximity-to-onset discharge monitoring and (ii) database analysis showing an (expected) observational global trend that the general hazard increases as a plasma performance metric increases. An important connection between the hazard function and its use as a conditional probability generator is reviewed in the Appendix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

Spatial scale of motion segmentation from speed cues

For the accurate perception of multiple, potentially overlapping, surfaces or objects, the visual system must distinguish different local motion vectors and selectively integrate similar motion vectors over space to segment the retinal image properly. We recently showed that large differences in speed are required to yield a percept of motion transparency. In the present study, to investigate the spatial scale of motion segmentation from speed cues alone, we measured the speed-segmentation threshold (the minimum speed difference required for 75% performance accuracy) for 'corrugated' random-dot patterns, i.e. patterns in which dots with two different speeds were alternately placed in adjacent bars of variable width. In a first experiment, we found that, at large bar widths, a smaller speed difference was required to segment and perceive the corrugated pattern of moving dots, while at small bar-widths, a larger speed difference was required to segment the two speeds and perceive two transparent surfaces of moving dots. Both the perceptual and segmentation performance transitions occurred at a bar width of around 0.4 degrees. In a second experiment, speed-segmentation thresholds were found to increase sharply when dots with different speeds were paired within a local pooling area. The critical pairing distance was about 0.2 degrees in the fovea and increased linearly with stimulus eccentricity. However, across the range of eccentricities tested (up to 15 degrees ), the critical pairing distance did not change much and remained close to the receptive field size of neurons within the primate primary visual cortex. In a third experiment, increasing dot density changed the relationship between speed-segmentation thresholds and bar width. Thresholds decreased for large bar widths, but increased for small bar widths. All of these results are well fit by a simple stochastic model, which estimates the probabilities of having identical or different motion vectors within a local pooling area whose size is the same as that of primate V1 neurons. Altogether, these results demonstrate that speed-based segmentation can function well, even at small spatial scales (i.e. high-spatial frequencies of spatial corrugation) and thereby emphasizes the critical role of a local pooling process early in the cortical motion-processing pathway.

NASA Discipline Neuroscience↗
NASA NTRS1986

Transition probability of the Al II 2669 intersystem line

Time-resolved observations of the spin-changing, or 'intersystem' emission at 2669.157 A obtained by the ion storage technique are used to measure the transition probability of the 3s2 1S0 - 3s3p 3P1 exp 0 line in Al II. A laser-generated plasma was used as the source of the metastable Al(+) ions. The A-value result obtained for the intersystem transition is 3.33 + or - 0.23 x 10 to the 3rd/sec at the 90-percent confidence level; this value is used to derive two line-intensity ratios which involve the intersystem line as a function of electron density and temperature.

Johnson, B. C.↗
NASA NTRS1993

Sharp edges to neutral hydrogen disks in galaxies and the extragalactic radiation field

It is shown that the very sharp truncation of the neutral hydrogen distribution seen in NGC 3198 (and probably M33) is well modeled as the result of ionization of the atomic gas by the extragalactic radiation field. Below a critical column density of about a few times 10 exp 19/sq cm the gas is dominantly ionized and undetectable in the 21-cm line. It is inferred from the photoionization models that the total disk gas distribution in NGC 3198 is actually fairly axisymmetric. The critical column density for ionization is not a strong function of galaxy mass or mass distribution; thus, all galaxies should show a cutoff at approximately the same column density. Specific models of 3198 suggest that the extragalactic ionizing photon flux is 5000-10,000 photons/sq cm s.

Maloney, Philip↗
NASA NTRS1999

A Distant, X-Ray Luminous Cluster of Galaxies at Redshift 0.83

We have observed the most distant (= 0.829) cluster of galaxies in the Einstein Extended Medium Sensitivity Survey (EMSS), with the ASCA and ROSAT satellites. We find an X-ray temperature of 12.3(sup 3.1, sub 2.2) keV for this cluster, and the ROSAT map reveals significant substructure. The high temperature of MS1054-0321 is consistent with both its approximate velocity dispersion, based on the redshifts of 12 cluster members we have obtained at the Keck and the Canada-France-Hawaii telescopes, and with its weak lensing signature. The X-ray temperature of this cluster implies a virial mass approximately 7.4 x 10(exp 14) /h solar mass, if the mean matter density in the universe equals the critical value (OMEGA(sub 0) = 1), or larger if OMEGA(sub 0) < 1. Finding such a hot, massive cluster in the EMSS is extremely improbable if clusters grew from Gaussian perturbations in an OMEGA(sub 0) = 1 universe. Combining the assumptions that OMEGA(sub 0) = 1 and that the initial perturbations were Gaussian with the observed X-ray temperature function at low redshift, we show that this probability of this cluster occurring in the volume sampled by the EMSS is less than a few times 10(exp -5). Nor is MS1054-0321 the only hot cluster at high redshift; the only two other z > 0.5 EMSS clusters already observed with ASCA also have temperatures exceeding 8 keV. Assuming again that the initial perturbations were Gaussian and OMEGA(sub 0) = 1, we find that each one is improbable at the < 10(exp -2) level. These observations, along with the fact that these luminosities and temperatures of the high-z clusters all agree with the low-z L(sub x) - T(sub x) relation, argue strongly that OMEGA(sub 0) < 1. Otherwise, the initial perturbations must be non-Gaussian, if these clusters' temperatures do indeed reflect their gravitational potentials.

Donahue, Megan↗
NASA NTRS2014

Ballistic Performance Model of Crater Formation in Monolithic, Porous Thermal Protection Systems

Porous monolithic ablative systems insulate atmospheric reentry vehicles from reentry plasmas generated by atmospheric braking from orbital and exo-orbital velocities. Due to the necessity that these materials create a temperature gradient up to several thousand Kelvin over their thickness, it is important that these materials are near their pristine state prior to reentry. These materials may also be on exposed surfaces to space environment threats like orbital debris and meteoroids leaving a probability that these exposed surfaces will be below their prescribed values. Owing to the typical small size of impact craters in these materials, the local flow fields over these craters and the ablative process afford some margin in thermal protection designs for these locally reduced performance values. In this work, tests to develop ballistic performance models for thermal protection materials typical of those being used on Orion are discussed. A density profile as a function of depth of a typical monolithic ablator and substructure system is shown in Figure 1a.

Miller, J. E.↗
DOE OSTI2025

Superspin renormalization and slow relaxation in random spin systems

We develop an excited-state real-space renormalization group (RSRG-X) formalism to describe the dynamics of conserved densities in randomly interacting spin-12 systems. Our formalism is suitable for systems with U(1) and Z2 symmetries, and we apply it to chains of randomly positioned spins with dipolar XX+YY interactions, as arise in Rydberg quantum simulators and other platforms. The formalism generates a sequence of effective Hamiltonians that provide approximate descriptions for dynamics on successively smaller energy scales. These effective Hamiltonians involve “superspins”: two-level collective degrees of freedom constructed from (anti)aligned microscopic spins. Conserved densities can then be understood as relaxing via coherent collective spin flips. For the well-studied simpler case of randomly interacting nearest-neighbor XX+YY chains, the superspins reduce to single spins. Our formalism also leads to a numerical method capable of simulating the dynamics up to an otherwise inaccessible combination of large system size and late time. Focusing on disorder-averaged infinite-temperature autocorrelation functions, in particular the spin survival probability Sp¯(t), we demonstrate quantitative agreement between our algorithm and exact diagonalization (ED) at low but nonzero frequencies. Such agreement holds for chains with nearest-neighbor, next-nearest-neighbor, and long-range dipolar interactions. Our results indicate decay of Sp¯(t) slower than any power law and feature no significant deviation from the ∼1/ln2(t) asymptote expected from the infinite-randomness fixed-point of the nearest-neighbor model. We also apply the RSRG-X formalism to two-dimensional long-range systems of moderate size and find slow late-time decay of Sp¯(t).

Zhao, Yi J↗
DOE OSTI2020

Bottom-Up Meets Top-Down: The Crossroads of Multiscale Chromatin Modeling

Chromatin can be viewed as a hierarchically structured fiber that regulates gene expression. It consists of a complex network of DNA and proteins whose characteristic dynamical modes facilitate compaction and rearrangement in the cell nucleus. These modes stem from chromatin's fundamental unit, the nucleosome, and their effects are propagated across length scales. Understanding the effects of nucleosome dynamics on the chromatin fiber, primarily through post-translational modifications that occur on the histones, is of central importance to epigenetics. Within the last decade, imaging and chromosome conformation capture techniques have revealed a number of structural and statistical features of the packaged chromatin fiber at a hitherto unavailable level of resolution. Such experiments have led to increased efforts to develop polymer models that aim to reproduce, explain, and predict the contact probability scaling and density heterogeneity. At nanometer scales, available models have focused on the role of the nucleosome and epigenetic marks on local chromatin structure. At micrometer scales, existing models have sought to explain scaling laws and density heterogeneity. Less work, however, has been done to reconcile these two approaches: bottom-up and top-down models of chromatin. Here, we highlight the multiscale simulation models that are driving toward an understanding of chromatin structure and function, from the nanometer to the micron scale, and we highlight areas of opportunity and some of the prospects for new frameworks that bridge these two scales. Taken together, experimental and modeling advances over the last few years have established a robust platform for the study of chromatin fiber structure and dynamics, which will be of considerable use to the chromatin community in developing an understanding of the interplay between epigenomic regulation and molecular structure.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Assessing density functionals for describing methane dissociative chemisorption on Pt(110)-(2×1) surface

In this work, we explore the suitability of several density functionals with the generalized gradient approximation (GGA) and beyond for describing the dissociative chemisorption of methane on the reconstructed Pt(110)-(2×1) surface. The bulk and surface structures of the metal, methane adsorption energy, and dissociation barrier are used to assess the functionals. A van der Waals corrected GGA functional (optPBE-vdW) and a metaGGA functional with van der Waals correction (MS PBEl-rVV10) are selected for ab initio molecular dynamics calculations of the sticking probability. Our results suggest that the use of these two functionals may lead to a better agreement with existing experimental results, thus serving as a good starting point for future development of reliable machine-learned potential energy surfaces for the dissociation of methane on the Pt(110)-(2×1) surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Measures of complexity and entanglement in many-fermion systems

There is no unique and widely accepted definition of the complexity measure (CM) of a many-fermion wave function in the presence of interactions. The simplest many-fermion wave function is a Slater determinant. In shell-model or configuration interaction (CI) and other related methods, the state is represented as a superposition of a large number of Slater determinants, which in the case of CI calculations reaches about 20 billion terms [Johnson, arXiv:1809.07869]. Although in practice this number has been used as a CM for decades, it is ill defined: it is not unique, and it depends on the particular type and the number of single-particle wave functions used to construct the Slater determinants. Further, the canonical wave functions and/or natural orbitals [Löwdin, Adv. Phys. 5, 1 (1956); Löwdin and Shull, Phys. Rev. 101, 1730 (1956); Bardeen et al., Phys. Rev. 108, 1175 (1957); N. N. Bogoljubov, Il Nuovo Cimento 7, 794 (1958); Valatin, Il Nuovo Cimento 7, 843 (1958); de Gennes, Superconductivity of Metals and Alloys (CRC Press, Boca Raton, FL, 1999); Ring and Schuck, The Nuclear Many-Body Problem, 1st ed. (Springer-Verlag, Berlin, 2004)] and their corresponding occupation probabilities are intrinsic properties of any many-body wave function, irrespective of the representation, and they provide a unique solution to characterize the CM. The non-negative orbital entanglement entropy, which vanishes for a Slater determinant, provides the simplest CM, while a more complete measure of complexity is the entanglement spectrum. We illustrate these aspects in the case of a complex nonequilibrium time-dependent process, induced nuclear fission described within a real-time density functional theory framework extended to superfluid systems, which can describe simultaneously the long-range and the short-range correlations between fermions. The orbital entanglement entropy of the fissioning nucleus illustrates the localization mechanism of the many-body wave function in Fock and/or Hilbert space. The (minimal) number of Slater determinants required to represent such a complex many-body wave function with a well-defined number of particles in the case presented here is about 10 500 . The realistic case of the highly nonequilibrium nuclear fission process illustrated here is equivalent to a system of 23.328×10 9 interacting quantum spin-1/2 particles, a very large system for the study of quantum entanglement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1998

A Very Hot, High Redshift Cluster of Galaxies: More Trouble for Omega(0) = 1

We have observed the most distant (= 0.829) cluster of galaxies in the Einstein Extended Medium Sensitivity Survey (EMSS), with the ASCA and ROSAT satellites. We find an X-ray temperature of 12.3 (sup +3.1) (sub -2.2)keV for this cluster, and the ROSAT map reveals significant substructure. The high temperature of MS1054-0321 is consistent with both its approximate velocity dispersion, based on the redshifts of 12 cluster members we have obtained at the Keck and the Canada-France-Hawaii telescopes, and with its weak lensing signature. The X-ray temperature of this cluster implies a virial mass approx. 7.4 x 10 (sup 14) h (sup -1) M (circle dot), if the mean matter density in the universe equals the critical value (OMEGA (sub 0) = 1), or larger if OMEGA (sub 0) is less than 1. Finding such a hot, massive cluster in the EMSS is extremely improbable if clusters grew from Gaussian perturbations in an OMEGA (sub 0) = 1 universe. Combining the assumptions that OMEGA (sub 0) = 1 and that the initial perturbations were Gaussian with the observed X-ray temperature function at low redshift, we show that this probability of this cluster occurring in the volume sampled by the EMSS is less than a few times 10 (sup -5). Nor is MS1054-0321 the only hot cluster at high redshift; the only two other z greater than 0.5 EMSS clusters already observed with ASCA also have temperatures exceeding 8 keV. Assuming again that the initial perturbations were Gaussian and OMEGA (sub 0) = 1, we find that each one is improbable at the less than 10 (sup -2) level. These observations, along with the fact that these luminosities and temperatures of the high-z clusters all agree with the low-z L (sub X) - T (sub X) relation, argue strongly that OMEGA (sub 0) less than 1. Otherwise, the initial perturbations must be non-Gaussian, if these clusters' temperatures do indeed reflect their gravitational potentials.

Donahue, Megan↗
DOE OSTI2021

Single-Molecule Kinetics of Styrene Hydrogenation on Silica-Supported Vanadium: The Role of Disorder for Single-Atom Catalysts

We report a theoretical approach for the study of supported atom catalysis is developed based on recent advances in the study of single-molecule kinetics. This view is particularly useful in exhibiting the role of disorder in single-atom and single-site catalysts on amorphous supports. The distribution of passage times (or waiting times) through a complex catalytic network originating from a set of coupled active sites is described by a probability distribution function (PDF), f(t), that reflects the local environment of the reaction center. An efficient algorithm is developed based on the linear algebra of the Markov transition matrix that produces f(t) or its moments. The kinetics of the hydrogenation reaction of styrene on an organovanadium(III) catalyst supported on amorphous silica is studied. A kinetic model consisting of three intertwined catalytic cycles emanating from three chemically distinct active sites is proposed to describe the chemistry. Density functional theory (DFT) calculations are employed to determine the free energy barriers of the reactions, which are used to construct the rate coefficient matrix. The disorder induced by the amorphous support material is divided into a low-dimensional short-range component reflecting the covalent structures near the reaction center and a weaker long-range component modeling the bulk randomness. The results are computed and analyzed for a wide range of concentration values and disorder scenarios. The unusual structure in the f(t) PDF is found to occur for certain cases that reveal the contribution of multiple catalytic pathways acting in concert.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1984

The cyanopolyynes as a chemical clock for molecular clouds

A new method is proposed for determining the ages of molecular clouds. The method utilizes the properties of the long-chain organic molecules known as the cyanopolyynes (HC/2k + 1/N, k = 0, 1, 2, 3 ...), which are found in a variety of clouds. The observed regular abundance decline of these molecules as a function of chain length suggests that they are formed sequentially. If so, the age of the cloud can be read off as the time to grow the longest chain present. Although the creation process for the chains is still unknown, the age can be obtained from knowledge of the chain destruction mechanism. Destruction is probably due to adsorption onto the surfaces of interstellar grains. Using the known properties of grains, the age can be obtained from the cloud density and the abundance ratios of the cyanopolyynes. As a first application, minimum ages for the four dark clouds B335, TMC-1, TMC-2, and L183 are obtained. These minimum ages are 1.3 x 10 to the 6th yr, 9.7 x 10 to the 5th yr, 3.4 x 10 to the 5th yr, and 3.3 x 10 to the 5th yr, respectively. In each case, the ages are greater than or equal to the cloud's free-fall collapse time. These four clouds are therefore either in a state of hydrostatic balance or have only recently begun to collapse, following an earlier period of hydrostatic support.

Stahler, S. W.↗
DOE OSTI2023

Molecular Dynamics Simulation of Hygroscopic Aging Effects in Epoxy Polymer

The automobile industry is incorporating more lightweight content in car designs to boost fuel-economy. New structural adhesives are needed to mitigate the corrosion and thermal expansion issues associated with joining dissimilar lightweight materials, but adhesive developers lack a fundamental understanding of the chemistry that occurs in the adhesive as the joint ages. In this study, we developed structural adhesive molecular models and applied classical molecular dynamics simulations and density functional theory calculations to gain molecular insights into the influence of water molecules on the properties of epoxy-based adhesives (DGEBA + Jeffamine (JD230)). The simulations were complemented by experimental synthesis and characterization. Our work underscores the impact of water molecules on the local structure of the epoxy network as well as resulting mechanical properties. Water molecules were mainly coordinated with hydroxyls, primary amines and secondary amines, but also weakly coordinated with ether linkages, which were found most probable to be labile. Simulated stress–strain data indicates that increasing the water content deteriorates the mechanical properties. The Young’s modulus decreased by ~ 30% when the water content increased to 3 wt%. We conclude, this integration of molecular-level chemical insights with mechanical property simulations of the hydrated epoxy system and experimental validation holds the promise to advance lightweight joint technologies.

36 MATERIALS SCIENCE↗