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 397 records · Page 22

Microscopic Simulation and Macroscopic Modeling for Thermal and Chemical Non-Equilibrium

This paper deals with the accurate microscopic simulation and macroscopic modeling of extreme non-equilibrium phenomena, such as encountered during hypersonic entry into a planetary atmosphere. The state-to-state microscopic equations involving internal excitation, de-excitation, dissociation, and recombination of nitrogen molecules due to collisions with nitrogen atoms are solved time-accurately. Strategies to increase the numerical efficiency are discussed. The problem is then modeled using a few macroscopic variables. The model is based on reconstructions of the state distribution function using the maximum entropy principle. The internal energy space is subdivided into multiple groups in order to better describe the non-equilibrium gases. The method of weighted residuals is applied to the microscopic equations to obtain macroscopic moment equations and rate coefficients. The modeling is completely physics-based, and its accuracy depends only on the assumed expression of the state distribution function and the number of groups used. The model makes no assumption at the microscopic level, and all possible collisional and radiative processes are allowed. The model is applicable to both atoms and molecules and their ions. Several limiting cases are presented to show that the model recovers the classical twotemperature models if all states are in one group and the model reduces to the microscopic equations if each group contains only one state. Numerical examples and model validations are carried out for both the uniform and linear distributions. Results show that the original over nine thousand microscopic equations can be reduced to 2 macroscopic equations using 1 to 5 groups with excellent agreement. The computer time is decreased from 18 hours to less than 1 second.

Modeling↗

The stability of perfect elliptic disks. 1: The maximum streaming case

Self-consistent distribution functions are constructed for two-dimensional perfect elliptic disks (for which the potential is exactly integrable) in the limit of maximum streaming; these are tested for stability by N-body integration. To obtain a discrete representation for each model, simulated annealing is used to choose a set of orbits which sample the distribution function and reproduce the required density profile while carrying the greatest possible amount of angular momentum. A quiet start technique is developed to place particles on these orbits uniformly in action-angle space, making the initial conditions as smooth as possible. The roundest models exhibit spiral instabilities similar to those of cold axisymmetric disks; the most elongated models show bending instabilities like those seen in prolate systems. Between these extremes, there is a range of axial ratios 0.25 approximately less than b/a approximately less than 0.6 within which these models appear to be stable. All the methods developed in this investigation can easily be extended to integrable potentials in three dimensions.

Levine, Stephen E.↗

The electron Boltzmann equation in a plasma generated by fission fragments

A Boltzmann equation formulation is presented for the determination of the electron distribution function in a plasma generated by fission fragments. The formulation takes into consideration ambipolar diffusion, elastic and inelastic collisions, recombination and ionization, and allows for the fact that the primary electrons are not monoenergetic. Calculations for He in a tube coated with fissionable material show that, over a wide pressure and neutron flux range, the distribution function is non-Maxwellian, but the electrons are essentially thermal. Moreover, about a third of the energy of the primary electrons is transferred into the inelastic levels of He. This fraction of energy transfer is almost independent of pressure and neutron flux but increases sharply in the presence of a sustainer electric field.

Hassan, H. A.↗

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↗

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.↗

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.↗

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

NASA is concerned with the health risks to astronauts, particularly those risks related to radiation exposure. 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 (A) peripheral leukocyte distribution; (B) plasma cytokine levels; (C) cytokine production profiles following whole blood stimulation of either T cells or monocytes.

Crucian, Brian E.↗

Eight-moment approximation solar wind models

Heat conduction from the corona is important in the solar wind energy budget. Until now all hydrodynamic solar wind models have been using the collisionally dominated gas approximation for the heat conductive flux. Observations of the solar wind show particle distribution functions which deviate significantly from a Maxwellian, and it is clear that the solar wind plasma is far from collisionally dominated. We have developed a numerical model for the solar wind which solves the full equation for the heat conductive flux together with the conservation equations for mass, momentum, and energy. The equations are obtained by taking moments of the Boltzmann equation, using an 8-moment approximation for the distribution function. For low-density solar winds the 8-moment approximation models give results which differ significantly from the results obtained in models assuming the gas to be collisionally dominated. The two models give more or less the same results in high density solar winds.

Olsen, Espen Lyngdal↗