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 847 records · Page 47

Tearing at the dayside magnetopause

Initial observations of the dayside magnetopause thickness using magnetometers on board ISEE 1 (International Sun-Earth Explorer) and ISEE 2 indicate that this boundary may only be a few ion gyro-radii thick during periods of observed southward interplanetary magnetic field (IMF). Such a thin current layer can destabilize the collisionless tearing mode, a probable first step in the reconnection sequence. The linear growth rate is calculated assuming typical magnetopause parameters and it is found that it is a sensitive function of IMF orientation, magnetopause thickness, and electron number density. It is also shown that the finite size of the dayside magnetopause necessitates the existence of a two-dimensional tearing mode wave vector spectrum, a consequence of the inability of a thermal electron to maintain Landau resonance with the wave for a growth period. Implications regarding reconnection are discussed.

Quest, K. B.↗

A complete southern sky redshift survey

We discuss a magnitude-limited redshift survey covering a region around the south Galactic pole (SSRS2). The survey includes 3592 galaxies and covers a region of 1.13 sr to a limiting m(sub B(0)) = 15.5, directly comparable to CfA2. Together, CfA2 and SSR2 cover more than a third of the sky. The northern and southern surveys are remarkably similar. Both contain voids with diameters as large as 5000 km/s. The southern survey contains the Southern Wall, similar to the northern Great Wall. The distributions of velocity dispersions for systems extracted from the SSRS2 and CfA2 surveys are also remarkably similar. For the SSRS2, an inhomogeneity-independent technique yields Schechter luminosity function parameters M(sup *) = -19.50 and alpha = -1.20 (H(sub O) = 100 h km/s/Mpc, h = 1 unless otherwise specified). The steep faint-end slope is probably attributable to bluer galaxies. By examining the normalized density fluctuations in the SSRS2 and CfA2 surveys, we suggest that the combined sample is not yet large enough to be 'fair.' There are large fluctuations in shells at 10,000 km/s.

Da Costa, L. Nicolaci↗

Effective spin-mixing conductance of topological-insulator/ferromagnet and heavy-metal/ferromagnet spin-orbit-coupled interfaces: A first-principles Floquet-nonequilibrium Green function approach

The spin-mixing conductance (SMC) is a key quantity determining efficiency of spin transport across interfaces. Thus, knowledge of its precise value is required for accurate measurement of parameters quantifying numerous effects in spintronics, such as spin-orbit torque, spin Hall magnetoresistance, spin Hall effect, and spin pumping. However, the standard expression for SMC, provided by the scattering theory in terms of the reflection probability amplitudes, is inapplicable when strong spin-orbit coupling (SOC) is present directly at the interface. This is the precisely the case of topological-insulator/ferromagnet and heavy-metal/ferromagnet interfaces of great contemporary interest. We introduce an approach where first-principles Hamiltonian of these interfaces, obtained from noncollinear density functional theory (ncDFT) calculations, is combined with charge-conserving Floquet-nonequilibrium-Green-function formalism to compute directly the pumped spin current $I^{S_z}_L$ into semi-infinite left lead of two-terminal heterostructures Cu/X/Co/Cu or Y/Co/Cu—where X = Bi 2 Se 3 and Y = Pt or W—due to microwave-driven steadily precessing magnetization of the Co layer. This allows us to extract an effective SMC as a prefactor in $I^{S_z}_L$ versus precession cone angle θ dependence, as long as it remains the same, $I^{S_z}$ ∝ sin 2 θ , as in the case where SOC is absent. By comparing calculations where SOC in switched off versus switched on in ncDFT calculations, we find that SOC consistently reduces the pumped spin current and therefore the effective SMC.

36 MATERIALS SCIENCE↗

X-ray clusters from a high-resolution hydrodynamic PPM simulation of the cold dark matter universe

A new three-dimensional hydrodynamic code based on the piecewise parabolic method (PPM) is utilized to compute the distribution of hot gas in the standard Cosmic Background Explorer (COBE)-normalized cold dark matter (CDM) universe. Utilizing periodic boundary conditions, a box with size 85 h(exp-1) Mpc, having cell size 0.31 h(exp-1) Mpc, is followed in a simulation with 270(exp 3)=10(exp 7.3) cells. Adopting standard parameters determined from COBE and light-element nucleosynthesis, Sigma(sub 8)=1.05, Omega(sub b)=0.06, we find the X-ray-emitting clusters, compute the luminosity function at several wavelengths, the temperature distribution, and estimated sizes, as well as the evolution of these quantities with redshift. The results, which are compared with those obtained in the preceding paper (Kang et al. 1994a), may be used in conjuction with ROSAT and other observational data sets. Overall, the results of the two computations are qualitatively very similar with regard to the trends of cluster properties, i.e., how the number density, radius, and temeprature depend on luminosity and redshift. The total luminosity from clusters is approximately a factor of 2 higher using the PPM code (as compared to the 'total variation diminishing' (TVD) code used in the previous paper) with the number of bright clusters higher by a similar factor. The primary conclusions of the prior paper, with regard to the power spectrum of the primeval density perturbations, are strengthened: the standard CDM model, normalized to the COBE microwave detection, predicts too many bright X-ray emitting clusters, by a factor probably in excess of 5. The comparison between observations and theoretical predictions for the evolution of cluster properties, luminosity functions, and size and temperature distributions should provide an important discriminator among competing scenarios for the development of structure in the universe.

Bryan, Greg L.↗

Characterizing a supernova’s standing accretion shock instability with neutrinos and gravitational waves

Here, we perform a novel multi-messenger analysis for the identification and parameter estimation of the Standing Accretion Shock Instability (SASI) in a core collapse supernova with neutrino and gravitational wave (GW) signals. In the neutrino channel, this method performs a likelihood ratio test for the presence of SASI in the frequency domain. For gravitational wave signals we process an event with a modified constrained likelihood method. Using simulated supernova signals, the properties of the Hyper-Kamiokande neutrino detector, and O3 LIGO Interferometric data, we produce the two-dimensional probability density distribution (PDF) of the SASI activity indicator and calculate the probability of detection P D as well as the false identification probability P FI . We discuss the probability to establish the presence of the SASI as a function of the source distance in each observational channel, as well as jointly. Compared to a single-messenger approach, the joint analysis results in P D (at P FI = 0.1) of SASI activities that is larger by up to ≈ 40% for a distance to the supernova of 5 kpc. We also discuss how accurately the frequency and duration of the SASI activity can be estimated in each channel separately. Our methodology is suitable for implementation in a realistic data analysis and a multi-messenger setting.

79 ASTRONOMY AND ASTROPHYSICS↗

Towards Probablistic Assessment of Hypobaric Decompression Sickness Treatment

INTRODUCTION: Pressure, oxygen (O2), and time are the pillars to effective treatment of decompression sickness (DCS). The NASA DCS Treatment Model links a decrease in computed bubble volume to the resolution of a symptom. The decrease in volume is realized in two stages: a) during the Boyle's Law compression and b) during subsequent dissolution of the gas phase by the O2 window. METHODS: The cumulative distribution of 154 symptoms that resolved during repressurization was described with a log-logistic density function of pressure difference (deltaP as psid) associated with symptom resolution and two other explanatory variables. The 154 symptoms originated from 119 cases of DCS during 969 exposures in 47 different altitude tests. RESULTS: The probability of symptom resolution [P(symptom resolution)] = 1 / (1+exp(- (ln(deltaP) - 1.682 + 1.089×AMB - 0.00395×SYMPTOM TIME) / 0.633)), where AMB is 1 when the subject ambulated as part of the altitude exposure or else 0 and SYMPTOM TIME is the elapsed time in min from start of the altitude exposure to recognition of a DCS symptom. The P(symptom resolution) was estimated from computed deltaP from the Tissue Bubble Dynamics Model based on the "effective" Boyle's Law change: P2 - P1 (deltaP, psid) = P1×V1/V2 - P1, where V1 is the computed volume of a spherical bubble in a unit volume of tissue at low pressure P1 and V2 is computed volume after a change to a higher pressure P2. V2 continues to decrease through time at P2, at a faster rate if 100% ground level O2 was breathed. The computed deltaP is the effective treatment pressure at any point in time as if the entire deltaP was just from Boyle's Law compression. DISCUSSION: Given the low probability of DCS during extravehicular activity and the prompt treatment of a symptom with options through the model it is likely that the symptom and gas phase will resolve with minimum resources and minimal impact on astronaut health, safety, and productivity.

Conkin, J.↗

Investigation of Transient Gas Phase Column Density Due to Droplet Evaporation

Analytical expressions are developed for estimating column density near a rapidly evaporating droplet along general paths. The influence for instantaneous evaporation is created first as a limiting case, where the peak value occurs at the time it takes a wave of vapor to reach the closest point along the optical path traveling at its most probable thermal speed. Next the case for finite-period evaporation is evaluated for constant conditions. Compared to the instantaneous case, peak column density occurs shortly after droplet extinction at a lower intensity. A new mathematical function is discovered that solves the integrals associated with this case. Finally, ways to account for droplet motion and changes in evaporation rate with size and temperature are discussed.

transient↗

Investigation of Transient Gas Phase Column Density Due to Droplet Evaporation

Analytical expressions are developed for estimating column density near a rapidly evaporating droplet along general paths. The influence for instantaneous evaporation is created first as a limiting case, where the peak value occurs at the time it takes a wave of vapor to reach the closest point along the optical path traveling at its most probable thermal speed. Next the case for finite-period evaporation is evaluated for constant conditions. Compared to the instantaneous case, peak column density occurs shortly after droplet extinction at a lower intensity. A new mathematical function is discovered that solves the integrals associated with this case. Finally, ways to account for droplet motion and changes in evaporation rate with size and temperature are discussed.

Column Density↗

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Photometric redshifts for the S-PLUS Survey: Is machine learning up to the task?

The Southern Photometric Local Universe Survey (S-PLUS) is a novel project that aims to map the Southern Hemisphere using a twelve filter system, comprising five broad-band SDSS-like filters and seven narrow-band filters optimized for important stellar features in the local universe. In this paper we use the photometry and morphological information from the first S-PLUS data release (S-PLUS DR1) cross-matched to unWISE data and spectroscopic redshifts from Sloan Digital Sky Survey DR15. We explore three different machine learning methods (Gaussian Processes with GPz and two Deep Learning models made with TensorFlow) and compare them with the currently used template-fitting method in the S-PLUS DR1 to address whether machine learning methods can take advantage of the twelve filter system for photometric redshift prediction. Using tests for accuracy for both single-point estimates such as the calculation of the scatter, bias, and outlier fraction, and probability distribution functions (PDFs) such as the Probability Integral Transform (PIT), the Continuous Ranked Probability Score (CRPS) and the Odds distribution, we conclude that a deep-learning method using a combination of a Bayesian Neural Network and a Mixture Density Network offers the most accurate photometric redshifts for the current test sample. In conclusion, it achieves single-point photometric redshifts with scatter (σ NMAD ) of 0.023, normalized bias of -0.001, and outlier fraction of 0.64% for galaxies with r_auto magnitudes between 16 and 21.

79 ASTRONOMY AND ASTROPHYSICS↗

Some effects of vibration and rotation on the drift of gyroscopic instruments

It is the purpose of this paper to describe and indicate possible analytical approaches to some of the most significant mechanisms by which vibration induced drifts may occur, and to indicate, where possible, some means of alleviating them. In all cases we will not consider gimbal bearing friction, mass unbalance or other sources of steady extraneous torques on the gyro rotor. The mechanisms to be discussed are: Non-Newtonian viscosity in single degree of freedom (SDF) integrating-rate gyros, nonisoelasticity in SDF gyros, noncommutative effects for SDF stabilized platform, gimbal inertia effects, 2 DF gyros, and spin axis torque coupling in 2 DF gyros. High-frequency elastic resonances are not discussed explicitly in detail, because they are a more familiar problem and one which depends strongly on the details of each particular gyro design.

Viscous Torque↗

Investigation of Transient Gas Phase Column Density Due to Droplet Evaporation

A NASA robotic refueling mission experiment on board the International Space Station (ISS) was designed to repeatedly transfer simulated cryogenic propellant between two dewars. After each metered transfer, the cryogen was to be vented overboard via sublimation or evaporation. ISS payload providers must conduct analyses to demonstrate that any planned gaseous venting generate no more than a certain level of material that may interfere with optical measurements from other payloads that may be located nearby. This requirement is expressed in terms of a maximum column number density (CND). For the experiment under consideration, fluid droplets likely would accompany any such vapor. Earlier related studies led to development of analytical expressions for column density along general paths for a number of practical cases under steady conditions, including droplet evaporation. These expressions are not appropriate for rapid evaporation of small droplets however, since the droplets cannot sustain mass flow rates consistent with steady conditions. In this work, analytical expressions are developed for estimating column density near a rapidly evaporating droplet along general paths. The influence for instantaneous evaporation is created first as a limiting case, where it is found the peak value occurs at the time it takes a wave of vapor to reach the closest point along the optical path traveling at its most probable thermal speed. Next, the case for finite-period evaporation is evaluated for constant conditions. Compared to the instantaneous case, peak column density occurs shortly after droplet extinction, but at a lower intensity. A new mathematical function is discovered that solves the integrals associated with this case. Finally, ways to account for droplet motion and changes in evaporation rate with size and temperature are discussed in this framework.

Woronowicz, M. S.↗

Oxygen vacancy-induced ferroelectric effect in (111) strontium titanate single crystals controls photoelectrochemical water oxidation

Ferroelectric materials, such as tetragonal BaTiO 3 , have a permanent electric polarization that can be controlled with an external electric field, however, a ferroelectric polarization in cubic SrTiO 3 is forbidden by the higher symmetry of the lattice. Here we demonstrate that hydrogen annealed SrTiO 3−x single crystals can be polarized electrically, and that the polarization controls the activity for photoelectrochemical water oxidation, a pathway to solar hydrogen fuel. Specifically, it is observed that the anodic water oxidation photocurrent increases from 0.99 to 2.22 mA cm −2 at 1.23 V RHE (60 mW cm −2 , UV illumination) or decreases to 0.50 mA cm −2 after electric polarization of hydrogen-annealed (111) SrTiO 3−x single crystals in forward or reverse direction. The polarization also modifies the surface photovoltage signal of the material and its flat band potential, based on Mott–Schottky measurements. These observations are attributed to the formation of an electric dipole at the (111) SrTiO 3−x surface, which alters the potential drop across the depletion layer at the solid–liquid junction, and with it the electron transfer barrier. Density functional theory calculations confirm that an electric dipole can result from the movement of oxygen vacancies between the surface or sub-surface layers of SrTiO 3−x . The filling of these surface oxygen vacancies is the probable cause for the observed disappearance of the electric polarization after 24 h storage in air and 48 h in argon. Overall, this work establishes a new surface-based ferroelectric effect in SrTiO 3−x and its use for solar energy conversion during photoelectrochemical water oxidation. Because oxygen vacancy defects are common, similar electric polarization effects are to be expected in other metal oxides.

Assavachin, Samutr [University of California, Davi↗

Beyond the proton drip line: Bayesian analysis of proton-emitting nuclei

The limits of the nuclear landscape are determined by nuclear binding energies. Beyond the proton drip lines, where the separation energy becomes negative, there is not enough binding energy to prevent protons from escaping the nucleus. Predicting properties of unstable nuclear states in the vast territory of proton emitters poses an appreciable challenge for nuclear theory as it often involves far extrapolations. In addition, significant discrepancies between nuclear models in the proton-rich territory call for quantified predictions. With the help of Bayesian methodology, we mix a family of nuclear mass models corrected with statistical emulators trained on the experimental mass measurements, in the proton-rich region of the nuclear chart. Separation energies were computed within nuclear density functional theory using several Skyrme and Gogny energy density functionals. We also considered mass predictions based on two models used in astrophysical studies. Quantified predictions were obtained for each model using Bayesian Gaussian processes trained on separation-energy residuals and combined via Bayesian model averaging. Here, we obtained a good agreement between averaged predictions of statistically corrected models and experiment. In particular, we quantified model results for one- and two-proton separation energies and derived probabilities of proton emission. This information enabled us to produce a quantified landscape of proton-rich nuclei. The most promising candidates for two-proton decay studies have been identified. Finally, the methodology used in this work has broad applications to model-based extrapolations of various nuclear observables. It also provides a reliable uncertainty quantification of theoretical predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

VLA and low-frequency VLBI observations of the radio source 0503 + 467 - Austere constraints on interstellar scattering in two media

The radio source 0503 + 467 lies near the Galactic plane (l = 161.0 deg, b = 3.7 deg) and at the edge of the supernova remnant (SNR) HB 9. The VLA observations show that it has a spectrum typical of a compact extragalactic radio source. The resultant small angular size of the source makes it an excellent probe of turbulence in two media: the diffuse, or 'type A', component of interstellar turbulence and a hypothesized region of hydromagnetic turbulence upstream of the supernova remnant. An eight-station VLBI experiment at 326 MHz indicates that the source is less than about 20 milliarcseconds (mas) in angular diameter. A value of 16 mas is most appropriate as an upper limit to the interstellar scattering contribution to the measured angular size. The implications of this upper limit are twofold. First, the galactocentric radial scale to the type-A turbulence is probably less than or equal to about 6 kpc. Second, no evidence is seen for shock-associated turbulence upstream of HB 9. The measurements make it possible to constrain a parameter which is a function of the rms density fluctuation in the upstream region, the outer scale to the density turbulence, and the thickness of SNR foreshock region.

Spangler, Steven R.↗

Data for 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.

AI/ML↗

One-bond 13 C– 13 C spin-coupling constants in saccharides: a comparison of experimental and calculated values by density functional theory using solid-state 13 C NMR and X-ray crystallography

Methyl aldohexopyranosides were 13 C-labeled at contiguous carbons, crystallized, and studied by single-crystal X-ray crystallography and solid-state 13 C nuclear magnetic resonance (NMR) spectroscopy to examine the degree to which density functional theory (DFT) can calculate one-bond 13C–13C spin-coupling constants ( 1 J CC ) in saccharides with sufficient accuracy to permit their use in MA'AT analysis, a newly-reported hybrid DFT/NMR method that provides probability distributions of molecular torsion angles in solution (Zhang et al., J. Phys. Chem. B, 2017, 121, 3042–3058; Meredith et al., J. Chem. Inf. Model., 2022, 62, 3135–3141). Experimental 1 J CC values in crystalline samples of the doubly 13 C-labeled compounds were measured by solid-state 13 C NMR and compared to those calculated from five different DFT models: (1) 1 J CC values calculated from single structures identical to those observed in crystalline samples by X-ray crystallography (all atom refinement); (2) 1 J CC values calculated from the single structures in (1) but after Hirshfeld atom refinement (HAR); (3) 1 J CC values calculated from the single structures in (1) after DFT-optimization of hydrogen atoms only; and (4 and 5) 1 J CC values calculated in rotamers of torsion angle θ 2 (C1–C2–O2–O2H) or ω(C4–C5–C6–O6) from which either specific or generalized parameterized equations were obtained and used to calculate 1 J CC values in the specific θ 2 or ω rotamers observed in crystalline samples. Good qualitative agreement was observed between calculated 1 J CC values and those measured by solid-state 13 C NMR regardless of the DFT model, but in no cases were calculated 1 J CC values quantitative, differing (over-estimated) on average by 4–5% from experimental values. These findings, and those reported recently from solution NMR studies (Tetrault et al., J. Phys. Chem. B 2022, 126, 9506–9515), indicate that improvements in DFT calculations are needed before calculated 1 J CC values can be used directly as reliable constraints in MA'AT analyses of saccharides in solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetic simulation of a 50 mTorr capacitively coupled argon discharge over a range of frequencies and comparison to experiments

The effect of driving frequency in the range of 13.56–73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct “temperatures” for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2–10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25% of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model, and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗