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 919 records · Page 51

Extreme Mean and Its Applications

Extreme value statistics obtained from normally distributed data are considered. An extreme mean is defined as the mean of p-th probability truncated normal distribution. An unbiased estimate of this extreme mean and its large sample distribution are derived. The distribution of this estimate even for very large samples is found to be nonnormal. Further, as the sample size increases, the variance of the unbiased estimate converges to the Cramer-Rao lower bound. The computer program used to obtain the density and distribution functions of the standardized unbiased estimate, and the confidence intervals of the extreme mean for any data are included for ready application. An example is included to demonstrate the usefulness of extreme mean application.

Swaroop, R.↗

On the alleged collisional origin of the Kirkwood Gaps

This paper examines two proposed mechanisms whereby asteroidal collisions and close approaches may have given rise to the Kirkwood Gaps. The first hypothesis is that asteroids in near-resonant orbits have markedly increased collision probabilities and so are preferentially destroyed, or suffer decay in population density, within the resonance zones. A simple order-of-magnitude analysis shows that this hypothesis is untenable since it leads to conclusions which are either unrealistic or not in accord with present understanding of asteroidal physics. The second hypothesis is the Brouwer-Jefferys theory that collisions would smooth an asteroidal distribution function, as a function of Jacobi constant, thus forming resonance gaps. This hypothesis is examined by direct numerical integration of 50 asteroid orbits near the 2:1 resonance, with collisions simulated by random variables. No tendency to form a gap was observed.

Heppenheimer, T. A.↗

Negative Ion Chemistry in the Coma of Comet 1P/Halley

Negative ions (anions) were identified in the coma of comet 1P/Halley from in-situ measurements performed by the Giotto spacecraft in 1986. These anions were detected with masses in the range 7-110 amu, but with insufficient mass resolution to permit unambiguous identification. We present details of a new chemical-hydrodynamic model for the coma of comet Halley that includes - for the first time - atomic and molecular anions, in addition to a comprehensive hydrocarbon chemistry. Anion number densities arc calculated as a function of radius in the coma, and compared with the Giotto results. Important anion production mechanisms arc found to include radiative electron attachment, polar photodissociation, dissociative electron attachment, and proton transfer. The polyyne anions C4H(-) and C6H(-) arc found to be likely candidates to explain the Giotto anion mass spectrum in the range 49-73 amu. Thc CN(-) anion probably makes a significant contribution to the mass spectrum at 26 amu. Larger carbon-chain anions such as C8H(1) can explain the peak near 100 amu provided there is a source of large carbon-chain-bearing molecules from the cometary nucleus.

Cordiner, M. A.↗

Self-similarity of k -nearest neighbour distributions in scale-free simulations

ABSTRACT We use the k-nearest neighbour probability distribution function (kNN-PDF; Banerjee & Abel 2021a) to assess convergence in a scale-free N-body simulation. Compared to our previous two-point analysis, the kNN-PDF allows us to quantify our results in the language of haloes and numbers of particles, while also incorporating non-Gaussian information. We find good convergence for 32 particles and greater at densities typical of haloes, while 16 particles and fewer appear unconverged. Halving the softening length extends convergence to higher densities, but not to fewer particles. Our analysis is less sensitive to voids, but we analyse a limited range of underdensities and find evidence for convergence at 16 particles and greater even in sparse voids.

79 ASTRONOMY AND ASTROPHYSICS↗

A microlensing model for QSO 2237 + 0305

Models of microlensing are presented for the four images of QSO 2237 + 0305, with the surface mass density and shear for each image based on the macrolensing model developed by Schneider et al. (1988). The computations were done with a new code that makes use of the hierarchical tree method. The high efficiency of this approach makes it possible to include up to 60,000 lensing stars with different masses. Examples are given of amplification maps, light curves, correlation functions, and amplification probabilities. All these are very strongly affected by the shear due to macrolensing. The observed changes in the luminosity of the four images of QSO 2237 + 0305 are most likely induced by microlensing.

Wambsganss, J.↗

A conservative Galerkin solver for the quasilinear diffusion model in magnetized plasmas

We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron impact excitation of highly charged sodium-like ions

Optical transition probabilities and electron collision strengths for Ca X, Fe XVI, Zn XX, Kr XXVI and Mo XXXII are calculated for transitions between n equal to 3 and n equal to 4 levels. The calculations neglect relativistic effects on the radial functions. A semi-empirical approach provides wave functions of the excited states; a distorted wave function without exchange is employed to obtain the excitation cross sections. The density dependence of the relative intensities of certain emission lines in the sodium isoelectronic sequence is also discussed.

Blaha, M.↗

An unbiased survey of field star X-ray emission

To determine the X-ray luminosity function of normal stars, 1700 stars brighter than 10th magnitude were surveyed by the imaging proportional counter aboard the Einstein Observatory. Seventy star positions were found to contain excess X-ray counting rates. The number of stars per square degree in a number of magnitude intervals was calculated as a function of spectral type and luminosity class, and the total number of stars for each spectral type brighter than magnitude 9.5 derived in this manner was compared with the 1700-star sample. The agreement is good, as is that between the surface density of soft X-ray sources and the number of stellar emitters predicted from the field star survey. It is concluded that stars probably do not contribute significantly to the diffuse soft X-ray background. The findings are consistent with the notion that stellar age and/or rotation velocity are important determinants of stellar X-ray emission level.

Helfand, D. J.↗

Observations of the earth's bow shock under high Mach number/high plasma beta solar wind conditions

Using 221 observations of the earth's bow shock by ISEE-1, the magnetic field characteristics of the shock are investigated as a function of Mach number and plasma beta. It is found that immediately behind the overshoot of the shock large amplitude waves develop preferentially with increasing Alfven Mach number and/or plasma beta. Their probability of occurrence is very high in the parallel regime, and declines to near zero probability when theta(BN) reaches 75 deg. It is estimated that approximately 10 percent of the total solar wind energy density is dissipated by these waves.

Winterhalter, Daniel↗

Anomalies and Local Structure of Liquid Water from Boiling to the Supercooled Regime as Predicted by the Many-Body MB-pol Model

For the past 50 years, researchers have sought molecular models that can accurately reproduce water’s microscopic structure and thermophysical properties across broad ranges of its complex phase diagram. Herein, molecular dynamics simulations with the many-body MB-pol model are performed to monitor the thermodynamic response functions and local structure of liquid water from the boiling point down to deeply supercooled temperatures at ambient pressure. The isothermal compressibility and isobaric heat capacity show maxima near 223 K, in excellent agreement with recent experiments, and the liquid density exhibits a minimum at ~208 K. A local tetrahedral arrangement, where each water molecule accepts and donates two hydrogen bonds, is found to be the most probable hydrogen-bonding topology at all temperatures. Further, this work suggests that MB-pol may provide predictive capability for studies of liquid water’s physical properties across broad ranges of thermodynamic states, including the so-called water’s “no man’s land” which is difficult to probe experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Statistical Standard Siren Measurement of the Hubble Constant from the LIGO/Virgo Gravitational Wave Compact Object Merger GW190814 and Dark Energy Survey Galaxies

We present a measurement of the Hubble constant H 0 using the gravitational wave (GW) event GW190814, which resulted from the coalescence of a 23 M ⊙ black hole with a 2.6 M ⊙ compact object, as a standard siren. No compelling electromagnetic counterpart has been identified for this event; thus our analysis accounts for thousands of potential host galaxies within a statistical framework. The redshift information is obtained from the photometric redshift (photo- z ) catalog from the Dark Energy Survey. The luminosity distance is provided by the LIGO/Virgo gravitational wave sky map. Since this GW event has the second-smallest localization volume after GW170817, GW190814 is likely to provide the best constraint on cosmology from a single standard siren without identifying an electromagnetic counterpart. Furthermore, our analysis uses photo- z probability distribution functions and corrects for photo- z biases. We also reanalyze the binary black hole GW170814 within this updated framework. We explore how our findings impact the H 0 constraints from GW170817, the only GW merger associated with a unique host galaxy. From a combination of GW190814, GW170814, and GW170817, our analysis yields H 0 =72.0$_{+12}^{–8.2}$ km s -1 Mpc -1 (68% highest-density interval, HDI) for a prior in H 0 uniform between [20 and 140] km s -1 Mpc -1 . The addition of GW190814 and GW170814 to GW170817 improves the 68% HDI from GW170817 alone by ~18%, showing how well-localized mergers without counterparts can provide a significant contribution to standard siren measurements, provided that a complete galaxy catalog is available at the location of the event.

79 ASTRONOMY AND ASTROPHYSICS↗

An altitude-dependent spacecraft charging model

A model for the altitude dependence of the hot plasma parameters responsible for the electrostatic charging of spacecraft was developed. Based upon plasma orbit theory, the directed velocity is a function of the ambient magnetic field flux density. A consequence of this approach is that while the thermal velocity distributions (assumed to be Maxwellian) of the plasma particles are independent of the magnetic field strength (and hence altitude), the particle densities increase with magnetic field strength. Thus, according to this model, while the equilibrium voltage is independent of altitude, the charging current density increases with decreasing altitude. However, the probability of such spacecraft charging decreases with decreasing altitude.

Haffner, J. W.↗

Variation of Lunar Sodium During Passage of the Moon through the Earth's Magnetotail

We measured sodium emission above the lunar equator over a range of lunar altitudes from 100 to 4000 km. The measurements were repeated approximately every 24 hours from June 7 to 16, 1998, covering the period during which the Moon passed through the Earth's magnetotail. Sodium temperatures derived from the altitude dependence of emission intensity ranged from 1200 to 2900 K. This result supports the view that photodesorption is a primary source of sodium in the exosphere since the most probable temperature of sodium form this source is in this range. Passage of the Moon through the Earth's magnetotail (where solar wind is essentially absent) affected the sodium density, such that it was higher before the Moon entered the Earth's magnetotail than after the Moon left it. This suggests that the solar wind plays a role in production of lunar sodium. We propose that its function is to mobilize sodium and bring it to the surface, where photodesorption can eject it into the exosphere. A two-step process such as this could help to explain the latitude dependence of sodium density, which varies as the second or higher power of cosine latitude.

Potter, Andrew E.↗

Extending Newtonian Dynamics to Include Stochastic Processes

A paper presents further results of continuing research reported in several previous NASA Tech Briefs articles, the two most recent being Stochastic Representations of Chaos Using Terminal Attractors (NPO-41519), [Vol. 30, No. 5 (May 2006), page 57] and Physical Principle for Generation of Randomness (NPO-43822) [Vol. 33, No. 5 (May 2009), page 56]. This research focuses upon a mathematical formalism for describing post-instability motions of a dynamical system characterized by exponential divergences of trajectories leading to chaos (including turbulence as a form of chaos). The formalism involves fictitious control forces that couple the equations of motion of the system with a Liouville equation that describes the evolution of the probability density of errors in initial conditions. These stabilizing forces create a powerful terminal attractor in probability space that corresponds to occurrence of a target trajectory with probability one. The effect in configuration space (ordinary three-dimensional space as commonly perceived) is to suppress exponential divergences of neighboring trajectories without affecting the target trajectory. As a result, the post-instability motion is represented by a set of functions describing the evolution of such statistical quantities as expectations and higher moments, and this representation is stable.

Zak, Michail↗

New method for estimating low-earth-orbit collision probabilities

An unconventional but general method is described for estimating the probability of collision between an earth-orbiting spacecraft and orbital debris. This method uses a Monte Caralo simulation of the orbital motion of the target spacecraft and each discrete debris object to generate an empirical set of distances, each distance representing the separation between the spacecraft and the nearest debris object at random times. Using concepts from the asymptotic theory of extreme order statistics, an analytical density function is fitted to this set of minimum distances. From this function, it is possible to generate realistic collision estimates for the spacecraft.

Vedder, John D.↗

ASCA measurements of the grain-scattered X-ray halos of eclipsing massive X-ray binaries: Vela X-1 and Centaurus X-3

We have measured the decaying dust-scattered X-ray halo of Cen X-3 during its binary eclipse with the ASCA solid-state imaging spectrometer (SIS). The surface brightness profile (SBP) of the image in the low-energy band (0.5-3 keV) lies substantially above the point-spread function (PSF) of the X-ray telescope, while the SBP in the high-energy band (5-10 keV) exhibits no significant deviation. By contrast, the SBPs of Vela X-1 during its eclipse are consistent with the PSF in both the low- and high-energy bands -- strong evidence that a dust halo is indeed present in Cen X-3. Accordingly, we modeled the SBP of Cen X-3 taken from six consecutive time segments under the principal assumptions that the dust is distributed uniformly along a segment of the line of sight, the grains have a power-law size distribution, and the low-energy source flux was the same function of orbital phase before as during our observation. The best-fit set of parameters included a grain density value of 1.3 g/cu cm, substanially less than the density of 'astronomical silicate.' This result supports the idea that interstellar grains are 'fluffy' aggregates of smaller solid particles. We attribute the failure to detect a halo of Vela X-1 during its eclipse phase to extended strong circumsource absorption that probably occurred before the eclipse and allowed the halo to decay away before the observation began.

Woo, Jonathan W.↗

Ground-state properties of the Hubbard model in one and two dimensions from the Gutzwiller conjugate gradient minimization theory

We introduce Gutzwiller conjugate gradient minimization (GCGM) theory, an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. GCGM uses the Gutzwiller wave function but does not use the commonly adopted Gutzwiller approximation (GA), which is a major source of inaccuracy. Instead, the theory uses an approximation that is based on the occupation probability of the on-site configurations, rather than approximations that decouple the site-site correlations as used in the GA. We test the theory in the one-dimensional and two-dimensional Hubbard models at various electron densities and find that GCGM reproduces energies and double occupancies in reasonable agreement with benchmark data at a very small computational cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Limb-brightening observations from the OSO-7 satellite. II - Comparison of Abel-inverted intensities of Fe XIV and Fe XIII EUV emission lines with predictions

Intensities of Fe XIV and Fe XIII EUV emission lines obtained at coronal locations beyond the limb by the Goddard spectroheliograph on the OSO 7 satellite have been corrected for the wavelength dependence of the instrument's sensitivity and have been Abel-inverted to provide a valid comparison with theoretical predictions for each ion. Details of the Abel-inversion procedure are given, including explicit formulas for application of Bracewell's (1956) method. The intensity ratios of pairs of lines originating from a common level are compared with expected theoretical transition probability ratios over a range of heliocentric distance; deviations in some cases yield information about adjacent unclassified lines. Comparison of the observations with predictions for Fe XIV and Fe XIII shows generally good agreement, with a few interesting discrepancies that may imply a corresponding need for more accurate collisional excitation cross sections. The same comparison yields the variation of electron density with heliocentric radius for each ion separately; the two density functions are found to agree within a factor of three.

Kastner, S. O.↗