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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Identification of hot water end-use process of electric water heaters from energy measurements

This paper presents an algorithm for the identification of parameters for a stochastic hot water end-use process that drives a homogeneous population of thermostatically controlled electric water heaters (EWH). Usually, only metered interval consumption data (kWh) is collected and the hot water end-use process is unobservable to utility and aggregators. However, the availability of EWHs for demand response (DR) is closely coupled with the hot water end-use process. In this context, the hot water end-use process is modeled as a two-state Markov chain (Use / No use), which causes the thermostatic ON-OFF switching process to behave as a Markov renewal process (MRP). A set of first passage-time problems is developed to obtain the moments of the transition probability densities of the MRP. These problems are addressed by establishing a system of coupled partial differential equations characterizing the temperature evolution of the EWH population. A key quantity in the methodology for estimating the parameters is the total time an EWH is ON within a period of interest. It is referred to as the total busy time. Total busy time in this approach is a random variable for which analytical expressions of the moments are developed as a function of the metered window length. The latter expressions become the basis of a hot water demand model identification algorithm which is validated using agent-based simulations of EWHs.

42 ENGINEERING↗
DOE OSTI2025

Uncovering the True Active Sites in Ni–N–C Catalysts for CO 2 Electroreduction

Understanding and designing active sites in single-atom catalysts (SACs) requires going beyond static models to capture their dynamic evolution under realistic electrochemical conditions. Here, in this work, we develop an integrated theoretical framework that accounts for operational conditions, by combining grand canonical density functional theory (GC-DFT) with machine-learning-accelerated sampling, to uncover structure–activity–stability relationships in Ni–N–C SACs for the CO 2 reduction reaction (CO 2 RR). A library of NiN x C 4–x (x = 0–4) motifs─representing coordination defects likely formed during high-temperature synthesis─was systematically evaluated. Under working conditions, these sites were found to undergo hydrogenation, and NiN 3 C 1_ H 1 was identified as the most probable active site. At reducing potentials, hydrogen adsorbs spontaneously at C–Ni bridge sites rather than Ni top sites, while subsurface hydrogen facilitates bent CO 2 adsorption crucial for activation. High CO 2 RR selectivity toward CO arises from site separation: Ni centers drive CO2RR, while the hydrogen evolution reaction (HER) occurs at the C–Ni bridge or N sites and from thermodynamic suppression of HER at moderate hydrogen coverage. At more negative potentials, a shift in the CO 2 RR rate-determining process (RDP) and Ni out-of-surface displacement induced by coadsorption of H and H 2 O jointly reduce activity and selectivity. Thus, both the high CO2RR selectivity of Ni–N–C catalysts and its reversal with more negative potentials can be rationalized by accounting for hydrogenated surfaces. This highlights the necessity of modeling realistic; in situ conditions. This framework provides generalizable insights into the dynamic behavior of active sites in SACs, offering guidance for the rational design of active and robust catalysts for a wide range of electrochemical reactions.

25 ENERGY STORAGE↗
DOE OSTI2024

On the statistical theory of self-gravitating collisionless dark matter flow: Scale and redshift variation of velocity and density distributions

The statistics of velocity and density fields are crucial for cosmic structure formation and evolution. Here, this paper extends our previous work on the two-point second-order statistics for the velocity field [Phys. Fluids 35, 077105 (2023)] to one-point probability distributions for both density and velocity fields. The scale and redshift variation of density and velocity distributions are studied by a halo-based non-projection approach. First, all particles are divided into halo and out-of-halo particles so that the redshift variation can be studied via generalized kurtosis of distributions for halo and out-of-halo particles, respectively. Second, without projecting particle fields onto a structured grid, the scale variation is analyzed by identifying all particle pairs on different scales $r$. We demonstrate that: (i) Delaunay tessellation can be used to reconstruct the density field. The density correlation, spectrum, and dispersion functions were obtained, modeled, and compared with the N-body simulation; (ii) the velocity distributions are symmetric on both small and large scales and are non-symmetric with a negative skewness on intermediate scales due to the inverse energy cascade on small scales with a constant rate $\varepsilon_u$; (iii) On small scales, the even order moments of pairwise velocity $\Delta u_L$ follow a two-thirds law $\propto{(-\varepsilon_ur)}^{2/3}$, while the odd order moments follow a linear scaling $\langle(\Delta u_L)^{2n+1}\rangle=(2n+1)\langle(\Delta u_L)^{2n}\rangle\langle\Delta u_L\rangle\propto{r}$; (iv) The scale variation of the velocity distributions was studied for longitudinal velocities $u_L$ or $u_L^{'}$, pairwise velocity (velocity difference) $\Delta u_L$=$u_L^{'}$-$u_L$ and velocity sum $\Sigma u_L$=$u^{'}_L$+$u_L$. Fully developed velocity fields are never Gaussian on any scale, despite that they can initially be Gaussian; (v) On small scales, $u_L$ and $\Sigma u_L$ can be modeled by a $X$ distribution to maximize the entropy of the system. The distribution of $\Delta u_L$ can be different; (vi) On large scales, $\Delta u_L$ and $\Sigma u_L$ can be modeled by a logistic or a $X$ distribution, while $u_L$ has a different distribution; (vii) the redshift variation of the velocity distributions follows the evolution of the $X$ distribution involving a shape parameter $\alpha(z)$ decreasing with time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Optimizing Earthquake Nowcasting With Machine Learning: The Role of Strain Hardening in the Earthquake Cycle

Abstract Nowcasting is a term originating from economics, finance, and meteorology. It refers to the process of determining the uncertain state of the economy, markets or the weather at the current time by indirect means. In this paper, we describe a simple two‐parameter data analysis that reveals hidden order in otherwise seemingly chaotic earthquake seismicity. One of these parameters relates to a mechanism of seismic quiescence arising from the physics of strain‐hardening of the crust prior to major events. We observe an earthquake cycle associated with major earthquakes in California, similar to what has long been postulated. An estimate of the earthquake hazard revealed by this state variable time series can be optimized by the use of machine learning in the form of the Receiver Operating Characteristic skill score. The ROC skill is used here as a loss function in a supervised learning mode. Our analysis is conducted in the region of 5° × 5° in latitude‐longitude centered on Los Angeles, a region which we used in previous papers to build similar time series using more involved methods (Rundle & Donnellan, 2020, https://doi.org/10.1029/2020EA001097 ; Rundle, Donnellan et al., 2021, https://doi.org/10.1029/2021EA001757 ; Rundle, Stein et al., 2021, https://doi.org/10.1088/1361-6633/abf893 ). Here we show that not only does the state variable time series have forecast skill, the associated spatial probability densities have skill as well. In addition, use of the standard ROC and Precision (PPV) metrics allow probabilities of current earthquake hazard to be defined in a simple, straightforward, and rigorous way.

58 GEOSCIENCES↗
DOE OSTI2024

Targeted Adaptive Design

Modern advanced manufacturing and advanced materials design often require searches of relatively high-dimensional process control parameter spaces for settings that result in optimal structure, property, and performance parameters. The mapping from the former to the latter must be determined from noisy experiments or from expensive simulations. Here, we abstract this problem to a mathematical framework in which an unknown function from a control space to a design space must be ascertained by means of expensive noisy measurements, which locate control settings generating desired design features within specified tolerances, with quantified uncertainty. We describe targeted adaptive design (TAD), a new algorithm that performs this sampling task efficiently. TAD creates a Gaussian process surrogate model of the unknown mapping at each iterative stage, proposing a new batch of control settings to sample experimentally and optimizing the updated expected log-predictive probability density of the target design. TAD either stops upon locating a solution with uncertainties that fit inside the tolerance box or uses a measure of expected future information to determine that the search space has been exhausted with no solution. TAD thus embodies the exploration-exploitation tension in a manner that recalls, but is essentially different from, Bayesian optimization and optimal experimental design.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

A BBGKY framework for fluid turbulence

A framework is presented for a systematic kinetic theory of turbulence originating from the Liouville equation for the Fourier coefficients of fluid variables. The real and imaginary parts of these Fourier coefficients play the role that particle coordinates (positions and momenta) play in the BBGKY theory. The basic relations of the problem are the incompressible Navier-Stokes equations in two dimensions with zero viscosity, with the probability distributions of Fourier coefficients rather than moments being the basic variables of the theory. A kinetic equation is derived and shown to possess a number of requirements that any reasonable kinetic equation must have: conservation laws, positive-definite spectral densities, and an H-theorem. The major lack in the theory is any reliable information on the relaxation predicted by the complicated linear operator H. Closure of the hierarchy is achieved by the hypothesis that the five-coefficient correlation function is negligible. Problems associated with inclusion of viscosity and external driving forces are discussed.

Montgomery, D.↗
NASA NTRS1981

Evaluation and universal curves of the photon loss probability

In the transfer of radiation in a plasma, an important parameter (which measures the coupling between the radiation field and the electron gas) is the photon loss probability per scattering epsilon, defined as the ratio of collisional de-excitation to total de-excitation. The evaluation of epsilon requires normally knowledge of specific collisional excitation rates and radiative transition rate coefficients. In the considered investigation a general formula for epsilon is obtained which depends only on electron density, temperature, and line wavelength. Use is made of a simple Bethe-Coulomb expression recently derived by Kastner (1980) for the collisional excitation rate coefficient. Epsilon is obtained as a function of the three variables, and universal curves of it can be constructed with one of the variables as parameter. As an example of the application of the considered procedure, the variation of epsilon with height h in the solar atmosphere is shown in a graph for several solar lines.

Kastner, S. O.↗
NASA NTRS1991

Investigation of loss processes of Tm and Tm,Ho in YAG

The loss of excitation from various manifolds of Tm and Tm,Ho in YAG as a function of temperature and concentration is studied. Two probable loss mechanisms - a Tm up-conversion and a Ho up-conversion - are identified. A 785-nm CW diode laser with 400-nW peak power was focused to a small spot on the sample. The emission from the sample observed at 90 deg was monitored through a monochromator with slits open to 3 mm. Intensity of emission was measured by varying the power of the excitation source using a set of neutral density filters. Power is reported as the percentage of the peak power, and the intensity curves were normalized below 20 percent of transmission. The fact that there is emission above the pump energy indicates an up-conversion from excited manifolds. Nonlinear changes in the intensity of the emission from the Tm 3F4 manifold with the pump power reveals a loss of excitation from this manifold. The linear dependence of the 5I7 manifold emission with pump power at low Tm and high Ho concentrations and the gain of energy in the 5I6 manifold of Ho indicate that the 5I7 manifold loss is due to the coupling of Tm and Ho ions.

Armagan, G.↗
DOE OSTI2021

First-principles studies of the concentration-dependent tritium diffusion in the zirconium hydrides with and without Sn impurity

We performed first-principles calculations to study the 3 H diffusion in the zirconium hydrides with or without Sn impurity. Our results show that the formation of Zr 3 H x becomes preferable as the 3 H concentration increases. Based on the most stable configurations of Zr 3 H x (x = 0.5, 1.0, 1.5, 2.0), we studied the interstitial 3 H diffusion with or without Sn impurity. The results show that the interstitial 3 H diffusion becomes less probable in the hydrides with a higher 3 H concentration. When introducing a substitutional Sn impurity, the diffusion barriers increase. Therefore, the substitutional Sn in the diffusion pathway hampers the 3 H diffusion in Zr hydrides.

36 MATERIALS SCIENCE↗
NASA NTRS1975

Primary gamma-rays

Within our Galaxy, cosmic rays can reveal their presence in interstellar space and probably in source regions by their interactions with interstellar matter which lead to gamma-rays with a very characteristic energy spectrum. From the study of the intensity of the high energy gamma radiation as a function of galactic longitude, it is already clear that cosmic rays are almost certainly not uniformly distributed in the Galaxy and are not concentrated in the center of the Galaxy. The galactic cosmic rays appear to be tied to galactic structural features, presumably by the galactic magnetic fields which are in turn held by the matter in the arm segments and the clouds. On the extra-galactic scale, it is now possible to say that cosmic rays are probably not at the density seen near the earth.

Fichtel, C. E.↗
NASA NTRS1999

Effect of a 17 day spaceflight on contractile properties of human soleus muscle fibres

1. Soleus biopsies were obtained from four male astronauts 45 days before and within 2 h after a 17 day spaceflight. 2. For all astronauts, single chemically skinned post-flight fibres expressing only type I myosin heavy chain (MHC) developed less average peak Ca2+ activated force (Po) during fixed-end contractions (0.78 +/- 0. 02 vs. 0.99 +/- 0.03 mN) and shortened at a greater mean velocity during unloaded contractions (Vo) (0.83 +/- 0.02 vs. 0.64 +/- 0.02 fibre lengths s-1) than pre-flight type I fibres. 3. The flight-induced decline in absolute Po was attributed to reductions in fibre diameter and/or Po per fibre cross-sectional area. Fibres from the astronaut who experienced the greatest relative loss of peak force also displayed a reduction in Ca2+ sensitivity. 4. The elevated Vo of the post-flight slow type I fibres could not be explained by alterations in myosin heavy or light chain composition. One alternative possibility is that the elevated Vo resulted from an increased myofilament lattice spacing. This hypothesis was supported by electron micrographic analysis demonstrating a reduction in thin filament density post-flight. 5. Post-flight fibres shortened at 30 % higher velocities than pre-flight fibres at external loads associated with peak power output. This increase in shortening velocity either reduced (2 astronauts) or prevented (2 astronauts) a post-flight loss in fibre absolute peak power (microN (fibre length) s-1). 6. The changes in soleus fibre diameter and function following spaceflight were similar to those observed after 17 days of bed rest. Although in-flight exercise countermeasures probably reduced the effects of microgravity, the results support the idea that ground-based bed rest can serve as a model of human spaceflight. 7. In conclusion, 17 days of spaceflight decreased force and increased shortening velocity of single Ca2+-activated muscle cells expressing type I MHC. The increase in shortening velocity greatly reduced the impact that impaired force production had on absolute peak power.

Clinical Trial↗
DOE OSTI2020

STRIDES: Spectroscopic and photometric characterization of the environment and effects of mass along the line of sight to the gravitational lenses DES J0408–5354 and WGD 2038–4008

In time-delay cosmography, three of the key ingredients are (1) determining the velocity dispersion of the lensing galaxy, (2) identifying galaxies and groups along the line of sight with sufficient proximity and mass to be included in the mass model, and (3) estimating the external convergence κ_ext from less massive structures that are not included in the mass model. We present results on all three of these ingredients for two time-delay lensed quad quasar systems, DES J0408–5354 and WGD 2038–4008 . We use the Gemini, Magellan, and VLT telescopes to obtain spectra to both measure the stellar velocity dispersions of the main lensing galaxies and to identify the line-of-sight galaxies in these systems. Next, we identify 10 groups in DES J0408–5354 and two groups in WGD 2038–4008 using a group-finding algorithm. We then identify the most significant galaxy and galaxy-group perturbers using the ‘flexion shift’ criterion. We determine the probability distribution function of the external convergence κ_ext for both of these systems based on our spectroscopy and on the DES-only multiband wide-field observations. Using weighted galaxy counts, calibrated based on the Millennium Simulation, we find that DES J0408–5354 is located in a significantly underdense environment, leading to a tight (width |$\sim 3{{\ \rm per\ cent}}$|⁠), negative-value κ_ext distribution. On the other hand, WGD 2038–4008 is located in an environment of close to unit density, and its low source redshift results in a much tighter κ_ext of |$\sim 1{{\ \rm per\ cent}}$|⁠, as long as no external shear constraints are imposed.

WGD 2038-4008↗
NASA NTRS1987

A probabilistic model of brittle crack formation

Probability of a brittle crack formation in an elastic solid with fluctuating strength is considered. A set Omega of all possible crack trajectories reflecting the fluctuation of the strength field is introduced. The probability P(X) that crack penetration depth exceeds X is expressed as a functional integral over Omega of a conditional probability of the same event taking place along a particular path. Various techniques are considered to evaluate the integral. Under rather nonrestrictive assumptions, the integral is reduced to solving a diffusion-type equation. A new characteristic of fracture process, 'crack diffusion coefficient', is introduced. An illustrative example is then considered where the integration is reduced to solving an ordinary differential equation. The effect of the crack diffusion coefficient and of the magnitude of strength fluctuations on probability density of crack penetration depth is presented. Practical implications of the proposed model are discussed.

Chudnovsky, A.↗
NASA NTRS1993

Statistics of highly magnified gravitational images in clusters of galaxies. I - Implications for the clusters. II - Implications for the sources

A study is presently conducted of the probability distribution of observing images of background galaxies that are gravitationally lensed in the form of elongated arcs by clusters of galaxies, as a function of arc width, curvature radius, and the elongation's orientation. An upper limit is obtained on the size of the cores which is smaller than cluster core radii obtained for the gas from the X-ray surface brightness. It is argued that lensing clusters tend to be elongated along the line of sight; projection effects will enhance surface density. In the second part of this work, it is shown that the number of arcs expected in lensing clusters is primarily dependent on the angular size of the sources, since lower magnification is required to generate an arc of fixed length for a larger source.

Miralda-Escude, Jordi↗
NASA NTRSDate not supplied

Dose, LET, time and strain dependence of radiation-induced 53BP1 foci in 15 mouse strains ex vivo and associations to in vivo radiation susceptibility

We present a comparative analysis on the repair of radiation-induced DNA damage ex vivo in 15 strains of mice, including 5 inbred reference strains and 10 collaborative-cross strains, of both sexes. Non-immortalized primary skin fibroblasts derived from 76 mice were subjected to both low- and high-LET radiation (0.1, 1 and 4 Gy of X rays; 1.1 and 3 particles/100μm2 of 350 MeV/n 40Ar and 600 MeV/n 56Fe). Automated image quantification of 53BP1 radiation-induced foci (RIF) during the first 4-48 h post-irradiation was performed as a function of dose and LET. Similarly to what we had previously reported for immortalized human cell lines [1], we observed a saturation of RIF number with dose at 4h post-irradiation, with more RIF/Gy for lower LET (X rays and 40Ar) compared to 56Fe. However at later time points (24h and above), the trend was inverted with more RIF/Gy for higher LET. Our data suggest that multiple DSBs cluster into RIF: as the linear density of DSBs increases with LET, so does the probability of having more DSBs per RIF, which makes it more difficult for cells to fully resolve high-LET-induced RIF, explaining the hypersensitivity to high-LET radiation despite a low number of RIF. Taking into account the amount of clustering at a given dose and LET, but also the kinetics of DNA damage repair, we introduced a novel mathematical formalism to evaluate the number of remaining RIF over time. We showed that the newly introduced kinetic metrics can be used as surrogate biomarkers for in vivo radiation toxicity, with potential applications in radiotherapy and human space exploration. In particular, we observed an association between the repairable fraction of RIF measured in vitro and survival levels of immune cells collected from irradiated mice. Moreover, the speed of DNA damage repair correlated with spontaneous cancer incidence data collected from the Mouse Tumor Biology database, suggesting a relationship between the efficiency of DSB repair after irradiation and cancer risk. In addition to the efficacy of repair and persistent RIF levels, even the amount of spontaneous foci without irradiation was shown to be strain dependent, indicating that these phenotypes are at least partially driven by genetics, and supporting their potential as indicators of individual radiation sensitivity. [1] Neumaier, T., et al., PNAS, 2012 (8) 109:443

Radiation, DNA damage, repair kinetics↗
DOE OSTI2020

An extended moments model of quantum efficiency for metals and semiconductors

The complexity of photocathode designs and detector materials, and the need to model their performance for short pulse durations, response to high frequency photons, presence of coatings and/or thinness of the absorptive layer, necessitates modifications to three-step and moments models of photoemission that are used in simulation codes. Here, methods to include input from computationally intensive approaches, such as density functional theory to model optical properties and transfer matrix approaches to treat emission from the surface or transport past coatings, by means of parametric models are demonstrated. First, a technique to accurately represent optical behavior so as to model reflectivity and penetration depth is given. Second, modifications to bulk models arising from the usage of thin film architectures, and a means to rapidly calculate them, is provided. Third, a parameterization to model the impact of wells associated with coatings and surface layers on the transmission probably is given. In all cases, the methods are computationally efficient and designed to allow for including input from numerically intensive approaches that would otherwise be unavailable for simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Spatial analysis of cell patterning to aid genetic and phenotypic understanding of grass stomatal density: A case study in maize

Biological processes involve complex hierarchies where composite traits result from multiple component traits. However, holistically understanding of how sets of component traits interact to underpin genotype-to-phenotype relationships is generally lacking. Stomatal density (SD) is a tractable model system for exploring how high-throughput phenotyping (HTP) data could be exploited by a new spatial analysis approach to better understand a developmentally and functionally important trait. SD is a composite trait, resulting from various components related to cell identity and size, which are themselves governed by a series of spatio-developmental processes. Data from 192 recombinant inbred lines of maize [Zea mays (L.)] were analyzed by a new stomatal patterning phenotype (SPP) to (1) describe the average spatial probability distribution of the nearest neighboring stomata; (2) derive a core set of component traits related to cell size, cell packing, and positional probabilities; (3) build a structural equation model of component traits underlying SD; and (4) identify stomatal patterning quantitative trait loci (QTL). The core set of SPP-derived traits explained 74% of the variation in SD. Analyzing SPP component traits allowed some loci previously identified as generic SD QTL to be recognized as specific to lateral versus longitudinal elements of stomatal patterning. Therefore, this study highlights how novel insights can be gained by decomposing a composite trait (e.g., SD) into a set of component traits that were present in HTP data but not previously exploited.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2007

Code Assesses Risks Posed by Meteoroids and Orbital Debris

BUMPER II version 1.92e is a computer code for assessing the risk of damage from impacts of micrometeoroids and orbital debris on the International Space Station (ISS), including those parts of the ISS covered by shielding that affords partial protection against such impacts. (Other versions of BUMPER II have been written for other spacecraft.) Bumper II quantifies the probability of penetration of shielding and the damage to spacecraft equipment as functions of the size, shape, and orientation of the spacecraft; the parameters of its orbit; failure criteria that quantify impact damage at the threshold of failure for each spacecraft surface; and the impact-damage resistance of each spacecraft surface as defined by "ballistic limit equations" that return the size of a failure-causing particle as a function of target parameters (including materials, configurations, thicknesses, and gap distances) and impact conditions (impact velocity and the density and shape of the impactor). BUMPER II version 1.92e contains several dozen ballistic limit equations that are based on results from thousands of hypervelocity impact tests conducted by NASA on ISS shielding and other hardware, and on results from numerical simulations of impacts.

Christiansen, Eric L.↗