Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution 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 685 records · Page 38

Solution of an integral equation encountered in studies on radiative transfer in completely absorbing leaf canopies

The leaf normal distribution function is related to the extinction coefficient of direct solar radiation by a Fredholm integral equation of the first kind. The integral equation is solved by the method of constrained least squares. The solutions are stable and correct only in extreme cases and, in general, it is not possible to recover an arbitrary leaf normal distribution function.

Myneni, Ranga B.↗

A predicted energetic ion signature of a neutral line in the geomagnetic tail

Under certain assumptions one can predict the formation of a magnetic neutral line on theoretical grounds, generally as a result of an unstable current sheet. Neutral points also seen to form spontaneously in a wide variety of plasma simulatioons, from MHD to fully kinetic. Hence, such magnetic null regions are central to some models of magnetotail energization and magnetic substorms. Yet the observational evidence for such regions remains controversial. This paper presents a new signature of a neutral line: the field-aligned ion beam characteristic of current sheet acceleration becomes dispersed into a 'ridge' in the velocity space distribution function as a neutral line is approached. The results were obtained using single-particle simulations in a fixed magnetic field model. The origin of the ridge is shown to be due to initially earthward going ions which pass through the neutral line region. These ions come from a smaller part of the initially antiearthward flowing distribution, generating the depletion of the distribution function below the ridge.

Martin, R. F., Jr.↗

Combined Effects of Gamma Radiation and High Dietary Iron on Peripheral Leukocyte Distribution and Function

Both radiation and increased iron stores can independently increase oxidative damage, resulting in protein, lipid and DNA oxidation. Oxidative stress increases the risk of many health problems including cancer, cataracts, and heart disease. This study, a subset of a larger interdisciplinary investigation of the combined effect of iron overload on sensitivity to radiation injury, monitored immune parameters in the peripheral blood of rats subjected to gamma radiation, high dietary iron or both. Specific immune measures consisted of: (1) peripheral leukocyte distribution, (2) plasma cytokine levels and (3) cytokine production profiles following whole blood mitogenic stimulation

Crucian, Brian E.↗

Higher order Compton-Getting anisotropies.

The Compton-Getting effect relates the distribution function of a particle population in a given frame of reference to that in another frame, moving at constant velocity V relative to the first. A study is made of the effects that arise when, contrary to previous assumptions, second-order terms in (V/v), where v is the particle velocity, are retained, and when the given distribution function is not isotropic. The importance of the extra terms is then examined in the specific case of low-energy solar protons in interplanetary space.

Balogh, A.↗

The nature of cometary dust as determined from infrared observations

The infrared measurements of comets, the compositional information available from interplanetary dust particles (IDPs), and the recent results of flybys to Comet Halley can help in restricting the nature and composition of cometary dust models (c.f., Proceedings of the 20th ESLAB Symposium on Exploration of Halley's Comet, 1986). Researchers tried to incorporate some of these results into a coherent model to account for the observed cometary infrared emission. The presence of 10 and 3.4 micron features in Comet Halley (c.f. Bregman et al. 1987; Wickramasinghe and Allen 1986) indicated the presence of at least two components in the grain material, namely silicates and some form of amorphous carbon. These two components could reside in separate grains or may be parts of composite particles. Both these cases have been considered (see Krishna Swamy el a. 1988a, 1988b). In the absence of refractive index data for cometary analogs, the authors used the optical constants of olivine-rich lunar material 12009.48 (Perry et al. 1972) for the infrared region and that of alpha:C-H film for amorphous carbon (angus et al. 1986). For the visible region, a value of m = 1.38-0.39i was used for the silicates, and values published by Arakawa et al. (1985) were used for the amorphous carbon. These materials should give a representative behavior of the expected results. The model results were compared to observational data. The strength of the 3.4 micron and 10 micron features relative to the adjacent continuum, as well as the slope of the continuum between 2500 and 1250 cm(exp -1) (4 to 8 microns), were used as criteria for comparison. Model calculations with alpha approx. equals -3.5, and also the size distribution function inferred for Comet Halley, with a mass fraction (X) of silicate to amorphous carbon grains of about 40 to 1 can fit the data. A good match is obtained for the infrared spectra of Comets Halley and West from a 40 to 1 mixture of silicate and amorphous carbon grains with a a(exp -3.5) size distribution function. The results are consistent with compositional constraints provided by interplanetary dust particles (IPDs) and Halley flyby data. The variation of grain temperature with heliocentric distance appears to account for the major changes observed in cometary spectra.

Swamy, K. S. Krishna↗

Statistical description of coalescing magnetic islands via magnetic reconnection

The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in two dimensions. We solve our island kinetic equation numerically for three different types of initial distribution: Dirac delta, Gaussian and power-law distributions. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as $\tilde {t}^{-1}$ , the number of islands decreases as $\tilde {t}^{-1}$ and the averaged area of islands grows as $\tilde {t}$ , where $\tilde {t}$ is the time normalised to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied. Finally, we discuss the underlying connection of our island-merger models to the (self-similar) decay of magnetohydrodynamic turbulence.

Physics↗

Augmenting machine learning of energy landscapes with local structural information

We present a machine learning approach for accurately predicting formation energies of binary compounds in the context of crystal structure predictions. The success of any machine learning model depends significantly on the choice of representation used to encode the relevant physical information into machine-learnable data. We test different representation schemes based on partial radial and angular distribution functions (RDF+ADF) on Al–Ni and Cd–Te structures generated using our genetic algorithm for structure prediction. We observe a remarkable improvement in predictive accuracy upon transitioning from global to atom-centered representations, resulting in a threefold decrease in prediction errors. We show that a support vector regression model using a combination of atomic radial and angular distribution functions performs best at the formation energy prediction task, providing small root mean squared errors of 3.9 meV/atom and 10.9 meV/atom for Al–Ni and Cd–Te, respectively. We test the performance of our models against common traditional descriptors and find that RDF- and ADF-based representations significantly outperform many of those in the prediction of formation energies. The high accuracy of predictions makes our machine learning models great candidates for the exploration of energy landscapes.

Honrao, Shreyas J. (ORCID:0000000166155257)↗

Distinguishing Quarks and Gluons in Pion and Kaon PDFs

The leading-twist parton distribution functions of the pion and kwon are calculated for the first time using a rainbow-ladder truncation of QCD's Dyson-Schwinger equations (DSEs) that self-consistently sums all planar diagrams. The nonperturbative gluon dressing of the quarks is thereby correctly accounted for, which in practice means solving the inhomogeneous Bethe-Salpeter equation (BSE) for the quark operator that defines the spin-independent quark distribution functions. An immediate consequence of using this dressed vertex is that gluons carry 35% of the pion's and 30% of the kaon's light-cone momentum, with the remaining momentum carried by the quarks. The scale associated with these DSE results is μ 0 = 0.78 GeV. The gluon effects generated by the inhomogeneous BSE are inherently nonperturbative and cannot be mimicked by the perturbative QCD evolution equations. A key consequence of this gluon dressing is that the valence quarks have reduced support at low-to-intermediate x, where the gluons dominate, and increased support at large x. As a result, our DSE calculation of the pion's valence quark distribution is in excellent agreement with the Conway et al. pion-induced Drell-Yan data, but nevertheless exhibits the q π (x) similar or equal to (1 - x) 2 behavior as x -> 1 predicted by perturbative QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fast resonance decays in nuclear collisions

Here, we present a new method to calculate the final decay spectrum of direct resonance decays directly from hydrodynamic fields on a freeze-out surface. The method is based on identifying components of the final particle spectrum that transform in an irreducible way under rotations in the fluid-restframe. Corresponding distribution functions can be pre-computed including all resonance decays. Just a few of easily tabulated scalar functions then determine the Lorentz invariant decay spectrum from each space-time point, and simple integrals of these scalar functions over the freeze-out surface determine the final decay products. This by-passes numerically costly event-by-event calculations of the intermediate resonances. The method is of considerable practical use for making realistic data to model comparisons of the identified particle yields and flow harmonics, and for studying the viscous corrections to the freeze-out distribution function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Anisotropies in the fluxes of Pioneer 10 protons

One-hour-averaged fluxes of 1.8- to 2.15-MeV protons observed by the LET2 detector on Pioneer 10 on the inbound trajectory showed anisotropies attributable to corotation of Jupiter's magnetodisc only when Pioneer was near the dipole equator. Most of the time the anisotropy greatly exceeded the value expected from corotation. Gradients in the distribution function can be used to account for this excess anisotropy, but the amount of gradient required is unacceptably large by 1-2 orders of magnitude. If they were taken as real, these gradients would predict almost complete disappearance of these protons from Jupiter's magnetosphere in a matter of hours. The remedy is to introduce into the model of the distribution function proton flow along field lines away from the equator into both the southern and the northern hemisphere. The parallel flux at the southernmost latitudes reached by Pioneer can reach 25% of the product of proton density and velocity, i.e., 25% of the maximum possible.

Northrop, T. G.↗

Application of Powder Diffraction Methods to the Analysis of Short- and Long-Range Atomic Order in Nanocrystalline Diamond and SiC: The Concept of the Apparent Lattice Parameter (alp)

Two methods of the analysis of powder diffraction patterns of diamond and SiC nanocrystals are presented: (a) examination of changes of the lattice parameters with diffraction vector Q ('apparent lattice parameter', alp) which refers to Bragg scattering, and (b), examination of changes of inter-atomic distances based on the analysis of the atomic Pair Distribution Function, PDF. Application of these methods was studied based on the theoretical diffraction patterns computed for models of nanocrystals having (i) a perfect crystal lattice, and (ii), a core-shell structure, i.e. constituting a two-phase system. The models are defined by the lattice parameter of the grain core, thickness of the surface shell, and the magnitude and distribution of the strain field in the shell. X-ray and neutron experimental diffraction data of nanocrystalline SiC and diamond powders of the grain diameter from 4 nm up to micrometers were used. The effects of the internal pressure and strain at the grain surface on the structure are discussed based on the experimentally determined dependence of the alp values on the Q-vector, and changes of the interatomic distances with the grain size determined experimentally by the atomic Pair Distribution Function (PDF) analysis. The experimental results lend a strong support to the concept of a two-phase, core and the surface shell structure of nanocrystalline diamond and SiC.

Palosz, B.↗

Interaction with the lower ionosphere of electromagnetic pulses from lightning - Heating, attachment, and ionization

A Boltzmann formulation of the electron distribution function and Maxwell's equations for the EM fields are used to simulate the interaction of lightning radiated EM pulses with the lower ionosphere. Ionization and dissociative attachment induced by the heated electrons cause significant changes in the local electron density, N(e). Due to 'slow' field changes of typical lightning EM pulses over time scales of tens of microsec, the distribution function follows the quasi-equilibrium solution of the Boltzmann equation in the altitude range of interest (70 to 100 km). The EM pulse is simulated as a planar 100 microsec long single period oscillation of a 10 kHz wave injected at 70 km. Under nighttime conditions, individual pulses of intensity 10-20 V/m (normalized to 100 km horizontal distance) produce changes in N(e) of 1-30 percent while a sequence of pulses leads to strong modification of N(e) at altitudes less than 95 km. The N(e) changes produce a 'sharpening' of the lower ionospheric boundary by causing a reduction in electron density at 75-85 km (due to attachment) and a substantial increase at 85-95 km (due to ionization) (e.g., the scale height decreases by a factor of about 2 at about 85 km for a single 20 V/m EM pulse). No substantial N(e) changes occur during daytime.

Taranenko, Y. N.↗

Probabilistic micromechanics and macromechanics of polymer matrix composites

A probabilistic evaluation of an eight ply graphite-epoxy quasi-isotropic laminate was completed using the Integrated Composite Analyzer (ICAN) in conjunction with Monte Carlo simulation and Fast Probability Integration (FPI) techniques. Probabilistic input included fiber and matrix properties, fiber misalignment, fiber volume ratio, void volume ratio, ply thickness and ply layup angle. Cumulative distribution functions (CDFs) for select laminate properties are given. To reduce the number of simulations, a Fast Probability Integration (FPI) technique was used to generate CDFs for the select properties in the absence of fiber misalignment. These CDFs were compared to a second Monte Carlo simulation done without fiber misalignment effects. It was found that FPI requires fewer simulations to obtain the cumulative distribution functions as opposed to Monte Carlo simulation techniques. Furthermore, FPI provides valuable information regarding the sensitivities of composite properties to the constituent properties, fiber volume ratio and void volume ratio.

Mase, G. T.↗

Probabilistic micromechanics and macromechanics of polymer matrix composites

A probabilistic evaluation of an eight ply graphite-epoxy quasi-isotropic laminate was completed using the Integrated Composite Analyzer (ICAN) in conjunction with Monte Carlo simulation and Fast Probability Integration (FPI) techniques. Probabilistic input included fiber and matrix properties, fiber misalignment, fiber volume ratio, void volume ratio, ply thickness and ply layup angle. Cumulative distribution functions (CDFs) for select laminate properties are given. To reduce the number of simulations, a Fast Probability Integration (FPI) technique was used to generate CDFs for the select properties in the absence of fiber misalignment. These CDFs were compared to a second Monte Carlo simulation done without fiber misalignment effects. It was found that FPI requires fewer simulations to obtain the cumulative distribution functions as opposed to Monte Carlo simulation techniques. Furthermore, FPI provides valuable information regarding the sensitivities of composite properties to the constituent properties, fiber volume ratio and void volume ratio.

Mase, G. T.↗

Location and characteristics of the reconnection X-line deduced from low-altitude satellite and radar observations

We present an analysis of a cusp ion step observed between two poleward-moving events of enhanced ionospheric electron temperature. From the computed variation of the reconnection rate and the onset times of the associated ionospheric events, the distance between the satellite and the X-line can be estimated, but with a large uncertainty due to that in the determination of the low-energy cut-off of the ion velocity distribution function, f(E). Nevertheless, analysis of the time series f(t) shows the reconnection site to be on the dayside magnetopause, consistent with the pulsating cusp model, and the best estimate of the X-line location is 13 R(E) from the satellite. The ion precipitation is used to reconstruct the field-parallel part of the Cowley-D ion distribution function injected into the open low latitude boundary layer (LLBL) in the vicinity of the X-line. From this the Alfven speed, plasma density, magnetic field, parallel ion temperature, and flow velocity of the magnetosheath near the X-line can be derived.

Lockwood, M.↗

Cometary pick-up ions observed near Giacobini-Zinner

The number and energy density of cometary water-group ions observed near Comet Giacobini-Zinner are derived using the rest-frame distribution functions. The data reveal that density profiles of inbound and outbound passes and their shape correlate with pick-up ion production model predictions. The lose rate and production rate of water-group cometary molecules calculated from predicted and measured density profiles are 2 x 10 to the -6th/sec and 2.6 x 10 to the 28th/sec respectively. The shapes of the distribution functions are examined to study the solar wind/cometary ions interaction process.

Gloeckler, G.↗

A small-box approach to the local crystal structure of Y 3 NbO 7

The local crystal structure of Y 3 NbO 7 was probed using neutron total scattering. Five different structural models were fit to the experimental pair distribution function of this material up to 15 Å. Fits of defect fluorite (Fm$\overline{3}$m), pyrochlore (Fd$\overline{3}$m), and orthorhombic C weberite models (C222 1 ) were attempted based on prior structural investigations of Y 3 NbO 7 and related fluorite-based materials. Orthorhombic P weberite (P222 1 ) and monoclinic weberite models (P112 1 ) were derived via symmetry reduction of the orthorhombic C weberite structure. All models except monoclinic weberite failed to reproduce the experimental pair distribution function. Misfits were clearly visible for peaks arising from atom pairs within the first coordination shell (defect fluorite and pyrochlore) and beyond (orthorhombic C and P weberite). By contrast, monoclinic weberite (γ ≈ 91.8°) provided an excellent fit, with no major misfits observed. This unit cell features distorted YO 8 dodecahedra, YO 7 monocapped octahedra, and NbO 6 octahedra as building blocks. Comparison between the topologies of the cationic and anionic substructures in monoclinic weberite and its parent orthorhombic C structure revealed displacements of the metal atoms along the [100] direction and a less compact spatial distribution of oxygen atoms along the a axis. In conclusion, the monoclinic weberite model presented herein may serve as a starting point to revisit structural modeling of fluorite-based materials with challenging local structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Magnetopause modeling - Flux transfer events and magnetosheath quasi-trapped distributions

Three-dimensional distribution functions for energetic ions are studied numerically in the magnetosphere, through the magnetopause, and in the magnetosheath using a simple one-dimensional quasi-static model and ISEE 1 magnetopause crossing data for November 10, 1977. Quasi-trapped populations in the magnetosheath observed near flux transfer events (FTEs) are investigated, and it is shown that the population in the sheath appears to sandwich the FTE distributions. These quasi-trapped distributions are due to slow, large pitch angle, outward moving particles left behind by the outward rush of the ions more field-aligned at the time the flux was opened. It is found that sheath convective flows can map along the connected flux tube without drastically changing the distribution function, and results suggest that localized tangential fields above the upper limit may exist.

Speiser, T. W.↗