Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DISTRIBUTION FUNCTIONS”

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 505 records · Page 28

Angular momentum eigenstates of the isotropic 3-D harmonic oscillator: Phase-space distributions and coalescence probabilities

The isotropic 3-dimensional harmonic oscillator potential can serve as an approximate description of many systems in atomic, solid state, nuclear, and particle physics. In particular, the question of 2 particles binding (or coalescing) into angular momentum eigenstates in such a potential has interesting applications. Here we compute the probabilities for coalescence of two distinguishable, non-relativistic particles into such a bound state, where the initial particles are represented by generic wave packets of given average positions and momenta. We use a phase-space formulation and hence need the Wigner distribution functions of angular momentum eigenstates in isotropic 3-dimensional harmonic oscillators. These distribution functions have been discussed in the literature before but we utilize an alternative approach to obtain these functions. Along the way, we derive a general formula that expands angular momentum eigenstates in terms of products of 1-dimensional harmonic oscillator eigenstates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations

In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.

36 MATERIALS SCIENCE↗

Photon-photon absorption and the uniqueness of the spectra of active galactic nuclei

The effects of the feedback of e(+)-e(-) pair reinjection in a plasma due to photon-photon absorption of its own radiation was examined. Under the assumption of continuous electron injection with a power law spectrum E to the minus gamma power and Compton losses only, it is shown that for gamma 2 the steady state electron distribution function has a unique form independent of the primary injection spectrum. This electron distribution function can, by synchrotron emission, reproduce the general characteristics of the observed radio to optical active galactic nuclei spectra. Inverse Compton scattering of the synchrotron photons by the same electron distribution can account for their X-ray spectra, and also implies gamma ray emission from these objects. This result is invoked to account for the similarity of these spectra, and it is consistent with observations of the diffuse gamma ray background.

Kazanas, D.↗

Photon-photon absorption and the uniqueness of the spectra of active galactic nuclei

The effects of the feedback of e(+)-e(-) pair reinjection in a plasma due to photon-photon absorption of its own radiation was examined. Under the assumption of continuous electron injection with a power law spectrum E to the minus gamma power and Compton losses only, it is shown that for gamma 2 the steady state electron distribution function has a unique form independent of the primary injection spectrum. This electron distribution function can, by synchrotron emission, reproduce the general characteristics of the observed radio to optical active galactic nuclei spectra. Inverse Compton scattering of the synchrotron photons by the same electron distribution can account for their X-ray spectra, and also implies gamma ray emission from these objects. This result is invoked to account for the similarity of these spectra, and it is consistent with observations of the diffuse gamma ray background.

Kazanas, D.↗

Insights on the structure and properties of sodium iron phosphate glasses from molecular dynamics simulations

Iron phosphate glasses are promising nuclear waste forms while more detailed understanding of their structures and structure-property relations are still needed to better design waste glass compositions. In this work we report studies of three series of sodium iron phosphate (NFP) glasses: 60P 2 O 5 -(40-x)Fe 2 O 3 -xNa 2 O (x = 0→40), (100–2x)P 2 O 5 -xFe 2 O 3 -xNa 2 O (x = 5→17.5) and one with different iron redox ratio, to understand the composition as well as the iron redox effects on the structure and properties of these glasses using molecular dynamics simulations with effective two-body and three-body potentials. Structural analyses, including pair distribution function, bond angle distribution, Q n distribution, and polyhedral connectivity, were performed to obtain in-depth information on short-range and medium-range structural features. The P-O pair distributions showed a first peak splitting with phosphorus-bridging and non-bridging oxygen contributions. This and the average P-O and other cation-oxygen bond distances are in excellent agreement with experiments. The coordination number of P 5+ remained four while that of Fe 3+ increased from 4.30 to 4.72 with decreasing Fe/Na ratio. Polyhedral linkage analysis showed [PO 4 ] units linked with [PO 4 ] and [FeO x ] through corner-sharing while the [PO 4 ]-[FeO x ] linkages become dominant for compositions with Fe 2 O 3 larger than 15 mol%. The effect of iron redox ratio on the structure of NFP glasses was also studied and it was found that bond lengths and coordination numbers were not strongly affected, while the reduction of iron introduced higher network distortions, as evident by O-P-O bond angle and Q n distribution. The glass transition temperature (T g ) showed a monotonic increase with Fe 2 O 3 in the first series, in good agreement with experiments, while those of the second series showed a maximum at P 2 O 5 = 82 mol%. Here, calculated elastic moduli were found to increase with Fe 2 O 3 in the first glass series, which was be explained by the increase of network connectivity, while those of the second series decrease with Fe 2 O 3 due to decrease of P 2 O 5 .

36 MATERIALS SCIENCE↗

An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system

The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Estimating Electron Temperature and Density Using Van Allen Probe Data: Typical Behavior of Energetic Electrons in the Inner Magnetosphere

Abstract The Earth's inner magnetosphere contains multiple electron populations influenced by different factors. The cold electrons of the plasmasphere, warm plasma that contributes to the ring current, and the relativistic plasma of the radiation belts often seem to behave independently. Using omni‐directional flux and energy measurements from the HOPE and Magnetic Electron Ion Spectrometer instruments aboard the Van Allen Probes, we provide a detailed density and temperature description of the inner magnetosphere, offering a comprehensive statistical analysis of the entire Van Allen Probe era. While number density and temperature data at geosynchronous orbit are available, this study focuses on the warm plasma in the inner magnetosphere . Values of density and temperature are extracted by fitting energy and phase space density to obtain the distribution function. The fitted distributions are related to the zeroth and second moments to estimate the number density and temperature. Analysis has indicated that a two Maxwellian fit is sufficient over a wide range of and that there are two independent plasma populations. The more energetic population has a median number density of approximately and a temperature of around 130 keV, with a temperature peak observed between L * = 4 and L * = 4.5. This population is relatively uniform in magnetic local time (MLT). In contrast, the less energetic warm electron population has a median number density of about and a temperature of 7.4 keV. Strong statistical trends in density and temperature across both L * and MLT are presented, along with potential sources driving these variations.

58 GEOSCIENCES↗

Temperature-induced gradual polymorphic structural crossover in liquid indium, tin, and antimony

Temperature-induced polymorphic structural crossover in metallic liquid indium (In), tin (Sn), and antimony (Sb) is investigated by in situ high-energy x-ray diffraction and ab initio molecular dynamics simulations. The results demonstrate the existence of a temperature-induced reversible and gradual structural crossover in these three liquids, reflecting from both “static” structures, e.g., peak positions of the structure factor and pair distribution function, bond angle distribution, number of tetrahedra, cluster connection, free volume, and dynamical behaviors, including the relaxation time, diffusion coefficient, and non-Gaussian parameter. Further, the main difference from liquid In, Sn, to Sb is that not only the large-sized free volume but also the anharmonicity and the dynamical heterogeneity increase significantly with the appearance of strong covalent bonds. These findings will deepen the understanding of liquid structure and the liquid-liquid transition in metallic liquids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurements of jet cross-section ratios in 13 TeV proton-proton collisions with ATLAS

Measurements of jet cross-section ratios between inclusive bins of jet multiplicity are performed in 140 fb -1 of proton-proton collisions with √s = 13 TeV center-of-mass energy, recorded with the ATLAS detector at CERN’s Large Hadron Collider. These ratios are constructed from double-differential cross-section measurements that are made in bins of jet multiplicity and other observables that are sensitive the energy scale and angular distribution of radiation due to the strong interaction in the final state. Additionally, the scalar sum of the two leading jets’ transverse momenta is measured triple differentially, in bins of the third jet’s transverse momentum and of jet multiplicity. These measurements are unfolded to account for acceptance and detector-related effects. The measured distributions are used to construct ratios of the inclusive jet-multiplicity bins, which have been shown to be sensitive to the strong coupling Formula Presented while being less sensitive than other observables to systematic uncertainties and parton distribution functions. The measured distributions are compared with state-of-the-art QCD calculations, including next-to-next-to-leading-order predictions for two- and three-jet events. These predictions are generally found to model the data well and perform best in bins with a modest requirement on the third jet’s transverse momentum. Significant differences between data and Monte Carlo predictions are observed in events with large rapidity gaps and invariant masses of the leading jet pair. Studies leading to reduced jet energy scale uncertainties significantly improve the precision of this work and are documented herein.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On the statistics of correlator outputs

Characteristic function, probability density function, and probability distribution function computations for correlation detector output, and correlator output statistics

Wintz, P. A.↗

On the probability of cycle-slipping in first-order phase-locked loops.

The first-passage time boundary value problem for first-order phase-locked loops (PLL) is analyzed, and spectral representations are developed for the probability density function (pdf), the distribution function, and the moments of the first time to passage (or cycle-slip). For the sinusoidal PLL, an asymptotic formula, that is surprisingly accurate even at low loop SNR's and large frequency offsets, is obtained for the pdf of the time to cycle-slip, in terms of the mean time to slip.

La Frieda, J. R.↗

Properties of the Io plasma torus inferred from Voyager EUV data

A physical model for the Io plasma torus is constructed to explain the EUV radiative emission observed by the Voyager UV spectrometer. Electron impact excitation rate coefficients for electronic transitions of S III, S IV, O II and O III are calculated by the method of distorted waves (Davis, Kepple, and Blaha, 1976); these coefficients account for the asymmetric shape of the 686 A feature. It is concluded that the electron gas must have a distribution function with a non-Maxwellian tail. An approximate representation of the distribution function as two temperature components requires a cold component of 3.5-4 eV and density of 2000 per cu cm and a hot component of about 100 eV and density of 50-100 per cu cm to satisfy observational constraints.

Strobel, D. F.↗

Martian Crater Size Distributions and Terrain Age

The crater size/frequency distributions of large ( 8 km) craters on the Moon and terrestrial planets display two very different curves representing two crater populations. The heavily cratered regions of the Moon, Mercury, and Mars show the same highly structured curve which cannot be represented by a single slope distribution function. In contrast, the lunar post mare crater population has a size/frequency distribution which differs significantly from that in the highlands over the same diameter range, and can be represented by a single-slope distribution function of -2.8 differential. On areas of martian lightly cratered northern plains, the crater population is essentially identical to that of the post mare population. This indicates that the same two families of impacting objects were responsible for the cratering records on both Moon and Mars. The thickness of mantling material varies among the various plains units, and can be calculated from the depth/diameter scaling relations for martian craters.

Nadine G Barlow↗

A multiple scattering theory for EM wave propagation in a dense random medium

For a dense medium of randomly distributed scatterers an integral formulation for the total coherent field has been developed. This formulation accounts for the multiple scattering of electromagnetic waves including both the twoand three-particle terms. It is shown that under the Markovian assumption the total coherent field and the effective field have the same effective wave number. As an illustration of this theory, the effective wave number and the extinction coefficient are derived in terms of the polarizability tensor and the pair distribution function for randomly distributed small spherical scatterers. It is found that the contribution of the three-particle term increases with the particle size, the volume fraction, the frequency and the permittivity of the particle. This increase is more significant with frequency and particle size than with other parameters.

Karam, M. A.↗

Particle distributions in a two-dimensional reconnection field geometry

In order to understand space phase density variations near the magnetopause and in the geomagnetic tail, particle trajectories in a two-dimensional steady state reconnection field geometry are studied by tracing particles from a spacecraft position back in time. The distribution function at the spacecraft is calculated using Liouville's theorem and known distribution functions in the source regions. It is found, in particular, that multiple ion beams in the plasma sheet boundary layer are a natural consequence of reconnection. In the center of the central plasma sheet, the omnidirectional flux of keV ions should be at a minimum. The depth of this minimum is very large near the neutral point but decreases as the distance from the neutral point increases.

Gurran, D. B.↗

A self-consistent model for the particles and fields upstream of an outgassing comet. I - Gyrotropic and isotropic ion distributions

The magnetic field and solar wind flow parameters in the vicinity of a weakly outgassing comet are determined using a self-consistent model which treats the cometary ions kinetically. Two different assumptions are made concerning the cometary ion distribution function; the pickup cometary ions from either a velocity space gyrotropic ring distribution or a velocity space isotropic shell distribution in the solar wind frame of reference. The individual currents due to the solar wind and cometary ions, which determine the magnetic field, are calculated in each case. Using theoretically determined parameters which reflect the conditions expected to comet Kopff (a typical short-period comet and one possible target for the future CRAF/Cassini mission) for various heliocentric distances, the way in which global field and flow characteristics are influenced by the nature of the cometary ion distribution function is determined.

Flammer, K. R.↗

Probabilistic material degradation model for aerospace materials subjected to high temperature, mechanical and thermal fatigue, and creep

A probabilistic general material strength degradation model has been developed for structural components of aerospace propulsion systems subjected to diverse random effects. The model has been implemented in two FORTRAN programs, PROMISS (Probabilistic Material Strength Simulator) and PROMISC (Probabilistic Material Strength Calibrator). PROMISS calculates the random lifetime strength of an aerospace propulsion component due to as many as eighteen diverse random effects. Results are presented in the form of probability density functions and cumulative distribution functions of lifetime strength. PROMISC calibrates the model by calculating the values of empirical material constants.

Boyce, L.↗

Wave generation by contaminant ions

Investigations dealing with the generation of waves by contaminant ions are reported. The studies included the properties of the velocity distribution function of such ions. It was found that it takes about one ion-cyclotron period for the distribution function to transform from a beam to a ring distribution. A linear instability analysis was performed to examine the possibility of wave excitation by an ion beam streaming perpendicular to the ambient magnetic field in an ionospheric type of plasma. A simulation code was developed to study the nonlinear behavior of the waves excited by beam and ring types of ion distributions. The code treats both electrons and ions as particles kinetically to also factor in the possible involvement of lower-hybrid waves.

Singh, N.↗